User Manual Privacy Policy Disclaimer Contact us
  Advanced SearchBrowse




Journal Article

Conservative transport schemes for spherical geodesic grids: High-order flux operators for ODE-based time integration


Gassmann,  A.
Director’s Research Group OES, The Ocean in the Earth System, MPI for Meteorology, Max Planck Society;

There are no locators available
Fulltext (public)
There are no public fulltexts available
Supplementary Material (public)
There is no public supplementary material available

Skamarock, W., & Gassmann, A. (2011). Conservative transport schemes for spherical geodesic grids: High-order flux operators for ODE-based time integration. Monthly Weather Review, 139, 2962-2975. doi:10.1175/MWR-D-10-05056.1.

Cite as: http://hdl.handle.net/11858/00-001M-0000-0012-2534-F
Higher-order finite-volume flux operators for transport algorithms used within Runge-Kutta time integration schemes on irregular Voronoi (hexagonal) meshes are proposed and tested. These operators are generalizations of third- and fourth-order operators currently used in atmospheric models employing regular, orthogonal rectangular meshes. Two-dimensional least squares fit polynomials are used to evaluate the higher-order spatial derivatives needed to cancel the leading-order truncation error terms of the standard second-order centered formulation. Positive definite or monotonic behavior is achieved by applying an appropriate limiter during the final Runge-Kutta stage within a given time step. The third- and fourth-order formulations are evaluated using standard transport tests on the sphere. The new schemes are more accurate and significantly more efficient than the standard second-order scheme and other schemes in the literature examined by the authors. The third-order formulation is equivalent to the fourth-order formulation plus an additional diffusion term that is proportional to the Courant number. An optimal value for the coefficient scaling this diffusion term is chosen based on qualitative evaluation of the test results. Improvements using the higher-order scheme in place of the traditional second-order centered approach are illustrated within 3D unstable baroclinic wave simulations produced using two global non-hydrostatic models employing spherical Voronoi meshes.