Search

We propose a generalized Eulerian-Lagrangian (GEL) discontinuous Galerkin (DG) method. The method is a generalization of the Eulerian-Lagrangian (EL) DG method for transport problems proposed in Cai et al. (2021) [5] , which tracks solution along approximations to characteristics in the DG framework, allowing extra large time stepping size with stability. The newly proposed GEL DG method in this paper is motivated for solving linear hyperbolic systems with variable coefficients, where the velocity field for adjoint problems of the test functions is frozen to constant. In this paper, in a simplified scalar setting, we propose the GEL DG methodology by freezing the velocity field of adjoint problems, and by formulating the semi-discrete scheme over the space-time region partitioned by linear lines approximating characteristics. The fully-discrete schemes are obtained by method-of-lines Runge-Kutta methods. We further design flux limiters for the schemes to satisfy the discrete geometric conservation law (DGCL) and maximum principle preserving (MPP) properties. Numerical results on 1D and 2D linear transport problems are presented to demonstrate great properties of the GEL DG method. These include the high order spatial and temporal accuracy, stability with extra large time stepping size, and satisfaction of DGCL and MPP properties. • We develop a generalized Eulerian Lagrangian (EL) discontinuous Galerkin (DG) method. • We establish connection of the generalized EL DG method with the EL DG method. • Numerical results on linear transport problems. [ABSTRACT FROM AUTHOR]

We propose a new Eulerian-Lagrangian (EL) discontinuous Galerkin (DG) method formulated by introducing a modified adjoint problem for the test function and by performing the integration of PDE over a space-time region partitioned by time-dependent linear functions approximating characteristics. The error incurred in characteristics approximation in the modified adjoint problem can then be taken into account by a new flux term, and can be integrated by method-of-line Runge-Kutta (RK) methods. The ELDG framework is designed as a generalization of the semi-Lagrangian (SL) DG method and classical Eulerian RK DG method for linear advection problems. It takes advantages of both formulations. In the EL DG framework, characteristics are approximated by a linear function in time, thus shapes of upstream cells are quadrilaterals in general two-dimensional problems. No quadratic-curved quadrilaterals are needed to design higher than second order schemes as in the SL DG scheme. On the other hand, the time step constraint from a classical Eulerian RK DG method is greatly mitigated, as it is evident from our theoretical and numerical investigations. Connection of the proposed EL DG method with the arbitrary Lagrangian-Eulerian (ALE) DG is observed. 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. • The novel design of the EL DG framework that accommodates both the SL DG method and classical Eulerian RK DG method. • Characteristics are approximated by a linear function in time, thus shapes of upstream cells are quadrilaterals in general two-dimensional problems. • The time step constraint from a classical Eulerian RK DG method is greatly mitigated. [ABSTRACT FROM AUTHOR]

We propose a new Eulerian-Lagrangian Runge-Kutta finite volume method for numerically solving convection and convection-diffusion equations. Eulerian-Lagrangian and semi-Lagrangian methods have grown in popularity mostly due to their ability to allow large time steps. Our proposed scheme is formulated by integrating the PDE on a space-time region partitioned by approximations of the characteristics determined from the Rankine-Hugoniot jump condition; and then rewriting the time-integral form into a time differential form to allow application of Runge-Kutta (RK) methods via the method-of-lines approach. The scheme can be viewed as a generalization of the standard Runge-Kutta finite volume (RK-FV) scheme for which the space-time region is partitioned by approximate characteristics with zero velocity. The high-order spatial reconstruction is achieved using the recently developed weighted essentially non-oscillatory schemes with adaptive order (WENO-AO); and the high-order temporal accuracy is achieved by explicit RK methods for convection equations and implicit-explicit (IMEX) RK methods for convection-diffusion equations. Our algorithm extends to higher dimensions via dimensional splitting. Numerical experiments demonstrate our algorithm's robustness, high-order accuracy, and ability to handle extra large time steps. • The novel design of the EL RK FV framework that accommodates both the SL FV method and classical Eulerian RK FV method. • Partition of space-time region that follows dynamics of characteristics. • The time step constraint from an Eulerian method is greatly mitigated. [ABSTRACT FROM AUTHOR]