An Eulerian-Lagrangian discontinuous Galerkin method for transport problems and its application to nonlinear dynamics

We propose a new Eulerian-Lagrangian (EL) discontinuous Galerkin (DG) method. The method is designed as a generalization of the semi-Lagrangian (SL) DG method for linear advection problems proposed in [J. Sci. Comput. 73: 514-542, 2017], which is formulated based on an adjoint problem and tracing upstream cells by tracking characteristics curves highly accurately. In the SLDG method, depending on the velocity field, upstream cells could be of arbitrary shape. Thus, a more sophisticated approximation to sides of the upstream cells is required to get high order approximation. For example, quadratic-curved (QC) quadrilaterals were proposed to approximate upstream cells for a third-order spatial accuracy in a swirling deformation example. In this paper, for linear advection problems, we propose a more general formulation, named the ELDG method. The scheme is formulated based on a {\em modified} adjoint problem for which the upstream cells are always quadrilaterals, which avoids the need to use QC quadrilaterals in the SLDG algorithm. The newly proposed ELDG method can be viewed as a new general framework, in which both the classical Eulerian Runge-Kutta DG formulation and the SL DG formulation can fit in. Numerical results on linear transport problems, as well as the nonlinear Vlasov and incompressible Euler dynamics using the exponential RK time integrators, are presented to demonstrate the effectiveness of the ELDG method.

[1]  Elena Celledoni,et al.  Commutator-free Lie group methods , 2003, Future Gener. Comput. Syst..

[2]  Tao Xiong,et al.  High Order Multi-dimensional Characteristics Tracing for the Incompressible Euler Equation and the Guiding-Center Vlasov Equation , 2018, J. Sci. Comput..

[3]  Jianxian Qiu,et al.  An h-Adaptive RKDG Method for the Two-Dimensional Incompressible Euler Equations and the Guiding Center Vlasov Model , 2017, J. Sci. Comput..

[4]  Eric Sonnendrücker,et al.  Conservative semi-Lagrangian schemes for Vlasov equations , 2010, J. Comput. Phys..

[5]  David C. Seal,et al.  A positivity-preserving high-order semi-Lagrangian discontinuous Galerkin scheme for the Vlasov-Poisson equations , 2010, J. Comput. Phys..

[6]  Chi-Wang Shu,et al.  TVB Runge-Kutta local projection discontinuous galerkin finite element method for conservation laws. II: General framework , 1989 .

[7]  Jing-Mei Qiu,et al.  A Conservative Semi-Lagrangian Discontinuous Galerkin Scheme on the Cubed Sphere , 2014 .

[8]  Todd Arbogast,et al.  An Eulerian–Lagrangian Weighted Essentially Nonoscillatory scheme for nonlinear conservation laws , 2017 .

[9]  Elena Celledoni,et al.  Semi-Lagrangian Runge-Kutta Exponential Integrators for Convection Dominated Problems , 2009, J. Sci. Comput..

[10]  Chi-Wang Shu,et al.  The Local Discontinuous Galerkin Method for Time-Dependent Convection-Diffusion Systems , 1998 .

[11]  P. Ciarlet,et al.  Mathematical elasticity, volume I: Three-dimensional elasticity , 1989 .

[12]  Julia S. Mullen,et al.  Filter-based stabilization of spectral element methods , 2001 .

[13]  Xiangxiong Zhang,et al.  On maximum-principle-satisfying high order schemes for scalar conservation laws , 2010, J. Comput. Phys..

[14]  Andrew M. Bradley,et al.  Conservative Multimoment Transport along Characteristics for Discontinuous Galerkin Methods , 2019, SIAM J. Sci. Comput..

[15]  Chuanju Xu,et al.  Stabilization Methods for Spectral Element Computations of Incompressible Flows , 2006, J. Sci. Comput..

[16]  Todd Arbogast,et al.  An Implicit Eulerian–Lagrangian WENO3 Scheme for Nonlinear Conservation Laws , 2018, J. Sci. Comput..

[17]  Richard E. Ewing,et al.  A Family of Eulerian-Lagrangian Localized Adjoint Methods for Multi-dimensional Advection-Reaction Equations , 1999 .

[18]  Jing-Mei Qiu,et al.  High Order Semi-Lagrangian Discontinuous Galerkin Method Coupled with Runge-Kutta Exponential Integrators for Nonlinear Vlasov Dynamics , 2021, J. Comput. Phys..

[19]  Lorenzo Pareschi,et al.  Modeling and Computational Methods for Kinetic Equations , 2012 .

[20]  Douglas N. Arnold,et al.  Unified Analysis of Discontinuous Galerkin Methods for Elliptic Problems , 2001, SIAM J. Numer. Anal..

[21]  Ilaria Perugia,et al.  An A Priori Error Analysis of the Local Discontinuous Galerkin Method for Elliptic Problems , 2000, SIAM J. Numer. Anal..

[22]  Dongmi Luo,et al.  A quasi-Lagrangian moving mesh discontinuous Galerkin method for hyperbolic conservation laws , 2018, J. Comput. Phys..

[23]  Jaime Peraire,et al.  Discontinuous Galerkin Solution of the Navier-Stokes Equations on Deformable Domains , 2007 .

[24]  T. F. Russell,et al.  An overview of research on Eulerian-Lagrangian localized adjoint methods (ELLAM). , 2002 .

[25]  Wei Guo,et al.  A High Order Conservative Semi-Lagrangian Discontinuous Galerkin Method for Two-Dimensional Transport Simulations , 2016, Journal of Scientific Computing.

[26]  Christian Klingenberg,et al.  Arbitrary Lagrangian-Eulerian discontinuous Galerkin method for conservation laws: Analysis and application in one dimension , 2016, Math. Comput..

[27]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[28]  Wei Guo,et al.  A high order semi-Lagrangian discontinuous Galerkin method for Vlasov-Poisson simulations without operator splitting , 2017, J. Comput. Phys..

[29]  Luming Wang,et al.  A high-order discontinuous Galerkin method with unstructured space–time meshes for two-dimensional compressible flows on domains with large deformations , 2015 .

[30]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[31]  Chi-Wang Shu,et al.  The Runge-Kutta local projection $P^1$-discontinuous-Galerkin finite element method for scalar conservation laws , 1988, ESAIM: Mathematical Modelling and Numerical Analysis.

[32]  Magdi Shoucri A two-level implicit scheme for the numerical solution of the linearized vorticity equation , 1981 .

[33]  Todd Arbogast,et al.  An Eulerian-Lagrangian WENO finite volume scheme for advection problems , 2012, J. Comput. Phys..

[34]  Yinhua Xia,et al.  Stability analysis and error estimates of arbitrary Lagrangian–Eulerian discontinuous Galerkin method coupled with Runge–Kutta time-marching for linear conservation laws , 2019, ESAIM: Mathematical Modelling and Numerical Analysis.

[35]  Jianxian Qiu,et al.  An h-Adaptive RKDG Method for the Vlasov–Poisson System , 2016, J. Sci. Comput..

[36]  Zhengfu Xu,et al.  High order maximum principle preserving semi-Lagrangian finite difference WENO schemes for the Vlasov equation , 2013, J. Comput. Phys..

[37]  ShuChi-Wang,et al.  Efficient implementation of essentially non-oscillatory shock-capturing schemes, II , 1989 .

[38]  Todd Arbogast,et al.  A semi-Lagrangian finite difference WENO scheme for scalar nonlinear conservation laws , 2016, J. Comput. Phys..

[39]  Chi-Wang Shu,et al.  A High-Order Discontinuous Galerkin Method for 2D Incompressible Flows , 2000 .

[40]  Chi-Wang Shu,et al.  Positivity preserving semi-Lagrangian discontinuous Galerkin formulation: Theoretical analysis and application to the Vlasov-Poisson system , 2011, J. Comput. Phys..

[41]  T. F. Russell,et al.  An Eulerian-Lagrangian localized adjoint method for the advection-diffusion equation , 1990 .