Optimal Disturbances in Three-Dimensional Natural Convection Flows

Author(s):  
Elia Merzari ◽  
Paul Fischer ◽  
W. David Pointer

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.

2017 ◽  
Vol 822 ◽  
pp. 813-847 ◽  
Author(s):  
Azan M. Sapardi ◽  
Wisam K. Hussam ◽  
Alban Pothérat ◽  
Gregory J. Sheard

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.


2012 ◽  
Vol 711 ◽  
pp. 27-39 ◽  
Author(s):  
Bo-Fu Wang ◽  
Dong-Jun Ma ◽  
Cheng Chen ◽  
De-Jun Sun

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.


Author(s):  
Christoph Glocker

Curve squealing of railway wheels occurs in narrow curves with a frequency of about 4 kHz and is caused by a self-excited stick-slip oscillation in the wheel-rail contact. A squeal model, consisting of the first modal forms of an elastic wheel and a three-dimensional hard Coulomb contact, is developed to analyze the stability of the stationary run through a curve via linearization about the steady sliding equilibrium. The results show that in particular the front inner of the four wheels of a certain bogie tends to squeal.


2002 ◽  
Vol 455 ◽  
pp. 1-19 ◽  
Author(s):  
CHO LIK CHAN ◽  
WEN-YAU CHEN ◽  
C. F. CHEN

The three-dimensional motion observed by Chen & Chen (1997) in the convection cells generated by sideways heating of a solute gradient is further examined by experiments and linear stability analysis. In the experiments, we obtained visualizations and PIV measurements of the velocity of the fluid motion in the longitudinal plane perpendicular to the imposed temperature gradient. The flow consists of a horizontal row of counter-rotating vortices within each convection cell. The magnitude of this secondary motion is approximately one-half that of the primary convection cell. Results of a linear stability analysis of a parallel double-diffusive flow model of the actual ow show that the instability is in the salt-finger mode under the experimental conditions. The perturbation streamlines in the longitudinal plane at onset consist of a horizontal row of counter-rotating vortices similar to those observed in the experiments.


2001 ◽  
Author(s):  
Davide Valtorta ◽  
Khaled E. Zaazaa ◽  
Ahmed A. Shabana ◽  
Jalil R. Sany

Abstract The lateral stability of railroad vehicles travelling on tangent tracks is one of the important problems that has been the subject of extensive research since the nineteenth century. Early detailed studies of this problem in the twentieth century are the work of Carter and Rocard on the stability of locomotives. The linear theory for the lateral stability analysis has been extensively used in the past and can give good results under certain operating conditions. In this paper, the results obtained using a linear stability analysis are compared with the results obtained using a general nonlinear multibody methodology. In the linear stability analysis, the sources of the instability are investigated using Liapunov’s linear theory and the eigenvalue analysis for a simple wheelset model on a tangent track. The effects of the stiffness of the primary and secondary suspensions on the stability results are investigated. The results obtained for the simple model using the linear approach are compared with the results obtained using a new nonlinear multibody based constrained wheel/rail contact formulation. This comparative numerical study can be used to validate the use of the constrained wheel/rail contact formulation in the study of lateral stability. Similar studies can be used in the future to define the limitations of the linear theory under general operating conditions.


Fluids ◽  
2021 ◽  
Vol 6 (10) ◽  
pp. 342
Author(s):  
Maksims Zigunovs ◽  
Andrei Kolyshkin ◽  
Ilmars Iltins

Linear stability analysis of a steady convective flow in a tall vertical annulus caused by nonlinear heat sources is conducted in the paper. Heat sources are generated as a result of a chemical reaction. The effect of radial cross-flow through permeable porous walls of the annulus is analyzed. The problem is relevant to biomass thermal conversion. The base flow solution is obtained by solving nonlinear boundary value problem. Linear stability analysis is performed, using collocation method. The calculations show that radial inward or outward flow has a stabilizing effect on the flow, while the increase in the Frank–Kamenetskii parameter (proportional to the intensity of the chemical reaction) destabilizes the flow. The increase in the Reynolds number based on the radial velocity leads to the appearance of the second minimum on the marginal stability curves. The rate of increase in the critical Grashof number with respect to the Reynolds number is different for inward and outward radial flows.


2018 ◽  
Vol 192 ◽  
pp. 892-905 ◽  
Author(s):  
Saikrishna Nadella ◽  
Abhishek Kumar Srivastava ◽  
Naresh Kumar Maheshwari

2011 ◽  
Vol 676 ◽  
pp. 110-144 ◽  
Author(s):  
P. BOHORQUEZ ◽  
E. SANMIGUEL-ROJAS ◽  
A. SEVILLA ◽  
J. I. JIMÉNEZ-GONZÁLEZ ◽  
C. MARTÍNEZ-BAZÁN

We investigate the stability properties and flow regimes of laminar wakes behind slender cylindrical bodies, of diameter D and length L, with a blunt trailing edge at zero angle of attack, combining experiments, direct numerical simulations and local/global linear stability analyses. It has been found that the flow field is steady and axisymmetric for Reynolds numbers below a critical value, Recs (L/D), which depends on the length-to-diameter ratio of the body, L/D. However, in the range of Reynolds numbers Recs(L/D) < Re < Reco(L/D), although the flow is still steady, it is no longer axisymmetric but exhibits planar symmetry. Finally, for Re > Reco, the flow becomes unsteady due to a second oscillatory bifurcation which preserves the reflectional symmetry. In addition, as the Reynolds number increases, we report a new flow regime, characterized by the presence of a secondary, low frequency oscillation while keeping the reflectional symmetry. The results reported indicate that a global linear stability analysis is adequate to predict the first bifurcation, thereby providing values of Recs nearly identical to those given by the corresponding numerical simulations. On the other hand, experiments and direct numerical simulations give similar values of Reco for the second, oscillatory bifurcation, which are however overestimated by the linear stability analysis due to the use of an axisymmetric base flow. It is also shown that both bifurcations can be stabilized by injecting a certain amount of fluid through the base of the body, quantified here as the bleed-to-free-stream velocity ratio, Cb = Wb/W∞.


Sign in / Sign up

Export Citation Format

Share Document