基于特征理论的二维可压缩流动的二阶拉氏算法

提出一种求解二维拉氏可压缩流体力学方程的中心型二阶精度有限体积方法.利用特征理论构造网格节点处的局部近似演化算子,算子用来求解网格节点处的速度及压力,利用这些物理量更新节点位置及计算网格界面通量.通过结合一定的重构方案,该方法达到时、空二阶精度,并且形式简单、计算量小,适用于结构网格与非结构网格.典型数值实验表明,本文格式具有良好的收敛性、对称性及鲁棒性,且能自然地求解多物质流动问题.     

热传导方程的一类无网格方法

构造求解热传导方程的一类无网格方法,只要选择好每个节点的适当的邻点集合,便可利用节点信息顺利进行计算.作为特殊情形,也可在各种结构或非结构网格上进行计算.在矩形或均匀平行四边形网格上进行计算时具有二阶精度,当在任意的不规则四边形或三角形网格上计算时仍然是守恒的和相容的,且至少具有一阶精度.作为数值试验,将该方法用于在不规则四边形网格上及四边形与三角形混合网格上求解二维非线性抛物型方程,并在不规则四边形网格上求解二维三温辐射热传导方程,均获得了较为理想的数值结果.     

求解多维欧拉方程的二阶非结构网格混合旋转Riemann求解器

将基于旋转近似Riemann求解器的二阶精度迎风型有限体积方法推广到非结构网格,采用基于网格中心的有限体积法,梯度的计算采用基于节点的方法引入更多的控制体模板,限制器的构造采用与非结构化网格相适应的形式.在求解Riemann问题时,沿具有一定物理意义的两个迎风方向,即控制体界面两侧速度差矢量方向及与之正交的方向.能够完全消除基于Riemann求解器的通量差分裂格式存在的激波不稳定或"红斑"现象.为减小计算量,采用HLL和Roe FDS混合旋转格式.     

三维非匹配网格上的扩散方程数值求解

将子网格剖分的支撑算子方法,拓展应用于三维非匹配网格上的扩散方程求解.算例表明该方法在正交非匹配网格上能够精确获得线性解;在一般非匹配网格上可以达到二阶精度;在求解曲面网格和节点不共面网格时,精度比平面近似的方法要高,也可以达到2阶精度,同时也适合求解含有物质界面的混合介质网格.     

非结构网格上Level Set方程的间断有限元解法

针对间断有限元弱形式难于求解可压缩流场中Level Set方程的问题提出间断有限元强形式,从而在统-框架下解决Level Set方程在可压缩与不可压缩流场中的求解问题.通过非结构网格上采用Legendre-Gauss-Lobatto节点构造基函数,在复杂区域上可以达到任意高阶的精度.将若干-、二、三维算例与已有文献或解析解比较,验证方法追踪自由界面的有效性.结果表明,该方法适合各种情形下Level Set方程求解,易于在复杂区域的非结构网格上实施,精度高、分辨率高且具有高质量守恒性,既能避免重新初始化过程又方便向高维扩展.     

求解可压缩流的高精度非结构网格WENO有限体积法

提出-种基于最小二乘重构和WENO限制器的非结构网格高精度有限体积方法.用中心网格的某些邻居网格建立重构多项式,给出-定的原则搜索和存储足够多的邻居网格以建立重构多项式,采用最小二乘法求解重构多项式的系数.用-种通用的方法控制重构邻居个数,以减少存储和计算,采用WENO限制器和旋转Riemann求解器以达到统-的高精度并且抑制守恒律方程求解中的非物理振荡.为检验上述算法,以基于节点的梯度重构,Bath and Jesperson限制器的二阶算法为基准,给出三阶和四阶格式与二阶格式以及高阶格式若干经典算例计算结果的对比和分析.     

带源参数的热传导反问题的无网格方法

利用无网格有限点法求带有源参数的一维热传导反问题，推导了相应的离散方程.与其他基于网格的方法相比，有限点法采用移动最小二乘法构造形函数，只需要节点信息，不需要划分网格，用配点法离散求解方程，可以直接施加边界条件，不需要在区域内部求积分，减小了计算量.用有限点法求解热传导反问题具有数值实现简单、计算量小、可以任意布置节点等优点.最后通过算例验证了该方法的有效性.     

定时间步长变坐标步长差分求解单相Stefan问题

对单相Stefan问题提出了一种定时间步长、变坐标步长的差分求解方法.在固定时间步长内,以计算得到的移动界面位置作为网格节点的坐标,前后界面位置之差为空间步长,逐步地自动形成网格的划分,计算这些节点处的温度,从而获得下一时刻的移动界面的位置.     

九点格式中网格边上切向流的计算——节点自适应加权格式

针对网格扭曲的不同情形,直接考虑网格边上切向流的离散.基于扩散方程法向流连续的条件,给出离散法向流的构造,导出扭曲网格上九点计算格式中网格边上离散切向流的表达式,从而推导出加权系数的计算公式,适应于各种扭曲的网格.数值结果表明,与九点格式中节点量简单加权的方法相比,基于网格边离散切向流的节点自适应加权九点格式的精度有明显改进,迭代求解次数减少,计算效率明显提高.     

多流管网格上九点格式的节点插值方法

多流管方法是二维多介质辐射流体力学数值模拟中一类常用的求解方法,它采用Lagrange-Euler混合型四边形网格,称为多流管网格。通常其网格品质高于一般的四边形网格。在这类网格上,可以利用网格特性对九点扩散格式中的节点插值方法进行改进。本文利用调和平均点和梯度离散构造的方法提出几种节点插值方法。并给出数值实验,说明现有应用程序中的节点插值方法损失精度,而新的节点插值方法能够使得九点格式在多流管网格上具有二阶精度。     

A new high-order discontinuous Galerkin spectral finite element method for Lagrangian gas dynamics in two-dimensions

This paper presents a new high-order cell-centered Lagrangian scheme for two-dimensional compressible flow. The scheme uses a fully Lagrangian form of the gas dynamics equations, which is a weakly hyperbolic system of conservation laws. The system of equations is discretized in the Lagrangian space by discontinuous Galerkin method using a spectral basis. The vertex velocities and the numerical fluxes through the cell interfaces are computed consistently in the Eulerian space by virtue of an improved nodal solver. The nodal solver uses the HLLC approximate Riemann solver to compute the velocities of the vertex. The time marching is implemented by a class of TVD Runge–Kutta type methods. A new HWENO (Hermite WENO) reconstruction algorithm is developed and used as limiters for RKDG methods to maintain compactness of RKDG methods. The scheme is conservative for the mass, momentum and total energy. It can maintain high-order accuracy both in space and time, obey the geometrical conservation law, and achieve at least second order accuracy on quadrilateral meshes. Results of some numerical tests are presented to demonstrate the accuracy and the robustness of the scheme.     

Cubic interaction vertex of higher-spin fields with external electromagnetic field

We fulfill the detailed analysis of coupling the charged bosonic higher-spin fields to external constant electromagnetic field in first order in external field strength. Cubic interaction vertex of arbitrary massive and massless bosonic higher-spin fields with external field is found. Construction is based on deformation of free Lagrangian and free gauge transformations by terms linear in electromagnetic field strength. In massive case a formulation with Stueckelberg fields is used. We begin with the most general form of deformations for Lagrangian and gauge transformations, admissible by Lorentz covariance and gauge invariance and containing some number of arbitrary coefficients, and require the gauge invariance of the deformed theory in first order in strength. It yields the equations for the coefficients which are exactly solved. As a result, the complete interacting Lagrangian of arbitrary bosonic higher-spin fields with constant electromagnetic field in first order in electromagnetic strength is obtained. Causality of massive spin-2 and spin-3 fields propagation in the corresponding electromagnetic background is proved.     

SOME DISCUSSIONS OF POSSIBLE MECHANISM OF Z0→l+l-γ

The possible mechanisms for anomalous decay Z⁰→l⁺l^-γ are discussed. The mechanism with anomaious ZZ_Y vertex is investipted indetail. It is pointed out that this vertex may appear via some effective Lagrangian which can be derived from an SU(2)×U (1) gauge invariant effective Lagrangian after the spontaneous symmetry breaking. The predictions based on this mechanism are consistent with recent experiments and can be verified experimentally. This mechanism is consistent with a composite model of Z⁰ with radius 10^-15 cm.     

Radiative corrections to low-energy ππ scattering

Radiative corrections to low-energy ππ scattering are calculated. A chiral SU(2) × SU(2) Lagrangian is used to describe low-energy scattering. This strong interaction Lagrangian is modified in the presence of the electromagnetic field in a gauge invariant way. It is found that photons do not probe the strong vertex. Corrections to the cross sections are calculated to O(α). All the ultra-violet divergences cancel, and the results are independent of cut-off parameters. Corrections to the Weinberg scattering lengths of the order of 1% are defined by removal of the Coulomb pole terms. The field theoretic calculation is compared with a potential theory approach. It is argued that the potential theory calculation of radiative corrections is more reliable.     

The quintic interaction vertex in light-cone gravity

We consider pure gravity in light-cone gauge and derive the complete quintic interaction vertex. Up to quartic order, the Kawai–Lewellen–Tye (KLT) relations can be made manifest at the level of the Einstein–Hilbert Lagrangian. The quintic interaction vertex represents an essential first step in further extending the off-shell validity of the KLT relations to higher order vertices.     

Face offsetting: A unified approach for explicit moving interfaces

Dynamic moving interfaces are central to many scientific, engineering, and graphics applications. In this paper, we introduce a novel method for moving surface meshes, called the face offsetting method, based on a generalized Huygens’ principle. Our method operates directly on a Lagrangian surface mesh, without requiring an Eulerian volume mesh. Unlike traditional Lagrangian methods, which move each vertex directly along an approximate normal or user-specified direction, our method propagates faces and then reconstructs vertices through an eigenvalue analysis locally at each vertex to resolve normal and tangential motion of the interface simultaneously. The method also includes techniques for ensuring the integrity of the surface as it evolves. Face offsetting provides a unified framework for various dynamic interface problems and delivers accurate physical solutions even in the presence of singularities and large curvatures. We present the theoretical foundation of our method, and also demonstrate its accuracy, efficiency, and flexibility for a number of benchmark problems and a real-world application.     

PCAC and the non-relativistic pion-nucleon vertex

The Foldy-Wouthuysen reduction technique is applied to the sigma model and pseudovector coupling pion-nucleon Lagrangians. In contrast to electrodynamics the requirement that the non-relativistic pion-nucleon Lagrangian have an approximate symmetry under the same global transformation as the relativistic theory does not uniquely fix the Foldy-Wouthuysen transformation. It is shown that the partially conserved axial-vector current relation is satisfied by the non-relativistic theory and the implications of this on the form of the effective pion-nucleon vertex are discussed.     

Multi-scale Godunov-type method for cell-centered discrete Lagrangian hydrodynamics

This work presents a multi-dimensional cell-centered unstructured finite volume scheme for the solution of multimaterial compressible fluid flows written in the Lagrangian formalism. This formulation is considered in the Arbitrary-Lagrangian–Eulerian (ALE) framework with the constraint that the mesh velocity and the fluid velocity coincide. The link between the vertex velocity and the fluid motion is obtained by a formulation of the momentum conservation on a class of multi-scale encased volumes around mesh vertices. The vertex velocity is derived with a nodal Riemann solver constructed in such a way that the mesh motion and the face fluxes are compatible. Finally, the resulting scheme conserves both momentum and total energy and, it satisfies a semi-discrete entropy inequality. The numerical results obtained for some classical 2D and 3D hydrodynamic test cases show the robustness and the accuracy of the proposed algorithm.     

The extended Nambu-Jona-Lasinio model with separable interaction: Low energy pion physics and pion-nucleon form factor

A Lagrangian formulation of the Nambu-Jona-Lasinio model with separable interaction is given. The electromagnetic interaction is introduced in a nonminimal way to the nonlocal quark current. Various choices of the vertex form factors characterizing the composite structure of mesons and baryon are investigated. We find that the physical observables depend very weakly on form factor shapes. We calculate the πNN form factor considering nucleon as a three-quark system.