A high order conservative semi-Lagrangian discontinuous Galerkin method for two-dimensional transport simulations

TL;DR

This paper introduces a high order conservative semi-Lagrangian discontinuous Galerkin method for 2D transport equations, achieving accuracy, stability, and mass conservation without splitting errors, and verified through benchmark tests.

Contribution

The paper presents a novel high order, non-splitting semi-Lagrangian DG method with quadratic-curved cells and positivity preservation for multi-dimensional transport.

Findings

Achieves formal third order accuracy.

Maintains mass conservation and stability.

Demonstrates efficiency on benchmark problems.

Abstract

In this paper, we develop a class of high order conservative semi-Lagrangian (SL) discontinuous Galerkin (DG) methods for solving multi-dimensional linear transport equations. The methods rely on a characteristic Galerkin weak formulation, leading to $L^2$ stable discretizations for linear problems. Unlike many existing SL methods, the high order accuracy and mass conservation of the proposed methods are realized in a non-splitting manner. Thus, the detrimental splitting error, which is known to significantly contaminate long term transport simulations, will be not incurred. One key ingredient in the scheme formulation, borrowed from CSLAM [Lauritzen, Nair, and Ullrich, 2010], is the use of Green's theorem which allows us to convert volume integrals into a set of line integrals. The resulting line integrals are much easier to approximate with high order accuracy, hence facilitating the implementation. Another novel ingredient is the construction of quadratic curves in approximating sides of upstream cell, leading to quadratic-curved quadrilateral upstream cells. Formal third order accuracy is obtained by such a construction. The desired positivity-preserving property is further attained by incorporating a high order bound-preserving filter. To assess the performance of the proposed methods, we test and compare the numerical schemes with a variety of configurations for solving several benchmark transport problems with both smooth and nonsmooth solutions. The efficiency and efficacy are numerically verified.