Stability analysis for Eulerian and semi-Lagrangian finite-element formulation of the advection-diffusion equation
Abstract
This paper analyzes the stability of the finite-element approximation to the linearized two-dimensional advection-diffusion equation. Bilinear basis functions on rectangular elements are considered. This is one of the two best schemes as was shown by Neta and Williams [1]. Time is discretized with the theta algorithms that yield the explicit (theta = 0), semi-implicit (theta = 1/2), and implicit (theta = 1) methods. This paper extends the results of Neta and Williams [1] for the advection equation. Giraldo and Neta [2] have numerically compared the Eulerian and semi-Lagrangian finite-element approximation for the advection-diffusion equation. This paper analyzes the finite element schemes used there. The stability analysis shows that the semi-Lagrangian method is unconditionally stable for all values of a while the Eulerian method is only unconditionally stable for 1/2 < theta < 1. This analysis also shows that the best methods are the semi-implicit ones (theta = 1/2). In essence this paper analytically compares a semi-implicit Eulerian method with a semi-implicit semi-Lagrangian method. It is concluded that (for small or no diffusion) the semi-implicit semi-Lagrangian method exhibits better amplitude, dispersion and group velocity errors than the semi-implicit Eulerian method thereby achieving better results. In the case the diffusion coefficient is large, the semi-Lagrangian loses its competitiveness. Published by Elsevier Science Ltd.
Description
The article of record as published may be located at http://dx.doi.org/10.1016/S0898-1221(99)00185-6
Rights
This publication is a work of the U.S. Government as defined in Title 17, United States Code, Section 101. Copyright protection is not available for this work in the United States.Collections
Related items
Showing items related by title, author, creator and subject.
-
A spectral element semi-Lagrangian (SESL) method for the spherical shallow water equations
Giraldo, F.X.; Perot, J. B.; Fischer, P. F. (2003);A spectral element semi-Lagrangian (SESL) method for the shallow water equations on the sphere is presented. The sphere is discretized using a hexahedral grid although any grid imaginable can be used as long as it is ... -
Hybrid Eulerian-Lagrangian semi-implicit time-integrators
Giraldo, F.X. (2006);Hybrid Eulerian-Lagrangian semi-implicit time-integrators (HELSI) are presented which use the standard semi-implicit formulation as their starting point. The advantage of such an approach is that existing models which ... -
Space-time transformation in flux-form semi-Lagrangian schemes
Fan, C.W.; Chu, Peter C. (2010);With a finite volume approach, a flux-form semi-Lagrangian (TFSL) scheme with space-time transformation was developed to provide stable and accurate algorithm in solving the advection-diffusion equation. Different from the ...