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.

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.

scholarly journals Stability analysis of the implicit finite-difference-based upwind lattice Boltzmann schemes

2019 ◽  
pp. 116-127
Author(s):  
Г.В. Кривовичев ◽  
М.П. Мащинская
Keyword(s):  
Stability Analysis ◽  
Finite Difference ◽  
Lattice Boltzmann ◽  
Kinetic Equations ◽  
Finite Difference Schemes ◽  
Courant Number ◽  
Von Neumann ◽  
Upwind Schemes ◽  
Implicit Finite Difference ◽  
The Stability

Статья посвящена анализу устойчивости неявных конечно-разностных схем для системы кинетических уравнений, применяемых для проведения гидродинамических расчетов в рамках метода решеточных уравнений Больцмана. Представлены семейства двухслойных и трехслойных схем с направленными разностями первого-четвертого порядков аппроксимации по пространственным переменным. Важной особенностью схем является то, что конвективные слагаемые аппроксимируются одной конечной разностью. Показано, что в выражении для аппроксимационной вязкости схем высоких порядков отсутствуют фиктивные слагаемые, что позволяет применять их во всем диапазоне значений времени релаксации. Анализ устойчивости проводится по линейному приближению с использованием метода Неймана. Получены приближенные условия устойчивости в виде неравенств на значения параметра Куранта. При расчетах показано, что площади областей устойчивости в пространстве параметров у двухслойных схем больше, чем у трехслойных. Исследованные схемы могут применяться при расчетах как непосредственно, так и в методах типа предиктор-корректор. The paper is devoted to the stability analysis of the implicit finite-difference schemes for the system of kinetic equations used for the hydrodynamic computations in the framework of the lattice Boltzmann method. The families of two- and three-layer upwind schemes of the first to fourth approximation orders on spatial variables are considered. An important feature of the presented schemes is that the convective terms are approximated by one finite difference. It is shown that, for the high-order schemes, in the expression for the current viscosity there are no fictitious terms, which makes it possible to perform computations in the whole range of relaxation time values. The stability analysis is based on the application of the von Neumann method to the linear approximations of the schemes. The stability conditions are obtained in the form of inequalities imposed on the Courant number values. It is also shown that the areas of stability domains for the two-layer schemes are greater than for the three-layer schemes in the parameter space. The considered schemes can be used as the fully implicit schemes in computational algorithms directly or in the predictor-corrector methods.

Transport Schemes in GungHo

2021 ◽  
Author(s):  
James Kent
Keyword(s):  
Finite Element ◽  
Finite Volume ◽  
Mixed Finite Element ◽  
Method Of Lines ◽  
Finite Volume Methods ◽  
Kutta Method ◽  
Finite Volume Schemes ◽  
Dynamical Core ◽  
Runge Kutta Method ◽  
The Stability

<p>GungHo is the mixed finite-element dynamical core under development by the Met Office. A key component of the dynamical core is the transport scheme, which advects density, temperature, moisture, and the winds, throughout the atmosphere. Transport in GungHo is performed by finite-volume methods, to ensure conservation of certain quantaties. There are a range of different finite-volume schemes being considered for transport, including the Runge-Kutta/method-of-lines and COSMIC/Lin-Rood schemes. Additional horizontal/vertical splitting approaches are also under consideration, to improve the stability aspects of the model. Here we discuss these transport options and present results from the GungHo framework, featuring both prescribed velocity advection tests and full dry dynamical core tests. </p>

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.

scholarly journals Linear stability of confined flow around a 180-degree sharp bend

2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard
Keyword(s):  
Reynolds Number ◽  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Critical Reynolds Number ◽  
Three Dimensional ◽  
Base Flow ◽  
Two Dimensional ◽  
Two Dimensional Flow ◽  
The Stability

This study seeks to characterise the breakdown of the steady two-dimensional solution in the flow around a 180-degree sharp bend to infinitesimal three-dimensional disturbances using a linear stability analysis. The stability analysis predicts that three-dimensional transition is via a synchronous instability of the steady flows. A highly accurate global linear stability analysis of the flow was conducted with Reynolds number $\mathit{Re}<1150$ and bend opening ratio (ratio of bend width to inlet height) $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 5$. This range of $\mathit{Re}$ and $\unicode[STIX]{x1D6FD}$ captures both steady-state two-dimensional flow solutions and the inception of unsteady two-dimensional flow. For $0.2\leqslant \unicode[STIX]{x1D6FD}\leqslant 1$, the two-dimensional base flow transitions from steady to unsteady at higher Reynolds number as $\unicode[STIX]{x1D6FD}$ increases. The stability analysis shows that at the onset of instability, the base flow becomes three-dimensionally unstable in two different modes, namely a spanwise oscillating mode for $\unicode[STIX]{x1D6FD}=0.2$ and a spanwise synchronous mode for $\unicode[STIX]{x1D6FD}\geqslant 0.3$. The critical Reynolds number and the spanwise wavelength of perturbations increase as $\unicode[STIX]{x1D6FD}$ increases. For $1<\unicode[STIX]{x1D6FD}\leqslant 2$ both the critical Reynolds number for onset of unsteadiness and the spanwise wavelength decrease as $\unicode[STIX]{x1D6FD}$ increases. Finally, for $2<\unicode[STIX]{x1D6FD}\leqslant 5$, the critical Reynolds number and spanwise wavelength remain almost constant. The linear stability analysis also shows that the base flow becomes unstable to different three-dimensional modes depending on the opening ratio. The modes are found to be localised near the reattachment point of the first recirculation bubble.

Logically rectangular finite volume methods with adaptive refinement on the sphere

2009 ◽  
Vol 367 (1907) ◽  
pp. 4483-4496 ◽  
Author(s):  
Marsha J. Berger ◽  
Donna A. Calhoun ◽  
Christiane Helzel ◽  
Randall J. LeVeque
Keyword(s):  
Finite Volume ◽  
Shallow Water Equations ◽  
Adaptive Mesh Refinement ◽  
Three Dimensional ◽  
Adaptive Mesh ◽  
Finite Volume Methods ◽  
Adaptive Refinement ◽  
Two Dimensional ◽  
Equilibrium Solutions ◽  
Courant Number

The logically rectangular finite volume grids for two-dimensional partial differential equations on a sphere and for three-dimensional problems in a spherical shell introduced recently have nearly uniform cell size, avoiding severe Courant number restrictions. We present recent results with adaptive mesh refinement using the G eo C law software and demonstrate well-balanced methods that exactly maintain equilibrium solutions, such as shallow water equations for an ocean at rest over arbitrary bathymetry.

Sign in / Sign up

Export Citation Format

Share Document