An Eulerian-Lagrangian discontinuous Galerkin method for transport problems and its application to nonlinear dynamics
share
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 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.
READ FULL TEXT