Editorial Type:
Article
Article Type:
Research Article

The Contour-Advective Semi-Lagrangian Algorithm for the Shallow Water Equations

David G. Dritschel Mathematics Institute, University of St. Andrews, Fife, United Kingdom

Search for other papers by David G. Dritschel in
Current site
Google Scholar
PubMed Close
,
Lorenzo M. Polvani Department of Applied Physics and Applied Mathematics, Columbia University, New York, New York

Search for other papers by Lorenzo M. Polvani in
Current site
Google Scholar
PubMed Close
, and
Ali R. Mohebalhojeh Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge, United Kingdom

Search for other papers by Ali R. Mohebalhojeh in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 1999
DOI:
https://doi.org/10.1175/1520-0493(1999)127<1551:TCASLA>2.0.CO;2
Page(s):
1551–1565
Restricted access

Abstract

A new method for integrating shallow water equations, the contour-advective semi-Lagrangian (CASL) algorithm, is presented. This is the first implementation of a contour method to a system of equations for which exact potential vorticity invertibility does not exist. The new CASL method fuses the recent contour-advection technique with the traditional pseudospectral (PS) method. The potential vorticity field, which typically develops steep gradients and evolves into thin filaments, is discretized by level sets separated by contours that are advected in a fully Lagrangian way. The height and divergence fields, which are intrinsically broader in scale, are treated in an Eulerian way: they are discretized on an fixed grid and time stepped with a PS scheme.

In fact, the CASL method is similar to the widely used semi-Lagrangian (SL) method in that material conservation of potential vorticity along particle trajectories is used to determine the potential vorticity at each time step from the previous one. The crucial difference is that, whereas in the CASL method the potential vorticity is merely advected, in the SL method the potential vorticity needs to be interpolated at each time step. This interpolation results in numerical diffusion in the SL method.

By directly comparing the CASL, SL, and PS methods, it is demonstrated that the implicit diffusion associated with potential vorticity interpolation in the SL method and the explicit diffusion required for numerical stability in the PS method seriously degrade the solution accuracy compared with the CASL method. Moreover, it is shown that the CASL method is much more efficient than the SL and PS methods since, for a given solution accuracy, a much coarser grid can be used and hence much faster computations can be performed.

Corresponding author address: Dr. Lorenzo M. Polvani, Department of Applied Physics and Applied Mathematics, Columbia University, 209 S W MUDD, Mail Code 4701, New York, NY 10027.

Email: polvani@columbia.edu

Abstract

A new method for integrating shallow water equations, the contour-advective semi-Lagrangian (CASL) algorithm, is presented. This is the first implementation of a contour method to a system of equations for which exact potential vorticity invertibility does not exist. The new CASL method fuses the recent contour-advection technique with the traditional pseudospectral (PS) method. The potential vorticity field, which typically develops steep gradients and evolves into thin filaments, is discretized by level sets separated by contours that are advected in a fully Lagrangian way. The height and divergence fields, which are intrinsically broader in scale, are treated in an Eulerian way: they are discretized on an fixed grid and time stepped with a PS scheme.

In fact, the CASL method is similar to the widely used semi-Lagrangian (SL) method in that material conservation of potential vorticity along particle trajectories is used to determine the potential vorticity at each time step from the previous one. The crucial difference is that, whereas in the CASL method the potential vorticity is merely advected, in the SL method the potential vorticity needs to be interpolated at each time step. This interpolation results in numerical diffusion in the SL method.

By directly comparing the CASL, SL, and PS methods, it is demonstrated that the implicit diffusion associated with potential vorticity interpolation in the SL method and the explicit diffusion required for numerical stability in the PS method seriously degrade the solution accuracy compared with the CASL method. Moreover, it is shown that the CASL method is much more efficient than the SL and PS methods since, for a given solution accuracy, a much coarser grid can be used and hence much faster computations can be performed.

Corresponding author address: Dr. Lorenzo M. Polvani, Department of Applied Physics and Applied Mathematics, Columbia University, 209 S W MUDD, Mail Code 4701, New York, NY 10027.

Email: polvani@columbia.edu

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Ambaum, M. H., and D. G. Dritschel, 1998: The inclusion of non-conservative forcing into a conservative, contour advection algorithm. Numerical Methods for Fluid Mechanics VI, M. J. Baines, Ed., ICFD, 99–110.

  • Asselin, R., 1972: Frequency filter for time integrations. Mon. Wea. Rev.,100, 487–490.

  • Bates, J. R., Y. Li, A. Brandt, S. F. McCormick, and J. Ruge, 1995:A global shallow-water numerical model based on the semi-Lagrangian advection of potential vorticity. Quart. J. Roy. Meteor. Soc.,121, 1981–2005.

  • Dritschel, D. G., 1989: Contour dynamics and contour surgery: Numerical algorithms for extended, high-resolution modelling of vortex dynamics in two-dimensional, inviscid, incompressible flows. Comp. Phys. Rep.,10, 77–146.

  • ——, and M. H. P. Ambaum, 1997: A contour-advective semi-Lagrangian algorithm for the simulation of fine-scale conservative fields. Quart. J. Roy. Meteor. Soc.,123, 1097–1130.

  • Edouard, S., B. Legras, F. Lefèvre, and R. Eymard, 1996: The effect of small-scale inhomogeneities on ozone depletion. Nature,384, 444–447.

  • Gravel, S., 1996: The semi-Lagrangian Method. Proc. Centre de Recherches Mathematiques, Montreal, PQ, Canada, 1–11.

  • Hoskins, B. J., M. E. McIntyre, and A. W. Robertson, 1985: On the use and significance of isentropic potential-vorticity maps. Quart. J. Roy. Meteor. Soc.,111, 877–946.

  • Jiménez, J., 1994: Hyperviscous vortices. J. Fluid Mech.,279, 169–176.

  • Macaskill, C., and B. M. Bewick, 1995: The use of hyperviscosity in simulations of geophysical fluid flows. Proc. 12th Australasian Fluid Mechanics Conf., Sydney, Australia, 415–418.

  • Mariotti, A., B. Legras, and D. G. Dritschel, 1994: Vortex stripping and the erosion of coherent structures in two-dimensional flows. Phys. Fluids,6, 3954–3962.

  • McIntyre, M. E., 1995: The stratospheric polar vortex and sub-vortex:Fluid dynamics and midlatitude ozone loss. Philos. Trans. Roy. Soc. London,352, 227–240.

  • Methven, J., and B. Hoskins, 1998: Spirals in potential vorticity. Part I: Measures of structure. J. Atmos. Sci.,55, 2053–2066.

  • Norton, W. A., 1988: Balance and potential vorticity inversion in atmospheric dynamics. Ph.D. thesis, University of Cambridge, 167 pp.

  • ——, 1994: Breaking Rossby waves in a model stratosphere diagnosed by a vortex-following coordinate system and a technique for advecting material contours. J. Atmos. Sci.,51, 654–673.

  • Plumb, R. A., and Coauthors, 1994: Intrusions into the lower stratospheric Arctic vortex during the winter of 1991–1992. J. Geophys. Res.,99, 1089–1106.

  • Ritchie, H., 1988: Application of the semi-Lagrangian method to a spectral model of the shallow water equations. Mon. Wea. Rev.,116, 1687–1598.

  • ——, C. Temperton, A. Simmons, M. Hortal, T. Davies, D. Dent, and M. Harmud, 1995: Implementation of the semi-Lagrangian method in a high resolution version of the ECMWF forecast model. Mon. Wea. Rev.,123, 489–514.

  • Robert, A. J., 1966: The integration of a low order spectral form of the primitive meteorological equations. J. Meteor. Soc. Japan,44, 237–245.

  • Staniforth, A., and J. Côté, 1991: Semi-Lagrangian integration schemes for atmospheric models—A review. Mon. Wea. Rev.,119, 2206–2223.

  • Tan, D. G. H., P. H. Haynes, A. R. MacKenzie, and J. A. Pyle, 1997:Effects of fluid-dynamical stirring and mixing on the deactivation of stratospheric chlorine. J. Geophys. Res.,103 (D1), 1585–1605.

  • Temperton, C., and A. Staniforth, 1987: An efficient two-time-level semi-Lagrangian semi-implicit integrating scheme. Quart. J. Roy. Meteor. Soc.,113, 1025–1040.

  • Waugh, D. W., and R. A. Plumb, 1994: Contour advection with surgery: A technique for investigating finescale structure in tracer transport. J. Atmos. Sci.,51, 530–540.

  • ——, and Coauthors, 1994: Transport of material out of the stratospheric Arctic vortex by Rossby wave breaking. J. Geophys. Res.,99, 1071–1088.

  • Williamson, D. L., J. B. Drake, J. J. Hack, R. Jakob, and P. N. Swarztrauber, 1992: A standard test set for numerical approximations to the shallow water equations in spherical geometry. J. Comput. Phys.,102, 211–224.

  • Wunsch, C., and D. Stammer, 1995: The global frequency-wavenumber spectrum of oceanic variability estimated from TOPEX/Poseidon altimetric measurements. J. Geophys. Res.,100 (C12), 24 895–24 910.

  • Yao, H. B., D. G. Dritschel, and N. J. Zabusky, 1995: High-gradient phenomena in 2D vortex interactions. Phys. Fluids,7, 539–548.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 190 47 5
PDF Downloads 172 53 8