scholarly journals Influence of confinement on a two-dimensional wake

2011 ◽  
Vol 688 ◽  
pp. 297-320 ◽  
Author(s):  
Luca Biancofiore ◽  
François Gallaire ◽  
Richard Pasquetti
Keyword(s):  
Stability Analysis ◽  
Velocity Ratio ◽  
Reynolds Numbers ◽  
Base Flow ◽  
Wake Flow ◽  
Two Dimensional ◽  
Dimensional Parameters ◽  
Spatio Temporal ◽  
Spatial Locations ◽  
Theoretical Results

AbstractThe spatio-temporal development of an incompressible two-dimensional viscous wake flow confined by two flat slipping plates is investigated by means of direct numerical simulation (DNS), using a spectral Chebyshev multi-domain method. The limit between unstable and stable configurations is determined with respect to several non-dimensional parameters: the confinement, the velocity ratio and two different Reynolds numbers, $100$ and $500$. The comparison of such limit curves with theoretical results obtained by Juniper (J. Fluid Mech., vol. 565, 2006, pp. 171–195) confirms the existence of a region at moderate confinement where the instability is maximal. Moreover, instabilities are also observed under sustained co-flow, in the form of a vacillating front. Using a direct computation of the two-dimensional base flow, we perform a local linear stability analysis for several velocity profiles prevailing at different spatial locations, so as to determine the local spatio-temporal nature of the flow: convectively unstable or absolutely unstable. Comparisons of the DNS and local stability analysis results are provided and discussed.

Stability and dynamics of the laminar wake past a slender blunt-based axisymmetric body

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
Keyword(s):  
Stability Analysis ◽  
Numerical Simulations ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Velocity Ratio ◽  
Reynolds Numbers ◽  
Base Flow ◽  
Direct Numerical Simulations ◽  
The Body ◽  
Stream Velocity

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∞.

Stability of the wakes of cylinders with triangular cross-sections

2018 ◽  
Vol 844 ◽  
pp. 721-745 ◽  
Author(s):  
Zhi Y. Ng ◽  
Tony Vo ◽  
Gregory J. Sheard
Keyword(s):  
Stability Analysis ◽  
Cross Sections ◽  
Reynolds Numbers ◽  
Base Flow ◽  
Half Period ◽  
Two Dimensional ◽  
Instability Modes ◽  
The Stability ◽  
Time Periodic ◽  
Mode A

The stability of the wakes of cylinders with triangular cross-sections at incidence is investigated using Floquet stability analysis to elucidate the effects of cylinder inclination on the dominant flow instability. The upper limit of the Reynolds numbers (scaled by the height projected by the cylinder in this study) at which the wake of the two-dimensional base flow is time periodic is$Re\approx 140$for most cylinder inclinations, exceeding which the flow becomes aperiodic, restricting the range of Reynolds numbers permitted for the stability analysis. Two different instability modes are predicted to manifest as the first-occurring mode at various cylinder inclinations – a regular mode possessing perturbation structures consistent with mode A dominates the wakes of cylinders at inclinations$\unicode[STIX]{x1D6FC}\lesssim 34.6^{\circ }$and$\unicode[STIX]{x1D6FC}\gtrsim 55.4^{\circ }$, with a subharmonic mode consistent with mode C emerging as the primary mode in the wakes of the cylinder at the intermediate range of inclinations. For all inclinations, the mode B branch is not detected within the range of Reynolds numbers examined. The peak instability growth rates corresponding to mode A for all cylinder inclinations describe a linear variation with$(Re-Re_{A})/Re_{A}$, where$Re_{A}$is the mode A transition Reynolds number, while those corresponding to mode C vary only approximately linearly. The generalized trend most pertinently shows mode C to develop more rapidly than mode A at inclinations which permit it. Examination of the near wake of the two-dimensional time-periodic base flow demonstrates the dependence of the development and intensity of mode C on imbalances in the flow solution over each shedding period, directly implying that the two-dimensional base flow solutions deviate from the half-period-flip map as the cylinder inclination is increased. The degree of asymmetry of the two-dimensional base flow relative to the ideal half-period-flip map is quantified using several measures. The results show distinctly different trends in these asymmetry measures between inclinations where modes A or C are dominant, agreeing with results from the stability analysis. The nature of the predicted instability modes at transition are also investigated by applying the Stuart–Landau equation, showing the transitions to be supercritical for all cylinder inclinations, with mode C being consistently more strongly supercritical than mode A.

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.

Transition and turbulence in horizontal convection: linear stability analysis

2017 ◽  
Vol 821 ◽  
pp. 31-58 ◽  
Author(s):  
Pierre-Yves Passaggia ◽  
Alberto Scotti ◽  
Brian White
Keyword(s):  
Stability Analysis ◽  
Rayleigh Number ◽  
Boundary Conditions ◽  
Stokes Equations ◽  
Critical Rayleigh Number ◽  
Base Flow ◽  
Two Dimensional ◽  
Slip Boundary Conditions ◽  
Slip Boundary ◽  
Horizontal Convection

The linear instability mechanisms of horizontal convection in a rectangular cavity forced by a horizontal buoyancy gradient along its surface are investigated using local and global stability analyses for a Prandtl number equal to unity. The results show that the stability of the base flow, a steady circulation characterized by a narrow descending plume and a broad upwelling region, depends on the Rayleigh number, $Ra$. For free-slip boundary conditions at a critical value of $Ra\approx 2\times 10^{7}$, the steady base flow becomes unstable to three-dimensional perturbations, characterized by counter-rotating vortices originating within the plume region. A Wentzel–Kramers–Brillouin (WKB) method applied along closed streamlines demonstrates that this instability is of a Rayleigh–Taylor type and can be used to accurately reconstruct the global instability mode. In the case of no-slip boundary conditions, the base flow also becomes unstable to a self-sustained two-dimensional instability whose critical Rayleigh number is $Ra\approx 1.7\times 10^{8}$. Beyond this critical $Ra$, two-dimensional equilibrium stationary states of the Navier–Stokes equations are computed using the selective frequency damping method. The two-dimensional onset of instability is shown to be characterized by a family of modes also originating within the plume. A local spatio-temporal stability analysis shows that the flow becomes absolutely unstable at the origin of the plume. Taken together, these results illustrate the mechanisms behind the onset of turbulence that has been observed in horizontal convection.

scholarly journals A Mathematical Model of the Transition from the Normal Hematopoiesis to the Chronic and Accelerated Acute Stages in Myeloid Leukemia

2020 ◽  
Author(s):  
Lorand Gabriel Parajdi ◽  
Radu Precup ◽  
Eduard Alexandru Bonci ◽  
Ciprian Tomuleasa
Keyword(s):  
Mathematical Model ◽  
Stem Cells ◽  
Bone Marrow ◽  
Stability Analysis ◽  
Myeloid Leukemia ◽  
Transition Process ◽  
Two Dimensional ◽  
The Stability ◽  
Theoretical Results

A mathematical model given by a two - dimensional differential system is introduced in order to understand the transition process from the normal hematopoiesis to the chronic and accelerated acute stages in chronic myeloid leukemia. A previous model of Dingli and Michor is refined by introducing a new parameter in order to differentiate the bone marrow microenvironment sensitivities of normal and mutant stem cells. In the light of the new parameter, the system now has three distinct equilibria corresponding to the normal hematopoietic state, to the chronic state, and to the accelerated acute phase of the disease. A characterization of the three hematopoietic states is obtained based on the stability analysis. Numerical simulations are included to illustrate the theoretical results.

The stability of the variable-density Kelvin–Helmholtz billow

2008 ◽  
Vol 612 ◽  
pp. 237-260 ◽  
Author(s):  
JÉRÔME FONTANE ◽  
LAURENT JOLY
Keyword(s):  
Stability Analysis ◽  
Mixing Layer ◽  
Three Dimensional ◽  
Variable Density ◽  
Base Flow ◽  
Shear Instability ◽  
Numerical Continuation ◽  
Small Scale ◽  
Two Dimensional ◽  
Two Fluids

We perform a three-dimensional stability analysis of the Kelvin–Helmholtz (KH) billow, developing in a shear layer between two fluids with different density. We begin with two-dimensional simulations of the temporally evolving mixing layer, yielding the unsteady base flow fields. The Reynolds number is 1500 while the Schmidt and Froude numbers are infinite. Then exponentially unstable modes are extracted from a linear stability analysis performed at the saturation of the primary mode kinetic energy. The spectrum of the least stable modes exhibits two main classes. The first class comprises three-dimensional core-centred and braid-centred modes already present in the homogeneous case. The baroclinic vorticity concentration in the braid lying on the light side of the KH billow turns the flow into a sharp vorticity ridge holding high shear levels. The hyperbolic modes benefit from the enhanced level of shear in the braid whereas elliptic modes remain quite insensitive to the modifications of the base flow. In the second class, we found typical two-dimensional modes resulting from a shear instability of the curved vorticity-enhanced braid. For a density contrast of 0.5, the wavelength of the two-dimensional instability is about ten times shorter than that of the primary wave. Its amplification rate competes well against those of the hyperbolic three-dimensional modes. The vorticity-enhanced braid thus becomes the preferred location for the development of secondary instabilities. This stands as the key feature of the transition of the variable-density mixing layer. We carry out a fully resolved numerical continuation of the nonlinear development of the two-dimensional braid-mode. Secondary roll-ups due to a small-scale Kelvin–Helmholtz mechanism are promoted by the underlying strain field and develop rapidly in the compression part of the braid. Originally analysed by Reinoud et al. (Phys. Fluids, vol. 12, 2000, p. 2489) from two-dimensional non-viscous numerical simulations, this instability is shown to substantially increase the mixing.

scholarly journals Absolute instability of an annular jet: local stability analysis

Meccanica ◽  
2020 ◽  
Vol 55 (11) ◽  
pp. 2179-2198
Author(s):  
A. Boguslawski ◽  
K. Wawrzak
Keyword(s):  
Stability Analysis ◽  
Shear Layer ◽  
Flow Velocity ◽  
Azimuthal Velocity ◽  
Base Flow ◽  
Back Flow ◽  
Annular Jet ◽  
Swirl Intensity ◽  
The Absolute ◽  
Spatio Temporal

Abstract The paper presents parametric studies of the first and second azimuthal absolute modes in annular non-swirling and swirling jets. The spatio-temporal linear stability analysis is applied to investigate an influence of governing parameters including axial velocity gradients in inner and outer shear layers, back-flow velocity, swirl number and shape of the azimuthal velocity. A new base flow is formulated allowing a flexible variation of the shape of axial and azimuthal velocity profiles. It is shown that the first helical absolute mode is governed mainly by the back-flow velocity and swirl intensity. A steepness of the inner shear layer can control the absolute mode frequency. The velocity gradient in the outer shear layer and the shape of the azimuthal velocity have rather limited impact on the absolute mode characteristics. Finally, it is shown that the second helical absolute mode can dominate the flow with a stronger swirl intensity.

The wake behind a cylinder rolling on a wall at varying rotation rates

2010 ◽  
Vol 648 ◽  
pp. 225-256 ◽  
Author(s):  
B. E. STEWART ◽  
M. C. THOMPSON ◽  
T. LEWEKE ◽  
K. HOURIGAN
Keyword(s):  
Reynolds Number ◽  
Three Dimensional ◽  
Reynolds Numbers ◽  
Wake Flow ◽  
Recirculation Region ◽  
Two Dimensional ◽  
Number Range ◽  
Rotation Rates ◽  
Transition Mechanism ◽  
Three Dimensional Simulations

A study investigating the flow around a cylinder rolling or sliding on a wall has been undertaken in two and three dimensions. The cylinder motion is specified from a set of five discrete rotation rates, ranging from prograde through to retrograde rolling. A Reynolds number range of 20–500 is considered. The effects of the nearby wall and the imposed body motion on the wake structure and dominant wake transitions have been determined. Prograde rolling is shown to destabilize the wake flow, while retrograde rotation delays the onset of unsteady flow to Reynolds numbers well above those observed for a cylinder in an unbounded flow.Two-dimensional simulations show the presence of two recirculation zones in the steady wake, the lengths of which increase approximately linearly with the Reynolds number. Values of the lift and drag coefficient are also reported for the steady flow regime. Results from a linear stability analysis show that the wake initially undergoes a regular bifurcation from a steady two-dimensional flow to a steady three-dimensional wake for all rotation rates. The critical Reynolds number Rec of transition and the spanwise wavelength of the dominant mode are shown to be highly dependent on, but smoothly varying with, the rotation rate of the cylinder. Varying the rotation from prograde to retrograde rolling acts to increase the value of Rec and decrease the preferred wavelength. The structure of the fully evolved wake mode is then established through three-dimensional simulations. In fact it is found that at Reynolds numbers only marginally (~5%) above critical, the three-dimensional simulations indicate that the saturated state becomes time dependent, although at least initially, this does not result in a significant change to the mode structure. It is only at higher Reynolds numbers that the wake undergoes a transition to vortex shedding.An analysis of the three-dimensional transition indicates that it is unlikely to be due to a centrifugal instability despite the superficial similarity to the flow over a backward-facing step, for which the transition mechanism has been speculated to be centrifugal. However, the attached elongated recirculation region and distribution of the spanwise perturbation vorticity field, and the similarity of these features with those of the flow through a partially blocked channel, suggest the possibility that the mechanism is elliptic in nature. Some analysis which supports this conjecture is undertaken.

scholarly journals Onset of global instability in the flow past a circular cylinder cascade

2010 ◽  
Vol 668 ◽  
pp. 304-334 ◽  
Author(s):  
V. B. L. BOPPANA ◽  
J. S. B. GAJJAR
Keyword(s):  
Stability Analysis ◽  
Reynolds Numbers ◽  
Oscillatory Motion ◽  
Two Dimensional ◽  
Global Instability ◽  
Global Stability Analysis ◽  
Generalized Eigenvalue ◽  
Flow Past ◽  
Flow Past A Cylinder ◽  
Critical Reynolds Numbers

The effect of blockage on the onset of instability in the two-dimensional uniform flow past a cascade of cylinders is investigated. The same techniques as those described in Gajjar & Azzam (J. Fluid Mech., vol. 520, 2004, p. 51) are used to tackle the generalized eigenvalue problem arising from a global stability analysis of the linearized disturbance equations. Results have been obtained for the various mode classes, and our results show that for the odd–even modes, which correspond to anti-phase oscillatory motion about the midplane between the cylinders and are the modes most extensively studied in the literature, the effect of blockage has a marginal influence on the critical Reynolds numbers for instability. This is in sharp contrast to results cited in many studies with a fully developed inlet flow past a cylinder placed between confining walls. We are also able to find other unstable modes and in particular for low blockage ratios, the odd–odd modes which correspond to the in-phase oscillatory motion about the midplane between the cylinders are the first to become unstable as compared with the odd–even modes, and with much lower frequencies.

scholarly journals A Discrete Linear Stability Analysis Framework for Two-Dimensional Laminar Flows

2021 ◽  
Vol 16 ◽  
pp. 34-42
Author(s):  
Nitin Kumar ◽  
Sunil Chamoli ◽  
Sachin Tejyan ◽  
Pawan Kumar Pant
Keyword(s):  
Reynolds Number ◽  
Stability Analysis ◽  
Rayleigh Number ◽  
Linear Stability ◽  
Linear Stability Analysis ◽  
Laminar Flows ◽  
Base Flow ◽  
Two Dimensional ◽  
Analysis Framework ◽  
Square Cylinder

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.

Sign in / Sign up

Export Citation Format

Share Document