scholarly journals Three-dimensional Floquet stability analysis of the wake of a circular cylinder

1996 ◽  
Vol 322 ◽  
pp. 215-241 ◽  
Author(s):  
Dwight Barkley ◽  
Ronald D. Henderson

Results are reported from a highly accurate, global numerical stability analysis of the periodic wake of a circular cylinder for Reynolds numbers between 140 and 300. The analysis shows that the two-dimensional wake becomes (absolutely) linearly unstable to three-dimensional perturbations at a critical Reynolds number of 188.5±1.0. The critical spanwise wavelength is 3.96 ± 0.02 diameters and the critical Floquet mode corresponds to a ‘Mode A’ instability. At Reynolds number 259 the two-dimensional wake becomes linearly unstable to a second branch of modes with wavelength 0.822 diameters at onset. Stability spectra and corresponding neutral stability curves are presented for Reynolds numbers up to 300.

Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three-dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio of in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations fail to capture completely the total expansion effect on the flow, which couples both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When the two-dimensional geometry was consistent with the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations was compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, helping to sustain laminar flow symmetry to higher Reynolds numbers in comparison to nominally two-dimensional double-expansion geometries. Lastly, and most importantly, when the logarithm of the critical Reynolds number from all these studies was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


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):  
A. Inasawa ◽  
K. Toda ◽  
M. Asai

Disturbance growth in the wake of a circular cylinder moving at a constant acceleration is examined experimentally. The cylinder is installed on a carriage moving in the still air. The results show that the critical Reynolds number for the onset of the global instability leading to a self-sustained wake oscillation increases with the magnitude of acceleration, while the Strouhal number of the growing disturbance at the critical Reynolds number is not strongly dependent on the magnitude of acceleration. It is also found that with increasing the acceleration, the Ka´rma´n vortex street remains two-dimensional even at the Reynolds numbers around 200 where the three-dimensional instability occurs to lead to the vortex dislocation in the case of cylinder moving at constant velocity or in the case of cylinder wake in the steady oncoming flow.


2014 ◽  
Vol 751 ◽  
pp. 570-600 ◽  
Author(s):  
Mark C. Thompson ◽  
Alexander Radi ◽  
Anirudh Rao ◽  
John Sheridan ◽  
Kerry Hourigan

AbstractWhile the wake of a circular cylinder and, to a lesser extent, the normal flat plate have been studied in considerable detail, the wakes of elliptic cylinders have not received similar attention. However, the wakes from the first two bodies have considerably different characteristics, in terms of three-dimensional transition modes, and near- and far-wake structure. This paper focuses on elliptic cylinders, which span these two disparate cases. The Strouhal number and drag coefficient variations with Reynolds number are documented for the two-dimensional shedding regime. There are considerable differences from the standard circular cylinder curve. The different three-dimensional transition modes are also examined using Floquet stability analysis based on computed two-dimensional periodic base flows. As the cylinder aspect ratio (major to minor axis) is decreased, mode A is no longer unstable for aspect ratios below 0.25, as the wake deviates further from the standard Bénard–von Kármán state. For still smaller aspect ratios, another three-dimensional quasi-periodic mode becomes unstable, leading to a different transition scenario. Interestingly, for the 0.25 aspect ratio case, mode A restabilises above a Reynolds number of approximately 125, allowing the wake to return to a two-dimensional state, at least in the near wake. For the flat plate, three-dimensional simulations show that the shift in the Strouhal number from the two-dimensional value is gradual with Reynolds number, unlike the situation for the circular cylinder wake once mode A shedding develops. Dynamic mode decomposition is used to characterise the spatially evolving character of the wake as it undergoes transition from the primary Bénard–von Kármán-like near wake into a two-layered wake, through to a secondary Bénard–von Kármán-like wake further downstream, which in turn develops an even longer wavelength unsteadiness. It is also used to examine the differences in the two- and three-dimensional near-wake state, showing the increasing distortion of the two-dimensional rollers as the Reynolds number is increased.


2017 ◽  
Vol 826 ◽  
pp. 32-59 ◽  
Author(s):  
F. Y. Houdroge ◽  
T. Leweke ◽  
K. Hourigan ◽  
M. C. Thompson

This paper presents the characteristics of the different stages in the evolution of the wake of a circular cylinder rolling without slipping along a wall at constant speed, acquired through numerical stability analysis and two- and three-dimensional numerical simulations. Reynolds numbers between 30 and 300 are considered. Of importance in this study is the transition to three-dimensionality from the underlying two-dimensional periodic flow and, in particular, the way that the associated transitions influence the fluid forces exerted on the cylinder and the development and the structure of the wake. It is found that the steady two-dimensional flow becomes unstable to three-dimensional perturbations at $Re_{c,3D}=37$, and that the transition to unsteady two-dimensional flow – or periodic vortex shedding – occurs at $Re_{c,2D}=88$, thus validating and refining the results of Stewart et al. (J. Fluid Mech. vol. 648, 2010, pp. 225–256). The main focus here is on Reynolds numbers beyond the transition to unsteady flow at $Re_{c,2D}=88$. From impulsive start up, the wake almost immediately undergoes transition to a periodic two-dimensional wake state, which, in turn, is three-dimensionally unstable. Thus, the previous three-dimensional stability analysis based on the two-dimensional steady flow provides only an element of the full story. Floquet analysis based on the periodic two-dimensional flow was undertaken and new three-dimensional instability modes were revealed. The results suggest that an impulsively started cylinder rolling along a surface at constant velocity for $Re\gtrsim 90$ will result in the rapid development of a periodic two-dimensional wake that will be maintained for a considerable time prior to the wake undergoing three-dimensional transition. Of interest, the mean lift and drag coefficients obtained from full three-dimensional simulations match predictions from two-dimensional simulations to within a few per cent.


2005 ◽  
Vol 128 (4) ◽  
pp. 671-679 ◽  
Author(s):  
Francine Battaglia ◽  
George Papadopoulos

The effect of three dimensionality on low Reynolds number flows past a symmetric sudden expansion in a channel was investigated. The geometric expansion ratio in the current study was 2:1 and the aspect ratio was 6:1. Both experimental velocity measurements and two- and three-dimensional simulations for the flow along the centerplane of the rectangular duct are presented for Reynolds numbers in the range of 150 to 600. Comparison of the two-dimensional simulations with the experiments revealed that the simulations failed to capture completely the total expansion effect on the flow, which is influenced by both geometric and hydrodynamic effects. To properly do so requires the definition of an effective expansion ratio, which is the ratio of the downstream and upstream hydraulic diameters and is therefore a function of both the expansion and aspect ratios. When two-dimensional simulations were performed using the effective expansion ratio, the new results agreed well with the three-dimensional simulations and the experiments. Furthermore, in the range of Reynolds numbers investigated, the laminar flow through the expansion underwent a symmetry-breaking bifurcation. The critical Reynolds number evaluated from the experiments and the simulations were compared to other values reported in the literature. Overall, side-wall proximity was found to enhance flow stability, thus sustaining laminar flow symmetry to higher Reynolds numbers. Last, and most important, when the logarithm of the critical Reynolds number was plotted against the reciprocal of the effective expansion ratio, a linear trend emerged that uniquely captured the bifurcation dynamics of all symmetric double-sided planar expansions.


2018 ◽  
Vol 841 ◽  
pp. 636-653
Author(s):  
Ting-Yueh Chang ◽  
Falin Chen ◽  
Min-Hsing Chang

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.


Author(s):  
Xuemin Ye ◽  
Weiping Yan ◽  
Chunxi Li

When liquid film is under evaporating or condensing conditions, the flow stability is clearly different to that under isothermal condition due to thermal non-equilibrium effect at interface, especially under lower Reynolds number. The universal linear temporal and spatial stability formulations of the two-dimensional surface waves on evaporating or isothermal or condensing liquid films are established in present paper with the collocation method based on the boundary layer theory and complete boundary conditions. The models include the effects of Reynolds number, thermocapillarity, inclination angle, liquid property, evaporation, isothermal or condensation. The effects of above factors are investigated with the neutral stability curves at different Reynolds numbers, and stabilities characteristics are fully indicated in theory for evaporating or condensing films.


2011 ◽  
Vol 681 ◽  
pp. 411-433 ◽  
Author(s):  
HEMANT K. CHAURASIA ◽  
MARK C. THOMPSON

A detailed numerical study of the separating and reattaching flow over a square leading-edge plate is presented, examining the instability modes governing transition from two- to three-dimensional flow. Under the influence of background noise, experiments show that the transition scenario typically is incompletely described by either global stability analysis or the transient growth of dominant optimal perturbation modes. Instead two-dimensional transition effectively can be triggered by the convective Kelvin–Helmholtz (KH) shear-layer instability; although it may be possible that this could be described alternatively in terms of higher-order optimal perturbation modes. At least in some experiments, observed transition occurs by either: (i) KH vortices shedding downstream directly and then almost immediately undergoing three-dimensional transition or (ii) at higher Reynolds numbers, larger vortical structures are shed that are also three-dimensionally unstable. These two paths lead to distinctly different three-dimensional arrangements of vortical flow structures. This paper focuses on the mechanisms underlying these three-dimensional transitions. Floquet analysis of weakly periodically forced flow, mimicking the observed two-dimensional quasi-periodic base flow, indicates that the two-dimensional vortex rollers shed from the recirculation region become globally three-dimensionally unstable at a Reynolds number of approximately 380. This transition Reynolds number and the predicted wavelength and flow symmetries match well with those of the experiments. The instability appears to be elliptical in nature with the perturbation field mainly restricted to the cores of the shed rollers and showing the spatial vorticity distribution expected for that instability type. Indeed an estimate of the theoretical predicted wavelength is also a good match to the prediction from Floquet analysis and theoretical estimates indicate the growth rate is positive. Fully three-dimensional simulations are also undertaken to explore the nonlinear development of the three-dimensional instability. These show the development of the characteristic upright hairpins observed in the experimental dye visualisations. The three-dimensional instability that manifests at lower Reynolds numbers is shown to be consistent with an elliptic instability of the KH shear-layer vortices in both symmetry and spanwise wavelength.


Author(s):  
Majid Nabavi ◽  
Luc Mongeau

In this study, two-dimensional laminar incompressible and turbulent compressible flow through the planar diffuser (gradual expansion) for different divergence half angles of the diffuser (θ), and different Reynolds numbers (Re) was numerically studied. The effects of θ on the critical Reynolds number at which the onset of asymmetric flow is observed, were investigated. In the laminar flow regime, it was observed that for every values of θ, there is a critical Re beyond which the flow is asymmetric. However, in the turbulent flow regime, for θ ≥ 20°, even at low Reynolds number the flow is asymmetric. Only for θ ≤ 10°, symmetric flow was observed below a critical Re.


Sign in / Sign up

Export Citation Format

Share Document