Optimal Explicit Runge-Kutta Time Stepping for Density-Based Finite-Volume Solver

2022 ◽  
Author(s):  
Siavash Hedayati Nasab ◽  
Jean-Sébastien Cagnone ◽  
Brian C. Vermeire
2012 ◽  
Vol 140 (4) ◽  
pp. 1257-1284 ◽  
Author(s):  
Paul Ullrich ◽  
Christiane Jablonowski

This paper presents a new approach for discretizing the nonhydrostatic Euler equations in Cartesian geometry using an operator-split time-stepping strategy and unstaggered upwind finite-volume model formulation. Following the method of lines, a spatial discretization of the governing equations leads to a set of coupled nonlinear ordinary differential equations. In general, explicit time-stepping methods cannot be applied directly to these equations because the large aspect ratio between the horizontal and vertical grid spacing leads to a stringent restriction on the time step to maintain numerical stability. Instead, an A-stable linearly implicit Rosenbrock method for evolving the vertical components of the equations coupled to a traditional explicit Runge–Kutta formula in the horizontal is proposed. Up to third-order temporal accuracy is achieved by carefully interleaving the explicit and linearly implicit steps. The time step for the resulting Runge–Kutta–Rosenbrock–type semi-implicit method is then restricted only by the grid spacing and wave speed in the horizontal. The high-order finite-volume model is tested against a series of atmospheric flow problems to verify accuracy and consistency. The results of these tests reveal that this method is accurate, stable, and applicable to a wide range of atmospheric flows and scales.


1991 ◽  
Vol 25 (3) ◽  
pp. 243-263 ◽  
Author(s):  
J. W. van der Burg ◽  
J. G. M. Kuerten ◽  
P. J. Zandbergen

1993 ◽  
Vol 115 (4) ◽  
pp. 781-790 ◽  
Author(s):  
G. A. Gerolymos

In the present work an algorithm for the numerical integration of the three-dimensional unsteady Euler equations in vibrating transonic compressor cascades is described. The equations are discretized in finite-volume formulation in a mobile grid using isoparametric brick elements. They are integrated in time using Runge-Kutta schemes. A thorough discussion of the boundary conditions used and of their influence on results is undertaken. The influence of grid refinement on computational results is examined. Unsteady convergence of results is discussed.


2019 ◽  
Vol 12 (2) ◽  
pp. 651-676 ◽  
Author(s):  
Christian Kühnlein ◽  
Willem Deconinck ◽  
Rupert Klein ◽  
Sylvie Malardel ◽  
Zbigniew P. Piotrowski ◽  
...  

Abstract. We present a nonhydrostatic finite-volume global atmospheric model formulation for numerical weather prediction with the Integrated Forecasting System (IFS) at ECMWF and compare it to the established operational spectral-transform formulation. The novel Finite-Volume Module of the IFS (henceforth IFS-FVM) integrates the fully compressible equations using semi-implicit time stepping and non-oscillatory forward-in-time (NFT) Eulerian advection, whereas the spectral-transform IFS solves the hydrostatic primitive equations (optionally the fully compressible equations) using a semi-implicit semi-Lagrangian scheme. The IFS-FVM complements the spectral-transform counterpart by means of the finite-volume discretization with a local low-volume communication footprint, fully conservative and monotone advective transport, all-scale deep-atmosphere fully compressible equations in a generalized height-based vertical coordinate, and flexible horizontal meshes. Nevertheless, both the finite-volume and spectral-transform formulations can share the same quasi-uniform horizontal grid with co-located arrangement of variables, geospherical longitude–latitude coordinates, and physics parameterizations, thereby facilitating their comparison, coexistence, and combination in the IFS. We highlight the advanced semi-implicit NFT finite-volume integration of the fully compressible equations of IFS-FVM considering comprehensive moist-precipitating dynamics with coupling to the IFS cloud parameterization by means of a generic interface. These developments – including a new horizontal–vertical split NFT MPDATA advective transport scheme, variable time stepping, effective preconditioning of the elliptic Helmholtz solver in the semi-implicit scheme, and a computationally efficient implementation of the median-dual finite-volume approach – provide a basis for the efficacy of IFS-FVM and its application in global numerical weather prediction. Here, numerical experiments focus on relevant dry and moist-precipitating baroclinic instability at various resolutions. We show that the presented semi-implicit NFT finite-volume integration scheme on co-located meshes of IFS-FVM can provide highly competitive solution quality and computational performance to the proven semi-implicit semi-Lagrangian integration scheme of the spectral-transform IFS.


Sign in / Sign up

Export Citation Format

Share Document