Abstract:We conduct a thorough study of different forms of horizontally explicit and vertically implicit (HEVI) time-integration strategies for the compressible Euler equations on spherical domains typical of nonhydrostatic global atmospheric applications. We compare the computational time and complexity of two nonlinear variants (NHEVI-GMRES and NHEVI-LU) and a linear variant (LHEVI). We report on the performance of these three variants for a number of additive Runge-Kutta Methods ranging in order of accuracy from second through fifth, and confirm the expected order of accuracy of the HEVI methods for each time-integrator. To gauge the maximum usable time-step of each HEVI method, we run simulations of a nonhydrostatic baroclinic instability for 100 days and then use this time-step to compare the time-to-solution of each method. The results show that NHEVI-LU is 2x faster than NHEVI-GMRES, and LHEVI is 5x faster than NHEVI-LU, for the idealized cases tested. The baroclinic instability and inertia-gravity wave simulations indicate that the optimal choice of time-integrator is LHEVI with either second or third order schemes, as both schemes yield similar time to solution and relative L2 error at their maximum usable time-steps. In the future, we will report on whether these results hold for more complex problems using, e.g., real atmospheric data and/or a higher model top typical of space weather applications.

What problem does this paper attempt to address?

This paper attempts to solve the problems of efficiency and stability of time - integration methods in non - hydrostatic atmospheric models. Specifically, the research focuses on the performance of the horizontally explicit - vertically implicit (HEVI) time - integration strategy when solving the compressible Euler equations. The following are the main problems and objectives of this paper: 1. **Improve computational efficiency**: Traditional explicit time - integration methods are very computationally expensive because they require very small time steps to maintain numerical stability. By adopting the HEVI method, larger time steps can be used while maintaining stability, thereby significantly improving computational efficiency. 2. **Optimize time - integration methods**: The paper compares different types of HEVI methods, including two nonlinear variants (NHEVI - GMRES and NHEVI - LU) and one linear variant (LHEVI). Through these comparisons, the researchers hope to find the optimal time - integration scheme to balance computational speed and accuracy. 3. **Verify the applicability of the method**: To evaluate the actual performance of these HEVI methods, the researchers conducted multiple benchmark tests, including simulating non - hydrostatic baroclinic instability and inertial - gravity waves. Through these tests, the maximum available time step for each method can be determined, and their solution times can be compared. 4. **Provide references for future research**: The research also explores whether these results are applicable to more complex problems, such as space - weather applications using real - atmospheric data or higher - model vertices. This provides directions for further research. ### Specific problems of the research - **How to select an appropriate time - integration method**: By comparing the performance of different HEVI methods, the researchers hope to find the time - integration method most suitable for non - hydrostatic atmospheric models. - **How to improve computational efficiency**: The researchers hope to improve the existing HEVI methods so that they can improve computational efficiency without reducing stability and accuracy. - **Verify the stability and accuracy of the method**: Through long - term numerical simulations (such as a 100 - day non - hydrostatic baroclinic instability simulation), verify the stability and accuracy of different HEVI methods. ### Main conclusions According to the results of the paper, the LHEVI method performs best in idealized test cases, especially when using second - order or third - order schemes. It not only has similar solution times and relative L2 errors, but is also faster than other methods. Specifically: - NHEVI - LU is 2 times faster than NHEVI - GMRES. - LHEVI is 5 times faster than NHEVI - LU. These results indicate that for non - hydrostatic atmospheric models, the LHEVI method may be the optimal choice. Future work will further verify whether these results are applicable to more complex practical problems.

A Performance Study of Horizontally Explicit Vertically Implicit (HEVI) Time-Integrators for Non-Hydrostatic Atmospheric Models

An energetically balanced, quasi-Newton integrator for non-hydrostatic vertical atmospheric dynamics

A Comparison of Rosenbrock-Wanner and Crank-Nicolson Time Integrators for Atmospheric Modelling

Applicability Study of Euler–Lagrange Integration Scheme in Constructing Small-Scale Atmospheric Dynamics Models

A semi-implicit compressible model for atmospheric flows with seamless access to soundproof and hydrostatic dynamics

Parallel-in-time integration of the shallow water equations on the rotating sphere using Parareal and MGRIT

A comparative computational study of different formulations of the compressible Euler equations for mesoscale atmospheric flows in a finite volume framework

A Global Non-Hydrostatic Atmospheric Model with a Mass and Energy Conserving Vertically-Implicit-Correction (VIC) Scheme

High-Order Dynamic Integration Method (HODIM) for Modeling Turbulent Fluid Dynamics

Comparison of the performance of high-order schemes based on the gas-kinetic and HLLC fluxes

Finite-element schemes for extended integrations of atmospheric models

Exponential time integration for 3D compressible atmospheric models

Filter stabilization for the mildly compressible Euler equations with application to atmosphere dynamics simulations

Current and Emerging Time-Integration Strategies in Global Numerical Weather and Climate Prediction

Adaptively Implicit Advection for Atmospheric Flows

On the performances of standard and kinetic energy preserving time-integration methods for incompressible-flow simulations

Implicit high order discontinuous Galerkin method for compressible laminar and turbulent flow simulation

Stabilized bases for high-order, interpolation semi-Lagrangian, element-based tracer transport

Massively parallel solvers for elliptic partial differential equations in numerical weather and climate prediction

A new family of semi-implicit Finite Volume / Virtual Element methods for incompressible flows on unstructured meshes