An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables
Abstract
We derive a formulation of the nonhydrostatic equations in spherical geometry with a Lorenz staggered vertical discretization. The combination conserves a discrete energy in exact time integration when coupled with a mimetic horizontal discretization. The formulation is a version of Dubos and Tort (2014, https://doi.org/10.1175/MWR-D-14-00069.1) rewritten in terms of primitive variables. It is valid for terrain following mass or height coordinates and for both Eulerian or vertically Lagrangian discretizations. The discretization relies on an extension to Simmons and Burridge (1981, https://doi.org/10.1175/1520-0493(1981)109<0758:AEAAMC>2.0.CO;2) vertical differencing, which we show obeys a discrete derivative product rule. This product rule allows us to simplify the treatment of the vertical transport terms. Energy conservation is obtained via a term-by-term balance in the kinetic, internal, and potential energy budgets, ensuring an energy-consistent discretization up to time truncation error with no spurious sources of energy. We demonstrate convergence with respect to time truncation error in a spectral element code with a horizontal explicit vertically implicit implicit-explicit time stepping algorithm.
- Authors:
-
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Univ. of California, Davis, CA (United States)
- NVIDIA, Santa Clara, CA (United States)
- Univ. Grenoble Alpes (France)
- Publication Date:
- Research Org.:
- Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
- Sponsoring Org.:
- USDOE Office of Science (SC), Biological and Environmental Research (BER); USDOE National Nuclear Security Administration (NNSA); French National Research Agency
- OSTI Identifier:
- 1631000
- Report Number(s):
- SAND2020-5242J
Journal ID: ISSN 1942-2466; 686193
- Grant/Contract Number:
- AC04-94AL85000; NA0003525; ANR-14-CE23-0010
- Resource Type:
- Accepted Manuscript
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Additional Journal Information:
- Journal Volume: 12; Journal Issue: 1; Journal ID: ISSN 1942-2466
- Publisher:
- American Geophysical Union (AGU)
- Country of Publication:
- United States
- Language:
- English
- Subject:
- 97 MATHEMATICS AND COMPUTING; nonhydrostatic; hamiltonian; dynamical core; energy conservation; mimetric
Citation Formats
Taylor, Mark A., Guba, Oksana, Steyer, Andrew, Ullrich, Paul A., Hall, David M., and Eldrid, Christopher. An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables. United States: N. p., 2019.
Web. doi:10.1029/2019MS001783.
Taylor, Mark A., Guba, Oksana, Steyer, Andrew, Ullrich, Paul A., Hall, David M., & Eldrid, Christopher. An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables. United States. https://doi.org/10.1029/2019MS001783
Taylor, Mark A., Guba, Oksana, Steyer, Andrew, Ullrich, Paul A., Hall, David M., and Eldrid, Christopher. Sat .
"An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables". United States. https://doi.org/10.1029/2019MS001783. https://www.osti.gov/servlets/purl/1631000.
@article{osti_1631000,
title = {An Energy Consistent Discretization of the Nonhydrostatic Equations in Primitive Variables},
author = {Taylor, Mark A. and Guba, Oksana and Steyer, Andrew and Ullrich, Paul A. and Hall, David M. and Eldrid, Christopher},
abstractNote = {We derive a formulation of the nonhydrostatic equations in spherical geometry with a Lorenz staggered vertical discretization. The combination conserves a discrete energy in exact time integration when coupled with a mimetic horizontal discretization. The formulation is a version of Dubos and Tort (2014, https://doi.org/10.1175/MWR-D-14-00069.1) rewritten in terms of primitive variables. It is valid for terrain following mass or height coordinates and for both Eulerian or vertically Lagrangian discretizations. The discretization relies on an extension to Simmons and Burridge (1981, https://doi.org/10.1175/1520-0493(1981)109<0758:AEAAMC>2.0.CO;2) vertical differencing, which we show obeys a discrete derivative product rule. This product rule allows us to simplify the treatment of the vertical transport terms. Energy conservation is obtained via a term-by-term balance in the kinetic, internal, and potential energy budgets, ensuring an energy-consistent discretization up to time truncation error with no spurious sources of energy. We demonstrate convergence with respect to time truncation error in a spectral element code with a horizontal explicit vertically implicit implicit-explicit time stepping algorithm.},
doi = {10.1029/2019MS001783},
journal = {Journal of Advances in Modeling Earth Systems},
number = 1,
volume = 12,
place = {United States},
year = {2019},
month = {12}
}
Web of Science
Figures / Tables:

Works referenced in this record:
Towards a consistent numerical compressible non-hydrostatic model using generalized Hamiltonian tools
journal, July 2008
- Gassmann, Almut; Herzog, Hans-Joachim
- Quarterly Journal of the Royal Meteorological Society, Vol. 134, Issue 635
Optimization-based limiters for the spectral element method
journal, June 2014
- Guba, Oksana; Taylor, Mark; St-Cyr, Amik
- Journal of Computational Physics, Vol. 267
A General Method for Conserving Energy and Potential Enstrophy in Shallow-Water Models
journal, February 2007
- Salmon, Rick
- Journal of the Atmospheric Sciences, Vol. 64, Issue 2
The DOE E3SM Coupled Model Version 1: Overview and Evaluation at Standard Resolution
journal, July 2019
- Golaz, Jean‐Christophe; Caldwell, Peter M.; Van Roekel, Luke P.
- Journal of Advances in Modeling Earth Systems, Vol. 11, Issue 7
Poisson-Bracket Approach to the Construction of Energy- and Potential-Enstrophy-Conserving Algorithms for the Shallow-Water Equations
journal, August 2004
- Salmon, Rick
- Journal of the Atmospheric Sciences, Vol. 61, Issue 16
Various Vertical Coordinate Systems Used for Numerical Weather Prediction
journal, July 1974
- Kasahara, Akira
- Monthly Weather Review, Vol. 102, Issue 7
Energy and Numerical Weather Prediction
journal, November 1960
- Lorenz, Edward N.
- Tellus, Vol. 12, Issue 4
Dispersion analysis of the spectral element method
journal, February 2012
- Melvin, Thomas; Staniforth, Andrew; Thuburn, John
- Quarterly Journal of the Royal Meteorological Society, Vol. 138, Issue 668
Support-Operator Finite-Difference Algorithms for General Elliptic Problems
journal, April 1995
- Shashkov, Mikhail; Steinberg, Stanly
- Journal of Computational Physics, Vol. 118, Issue 1
Direct Discretization of Planar Div-Curl Problems
journal, February 1992
- Nicolaides, R. A.
- SIAM Journal on Numerical Analysis, Vol. 29, Issue 1
An Overview of the Atmospheric Component of the Energy Exascale Earth System Model
journal, August 2019
- Rasch, P. J.; Xie, S.; Ma, P. ‐L.
- Journal of Advances in Modeling Earth Systems, Vol. 11, Issue 8
A global hexagonal C-grid non-hydrostatic dynamical core (ICON-IAP) designed for energetic consistency
journal, June 2012
- Gassmann, Almut
- Quarterly Journal of the Royal Meteorological Society, Vol. 139, Issue 670
The Parabolic Spline Method (PSM) for conservative transport problems
journal, January 2006
- Zerroukat, M.; Wood, N.; Staniforth, A.
- International Journal for Numerical Methods in Fluids, Vol. 51, Issue 11
Dispersive behaviour of high order finite element schemes for the one-way wave equation
journal, February 2014
- Ainsworth, Mark
- Journal of Computational Physics, Vol. 259
The “Cubed Sphere”: A New Method for the Solution of Partial Differential Equations in Spherical Geometry
journal, March 1996
- Ronchi, C.; Iacono, R.; Paolucci, P. S.
- Journal of Computational Physics, Vol. 124, Issue 1
A primal–dual mimetic finite element scheme for the rotating shallow water equations on polygonal spherical meshes
journal, June 2015
- Thuburn, John; Cotter, Colin J.
- Journal of Computational Physics, Vol. 290
CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model
journal, November 2011
- Dennis, John M.; Edwards, Jim; Evans, Katherine J.
- The International Journal of High Performance Computing Applications, Vol. 26, Issue 1
Idealized tropical cyclone simulations of intermediate complexity: A test case for AGCMs
journal, February 2012
- Reed, Kevin A.; Jablonowski, Christiane
- Journal of Advances in Modeling Earth Systems, Vol. 4, Issue 2
The spectral element method on variable resolution grids: evaluating grid sensitivity and resolution-aware numerical viscosity
journal, January 2014
- Guba, O.; Taylor, M. A.; Ullrich, P. A.
- Geoscientific Model Development Discussions, Vol. 7, Issue 3
The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable
journal, January 1992
- Laprise, René
- Monthly Weather Review, Vol. 120, Issue 1
Runge–Kutta IMEX schemes for the Horizontally Explicit/Vertically Implicit (HEVI) solution of wave equations
journal, November 2013
- Weller, Hilary; Lock, Sarah-Jane; Wood, Nigel
- Journal of Computational Physics, Vol. 252
An Energy and Angular-Momentum Conserving Vertical Finite-Difference Scheme and Hybrid Vertical Coordinates
journal, April 1981
- Simmons, A. J.; Burridge, D. M.
- Monthly Weather Review, Vol. 109, Issue 4
A high-order staggered finite-element vertical discretization for non-hydrostatic atmospheric models
journal, January 2016
- Guerra, Jorge E.; Ullrich, Paul A.
- Geoscientific Model Development, Vol. 9, Issue 5
High-Resolution Mesh Convergence Properties and Parallel Efficiency of a Spectral Element Atmospheric Dynamical Core
journal, August 2005
- Dennis, John; Fournier, Aimé; Spotz, William F.
- The International Journal of High Performance Computing Applications, Vol. 19, Issue 3
On the use of a coordinate transformation for the solution of the Navier-Stokes equations
journal, February 1975
- Gal-Chen, Tzvi; Somerville, Richard C. J.
- Journal of Computational Physics, Vol. 17, Issue 2
Implicit-Explicit Formulations of a Three-Dimensional Nonhydrostatic Unified Model of the Atmosphere (NUMA)
journal, January 2013
- Giraldo, F. X.; Kelly, J. F.; Constantinescu, E. M.
- SIAM Journal on Scientific Computing, Vol. 35, Issue 5
DCMIP2016: the splitting supercell test case
journal, January 2019
- Zarzycki, Colin M.; Jablonowski, Christiane; Kent, James
- Geoscientific Model Development, Vol. 12, Issue 3
A unified approach to energy conservation and potential vorticity dynamics for arbitrarily-structured C-grids
journal, May 2010
- Ringler, T. D.; Thuburn, J.; Klemp, J. B.
- Journal of Computational Physics, Vol. 229, Issue 9
Operator-Split Runge–Kutta–Rosenbrock Methods for Nonhydrostatic Atmospheric Models
journal, April 2012
- Ullrich, Paul; Jablonowski, Christiane
- Monthly Weather Review, Vol. 140, Issue 4
Implicit–explicit (IMEX) Runge–Kutta methods for non-hydrostatic atmospheric models
journal, January 2018
- Gardner, David J.; Guerra, Jorge E.; Hamon, François P.
- Geoscientific Model Development, Vol. 11, Issue 4
Mixed finite elements for numerical weather prediction
journal, August 2012
- Cotter, C. J.; Shipton, J.
- Journal of Computational Physics, Vol. 231, Issue 21
Equations of Atmospheric Motion in Non-Eulerian Vertical Coordinates: Vector-Invariant Form and Quasi-Hamiltonian Formulation
journal, October 2014
- Dubos, Thomas; Tort, Marine
- Monthly Weather Review, Vol. 142, Issue 10
A Terrain-Following Coordinate with Smoothed Coordinate Surfaces
journal, July 2011
- Klemp, J. B.
- Monthly Weather Review, Vol. 139, Issue 7
Numerical analyses of Runge-Kutta implicit-explicit schemes for horizontally explicit, vertically implicit solutions of atmospheric models: Numerical Analyses of RK HEVI Solutions
journal, February 2014
- Lock, S. -J.; Wood, N.; Weller, H.
- Quarterly Journal of the Royal Meteorological Society, Vol. 140, Issue 682
A New Terrain-Following Vertical Coordinate Formulation for Atmospheric Prediction Models
journal, October 2002
- Schär, Christoph; Leuenberger, Daniel; Fuhrer, Oliver
- Monthly Weather Review, Vol. 130, Issue 10
DCMIP2016: a review of non-hydrostatic dynamical core design and intercomparison of participating models
journal, January 2017
- Ullrich, Paul A.; Jablonowski, Christiane; Kent, James
- Geoscientific Model Development, Vol. 10, Issue 12
A compatible and conservative spectral element method on unstructured grids
journal, August 2010
- Taylor, Mark A.; Fournier, Aimé
- Journal of Computational Physics, Vol. 229, Issue 17
A Generalization of the SLEVE Vertical Coordinate
journal, September 2010
- Leuenberger, Daniel; Koller, Marcel; Fuhrer, Oliver
- Monthly Weather Review, Vol. 138, Issue 9
Impact and importance of hyperdiffusion on the spectral element method: A linear dispersion analysis
journal, December 2018
- Ullrich, Paul A.; Reynolds, Daniel R.; Guerra, Jorge E.
- Journal of Computational Physics, Vol. 375
Conservative Scheme for the Compressible Nonhydrostatic Models with the Horizontally Explicit and Vertically Implicit Time Integration Scheme
journal, May 2002
- Satoh, Masaki
- Monthly Weather Review, Vol. 130, Issue 5
DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility
journal, January 2015
- Dubos, T.; Dubey, S.; Tort, M.
- Geoscientific Model Development, Vol. 8, Issue 10
A proposed baroclinic wave test case for deep- and shallow-atmosphere dynamical cores: A Baroclinic Wave Test Case for Deep and Shallow-atmosphere Dynamical Cores
journal, December 2013
- Ullrich, Paul A.; Melvin, Thomas; Jablonowski, Christiane
- Quarterly Journal of the Royal Meteorological Society, Vol. 140, Issue 682
SUNDIALS: Suite of nonlinear and differential/algebraic equation solvers
journal, September 2005
- Hindmarsh, Alan C.; Brown, Peter N.; Grant, Keith E.
- ACM Transactions on Mathematical Software, Vol. 31, Issue 3
Numerical representation of geostrophic modes on arbitrarily structured C-grids
journal, December 2009
- Thuburn, J.; Ringler, T. D.; Skamarock, W. C.
- Journal of Computational Physics, Vol. 228, Issue 22
Discretization of generalized Coriolis and friction terms on the deformed hexagonal C-grid
journal, October 2018
- Gassmann, Almut
- Quarterly Journal of the Royal Meteorological Society, Vol. 144, Issue 716
Figures / Tables found in this record: