Design Criteria for Finite-Difference Models for Eddy Diffusion With Winds That Guarantee Stability, Mass Conservation, and Nonnegative Masses

I. M. BASSETT School of Physics, Sydney University, Sydney, New South Wales, Australia

Search for other papers by I. M. BASSETT in
Current site
Google Scholar
PubMed
Close
,
R. G. L. HEWITT School of Physics, Sydney University, Sydney, New South Wales, Australia

Search for other papers by R. G. L. HEWITT in
Current site
Google Scholar
PubMed
Close
, and
BRIAN MARTIN School of Physics, Sydney University, Sydney, New South Wales, Australia

Search for other papers by BRIAN MARTIN in
Current site
Google Scholar
PubMed
Close
Restricted access

Abstract

A technique for following the circulation of a tracer in a turbulent fluid is developed from the integral form of the mass continuity equation. Numerical methods based on this technique are shown to be stable, to ensure that the total tracer mass is conserved, and that the mass in any region is always nonnegative. As an illustration of the utility of the technique, a numerical method is developed for a two-dimensional model of the stratosphere.

Abstract

A technique for following the circulation of a tracer in a turbulent fluid is developed from the integral form of the mass continuity equation. Numerical methods based on this technique are shown to be stable, to ensure that the total tracer mass is conserved, and that the mass in any region is always nonnegative. As an illustration of the utility of the technique, a numerical method is developed for a two-dimensional model of the stratosphere.

Save