Total energy and potential enstrophy conserving schemes for the shallow water equations using Hamiltonian methods – Part 1: Derivation and properties
Abstract
The shallow water equations provide a useful analogue of the fully compressible Euler equations since they have similar characteristics: conservation laws, inertia-gravity and Rossby waves, and a (quasi-) balanced state. In order to obtain realistic simulation results, it is desirable that numerical models have discrete analogues of these properties. Two prototypical examples of such schemes are the 1981 Arakawa and Lamb (AL81) C-grid total energy and potential enstrophy conserving scheme, and the 2007 Salmon (S07) Z-grid total energy and potential enstrophy conserving scheme. Unfortunately, the AL81 scheme is restricted to logically square, orthogonal grids, and the S07 scheme is restricted to uniform square grids. The current work extends the AL81 scheme to arbitrary non-orthogonal polygonal grids and the S07 scheme to arbitrary orthogonal spherical polygonal grids in a manner that allows for both total energy and potential enstrophy conservation, by combining Hamiltonian methods (work done by Salmon, Gassmann, Dubos, and others) and discrete exterior calculus (Thuburn, Cotter, Dubos, Ringler, Skamarock, Klemp, and others). Lastly, detailed results of the schemes applied to standard test cases are deferred to part 2 of this series of papers.
- Authors:
- Publication Date:
- Research Org.:
- Krell Institute, Ames, IA (United States)
- Sponsoring Org.:
- USDOE; National Science Foundation (NSF)
- OSTI Identifier:
- 1344140
- Alternate Identifier(s):
- OSTI ID: 1366511
- Grant/Contract Number:
- FG02-97ER25308
- Resource Type:
- Published Article
- Journal Name:
- Geoscientific Model Development (Online)
- Additional Journal Information:
- Journal Name: Geoscientific Model Development (Online) Journal Volume: 10 Journal Issue: 2; Journal ID: ISSN 1991-9603
- Publisher:
- Copernicus Publications, EGU
- Country of Publication:
- Germany
- Language:
- English
- Subject:
- 58 GEOSCIENCES
Citation Formats
Eldred, Christopher, and Randall, David. Total energy and potential enstrophy conserving schemes for the shallow water equations using Hamiltonian methods – Part 1: Derivation and properties. Germany: N. p., 2017.
Web. doi:10.5194/gmd-10-791-2017.
Eldred, Christopher, & Randall, David. Total energy and potential enstrophy conserving schemes for the shallow water equations using Hamiltonian methods – Part 1: Derivation and properties. Germany. https://doi.org/10.5194/gmd-10-791-2017
Eldred, Christopher, and Randall, David. Fri .
"Total energy and potential enstrophy conserving schemes for the shallow water equations using Hamiltonian methods – Part 1: Derivation and properties". Germany. https://doi.org/10.5194/gmd-10-791-2017.
@article{osti_1344140,
title = {Total energy and potential enstrophy conserving schemes for the shallow water equations using Hamiltonian methods – Part 1: Derivation and properties},
author = {Eldred, Christopher and Randall, David},
abstractNote = {The shallow water equations provide a useful analogue of the fully compressible Euler equations since they have similar characteristics: conservation laws, inertia-gravity and Rossby waves, and a (quasi-) balanced state. In order to obtain realistic simulation results, it is desirable that numerical models have discrete analogues of these properties. Two prototypical examples of such schemes are the 1981 Arakawa and Lamb (AL81) C-grid total energy and potential enstrophy conserving scheme, and the 2007 Salmon (S07) Z-grid total energy and potential enstrophy conserving scheme. Unfortunately, the AL81 scheme is restricted to logically square, orthogonal grids, and the S07 scheme is restricted to uniform square grids. The current work extends the AL81 scheme to arbitrary non-orthogonal polygonal grids and the S07 scheme to arbitrary orthogonal spherical polygonal grids in a manner that allows for both total energy and potential enstrophy conservation, by combining Hamiltonian methods (work done by Salmon, Gassmann, Dubos, and others) and discrete exterior calculus (Thuburn, Cotter, Dubos, Ringler, Skamarock, Klemp, and others). Lastly, detailed results of the schemes applied to standard test cases are deferred to part 2 of this series of papers.},
doi = {10.5194/gmd-10-791-2017},
journal = {Geoscientific Model Development (Online)},
number = 2,
volume = 10,
place = {Germany},
year = {Fri Feb 17 00:00:00 EST 2017},
month = {Fri Feb 17 00:00:00 EST 2017}
}
https://doi.org/10.5194/gmd-10-791-2017
Web of Science
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
A Framework for Mimetic Discretization of the Rotating Shallow-Water Equations on Arbitrary Polygonal Grids
journal, January 2012
- Thuburn, J.; Cotter, C. J.
- SIAM Journal on Scientific Computing, Vol. 34, Issue 3
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
Hamiltonian Fluid Mechanics
journal, January 1988
- Salmon, R.
- Annual Review of Fluid Mechanics, Vol. 20, Issue 1
Computational Modes and Grid Imprinting on Five Quasi-Uniform Spherical C Grids
journal, August 2012
- Weller, Hilary; Thuburn, John; Cotter, Colin J.
- Monthly Weather Review, Vol. 140, Issue 8
Numerical instabilities of vector-invariant momentum equations on rectangular C-grids: Vector-Invariant Momentum Equation Instabilities
journal, December 2016
- Bell, Michael J.; Peixoto, Pedro S.; Thuburn, John
- Quarterly Journal of the Royal Meteorological Society, Vol. 143, Issue 702
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
Why starting from differential equations for computational physics?
journal, January 2014
- Tonti, Enzo
- Journal of Computational Physics, Vol. 257
An initial-value problem for testing numerical models of the global shallow-water equations
journal, October 2004
- Galewsky, Joseph; Scott, Richard K.; Polvani, Lorenzo M.
- Tellus A, Vol. 56, Issue 5
Horizontal grids for global weather and climate prediction models: a review
journal, November 2011
- Staniforth, Andrew; Thuburn, John
- Quarterly Journal of the Royal Meteorological Society, Vol. 138, Issue 662
A mimetic, semi-implicit, forward-in-time, finite volume shallow water model: comparison of hexagonal–icosahedral and cubed-sphere grids
journal, January 2014
- Thuburn, J.; Cotter, C. J.; Dubos, T.
- Geoscientific Model Development, Vol. 7, Issue 3
Energy-conserving finite-difference schemes for quasi-hydrostatic equations: Energy-Conserving Schemes for QHEs
journal, August 2015
- Tort, Marine; Dubos, Thomas; Melvin, Thomas
- Quarterly Journal of the Royal Meteorological Society, Vol. 141, Issue 693
Energy- and enstrophy-conserving schemes for the shallow-water equations, based on mimetic finite elements: Energy- and enstrophy-conserving schemes
journal, February 2014
- McRae, Andrew T. T.; Cotter, Colin J.
- Quarterly Journal of the Royal Meteorological Society, Vol. 140, Issue 684
A finite element exterior calculus framework for the rotating shallow-water equations
journal, January 2014
- Cotter, C. J.; Thuburn, J.
- Journal of Computational Physics, Vol. 257
Optimized Icosahedral Grids: Performance of Finite-Difference Operators and Multigrid Solver
journal, December 2013
- Heikes, Ross P.; Randall, David A.; Konor, Celal S.
- Monthly Weather Review, Vol. 141, Issue 12
Some conservation issues for the dynamical cores of NWP and climate models
journal, March 2008
- Thuburn, J.
- Journal of Computational Physics, Vol. 227, Issue 7
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
An internal symmetric computational instability
journal, April 1983
- Hollingsworth, A.; Kållberg, P.; Renner, V.
- Quarterly Journal of the Royal Meteorological Society, Vol. 109, Issue 460
Controlling the Computational Modes of the Arbitrarily Structured C Grid
journal, October 2012
- Weller, Hilary
- Monthly Weather Review, Vol. 140, Issue 10
An initial-value problem for testing numerical models of the global shallow-water equations
journal, January 2004
- Galewsky, Joseph; Scott, Richard K.; Polvani, Lorenzo M.
- Tellus A: Dynamic Meteorology and Oceanography, Vol. 56, Issue 5
Statistical mechanics of Arakawa’s discretizations
journal, December 2007
- Dubinkina, Svetlana; Frank, Jason
- Journal of Computational Physics, Vol. 227, Issue 2
A general method for conserving quantities related to potential vorticity in numerical models
journal, June 2005
- Salmon, Rick
- Nonlinearity, Vol. 18, Issue 5
A unified theory of available potential energy 1
journal, March 1993
- Shepherd, Theodore G.
- Atmosphere-Ocean, Vol. 31, Issue 1
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
Differential forms for scientists and engineers
journal, January 2014
- Blair Perot, J.; Zusi, Christopher J.
- Journal of Computational Physics, Vol. 257
A two-dimensional mixed finite-element pair on rectangles: 2D Mixed Finite-Element Pair
journal, July 2013
- Melvin, Thomas; Staniforth, Andrew; Cotter, Colin
- Quarterly Journal of the Royal Meteorological Society, Vol. 140, Issue 680
Energy–Vorticity Theory of Ideal Fluid Mechanics
journal, July 2009
- Névir, Peter; Sommer, Matthias
- Journal of the Atmospheric Sciences, Vol. 66, Issue 7
Analysis of a mixed finite-element pair proposed for an atmospheric dynamical core
journal, January 2013
- Staniforth, Andrew; Melvin, Thomas; Cotter, Colin
- Quarterly Journal of the Royal Meteorological Society, Vol. 139, Issue 674
Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid
journal, January 2014
- Weller, H.
- Geoscientific Model Development, Vol. 7, Issue 3
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
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
A conservative scheme for the shallow-water system on a staggered geodesic grid based on a Nambu representation
journal, January 2009
- Sommer, Matthias; Névir, Peter
- Quarterly Journal of the Royal Meteorological Society, Vol. 135, Issue 639