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

Title: 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}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record
https://doi.org/10.5194/gmd-10-791-2017

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

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

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

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

Hamiltonian Fluid Mechanics
journal, January 1988


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
  • DOI: 10.1175/MWR-D-11-00193.1

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

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

Why starting from differential equations for computational physics?
journal, January 2014


An initial-value problem for testing numerical models of the global shallow-water equations
journal, October 2004


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

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
  • DOI: 10.5194/gmd-7-909-2014

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

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

A finite element exterior calculus framework for the rotating shallow-water equations
journal, January 2014


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
  • DOI: 10.1175/MWR-D-12-00236.1

Some conservation issues for the dynamical cores of NWP and climate models
journal, March 2008


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

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

Controlling the Computational Modes of the Arbitrarily Structured C Grid
journal, October 2012


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
  • DOI: 10.3402/tellusa.v56i5.14436

Statistical mechanics of Arakawa’s discretizations
journal, December 2007


A unified theory of available potential energy 1
journal, March 1993


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


Differential forms for scientists and engineers
journal, January 2014


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

Energy–Vorticity Theory of Ideal Fluid Mechanics
journal, July 2009

  • Névir, Peter; Sommer, Matthias
  • Journal of the Atmospheric Sciences, Vol. 66, Issue 7
  • DOI: 10.1175/2008JAS2897.1

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

Non-orthogonal version of the arbitrary polygonal C-grid and a new diamond grid
journal, January 2014


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


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

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