scholarly journals How we compute N matters to estimates of mixing in stratified flows

2017 ◽  
Vol 831 ◽  
Author(s):  
Robert S. Arthur ◽  
Subhas K. Venayagamoorthy ◽  
Jeffrey R. Koseff ◽  
Oliver B. Fringer

Most commonly used models for turbulent mixing in the ocean rely on a background stratification against which turbulence must work to stir the fluid. While this background stratification is typically well defined in idealized numerical models, it is more difficult to capture in observations. Here, a potential discrepancy in ocean mixing estimates due to the chosen calculation of the background stratification is explored using direct numerical simulation data of breaking internal waves on slopes. Two different methods for computing the buoyancy frequency $N$, one based on a three-dimensionally sorted density field (often used in numerical models) and the other based on locally sorted vertical density profiles (often used in the field), are used to quantify the effect of $N$ on turbulence quantities. It is shown that how $N$ is calculated changes not only the flux Richardson number $R_{f}$, which is often used to parameterize turbulent mixing, but also the turbulence activity number or the Gibson number $Gi$, leading to potential errors in estimates of the mixing efficiency using $Gi$-based parameterizations.

2008 ◽  
Vol 598 ◽  
pp. 67-80 ◽  
Author(s):  
HANS VAN HAREN

The large-scale vertical density stratification, represented by buoyancy frequency N, is generally very stable in the upper half of the ocean, and relatively weak in the lower half. However, closer inspection of density profiles demonstrates steps rather than a smooth increase with depth. As is demonstrated here using Richardson number, geostrophic balance and slantwise convective mixing arguments, these layers have a limited set of minimum, weak stratification, N-values Nmin indicating the transition between stably stratified and convective ‘homogeneous’ layers. Adopting the viewpoint that the transition occurs for neutral stability in the direction of Earth's rotation Ω instead of gravity g, three discrete states are hypothesized for mid-latitudes: (i) Nmin = 2fh under linear stability conditions, (ii) Nmin = fh(|ϕ| < 45°) and (iii) Nmin = 4fh, both under nonlinear stability, where horizontal component fh = 2Ω cos ϕ at latitude ϕ. The Nmin are not in terms of inertial frequency f = 2Ω sin ϕ, because the effect of fh is the tilting of vortex tubes away from the local vertical in the direction of Ω. The above explains very well deep-ocean North-Atlantic and Mediterranean observations on transitions in conductivity-temperature with depth profiles, inertial polarization and near-inertial shear. The latter peaks at sub-inertial 0.97f, which is associated with the lower inertio-gravity wave limit for Nmin = 4fh, thereby stressing the importance of fh for the dominant physics associated with mixing in the ocean.


2017 ◽  
Vol 837 ◽  
pp. 129-146 ◽  
Author(s):  
Jason Olsthoorn ◽  
Stuart B. Dalziel

The study of vortex-ring-induced mixing has been significant for understanding stratified turbulent mixing in the absence of a mean flow. Renewed interest in this topic has prompted the development of a one-dimensional model for the evolution of a stratified system in the context of isolated mixing events. This model is compared to numerical simulations and physical experiments of vortex rings interacting with a stratification. Qualitative agreement between the evolution of the density profiles is observed, along with close quantitative agreement of the mixing efficiency. This model highlights the key dynamical features of such isolated mixing events.


2014 ◽  
Vol 756 ◽  
pp. 1027-1057 ◽  
Author(s):  
Megan S. Davies Wykes ◽  
Stuart B. Dalziel

AbstractBoussinesq salt-water laboratory experiments of Rayleigh–Taylor instability (RTI) can achieve mixing efficiencies greater than 0.75 when the unstable interface is confined between two stable stratifications. This is much greater than that found when RTI occurs between two homogeneous layers when the mixing efficiency has been found to approach 0.5. Here, the mixing efficiency is defined as the ratio of energy used in mixing compared with the energy available for mixing. If the initial and final states are quiescent then the mixing efficiency can be calculated from experiments by comparison of the corresponding density profiles. Varying the functional form of the confining stratifications has a strong effect on the mixing efficiency. We derive a buoyancy-diffusion model for the rate of growth of the turbulent mixing region, $\def \xmlpi #1{}\def \mathsfbi #1{\boldsymbol {\mathsf {#1}}}\let \le =\leqslant \let \leq =\leqslant \let \ge =\geqslant \let \geq =\geqslant \def \Pr {\mathit {Pr}}\def \Fr {\mathit {Fr}}\def \Rey {\mathit {Re}}\dot{h} = 2 \sqrt{\alpha A g h}$ (where $A = A(h)$ is the Atwood number across the mixing region when it extends a height $h$, $g$ is acceleration due to gravity and $\alpha $ is a constant). This model shows good agreement with experiments when the value of the constant $\alpha $ is set to 0.07, the value found in experiments of RTI between two homogeneous layers (where the height of the turbulent mixing region increases as $h =\alpha A g t^2$, an expression which is equivalent to that derived for $\dot{h}$).


2019 ◽  
Vol 867 ◽  
pp. 323-333 ◽  
Author(s):  
Amrapalli Garanaik ◽  
Subhas K. Venayagamoorthy

Scaling arguments are presented to quantify the widely used diapycnal (irreversible) mixing coefficient $\unicode[STIX]{x1D6E4}=\unicode[STIX]{x1D716}_{PE}/\unicode[STIX]{x1D716}$ in stratified flows as a function of the turbulent Froude number $Fr=\unicode[STIX]{x1D716}/Nk$. Here, $N$ is the buoyancy frequency, $k$ is the turbulent kinetic energy, $\unicode[STIX]{x1D716}$ is the rate of dissipation of turbulent kinetic energy and $\unicode[STIX]{x1D716}_{PE}$ is the rate of dissipation of turbulent potential energy. We show that for $Fr\gg 1$, $\unicode[STIX]{x1D6E4}\propto Fr^{-2}$, for $Fr\sim \mathit{O}(1)$, $\unicode[STIX]{x1D6E4}\propto Fr^{-1}$ and for $Fr\ll 1$, $\unicode[STIX]{x1D6E4}\propto Fr^{0}$. These scaling results are tested using high-resolution direct numerical simulation (DNS) data from three different studies and are found to hold reasonably well across a wide range of $Fr$ that encompasses weakly stratified to strongly stratified flow conditions. Given that the $Fr$ cannot be readily computed from direct field measurements, we propose a practical approach that can be used to infer the $Fr$ from readily measurable quantities in the field. Scaling analyses show that $Fr\propto (L_{T}/L_{O})^{-2}$ for $L_{T}/L_{O}>O(1)$, $Fr\propto (L_{T}/L_{O})^{-1}$ for $L_{T}/L_{O}\sim O(1)$, and $Fr\propto (L_{T}/L_{O})^{-2/3}$ for $L_{T}/L_{O}<O(1)$, where $L_{T}$ is the Thorpe length scale and $L_{O}$ is the Ozmidov length scale. These formulations are also tested with DNS data to highlight their validity. These novel findings could prove to be a significant breakthrough not only in providing a unifying (and practically useful) parameterization for the mixing efficiency in stably stratified turbulence but also for inferring the dynamic state of turbulence in geophysical flows.


2015 ◽  
Vol 781 ◽  
pp. 113-126 ◽  
Author(s):  
Jason Olsthoorn ◽  
Stuart B. Dalziel

There is tantalizing evidence that some mechanically driven stratified flows tend towards a state of constant mixing efficiency. We provide insight into the energy balance leading to the constant mixing efficiency and isolate the responsible mechanism. The work presented demonstrates an important mixing efficiency regime for periodically forced externally driven stratified flows. Externally forced stratified turbulent mixing is often characterized by the associated eddies within the flow, which are the dominant mixing mechanism (Turner, J. Fluid Mech., vol. 173, 1986, pp. 431–471). Here, we study mixing induced by vortex rings in order to characterize the mixing induced by an individual eddy. By generating a long sequence of independent vortex-ring mixing events in a density-stratified fluid with a sharp interface, we determine the mixing efficiency of each ring. After an initial adjustment phase, we find that the mixing efficiency of each vortex ring is independent of the Richardson number. By studying the mixing mechanism here, we demonstrate consistent features of a volumetrically confined, periodically forced external mixing regime.


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.


2006 ◽  
Vol 128 (4) ◽  
pp. 874-879 ◽  
Author(s):  
Roberto C. Aguirre ◽  
Jennifer C. Nathman ◽  
Haris C. Catrakis

Flow geometry effects are examined on the turbulent mixing efficiency quantified as the mixture fraction. Two different flow geometries are compared at similar Reynolds numbers, Schmidt numbers, and growth rates, with fully developed turbulence conditions. The two geometries are the round jet and the single-stream planar shear layer. At the flow conditions examined, the jet exhibits an ensemble-averaged mixing efficiency which is approximately double the value for the shear layer. This substantial difference is explained fluid mechanically in terms of the distinct large-scale entrainment and mixing-initiation environments and is therefore directly due to flow geometry effects.


2009 ◽  
Vol 66 (4) ◽  
pp. 1063-1068 ◽  
Author(s):  
Ross Tulloch ◽  
K. Shafer Smith

Abstract The quasigeostrophic equations consist of the advection of linearized potential vorticity coupled with advection of temperature at the bounding upper and lower surfaces. Numerical models of quasigeostrophic flow often employ greater (scaled) resolution in the horizontal than in the vertical (the two-layer model is an extreme example). In the interior, this has the effect of suppressing interactions between layers at horizontal scales that are small compared to Nδz/f (where δz is the vertical resolution, N the buoyancy frequency, and f the Coriolis parameter). The nature of the turbulent cascade in the interior is, however, not fundamentally altered because the downscale cascade of potential enstrophy in quasigeostrophic turbulence and the downscale cascade of enstrophy in two-dimensional turbulence (occurring layerwise) both yield energy spectra with slopes of −3. It is shown here that a similar restriction on the vertical resolution applies to the representation of horizontal motions at the surfaces, but the penalty for underresolving in the vertical is complete suppression of the surface temperature cascade at small scales and a corresponding artificial steepening of the surface energy spectrum. This effect is demonstrated in the nonlinear Eady model, using a finite-difference representation in comparison with a model that explicitly advects temperature at the upper and lower surfaces. Theoretical predictions for the spectrum of turbulence in the nonlinear Eady model are reviewed and compared to the simulated flows, showing that the latter model yields an accurate representation of the cascade dynamics. To accurately represent dynamics at horizontal wavenumber K in the vertically finite-differenced model, it is found that the vertical grid spacing must satisfy δz ≲ 0.3f/(NK); at wavenumbers K &gt; 0.3f/(Nδz), the spectrum of temperature variance rolls off rapidly.


Ocean Science ◽  
2011 ◽  
Vol 7 (2) ◽  
pp. 219-229 ◽  
Author(s):  
L. de Montera ◽  
M. Jouini ◽  
S. Verrier ◽  
S. Thiria ◽  
M. Crepon

Abstract. Phytoplankton patchiness has been investigated with multifractal analysis techniques. We analyzed oceanic chlorophyll maps, measured by the SeaWiFS orbiting sensor, which are considered to be good proxies for phytoplankton. The study area is the Senegalo-Mauritanian upwelling region, because it has a low cloud cover and high chlorophyll concentrations. Multifractal properties are observed, from the sub-mesoscale up to the mesoscale, and are found to be consistent with the Corssin-Obukhov scale law of passive scalars. This result indicates that, in this specific region and within this scale range, turbulent mixing would be the dominant effect leading to the observed variability of phytoplankton fields. Finally, it is shown that multifractal patchiness can be responsible for significant biases in the nonlinear source and sink terms involved in biogeochemical numerical models.


Icarus ◽  
1980 ◽  
Vol 41 (3) ◽  
pp. 410-422 ◽  
Author(s):  
André Marten ◽  
Régis Courtin ◽  
Daniel Gautier ◽  
Anne Lacombe

Sign in / Sign up

Export Citation Format

Share Document