A COUPLING METHOD OF DIFFERENCE WITH HIGH ORDER ACCURACY AND BOUNDARY INTEGRAL EQUATION FOR EVOLUTIONARY EQUATION AND ITS ERROR ESTIMATES

In the present paper, a new numerical method for solving initial-boundary value problems of evolutionary equations is proposed and studied, combining difference method with high accuracy with boundary integral equation method. The numerical approximate schemes for both problems on a bounded or unbounded domain in R3 are proposed and their prior error estimates are obtained.

High Order Finite Difference WENO Methods for Shallow Water Equations on Curvilinear Meshes

A high order finite difference numerical scheme is developed for the shallow water equations on curvilinear meshes based on an alternative flux formulation of the weighted essentially non-oscillatory(WENO)scheme.The exact C-property is investigated,and comparison with the standard finite difference WENO scheme is made.Theoretical derivation and numerical results show that the proposed finite difference WENO scheme can maintain the exact C-property on both stationarily and dynamically generalized coordinate systems.The Harten-Lax-van Leer type flux is developed on general curvilinear meshes in two dimensions and verified on a number of benchmark problems,indicating smaller errors compared with the Lax-Friedrichs solver.In addition,we propose a positivity-preserving limiter on stationary meshes such that the scheme can preserve the non-negativity of the water height without loss of mass conservation.

Hermite WENO-based limiters for high order discontinuous Galerkin method on unstructured grids

A novel class of weighted essentially nonoscillatory （WENO） schemes based on Hermite polynomi- als, termed as HWENO schemes, is developed and applied as limiters for high order discontinuous Galerkin （DG） method on triangular grids. The developed HWENO methodology utilizes high-order derivative information to keep WENO re- construction stencils in the von Neumann neighborhood. A simple and efficient technique is also proposed to enhance the smoothness of the existing stencils, making higher-order scheme stable and simplifying the reconstruction process at the same time. The resulting HWENO-based limiters are as compact as the underlying DG schemes and therefore easy to implement. Numerical results for a wide range of flow conditions demonstrate that for DG schemes of up to fourth order of accuracy, the designed HWENO limiters can simul- taneously obtain uniform high order accuracy and sharp, es- sentially non-oscillatory shock transition.

High-order discontinuous Galerkin solver on hybrid anisotropic meshes for laminar and turbulent simulations

Efficient and robust solution strategies are developed for discontinuous Galerkin （DG） discretization of the Navier-Stokes （NS） and Reynolds-averaged NS （RANS） equations on structured/unstructured hybrid meshes. A novel line-implicit scheme is devised and implemented to reduce the memory gain and improve the computational eificiency for highly anisotropic meshes. A simple and effective technique to use the mod- ified Baldwin-Lomax （BL） model on the unstructured meshes for the DC methods is proposed. The compact Hermite weighted essentially non-oscillatory （HWENO） limiters are also investigated for the hybrid meshes to treat solution discontinuities. A variety of compressible viscous flows are performed to examine the capability of the present high- order DG solver. Numerical results indicate that the designed line-implicit algorithms exhibit weak dependence on the cell aspect-ratio as well as the discretization order. The accuracy and robustness of the proposed approaches are demonstrated by capturing com- plex flow structures and giving reliable predictions of benchmark turbulent problems.

A high order boundary scheme to simulate complex moving rigid body under impingement of shock wave

In the paper, we study a high order numerical boundary scheme for solving the complex moving boundary problem on a fixed Cartesian mesh, and numerically investigate the moving rigid body with the complex boundary under the impingement of an inviscid shock wave. Based on the high order inverse Lax-Wendroff(ILW) procedure developed in the previous work(TAN, S. and SHU, C. W. A high order moving boundary treatment for compressible inviscid flows. Journal of Computational Physics, 230(15),6023–6036(2011)), in which the authors only considered the translation of the rigid body,we consider both translation and rotation of the body in this paper. In particular, we reformulate the material derivative on the moving boundary with no-penetration condition, and the newly obtained formula plays a key role in the proposed algorithm. Several numerical examples, including cylinder, elliptic cylinder, and NACA0012 airfoil, are given to indicate the effectiveness and robustness of the present method.

High Order Semi-implicit Multistep Methods for Time-Dependent Partial Differential Equations

We consider the construction of semi-implicit linear multistep methods that can be applied to time-dependent PDEs where the separation of scales in additive form,typically used in implicit-explicit(IMEX)methods,is not possible.As shown in Boscarino et al.(J.Sci.Comput.68:975-1001,2016)for Runge-Kutta methods,these semi-implicit techniques give a great flexibility,and allow,in many cases,the construction of simple linearly implicit schemes with no need of iterative solvers.In this work,we develop a general setting for the construction of high order semi-implicit linear multistep methods and analyze their stability properties for a prototype lineal'advection-diffusion equation and in the setting of strong stability preserving(SSP)methods.Our findings are demonstrated on several examples,including nonlinear reaction-diffusion and convection-diffusion problems.

A Provable Positivity-Preserving Local Discontinuous Galerkin Method for the Viscous and Resistive MHD Equations

In this paper,we construct a high-order discontinuous Galerkin(DG)method which can preserve the positivity of the density and the pressure for the viscous and resistive magnetohydrodynamics(VRMHD).To control the divergence error in the magnetic field,both the local divergence-free basis and the Godunov source term would be employed for the multi-dimensional VRMHD.Rigorous theoretical analyses are presented for one-dimensional and multi-dimensional DG schemes,respectively,showing that the scheme can maintain the positivity-preserving(PP)property under some CFL conditions when combined with the strong-stability-preserving time discretization.Then,general frameworks are established to construct the PP limiter for arbitrary order of accuracy DG schemes.Numerical tests demonstrate the effectiveness of the proposed schemes.

A High Order Positivity-Preserving Discontinuous Galerkin Remapping Method Based on a Moving Mesh Solver for ALE Simulation of the Compressible Fluid Flow

The arbitrary Lagrangian-Eulerian(ALE)method is widely used in the field of compressible multi-material and multi-phase flow problems.In order to implement the indirect ALE approach for the simulation of compressible flow in the context of high order discontinuous Galerkin(DG)discretizations,we present a high order positivity-preserving DG remapping method based on a moving mesh solver in this paper.This remapping method is based on the ALE-DG method developed by Klingenberg et al.[17,18]to solve the trivial equation∂u/∂t=0 on a moving mesh,which is the old mesh before remapping at t=0 and is the new mesh after remapping at t=T.An appropriate selection of the final pseudo-time T can always satisfy the relatively mild smoothness requirement(Lipschitz continuity)on the mesh movement velocity,which guarantees the high order accuracy of the remapping procedure.We use a multi-resolution weighted essentially non-oscillatory(WENO)limiter which can keep the essentially non-oscillatory property near strong discontinuities while maintaining high order accuracy in smooth regions.We further employ an effective linear scaling limiter to preserve the positivity of the relevant physical variables without sacrificing conservation and the original high order accuracy.Numerical experiments are provided to illustrate the high order accuracy,essentially non-oscillatory performance and positivity-preserving of our remapping algorithm.In addition,the performance of the ALE simulation based on the DG framework with our remapping algorithm is examined in one-and two-dimensional Euler equations.

Local Discontinuous Galerkin Methods for High-Order Time-Dependent Partial Differential Equations

Discontinuous Galerkin (DG) methods are a class of finite element methodsusing discontinuous basis functions, which are usually chosen as piecewise polynomi-als. Since the basis functions can be discontinuous, these methods have the flexibilitywhich is not shared by typical finite element methods, such as the allowance of ar-bitrary triangulation with hanging nodes, less restriction in changing the polynomialdegrees in each element independent of that in the neighbors (p adaptivity), and localdata structure and the resulting high parallel efficiency. In this paper, we give a generalreview of the local DG (LDG) methods for solving high-order time-dependent partialdifferential equations (PDEs). The important ingredient of the design of LDG schemes,namely the adequate choice of numerical fluxes, is highlighted. Some of the applica-tions of the LDG methods for high-order time-dependent PDEs are also be discussed.

High-Order Accurate Entropy Stable Finite Difference Schemes for One- and Two-Dimensional Special Relativistic Hydrodynamics

This paper develops the high-order accurate entropy stable finite difference schemes for one-and two-dimensional special relativistic hydrodynamic equations.The schemes are built on the entropy conservative flux and the weighted essentially non-oscillatory(WENO)technique as well as explicit Runge-Kutta time discretization.The key is to technically construct the affordable entropy conservative flux of the semi-discrete second-order accurate entropy conservative schemes satisfying the semi-discrete entropy equality for the found convex entropy pair.As soon as the entropy conservative flux is derived,the dissipation term can be added to give the semidiscrete entropy stable schemes satisfying the semi-discrete entropy inequality with the given convex entropy function.The WENO reconstruction for the scaled entropy variables and the high-order explicit Runge-Kutta time discretization are implemented to obtain the fully-discrete high-order entropy stable schemes.Several numerical tests are conducted to validate the accuracy and the ability to capture discontinuities of our entropy stable schemes.

A New Approach of High OrderWell-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a Class of Hyperbolic Systems with Source Terms

Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms.In our earlier work[31–33],we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms.In this paper,we present a different approach to the same purpose:designing high order well-balanced finite volume weighted essentially non-oscillatory(WENO)schemes and RungeKutta discontinuous Galerkin(RKDG)finite element methods.We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly,if a small modification on either the initial condition or the flux is provided.The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method.The same idea can be applied to the finite volume WENO schemes.We will first describe the algorithms and prove the well balanced property for the shallow water equations,and then show that the result can be generalized to a class of other balance laws.We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions,the non-oscillatory property for general solutions with discontinuities,and the genuine high order accuracy in smooth regions.

Numerical Simulation of Time-Harmonic Waves in Inhomogeneous Media using Compact High Order Schemes

In many problems,one wishes to solve the Helmholtz equation with variable coefficients within the Laplacian-like term and use a high order accurate method(e.g.,fourth order accurate)to alleviate the points-per-wavelength constraint by reducing the dispersion errors.The variation of coefficients in the equation may be due to an inhomogeneous medium and/or non-Cartesian coordinates.This renders existing fourth order finite difference methods inapplicable.We develop a new compact scheme that is provably fourth order accurate even for these problems.We present numerical results that corroborate the fourth order convergence rate for several model problems.

High Order Finite Difference WENO Methods with Unequal-Sized Sub-Stencils for the Degasperis-Procesi Type Equations

In this paper,we develop twofinite difference weighted essentially non-oscillatory(WENO)schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi(DP)andµ-Degasperis-Procesi(µDP)equations,which contain nonlinear high order derivatives,and possibly peakon solutions or shock waves.By introducing auxiliary variable(s),we rewrite the DP equation as a hyperbolic-elliptic system,and theµDP equation as afirst order system.Then we choose a linearfinite difference scheme with suitable order of accuracy for the auxiliary variable(s),and twofinite difference WENO schemes with unequal-sized sub-stencils for the primal variable.One WENO scheme uses one large stencil and several smaller stencils,and the other WENO scheme is based on the multi-resolution framework which uses a se-ries of unequal-sized hierarchical central stencils.Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil,both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights.Another advantage is that thefinal reconstructed polynomial on the target cell is a polynomial of the same de-gree as the polynomial over the big stencil,while the classicalfinite difference WENO reconstruction can only be obtained for specific points inside the target interval.Nu-merical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.

Development and Comparison of Numerical Fluxes for LWDG Methods

The discontinuous Galerkin （DO） or local discontinuous Galerkin （LDG） method is a spatial discretization procedure for convection-diffusion equations, which employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes and limiters. The Lax- Wendroff time discretization procedure is an altemative method for time discretization to the popular total variation diminishing （TVD） Runge-Kutta time discretizations. In this paper, we develop fluxes for the method of DG with Lax-Wendroff time discretization procedure （LWDG） based on different numerical fluxes for finite volume or finite difference schemes, including the first-order monotone fluxes such as the Lax-Friedfichs flux, Godunov flux, the Engquist-Osher flux etc. and the second-order TVD fluxes. We systematically investigate the performance of the LWDG methods based on these different numerical fluxes for convection terms with the objective of obtaining better performance by choosing suitable numerical fluxes. The detailed numerical study is mainly performed for the one-dimensional system case, addressing the issues of CPU cost, accuracy, non-oscillatory property, and resolution of discontinuities. Numerical tests are also performed for two dimensional systems.

A STOPPING CRITERION FOR HIGHER-ORDER SWEEPING SCHEMES FOR STATIC HAMILTON-JACOBI EQUATIONS

We propose an effective stopping criterion for higher-order fast sweeping schemes for static Hamilton-Jacobi equations based on ratios of three consecutive iterations. To design the new stopping criterion we analyze the convergence of the first-order Lax-Friedrichs sweeping scheme by using the theory of nonlinear iteration. In addition, we propose a fifth-order Weighted PowerENO sweeping scheme for static Hamilton-Jacobi equations with convex Hamiltonians and present numerical examples that validate the effectiveness of the new stopping criterion.

An Efficient High Order Well-Balanced Finite Difference WENO Scheme for the Blood Flow Model

The blood flow model admits the steady state,in which the flux gradient is non-zero and is exactly balanced by the source term.In this paper,we present a high order well-balanced finite difference weighted essentially non-oscillatory(WENO)scheme,which exactly preserves the steady state.In order to maintain the wellbalanced property,we propose to reformulate the equation and apply a novel source term approximation.Extensive numerical experiments are carried out to verify the performances of the current scheme such as the maintenance of well-balanced property,the ability to capture the perturbations of such steady state and the genuine high order accuracy for smooth solutions.

High Order Accurate Direct Arbitrary-Lagrangian-Eulerian ADER-MOOD Finite Volume Schemes for Non-Conservative Hyperbolic Systems with Stiff Source Terms

In this paper we present a 2D/3D high order accurate finite volume scheme in the context of direct Arbitrary-Lagrangian-Eulerian algorithms for general hyperbolic systems of partial differential equations with non-conservative products and stiff source terms.This scheme is constructed with a single stencil polynomial reconstruction operator,a one-step space-time ADER integration which is suitably designed for dealing even with stiff sources,a nodal solver with relaxation to determine the mesh motion,a path-conservative integration technique for the treatment of non-conservative products and an a posteriori stabilization procedure derived from the so-called Multidimensional Optimal Order Detection(MOOD)paradigm.In this work we consider the seven equation Baer-Nunziato model of compressible multi-phase flows as a representative model involving non-conservative products as well as relaxation source terms which are allowed to become stiff.The new scheme is validated against a set of test cases on 2D/3D unstructured moving meshes on parallel machines and the high order of accuracy achieved by the method is demonstrated by performing a numerical convergence study.Classical Riemann problems and explosion problems with exact solutions are simulated in 2D and 3D.The overall numerical code is also profiled to provide an estimate of the computational cost required by each component of the whole algorithm.

On the Monotonicity of Q^(3) Spectral Element Method for Laplacian

The monotonicity of discrete Laplacian, i.e., inverse positivity of stiffness matrix, implies discrete maximum principle, which is in general not true for high order accurate schemes on unstructured meshes. On the other hand,it is possible to construct high order accurate monotone schemes on structured meshes. All previously known high order accurate inverse positive schemes are or can be regarded as fourth order accurate finite difference schemes, which is either an M-matrix or a product of two M-matrices. For the Q3spectral element method for the two-dimensional Laplacian, we prove its stiffness matrix is a product of four M-matrices thus it is unconditionally monotone. Such a scheme can be regarded as a fifth order accurate finite difference scheme. Numerical tests suggest that the unconditional monotonicity of Q^(k) spectral element methods will be lost for k ≥ 9 in two dimensions, and for k ≥ 4 in three dimensions. In other words, for obtaining a high order monotone scheme, only Q^(2) and Q^(3) spectral element methods can be unconditionally monotone in three dimensions.

A Third Order Conservative Lagrangian Type Scheme on Curvilinear Meshes for the Compressible Euler Equations

Based on the high order essentially non-oscillatory(ENO)Lagrangian type scheme on quadrilateral meshes presented in our earlier work[3],in this paper we develop a third order conservative Lagrangian type scheme on curvilinear meshes for solving the Euler equations of compressible gas dynamics.The main purpose of this work is to demonstrate our claim in[3]that the accuracy degeneracy phenomenon observed for the high order Lagrangian type scheme is due to the error from the quadrilateral mesh with straight-line edges,which restricts the accuracy of the resulting scheme to at most second order.The accuracy test given in this paper shows that the third order Lagrangian type scheme can actually obtain uniformly third order accuracy even on distorted meshes by using curvilinear meshes.Numerical examples are also presented to verify the performance of the third order scheme on curvilinear meshes in terms of resolution for discontinuities and non-oscillatory properties.

Genuinely Multidimensional Physical-Constraints-Preserving FiniteVolume Schemes for the Special Relativistic Hydrodynamics

This paper develops the genuinely multidimensional HLL Riemann solver for the two-dimensional special relativistic hydrodynamic equations on Cartesian meshes and studies its physical-constraint-preserving(PCP)property.Based on the resulting HLL solver,the first-and high-order accurate PCP finite volume schemes are proposed.In the high-order scheme,the WENO reconstruction,the third-order accurate strong-stability-preserving time discretizations and the PCP flux limiter are used.Several numerical results are given to demonstrate the accuracy,performance and resolution of the shock waves and the genuinely multi-dimensional wave structures etc.of our PCP finite volume schemes.