Bound-Preserving Discontinuous Galerkin Methods with Modified Patankar Time Integrations for Chemical Reacting Flows

In this paper,we develop bound-preserving discontinuous Galerkin(DG)methods for chemical reactive flows.There are several difficulties in constructing suitable numerical schemes.First of all,the density and internal energy are positive,and the mass fraction of each species is between 0 and 1.Second,due to the rapid reaction rate,the system may contain stiff sources,and the strong-stability-preserving explicit Runge-Kutta method may result in limited time-step sizes.To obtain physically relevant numerical approximations,we apply the bound-preserving technique to the DG methods.Though traditional positivity-preserving techniques can successfully yield positive density,internal energy,and mass fractions,they may not enforce the upper bound 1 of the mass fractions.To solve this problem,we need to(i)make sure the numerical fluxes in the equations of the mass fractions are consistent with that in the equation of the density;(ii)choose conservative time integrations,such that the summation of the mass fractions is preserved.With the above two conditions,the positive mass fractions have summation 1,and then,they are all between 0 and 1.For time discretization,we apply the modified Runge-Kutta/multi-step Patankar methods,which are explicit for the flux while implicit for the source.Such methods can handle stiff sources with relatively large time steps,preserve the positivity of the target variables,and keep the summation of the mass fractions to be 1.Finally,it is not straightforward to combine the bound-preserving DG methods and the Patankar time integrations.The positivity-preserving technique for DG methods requires positive numerical approximations at the cell interfaces,while Patankar methods can keep the positivity of the pre-selected point values of the target variables.To match the degree of freedom,we use polynomials on rectangular meshes for problems in two space dimensions.To evolve in time,we first read the polynomials at the Gaussian points.Then,suitable slope limiters can be applied to enforce the positivity of the solutions at those points,which can be preserved by the Patankar methods,leading to positive updated numerical cell averages.In addition,we use another slope limiter to get positive solutions used for the bound-preserving technique for the flux.Numerical examples are given to demonstrate the good performance of the proposed schemes.

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.

Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

In this paper,numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin(DG)methods(Dong et al.in J Sci Comput 66:321–345,2016;Dong and Wang in J Comput Appl Math 380:1–11,2020)for a one-dimensional stationary Schrödinger equation.Previous work showed that penalty parameters were required to be positive in error analysis,but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes.In this work,by performing extensive numerical experiments,we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods,and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

A Local Macroscopic Conservative(LoMaC)Low Rank Tensor Method with the Discontinuous Galerkin Method for the Vlasov Dynamics

In this paper,we propose a novel Local Macroscopic Conservative(LoMaC)low rank tensor method with discontinuous Galerkin(DG)discretization for the physical and phase spaces for simulating the Vlasov-Poisson(VP)system.The LoMaC property refers to the exact local conservation of macroscopic mass,momentum,and energy at the discrete level.The recently developed LoMaC low rank tensor algorithm(arXiv:2207.00518)simultaneously evolves the macroscopic conservation laws of mass,momentum,and energy using the kinetic flux vector splitting;then the LoMaC property is realized by projecting the low rank kinetic solution onto a subspace that shares the same macroscopic observables.This paper is a generalization of our previous work,but with DG discretization to take advantage of its compactness and flexibility in handling boundary conditions and its superior accuracy in the long term.The algorithm is developed in a similar fashion as that for a finite difference scheme,by observing that the DG method can be viewed equivalently in a nodal fashion.With the nodal DG method,assuming a tensorized computational grid,one will be able to(i)derive differentiation matrices for different nodal points based on a DG upwind discretization of transport terms,and(ii)define a weighted inner product space based on the nodal DG grid points.The algorithm can be extended to the high dimensional problems by hierarchical Tucker(HT)decomposition of solution tensors and a corresponding conservative projection algorithm.In a similar spirit,the algorithm can be extended to DG methods on nodal points of an unstructured mesh,or to other types of discretization,e.g.,the spectral method in velocity direction.Extensive numerical results are performed to showcase the efficacy of the method.

Stability and Time-Step Constraints of Implicit-Explicit Runge-Kutta Methods for the Linearized Korteweg-de Vries Equation

This paper provides a study on the stability and time-step constraints of solving the linearized Korteweg-de Vries(KdV)equation,using implicit-explicit(IMEX)Runge-Kutta(RK)time integration methods combined with either finite difference(FD)or local discontinuous Galerkin(DG)spatial discretization.We analyze the stability of the fully discrete scheme,on a uniform mesh with periodic boundary conditions,using the Fourier method.For the linearized KdV equation,the IMEX schemes are stable under the standard Courant-Friedrichs-Lewy(CFL)conditionτ≤λh.Here,λis the CFL number,τis the time-step size,and h is the spatial mesh size.We study several IMEX schemes and characterize their CFL number as a function ofθ=d/h^(2)with d being the dispersion coefficient,which leads to several interesting observations.We also investigate the asymptotic behaviors of the CFL number for sufficiently refined meshes and derive the necessary conditions for the asymptotic stability of the IMEX-RK methods.Some numerical experiments are provided in the paper to illustrate the performance of IMEX methods under different time-step constraints.

THE UNIFORM CONVERGENCE OF A DG METHOD FOR A SINGULARLY PERTURBED VOLTERRA INTEGRO-DIFFERENTIAL EQUATION

The purpose of this work is to implement a discontinuous Galerkin(DG)method with a one-sided flux for a singularly perturbed Volterra integro-differential equation(VIDE)with a smooth kernel.First,the regularity property and a decomposition of the exact solution of the singularly perturbed VIDE with the initial condition are provided.Then the existence and uniqueness of the DG solution are proven.Then some appropriate projection-type interpolation operators and their corresponding approximation properties are established.Based on the decomposition of the exact solution and the approximation properties of the projection type interpolants,the DG method achieves the uniform convergence in the L2 norm with respect to the singular perturbation parameter e when the space of polynomials with degree p is used.A numerical experiment validates the theoretical results.Furthermore,an ultra-convergence order 2p+1 at the nodes for the one-sided flux,uniform with respect to the singular perturbation parameter e,is observed numerically.

SPACE-TIME DISCONTINUOUS GALERKIN METHOD FOR MAXWELL EQUATIONS IN DISPERSIVE MEDIA

In this paper, a unified model for time-dependent Maxwell equations in dispersive media is considered. The space-time DG method developed in [29] is applied to solve the un-derlying problem. Unconditional L2-stability and error estimate of order O？τr＋1＋hk＋1/2？ are obtained when polynomials of degree at most r and k are used for the temporal dis-cretization and spatial discretization respectively. 2-D and 3-D numerical examples are given to validate the theoretical results. Moreover, numerical results show an ultra-convergence of order 2r＋1 in temporal variable t.

An h-Adaptivity DG Method on Locally Curved Tetrahedral Mesh for Solving Compressible Flows

For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions where the physical variables vary violently(for example,near the shock waves or in the boundary layers)and larger elements are expected for the regions where the solution is smooth.h-adaptive mesh has been widely used for complex flows.However,there are two difficulties when employing h-adaptivity for high-order discontinuous Galerkin(DG)methods.First,locally curved elements are required to precisely match the solid boundary,which significantly increases the difficulty to conduct the"refining"and"coarsening"operations since the curved information has to be maintained.Second,h-adaptivity could break the partition balancing,which would significantly affect the efficiency of parallel computing.In this paper,a robust and automatic h-adaptive method is developed for high-order DG methods on locally curved tetrahedral mesh,for which the curved geometries are maintained during the h-adaptivity.Furthermore,the reallocating and rebalancing of the computational loads on parallel clusters are conducted to maintain the parallel efficiency.Numerical results indicate that the introduced h-adaptive method is able to generate more reasonable mesh according to the structure of flow-fields.

A Uniformly Robust Staggered DG Method for the Unsteady Darcy-Forchheimer-Brinkman Problem

In this paper,we propose and analyze a uniformly robust staggered DG method for the unsteady Darcy-Forchheimer-Brinkman problem.Our formulation is based on velocity gradient-velocity-pressure and the resulting scheme can be flexibly applied to fairly general polygonal meshes.We relax the tangential continuity for velocity,which is the key ingredi-ent in achieving the uniform robustness.We present well-posedness and error analysis for both the semi-discrete scheme and the fully discrete scheme,and the theories indicate that the error estimates for velocity are independent of pressure.Several numerical experiments are presented to confirm the theoretical findings.

Low Dissipation Simulation for Vortex Flowfield of Rotor in Hover Based upon Discontinuous Galerkin Method

The discontinuous Galerkin(DG) method is established and innovatively conducted on accurately simulating the evolution of blade-tip vortex and the aerodynamic characteristics of helicopter rotor. Firstly,the Reynolds-Averaged Navier-Stokes(RANS)equations in rotating reference frame are employed,and the embedded grid system is developed with the finite volume method(FVM)and the DG method conducted on the blade grid and background grid respectively. Besides,the Harten-Lax-Van Leer contact(HLLC)scheme with high-resolution and low-dissipation is employed for spatial discretization,and the explicit third-order Runge-Kutta scheme is used to accomplish the temporal discretization. Secondly,the aerodynamic characteristics and the evolution of blade-tip vortex for Caradonna-Tung rotor are simulated by the established CFD method,and the numerical results are in good agreement with experimental data,which well validates the accuracy of the DG method and shows the advantages of DG method on capturing the detailed blade-tip vortex compared with the FVM method. Finally,the evolution of tip vortex at different blade tip Mach numbers and collective pitches is discussed.

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.

High-Order Discontinuous Galerkin Solution of Compressible Flows with a Hybrid Lattice Boltzmann Flux

A discontinuous Galerkin(DG)-based lattice Boltzmann method is employed to solve the Euler and Navier-Stokes equations.Instead of adopting the widely used local Lax-Friedrichs flux and Roe Flux etc.,a hybrid lattice Boltzmann flux solver(LBFS)is employed to evaluate the inviscid flux across the cell interfaces.The main advantage of the hybrid LBFS is its flexibility for capturing both strong shocks and thin boundary layers through introducing a function which varies from zero to one to control the artificial viscosity.Numerical results indicate that the hybrid lattice Boltzmann flux solver behaves very well combining with the high-order DG method when simulating both inviscid and viscous flows.

MODELING DAM-BREAK FLOOD OVER NATURAL RIVERS USING DISCONTINUOUS GALERKIN METHOD

A well-balanced numerical model is presented for two-dimensional, depth-averaged, shallow water flows based on the Discontinuous Galerkin （DG） method. The model is applied to simulate dam-break flood in natural rivers with wet/dry bed and complex topography. To eliminate numerical imbalance, the pressure force and bed slope terms are combined in the shallow water flow equations. For partially wet/dry elements, a treatment of the source term that preserves the well-balanced property is presented. A treatment for modeling flow over initially dry bed is presented. Numerical results show that the time step used is related to the dry bed criterion. The intercell numerical flux in the DG method is computed by the Harten-Lax-van Contact （HLLC） approximate Riemann solver. A two-dimensional slope limiting procedure is employed to prevent spurious oscillation. The robustness and accuracy of the model are demonstrated through several test cases, including dam-break flow in a channel with three bumps, laboratory dam-break tests over a triangular bump and an L-shape bend, dam-break flood in the Paute River, and the Malpasset dam-break case. Numerical results show that the model is robust and accurate to simulate dam-break flood over natural rivers with complex geometry and wet/dry beds.

A DISCONTINUOUS GALERKIN METHOD FOR THE FOURTH-ORDER CURL PROBLEM

In this paper, we present a discontinuous Galerkin （DG） method based on the N@d@lec finite element space for solving a fourth-order curl equation arising from a magnetohy- drodynamics model on a 3-dimensional bounded Lipschitz polyhedron. We show that the method has an optimal error estimate for a model problem involving a fourth-order curl operator. Furthermore, some numerical results in 2 dimensions are presented to verify the theoretical results.

A Well-Balanced Positivity-Preserving Quasi-Lagrange Moving Mesh DG Method for the Shallow Water Equations

A high-order, well-balanced, positivity-preserving quasi-Lagrange movingmesh DG method is presented for the shallow water equations with non-flat bottomtopography. The well-balance property is crucial to the ability of a scheme to simulate perturbation waves over the lake-at-rest steady state such as waves on a lake ortsunami waves in the deep ocean. The method combines a quasi-Lagrange movingmesh DG method, a hydrostatic reconstruction technique, and a change of unknownvariables. The strategies in the use of slope limiting, positivity-preservation limiting,and change of variables to ensure the well-balance and positivity-preserving properties are discussed. Compared to rezoning-type methods, the current method treatsmesh movement continuously in time and has the advantages that it does not need tointerpolate flow variables from the old mesh to the new one and places no constraintfor the choice of a update scheme for the bottom topography on the new mesh. A selection of one- and two-dimensional examples are presented to demonstrate the wellbalance property, positivity preservation, and high-order accuracy of the method andits ability to adapt the mesh according to features in the flow and bottom topography.

Time stepping in discontinuous Galerkin method

The time discretization in the Discontinuous Galerkin （DG） scheme has been traditionally based on the Total Variation Diminishing （TVD） second-order Runge-Kutta （RK2） scheme. Computational efficiency and accuracy with the Euler Forward （EF） and the TVD second-order RK2 time stepping schemes in the DG method are investigated in this work. Numerical tests are condu- cted with the scalar Burgers equation, 1-D and 2-D shallow water flow equations. The maximum Courant number or time step size required for stability for the EF scheme and RK2 scheme with different slope limiters are compared. Numerical results show that the slope limiters affect the stability requirement in the DG method. The RK2 scheme is generally more diffusive than the EF scheme, and the RK2 scheme allows larger time step sizes. The EF scheme is found to be more efficient and accurate than the RK2 scheme in the DG method in computation.

DISCONTINUOUS GALERKIN METHODS AND THEIR ADAPTIVITY FOR THE TEMPERED FRACTION AL(CONVECTION)DIFFUSION EQUATIONS

This paper focuses on the adaptive discontinuous Galerkin(DG)methods for the tempered fractional(convection)diffusion equations.The DG schemes with interior penalty for the diffusion term and numerical flux for the convection term are used to solve the equations,and the detailed stability and convergence analyses are provided.Based on the derived posteriori error estimates,the local error indicator is designed.The theoretical results and the effectiveness of the adaptive DG methods are,respectively,verified and displayed by the extensive numerical experiments.The strategy of designing adaptive schemes presented in this paper works for the general PDEs with fractional operators.

Implicit DG Method for Time Domain Maxwell’s Equations Involving Metamaterials

An implicit discontinuous Galerkin method is introduced to solve the timedomain Maxwell’s equations in metamaterials.The Maxwell’s equations in metamaterials are represented by integral-differential equations.Our scheme is based on discontinuous Galerkin method in spatial domain and Crank-Nicolson method in temporal domain.The fully discrete numerical scheme is proved to be unconditionally stable.When polynomial of degree at most p is used for spatial approximation,our scheme is verified to converge at a rate of O(τ^(2)+h^(p)+1/2).Numerical results in both 2D and 3D are provided to validate our theoretical prediction.

Discontinuous Galerkin Finite Element Method for a Nonlinear Boundary Value Problem

In this paper, we investigate the a priori and a posteriori error estimates for the discontinuous Galerkin finite element approximation to a regularization version of the variational inequality of the second kind. We show the optimal error estimates in the DG-norm （stronger than the H1 norm） and the L2 norm, respectively. Furthermore, some residual-based a posteriori error estimators are established which provide global upper bounds and local lower bounds on the discretization error. These a posteriori analysis results can be applied to develop the adaptive DG methods.

Error Estimates for the Iterative Discontinuous Galerkin Method to the Nonlinear Poisson-Boltzmann Equation

This paper is devoted to the error estimate for the iterative discontinuous Galerkin(IDG)method introduced in[P.Yin,Y.Huang and H.Liu.Commun.Comput.Phys.16:491-515,2014]to the nonlinear Poisson-Boltzmann equation.The total error includes both the iteration error and the discretization error of the direct DG method to linear elliptic equations.For the DDG method,the energy error is obtained by a constructive approach through an explicit global projection satisfying interface conditions dictated by the choice of numerical fluxes.The L^(2) error of order O(h^(m+1))for polynomials of degree m is further recovered.The bounding constant is also shown to be independent of the iteration times.Numerical tests are given to validate the established convergence theory.