Insights on a sign-preserving numerical method for the advection-diffusion equation

Abstract

In this paper we explore theoretically and numerically the application of the advection transport algorithm introduced by Smolarkiewicz to the one dimensional unsteady advection diffusion equation. The scheme consists of a sequence of upwind iterations, where the initial iteration is the first order accurate upwind scheme, while the subsequent iterations are designed to compensate for the truncation error of preceding step. Two versions of the method are discussed. One, the classical version of the method, regards the second order terms of the truncation error and the other considers additionally the third order terms. Stability and convergence are discussed and the theoretical considerations are illustrated through numerical tests. The numerical tests will also indicate in which situations is advantageous to use the numerical methods presented.