AFullyMultidimensionalPositiveDefiniteAdvection.pdfVIP

JOURNAL OF COMPUTATIONAL PYSICS 54, 325-362 (1984) A Fully Multidimensional Positive Definite Advection Transport Algorithm with Small Implicit Diffusion PIOTR K. SMOLARKIEWICZ* National Centerfor Atmospheric Research,* Boulder, Colorado 80307 Received July 8, 1983 The idea of the simple positive definite advection scheme presented previously in Monthly Weather Review (1 I I (1983), 479) is improved for an optional multidimensional case and is presented in a generalized format. The accuracy of the scheme is discussed and a review of existing options is presented and illustrated through numerical tests. 1. INTR~~LJCTI~N In numerical modeling of physical phenomena it is often necessaryto solve the advective transport equation for positive definite scalar functions. Numerical schemes of second- or higher-order accuracy can produce negative values in the solution due to the dispersive ripples. Lower-order schemes, such as the donor cell or Lax-Friedrichs, or higher-order schemeswith zeroth-order diffusion added produce no ripples but suffer from excessiveimplicit diffusion. In the last ten years a possible resolution of this dilemma has been developed in the form of hybrid schemes, in which the advective fluxes are given as a weighted average of a first-order positive definite scheme’sfluxes and a higher-order scheme’sfluxes. The difference in deter- mination of the weights in the calculation of the average advective fluxes h