Superconvergence Analysis of the Runge-Kutta Discontinuous Galerkin Method with Upwind-Biased Numerical Flux for Two-Dimensional Linear Hyperbolic Equation

In this paper,we shall establish the superconvergence properties of the Runge-Kutta dis-continuous Galerkin method for solving two-dimensional linear constant hyperbolic equa-tion,where the upwind-biased numerical flux is used.By suitably defining the correction function and deeply understanding the mechanisms when the spatial derivatives and the correction manipulations are carried out along the same or different directions,we obtain the superconvergence results on the node averages,the numerical fluxes,the cell averages,the solution and the spatial derivatives.The superconvergence properties in space are pre-served as the semi-discrete method,and time discretization solely produces an optimal order error in time.Some numerical experiments also are given.

Direct discontinuous Galerkin method for the generalized Burgers-Fisher equation

In this study, we use the direct discontinuous Galerkin method to solve the generalized Burgers-Fisher equation. The method is based on the direct weak formulation of the Burgers-Fisher equation. The two adjacent cells are jointed by a numerical flux that includes the convection numerical flux and the diffusion numerical flux. We solve the ordinary differential equations arising in the direct Galerkin method by using the strong stability preserving Runge^Kutta method. Numerical results are compared with the exact solution and the other results to show the accuracy and reliability 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.

Central Discontinuous Galerkin Method for the Navier-Stokes Equations

Central discontinuous Galerkin（CDG）method is used to solve the Navier-Stokes equations for viscous flow in this paper.The CDG method involves two pieces of approximate solutions defined on overlapping meshes.Taking advantages of the redundant representation of the solution on the overlapping meshes,the cell interface of one computational mesh is right inside the staggered mesh,hence approximate Riemann solvers are not needed at cell interfaces.Third order total variation diminishing（TVD）Runge-Kutta（RK）methods are applied in time discretization.Numerical examples for 1D and2 D viscous flow simulations are presented to validate the accuracy and robustness of the CDG method.

POSITIVITY-PRESERVING LOCAL DISCONTINUOUS GALERKIN METHOD FOR PATTERN FORMATION DYNAMICAL MODEL IN POLYMERIZING ACTIN FLOCKS

In this paper,we apply local discontinuous Galerkin(LDG)methods for pattern formation dynamical model in polymerizing actin focks.There are two main dificulties in designing effective numerical solvers.First of all,the density function is non-negative,and zero is an unstable equilibrium solution.Therefore,negative density values may yield blow-up solutions.To obtain positive numerical approximations,we apply the positivitypreserving(PP)techniques.Secondly,the model may contain stif source.The most commonly used time integration for the PP technique is the strong-stability-preserving Runge-Kutta method.However,for problems with stiff source,such time discretizations may require strictly limited time step sizes,leading to large computational cost.Moreover,the stiff source any trigger spurious filament polarization,leading to wrong numerical approximations on coarse meshes.In this paper,we combine the PP LDG methods with the semi-implicit Runge-Kutta methods.Numerical experiments demonstrate that the proposed method can yield accurate numerical approximations with relatively large time steps.

High-Order Accurate Runge-Kutta (Local) Discontinuous Galerkin Methods for One- and Two-Dimensional Fractional Diffusion Equations

As the generalization of the integer order partial differential equations(PDE),the fractional order PDEs are drawing more and more attention for their applications in fluid flow,finance and other areas.This paper presents high-order accurate Runge-Kutta local discontinuous Galerkin(DG)methods for one-and two-dimensional fractional diffusion equations containing derivatives of fractional order in space.The Caputo derivative is chosen as the representation of spatial derivative,because it may represent the fractional derivative by an integral operator.Some numerical examples show that the convergence orders of the proposed local Pk–DG methods are O(hk+1)both in one and two dimensions,where Pk denotes the space of the real-valued polynomials with degree at most k.

Simulations of Compressible Two-Medium Flow by Runge-Kutta Discontinuous Galerkin Methods with the Ghost Fluid Method

The original ghost fluid method (GFM) developed in [13] and the modifiedGFM (MGFM) in [26] have provided a simple and yet flexible way to treat twomediumflow problems. The original GFM and MGFM make the material interface"invisible" during computations and the calculations are carried out as for a singlemedium such that its extension to multi-dimensions becomes fairly straightforward.The Runge-Kutta discontinuous Galerkin (RKDG) method for solving hyperbolic conservationlaws is a high order accurate finite element method employing the usefulfeatures from high resolution finite volume schemes, such as the exact or approximateRiemann solvers, TVD Runge-Kutta time discretizations, and limiters. In this paper,we investigate using RKDG finite element methods for two-medium flow simulationsin one and two dimensions in which the moving material interfaces is treated via nonconservativemethods based on the original GFM and MGFM. Numerical results forboth gas-gas and gas-water flows are provided to show the characteristic behaviors ofthese combinations.

Runge-Kutta Discontinuous Galerkin Method Using WENO-Type Limiters:Three-Dimensional

This paper further considers weighted essentially non-oscillatory(WENO)and Hermite weighted essentially non-oscillatory(HWENO)finite volume methods as limiters for Runge-Kutta discontinuous Galerkin(RKDG)methods to solve problems involving nonlinear hyperbolic conservation laws.The application discussed here is the solution of 3-D problems on unstructured meshes.Our numerical tests again demonstrate this is a robust and high order limiting procedure,which simultaneously achieves high order accuracy and sharp non-oscillatory shock transitions.

Stability analysis and a priori error estimate of explicit Runge-Kutta discontinuous Galerkin methods for correlated random walk with density-dependent turning rates

In this paper,we analyze the explicit Runge-Kutta discontinuous Galerkin(RKDG)methods for the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology.The RKDG methods use a third order explicit total-variation-diminishing Runge-Kutta(TVDRK3)time discretization and upwinding numerical fluxes.By using the energy method,under a standard CourantFriedrichs-Lewy(CFL)condition,we obtain L2stability for general solutions and a priori error estimates when the solutions are smooth enough.The theoretical results are proved for piecewise polynomials with any degree k 1.Finally,since the solutions to this system are non-negative,we discuss a positivity-preserving limiter to preserve positivity without compromising accuracy.Numerical results are provided to demonstrate these RKDG methods.

An h-Adaptive Runge-Kutta Discontinuous Galerkin Method for Hamilton-Jacobi Equations

In[35,36],we presented an h-adaptive Runge-Kutta discontinuous Galerkin method using troubled-cell indicators for solving hyperbolic conservation laws.A tree data structure(binary tree in one dimension and quadtree in two dimensions)is used to aid storage and neighbor finding.Mesh adaptation is achieved by refining the troubled cells and coarsening the untroubled"children".Extensive numerical tests indicate that the proposed h-adaptive method is capable of saving the computational cost and enhancing the resolution near the discontinuities.In this paper,we apply this h-adaptive method to solve Hamilton-Jacobi equations,with an objective of enhancing the resolution near the discontinuities of the solution derivatives.One-and two-dimensional numerical examples are shown to illustrate the capability of the method.

A Well-Balanced Runge-Kutta Discontinuous Galerkin Method for Multilayer ShallowWater Equations with Non-Flat Bottom Topography

A well-balanced Runge-Kutta discontinuous Galerkin method is presented for the numerical solution of multilayer shallow water equations with mass exchange and non-flat bottom topography.The governing equations are reformulated as a non-linear system of conservation laws with differential source forces and reaction terms.Coupling between theflow layers is accounted for in the system using a set of ex-change relations.The considered well-balanced Runge-Kutta discontinuous Galerkin method is a locally conservativefinite element method whose approximate solutions are discontinuous across the inter-element boundaries.The well-balanced property is achieved using a special discretization of source terms that depends on the nature of hydrostatic solutions along with the Gauss-Lobatto-Legendre nodes for the quadra-ture used in the approximation of source terms.The method can also be viewed as a high-order version of upwindfinite volume solvers and it offers attractive features for the numerical solution of conservation laws for which standardfinite element methods fail.To deal with the source terms we also implement a high-order splitting operator for the time integration.The accuracy of the proposed Runge-Kutta discontinuous Galerkin method is examined for several examples of multilayer free-surfaceflows over bothflat and non-flat beds.The performance of the method is also demonstrated by comparing the results obtained using the proposed method to those obtained using the incompressible hydrostatic Navier-Stokes equations and a well-established kinetic method.The proposed method is also applied to solve a recirculationflow problem in the Strait of Gibraltar.

Runge-Kutta Discontinuous Galerkin Method with Front Tracking Method for Solving the Compressible Two-Medium Flow on Unstructured Meshes

In this paper,we extend using the Runge-Kutta discontinuous Galerkin method together with the front tracking method to simulate the compressible twomedium flow on unstructured meshes.A Riemann problem is constructed in the normal direction in the material interfacial region,with the goal of obtaining a compact,robust and efficient procedure to track the explicit sharp interface precisely.Extensive numerical tests including the gas-gas and gas-liquid flows are provided to show the proposed methodologies possess the capability of enhancing the resolutions nearby the discontinuities inside of the single medium flow and the interfacial vicinities of the two-medium flow in many occasions.

Enforcing Strong Stability of Explicit Runge-Kutta Methods with Superviscosity

A time discretization method is called strongly stable(or monotone),if the norm of its numerical solution is nonincreasing.Although this property is desirable in various of contexts,many explicit Runge-Kutta(RK)methods may fail to preserve it.In this paper,we enforce strong stability by modifying the method with superviscosity,which is a numerical technique commonly used in spectral methods.Our main focus is on strong stability under the inner-product norm for linear problems with possibly non-normal operators.We propose two approaches for stabilization:the modified method and the filtering method.The modified method is achieved by modifying the semi-negative operator with a high order superviscosity term;the filtering method is to post-process the solution by solving a diffusive or dispersive problem with small superviscosity.For linear problems,most explicit RK methods can be stabilized with either approach without accuracy degeneration.Furthermore,we prove a sharp bound(up to an equal sign)on diffusive superviscosity for ensuring strong stability.For nonlinear problems,a filtering method is investigated.Numerical examples with linear non-normal ordinary differential equation systems and for discontinuous Galerkin approximations of conservation laws are performed to validate our analysis and to test the performance.

ERROR ESTIMATE ON A FULLY DISCRETE LOCAL DISCONTINUOUS GALERKIN METHOD FOR LINEAR CONVECTION-DIFFUSION PROBLEM

In this paper we present the error estimate for the fully discrete local discontinuous Galerkin algorithm to solve the linear convection-diffusion equation with Dirichlet boundary condition in one dimension. The time is advanced by the third order explicit total variation diminishing Runge-Kutta method under the reasonable temporal-spatial condition as general. The optimal error estimate in both space and time is obtained by aid of the energy technique, if we set the numerical flux and the intermediate boundary condition properly.

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.

LOCAL ANALYSIS OF THE FULLY DISCRETE LOCAL DISCONTINUOUS GALERKIN METHOD FOR THE TIME-DEPENDENT SINGULARLY PERTURBED PROBLEM

In this paper we consider the fully discrete local discontinuous Galerkin method, where the third order explicit Runge-Kutta time marching is coupled. For the one-dimensional time-dependent singularly perturbed problem with a boundary layer, we shall prove that the resulted scheme is not only of good behavior at the local stability, but also has the double-optimal local error estimate. It is to say, the convergence rate is optimal in both space and time, and the width of the cut-off subdomain is also nearly optimal, if the boundary condition at each intermediate stage is given in a proper way. Numerical experiments are also given.

Nodal Discontinuous Galerkin Method for Time-Domain Lorentz Model Equations in Meta-Materials

In this paper,Nodal discontinuous Galerkin method is presented to approxi-mate Time-domain Lorentz model equations in meta-materials.The upwind flux is cho-sen in spatial discrete scheme.Low-storage five-stage fourth-order explicit Runge-Kutta method is employed in time discrete scheme.An error estimate of accuracy O(τ^(4)+h^(n))is proved under the L^(2)-norm,specially O(τ^(4)+h^(n+1))can be obtained.Numerical exper-iments for transverse electric(TE)case and transverse magnetic(TM)case are demon-strated to verify the stability and the efficiency of the method in low and higher wave frequency.

Numerical Diffusion Control of a Space-Time Discontinuous Galerkin Method

Variations on space-time Discontinuous Galerkin(STDG)discretization associated to Runge-Kutta schemes are developed.These new schemes while keeping the original scheme order can improve accuracy and stability.Numerical analysis is made on academic test cases and efficiency of these schemes are shown on propagating pressure waves.

Adaptive local discontinuous Galerkin methods with semi-implicit time discretizations for the Navier-Stokes equations

In this paper,we present a mesh adaptation algorithm for the unsteady compressible Navier-Stokes equations under the framework of local discontinuous Galerkin methods coupled with implicit-explicit Runge-Kutta or spectral deferred correction time discretization methods.In both of the two high order semi-implicit time integration methods,the convective flux is treated explicitly and the viscous and heat fluxes are treated implicitly.The remarkable benefits of such semi-implicit temporal discretizations are that they can not only overcome the stringent time step restriction compared with time explicit methods,but also avoid the construction of the large Jacobian matrix as is done for fully implicit methods,thus are relatively easy to implement.To save computing time as well as capture the flow structures of interest accurately,a local mesh refinement(h-adaptive)technique,in which we present detailed criteria for selecting candidate elements and complete strategies to refine and coarsen them,is also applied for the Navier-Stokes equations.Numerical experiments are provided to illustrate the high order accuracy,efficiency and capabilities of the semi-implicit schemes in combination with adaptive local discontinuous Galerkin methods for the Navier-Stokes equations.

小当量柱型装药水下近场爆炸固支单层方形钢板毁伤特性研究

为研究小当量柱形装药水下近场爆炸固支单层方形钢板毁伤特性,采用高精度龙格库塔间断迦辽金(Runge-Kutta discontinuous Galerkin,RKDG)方法求解柱形装药水下近场爆炸载荷,并与水下爆炸经验公式进行对比,吻合良好。将计算得到的爆炸载荷加载到LS-DYNA非线性有限元求解器中,得到固支单层方形钢板的毁伤特性,与试验结果进行对比,证明模型的精确性,表明龙格库塔间断迦辽金-有限元法(Runge-Kutta discontinuous Galerkin-finite element method,RKDG-FEM)耦合计算模型能够精确模拟柱形装药在近场爆炸条件下固支单层方形钢板的响应特性。随后研究了固支单层方形钢板在板厚4 mm、药量5~30 g和板厚3~8 mm、药量20 g的毁伤模式。研究发现,固支单层方形钢板的毁伤模式主要为塑性大变形、花瓣型破口和冲塞型圆形破口等。并拟合了考虑药量、板厚因素的破口尺寸估算公式,可为小当量柱形装药水下近场爆炸单层板的毁伤研究提供参考。