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

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.

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.


2021 ◽  
Vol 16 ◽  
pp. 34-42
Author(s):  
Nitin Kumar ◽  
Sunil Chamoli ◽  
Sachin Tejyan ◽  
Pawan Kumar Pant

A discrete linear stability analysis framework for two-dimensional laminar flows is presented. Using two case studies involving analysis of thermal and laminar flows, the stability of flows in the discrete numerical sense is addressed. The two-dimensional base flow for various values of the controlling parameter (Reynolds number for flow past a square cylinder and Rayleigh number for double-glazing problem) is computed numerically by using the lattice Boltzmann method. The governing equations, discretized using the finitedifference method in two-dimensions and are subsequently written in the form of perturbed equations with twodimensional disturbances. These equations are linearized around the base flow and form a set of partial differential equations that govern the evolution of the perturbations. The eigenvalues, stability of the base flow and the points of bifurcations are determined using normal mode analysis. The eigenvalue spectrum predicts that the critical Reynolds number is 52 and the critical Rayleigh number is 6 1.88×10 for the square cylinder and double-glazing problem, respectively, The results are consistent with the previous numerical and experimental observations.


2021 ◽  
Vol 16 ◽  
pp. 109-119
Author(s):  
Nitin Kumar ◽  
Sachin Tejyan ◽  
Sunil Chamoli ◽  
Pawan Kumar Pant

The present study focuses on the development of a numerical framework for predicting the onset of vortex sheading due to flow past a square cylinder. For this a discrete linear stability analysis framework for two-dimensional laminar flows have used. Initially the frame work is validating by using the analysis of thermal stability of flows in the discrete numerical sense. The two-dimensional base flow for various values of the controlling parameter (Reynolds number for flow past a square cylinder and Rayleigh number for double-glazing problem) is computed numerically by using the lattice Boltzmann method. The governing equations, discretized using the finite-difference method in two-dimensions and are subsequently written in the form of perturbed equations with two-dimensional disturbances. These equations are linearized around the base flow and form a set of partial differential equations that govern the evolution of the perturbations. The eigenvalues, stability of the base flow and the points of bifurcations are determined using normal mode analysis. The eigenvalue spectrum predicts that the critical Reynolds number is 52 for the flow past a square cylinder. The results are consistent with the previous numerical and experimental observations.


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.


Fluids ◽  
2020 ◽  
Vol 5 (4) ◽  
pp. 212
Author(s):  
Miles Owen ◽  
Abdelkader Frendi

The results from a temporal linear stability analysis of a subsonic boundary layer over a flat plate with a straight and wavy leading edge are presented in this paper for a swept and un-swept plate. For the wavy leading-edge case, an extensive study on the effects of the amplitude and wavelength of the waviness was performed. Our results show that the wavy leading edge increases the critical Reynolds number for both swept and un-swept plates. For the un-swept plate, increasing the leading-edge amplitude increased the critical Reynolds number, while changing the leading-edge wavelength had no effect on the mean flow and hence the flow stability. For the swept plate, a local analysis at the leading-edge peak showed that increasing the leading-edge amplitude increased the critical Reynolds number asymptotically, while the leading-edge wavelength required optimization. A global analysis was subsequently performed across the span of the swept plate, where smaller leading-edge wavelengths produced relatively constant critical Reynolds number profiles that were larger than those of the straight leading edge, while larger leading-edge wavelengths produced oscillating critical Reynolds number profiles. It was also found that the most amplified wavenumber was not affected by the wavy leading-edge geometry and hence independent of the waviness.


Author(s):  
Paul Ziadé ◽  
Pierre E. Sullivan

Large-eddy simulation and linear stability analysis were performed on a NACA 0025 airfoil at a chord Reynolds number of 105 and four angles of attack. The computations showed that the initial vortex roll-up quickly breaks down to three-dimensional turbulence. Flow separation was observed at all angles, whereas only the lowest angle of attack formed a laminar separation bubble due to flow transition occuring close to the airfoil surface. A Chebyshev collocation method was employed to solve the viscous and inviscid stability equations. Linear stability analysis demonstrated that high-frequency disturbances occur in the laminar separation bubble case, whereas lower frequencies are present for the fully separated angles of attack. The maximum disturbance growth rates were dampened with the addition of viscosity but negligible change in peak frequency was noted.


2015 ◽  
Vol 788 ◽  
pp. 129-146 ◽  
Author(s):  
Jānis Priede ◽  
Thomas Arlt ◽  
Leo Bühler

This study is concerned with the numerical linear stability analysis of liquid-metal flow in a square duct with thin electrically conducting walls subject to a uniform transverse magnetic field. We derive an asymptotic solution for the base flow that is valid for not only high but also moderate magnetic fields. This solution shows that, for low wall conductance ratios $c\ll 1$, an extremely strong magnetic field with Hartmann number $\mathit{Ha}\sim c^{-4}$ is required to attain the asymptotic flow regime considered in previous studies. We use a vector streamfunction–vorticity formulation and a Chebyshev collocation method to solve the eigenvalue problem for three-dimensional small-amplitude perturbations in ducts with realistic wall conductance ratios $c=1$, 0.1 and 0.01 and Hartmann numbers up to $10^{4}$. As for similar flows, instability in a sufficiently strong magnetic field is found to occur in the sidewall jets with characteristic thickness ${\it\delta}\sim \mathit{Ha}^{-1/2}$. This results in the critical Reynolds number and wavenumber increasing asymptotically with the magnetic field as $\mathit{Re}_{c}\sim 110\mathit{Ha}^{1/2}$ and $k_{c}\sim 0.5\mathit{Ha}^{1/2}$. The respective critical Reynolds number based on the total volume flux in a square duct with $c\ll 1$ is $\overline{\mathit{Re}}_{c}\approx 520$. Although this value is somewhat larger than $\overline{\mathit{Re}}_{c}\approx 313$ found by Ting et al. (Intl J. Engng Sci., vol. 29 (8), 1991, pp. 939–948) for the asymptotic sidewall jet profile, it still appears significantly lower than the Reynolds numbers at which turbulence is observed in experiments as well as in direct numerical simulations of this type of flow.


2007 ◽  
Vol 582 ◽  
pp. 319-340 ◽  
Author(s):  
M. D. GRIFFITH ◽  
M. C. THOMPSON ◽  
T. LEWEKE ◽  
K. HOURIGAN ◽  
W. P. ANDERSON

The two-dimensional flow through a constricted channel is studied. A semi-circular bump is located on one side of the channel and the extent of blockage is varied by adjusting the radius of the bump. The blockage is varied between 0.05 and 0.9 of the channel width and the upstream Reynolds number between 25 and 3000. The geometry presents a simplified blockage specified by a single parameter, serving as a starting point for investigations of other more complex blockage geometries. For blockage ratios in excess of 0.4, the variation of reattachment length with Reynolds number collapses to within approximately 15%, while at lower ratios the behaviour differs. For the constrained two-dimensional flow, various phenomena are identified, such as multiple mini-recirculations contained within the main recirculation bubble and vortex shedding at higher Reynolds numbers. The stability of the flow to three-dimensional perturbations is analysed, revealing a transition to a three-dimensional state at a critical Reynolds number which decreases with higher blockage ratios. Separation lengths and the onset and structure of three-dimensional instability observed from the geometry of blockage ratio 0.5 resemble results taken from backward-facing step investigations. The question of the underlying mechanism behind the instability being either centrifugal or elliptic in nature and operating within the initial recirculation zone is analytically tested.


Sign in / Sign up

Export Citation Format

Share Document