Editorial Type:
Article
Article Type:
Research Article

A Three-Dimensional, Adaptive, Godunov-Type Model for Global Atmospheric Flows

M. E. Hubbard DAMTP, University of Cambridge, Cambridge, United Kingdom

Search for other papers by M. E. Hubbard in
Current site
Google Scholar
PubMed Close
and
N. Nikiforakis DAMTP, University of Cambridge, Cambridge, United Kingdom

Search for other papers by N. Nikiforakis in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Aug 2003
DOI:
https://doi.org/10.1175//2568.1
Page(s):
1848–1864
Restricted access

Abstract

In this paper a Godunov-type methodology is applied to three-dimensional global atmospheric modeling. Numerical issues are addressed regarding the formulation of the tracer advection problem, the application of dimensional splitting, and the implementation of a Godunov-type scheme, based on the WAF approach, on spherical geometries. Particular attention is paid to addressing the problems that arise because of the convergence of the grid lines toward the Poles. A three-dimensional model is then built on the sphere that is based on a uniform longitude–latitude–height grid. This provides the framework within which an adaptive mesh refinement (AMR) algorithm is applied, to enhance the efficiency and accuracy with which results are obtained. These methods are not commonly used in the area of atmospheric modeling, but AMR in particular is commonly used with great success in other areas of computational fluid dynamics. The model is initially validated using a series of idealized case studies that have exact solutions, but is then developed into an offline model of tracer advection, forced by data from meteorological analyses, suitable to study the evolution of trace chemical species in the atmosphere.

Current affiliation: School of Computing, University of Leeds, Leeds, United Kingdom

Corresponding author address: Dr. M. E. Hubbard, School of Computing, University of Leeds, Leeds LS2 9JT, United Kingdom. Email: meh@comp.leeds.ac.uk

Abstract

In this paper a Godunov-type methodology is applied to three-dimensional global atmospheric modeling. Numerical issues are addressed regarding the formulation of the tracer advection problem, the application of dimensional splitting, and the implementation of a Godunov-type scheme, based on the WAF approach, on spherical geometries. Particular attention is paid to addressing the problems that arise because of the convergence of the grid lines toward the Poles. A three-dimensional model is then built on the sphere that is based on a uniform longitude–latitude–height grid. This provides the framework within which an adaptive mesh refinement (AMR) algorithm is applied, to enhance the efficiency and accuracy with which results are obtained. These methods are not commonly used in the area of atmospheric modeling, but AMR in particular is commonly used with great success in other areas of computational fluid dynamics. The model is initially validated using a series of idealized case studies that have exact solutions, but is then developed into an offline model of tracer advection, forced by data from meteorological analyses, suitable to study the evolution of trace chemical species in the atmosphere.

Current affiliation: School of Computing, University of Leeds, Leeds, United Kingdom

Corresponding author address: Dr. M. E. Hubbard, School of Computing, University of Leeds, Leeds LS2 9JT, United Kingdom. Email: meh@comp.leeds.ac.uk

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Behrens, J., 1996: An adaptive semi-Lagrangian advection scheme and its parallelization. Mon. Wea. Rev., 124 , 2386–2395.

  • Behrens, J., 1998: Atmospheric and ocean modeling with an adaptive finite element solver for the shallow-water equations. Appl. Numer. Math., 26 , 217–226.

    • Search Google Scholar
    • Export Citation

  • Bell, J., M. Berger, J. Saltzman, and M. Welcome, 1994: Three-dimensional adaptive mesh refinement for hyperbolic conservation laws. SIAM J. Sci. Comput., 15 , 127–138.

    • Search Google Scholar
    • Export Citation

  • Berger, M. J., cited 2002: Adaptive mesh refinement software for hyperbolic conservation laws. [Available online at http://cs.nyu.edu/cs/faculty/berger/amrsoftware.html.].

    • Search Google Scholar
    • Export Citation

  • Berger, M. J., and J. Oliger, 1984: Adaptive mesh refinement for hyperbolic partial differential equations. J. Comput. Phys., 53 , 482–512.

    • Search Google Scholar
    • Export Citation

  • Berger, M. J., and P. Colella, 1989: Local adaptive mesh refinement for shock hydrodynamics. J. Comput. Phys., 82 , 67–84.

  • Chipperfield, M. P., M. L. Santee, L. Froidevaux, G. L. Manney, W. G. Read, J. W. Waters, A. E. Roche, and J. M. Russell, 1996: Isentropic 3D chemical transport model and comparison with UARS data in southern polar vortex, September 1992. J. Geophys. Res., 101 (D13) 18861–18881.

    • Search Google Scholar
    • Export Citation

  • Colella, P., 1990: Multidimensional upwind methods for hyperbolic conservation laws. J. Comput. Phys., 87 , 171–200.

  • Doswell, C. A., 1984: A kinematic analysis of frontogenesis associated with a nondivergent vortex. J. Atmos. Sci., 41 , 1242–1248.

  • Easter, R. C., 1993: Two modified versions of Bott's positive-definite numerical advection scheme. Mon. Wea. Rev., 121 , 297–304.

  • Giraldo, F. X., 2000: The Lagrange–Galerkin method for the two-dimensional shallow water equations on adaptive grids. Int. J. Numer. Methods Fluids, 33 , 789–832.

    • Search Google Scholar
    • Export Citation

  • Leonard, B. P., 1979: A stable and accurate convection modeling procedure based on quadratic upstream interpolation. Comput. Methods Appl. Mech. Eng., 19 , 59–98.

    • Search Google Scholar
    • Export Citation

  • Leonard, B. P., and H. S. Niknafs, 1991: Sharp monotonic resolution of discontinuities without clipping of narrow extrema. Comput. Fluids, 19 , 141–154.

    • Search Google Scholar
    • Export Citation

  • Leonard, B. P., M. K. MacVean, and A. P. Lock, 1995: The flux integral method for multidimensional convection and diffusion. Appl. Math. Model, 19 , 333–342.

    • Search Google Scholar
    • Export Citation

  • Leonard, B. P., A. P. Lock, and M. K. MacVean, 1996: Conservative explicit unrestricted-time-step multidimensional constancy-preserving advection schemes. Mon. Wea. Rev., 124 , 2588–2606.

    • Search Google Scholar
    • Export Citation

  • Li, Y., and J. S. Chang, 1996: A mass-conservative, positive-definite, and efficient Eulerian advection scheme in spherical geometry and on a nonuniform grid system. J. Appl. Meteor., 35 , 1897–1913.

    • Search Google Scholar
    • Export Citation

  • Lin, S-J., and R. B. Rood, 1996: Multidimensional flux-form semi-Lagrangian transport schemes. Mon. Wea. Rev., 124 , 2046–2070.

  • Nair, R., J. Côté, and A. Staniforth, 1999: Cascade interpolation for semi-Lagrangian advection over the sphere. Quart. J. Roy. Meteor. Soc., 125 , 1445–1468.

    • Search Google Scholar
    • Export Citation

  • Petshek, A. G., and L. D. Libersky, 1975: Stability, accuracy and improvement of Crowley advection schemes. Mon. Wea. Rev., 103 , 1104–1109.

    • Search Google Scholar
    • Export Citation

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

    • Search Google Scholar
    • Export Citation

  • Prather, M. J., 1986: Numerical advection by conservation of second-order moments. J. Geophys. Res., 91 , 6671–6681.

  • Quirk, J. J., 1991: An adaptive grid algorithm for computational shock hydrodynamics. Ph.D. thesis, College of Aeronautics, Cranfield Institute of Technology, Cranfield, Bedfordshire, United Kingdom, 172 pp.

    • Search Google Scholar
    • Export Citation

  • Quirk, J. J., 1996: A parallel adaptive grid algorithm for computational shock hydrodynamics. Appl. Numer. Math., 20 , 427–453.

  • Quirk, J. J., and S. Karni, 1996: On the dynamics of shock-bubble interaction. J. Fluid Mech., 318 , 129–163.

  • Rasch, P. J., 1994: Conservative shape-preserving two-dimensional transport on a spherical reduced grid. Mon. Wea. Rev., 122 , 1337–1350.

    • Search Google Scholar
    • Export Citation

  • Rendleman, C. A., V. E. Beckner, M. Lijewski, W. Y. Crutchfield, and J. B. Bell, 2000: Parallelization of structured, hierarchical adaptive mesh refinement algorithms. Comput. Visualization Sci., 3 , 147–157.

    • Search Google Scholar
    • Export Citation

  • Skamarock, W. C., and J. B. Klemp, 1993: Adaptive grid refinement for two-dimensional and three-dimensional nonhydrostatic atmospheric flow. Mon. Wea. Rev., 121 , 788–804.

    • Search Google Scholar
    • Export Citation

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

    • Search Google Scholar
    • Export Citation

  • Strang, G., 1968: On the construction and comparison of difference schemes. SIAM J. Numer. Anal., 5 , 506–517.

  • Sweby, P. K., 1984: High resolution schemes using flux limiters for hyperbolic conservation laws. SIAM J. Numer. Anal., 21 , 995–1011.

  • Toro, E. F., 1989: A weighted average flux method for hyperbolic conservation laws. Proc. Roy. Soc. London, 423 , 401–418.

  • van Leer, B., 1979: Towards the ultimate conservative difference scheme V: A second order sequel to Godunov's method. J. Comput. Phys., 32 , 101–136.

    • Search Google Scholar
    • Export Citation

  • Williamson, D. L., and P. J. Rasch, 1989: Two-dimensional semi-Lagrangian transport with shape-preserving interpolation. Mon. Wea. Rev., 117 , 102–129.

    • Search Google Scholar
    • Export Citation

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

    • Search Google Scholar
    • Export Citation
All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 213 79 5
PDF Downloads 171 58 3