Flow instabilities in the wake of a rounded square cylinder

2016 ◽  
Vol 793 ◽  
pp. 915-932 ◽  
Author(s):  
Doohyun Park ◽  
Kyung-Soo Yang

Instabilities in the flow past a rounded square cylinder have been numerically studied in order to clarify the effects of rounding the sharp edges of a square cylinder on the primary and secondary instabilities associated with the flow. Rounding the edges was done by inscribing a quarter circle of radius $r$ in each edge of a square cylinder of height $d$. Nine cases of rounding were considered, ranging from a square cylinder ($r/d=0$) to a circular cylinder ($r/d=0.5$) with an increment of 0.0625. Each cross-section was numerically implemented in a Cartesian grid system by using an immersed boundary method. The key parameters are the Reynolds number (Re) and the edge-radius ratio ($r/d$). For low Re, the flow is steady and symmetric with respect to the centreline. Over the first critical Reynolds number ($Re_{c1}$), the flow undergoes a Hopf bifurcation to a time-periodic flow, termed the primary instability. As Re further increases, the onset of the secondary instability of three-dimensional (3D) nature is detected beyond the second critical Reynolds number ($Re_{c2}$). Rounding the sharp edges of a square cylinder significantly affects the flow topology, leading to noticeable changes in both instabilities. By employing the Stuart–Landau equation, we investigated the criticality of the primary instability depending upon $r/d$. The onset of the 3D secondary instability was detected by using Floquet stability analysis. The temporal and spatial characteristics of the dominant modes (A, B, QP) were described. The neutral stability curves of each mode were computed depending upon $r/d$.

1990 ◽  
Vol 210 ◽  
pp. 537-563 ◽  
Author(s):  
O. John E. Matsson ◽  
P. Henrik Alfredsson

In a curved channel streamwise vortices, often called Dean vortices, may develop above a critical Reynolds number owing to centrifugal effects. Similar vortices can occur in a rotating plane channel due to Coriolis effects if the axis of rotation is normal to the mean flow velocity and parallel to the walls. In this paper the flow in a curved rotating channel is considered. It is shown from linear stability theory that there is a region for which centrifugal effects and Coriolis effects almost cancel each other, which increases the critical Reynolds number substantially. The flow visualization experiments carried out show that a complete cancellation of Dean vortices can be obtained for low Reynolds number. The rotation rate for which this occurs is in close agreement with predictions from linear stability theory. For curved channel flow a secondary instability of travelling wave type is found at a Reynolds number about three times higher than the critical one for the primary instability. It is shown that rotation can completely cancel the secondary instability.


2017 ◽  
Vol 821 ◽  
Author(s):  
M. H. Allouche ◽  
V. Botton ◽  
S. Millet ◽  
D. Henry ◽  
S. Dagois-Bohy ◽  
...  

The main objective of this work is to study experimentally the primary instability of non-Newtonian film flows down an inclined plane. We focus on low-concentration shear-thinning aqueous solutions obeying the Carreau law. The experimental study essentially consists of measuring wavelengths in marginal conditions, which yields the primary stability threshold for a given slope. The experimental results for neutral curves presented in the $(Re,f_{c})$ and $(Re,k)$ planes (where $f_{c}$ is the driving frequency, $k$ is the wavenumber and $Re$ is the Reynolds number) are in good agreement with the numerical results obtained by a resolution of the generalized Orr–Sommerfeld equation. The long-wave asymptotic extension of our results is consistent with former theoretical predictions of the critical Reynolds number. This is the first experimental evidence of the destabilizing effect of the shear-thinning behaviour in comparison with the Newtonian case: the critical Reynolds number is smaller, and the ratio between the critical wave celerity and the flow velocity at the free surface is larger.


1998 ◽  
Vol 357 ◽  
pp. 123-140 ◽  
Author(s):  
V. KUMARAN

The stability of the flow of a fluid in a flexible tube is analysed over a range of Reynolds numbers 1<Re<104 using a linear stability analysis. The system consists of a Hagen–Poiseuille flow of a Newtonian fluid of density ρ, viscosity η and maximum velocity V through a tube of radius R which is surrounded by an incompressible viscoelastic solid of density ρ, shear modulus G and viscosity ηs in the region R<r<HR. In the intermediate Reynolds number regime, the stability depends on the Reynolds number Re=ρVR/η, a dimensionless parameter [sum ]=ρGR2/η2, the ratio of viscosities ηr= ηs/η, the ratio of radii H and the wavenumber of the perturbations k. The neutral stability curves are obtained by numerical continuation using the analytical solutions obtained in the zero Reynolds number limit as the starting guess. For ηr=0, the flow becomes unstable when the Reynolds number exceeds a critical value Rec, and the critical Reynolds number increases with an increase in [sum ]. In the limit of high Reynolds number, it is found that Rec∝[sum ]α, where α varies between 0.7 and 0.75 for H between 1.1 and 10.0. An analysis of the flow structure indicates that the viscous stresses are confined to a boundary layer of thickness Re−1/3 for Re[Gt ]1, and the shear stress, scaled by ηV/R, increases as Re1/3. However, no simple scaling law is observed for the normal stress even at 103<Re<105, and consequently the critical Reynolds number also does not follow a simple scaling relation. The effect of variation of ηr on the stability is analysed, and it is found that a variation in ηr could qualitatively alter the stability characteristics. At relatively low values of [sum ] (about 102), the system could become unstable at all values of ηr, but at relatively high values of [sum ] (greater than about 104), an instability is observed only when the viscosity ratio is below a maximum value η*rm.


2014 ◽  
Vol 753 ◽  
pp. 131-149 ◽  
Author(s):  
Arghya Samanta

AbstractThe study of a film falling down an inclined plane is revisited in the presence of imposed shear stress. Earlier studies regarding this topic (Smith, J. Fluid Mech., vol. 217, 1990, pp. 469–485; Wei, Phys. Fluids, vol. 17, 2005a, 012103), developed on the basis of a low Reynolds number, are extended up to moderate values of the Reynolds number. The mechanism of the primary instability is provided under the framework of a two-wave structure, which is normally a combination of kinematic and dynamic waves. In general, the primary instability appears when the kinematic wave speed exceeds the speed of dynamic waves. An equality criterion between their speeds yields the neutral stability condition. Similarly, it is revealed that the nonlinear travelling wave solutions also depend on the kinematic and dynamic wave speeds, and an equality criterion between the speeds leads to an analytical expression for the speed of a family of travelling waves as a function of the Froude number. This new analytical result is compared with numerical prediction, and an excellent agreement is achieved. Direct numerical simulations of the low-dimensional model have been performed in order to analyse the spatiotemporal behaviour of nonlinear waves by applying a constant shear stress in the upstream and downstream directions. It is noticed that the presence of imposed shear stress in the upstream (downstream) direction makes the evolution of spatially growing waves weaker (stronger).


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.


This paper considers the hydrodynamic stability of flows in which there are no solid boundaries in the field of flow. The method used is an extension of that initiated by McKoen (1957), in which the fourth derivative, 0 iv , is assumed to be significant only near to the singular layer, but otherwise the complete fourth-order Orr—Sommerfeld equation is considered. An alternative derivation is given for McKoen’s integral form of the boundary condition for an antisymmetrical perturbation. In this integral it is necessary to approximate for (j) but not for any of its derivatives. It is shown that the present method will always lead to a neutral stability curve of wave number against Reynolds number, having two branches as R ->oo and hence a least critical R . The case of the plane laminar jet is considered, and a critical Reynolds number of 4 is obtained, which does not compare unreasonably with experiment in which unsteadiness is first detected at a Reynolds number of about 10. The lower branch of the neutral curve is found to be almost coincident with the R -axis.


1970 ◽  
Vol 43 (4) ◽  
pp. 813-818 ◽  
Author(s):  
M. D. J. Barry ◽  
M. A. S. Ross

Numerical analysis has been used to find the neutral stability curve for the flat plate boundary layer in zero pressure gradient when the main terms representing the growth of boundary-layer thickness are either included or excluded. The boundary layer is found to be slightly less stable when the extra terms are included. The calculations give a critical Reynolds number of 500.


2009 ◽  
Vol 636 ◽  
pp. 217-277 ◽  
Author(s):  
D. GIANNAKIS ◽  
R. ROSNER ◽  
P. F. FISCHER

We study the linear stability of the flow of a viscous electrically conducting capillary fluid on a planar fixed plate in the presence of gravity and a uniform magnetic field, assuming that the plate is either a perfect electrical insulator or a perfect conductor. We first confirm that the Squire transformation for magnetohydrodynamics is compatible with the stress and insulating boundary conditions at the free surface but argue that unless the flow is driven at fixed Galilei and capillary numbers, respectively parameterizing gravity and surface tension, the critical mode is not necessarily two-dimensional. We then investigate numerically how a flow-normal magnetic field and the associated Hartmann steady state affect the soft and hard instability modes of free-surface flow, working in the low-magnetic-Prandtl-number regime of conducting laboratory fluids (Pm ≤ 10−4). Because it is a critical-layer instability (moderately modified by the presence of the free surface), the hard mode exhibits similar behaviour as the even unstable mode in channel Hartmann flow, in terms of both the weak influence of Pm on its neutral-stability curve and the dependence of its critical Reynolds number Rec on the Hartmann number Ha. In contrast, the structure of the soft mode's growth-rate contours in the (Re, α) plane, where α is the wavenumber, differs markedly between problems with small, but non-zero, Pm and their counterparts in the inductionless limit, Pm ↘ 0. As derived from large-wavelength approximations and confirmed numerically, the soft mode's critical Reynolds number grows exponentially with Ha in inductionless problems. However, when Pm is non-zero the Lorentz force originating from the steady-state current leads to a modification of Rec(Ha) to either a sub-linearly increasing or a decreasing function of Ha, respectively for problems with insulating or perfectly conducting walls. In insulating-wall problems we also observe pairs of counter-propagating Alfvén waves, the upstream-propagating wave undergoing an instability driven by energy transferred from the steady-state shear to both of the velocity and magnetic degrees of freedom. Movies are available with the online version of the paper.


1989 ◽  
Vol 202 ◽  
pp. 543-557 ◽  
Author(s):  
P. Henrik Alfredsson ◽  
Håkan Persson

A flow visualization study of instabilities caused by Coriolis effects in plane rotating Poiseuille flow has been carried out. The primary instability takes the form of regularly spaced roll cells aligned in the flow direction. They may occur at Reynolds numbers as low as 100, i.e. almost two orders of magnitude lower than the critical Reynolds number for Tollmien-Schlichting waves in channel flow without rotation. The development of such roll cells was studied as a function of both the Reynolds number and the rotation rate and their properties compared with results from linear spatial stability theory. The theoretically obtained most unstable wavenumber agrees fairly well with the experimentally observed value. At high Reynolds number a secondary instability sets in, which is seen as a twisting of the roll cells. A wavytype disturbance is also seen at this stage which, if the rotational speed is increased, develops into large-scale ‘turbulence’ containing imbedded roll cells.


Sign in / Sign up

Export Citation Format

Share Document