skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Turbulence, entrainment and low-order description of a transitional variable-density jet
Geophysical flows occur over a large range of scales, with Reynolds numbers and Richardson numbers varying over several orders of magnitude. For this study, jets of different densities were ejected vertically into a large ambient region, considering conditions relevant to some geophysical phenomena. Using particle image velocimetry, the velocity fields were measured for three different gases exhausting into air – specifically helium, air and argon. Measurements focused on both the jet core and the entrained ambient. Experiments considered relatively low Reynolds numbers from approximately 1500 to 10 000 with Richardson numbers near 0.001 in magnitude. These included a variety of flow responses, notably a nearly laminar jet, turbulent jets and a transitioning jet in between. Several features were studied, including the jet development, the local entrainment ratio, the turbulent Reynolds stresses and the eddy strength. Compared to a fully turbulent jet, the transitioning jet showed up to 50 % higher local entrainment and more significant turbulent fluctuations. For this condition, the eddies were non-axisymmetric and larger than the exit radius. For turbulent jets, the eddies were initially smaller and axisymmetric while growing with the shear layer. At lower turbulent Reynolds number, the turbulent stresses were more than 50 % higher than at higher turbulent Reynolds number. In either case, the low-density jet developed faster than a comparable non-buoyant jet. Quadrant analysis and proper orthogonal decomposition were also utilized for insight into the entrainment of the jet, as well as to assess the energy distribution with respect to the number of eigenmodes. Reynolds shear stresses were dominant in Q1 and Q3 and exhibited negligible contributions from the remaining two quadrants. Both analysis techniques showed that the development of stresses downstream was dependent on the Reynolds number while the spanwise location of the stresses depended on the Richardson number.  more » « less
Award ID(s):
1346580
PAR ID:
10414303
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
836
ISSN:
0022-1120
Page Range / eLocation ID:
1009 to 1049
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The interplay between viscoelasticity and inertia in dilute polymer solutions at high deformation rates can result in inertioelastic instabilities. The nonlinear evolution of these instabilities generates a state of turbulence with significantly different spatiotemporal features compared to Newtonian turbulence, termed elastoinertial turbulence (EIT). We ex- plore EIT by studying the dynamics of a submerged planar jet of a dilute aqueous polymer solution injected into a quiescent tank of water using a combination of schlieren imaging and laser Doppler velocimetry (LDV). We show how fluid elasticity has a nonmonotonic effect on the jet stability depending on its magnitude, creating two distinct regimes in which elastic effects can either destabilize or stabilize the jet. In agreement with linear stability analyses of viscoelastic jets, an inertioelastic shear-layer instability emerges near the edge of the jet for small levels of elasticity, independent of bulk undulations in the fluid column. The growth of this disturbance mode destabilizes the flow, resulting in a turbulence transition at lower Reynolds numbers and closer to the nozzle compared to the conditions required for the transition to turbulence in a Newtonian jet. Increasing the fluid elasticity merges the shear-layer instability into a bulk instability of the jet column. In this regime, elastic tensile stresses generated in the shear layer act as an “elastic membrane” that partially stabilizes the flow, retarding the transition to turbulence to higher levels of inertia and greater distances from the nozzle. In the fully turbulent state far from the nozzle, planar viscoelastic jets exhibit unique spatiotemporal features associated with EIT. The time-averaged angle of jet spreading, an Eulerian measure of the degree of entrainment, and the centerline velocity of the jets both evolve self-similarly with distance from the nozzle. The autocovariance of the schlieren images in the fully turbulent region of the jets shows coherent structures that are elongated in the streamwise direction, consistent with the suppression of streamwise vortices by elastic stresses. These coherent structures give a higher spectral energy to small frequency modes in EIT characterized by LDV measurements of the velocity fluctuations at the jet centerline. Finally, our LDV measurements reveal a frequency spectrum characterized by a −3 power-law exponent, different from the well-known −5/3 power-law exponent characteristic of Newtonian turbulence. 
    more » « less
  2. Direct Numerical Simulation (DNS) of compressible spatially-developing turbulent boundary layers (SDTBL) is performed at a Mach number of 2.5 and low/high Reynolds numbers over isothermal Zero-Pressure Gradient (ZPG) flat plates. Turbulent inflow information is generated via a dynamic rescaling-recycling approach (J. Fluid Mech., 670, pp. 581-605, 2011), which avoids the use of empirical correlations in the computation of inlet turbulent scales. The range of the low Reynolds number case is approximately 400-800, based on the momentum thickness, freestream velocity and wall viscosity. DNS at higher Reynolds numbers (~3,000, about four-fold larger) is also carried out with the purpose of analyzing the effect of Reynolds number on the transport phenomena in the supersonic regime. Additionally, low/high order flow statistics are compared with DNS of an incompressible isothermal ZPG boundary layer at similar low Reynolds numbers and the temperature regarded as a passive scalar. Peaks of turbulence intensities move closer to the wall as the Reynolds number increases in the supersonic flat plate. Furthermore, Reynolds shear stresses depict a much larger "plateau" (constant shear layer) at the highest Reynolds number considered in present study. 
    more » « less
  3. Theoretical predictions and numerical simulations are used to determine the transition to bubble and conical vortex breakdown in low-Mach-number laminar axisymmetric variable-density swirling jets. A critical value of the swirl number $$S$$ for the onset of the bubble ( $$S^*_B$$ ) and the cone ( $$S^*_C$$ ) is determined as the jet-to-ambient density ratio $$\varLambda$$ is varied, with the temperature dependence of the gas density and viscosity appropriate to that of air. The criterion of failure of the slender quasi-cylindrical approximation predicts $$S^*_B$$ that decreases with increasing values of $$\varLambda$$ for a jet in solid-body rotation emerging sharply into a quiescent atmosphere. In addition, a new criterion for the onset of conical breakdown is derived from divergence of the initial value of the radial spreading rate of the jet occurring at $$S^*_C$$ , found to be independent of $$\varLambda$$ , in an asymptotic analysis for small distances from the inlet plane. To maintain stable flow in the unsteady numerical simulations, an effective Reynolds number $$Re_{eff}$$ , defined employing the geometric mean of the viscosity in the jet and ambient atmosphere, is fixed at $$Re_{eff}=200$$ for all $$\varLambda$$ . Similar to the theoretical predictions, numerical calculations of $$S^*_B$$ decrease monotonically as $$\varLambda$$ is increased. The critical swirl numbers for the cone, $$S^*_C$$ , are found to depend strongly on viscous effects; for $$\varLambda =1/5$$ , the low jet Reynolds number (51) at $$Re_{eff}=200$$ delays the transition to the cone, while for $$\varLambda =5$$ at $$Re_{eff}=200$$ , the large increase in kinematic viscosity in the external fluid produces a similar trend, significantly increasing $$S^*_C$$ . 
    more » « less
  4. Liquid jets in surrounding air face capillary and shear forces which eventually disintegrate the jet into droplets or spray. The instabilities developed in the flow inevitably break down an initial laminar (coherent) jet into a turbulent one. In the manufacturing process called grinding, one of the oldest approaches of shaping metals and other materials, liquid coolant jets are frequently used. A non-coherent or turbulent jet has a reduced flow rate due to cavitation, air entrapment and atomization of the fluid particles. The jet spread does not allow the coolant jet to effectively breach the high-speed rotating air layer, created by entrainment of air along the surface of rapidly rotating grinding wheel. The coherent, nearly columnal jet should be sufficiently long to maintain its initial velocity to penetrate the layer of air rotating with the grinding wheel. Thus, in many critical grinding applications, it is advised to use a coherent jet instead of a spray to eradicate defects of ground surface. In this study, we present simulations of liquid jet flows to see how the jet develops and breaks due to surface tension and shear forces. Creating an accurate model to predict liquid jet characteristics, especially for high-speed applications such as grinding wheel cooling would require wellresolving numerical grids and turbulence model selection. The problem being multi-phased with a density ratio of coolant-to-air being order of 1000 adds to the computational complexity. The presented numerical model and results are different compared to the previous simulations of liquid jets as the characteristics of jet disintegration are explored under conditions that closely resemble a grinding cooling application. Finite volume discretization of the flow domain and calculation of flow field characteristics were done by commercial software ANSYS Mesh and ANSYS Fluent modules, respectively. The numerical calculation and visualization of disintegration of free jet and the jet impinging into grinding wheel will be presented. 
    more » « less
  5. Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations. 
    more » « less