On stable and explicit numerical methods for the advection-diffusion equation

In this paper two stable and explicit numerical methods to integrate the one-dimensional (1D) advection-diffusion equation are presented. These schemes are stable by design and follow the main general concept behind the semi-Lagrangian method by constructing a virtual grid where the explicit method becomes stable. It is shown that the new schemes compare well with analytic solutions and are often more accurate than implicit schemes. In particular, the diffusion-only case is explored in some detail. The error produced by the stable and explicit method is a function of the ratio between the standard deviation @s"0 of the initial Gaussian state and the characteristic virtual grid distance @DS. Larger values of this ratio lead to very accurate results when compared to implicit methods, while lower values lead to less accuracy. It is shown that the @s"0/@DS ratio is also significant in the advection-diffusion problem: it determines the maximum error generated by new methods, obtained with a certain combination of the advection and diffusion values. In addition, the error becomes smaller when the problem becomes more advective or more diffusive.

[1]  João Teixeira.  Stable schemes for partial differential equations: the one-dimensional reaction-diffusion equation , 2004, Math. Comput. Simul..

[2]  A. Simmons,et al.  Implementation of the Semi-Lagrangian Method in a High-Resolution Version of the ECMWF Forecast Model , 1995 .

[3]  Guriĭ Ivanovich Marchuk,et al.  Numerical Methods in Weather Prediction , 1974 .

[4]  A. Robert.  A Semi-Lagrangian and Semi-Implicit Numerical Integration Scheme for the Primitive Meteorological Equations , 1982 .

[5]  G. J. Haltiner.  Numerical Prediction and Dynamic Meteorology , 1980 .

[6]  André Robert,et al.  A stable numerical integration scheme for the primitive meteorological equations , 1981 .

[7]  M. van Loon.  Tests on semi-Lagrangian transport and interpolation , 1992 .

[8]  A. Staniforth,et al.  Semi-Lagrangian integration schemes for atmospheric models - A review , 1991 .

[9]  C. Fletcher.  Computational techniques for fluid dynamics , 1992 .

[10]  C. Fletcher.  Computational techniques for fluid dynamics. Volume 1 - Fundamental and general techniques. Volume 2 - Specific techniques for different flow categories , 1988 .

[11]  Lance M. Leslie,et al.  An Efficient Interpolation Procedure for High-Order Three-Dimensional Semi-Lagrangian Models , 1991 .

[12]  David L. Williamson,et al.  A comparison of semi-lagrangian and Eulerian polar climate simulations , 1998 .

[13]  P. Smolarkiewicz,et al.  A class of semi-Lagrangian approximations for fluids. , 1992 .

[14]  Claude Girard,et al.  Stable Schemes for Nonlinear Vertical Diffusion in Atmospheric Circulation Models , 1990 .

[15]  Joao Teixeira.  Stable Schemes for Partial Differential Equations , 1999 .

[16]  João Paulo Teixeira,et al.  An eddy‐diffusivity/mass‐flux parametrization for dry and shallow cumulus convection , 2004 .

[17]  J. Doyle,et al.  A New Mixing-Length Formulation for the Parameterization of Dry Convection: Implementation and Evaluation in a Mesoscale Model , 2004 .

[18]  A. P. Siebesma.  An Advection-Diffusion scheme for the convective boundary layer: description and 1d-results , 2000 .