The stability of stationary solution for outflow problem on the navier-stokes-poisson system

2016 ◽  
Vol 36 (4) ◽  
pp. 1098-1116 ◽  
Author(s):  
Mina JIANG ◽  
Suhua LAI ◽  
Haiyan YIN ◽  
Changjiang ZHU
2013 ◽  
Vol 136 (2) ◽  
Author(s):  
Graham Ashcroft ◽  
Christian Frey ◽  
Kathrin Heitkamp ◽  
Christian Weckmüller

This is the first part of a series of two papers on unsteady computational fluid dynamics (CFD) methods for the numerical simulation of aerodynamic noise generation and propagation. In this part, the stability, accuracy, and efficiency of implicit Runge–Kutta schemes for the temporal integration of the compressible Navier–Stokes equations are investigated in the context of a CFD code for turbomachinery applications. Using two model academic problems, the properties of two explicit first stage, singly diagonally implicit Runge–Kutta (ESDIRK) schemes of second- and third-order accuracy are quantified and compared with more conventional second-order multistep methods. Finally, to assess the ESDIRK schemes in the context of an industrially relevant configuration, the schemes are applied to predict the tonal noise generation and transmission in a modern high bypass ratio fan stage and comparisons with the corresponding experimental data are provided.


1976 ◽  
Vol 78 (2) ◽  
pp. 355-383 ◽  
Author(s):  
H. Fasel

The stability of incompressible boundary-layer flows on a semi-infinite flat plate and the growth of disturbances in such flows are investigated by numerical integration of the complete Navier–;Stokes equations for laminar two-dimensional flows. Forced time-dependent disturbances are introduced into the flow field and the reaction of the flow to such disturbances is studied by directly solving the Navier–Stokes equations using a finite-difference method. An implicit finitedifference scheme was developed for the calculation of the extremely unsteady flow fields which arose from the forced time-dependent disturbances. The problem of the numerical stability of the method called for special attention in order to avoid possible distortions of the results caused by the interaction of unstable numerical oscillations with physically meaningful perturbations. A demonstration of the suitability of the numerical method for the investigation of stability and the initial growth of disturbances is presented for small periodic perturbations. For this particular case the numerical results can be compared with linear stability theory and experimental measurements. In this paper a number of numerical calculations for small periodic disturbances are discussed in detail. The results are generally in fairly close agreement with linear stability theory or experimental measurements.


2013 ◽  
Vol 2013 ◽  
pp. 1-6
Author(s):  
Wen-Juan Wang ◽  
Yan Jia

We study the stability issue of the generalized 3D Navier-Stokes equations. It is shown that if the weak solutionuof the Navier-Stokes equations lies in the regular class∇u∈Lp(0,∞;Bq,∞0(ℝ3)),(2α/p)+(3/q)=2α,2<q<∞,0<α<1, then every weak solutionv(x,t)of the perturbed system converges asymptotically tou(x,t)asvt-utL2→0,t→∞.


Author(s):  
Fayçal Hammami ◽  
Nader Ben Cheikh ◽  
Brahim Ben Beya

This paper deals with the numerical study of bifurcations in a two-sided lid driven cavity flow. The flow is generated by moving the upper wall to the right while moving the left wall downwards. Numerical simulations are performed by solving the unsteady two dimensional Navier-Stokes equations using the finite volume method and multigrid acceleration. In this problem, the ratio of the height to the width of the cavity are ranged from H/L = 0.25 to 1.5. The code for this cavity is presented using rectangular cavity with the grids 144 × 36, 144 × 72, 144 × 104, 144 × 136, 144 × 176 and 144 × 216. Numerous comparisons with the results available in the literature are given. Very good agreements are found between current numerical results and published numerical results. Various velocity ratios ranged in 0.01≤ α ≤ 0.99 at a fixed aspect ratios (A = 0.5, 0.75, 1.25 and 1.5) were considered. It is observed that the transition to the unsteady regime follows the classical scheme of a Hopf bifurcation. The stability analysis depending on the aspect ratio, velocity ratios α and the Reynolds number when transition phenomenon occurs is considered in this paper.


2011 ◽  
Vol 1 (3) ◽  
pp. 215-234 ◽  
Author(s):  
Lizhen Chen ◽  
Jie Shen ◽  
Chuanju Xu

AbstractWe propose and analyze spectral direction splitting schemes for the incompressible Navier-Stokes equations. The schemes combine a Legendre-spectral method for the spatial discretization and a pressure-stabilization/direction splitting scheme for the temporal discretization, leading to a sequence of one-dimensional elliptic equations at each time step while preserving the same order of accuracy as the usual pressure-stabilization schemes. We prove that these schemes are unconditionally stable, and present numerical results which demonstrate the stability, accuracy, and efficiency of the proposed methods.


Author(s):  
Sakir Amiroudine

The case of a supercritical fluid heated from below (Rayleigh-Bénard) in a rectangular cavity is first presented. The stability of the two boundary layers (hot and cold) is analyzed by numerically solving the Navier-Stokes equations with a van der Waals gas and stability diagrams are derived. The very large compressibility and the very low heat diffusivity of near critical pure fluids induce very large density gradients which lead to a Rayleigh–Taylor-like gravitational instability of the heat diffusion layer and results in terms of growth rates and wave numbers are presented. Depending on the relative direction of the interface or the boundary layer with respect to vibration, vibrational forces can destabilize a thermal boundary layer, resulting in parametric/Rayleigh vibrational instabilities. This has recently been achieved by using a numerical model which does not require any equation of state and directly calculates properties from NIST data base, for instance.


Author(s):  
Jean-Yves Chemin ◽  
Benoit Desjardins ◽  
Isabelle Gallagher ◽  
Emmanuel Grenier

Let us now detail the stability properties of an Ekman layer introduced in Part I, page 11. First we will recall how to compute the critical Reynolds number. Then we will describe briefly what happens at larger Reynolds numbers. The first step in the study of the stability of the Ekman layer is to consider the linear stability of a pure Ekman spiral of the form where U∞ is the velocity away from the layer and ζ is the rescaled vertical component ζ = x3/√εν. The corresponding Reynolds number is Let us consider the Navier–Stokes–Coriolis equations, linearized around uE The problem is now to study the (linear) stability of the 0 solution of the system (LNSCε). If u=0 is stable we say that uE is linearly stable, if not we say that it is linearly unstable. Numerical results show that u=0 is stable if and only if Re<Rec where Rec can be evaluated numerically. Up to now there is no mathematical proof of this fact, and it is only possible to prove that 0 is linearly stable for Re<Re1 and unstable for Re>Re2 with Re1<Rec<Re2, Re1 being obtained by energy estimates and Re2 by a perturbative analysis of the case Re=∞. We would like to emphasize that the numerical results are very reliable and can be considered as definitive results, since as we will see below, the stability analysis can be reduced to the study of a system of ordinary differential equations posed on the half-space, with boundary conditions on both ends, a system which can be studied arbitrarily precisely, even on desktop computers (first computations were done in the 1960s by Lilly).


Sign in / Sign up

Export Citation Format

Share Document