High-Order Decoupled and Bound Preserving Local Discontinuous Galerkin Methods for a Class of Chemotaxis Models

In this paper,we explore bound preserving and high-order accurate local discontinuous Galerkin(LDG)schemes to solve a class of chemotaxis models,including the classical Keller-Segel(KS)model and two other density-dependent problems.We use the convex splitting method,the variant energy quadratization method,and the scalar auxiliary variable method coupled with the LDG method to construct first-order temporal accurate schemes based on the gradient flow structure of the models.These semi-implicit schemes are decoupled,energy stable,and can be extended to high accuracy schemes using the semi-implicit spectral deferred correction method.Many bound preserving DG discretizations are only worked on explicit time integration methods and are difficult to get high-order accuracy.To overcome these difficulties,we use the Lagrange multipliers to enforce the implicit or semi-implicit LDG schemes to satisfy the bound constraints at each time step.This bound preserving limiter results in the Karush-Kuhn-Tucker condition,which can be solved by an efficient active set semi-smooth Newton method.Various numerical experiments illustrate the high-order accuracy and the effect of bound preserving.

A Sub-element Adaptive Shock Capturing Approach for Discontinuous Galerkin Methods

In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.

A multithreaded parallel upwind sweep algorithm for the S_(N) transport equations discretized with discontinuous finite elements

The complex structure and strong heterogeneity of advanced nuclear reactor systems pose challenges for high-fidelity neutron-shielding calculations. Unstructured meshes exhibit strong geometric adaptability and can overcome the deficiencies of conventionally structured meshes in complex geometry modeling. A multithreaded parallel upwind sweep algorithm for S_(N) transport was proposed to achieve a more accurate geometric description and improve the computational efficiency. The spatial variables were discretized using the standard discontinuous Galerkin finite-element method. The angular flux transmission between neighboring meshes was handled using an upwind scheme. In addition, a combination of a mesh transport sweep and angular iterations was realized using a multithreaded parallel technique. The algorithm was implemented in the 2D/3D S_(N) transport code ThorSNIPE, and numerical evaluations were conducted using three typical benchmark problems:IAEA, Kobayashi-3i, and VENUS-3. These numerical results indicate that the multithreaded parallel upwind sweep algorithm can achieve high computational efficiency. ThorSNIPE, with a multithreaded parallel upwind sweep algorithm, has good reliability, stability, and high efficiency, making it suitable for complex shielding calculations.

Meshless Local Discontinuous Petrov-Galerkin Method with Application to Blasting Problems

A meshless local discontinuous Petrov-Galerkin (MLDPG) method based on the local symmetric weak form (LSWF) is presented with the application to blasting problems. The derivation is similar to that of mesh-based Runge-Kutta Discontinuous Galerkin (RKDG) method. The solutions are reproduced in a set of overlapped spherical sub-domains, and the test functions are employed from a partition of unity of the local basis functions. There is no need of any traditional non-overlapping mesh either for local approximation purpose or for Galerkin integration purpose in the presented method. The resulting MLDPG method is a meshless, stable, high-order accurate and highly parallelizable scheme which inherits both the advantages of RKDG and meshless method (MM), and it can handle the problems with extremely complicated physics and geometries easily. Three numerical examples of the one-dimensional Sod shock-tube problem, the blast-wave problem and the Woodward-Colella interacting shock wave problem are given. All the numerical results are in good agreement with the closed solutions. The higher-order MLDPG schemes can reproduce more accurate solution than the lower-order schemes.

A Local Discontinuous Galerkin Method for Two-Dimensional Time Fractional Diffusion Equations

For two-dimensional(2D)time fractional diffusion equations,we construct a numerical method based on a local discontinuous Galerkin(LDG)method in space and a finite differ-ence scheme in time.We investigate the numerical stability and convergence of the method for both rectangular and triangular meshes and show that the method is unconditionally stable.Numerical results indicate the effectiveness and accuracy of the method and con-firm the analysis.

A Posteriori Error Estimates for Finite Element Methods for Systems of Nonlinear,Dispersive Equations

The present study regards the numerical approximation of solutions of systems of Korteweg-de Vries type,coupled through their nonlinear terms.In our previous work[9],we constructed conservative and dissipative finite element methods for these systems and presented a priori error estimates for the semidiscrete schemes.In this sequel,we present a posteriori error estimates for the semidiscrete and fully discrete approximations introduced in[9].The key tool employed to effect our analysis is the dispersive reconstruction devel-oped by Karakashian and Makridakis[20]for related discontinuous Galerkin methods.We conclude by providing a set of numerical experiments designed to validate the a posteriori theory and explore the effectivity of the resulting error indicators.

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.

Discontinuous-Galerkin-Based Analysis of Traffic Flow Model Connected with Multi-Agent Traffic Model

As the number of automobiles continues to increase year after year,the associated problem of traffic congestion has become a serious societal issue.Initiatives to mitigate this problem have considered methods for optimizing traffic volumes in wide-area road networks,and traffic-flow simulation has become a focus of interest as a technique for advance characterization of such strategies.Classes of models commonly used for traffic-flow simulations include microscopic models based on discrete vehicle representations,macroscopic models that describe entire traffic-flow systems in terms of average vehicle densities and velocities,and mesoscopic models and hybrid(or multiscale)models incorporating both microscopic and macroscopic features.Because traffic-flow simulations are designed to model traffic systems under a variety of conditions,their underlyingmodelsmust be capable of rapidly capturing the consequences of minor variations in operating environments.In other words,the computation speed of macroscopic models and the precise representation of microscopic models are needed simultaneously.Thus,in this study we propose a multiscale model that combines a microscopic model—for detailed analysis of subregions containing traffic congestion bottlenecks or other localized phenomena of interest-with a macroscopic model enabling simulation of wide target areas at a modest computational cost.In addition,to ensure analytical stability with robustness in the presence of discontinuities,we discretize our macroscopic model using a discontinuous Galerkin finite element method(DGFEM),while to conjoin microscopic and macroscopic models,we use a generating/absorbing sponge layer,a technique widely used for numerical analysis of long-wavelength phenomena in shallow water,to enable traffic-flow simulations with stable input and output regions.

A combined mixed finite element method and local discontinuous Galerkin method for miscible displacement problem in porous media

A combined method consisting of the mixed finite element method for flow and the local discontinuous Galerkin method for transport is introduced for the one-dimensional coupled system of incompressible miscible displacement problem. Optimal error estimates in L∞(0,T;L2) for concentration c,in L2(0,T;L2)for cxand L∞(0,T;L2) for velocity u are derived. The main technical difficulties in the analysis include the treatment of the inter-element jump terms which arise from the discontinuous nature of the numerical method,the nonlinearity,and the coupling of the models. Numerical experiments are performed to verify the theoretical results. Finally,we apply this method to the one-dimensional compressible miscible displacement problem and give the numerical experiments to confirm the efficiency of the scheme.

FINITE ELEMENT AND DISCONTINUOUS GALERKIN METHOD FOR STOCHASTIC HELMHOLTZ EQUATION IN TWO-AND THREE-DIMENSIONS

In this paper, we consider the finite element method and discontinuous Galerkin method for the stochastic Helmholtz equation in R^d （d = 2, 3）. Convergence analysis and error estimates are presented for the numerical solutions. The effects of the noises on the accuracy of the approximations are illustrated. Numerical experiments are carried out to verify our theoretical results.

A Discontinuous Galerkin Finite Element Method without Interior Penalty Terms

A conforming discontinuous Galerkinfinite element method was introduced by Ye and Zhang,on simplicial meshes and on polytopal meshes,which has theflexibility of using discontinuous approximation and an ultra simple formulation.The main goal of this paper is to improve the above discontinuous Galerkinfinite element method so that it can handle nonhomogeneous Dirichlet boundary conditions effectively.In addition,the method has been generalized in terms of approximation of the weak gradient.Error estimates of optimal order are established for the correspond-ing discontinuousfinite element approximation in both a discrete H1 norm and the L2 norm.Numerical results are presented to confirm the theory.

Time discontinuous Galerkin space-time finite element method for nonlinear Sobolev equations

This article presents a complete discretization of a nonlinear Sobolev equation using space-time discontinuous Galerkin method that is discontinuous in time and continuous in space. The scheme is formulated by introducing the equivalent integral equation of the primal equation. The proposed scheme does not explicitly include the jump terms in time, which represent the discontinuity characteristics of approximate solution. And then the complexity of the theoretical analysis is reduced. The existence and uniqueness of the approximate solution and the stability of the scheme are proved. The optimalorder error estimates in L2 （H1） and L2 （L2） norms are derived. These estimates are valid under weak restrictions on the space-time mesh, namely, without the condition kn ≥ ch2, which is necessary in traditional space-time discontinuous Galerkin methods. Numerical experiments are presented to verify the theoretical results.

A Posteriori Error Estimates of a Combined Mixed Finite Element and Discontinuous Galerkin Method for a Kind of Compressible Miscible Displacement Problems

Akind of compressiblemiscible displacement problemswhich includemolecular diffusion and dispersion in porous media are investigated.The mixed finite element method is applied to the flow equation,and the transport one is solved by the symmetric interior penalty discontinuous Galerkin method.Based on a duality argument,employing projection estimates and approximation properties,a posteriori residual-type hp error estimates for the coupled system are presented,which is often used for guiding adaptivity.Comparing with the error analysis carried out by Yang(Int.J.Numer.Meth.Fluids,65(7)(2011),pp.781-797),the current work is more complicated and challenging.

A DISCONTINUOUS GALERKIN METHOD COMBINED WITH MIXED FINITE ELEMENT FOR SEAWATER INTRUSION PROBLEM

Seawater intrusion problem is considered in this paper.Its mathematical model is anonlinear coupled system of partial differential equations with initial boundary problem.It consistsof the water head equation and the salt concentration equation.A combined method is developedto approximate the water head equation by mixed finite element method and concentration equationby discontinuous Galerkin method.The scheme is continuous in time and optimal order estimates inH^1-norm and L^2-norm are derived for the errors.

Optimal Superconvergence of Energy Conserving Local Discontinuous Galerkin Methods for Wave Equations

This paper is concerned with numerical solutions of the LDG method for 1D wave equations.Superconvergence and energy conserving properties have been studied.We first study the superconvergence phenomenon for linear problems when alternating fluxes are used.We prove that,under some proper initial discretization,the numerical trace of the LDG approximation at nodes,as well as the cell average,converge with an order 2k+1.In addition,we establish k+2-th order and k+1-th order superconvergence rates for the function value error and the derivative error at Radau points,respectively.As a byproduct,we prove that the LDG solution is superconvergent with an order k+2 towards the Radau projection of the exact solution.Numerical experiments demonstrate that in most cases,our error estimates are optimal,i.e.,the error bounds are sharp.In the second part,we propose a fully discrete numerical scheme that conserves the discrete energy.Due to the energy conserving property,after long time integration,our method still stays accurate when applied to nonlinear Klein-Gordon and Sine-Gordon equations.

L1/Local Discontinuous Galerkin Method for the Time-Fractional Stokes Equation

In this paper,L1/local discontinuous Galerkin method seeking the numerical solution to the time-fractional Stokes equation is displayed,where the timefractional derivative is in the sense of Caputo with derivative orderα∈(0,1).Although the time-fractional derivative is used,its solution may be smooth since such examples can be easily constructed.In this case,we use the uniform L1 scheme to approach the temporal derivative and use the local discontinuous Galerkin(LDG)method to approximate the spatial derivative.If the solution has a certain weak regularity at the initial time,we use the non-uniform L1 scheme to discretize the time derivative and still apply LDG method to discretizing the spatial derivative.The numerical stability and error analysis for both situations are studied.Numerical experiments are also presented which support the theoretical analysis.

起伏地表弹性波传播的间断Galerkin有限元数值模拟方法

间断Galerkin有限元法(DG-FEM)作为一种有效的高阶有限元法受到了国内外学者的广泛关注.本文基于任意高阶间断Galerkin有限元法对弹性波方程进行空间离散,并将离散后所得的非齐次线性常微分方程系统齐次化,最后结合针对齐次问题的强稳定性保持龙格库塔(SSP Runge-Kutta)算法,将DG-FEM推广至时间任意高阶精度.另外,借鉴近最佳匹配层(NPML)的思想,基于复频移(CFS)拉伸坐标变换推导了一种新的PML吸收边界条件(简称为CFS-NPML),该CFS-NPML能够与DG-FEM算法很好地结合,形成有效的起伏地表地震波传播数值模拟技术.数值试验结果表明,DG-FEM具有高阶精度,可以适应任意复杂起伏地表和复杂构造情况下的弹性波传播数值模拟.同时,CFS-NPML对包括面波等震相的人为边界反射都具有良好的吸收效果.

任意单元间断Galerkin有限元计算方法研究

基于龙格库塔间断Galerkin(RKDG)有限元法的构造思想,通过局部坐标变换,发展了非正交单元DG有限元计算方法;借鉴非结构网格有限体积隐式计算方法,发展了适应于DG有限元方法的隐式计算方法;借鉴一维和二维(三角形单元)DG有限元限制器构造方法,提出了非正交三棱柱单元限制器方法;利用上述方法数值模拟了球头和双椭球的高超声速粘性绕流,得到了清晰的流场结构,并得到了较好的压力和热流分布,表明该方法在复杂高超声速流动的数值模拟方面具有广阔的应用前景。

Navier-Stokes方程间断Galerkin有限元方法研究

通过引入全局提升算子和局部提升算子,发展了求解Navier-Stokes方程的间断Galerkin(discontinuousGalerkin,DG)有限元方法的一般框架,并在此框架下给出了几种典型黏性离散格式的具体表达形式.对局部提升算子的求解给出了详细的计算步骤.同时还给出了一种简单有效的计算方法来对物面边界进行高阶近似.为了能够对NS方程进行精度测试,采用对原始系统添加源项的方法构造精确解.二维Euler和NS系统的精度测试表明该方法达到了DG方法的理论精度.二维圆柱无黏绕流的计算结果表明关于物面边界的高阶近似方法能够保持DG方法原有的精度.卡门涡街数值模拟则进一步验证了该方法的正确性并且显示出DG方法较高的计算精度和分辨率.