scholarly journals Global stability of buoyant jets and plumes

2017 ◽  
Vol 835 ◽  
pp. 654-673 ◽  
Author(s):  
R. V. K. Chakravarthy ◽  
L. Lesshafft ◽  
P. Huerre

The linear global stability of laminar buoyant jets and plumes is investigated under the low-Mach-number approximation. For Richardson numbers in the range $10^{-4}\leqslant Ri\leqslant 10^{3}$ and density ratios $S=\unicode[STIX]{x1D70C}_{\infty }/\unicode[STIX]{x1D70C}_{jet}$ between 1.05 and 7, only axisymmetric perturbations are found to exhibit global instability, consistent with experimental observations in helium jets. By varying the Richardson number over seven decades, the effects of buoyancy on the base flow and on the instability dynamics are characterised, and distinct behaviour is observed in the low-$Ri$ (jet) and in the high-$Ri$ (plume) regimes. A sensitivity analysis indicates that different physical mechanisms are responsible for the global instability dynamics in both regimes. In buoyant jets at low Richardson number, the baroclinic torque enhances the basic shear instability, whereas buoyancy provides the dominant instability mechanism in plumes at high Richardson number. The onset of axisymmetric global instability in both regimes is consistent with the presence of absolute instability. While absolute instability also occurs for helical perturbations, it appears to be too weak or too localised to give rise to a global instability.

2001 ◽  
Vol 428 ◽  
pp. 349-386 ◽  
Author(s):  
E. J. STRANG ◽  
H. J. S. FERNANDO

The results of a laboratory experiment designed to study turbulent entrainment at sheared density interfaces are described. A stratified shear layer, across which a velocity difference ΔU and buoyancy difference Δb is imposed, separates a lighter upper turbulent layer of depth D from a quiescent, deep lower layer which is either homogeneous (two-layer case) or linearly stratified with a buoyancy frequency N (linearly stratified case). In the parameter ranges investigated the flow is mainly determined by two parameters: the bulk Richardson number RiB = ΔbD/ΔU2 and the frequency ratio fN = ND=ΔU.When RiB > 1.5, there is a growing significance of buoyancy effects upon the entrainment process; it is observed that interfacial instabilities locally mix heavy and light fluid layers, and thus facilitate the less energetic mixed-layer turbulent eddies in scouring the interface and lifting partially mixed fluid. The nature of the instability is dependent on RiB, or a related parameter, the local gradient Richardson number Rig = N2L/ (∂u/∂z)2, where NL is the local buoyancy frequency, u is the local streamwise velocity and z is the vertical coordinate. The transition from the Kelvin–Helmholtz (K-H) instability dominated regime to a second shear instability, namely growing Hölmböe waves, occurs through a transitional regime 3.2 < RiB < 5.8. The K-H activity completely subsided beyond RiB ∼ 5 or Rig ∼ 1. The transition period 3.2 < RiB < 5 was characterized by the presence of both K-H billows and wave-like features, interacting with each other while breaking and causing intense mixing. The flux Richardson number Rif or the mixing efficiency peaked during this transition period, with a maximum of Rif ∼ 0.4 at RiB ∼ 5 or Rig ∼ 1. The interface at 5 < RiB < 5.8 was dominated by ‘asymmetric’ interfacial waves, which gradually transitioned to (symmetric) Hölmböe waves at RiB > 5:8.Laser-induced fluorescence measurements of both the interfacial buoyancy flux and the entrainment rate showed a large disparity (as large as 50%) between the two-layer and the linearly stratified cases in the range 1.5 < RiB < 5. In particular, the buoyancy flux (and the entrainment rate) was higher when internal waves were not permitted to propagate into the deep layer, in which case more energy was available for interfacial mixing. When the lower layer was linearly stratified, the internal waves appeared to be excited by an ‘interfacial swelling’ phenomenon, characterized by the recurrence of groups or packets of K-H billows, their degeneration into turbulence and subsequent mixing, interfacial thickening and scouring of the thickened interface by turbulent eddies.Estimation of the turbulent kinetic energy (TKE) budget in the interfacial zone for the two-layer case based on the parameter α, where α = (−B + ε)/P, indicated an approximate balance (α ∼ 1) between the shear production P, buoyancy flux B and the dissipation rate ε, except in the range RiB < 5 where K-H driven mixing was active.


2018 ◽  
Vol 843 ◽  
pp. 575-600 ◽  
Author(s):  
Jean-Philippe Matas ◽  
Antoine Delon ◽  
Alain Cartellier

We study the destabilization of a round liquid jet by a fast annular gas stream. We measure the frequency of the shear instability waves for several geometries and air/water velocities. We then carry out a linear stability analysis, and show that there are three competing mechanisms for the destabilization: a convective instability, an absolute instability driven by surface tension and an absolute instability driven by confinement. We compare the predictions of this analysis with experimental results, and propose scaling laws for wave frequency in each regime. We finally introduce criteria to predict the boundaries between these three regimes.


2011 ◽  
Vol 685 ◽  
pp. 54-82 ◽  
Author(s):  
Hieu T. Pham ◽  
Sutanu Sarkar

AbstractThe fine-scale response of a subsurface stable stratified jet subject to the forcing of surface wind stress and surface cooling is investigated using direct numerical simulation. The initial velocity profile consists of a symmetric jet located below a surface layer driven by a constant wind stress. The initial density profile is well-mixed in the surface layer and linearly stratified in both upper and lower flanks of the jet. The minimum value of the gradient Richardson number in the upper flank of the jet exceeds the critical value of 0.25 for linear shear instability. Broadband finite-amplitude fluctuations are introduced to the surface layer to initiate the simulation. Turbulence is generated in the surface layer and deepens into the jet upper flank. Internal waves generated by the turbulent surface layer are observed to propagate downward across the jet. The momentum flux carried by the waves is significantly smaller than the Reynolds shear stress extracted from the background velocity. The wave energy flux is also smaller than the turbulence production by mean shear. Ejections of fluid parcels by horseshoe-like vortices cause intermittent patches of intense dissipation inside the jet upper flank where the background gradient Richardson number is larger than 0.25. Drag due to the wind stress is smaller than the drag caused by turbulent stress in the flow. Analysis of the mean and turbulent kinetic energy budgets suggests that the energy input by surface forcing is considerably smaller than the energy extracted from the initially imposed background shear in the surface layer.


2009 ◽  
Vol 622 ◽  
pp. 291-320 ◽  
Author(s):  
M. D. GRIFFITH ◽  
T. LEWEKE ◽  
M. C. THOMPSON ◽  
K. HOURIGAN

Pulsatile inlet flow through a circular tube with an axisymmetric blockage of varying size is studied both numerically and experimentally. The geometry consists of a long, straight tube and a blockage, semicircular in cross-section, serving as a simplified model of an arterial stenosis. The stenosis is characterized by a single parameter, the aim being to highlight fundamental behaviours of constricted pulsatile flows. The Reynolds number is varied between 50 and 700 and the stenosis degree by area between 0.20 and 0.90. Numerically, a spectral element code is used to obtain the axisymmetric base flow fields, while experimentally, results are obtained for a similar set of geometries, using water as the working fluid. For low Reynolds numbers, the flow is characterized by a vortex ring which forms directly downstream of the stenosis, for which the strength and downstream propagation velocity vary with the stenosis degree. Linear stability analysis is performed on the simulated axisymmetric base flows, revealing a range of absolute instability modes. Comparisons are drawn between the numerical linear stability analysis and the observed instability in the experimental flows. The observed flows are less stable than the numerical analysis predicts, with convective shear layer instability present in the experimental flows. Evidence is found of Kelvin–Helmholtz-type shear layer roll-ups; nonetheless, the possibility of the numerically predicted absolute instability modes acting in the experimental flow is left open.


2013 ◽  
Vol 43 (11) ◽  
pp. 2475-2489 ◽  
Author(s):  
Oliver M. Sun ◽  
Steven R. Jayne ◽  
Kurt L. Polzin ◽  
Bryan A. Rahter ◽  
Louis C. St. Laurent

Abstract Data from three midlatitude, month-long surveys are examined for evidence of enhanced vertical mixing associated with the transition layer (TL), here defined as the strongly stratified layer that exists between the well mixed layer and the thermocline below. In each survey, microstructure estimates of turbulent dissipation were collected concurrently with fine-structure stratification and shear. Survey-wide averages are formed in a “TL coordinate” zTL, which is referenced around the depth of maximum stratification for each profile. Averaged profiles show characteristic TL structures such as peaks in stratification N2 and shear variance S2, which fall off steeply above zTL = 0 and more gradually below. Turbulent dissipation rates ɛ are 5–10 times larger than those found in the upper thermocline (TC). The gradient Richardson number Ri = N2/S2 becomes unstable (Ri &lt; 0.25) within ~10 m of the TL upper boundary, suggesting that shear instability is active in the TL for zTL &gt; 0. Ri is stable for zTL ≤ 0. Turbulent dissipation is found to scale exponentially with depth for zTL ≤ 0, but the decay scales are different for the TL and upper TC: ɛ scales well with either N2 or S2. Owing to the strong correlation between S2 and N2, existing TC scalings of the form ɛ ~ |S|p|N|q overpredict variations in ɛ. The scale dependence of shear variance is not found to significantly affect the scalings of ɛ versus N2 and S2 for zTL ≤ 0. However, the onset of unstable Ri at the top of the TL is sensitively dependent to the resolution of the shears.


2016 ◽  
Vol 801 ◽  
pp. 508-553 ◽  
Author(s):  
Eunok Yim ◽  
Paul Billant ◽  
Claire Ménesguen

This paper investigates the stability of an axisymmetric pancake vortex with Gaussian angular velocity in radial and vertical directions in a continuously stratified-rotating fluid. The different instabilities are determined as a function of the Rossby number $Ro$, Froude number $F_{h}$, Reynolds number $Re$ and aspect ratio ${\it\alpha}$. Centrifugal instability is not significantly different from the case of a columnar vortex due to its short-wavelength nature: it is dominant when the absolute Rossby number $|Ro|$ is large and is stabilized for small and moderate $|Ro|$ when the generalized Rayleigh discriminant is positive everywhere. The Gent–McWilliams instability, also known as internal instability, is then dominant for the azimuthal wavenumber $m=1$ when the Burger number $Bu={\it\alpha}^{2}Ro^{2}/(4F_{h}^{2})$ is larger than unity. When $Bu\lesssim 0.7Ro+0.1$, the Gent–McWilliams instability changes into a mixed baroclinic–Gent–McWilliams instability. Shear instability for $m=2$ exists when $F_{h}/{\it\alpha}$ is below a threshold depending on $Ro$. This condition is shown to come from confinement effects along the vertical. Shear instability transforms into a mixed baroclinic–shear instability for small $Bu$. The main energy source for both baroclinic–shear and baroclinic–Gent–McWilliams instabilities is the potential energy of the base flow instead of the kinetic energy for shear and Gent–McWilliams instabilities. The growth rates of these four instabilities depend mostly on $F_{h}/{\it\alpha}$ and $Ro$. Baroclinic instability develops when $F_{h}/{\it\alpha}|1+1/Ro|\gtrsim 1.46$ in qualitative agreement with the analytical predictions for a bounded vortex with angular velocity slowly varying along the vertical.


2008 ◽  
Vol 612 ◽  
pp. 237-260 ◽  
Author(s):  
JÉRÔME FONTANE ◽  
LAURENT JOLY

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.


2013 ◽  
Vol 722 ◽  
pp. 618-644 ◽  
Author(s):  
G. Boyer ◽  
G. Casalis ◽  
J. L. Estivalèzes

AbstractThe present article aims at enhancing the computation of the global stability modes of the internal flow of solid rocket motors (SRMs) approximated by the Taylor–Culick solution. This modal approach suffers from the consequences of the non-normality of the global linearized incompressible Navier–Stokes operator, namely the lack of robustness of the eigenvalues that can lead to the computation of pseudo-modes rather than actual eigenmodes. In this respect, the effects of non-normality associated with strongly amplified eigenfunctions are highlighted on a simplified convective–diffusive stability problem with uniformly accelerated base state, the latter property being a typical characteristic of the Taylor–Culick flow. Non-convergence zones for the eigenvalues are exhibited for large Reynolds numbers and are related to the critical sensitivity to disturbances applied to one of the boundary conditions. For this reason, and according to experimental and numerical data related to the stability of simplified SRMs, a global stability analysis is performed assuming that the hydrodynamic fluctuations emerge from a geometrical defect applied at the sidewall. This comes to fix the upstream boundary condition at the abscissa of the sidewall disturbance. The resulting eigenmodes are shown to be discrete, numerically converged, well identified by a finite number of points of undefined phase of the velocity fluctuations. They marginally depend on Reynolds number variations, but are modified by changes on the boundaries location. As in the simplified problem, the inflow boundary condition is the most critical in terms of sensitivity to numerical errors, although not dramatic. Finally, the sensitivity analysis to infinitesimal base flow changes indicates that the variations applied close to the inflow boundary condition induce the largest move of the eigenvalues. In spite of the large non-normal effects induced by the large polynomial growth of the eigenfunctions, this paper shows that discrete instabilities may emerge from a wall defect, in contrast to configurations without such a geometrical perturbation whose dynamics may be rather driven by pseudo-modes.


2010 ◽  
Vol 660 ◽  
pp. 499-526 ◽  
Author(s):  
P. MELIGA ◽  
D. SIPP ◽  
J.-M. CHOMAZ

We study the linear dynamics of global eigenmodes in compressible axisymmetric wake flows, up to the high subsonic regime. We consider both an afterbody flow at zero angle of attack and a sphere, and find that the sequence of bifurcations destabilizing the axisymmetric steady flow is independent of the Mach number and reminiscent of that documented in the incompressible wake past a sphere and a disk (Natarajan & Acrivos, J. Fluid Mech., vol. 254, 1993, p. 323), hence suggesting that the onset of unsteadiness in this class of flows results from a global instability. We determine the boundary separating the stable and unstable domains in the (M, Re) plane, and show that an increase in the Mach number yields a stabilization of the afterbody flow, but a destabilization of the sphere flow. These compressible effects are further investigated by means of adjoint-based sensitivity analyses relying on the computation of gradients or sensitivity functions. Using this theoretical formalism, we show that they do not act through specific compressibility effects at the disturbance level but mainly through implicit base flow modifications, an effect that had not been taken into consideration by previous studies based on prescribed parallel base flow profiles. We propose a physical interpretation for the observed compressible effects, based on the competition between advection and production of disturbances, and provide evidence linking the stabilizing/destabilizing effect observed when varying the Mach number to a strengthening/weakening of the disturbance advection mechanism. We show, in particular, that the destabilizing effect of compressibility observed in the case of the sphere results from a significant increase of the backflow velocity in the whole recirculating bubble, which opposes the downstream advection of disturbances.


Author(s):  
Sukanta Bhattacharjee ◽  
Sumanta Acharya

A numerical study is conducted to understand the global instability of very low-density jets (as encountered in thermal plasmas). The simulations have been carried out for different density ratios, S = ρj /ρ∞, ranging from 0.5 to 0.01, and different Reynolds numbers. The results show that the global instability mode exists in the range of density ratios investigated. A strategy for reducing the instability by altering the density profile of the surrounding gas stream is explored. Specifically, a ramp density profile or a step shape with an outward offset was studied, and it was observed that there was a reduction in the instability amplitude with the modified density profiles. Such a lowering in the instability fluctuations can be beneficial in stabilizing the thermal plasma behavior.


Sign in / Sign up

Export Citation Format

Share Document