A semi-Lagrangian discontinuous Galerkin method for scalar advection by incompressible flows

A new, conservative semi-Lagrangian formulation is proposed for the discretization of the scalar advection equation in flux form. The approach combines the accuracy and conservation properties of the Discontinuous Galerkin (DG) method with the computational efficiency and robustness of Semi-Lagrangian (SL) techniques. Unconditional stability in the von Neumann sense is proved for the proposed discretization in the one-dimensional case. A monotonization technique is then introduced, based on the Flux Corrected Transport approach. This yields a multi-dimensional monotonic scheme for the piecewise constant component of the computed solution that is characterized by a smaller amount of numerical diffusion than standard DG methods. The accuracy and stability of the method are further demonstrated by two-dimensional tracer advection tests in the case of incompressible flows. The comparison with results obtained by standard SL and DG methods highlights several advantages of the new technique. (c) 2005 Elsevier Inc. All rights reserved.