DOE PAGES title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: 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:
 [1];  [1];  [1];  [2];  [3];  [4]
  1. Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)
  2. Univ. of California, Davis, CA (United States)
  3. NVIDIA, Santa Clara, CA (United States)
  4. 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 = {Sat Dec 28 00:00:00 EST 2019},
month = {Sat Dec 28 00:00:00 EST 2019}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 21 works
Citation information provided by
Web of Science

Figures / Tables:

Figure 1 Figure 1: The terrain following s coordinate levels and level indexing. There are n midpoint levels indexed by i = 1, 2, ... , n and n + 1 interface levels indexed by i + $^{1}_{2}$, i = 0, 1, ... , n.

Save / Share:

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
  • DOI: 10.1002/qj.297

Optimization-based limiters for the spectral element method
journal, June 2014


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
  • DOI: 10.1175/JAS3837.1

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
  • DOI: 10.1029/2018MS001603

Various Vertical Coordinate Systems Used for Numerical Weather Prediction
journal, July 1974


Energy and Numerical Weather Prediction
journal, November 1960


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
  • DOI: 10.1002/qj.1906

Support-Operator Finite-Difference Algorithms for General Elliptic Problems
journal, April 1995

  • Shashkov, Mikhail; Steinberg, Stanly
  • Journal of Computational Physics, Vol. 118, Issue 1
  • DOI: 10.1006/jcph.1995.1085

Direct Discretization of Planar Div-Curl Problems
journal, February 1992

  • Nicolaides, R. A.
  • SIAM Journal on Numerical Analysis, Vol. 29, Issue 1
  • DOI: 10.1137/0729003

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
  • DOI: 10.1029/2019MS001629

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
  • DOI: 10.1002/qj.1960

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
  • DOI: 10.1002/fld.1154

Dispersive behaviour of high order finite element schemes for the one-way wave equation
journal, February 2014


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
  • DOI: 10.1006/jcph.1996.0047

A primal–dual mimetic finite element scheme for the rotating shallow water equations on polygonal spherical meshes
journal, June 2015


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
  • DOI: 10.1177/1094342011428142

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
  • DOI: 10.1029/2011MS000099

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
  • DOI: 10.5194/gmdd-7-4081-2014

Runge–Kutta IMEX schemes for the Horizontally Explicit/Vertically Implicit (HEVI) solution of wave equations
journal, November 2013


A high-order staggered finite-element vertical discretization for non-hydrostatic atmospheric models
journal, January 2016


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
  • DOI: 10.1177/1094342005056108

On the use of a coordinate transformation for the solution of the Navier-Stokes equations
journal, February 1975


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
  • DOI: 10.1137/120876034

DCMIP2016: the splitting supercell test case
journal, January 2019

  • Zarzycki, Colin M.; Jablonowski, Christiane; Kent, James
  • Geoscientific Model Development, Vol. 12, Issue 3
  • DOI: 10.5194/gmd-12-879-2019

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
  • DOI: 10.1016/j.jcp.2009.12.007

Operator-Split Runge–Kutta–Rosenbrock Methods for Nonhydrostatic Atmospheric Models
journal, April 2012


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
  • DOI: 10.5194/gmd-11-1497-2018

Mixed finite elements for numerical weather prediction
journal, August 2012


Equations of Atmospheric Motion in Non-Eulerian Vertical Coordinates: Vector-Invariant Form and Quasi-Hamiltonian Formulation
journal, October 2014


A Terrain-Following Coordinate with Smoothed Coordinate Surfaces
journal, July 2011


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
  • DOI: 10.1002/qj.2246

A New Terrain-Following Vertical Coordinate Formulation for Atmospheric Prediction Models
journal, October 2002


Conservation of Mass and Energy for the Moist Atmospheric Primitive Equations on Unstructured Grids
book, January 2011


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
  • DOI: 10.5194/gmd-10-4477-2017

A compatible and conservative spectral element method on unstructured grids
journal, August 2010


A Generalization of the SLEVE Vertical Coordinate
journal, September 2010

  • Leuenberger, Daniel; Koller, Marcel; Fuhrer, Oliver
  • Monthly Weather Review, Vol. 138, Issue 9
  • DOI: 10.1175/2010MWR3307.1

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
  • DOI: 10.1016/j.jcp.2018.06.035

DYNAMICO-1.0, an icosahedral hydrostatic dynamical core designed for consistency and versatility
journal, January 2015


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
  • DOI: 10.1002/qj.2241

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
  • DOI: 10.1145/1089014.1089020

A “Vertically Lagrangian” Finite-Volume Dynamical Core for Global Models
journal, October 2004


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
  • DOI: 10.1016/j.jcp.2009.08.006

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
  • DOI: 10.1002/qj.3294

Figures/Tables have been extracted from DOE-funded journal article accepted manuscripts.