Self-organized criticality of turbulence in strongly stratified mixing layers

Hesam Salehipour,W. R. Peltier,C. P. Caulfield

DOI: https://doi.org/10.1017/jfm.2018.695

2018-09-10

Abstract:Motivated by the importance of stratified shear flows in geophysical and environmental circumstances, we characterize their energetics, mixing and spectral behavior through a series of direct numerical simulations of turbulence generated by Holmboe wave instability (HWI) under various initial conditions. We focus on circumstances where the stratification is sufficiently `strong' so that HWI is the dominant primary instability of the flow. Our numerical findings demonstrate the emergence of self-organised criticality (SOC) that is manifest as an adjustment of an appropriately defined gradient Richardson number, $Ri_g$, associated with the horizontally-averaged mean flow, in such a way that it is continuously attracted towards a critical value of $Ri_g \sim 1/4$. This self-organization occurs through a continuously reinforced localisation of the `scouring' motions (i.e. `avalanches') that are characteristic of the turbulence induced by the break down of Holmboe wave instabilities and are developed on the upper and lower flanks of the sharply localized density interface, embedded within a much more diffuse shear layer. These localised `avalanches' are also found to exhibit the expected scale invariant characteristics. From an energetics perspective, the emergence of SOC is expressed in the form of a long-lived turbulent flow that remains in a `quasi-equilibrium' state for an extended period of time. Most importantly, the irreversible mixing that results from such self-organised behavior appears to be characterized generically by a universal cumulative turbulent flux coefficient of $\Gamma_c \sim 0.2$ only for turbulent flows engendered by Holmboe wave instability. The existence of this self-organised critical state corroborates the original physical arguments associated with self-regulation of stratified turbulent flows as involving a `kind of equilibrium' as described by Turner (1973).

Fluid Dynamics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the self - organized criticality (SOC) of turbulence in strongly stratified shear flows. Specifically, through a series of direct numerical simulations (DNS), the authors studied the energy characteristics, mixing behavior, and spectral behavior of turbulence induced by Holmboe wave instability (HWI) under different initial conditions. The focus of the study is on the behavioral characteristics of these flows when the stratification is "strong" enough so that HWI becomes the main instability of the flow.

### Main problems
1. **Manifestation of self - organized criticality**: The study shows that in this case, the flow exhibits self - organized criticality, manifested as the gradient Richardson number \(R_{ig}\) adjusting to a critical value \(R_{ig}\sim1/4\). This self - organizing process is achieved through localized "erosion" motions (i.e., "avalanches") that are generated due to the breakdown of Holmboe wave instability and mainly occur at the upper and lower edges of the density interface.

2. **Energy balance state**: From an energy perspective, the emergence of self - organized criticality is manifested as a turbulent flow that maintains a "quasi - equilibrium" state for a long time.

3. **Universality of mixing efficiency**: The study indicates that the irreversible mixing resulting from this self - organized behavior has a universal cumulative turbulent flux coefficient \(\Gamma_c\sim0.2\), which is only applicable to turbulent flows induced by Holmboe wave instability.

### Key formulas
- Definition of the gradient Richardson number \(R_{ig}\):
\[
R_{ig}(z,t)=\frac{N^{2}}{S^{2}}=\frac{-\frac{g}{\rho_r}\frac{\partial\rho}{\partial z}}{\left(\frac{\partial u}{\partial z}\right)^{2}}
\]
where \(\rho(z,t)\) and \(u(z,t)\) represent the horizontally averaged density and velocity profiles respectively, \(g\) is the acceleration due to gravity, and \(\rho_r\) is the hydrostatic reference density.

- Definition of the bulk Richardson number \(R_{ib}\):
\[
R_{ib}=\frac{g\rho_0d}{\rho_rU_0^{2}}
\]
where \(\rho_0\) and \(U_0\) represent the characteristic density and velocity changes across the characteristic length \(d\) respectively.

- Definition of the turbulent dissipation rate \(\epsilon'(z,t)\):
\[
\epsilon'(z,t) = 2\nu s'_{ij}s'_{ij}
\]
where \(\nu\) is the kinematic viscosity, \(s'_{ij}=\frac{1}{2}\left(\frac{\partial u'_i}{\partial x_j}+\frac{\partial u'_j}{\partial x_i}\right)\) is the perturbation strain - rate tensor, and \(u'=(u',v',w')\) represents the velocity field deviating from the horizontally averaged mean.

- Definition of the total energy budget \(E_{ST}\):
\[
E_{ST}=K'+PA
\]
where \(K'\) is the perturbation kinetic energy and \(PA\) is the available potential energy.

- Definition of the available potential energy \(PA\):
\[
PA = P - P_B=R_{ib}\langle(\rho-\rho^*)z\rangle
\]

- Energy evolution equation:
\[
\frac{dE_{ST}}{dt}=P - M - D
\]