Open AccessLinear Stability Analysis of Runge–Kutta-Based Partial Time-Splitting Schemes for the Euler Equations
Author(s) -
Michael Baldauf
Publication year - 2010
Publication title -
monthly weather review
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.862
H-Index - 179
eISSN - 1520-0493
pISSN - 0027-0644
DOI - 10.1175/2010mwr3355.1
Subject(s) - euler equations , runge–kutta methods , backward euler method , stability (learning theory) , discretization , total variation diminishing , mathematics , euler's formula , truncation error , courant–friedrichs–lewy condition , divergence (linguistics) , eulerian path , semi implicit euler method , euler method , numerical analysis , mathematical analysis , computer science , machine learning , linguistics , philosophy , lagrangian
For atmospheric simulation models with resolutions from about 10 km to the subkilometer cloud-resolving scale, the complete nonhydrostatic compressible Euler equations are often used. An important integration technique for them is the time-splitting (or split explicit) method. This article presents a comprehensive numerical stability analysis of Runge–Kutta (RK)-based partial time-splitting schemes. To this purpose a linearized two-dimensional (2D) compressible Euler system containing advection (as the slow process), sound, and gravity wave terms (as fast processes) is considered. These processes are the most important ones in limiting stability. First, the detailed stability properties are discussed with regard to several off-centering weights for each fast process described by horizontally explicit, vertically implicit schemes. Then the stability properties of the temporally and spatially discretized three-stage RK scheme for the complete 2D Euler equations and their stabilization (e.g., by divergence damping) are discussed. The main goal is to find optimal values for all of the occurring numerical parameters to guarantee stability in operational model applications. Furthermore, formal orders of temporal truncation errors for the time-splitting schemes are calculated. With the same methodology, two alternatives to the three-stage RK method, a so-called RK3-TVD method, and a new four-stage, second-order RK method are inspected.