Linear Stability Analysis of Runge–Kutta-Based Partial Time-Splitting Schemes for the Euler Equations

Abstract 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 div...

[1]  F. Bretherton.  The propagation of groups of internal gravity waves in a shear flow , 1966 .

[2]  Ed Anderson,et al.  LAPACK Users' Guide , 1995 .

[3]  J. Klemp,et al.  The Simulation of Three-Dimensional Convective Storm Dynamics , 1978 .

[4]  Jimy Dudhia,et al.  Conservative Split-Explicit Time Integration Methods for the Compressible Nonhydrostatic Equations , 2007 .

[5]  D. Durran.  Numerical methods for wave equations in geophysical fluid dynamics , 1999 .

[6]  Louis J. Wicker,et al.  A Time-Splitting Scheme for the Elastic Equations Incorporating Second-Order Runge–Kutta Time Differencing , 1998 .

[7]  Chi-Wang Shu.  Total-variation-diminishing time discretizations , 1988 .

[8]  C. Epifanio,et al.  An Analysis of Klemp–Wilhelmson Schemes as Applied to Large-Scale Wave Modes , 2008 .

[9]  Louis J. Wicker,et al.  Time-Splitting Methods for Elastic Models Using Forward Time Schemes , 2002 .

[10]  Rüdiger Weiner,et al.  Explicit Two-Step Peer Methods for the Compressible Euler Equations , 2009 .

[11]  W. Skamarock,et al.  The stability of time-split numerical methods for the hydrostatic and the nonhydrostatic elastic equations , 1992 .

[12]  Randall J. LeVeque,et al.  Numerical methods based on additive splittings for hyperbolic partial differential equations , 1981 .

[13]  Almut Gassmann,et al.  A Consistent Time-Split Numerical Scheme Applied to the Nonhydrostatic Compressible Equations* , 2007 .

[14]  Michael Baldauf.  Stability analysis for linear discretisations of the advection equation with Runge-Kutta time integration , 2008, J. Comput. Phys..

[15]  William C. Skamarock,et al.  A time-split nonhydrostatic atmospheric model for weather research and forecasting applications , 2008, J. Comput. Phys..

[16]  Jack Dongarra,et al.  LAPACK Users' Guide, 3rd ed. , 1999 .

[17]  J. Steppeler,et al.  Meso-gamma scale forecasts using the nonhydrostatic model LM , 2003 .

[18]  James Demmel,et al.  LAPACK Users' Guide, Third Edition , 1999, Software, Environments and Tools.

[19]  Steven J. Ruuth,et al.  High-Order Strong-Stability-Preserving Runge-Kutta Methods with Downwind-Biased Spatial Discretizations , 2004, SIAM J. Numer. Anal..