U.S. Department of Energy Office of Biological and Environmental Research

BER Research Highlights

Discrete Conservation Properties for Shallow Water Flows using Mixed Mimetic Spectral Elements
Published: December 20, 2017
Posted: April 10, 2018

New methods for simulating fluid dynamics in the atmosphere and ocean.

The Science
This paper explores a new method for the solution of the rotating shallow water equations - a simpler analog often used to test approaches for full general circulation models. This new method exactly preserves conservation properties of the fluid flow with arbitrarily high accuracy and is a strong candidate for use in climate system simulations at long timescales.

The Impact
The conservation and error convergence properties of this new approach make it a strong candidate for modeling geophysical dynamics over long time integrations with high accuracy and minimal bias. In particular, it should be well suited to the modeling of both gravity waves on fast time scales and Rossby waves on slow time scales that are responsible for many aspects of climate and weather dynamics.

To simulate fluid flow on a computer, we must translate the continuous fluid equations into a discrete form that can be calculated on a computational mesh covering the sphere. Mimetic methods are a technical term for a discrete formulation that can exactly capture properties of the continuous mathematical equations. Spectral element methods are a choice of discrete model that uses a series of carefully chosen mathematical functions, called basis functions, to provide a locally accurate approximation to the solution in each mesh cell (element). In this work, mimetic and spectral element approaches are combined and applied to a set of simpler fluid equations known as a shallow water model. The combined mixed mimetic spectral element method introduces a set of recently developed basis functions for the spectral element approach that exactly satisfy properties of the underlying mathematical operators (i.e., mimetic). This allows for the exact conservation of both simpler quantities (e.g., mass) as well as higher moments like energy or enstrophy, subject to the truncation error of the time stepping scheme. The combination of very high accuracy associated with spectral element methods and the preservation of conservation properties needed to represent the Earth system in long timescale simulations appears to provide a much improved method for solving the fluid equations for use in large-scale Earth system models. This work has provided an initial implementation of these idea and has demonstrated much improved results. Future work will integrate these schemes in the full dynamical cores of Earth system models.

Contacts (BER PM)
Dorothy Koch
Earth System Modeling

(PI Contact)
Phil Jones
Los Alamos National Laboratory

Funding was provided by the Earth System Modeling program within the DOE Office of Science Office of Biological and Environmental Research.

Lee, D., A. Palha, and M. Gerritsma. "Discrete conservation properties for shallow water flows using mixed mimetic spectral elements." Journal of Computational Physics 357, 282-304 (2017). [DOI:10.1016/j.jcp.2017.12.022]

Related Links
Launching an Extreme-scale ACME Prototype - Transport (LEAP-T)

Topic Areas:

  • Research Area: Earth and Environmental Systems Modeling

Division: SC-23.1 Climate and Environmental Sciences Division, BER


BER supports basic research and scientific user facilities to advance DOE missions in energy and environment. More about BER

Recent Highlights

Aug 24, 2019
New Approach for Studying How Microbes Influence Their Environment
A diverse group of scientists suggests a common framework and targeting of known microbial processes [more...]

Aug 08, 2019
Nutrient-Hungry Peatland Microbes Reduce Carbon Loss Under Warmer Conditions
Enzyme production in peatlands reduces carbon lost to respiration under future high temperatures. [more...]

Aug 05, 2019
Amazon Forest Response to CO2 Fertilization Dependent on Plant Phosphorus Acquisition
AmazonFACE Model Intercomparison. The Science Plant growth is dependent on the availabi [more...]

Jul 29, 2019
A Slippery Slope: Soil Carbon Destabilization
Carbon gain or loss depends on the balance between competing biological, chemical, and physical reac [more...]

Jul 15, 2019
Field Evaluation of Gas Analyzers for Measuring Ecosystem Fluxes
How gas analyzer type and correction method impact measured fluxes. The Science A side- [more...]

List all highlights (possible long download time)