Linear stability analysis of cylindrical Rayleigh–Bénard convection

2012 ◽  
Vol 711 ◽  
pp. 27-39 ◽  
Author(s):  
Bo-Fu Wang ◽  
Dong-Jun Ma ◽  
Cheng Chen ◽  
De-Jun Sun
Keyword(s):  
Stability Analysis ◽  
Prandtl Number ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Axisymmetric Flow ◽  
Base Flow ◽  
Prandtl Numbers ◽  
The Stability ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

AbstractThe instabilities and transitions of flow in a vertical cylindrical cavity with heated bottom, cooled top and insulated sidewall are investigated by linear stability analysis. The stability boundaries for the axisymmetric flow are derived for Prandtl numbers from 0.02 to 1, for aspect ratio $A$ ($A= H/ R= \mathrm{height} / \mathrm{radius} $) equal to 1, 0.9, 0.8, 0.7, respectively. We found that there still exists stable non-trivial axisymmetric flow beyond the second bifurcation in certain ranges of Prandtl number for $A= 1$, $0. 9$ and 0.8, excluding the $A= 0. 7$ case. The finding for $A= 0. 7$ is that very frequent changes of critical mode (azimuthal Fourier mode) of the second bifurcation occur when the Prandtl number is changed, where five kinds of steady modes $m= 1, 2, 8, 9, 10$ and three kinds of oscillatory modes $m= 3, 4, 6$ are presented. These multiple modes indicate different flow structures triggered at the transitions. The instability mechanism of the flow is explained by kinetic energy transfer analysis, which shows that the radial or axial shear of base flow combined with buoyancy mechanism leads to the instability results.

Rayleigh–Bénard instability generated by a diffusion flame

2013 ◽  
Vol 736 ◽  
pp. 464-494 ◽  
Author(s):  
P. Pearce ◽  
J. Daou
Keyword(s):  
Stability Analysis ◽  
Rayleigh Number ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Diffusion Flame ◽  
Damkohler Number ◽  
Benard Convection ◽  
The Stability ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

AbstractWe investigate the Rayleigh–Bénard convection problem within the context of a diffusion flame formed in a horizontal channel where the fuel and oxidizer concentrations are prescribed at the porous walls. This problem seems to have received no attention in the literature. When formulated in the low-Mach-number approximation the model depends on two main non-dimensional parameters, the Rayleigh number and the Damköhler number, which govern gravitational strength and reaction speed respectively. In the steady state the system admits a planar diffusion flame solution; the aim is to find the critical Rayleigh number at which this solution becomes unstable to infinitesimal perturbations. In the Boussinesq approximation, a linear stability analysis reduces the system to a matrix equation with a solution comparable to that of the well-studied non-reactive case of Rayleigh–Bénard convection with a hot lower boundary. The planar Burke–Schumann diffusion flame, which has been previously considered unconditionally stable in studies disregarding gravity, is shown to become unstable when the Rayleigh number exceeds a critical value. A numerical treatment is performed to test the effects of compressibility and finite chemistry on the stability of the system. For weak values of the thermal expansion coefficient $\alpha $, the numerical results show strong agreement with those of the linear stability analysis. It is found that as $\alpha $ increases to a more realistic value the system becomes considerably more stable, and also exhibits hysteresis at the onset of instability. Finally, a reduction in the Damköhler number is found to decrease the stability of the system.

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.

Weak Turbulence in Small Prandtl Number Convection in Porous Media

2010 ◽  
Author(s):  
Peter Vadasz
Keyword(s):  
Porous Media ◽  
Stability Analysis ◽  
Prandtl Number ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Weak Turbulence ◽  
Numerical Computations ◽  
Prandtl Numbers ◽  
Convection In Porous Media

The dynamics of weak turbulence in small Prandtl number convection in porous media is substantially distinct than the corresponding dynamics for moderate and large Prandtl numbers. Linear stability analysis is performed and its results compared with numerical computations to reveal the underlying phenomena.

scholarly journals Bioconvection under uniform shear: linear stability analysis

2013 ◽  
Vol 738 ◽  
pp. 522-562 ◽  
Author(s):  
Yongyun Hwang ◽  
T. J. Pedley
Keyword(s):  
Stability Analysis ◽  
Shear Rate ◽  
Linear Stability ◽  
Couette Flow ◽  
Linear Stability Analysis ◽  
Plane Couette Flow ◽  
Shear Rates ◽  
Uniform Shear ◽  
Rayleigh Benard Convection ◽  
Rayleigh Benard

AbstractThe role of uniform shear in bioconvective instability in a shallow suspension of swimming gyrotactic cells is studied using linear stability analysis. The shear is introduced by applying a plane Couette flow, and it significantly disturbs gravitaxis of the cell. The unstably stratified basic state of the cell concentration is gradually relieved as the shear rate is increased, and it even becomes stably stratified at very large shear rates. Stability of the basic state is significantly changed. The instability at high wavenumbers is drastically damped out with the shear rate, while that at low wavenumbers is destabilized. However, at very large shear rates, the latter is also suppressed. The most unstable mode is found as a pair of streamwise uniform rolls aligned with the shear, analogous to Rayleigh–Bénard convection in plane Couette flow. To understand these findings, the physical mechanism of the bioconvective instability is reexamined with several sets of numerical experiments. It is shown that the bioconvective instability in a shallow suspension originates from three different physical processes: gravitational overturning, gyrotaxis of the cell and negative cross-diffusion flux. The first mechanism is found to rule the behaviour of low-wavenumber instability whereas the last two mechanisms are mainly associated with high-wavenumber instability. With the increase of the shear rate, the former is enhanced, thereby leading to destabilization at low wavenumbers, whereas the latter two mechanisms are significantly suppressed. For streamwise varying perturbations, shear with sufficiently large rates is also found to play a stabilizing role as in Rayleigh–Bénard convection. However, at small shear rates, it destabilizes these perturbations through the mechanism of overstability discussed by Hill, Pedley and Kessler (J. Fluid Mech., vol. 208, 1989, pp. 509–543). Finally, the present findings are compared with a recent experiment by Croze, Ashraf and Bees (Phys. Biol., vol. 7, 2010, 046001) and they are in qualitative agreement.

Optimal Disturbances in Three-Dimensional Natural Convection Flows

2012 ◽  
Author(s):  
Elia Merzari ◽  
Paul Fischer ◽  
W. David Pointer
Keyword(s):  
Stability Analysis ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Natural Circulation ◽  
Three Dimensional ◽  
Boussinesq Approximation ◽  
Spectral Element ◽  
Base Flow ◽  
The Stability ◽  
Optimal Disturbances

Buoyancy-driven systems are subject to several types of flow instabilities. To evaluate the performance of such systems it is becoming increasingly crucial to be able to predict the stability of a given base flow configuration. Traditional Modal Linear stability Analysis requires the solution of very large eigenvalue systems for three-dimensional flows, which make this problem difficult to tackle. An alternative to modal Linear stability Analysis is the use of adjoint solvers [1] in combination with a power iteration [2]. Such methodology allows for the identification of an optimal disturbance or forcing and has been recently used to evaluate the stability of several isothermal flow systems [2]. In this paper we examine the extension of the methodology to non-isothermal flows driven by buoyancy. The contribution of buoyancy in the momentum equation is modeled through the Boussinesq approximation. The method is implemented in the spectral element code Nek5000. The test case is the flow is a two-dimensional cavity with differential heating and conductive walls and the natural circulation flow in a toroidal thermosiphon.

Sign in / Sign up

Export Citation Format

Share Document