scholarly journals Dynamics of a Turbulent Buoyant Plume in a Stratified Fluid: An Idealized Model of Subglacial Discharge in Greenland Fjords

2017 ◽  
Vol 47 (10) ◽  
pp. 2611-2630 ◽  
Author(s):  
Ekaterina Ezhova ◽  
Claudia Cenedese ◽  
Luca Brandt

AbstractThis study reports the results of large-eddy simulations of an axisymmetric turbulent buoyant plume in a stratified fluid. The configuration used is an idealized model of the plume generated by a subglacial discharge at the base of a tidewater glacier with an ambient stratification typical of Greenland fjords. The plume is discharged from a round source of various diameters and characteristic stratifications for summer and winter are considered. The classical theory for the integral parameters of a turbulent plume in a homogeneous fluid gives accurate predictions in the weakly stratified lower layer up to the pycnocline, and the plume dynamics are not sensitive to changes in the source diameter. In winter, when the stratification is similar to an idealized two-layer case, turbulent entrainment and generation of internal waves by the plume top are in agreement with the theoretical and numerical results obtained for turbulent jets in a two-layer stratification. In summer, instead, the stratification is more complex and turbulent entrainment by the plume top is significantly reduced. The subsurface layer in summer is characterized by a strong density gradient and the oscillating plume generates internal waves that might serve as an indicator of submerged plumes not penetrating to the surface.

1993 ◽  
Vol 254 ◽  
pp. 113-126 ◽  
Author(s):  
Nikolay Kuznetsov

A horizontal channel of infinite length and depth and of constant width contains inviscid, incompressible, two-layer fluid under gravity. The upper layer has constant finite depth and is occupied by a fluid of constant density ρ*. The lower layer has infinite depth and is occupied by a fluid of constant density ρ > ρ*. The parameter ε = (ρ/ρ*)–1 is assumed to be small. The lower fluid is bounded internally by an immersed horizontal cylinder which extends right across the channel and has its generators normal to the sidewalls. The free, time-harmonic oscillations of fluid, which have finite kinetic and potential energy (such oscillations are called trapped modes), are investigated. Trapped modes in homogeneous fluid above submerged cylinders and other obstacles are well known. In the present paper it is shown that there are two sets of frequencies of trapped modes for the two-layer fluid. The frequencies of the first finite set are close to the frequencies of trapped modes in the homogeneous fluid (when ρ* = ρ). They correspond to the trapped modes of waves on the free surface of the upper fluid. The frequencies of the second finite set are proportional to ε, and hence, are small. These latter frequencies correspond to the trapped modes of internal waves on the interface between two fluids. To obtain these results the perturbation method for a quadratic operator family was applied. The quadratic operator family with bounded, symmetric, linear, integral operators in the space L2(−∞, +∞) arises as a result of two reductions of the original problem. The first reduction allows to consider the potential in the lower fluid only. The second reduction is the same as used by Ursell (1987).


2016 ◽  
Vol 46 (11) ◽  
pp. 3415-3437 ◽  
Author(s):  
Ekaterina Ezhova ◽  
Claudia Cenedese ◽  
Luca Brandt

AbstractThe behavior of an axisymmetric vertical turbulent jet in an unconfined stratified environment is studied by means of well-resolved, large-eddy simulations. The stratification is two uniform layers separated by a thermocline. This study considers two cases: when the thermocline thickness is small and on the order of the jet diameter at the thermocline entrance. The Froude number of the jet at the thermocline varies from 0.6 to 1.9, corresponding to the class of weak fountains. The mean jet penetration, stratified turbulent entrainment, jet oscillations, and the generation of internal waves are examined. The mean jet penetration is predicted well by a simple model based on the conservation of the source energy in the thermocline. The entrainment coefficient for the thin thermocline is consistent with the theoretical model for a two-layer stratification with a sharp interface, while for the thick thermocline entrainment is larger at low Froude numbers. The data reveal the presence of a secondary horizontal flow in the upper part of the thick thermocline, resulting in the entrainment of fluid from the thermocline rather than from the upper stratification layer. The spectra of the jet oscillations in the thermocline display two peaks, at the same frequencies for both stratifications at fixed Froude number. For the thick thermocline, internal waves are generated only at the lower frequency, since the higher peak exceeds the maximal buoyancy frequency. For the thin thermocline, conversely, the spectra of the internal waves show the two peaks at low Froude numbers, whereas only one peak at the lower frequency is observed at higher Froude numbers.


2014 ◽  
Vol 747 ◽  
pp. 605-634 ◽  
Author(s):  
G. M. Reznik

AbstractWe examine nonlinear geostrophic adjustment in a rapidly rotating (small Rossby number Ro) stably neutrally stratified (SNS) fluid consisting of a stratified upper layer with $N>f$ ($N$ is the buoyancy frequency, $f$ the Coriolis parameter) and a homogeneous lower layer, the density and other fields being continuous at the interface between the layers. The angular speed of rotation is non-parallel to gravity; the traditional and hydrostatic approximations are not used. The wave spectrum in the model consists of internal and gyroscopic waves. During the adjustment an arbitrary long-wave perturbation is split in a unique way into slow quasi-geostrophic (QG) and fast ageostrophic components with typical time scales $(Ro\, f)^{-1}$ and $f^{-1}$, respectively. The QG flow is governed by two coupled nonlinear equations of conservation of QG potential vorticity (PV) in the layers. The fast component is a sum of internal waves and inertial oscillations (long gyroscopic waves) confined to the homogeneous layer and modulated by an amplitude depending on coordinates and slow time. On times $t\sim (\, f\, Ro)^{-1}$ the slow component is not influenced by the fast one but the inertial oscillations amplitude is coupled to the QG flow and obeys an equation practically coinciding with that in the barotropic case (Reznik, J. Fluid Mech., vol. 743, 2014, pp. 585–605). A non-stationary boundary layer with large vertical gradients of horizontal velocity develops in the stratified layer near the interface to prevent penetration of the inertial oscillations into the stratified fluid; an analogous weaker boundary layer arises near the upper rigid lid. At large times the internal waves gradually decay because of dispersion and the resulting motion consists of the slow QG component and inertial oscillations confined to the barotropic lower layer.


AIAA Journal ◽  
2001 ◽  
Vol 39 ◽  
pp. 1469-1477
Author(s):  
W. Zhao ◽  
S. H. Frankel ◽  
L. Mongeau

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.


2008 ◽  
Vol 616 ◽  
pp. 327-356 ◽  
Author(s):  
BRIAN L. WHITE ◽  
KARL R. HELFRICH

A steady theory is presented for gravity currents propagating with constant speed into a stratified fluid with a general density profile. Solution curves for front speed versus height have an energy-conserving upper bound (the conjugate state) and a lower bound marked by the onset of upstream influence. The conjugate state is the largest-amplitude nonlinear internal wave supported by the ambient stratification, and in the limit of weak stratification approaches Benjamin's energy-conserving gravity current solution. When the front speed becomes critical with respect to linear long waves generated above the current, steady solutions cannot be calculated, implying upstream influence. For non-uniform stratification, the critical long-wave speed exceeds the ambient long-wave speed, and the critical-Froude-number condition appropriate for uniform stratification must be generalized. The theoretical results demonstrate a clear connection between internal waves and gravity currents. The steady theory is also compared with non-hydrostatic numerical solutions of the full lock release initial-value problem. Some solutions resemble classic gravity currents with no upstream disturbance, but others show long internal waves propagating ahead of the gravity current. Wave generation generally occurs when the stratification and current speed are such that the steady gravity current theory fails. Thus the steady theory is consistent with the occurrence of either wave-generating or steady gravity solutions to the dam-break problem. When the available potential energy of the dam is large enough, the numerical simulations approach the energy-conserving conjugate state. Existing laboratory experiments for intrusions and gravity currents produced by full-depth lock exchange flows over a range of stratification profiles show excellent agreement with the conjugate state solutions.


Author(s):  
Tobias Kukulka ◽  
Todd Thoman

AbstractDispersion processes in the ocean surface boundary layer (OSBL) determine marine material distributions such as those of plankton and pollutants. Sheared velocities drive shear dispersion, which is traditionally assumed to be due to mean horizontal currents that decrease from the surface. However, OSBL turbulence supports along-wind jets; located in near-surface convergence and downwelling regions, such turbulent jets contain strong local shear. Through wind-driven idealized and large eddy simulation (LES) models of the OSBL, this study examines the role of turbulent along-wind jets in dispersing material. In the idealized model, turbulent jets are generated by prescribed cellular flow with surface convergence and associated downwelling regions. Numeric and analytic model solutions reveal that horizontal jets substantially contribute to along-wind dispersion for sufficiently strong cellular flows and exceed contributions due to vertical mean shear for buoyant surface-trapped material. However, surface convergence regions also accumulate surface-trapped material, reducing shear dispersion by jets. Turbulence resolving LES results of a coastal depth-limited ocean agree qualitatively with the idealized model and reveal long-lived coherent jet structures that are necessary for effective jet dispersion. These coastal results indicate substantial jet contributions to along-wind dispersion. However, jet dispersion is likely less effective in the open ocean because jets are shorter lived, less organized, and distorted due to spiraling Ekman currents.


AIAA Journal ◽  
1985 ◽  
Vol 23 (5) ◽  
pp. 693-700 ◽  
Author(s):  
H. E. Gilreath ◽  
A. Brandt

Sign in / Sign up

Export Citation Format

Share Document