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

2016 ◽  
Vol 138 (10) ◽  
Author(s):  
Santosh Konangi ◽  
Nikhil K. Palakurthi ◽  
Urmila Ghia
Keyword(s):  
Stability Analysis ◽  
Euler Equations ◽  
Two Dimensional ◽  
Von Neumann ◽  
One Dimensional ◽  
Flow Equations ◽  
Solution Scheme ◽  
Von Neumann Stability ◽  
The Stability ◽  
Stability Limits

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.

A modified two-dimensional triangular lattice model under honk environment

2020 ◽  
Vol 31 (06) ◽  
pp. 2050089
Author(s):  
Cong Zhai ◽  
Weitiao Wu
Keyword(s):  
Stability Analysis ◽  
Hydrodynamic Model ◽  
High Dimensional ◽  
Two Dimensional ◽  
Lattice Hydrodynamic Model ◽  
Dimensional Lattice ◽  
One Dimensional ◽  
The Stability ◽  
Honk Effect ◽  
Real Traffic

The honk effect is not uncommon in the real traffic and may exert great influence on the stability of traffic flow. As opposed to the linear description of the traditional one-dimensional lattice hydrodynamic model, the high-dimensional lattice hydrodynamic model is a gridded analysis of the real traffic environment, which is a generalized form of the one-dimensional lattice model. Meanwhile, the high-dimensional traffic flow exposed to the open-ended environment is more likely to be affected by the honk effect. In this paper, we propose an extension of two-dimensional triangular lattice hydrodynamic model under honk environment. The stability condition is obtained via the linear stability analysis, which shows that the stability region in the phase diagram can be effectively enlarged under the honk effect. Modified Korteweg–de Vries equations are derived through the nonlinear stability analysis method. The kink–antikink solitary wave solution is obtained by solving the equation, which can be used to describe the propagation characteristics of density waves near the critical point. Finally, the simulation example verifies the correctness of the above theoretical analysis.

scholarly journals A Stability Analysis of Finite-Volume Advection Schemes Permitting Long Time Steps

2007 ◽  
Vol 135 (7) ◽  
pp. 2658-2673 ◽  
Author(s):  
Peter Hjort Lauritzen
Keyword(s):  
Stability Analysis ◽  
Finite Volume ◽  
Two Dimensional ◽  
Finite Volume Schemes ◽  
Courant Number ◽  
Von Neumann ◽  
Dispersion Properties ◽  
Advection Schemes ◽  
Long Time ◽  
The Stability

Abstract Finite-volume schemes developed in the meteorological community that permit long time steps are considered. These include Eulerian flux-form schemes as well as fully two-dimensional and cascade cell-integrated semi-Lagrangian (CISL) schemes. A one- and two-dimensional Von Neumann stability analysis of these finite-volume advection schemes is given. Contrary to previous analysis, no simplifications in terms of reducing the formal order of the schemes, which makes the analysis mathematically less complex, have been applied. An interscheme comparison of both dissipation and dispersion properties is given. The main finding is that the dissipation and dispersion properties of Eulerian flux-form schemes are sensitive to the choice of inner and outer operators applied in the scheme that can lead to increased numerical damping for large Courant numbers. This spurious dependence on the integer value of the Courant number disappears if the inner and outer operators are identical, in which case, under the assumptions used in the stability analysis, the Eulerian flux-form scheme becomes identical to the cascade scheme. To explain these properties a conceptual interpretation of the flux-based Eulerian schemes is provided. Of the two CISL schemes, the cascade scheme has superior stability properties.

One-dimensional model for the Rijke tube

1989 ◽  
Vol 202 ◽  
pp. 83-96 ◽  
Author(s):  
C. Nicoli ◽  
P. Pelcé
Keyword(s):  
Heat Transfer ◽  
Transfer Function ◽  
Dimensional Model ◽  
Unsteady Heat Transfer ◽  
Order Unity ◽  
Rijke Tube ◽  
Thermoacoustic Instability ◽  
One Dimensional ◽  
The Stability ◽  
Stability Limits

We develop a simple model in which longitudinal, compressible, unsteady heat transfer between heater and gas is computed in the small-Mach-number limit. This calculation is used to determine the transfer function of the heater, which plays an important role in the stability limits of the thermoacoustic instability of the Rijke tube. The transfer function is determined analytically in the limit of small expansion parameter γ, and numerically for γ of order unity. In the case ρμ/cp = constant, an analytical solution can be found.

Chaos in the Rulkov Neuron Model Based on Marotto’s Theorem

2021 ◽  
Vol 31 (15) ◽  
Author(s):  
Penghe Ge ◽  
Hongjun Cao
Keyword(s):  
Neuron Model ◽  
Initial Conditions ◽  
Stability Conditions ◽  
Two Dimensional ◽  
Bifurcation Diagrams ◽  
Fast Subsystem ◽  
One Dimensional ◽  
Sensitive Dependence ◽  
The Stability ◽  
Slow Subsystem

The existence of chaos in the Rulkov neuron model is proved based on Marotto’s theorem. Firstly, the stability conditions of the model are briefly renewed through analyzing the eigenvalues of the model, which are very important preconditions for the existence of a snap-back repeller. Secondly, the Rulkov neuron model is decomposed to a one-dimensional fast subsystem and a one-dimensional slow subsystem by the fast–slow dynamics technique, in which the fast subsystem has sensitive dependence on the initial conditions and its snap-back repeller and chaos can be verified by numerical methods, such as waveforms, Lyapunov exponents, and bifurcation diagrams. Thirdly, for the two-dimensional Rulkov neuron model, it is proved that there exists a snap-back repeller under two iterations by illustrating the existence of an intersection of three surfaces, which pave a new way to identify the existence of a snap-back repeller.

Three-dimensional stability analysis for a salt-finger convecting layer

2018 ◽  
Vol 841 ◽  
pp. 636-653
Author(s):  
Ting-Yueh Chang ◽  
Falin Chen ◽  
Min-Hsing Chang
Keyword(s):  
Stability Analysis ◽  
Three Dimensional ◽  
Longitudinal Mode ◽  
Transverse Mode ◽  
Solute Diffusion ◽  
Two Dimensional ◽  
Grashof Number ◽  
Significant Difference ◽  
The Stability ◽  
Mode A

A three-dimensional linear stability analysis is carried out for a convecting layer in which both the temperature and solute distributions are linear in the horizontal direction. The three-dimensional results show that, for $Le=3$ and 100, the most unstable mode occurs invariably as the longitudinal mode, a vortex roll with its axis perpendicular to the longitudinal plane, suggesting that the two-dimensional results are sufficient to illustrate the stability characteristics of the convecting layer. Two-dimensional results show that the stability boundaries of the transverse mode (a vortex roll with its axis perpendicular to the transverse plane) and the longitudinal modes are virtually overlapped in the regime dominated by thermal diffusion and the regime dominated by solute diffusion, while these two modes hold a significant difference in the regime the salt-finger instability prevails. More precisely, the instability area in terms of thermal Grashof number $Gr$ and solute Grashof number $Gs$ is larger for the longitudinal mode than the transverse mode, implying that, under any circumstance, the longitudinal mode is always more unstable than the transverse mode.

Sign in / Sign up

Export Citation Format

Share Document