Multiderivative Combined Dissipative Compact Scheme Satisfying Geometric Conservation Law III: Characteristic-Wise Hybrid Method

Based on newly developed weight-based smoothness detectors and non-linear interpolations designed to capture discontinuities for the multiderivative com-bined dissipative compact scheme(MDCS),hybrid linear and nonlinear interpolations are proposed to form hybrid MDCS.These detectors are derived from the weights used for the nonlinear interpolations and can provide suitable switches between the linear and the nonlinear schemes to realize the characteristics for the hybrid MDCS of capturing discontinuities and maintaining high resolution in the region without large discontinuities.To save computational cost,the nonlinear scheme with characteris-tic decomposition is only applied in the detected discontinuities region by specially designed hybrid strategy.Typical tests show that the hybrid MDCS is capable of cap-turing discontinuities and maintaining high resolution power for the smooth region at the same time.With the satisfaction of the geometric conservative law(GCL),the MDCS is further applied on curvilinear mesh to present its promising capability of handling pragmatic simulations.

High-order field theory and a weak Euler–Lagrange–Barut equation for classical relativistic particle-field systems

In both quantum and classical field systems,conservation laws such as the conservation of energy and momentum are widely regarded as fundamental properties.A broadly accepted approach to deriving conservation laws is built using Noether's method.However,this procedure is still unclear for relativistic particle-field systems where particles are regarded as classical world lines.In the present study,we establish a general manifestly covariant or geometric field theory for classical relativistic particle-field systems.In contrast to quantum systems,where particles are viewed as quantum fields,classical relativistic particle-field systems present specific challenges.These challenges arise from two sides.The first comes from the mass-shell constraint.To deal with the mass-shell constraint,the Euler–Lagrange–Barut(ELB)equation is used to determine the particle's world lines in the four-dimensional(4D)Minkowski space.Besides,the infinitesimal criterion,which is a differential equation in formal field theory,is reconstructed by an integro-differential form.The other difficulty is that fields and particles depend on heterogeneous manifolds.To overcome this challenge,we propose using a weak version of the ELB equation that allows us to connect local conservation laws and continuous symmetries in classical relativistic particle-field systems.By applying a weak ELB equation to classical relativistic particle-field systems,we can systematically derive local conservation laws by examining the underlying symmetries of the system.Our proposed approach provides a new perspective on understanding conservation laws in classical relativistic particle-field systems.

Developing Hybrid cell-edge and cell-node Dissipative Compact Scheme for Complex Geometry Flows

Developing high resolution finite difference scheme and enabling the use of this scheme on complex geometry are the aims of this study.High resolution has been achieved by Dissipative Compact Schemes(DCS),however,according to the recent research,applications of DCS on complex geometry may have serious problem for that the Geometric Conservation Law(GCL)is not satisfied,and this may cause numerical instability.To cope with this problem,a new scheme named Hybrid cell-edge and cell-node Dissipative Compact Scheme(HDCS)has been formulated.The formulation of the HDCS contains two steps.First,a new central compact scheme is formulated for the purpose of conveniently fulfilling the GCL,and then dissipation is added on the central scheme by high-order dissipative interpolation of cell-edge variables.The solutions of Euler and Navier-Stokes equations show that the HDCS can be applied successfully on complex geometry,while the DCS may suffer numerical instabilities.Moreover,high resolution of the HDCS may be observed in the test of scattering of acoustic waves by multiple cylinders.

High-Order and High Accurate CFD Methods and Their Applications for Complex Grid Problems

The purpose of this article is to summarize our recent progress in high-order and high accurate CFD methods for flow problems with complex grids as well as to discuss the engineering prospects in using these methods.Despite the rapid development of high-order algorithms in CFD,the applications of high-order and high accurate methods on complex configurations are still limited.One of the main reasons which hinder the widely applications of thesemethods is the complexity of grids.Many aspects which can be neglected for low-order schemes must be treated carefully for high-order ones when the configurations are complex.In order to implement highorder finite difference schemes on complex multi-block grids,the geometric conservation lawand block-interface conditions are discussed.A conservativemetricmethod is applied to calculate the grid derivatives,and a characteristic-based interface condition is employed to fulfil high-order multi-block computing.The fifth-order WCNS-E-5 proposed by Deng[9,10]is applied to simulate flows with complex grids,including a double-delta wing,a transonic airplane configuration,and a hypersonic X-38 configuration.The results in this paper and the references show pleasant prospects in engineering-oriented applications of high-order schemes.

A REMAPPING METHOD BASED ON MULTI-POINT FLUX CORNER TRANSPORT UPWIND ADVECTION ALGORITHM

A local remapping algorithm for scalar function on quadrilateral meshes is described. The remapper from a distorted grid to a rezoned grid is usually regarded as a conservative interpolation problem. The present paper introduces a pseudo time to transform the interpolation into an initial value problem on a moving grid, and construct a moving mesh method to solve it. The new feature of the algorithm is the introduction of multi- point information on each edge, which leads to the numerical flux consistent with grid node motion. During the procedure of deriving scheme, we illustrate a framework about how the algorithms on a rectangular mesh are easily generated to those on a moving mesh. The basic ideas include： （i） introducing coordinate transformation, which maps the irregular domain in physical space to a perfectly regular computational domain, and （ii） deriving finite volume methods in the physical domain, which can be viewed as a discretization of the transformed equation. The resulting scheme is second-order accurate, conservative and monotonicity preserving. Numerical examples are carried out to show the good performance of ore＂ schemes.