Article Type:
Research Article

Trajectory Calculations for Spherical Geodesic Grids in Cartesian Space

Francis X. Giraldo Naval Research Laboratory, Monterey, California

Search for other papers by Francis X. Giraldo in
Current site
Google Scholar
PubMed Close
Print Publication:
01 Jul 1999
DOI:
https://doi.org/10.1175/1520-0493(1999)127<1651:TCFSGG>2.0.CO;2
Page(s):
1651–1662
Restricted access

Abstract

This paper shows how to obtain accurate and efficient trajectory calculations for spherical geodesic grids in Cartesian space. Determination of the departure points is essential to characteristic-based methods that trace the value of a function to the foot of the characteristics and then either integrate or interpolate at this location. In this paper, the departure points are all computed in relation to the spherical geodesic grids that are composed of a disjoint set of unstructured equilateral triangles. Interpolating and noninterpolating trajectory calculation approaches are both illustrated and the accuracy of both methods are compared. The noninterpolating method of McGregor results in the most accurate trajectories. The challenge in using McGregor’s method on unstructured triangular grids lies in the computation of the derivatives required in the high-order terms of the Taylor series expansion. This paper extends McGregor’s method to unstructured triangular grids by describing an accurate and efficient method for constructing the derivatives in an element by element approach typical of finite element methods. An order of accuracy analysis reveals that these numerical derivatives are second-order accurate.

Corresponding author address: Francis X. Giraldo, Naval Research Laboratory, 7 Grace Hopper Ave., Monterey, CA 93943.

Email: giraldo@nrlmry.navy.mil

Abstract

This paper shows how to obtain accurate and efficient trajectory calculations for spherical geodesic grids in Cartesian space. Determination of the departure points is essential to characteristic-based methods that trace the value of a function to the foot of the characteristics and then either integrate or interpolate at this location. In this paper, the departure points are all computed in relation to the spherical geodesic grids that are composed of a disjoint set of unstructured equilateral triangles. Interpolating and noninterpolating trajectory calculation approaches are both illustrated and the accuracy of both methods are compared. The noninterpolating method of McGregor results in the most accurate trajectories. The challenge in using McGregor’s method on unstructured triangular grids lies in the computation of the derivatives required in the high-order terms of the Taylor series expansion. This paper extends McGregor’s method to unstructured triangular grids by describing an accurate and efficient method for constructing the derivatives in an element by element approach typical of finite element methods. An order of accuracy analysis reveals that these numerical derivatives are second-order accurate.

Corresponding author address: Francis X. Giraldo, Naval Research Laboratory, 7 Grace Hopper Ave., Monterey, CA 93943.

Email: giraldo@nrlmry.navy.mil

Save
Cover Monthly Weather Review
Monthly Weather Review

  • Côté, J., 1988: A Lagrange multiplier approach for the metric terms of semi-Lagrangian models on the sphere. Quart. J. Roy. Meteor. Soc.,114, 1347–1352.

  • Giraldo, F. X., 1997: Lagrange–Galerkin methods on spherical geodesic grids. J. Comput. Phys.,136, 197–217.

  • Le Roux, D. Y., C. A. Lin, and A. Staniforth, 1997: An accurate interpolating scheme for semi-Lagrangian advection on an unstructured mesh for ocean modelling. Tellus,49A, 119–138.

  • McDonald, A., and J. R. Bates, 1989: Semi-Lagrangian integration of a gridpoint shallow water model on the sphere. Mon. Wea. Rev.,117, 130–137.

  • McGregor, J. L., 1993: Economical determination of departure points for semi-Lagrangian models. Mon. Wea. Rev.,121, 221–230.

  • Ritchie, H., 1987: Semi-Lagrangian advection on a Gaussian grid. Mon. Wea. Rev.,115, 608–619.

  • Silvester, P., 1969: High-order polynomial triangular finite elements for potential problems. Int. J. Eng. Sci.,7, 849–861.

  • 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.

All Time Past Year Past 30 Days
Abstract Views 0 0 0
Full Text Views 319 88 9
PDF Downloads 185 51 6