von Neumann Stability Analysis of a Segregated Pressure-Based Solution Scheme for One-Dimensional and Two-Dimensional Flow Equations

Abstract

The goal of this paper is to derive the von Neumann stability conditions for the pressure-based solution scheme, semi-implicit method for pressure-linked equations (SIMPLE). The SIMPLE scheme lies at the heart of a class of computational fluid dynamics (CFD) algorithms built into several commercial and open-source CFD software packages. To the best of the authors' knowledge, no readily usable stability guidelines appear to be available for this popularly employed scheme. The Euler equations are examined, as the inclusion of viscosity in the Navier–Stokes (NS) equation serves to only soften the stability limits. First, the one-dimensional (1D) Euler equations are studied, and their stability properties are delineated. Next, a rigorous stability analysis is carried out for the two-dimensional (2D) Euler equations; the analysis of the 2D equations is considerably more challenging as compared to analysis of the 1D form of equations. The Euler equations are discretized using finite differences on a staggered grid, which is used to achieve equivalence to finite-volume discretization. Error amplification matrices are determined from the stability analysis, stable and unstable regimes are identified, and practical stability limits are predicted in terms of the maximum allowable Courant–Friedrichs–Lewy (CFL) number as a function of Mach number. The predictions are verified using the Riemann problem, and very good agreement is obtained between the analytically predicted and the “experimentally” observed CFL values. The successfully tested stability limits are presented in graphical form, as compared to complicated mathematical expressions often reported in published literature. Since our analysis accounts for the solution scheme along with the full system of flow equations, the conditions reported in this paper offer practical value over the conditions that arise from analysis of simplified 1D model equations.

Kama Şekilli Kanat Üzerindeki Süpersonik Akışı Çözmek için Sayısal Bir Algoritma

“…Through Von Neumann stability analysis, it can be seen that for high frequencies, this scheme is non-dissipative (Leveque, 1992;Strikwerda, 2004;von Neuman and Richtmyer, 1950;Konangi et al, 2016). As a consequence large oscillations will be apparent in regions around a shock wave, and artificial dissipation will have to be added to the scheme (Degani and Fox, 1994;Cavus, 2013).…”

“…As a consequence large oscillations will be apparent in regions around a shock wave, and artificial dissipation will have to be added to the scheme (Degani and Fox, 1994;Cavus, 2013). Additionally, the implicit scheme will require a solution to a tridiagonal matrix (Strikwerda, 2004;Laney, 1998;von Neuman and Richtmyer, 1950;Konangi et al, 2016).…”