A framework to evaluate IMEX schemes for atmospheric models

Guba, Oksana; Taylor, Mark A.; Bradley, Andrew M.; Bosler, Peter A.; Steyer, Andrew

doi:10.5194/gmd-13-6467-2020

Title: A framework to evaluate IMEX schemes for atmospheric models

Abstract

We present a new evaluation framework for implicit and explicit (IMEX) Runge–Kutta time-stepping schemes. The new framework uses a linearized nonhydrostatic system of normal modes. We utilize the framework to investigate the stability of IMEX methods and their dispersion and dissipation of gravity, Rossby, and acoustic waves. We test the new framework on a variety of IMEX schemes and use it to develop and analyze a set of second-order low-storage IMEX Runge–Kutta methods with a high Courant–Friedrichs–Lewy (CFL) number. We show that the new framework is more selective than the 2-D acoustic system previously used in the literature. Schemes that are stable for the 2-D acoustic system are not stable for the system of normal modes.

Authors:

^[1];

^[1]; Bradley, Andrew M. ^[1]; Bosler, Peter A. ^[1]; Steyer, Andrew ^[1]

Sandia National Lab. (SNL-NM), Albuquerque, NM (United States)

Publication Date:: Tue Dec 22 00:00:00 EST 2020

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)

OSTI Identifier:: 1765755

Report Number(s):: SAND-2020-5447J
Journal ID: ISSN 1991-9603; 686315

Grant/Contract Number:: AC04-94AL85000; NA0003525

Resource Type:: Accepted Manuscript

Journal Name:: Geoscientific Model Development (Online)

Additional Journal Information:: Journal Name: Geoscientific Model Development (Online); Journal Volume: 13; Journal Issue: 12; Journal ID: ISSN 1991-9603

Publisher:: European Geosciences Union

Country of Publication:: United States

Language:: English

Subject:: 54 ENVIRONMENTAL SCIENCES

Citation Formats


                    Guba, Oksana, Taylor, Mark A., Bradley, Andrew M., Bosler, Peter A., and Steyer, Andrew. A framework to evaluate IMEX schemes for atmospheric models.  United States: N. p., 2020. 
Web.  doi:10.5194/gmd-13-6467-2020.

Copy to clipboard


                    Guba, Oksana, Taylor, Mark A., Bradley, Andrew M., Bosler, Peter A., & Steyer, Andrew. A framework to evaluate IMEX schemes for atmospheric models.  United States.  https://doi.org/10.5194/gmd-13-6467-2020

Copy to clipboard


                    Guba, Oksana, Taylor, Mark A., Bradley, Andrew M., Bosler, Peter A., and Steyer, Andrew. Tue .  
"A framework to evaluate IMEX schemes for atmospheric models".  United States.  https://doi.org/10.5194/gmd-13-6467-2020.  https://www.osti.gov/servlets/purl/1765755.

Copy to clipboard


                    
@article{osti_1765755,

  title        = {A framework to evaluate IMEX schemes for atmospheric models},

  author       = {Guba, Oksana and Taylor, Mark A. and Bradley, Andrew M. and Bosler, Peter A. and Steyer, Andrew},

  abstractNote = {We present a new evaluation framework for implicit and explicit (IMEX) Runge–Kutta time-stepping schemes. The new framework uses a linearized nonhydrostatic system of normal modes. We utilize the framework to investigate the stability of IMEX methods and their dispersion and dissipation of gravity, Rossby, and acoustic waves. We test the new framework on a variety of IMEX schemes and use it to develop and analyze a set of second-order low-storage IMEX Runge–Kutta methods with a high Courant–Friedrichs–Lewy (CFL) number. We show that the new framework is more selective than the 2-D acoustic system previously used in the literature. Schemes that are stable for the 2-D acoustic system are not stable for the system of normal modes.},

  doi          = {10.5194/gmd-13-6467-2020},

  journal      = {Geoscientific Model Development (Online)},

  number       = 12,

  volume       = 13,

  place        = {United States},

  year         = {Tue Dec 22 00:00:00 EST 2020},

  month        = {Tue Dec 22 00:00:00 EST 2020}

}

Copy to clipboard

Journal Article:

Free Publicly Available Full Text

Accepted Manuscript (DOE)

Publisher's Version of Record

https://doi.org/10.5194/gmd-13-6467-2020

Other availability

Search WorldCat to find libraries that may hold this journal

Save / Share:

Export Metadata

Save to My Library

Works referenced in this record:

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

Implicit–Explicit Multistep Methods for Fast-Wave–Slow-Wave Problems
journal, April 2012

Durran, Dale R.; Blossey, Peter N.
Monthly Weather Review, Vol. 140, Issue 4
DOI: 10.1175/MWR-D-11-00088.1

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

One step integration methods with maximum stability regions
journal, April 1984

Kinnmark, Ingemar P. E.; Gray, William G.
Mathematics and Computers in Simulation, Vol. 26, Issue 2
DOI: 10.1016/0378-4754(84)90039-9

One step integration methods of third-fourth order accuracy with large hyperbolic stability limits
journal, June 1984

Kinnmark, Ingemar P. E.; Gray, William G.
Mathematics and Computers in Simulation, Vol. 26, Issue 3
DOI: 10.1016/0378-4754(84)90056-9

The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable
journal, January 1992

Laprise, René
Monthly Weather Review, Vol. 120, Issue 1
DOI: 10.1175/1520-0493(1992)120<0197:TEEOMW>2.0.CO;2

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

Lin, Shian-Jiann
Monthly Weather Review, Vol. 132, Issue 10
DOI: 10.1175/1520-0493(2004)132<2293:AVLFDC>2.0.CO;2

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

Algorithms for the Assignment and Transportation Problems
journal, March 1957

Munkres, James
Journal of the Society for Industrial and Applied Mathematics, Vol. 5, Issue 1
DOI: 10.1137/0105003

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

An Optimally Stable and Accurate Second-Order SSP Runge-Kutta IMEX Scheme for Atmospheric Applications
journal, January 2018

Rokhzadi, Arman; Mohammadian, Abdolmajid; Charron, Martin
Journal of Advances in Modeling Earth Systems, Vol. 10, Issue 1
DOI: 10.1002/2017MS001065

Vertical discretizations for compressible Euler equation atmospheric models giving optimal representation of normal modes
journal, March 2005

Thuburn, J.; Woollings, T. J.
Journal of Computational Physics, Vol. 203, Issue 2
DOI: 10.1016/j.jcp.2004.08.018

Normal modes of deep atmospheres. I: Spherical geometry
journal, July 2002

Thuburn, John; Wood, Nigel; Staniforth, Andrew
Quarterly Journal of the Royal Meteorological Society, Vol. 128, Issue 584
DOI: 10.1256/003590002320603403

Normal modes of deep atmospheres. II: f–F-plane geometry
journal, July 2002

Thuburn, John; Wood, Nigel; Staniforth, Andrew
Quarterly Journal of the Royal Meteorological Society, Vol. 128, Issue 584
DOI: 10.1256/003590002320603412

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

One step integration methods of third-fourth order accuracy with large hyperbolic stability limits
journal, June 1984

Kinnmark, Ingemar P. E.; Gray, William G.
Mathematics and Computers in Simulation, Vol. 26, Issue 3
DOI: 10.1016/0378-4754(84)90056-9

Algorithms for the Assignment and Transportation Problems
journal, March 1957

Munkres, James
Journal of the Society for Industrial and Applied Mathematics, Vol. 5, Issue 1
DOI: 10.1137/0105003

Similar Records in DOE PAGES and OSTI.GOV collections:

Acceleration of the IMplicit–EXplicit nonhydrostatic unified model of the atmosphere on manycore processors

Journal Article Abdi, Daniel S. ; Giraldo, Francis X. ; Constantinescu, Emil M. ; ... - International Journal of High Performance Computing Applications

We present the acceleration of an IMplicit-EXplicit (IMEX) nonhydrostatic atmospheric model on manycore processors such as graphic processing units (GPUs) and Intel's Many Integrated Core (MIC) architecture. IMEX time integration methods sidestep the constraint imposed by the Courant-Friedrichs-Lewy condition on explicit methods through corrective implicit solves within each time step. In this work, we implement and evaluate the performance of IMEX on manycore processors relative to explicit methods. Using 3D-IMEX at Courant number C = 15, we obtained a speedup of about 4x relative to an explicit time stepping method run with the maximum allowable C = 1. Moreover, themore »« less
Cited by 13
https://doi.org/10.1177/1094342017732395

Full Text Available
thornado-transport: IMEX schemes for two-moment neutrino transport respecting Fermi-Dirac statistics

Conference Chu, Ran ; Endeve, Eirik ; Hauck, Cory ; ...

We develop implicit-explicit (IMEX) schemes for neutrino transport in a background material in the context of a two-moment model that evolves the angular moments of a neutrino phase-space distribution function. Considering the upper and lower bounds that are introduced by Pauli's exclusion principle on the moments, an algebraic moment closure based on Fermi-Dirac statistics and a convex-invariant time integrator both are demanded. A finite-volume/first-order discontinuous Galerkin(DG) method is used to illustrate how an algebraic moment closure based on Fermi-Dirac statistics is needed to satisfy the bounds. Several algebraic closures are compared with these bounds in mind, and the Cernohorsky andmore »« less
https://doi.org/10.1088/1742-6596/1225/1/012013

Full Text Available
On an implicit ENO scheme

Journal Article Wilcoxson, M ; Manousiouthakis, V - Journal of Computational Physics

Essentially non-oscillatory (ENO) schemes, which have high order accuracy in regions where solutions are smooth and effectively capture shocks and discontinuities, have been developed for numerical solution of hyperbolic equations. These methods have traditionally been implemented using explicit Euler or explicit Runge-Kutta time marching schemes and consequently suffer from the Courant-Friedrichs-Lewy (CFL) time step restriction. This restriction is significant for steady-state problems and for transient problems with fronts whose speed is slow compared to the propagation speed in the rest of the domain. In this paper, a novel implicit time marching implementation of the ENO-Roe scheme is developed. Simulation resultsmore »« less
https://doi.org/10.1006/jcph.1994.1203
Efficient IMEX Runge-Kutta Methods for Nonhydrostatic Dynamics

Journal Article Steyer, Andrew ; Vogl, Christopher J. ; Taylor, Mark ; ... - SIAM Journal on Scientific Computing

We develop new implicit-explicit Runge-Kutta (IMEX RK) methods for integrating horizontally explicit, vertically implicit (HEVI) partitionings of nonhydrostatic atmosphere models (HEVI models). These new methods, termed IMKG methods, are IMEX RK methods whose explicit part has optimal or near-optimal stability on the imaginary axis and whose implicit part is I-stable. A specialized stability region is presented for characterizing the stability of IMKG and other IMEX methods integrating HEVI models. Subsequently, we formulate two families of IMEX RK methods to enable deriving IMKG methods with a high explicit stage count for integrating HEVI models with large, stable time-steps. Here, we thenmore »« less
Full Text Available
Implicit–explicit (IMEX) Runge–Kutta methods for non-hydrostatic atmospheric models

Journal Article Gardner, David J. ; Guerra, Jorge E. ; Hamon, François P. ; ... - Geoscientific Model Development (Online)

The efficient simulation of non-hydrostatic atmospheric dynamics requires time integration methods capable of overcoming the explicit stability constraints on time step size arising from acoustic waves. In this work, we investigate various implicit–explicit (IMEX) additive Runge–Kutta (ARK) methods for evolving acoustic waves implicitly to enable larger time step sizes in a global non-hydrostatic atmospheric model. The IMEX formulations considered include horizontally explicit – vertically implicit (HEVI) approaches as well as splittings that treat some horizontal dynamics implicitly. In each case, the impact of solving nonlinear systems in each implicit ARK stage in a linearly implicit fashion is also explored.The accuracymore »« less
Cited by 30
https://doi.org/10.5194/gmd-11-1497-2018

Similar Records

Title: A framework to evaluate IMEX schemes for atmospheric models

Abstract

Citation Formats

CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model journal, November 2011

Implicit–Explicit Multistep Methods for Fast-Wave–Slow-Wave Problems journal, April 2012

Implicit-Explicit Formulations of a Three-Dimensional Nonhydrostatic Unified Model of the Atmosphere (NUMA) journal, January 2013

One step integration methods with maximum stability regions journal, April 1984

One step integration methods of third-fourth order accuracy with large hyperbolic stability limits journal, June 1984

The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable journal, January 1992

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

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

Algorithms for the Assignment and Transportation Problems journal, March 1957

An Overview of the Atmospheric Component of the Energy Exascale Earth System Model journal, August 2019

An Optimally Stable and Accurate Second-Order SSP Runge-Kutta IMEX Scheme for Atmospheric Applications journal, January 2018

Vertical discretizations for compressible Euler equation atmospheric models giving optimal representation of normal modes journal, March 2005

Normal modes of deep atmospheres. I: Spherical geometry journal, July 2002

Normal modes of deep atmospheres. II: f–F-plane geometry journal, July 2002

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

One step integration methods of third-fourth order accuracy with large hyperbolic stability limits journal, June 1984

Algorithms for the Assignment and Transportation Problems journal, March 1957

CAM-SE: A scalable spectral element dynamical core for the Community Atmosphere Model
journal, November 2011

Implicit–Explicit Multistep Methods for Fast-Wave–Slow-Wave Problems
journal, April 2012

Implicit-Explicit Formulations of a Three-Dimensional Nonhydrostatic Unified Model of the Atmosphere (NUMA)
journal, January 2013

One step integration methods with maximum stability regions
journal, April 1984

One step integration methods of third-fourth order accuracy with large hyperbolic stability limits
journal, June 1984

The Euler Equations of Motion with Hydrostatic Pressure as an Independent Variable
journal, January 1992

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

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

Algorithms for the Assignment and Transportation Problems
journal, March 1957

An Overview of the Atmospheric Component of the Energy Exascale Earth System Model
journal, August 2019

An Optimally Stable and Accurate Second-Order SSP Runge-Kutta IMEX Scheme for Atmospheric Applications
journal, January 2018

Vertical discretizations for compressible Euler equation atmospheric models giving optimal representation of normal modes
journal, March 2005

Normal modes of deep atmospheres. I: Spherical geometry
journal, July 2002

Normal modes of deep atmospheres. II: f–F-plane geometry
journal, July 2002

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

One step integration methods of third-fourth order accuracy with large hyperbolic stability limits
journal, June 1984

Algorithms for the Assignment and Transportation Problems
journal, March 1957