Turbulent Potential Energy Research Articles

On dissipation timescales of the basic second-order moments: the effect on the energy and flux budget (EFB) turbulence closure for stably stratified turbulence

Abstract. The dissipation rates of the basic second-order moments are the key parameters playing a vital role in turbulence modelling and controlling turbulence energetics and spectra and turbulent fluxes of momentum and heat. In this paper, we use the results of direct numerical simulations (DNSs) to evaluate dissipation rates of the basic second-order moments and revise the energy and flux budget (EFB) turbulence closure theory for stably stratified turbulence. We delve into the theoretical implications of this approach and substantiate our closure hypotheses through DNS data. We also show why the concept of down-gradient turbulent transport becomes incomplete when applied to the vertical turbulent flux of potential temperature under stable stratification. We reveal essential feedback between the turbulent kinetic energy (TKE), the vertical turbulent flux of buoyancy, and the turbulent potential energy (TPE), which is responsible for maintaining shear-produced stably stratified turbulence for any Richardson number.

Understanding the effect of Prandtl number on momentum and scalar mixing rates in neutral and stably stratified flows using gradient field dynamics

Recently, direct numerical simulations (DNS) of stably stratified turbulence have shown that as the Prandtl number ( $Pr$ ) is increased from 1 to 7, the mean turbulent potential energy dissipation rate (TPE-DR) drops dramatically, while the mean turbulent kinetic energy dissipation rate (TKE-DR) increases significantly. Through an analysis of the equations governing the fluctuating velocity and density gradients we provide a mechanistic explanation for this surprising behaviour and test the predictions using DNS. We show that the mean density gradient gives rise to a mechanism that opposes the production of fluctuating density gradients, and this is connected to the emergence of ramp cliffs. The same term appears in the velocity gradient equation but with the opposite sign, and is the contribution from buoyancy. This term is ultimately the reason why the TPE-DR reduces while the TKE-DR increases with increasing $Pr$ . Our analysis also predicts that the effects of buoyancy on the smallest scales of the flow become stronger as $Pr$ is increased, and this is confirmed by our DNS data. A consequence of this is that the standard buoyancy Reynolds number does not correctly estimate the impact of buoyancy at the smallest scales when $Pr$ deviates from 1, and we derive a suitable alternative parameter. Finally, an analysis of the filtered gradient equations reveals that the mean density gradient term changes sign at sufficiently large scales, such that buoyancy acts as a source for velocity gradients at small scales, but as a sink at large scales.

Evaluation of subgrid-scale models in decaying rotating stratified turbulence

The results of large eddy simulation (LES) using four sub-grid scale (SGS) models, namely: constant coefficient Smagorinsky (CS), dynamic Smagorinsky (DS), a dynamic Clark (DC) model, and dynamic Lund-Novikov II-term model (DLN2M) for rotating stratified turbulence in the absence of forcing using large-scale isotropic initial conditions, are reported here. Three cases with varying ratios of Brunt-Väisälä frequency to the inertial wave frequency, N / f , have been chosen to evaluate the performance of SGS models. The Reynolds number and N / f are chosen as (a) Case 1: Re=3704, N / f = 5 , (b) Case 2: Re=6667, N / f = 40 and (c) Case 3: Re=6667, N / f = 138 . Various quantities including turbulent kinetic energy (tke), turbulent potential energy (tpe), total dissipation, potential and total energy spectra, and their fluxes, are analyzed. It is observed that the increase in the value of N / f leads to increased oscillations in the evolution of kinetic and potential energy. Our results suggest that the inclusion of the rotation rate tensor in SGS modeling leads to overall improvements in energy evolution predictions. The CS model produces the largest errors, while the DS and the DC models are more comparable to each other. Spectral analysis of total energy and their fluxes shows that the CS, DS and DC models can reasonably predict the large-scale motion ( κ < 10 ), while the DLN2M model improves the energy prediction in the wavenumber range ( 10 < κ < 64 ) in comparison to other SGS models.

Analysis of scale-dependent kinetic and potential energy in sheared, stably stratified turbulence

Budgets of turbulent kinetic energy (TKE) and turbulent potential energy (TPE) at different scales $\ell$ in sheared, stably stratified turbulence are analysed using a filtering approach. Competing effects in the flow are considered, along with the physical mechanisms governing the energy fluxes between scales, and the budgets are used to analyse data from direct numerical simulation at buoyancy Reynolds number $Re_b=O(100)$ . The mean TKE exceeds the TPE by an order of magnitude at the large scales, with the difference reducing as $\ell$ is decreased. At larger scales, buoyancy is never observed to be positive, with buoyancy always converting TKE to TPE. As $\ell$ is decreased, the probability of locally convecting regions increases, though it remains small at scales down to the Ozmidov scale. The TKE and TPE fluxes between scales are both downscale on average, and their instantaneous values are correlated positively, but not strongly so, and this occurs due to the different physical mechanisms that govern these fluxes. Moreover, the contributions to these fluxes arising from the sub-grid fields are shown to be significant, in addition to the filtered scale contributions associated with the processes of strain self-amplification, vortex stretching and density gradient amplification. Probability density functions (PDFs) of the $Q,R$ invariants of the filtered velocity gradient are considered and show that as $\ell$ increases, the sheared-drop shape of the PDF becomes less pronounced and the PDF becomes more symmetric about $R=0$ .

XXI-st century: a shift of paradigm in the planetary boundary layer theory

At the beginning of 21st century, Sergej Sergeevich Zilitinkevich initiated a shift of the fundamental paradigm in the planetary boundary layer theory, and wider, in statistical hydromechanics of stratified flows. He held the position of the chair in meteorology at Uppsala University, Sweden in 1998–2003, where he and his colleagues consistently introduced a concept of turbulent potential and total energy and not least reviewed fundamental understanding of turbulent mixing in stratified flows. Zilitinkevich continued as a professor-emeritus at University of Helsinki, Finland, in 2004–2021, where he attracted prestigious mega-grants both in Russia and European Union. As such, he created and led a virtual international research institute — the third institute in course of his carrier in science — for development of theoretical foundations and practical implications of the new paradigm for turbulent energies. A network of researchers and practitioners established by Zilitinkevich continue their joint research in the Pan-Eurasian Experiment frameworks.

Assessment of SGS Models for Large Eddy Simulation (LES) of a Stratified Taylor–Green Vortex

Eddy viscosity-based and regularization-based subgrid-scale (SGS) models are quantitatively assessed for large eddy simulations (LES) of stratified Taylor–Green vortex (TGV). LES calculations using the dynamic Smagorinsky, Leray- $$\alpha$$ , LANS- $$\alpha$$ and Clark- $$\alpha$$ models are conducted using a pseudo-spectral formulation on a $$128^{3}$$ grid and using $$\alpha =\frac{1}{160}$$ (for all the regularization models) for Froude number, $$Fr_0=\frac{{\mathcal {U}}}{{{\mathcal {N}}}{{\mathcal {L}}}}=1$$ , and Reynolds number, $$Re_{0}= \frac{\mathcal {U L} }{\nu }=1600$$ . Results are compared in detail with in-house direct numerical simulation (DNS) calculations. Validation of the formulation is carried out through grid-dependent studies using grids of $$64^3$$ , $$128^3$$ and $$256^3$$ , $$\alpha$$ (regularization parameter)-dependent studies using $$\alpha$$ values of $$\frac{1}{40}, \frac{1}{80}, \frac{1}{160}$$ and $$\frac{1}{320}$$ , and comparisons with previously published results Remmler and Hickel, 2012, of the same problem. Various quantities including turbulent kinetic energy (tke), turbulent potential energy (tpe), Q-criteria based vortical structures, potential & total energy spectra, and horizontal & vertical energy fluxes in the respective directions, are analyzed to understand the capability of the various LES models in predicting the development of stratified turbulence. Results showed that all the SGS model predictions were very similar to each other with Smagorinsky and Leray displaying the lowest and highest errors, respectively, in comparison to DNS. The effect of $$Fr_0$$ was also studied through additional LES calculations at $$Fr_0=0.5$$ and 0.16, again using the dynamic Smagorinsky and regularization models and comparing various parameters including tke, tpe and energy spectra to DNS. Decay rates of tke and tpe decreased and turbulence was significantly suppressed with decrease in $$Fr_0$$ . At these lower $$Fr_0$$ , all the SGS model predictions were almost identical to each other and compared extremely well with the DNS results.

Flow kinematics and air entrainment under laboratory spilling breaking waves

Laboratory measurements of velocity fields and void fraction under spilling breaking waves are presented. Modified particle image velocimetry was used to quantify the flow kinematics and turbulence while fibre optic reflectometry was used to quantify the breaking-induced air entrainment inside the aerated region of the spilling breakers. The measurements confirmed that the ratio of the local energy flux and the local energy density sharply increases and exceeds the threshold value of 0.85 near the onset of breaking. Based on the measured velocity fields, the maximum horizontal velocity reached at the onset of breaking, with being the phase speed of the primary breaking wave. The maximum horizontal velocity then reached at approximately one-quarter of a wave period after the onset of breaking. The results also confirmed that the wavelet-educed turbulence length scale estimates are comparable to the previously reported values with different wave parameters, suggesting that the dependence of the size of energy-containing eddies on the physical scales of the breaking waves is insignificant. The measured void fraction showed a similarity profile although the measurement locations span one wavelength. The mean kinetic energy, turbulent kinetic energy, potential energy and total energy were quantified with and without the void fraction being accounted for. Results show near 60 % and 40 % overestimates of the kinetic energy and the potential energy if the void fraction is not considered. After correcting the density variation due to air entrainment, the total energy dissipated following an exponential decay, with 43 % and 65 % energy being dissipated at one and two wavelengths downstream from the breaking point, respectively. The equipartition assumption was found to be applicable before and during the entire breaking process in the present spilling breakers.

Internal gravity waves in the energy and flux budget turbulence-closure theory for shear-free stably stratified flows.

We have advanced the energy and flux budget turbulence closure theory that takes into account a two-way coupling between internal gravity waves (IGWs) and the shear-free stably stratified turbulence. This theory is based on the budget equation for the total (kinetic plus potential) energy of IGWs, the budget equations for the kinetic and potential energies of fluid turbulence, and turbulent fluxes of potential temperature for waves and fluid flow. The waves emitted at a certain level propagate upward, and the losses of wave energy cause the production of turbulence energy. We demonstrate that due to the nonlinear effects more intensive waves produce more strong turbulence, and this, in turn, results in strong damping of IGWs. As a result, the penetration length of more intensive waves is shorter than that of less intensive IGWs. The anisotropy of the turbulence produced by less intensive IGWs is stronger than that caused by more intensive waves. The low-amplitude IGWs produce turbulence consisting up to 90% of turbulent potential energy. This resembles the properties of the observed high-altitude tropospheric strongly anisotropic (nearly two-dimensional) turbulence.

Energetics and mixing in buoyancy-driven near-bottom stratified flow

Turbulence and mixing in a near-bottom convectively driven flow are examined by numerical simulations of a model problem: a statically unstable disturbance at a slope with inclination $\unicode[STIX]{x1D6FD}$ in a stable background with buoyancy frequency $N$ . The influence of slope angle and initial disturbance amplitude are quantified in a parametric study. The flow evolution involves energy exchange between four energy reservoirs, namely the mean and turbulent components of kinetic energy (KE) and available potential energy (APE). In contrast to the zero-slope case where the mean flow is negligible, the presence of a slope leads to a current that oscillates with $\unicode[STIX]{x1D714}=N\sin \unicode[STIX]{x1D6FD}$ and qualitatively changes the subsequent evolution of the initial density disturbance. The frequency, $N\sin \unicode[STIX]{x1D6FD}$ , and the initial speed of the current are predicted using linear theory. The energy transfer in the sloping cases is dominated by an oscillatory exchange between mean APE and mean KE with a transfer to turbulence at specific phases. In all simulated cases, the positive buoyancy flux during episodes of convective instability at the zero-velocity phase is the dominant contributor to turbulent kinetic energy (TKE) although the shear production becomes increasingly important with increasing $\unicode[STIX]{x1D6FD}$ . Energy that initially resides wholly in mean available potential energy is lost through conversion to turbulence and the subsequent dissipation of TKE and turbulent available potential energy. A key result is that, in contrast to the explosive loss of energy during the initial convective instability in the non-sloping case, the sloping cases exhibit a more gradual energy loss that is sustained over a long time interval. The slope-parallel oscillation introduces a new flow time scale $T=2\unicode[STIX]{x03C0}/(N\sin \unicode[STIX]{x1D6FD})$ and, consequently, the fraction of initial APE that is converted to turbulence during convective instability progressively decreases with increasing $\unicode[STIX]{x1D6FD}$ . For moderate slopes with $\unicode[STIX]{x1D6FD}&lt;10^{\circ }$ , most of the net energy loss takes place during an initial, short ( $Nt\approx 20$ ) interval with periodic convective overturns. For steeper slopes, most of the energy loss takes place during a later, long ( $Nt&gt;100$ ) interval when both shear and convective instability occur, and the energy loss rate is approximately constant. The mixing efficiency during the initial period dominated by convectively driven turbulence is found to be substantially higher (exceeds 0.5) than the widely used value of 0.2. The mixing efficiency at long time in the present problem of a convective overturn at a boundary varies between 0.24 and 0.3.

On the inference of the state of turbulence and mixing efficiency in stably stratified flows

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}&gt;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}&lt;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.

Turbulent Mechanical Energy Budget in Stably Stratified Baroclinic Flows over Sloping Terrain

Analysis of second-moment budget equations in a slope-oriented coordinate frame exhibits the pathways of exchange between the potential energy of mean flow and the total turbulent mechanical energy. It is shown that this process is controlled by the inclination of the potential temperature gradient. Hence, this parameter should be considered in studies of turbulence in slope flows as well as the slope inclination. The concept of turbulent potential energy is generalized to include baroclinicity, and is used to explain the role of along-slope turbulent heat flux in energy conversions. A generalization of static stability criteria for baroclinic conditions is also proposed. In addition, the presence of feedback between the turbulent heat flux and the temperature variance in stably-stratified flows is identified, which implies the existence of oscillatory modes characterized by the Brunt–Väisäla frequency.

Statistical Modeling of Stratified Two-Phase Flow

A new numerical model for stratified two-phase flows with wavy interface is derived in this study. Assuming an equilibrium condition between turbulent kinetic energy, potential energy, and surface energy, the turbulent length scale in the inner region of a two-layer turbulence approach can be described by a statistical model to account for the influence of the waves. The average wave number, which is an input parameter to this model, is taken from wave spectra. They can be predicted from a Boltzmann statistic of turbulent kinetic energy. The new turbulence model is compared with the two-phase k–ϵ turbulence model. Time-averaged flow properties calculated by the new approach, such as velocity, turbulence, and void profiles, are shown to be in good agreement with experimental data.

Improving Wind-Ramp Forecasts in the Stable Boundary Layer

The viability of wind-energy generation is dependent on highly accurate numerical wind forecasts, which are impeded by inaccuracies in model representation of boundary-layer processes. This study revisits the basic theory of the Mellor, Yamada, Nakanishi, and Niino (MYNN) planetary boundary-layer parametrization scheme, focusing on the onset of wind-ramp events related to nocturnal low-level jets. Modifications to the MYNN scheme include: (1) calculation of new closure parameters that determine the relative effects of turbulent energy production, dissipation, and redistribution; (2) enhanced mixing in the stable boundary layer when the mean wind speed exceeds a specified threshold; (3) explicit accounting of turbulent potential energy in the energy budget. A mesoscale model is used to generate short-term (24 h) wind forecasts for a set of 15 cases from both the U.S.A. and Germany. Results show that the new set of closure parameters provides a marked forecast improvement only when used in conjunction with the new mixing length formulation and only for cases that are originally under- or over-forecast (10 of the 15 cases). For these cases, the mean absolute error (MAE) of wind forecasts at turbine-hub height is reduced on average by 17%. A reduction in MAE values on average by 26% is realized for these same cases when accounting for the turbulent potential energy together with the new mixing length. This last method results in an average reduction by at least 13% in MAE values across all 15 cases.

Phenomenology of two-dimensional stably stratified turbulence under large-scale forcing

ABSTRACTIn this paper, we characterise the scaling of energy spectra, and the interscale transfer of energy and enstrophy, for strongly, moderately and weakly stably stratified two-dimensional (2D) turbulence, restricted in a vertical plane, under large-scale random forcing. In the strongly stratified case, a large-scale vertically sheared horizontal flow (VSHF) coexists with small scale turbulence. The VSHF consists of internal gravity waves and the turbulent flow has a kinetic energy (KE) spectrum that follows an approximate k−3 scaling with zero KE flux and a robust positive enstrophy flux. The spectrum of the turbulent potential energy (PE) also approximately follows a k−3 power-law and its flux is directed to small scales. For moderate stratification, there is no VSHF and the KE of the turbulent flow exhibits Bolgiano–Obukhov scaling that transitions from a shallow k−11/5 form at large scales, to a steeper approximate k−3 scaling at small scales. The entire range of scales shows a strong forward enstrophy flux, and interestingly, large (small) scales show an inverse (forward) KE flux. The PE flux in this regime is directed to small scales, and the PE spectrum is characterised by an approximate k−1.64 scaling. Finally, for weak stratification, KE is transferred upscale and its spectrum closely follows a k−2.5 scaling, while PE exhibits a forward transfer and its spectrum shows an approximate k−1.6 power-law. For all stratification strengths, the total energy always flows from large to small scales and almost all the spectral indicies are well explained by accounting for the scale-dependent nature of the corresponding flux.

Closure Schemes for Stably Stratified Atmospheric Flows without Turbulence Cutoff

Abstract Two recently proposed turbulence closure schemes are compared against the conventional Mellor–Yamada (MY) model for stably stratified atmospheric flows. The Energy- and Flux-Budget (EFB) approach solves the budgets of turbulent momentum and heat fluxes and turbulent kinetic and potential energies. The Cospectral Budget (CSB) approach is formulated in wavenumber space and integrated across all turbulent scales to obtain flow variables in physical space. Unlike the MY model, which is subject to a “critical gradient Richardson number,” both EFB and CSB models allow turbulence to exist at any gradient Richardson number and predict a saturation of flux Richardson number () at sufficiently large . The CSB approach further predicts the value of and reveals a unique expression linking the Rotta and von Kármán constants. Hence, all constants in the CSB model are nontunable and stability independent. All models agree that the dimensionless sensible heat flux decays with increasing . However, the decay rate and subsequent cutoff in the MY model appear abrupt. The MY model further exhibits an abrupt cutoff in the turbulent stress normalized by vertical velocity variance, while the CSB and EFB models display increasing trends. The EFB model produces a rapid increase in the ratio of turbulent potential energy and vertical velocity variance as is approached, suggesting a strong self-preservation mechanism. Vertical anisotropy in the turbulent kinetic energy is parameterized in different ways in MY and EFB, but this consideration is not required in CSB. Differences between EFB and CSB model predictions originate from how the vertical anisotropy is specified in the EFB model.

On the flux Richardson number in stably stratified turbulence

The flux Richardson number $R_{f}$ (often referred to as the mixing efficiency) is a widely used parameter in stably stratified turbulence which is intended to provide a measure of the amount of turbulent kinetic energy $k$ that is irreversibly converted to background potential energy (which is by definition the minimum potential energy that a stratified fluid can attain that is not available for conversion back to kinetic energy) due to turbulent mixing. The flux Richardson number is traditionally defined as the ratio of the buoyancy flux $B$ to the production rate of turbulent kinetic energy $P$. An alternative generalized definition for $R_{f}$ was proposed by Ivey &amp; Imberger (J. Phys. Oceanogr., vol. 21, 1991, pp. 650–658), where the non-local transport terms as well as unsteady contributions were included as additional sources to the production rate of $k$. While this definition precludes the need to assume that turbulence is statistically stationary, it does not properly account for countergradient fluxes that are typically present in more strongly stratified flows. Hence, a third definition that more rigorously accounts for only the irreversible conversions of energy has been defined, where only the irreversible fluxes of buoyancy and production (i.e. the dissipation rates of $k$ and turbulent potential energy ($E_{PE}^{\prime }$)) are used. For stationary homogeneous shear flows, all of the three definitions produce identical results. However, because stationary and/or homogeneous flows are typically not found in realistic geophysical situations, clarification of the differences/similarities between these various definitions of $R_{f}$ is imperative. This is especially true given the critical role $R_{f}$ plays in inferring turbulent momentum and heat fluxes using indirect methods, as is commonly done in oceanography, and for turbulence closure parameterizations. To this end, a careful analysis of two existing direct numerical simulation (DNS) datasets of stably stratified homogeneous shear and channel flows was undertaken in the present study to compare and contrast these various definitions. We find that all three definitions are approximately equivalent when the gradient Richardson number $Ri_{g}\leqslant 1/4$. Here, $Ri_{g}=N^{2}/S^{2}$, where $N$ is the buoyancy frequency and $S$ is the mean shear rate, provides a measure of the stability of the flow. However, when $Ri_{g}&gt;1/4$, significant differences are noticeable between the various definitions. In addition, the irreversible formulation of $R_{f}$ based on the dissipation rates of $k$ and $E_{PE}^{\prime }$ is the only definition that is free from oscillations at higher gradient Richardson numbers. Both the traditional definition and the generalized definition of $R_{f}$ exhibit significant oscillations due to the persistence of linear internal wave motions and countergradient fluxes that result in reversible exchanges between $k$ and $E_{PE}^{\prime }$. Finally, we present a simple parameterization for the irreversible flux Richardson number $R_{f}^{\ast }$ based on $Ri_{g}$ that produces excellent agreement with the DNS results for $R_{f}^{\ast }$.

On the application of a turbulence closure modified model to the description of the density jump evolution in a stably stratified medium

The self-similar turbulent density jump evolution has been studied in the scope of a turbulence closure modernized theory which takes into account the anisotropy and mutual transformation of the turbulent fluctuation kinetic and potential energy for a stably stratified fluid. The numerical calculation, performed using the equations for the average density and kinetic and potential energies of turbulent fluctuations, indicates that the vertical profiles of the buoyancy frequency, turbulence scale, and kinetic and potential energies drastically change when the turbulence anisotropy is strong. The vertical profiles of the corresponding energy and spatial discontinuity parameters, calculated at a weaker anisotropy, indicate that similar drastic changes are absent and a qualitative agreement exists with the known analytical solution, which describes the density jump evolution in a freshwater basin and was obtained previously [5, 8] in the scope of a turbulence local-similarity hypothesis applied in combination with the budget equation for the turbulent fluctuation kinetic energy.

Effects of three-dimensionality on instability and turbulence in a frontal zone

Linear stability analysis and direct numerical simulation are used to investigate the evolution of a symmetrically unstable uniform frontal zone. Simulations in a three-dimensional computational domain capable of resolving near-symmetric currents develop strong nonlinearities without the emergence of pure symmetric instability. Linear stability analysis demonstrates that for $ft&gt;1$ ( $f$ is the Coriolis parameter and $t$ denotes time) the flow generates strongly asymmetric structures which become nearly symmetric when $ft\gg 1$. Unlike the currents generated during pure symmetric instability, near-symmetric instability generates currents that do not align with isopycnals. This greatly modifies their energetics and evolution, leading to regions of the flow that are unstable to gravitational instability and energized by the reservoir of available potential energy. A high-resolution simulation demonstrates the flow evolution from near-symmetric currents to secondary shear-convective instabilities and finally, through tertiary instabilities, to fully three-dimensional turbulence. The effect of this sequence of instabilities is quantified through velocity and vorticity statistics as well as budgets for turbulent kinetic and potential energy. It is not until $ft\sim 10$ that the energy source for fluctuations is primarily shear, in contrast to the purely symmetric instability which draws its energy exclusively from shear production.

Biases in Thorpe-Scale Estimates of Turbulence Dissipation. Part II: Energetics Arguments and Turbulence Simulations

AbstractThis paper uses the energetics framework developed by Scotti and White to provide a critical assessment of the widely used Thorpe-scale method, which is used to estimate dissipation and mixing rates in stratified turbulent flows from density measurements along vertical profiles. This study shows that the relevant displacement scale in general is not the rms value of the Thorpe displacement. Rather, the displacement field must be Reynolds decomposed to separate the mean from the turbulent component, and it is the turbulent component that ought to be used to diagnose mixing and dissipation. In general, the energetics of mixing in an overall stably stratified flow involves potentially complex exchanges among the available potential energy and kinetic energy associated with the mean and turbulent components of the flow. The author considers two limiting cases: shear-driven mixing, where mixing comes at the expense of the mean kinetic energy of the flow, and convective-driven mixing, which taps the available potential energy of the mean flow to drive mixing. In shear-driven flows, the rms of the Thorpe displacement, known as the Thorpe scale is shown to be equivalent to the turbulent component of the displacement. In this case, the Thorpe scale approximates the Ozmidov scale, or, which is the same, the Thorpe scale is the appropriate scale to diagnose mixing and dissipation. However, when mixing is driven by the available potential energy of the mean flow (convective-driven mixing), this study shows that the Thorpe scale is (much) larger than the Ozmidov scale. Using the rms of the Thorpe displacement overestimates dissipation and mixing, since the amount of turbulent available potential energy (measured by the turbulent displacement) is only a fraction of the total available potential energy (measured by the Thorpe scale). Corrective measures are discussed that can be used to diagnose mixing from knowledge of the Thorpe displacement. In a companion paper, Mater et al. analyze field data and show that the Thorpe scale can indeed be much larger than the Ozmidov scale.

A numerical study of a weakly stratified turbulent wake

Direct numerical simulation (DNS) of a time-developing turbulent wake evolving in a stably stratified background is presented. A large initial Froude number is chosen to allow the wake to become fully turbulent and axisymmetric before stratification affects the spreading rate of the mean defect. Turbulence statistics are formed by averaging over the homogeneous streamwise direction of a domain that is larger than earlier stratified-wake simulations in order to reduce the statistical uncertainty. The DNS results are used to cast light on the mechanisms that lead to the various states of this flow – namely the three-dimensional (essentially unstratified), non-equilibrium (or ‘wake-collapse’) and quasi-two-dimensional (or ‘two-component’) regimes, previously observed for wakes embedded in both weakly and strongly stratified backgrounds. For this relatively high-initial-Reynolds- and Froude-number simulation, we find that the signature reduction in the rate of decay of the maximum mean defect velocity during the wake-collapse regime is due to buoyancy-induced alterations of the turbulence structure, which weaken and redistribute the Reynolds shear stresses whose gradients appear in the streamwise mean momentum equation. The change in the rate of decay of the turbulence kinetic energy (TKE) observed during the wake-collapse regime (which occurs well after the mean velocity decay reduction begins) is not caused by transfer of turbulent gravitational potential energy to TKE, as has been previously suggested. The results instead reveal that the reduction in TKE decay – which for this flow, with its relatively weak internal waves, eventually leads to TKE growth, heralding the arrival of the two-component regime – is caused by an increase in the rate of TKE production associated with the wake structure becoming increasingly two-dimensional, such that the lateral Reynolds shear stress,$-\overline{u^{\prime }v^{\prime }}$, becomes dominant. The present results are also compared with those of previous simulations at different Froude and Reynolds numbers, and whose initial conditions contain different turbulence structures. This comparison confirms a strong degree of commonality in the late-wake behaviour, which lends support to the hypothesis that all wakes in stably stratified environments achieve a universal state in the final stages of decay.