Using numerical simulations, we investigate the near-field temporal variability of axisymmetric helium plumes as a function of inlet-based Richardson ( ${Ri}_0$ ) and Reynolds ( ${Re}_0$ ) numbers. Previous studies have shown that ${Ri}_0$ plays a leading-order role in determining the frequency at which large-scale vortices are produced (commonly called the ‘puffing’ frequency). By contrast, ${Re}_0$ dictates the strength of localized gradients, which are important during the transition from laminar to turbulent flow. The simulations presented here span a range of ${Ri}_0$ and ${Re}_0$ , and use adaptive mesh refinement to achieve high spatial resolutions. We find that as ${Re}_0$ increases for a given ${Ri}_0$ , the puffing motion undergoes a transition at a critical ${Re}_0$ , marking the onset of chaotic dynamics. Moreover, the critical ${Re}_0$ decreases as ${Ri}_0$ increases. When the puffing instability is non-chaotic, time series of different variables are well-correlated, exhibiting only modest changes in the dynamics (including period doubling and flapping). Once the flow becomes chaotic, denser ambient fluid penetrates the core of the plume, similar to penetrating ‘spikes’ formed by Rayleigh–Taylor instabilities, leading to only moderately correlated flow variables. These changes result in a non-trivial dependence of the puffing frequency on ${Re}_0$ . Specifically, at sufficiently low ${Re}_0$ , the puffing frequency falls below the prediction from Wimer et al. ( J. Fluid Mech. , vol. 895, 2020). As ${Re}_0$ increases beyond the critical ${Re}_0$ , the puffing frequency increases and then drops back down to the predicted scaling. The dependence of the puffing frequency on ${Re}_0$ provides a possible explanation for previously observed changes in the scaling of the puffing frequency for high ${Ri}_0$ .
more »
« less
Richardson and Reynolds number effects on the near field of buoyant plumes: flow statistics and fluxes
The near-field characteristics of highly buoyant plumes, commonly referred to as lazy plumes, remain relatively poorly understood across a range of flow conditions, particularly compared with our understanding of far-field characteristics. Here, we perform fully resolved three-dimensional numerical simulations of round helium plumes to characterize the effects of different inlet Richardson, ${Ri}_0$ , and Reynolds, ${Re}_0$ , numbers on first- and second-order statistical moments as well as average vertical fluxes in the near field. For sufficiently high ${Re}_0$ at a particular ${Ri}_0$ , heavy air can penetrate the core of the plume, reminiscent of spikes in the classical Rayleigh–Taylor instability. In the most turbulent simulation, this penetration becomes so strong that a recirculation zone forms along the centreline of the plume. Vertical fluxes are found to scale linearly with vertical distance from the plume inlet, consistent with experimental and numerical observations (Jiang & Luo, Flow Turbul. Combust. , vol. 64, 2000, pp. 43–69; Kaye & Hunt, Intl J. Heat Fluid Flow , vol. 30, 2009, pp. 1099–1105). We analytically derive this linear scaling from the governing equations by making a radial entrainment hypothesis whereby ambient fluid is entrained, on average, only in the radial direction at a finite distance from the inlet. Through this derivation, we identify physical mechanisms that can cause these relationships to remain only approximately valid for the present simulations. Lastly, we identify near-field power-law scaling relations for the flux magnitudes based on ${Ri}_0$ , and also examine vertical profiles of the non-dimensional Richardson number flux. Ultimately, insights from the present simulations are used to define near-, intermediate- and far-field regions in buoyant plumes.
more »
« less
- Award ID(s):
- 1847111
- PAR ID:
- 10425933
- Date Published:
- Journal Name:
- Journal of Fluid Mechanics
- Volume:
- 961
- ISSN:
- 0022-1120
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
We perform a direct numerical simulation (DNS) of interacting Kelvin–Helmholtz instabilities (KHI) that arise at a stratified shear layer where KH billow cores are misaligned or exhibit varying phases along their axes. Significant evidence of these dynamics in early laboratory shear-flow studies by Thorpe ( Geophys. Astrophys. Fluid Dyn. , vol. 34, 1985, pp. 175–199) and Thorpe ( J. Geophys. Res. , vol. 92, 1987, pp. 5231–5248), in observations of KH billow misalignments in tropospheric clouds (Thorpe, Q. J. R. Meteorol. Soc. , vol. 128, 2002, pp. 1529–1542) and in recent direct observations of such events in airglow and polar mesospheric cloud imaging in the upper mesosphere reveals that these dynamics are common. More importantly, the laboratory and mesospheric observations suggest that these dynamics lead to more rapid and more intense instabilities and turbulence than secondary convective instabilities in billow cores and secondary KHI in stratified braids between and around adjacent billows. To date, however, no simulations exploring the dynamics and energetics of interacting KH billows (apart from pairing) have been performed. Our DNS performed for Richardson number $Ri=0.10$ and Reynolds number $Re=5000$ demonstrates that KHI tubes and knots (i) comprise strong and complex vortex interactions accompanying misaligned KH billows, (ii) accelerate the transition to turbulence relative to secondary instabilities of individual KH billows, (iii) yield significantly stronger turbulence than secondary KHI in billow braids and secondary convective instabilities in KHI billow cores and (iv) expand the suite of secondary instabilities previously recognized to contribute to KHI dynamics and breakdown to turbulence in realistic geophysical environments.more » « less
-
We perform a direct numerical simulation (DNS) of interacting Kelvin–Helmholtz instabilities (KHI) that arise at a stratified shear layer where KH billow cores are misaligned or exhibit varying phases along their axes. Significant evidence of these dynamics in early laboratory shear-flow studies by Thorpe (Geophys. Astrophys. Fluid Dyn., vol. 34, 1985, pp. 175–199) and Thorpe (J. Geophys. Res., vol. 92, 1987, pp. 5231–5248), in observations of KH billow misalignments in tropospheric clouds (Thorpe, Q. J. R. Meteorol. Soc., vol. 128, 2002, pp. 1529–1542) and in recent direct observations of such events in airglow and polar mesospheric cloud imaging in the upper mesosphere reveals that these dynamics are common. More importantly, the laboratory and mesospheric observations suggest that these dynamics lead to more rapid and more intense instabilities and turbulence than secondary convective instabilities in billow cores and secondary KHI in stratified braids between and around adjacent billows. To date, however, no simulations exploring the dynamics and energetics of interacting KH billows (apart from pairing) have been performed. Our DNS performed for Richardson number Ri = 0.10 and Reynolds number Re = 5000 demonstrates that KHI tubes and knots (i) comprise strong and complex vortex interactions accompanying misaligned KH billows, (ii) accelerate the transition to turbulence relative to secondary instabilities of individual KH billows, (iii) yield significantly stronger turbulence than secondary KHI in billow braids and secondary convective instabilities in KHI billow cores and (iv) expand the suite of secondary instabilities previously recognized to contribute to KHI dynamics and breakdown to turbulence in realistic geophysical environments.more » « less
-
more » « less
This study presents observations of a buoyant plume off Winyah Bay, South Carolina, which was formed under conditions of high freshwater discharge and upwelling-favorable wind forcing. Analysis of observations demonstrates that the response of the anticyclonic bulge formed by tidally modulated estuarine outflow to the light upwelling-favorable wind is more complex than the previously studied far-field response. The latter can be described by a slab-like model with mixing concentrating at the offshore edge of a buoyant layer. The observed plume depth increased from ~3 m near the mouth to 6 m at the offshore edge, with plume depth changing in a steplike fashion rather than continuously. CTD profiles near these steps revealed overturning indicative of vigorous mixing. Estimates of the gradient Richardson number confirmed the likelihood of mixing/entrainment not only at the offshore edge of the plume but also in the proximity of the observed steps. We hypothesize that these steps represent tidal fronts that undergo geostrophic adjustment and are advected offshore by the superimposed Ekman drift. Scaling analysis suggests that mixing and entrainment at the observed interior fronts can be enhanced by superposition of geostrophic and wind-induced shear.
-
We numerically and theoretically investigate the Boussinesq Eady model, where a rapidly rotating density-stratified layer of fluid is subject to a meridional temperature gradient in thermal wind balance with a uniform vertically sheared zonal flow. Through a suite of numerical simulations, we show that the transport properties of the resulting turbulent flow are governed by quasigeostrophic (QG) dynamics in the rapidly rotating strongly stratified regime. The ‘vortex gas’ scaling predictions put forward in the context of the two-layer QG model carry over to this fully three-dimensional system: the functional dependence of the meridional flux on the control parameters is the same, the two adjustable parameters entering the theory taking slightly different values. In line with the QG prediction, the meridional heat flux is depth-independent. The vertical heat flux is such that turbulence transports buoyancy along isopycnals, except in narrow layers near the top and bottom boundaries, the thickness of which decreases as the diffusivities go to zero. The emergent (re)stratification is set by a simple balance between the vertical heat flux and diffusion along the vertical direction. Overall, this study demonstrates how the vortex-gas scaling theory can be adapted to quantitatively predict the magnitude and vertical structure of the meridional and vertical heat fluxes, and of the emergent stratification, without additional fitting parameters.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1017/jfm.2023.208
