Linear and nonlinear contributions to the localization and dynamics of internal gravity waves in a stably stratified turbulent channel flow are investigated using data from direct numerical simulations (DNS). The classification into linear and nonlinear mechanisms is based on the resolvent formulation of the Navier–Stokes equations, which interprets velocity and temperature fluctuations (flow response) as the result of a linear operator (resolvent) acting on the nonlinear advection terms (forcing). Spatial and spatio-temporal power spectral densities computed from DNS data demonstrate that the stratified flow response is localized in spectral space and in the channel core, while the nonlinear forcing is broadband and spans up to the entire channel height. The localization of the velocity and temperature fluctuations in wavenumber and frequency is captured by the leading singular value of the resolvent operator. The wall-normal localization on the other hand results from a combination of linear dynamics and nonlinear forcing, and the latter is further examined using the cross-spectral density (CSD) tensor. Wall-normal subsets of the forcing CSD lead to flow responses that reveal a three-layer structure. The middle one hosts the critical layer of the gravity wave, and is termed the outer layer since it is flanked by an inner layer at the wall and the core region at the channel centre. Forcing within this outer layer generates the majority of the flow response in the channel core. Furthermore, a decomposition of the forcing CSD into velocity and temperature demonstrates that each imprints distinct phase relations on their associated responses, which lead to destructive interference and localization of the gravity waves in the channel core.

The seminal Bolgiano–Obukhov (BO) theory established the fundamental framework for turbulent mixing and energy transfer in stably stratified fluids. However, the presence of BO scalings remains debatable despite their being observed in stably stratified atmospheric layers and convective turbulence. In this study, we performed precise temperature measurements with 51 high-resolution loggers above the seafloor for 46 h on the continental shelf of the northern South China Sea. The temperature observation exhibits three layers with increasing distance from the seafloor: the bottom mixed layer (BML), the mixing zone and the internal wave zone. A BO-like scaling $\alpha =-1.34\pm 0.10$ is observed in the temperature spectrum when the BML is in a weakly stable stratified ($N\sim 0.0018$ rad s$^{-1}$) and strongly sheared ($Ri\sim 0.0027$) condition, whereas in the unstably stratified convective turbulence of the BML, the scaling $\alpha =-1.76\pm 0.10$ clearly deviated from the BO theory but approached the classical $-$5/3 scaling in isotropic turbulence. This suggests that the convective turbulence is not the promise of BO scaling. In the mixing zone, where internal waves alternately interact with the BML, the scaling follows the Kolmogorov scaling. In the internal wave zone, the scaling $\alpha =-2.12 \pm 0.15$ is observed in the turbulence range and possible mechanisms are provided.

The flow field of a bluff body, a circular disk, that moves horizontally in a stratified environment is studied using large-eddy simulations. Five levels of stratification (body Froude numbers of ${{Fr}} = 0.5, 1, 1.5, 2$ and $5$) are simulated at Reynolds number of ${{Re}} = 5000$ and Prandtl number of $Pr =1$. A higher ${{Re}} = 50\,000$ database at ${{Fr}} = 2, 10$ and $Pr =1$ is also examined for comparison. The wavelengths and amplitudes of steady lee waves are compared with a linear-theory analysis. Excellent agreement is found over the entire range of ${{Fr}}$ if an ‘equivalent body’ that includes the separation region is employed for the linear theory. For asymptotically large distances, the velocity amplitude varies theoretically as ${{Fr}}^{-1}$ but a correction owing to the dependence of the separation zone on ${{Fr}}$ is needed. The wake waves propagate in a narrow band of angles with the vertical, and have a wavelength that increases with increasing ${{Fr}}$. The envelope of wake waves, demarcated using buoyancy variance, exhibits self-similar behaviour. The higher ${{Re}}$ results are consistent with the buoyancy effects exhibited at the lower ${{Re}}$. The wake wave energy is larger at ${{Re}} = 50\,000$. Nevertheless, independent of ${{Fr}}$ and ${{Re}}$, the ratio of the wake wave potential energy to the wake turbulent energy increases to approximately 0.6–0.7 in the non-equilibrium stage showing their energetic importance besides suggesting universality in this statistic. There is a crossover of energetic dominance of lee waves at ${{Fr}} <2$ to wake-wave dominance at ${{Fr}} \approx 5$.

Quantifying transport by strongly stratified turbulence in low Prandtl number ($Pr$) fluids is critically important for the development of better models for the structure and evolution of stellar and planetary interiors. Motivated by recent numerical simulations showing strongly anisotropic flows suggestive of a scale-separated dynamics, we perform a multiscale asymptotic analysis of the governing equations. We find that, in all cases, the resulting slow–fast systems take a quasilinear form. Our analysis also reveals the existence of several distinct dynamical regimes depending on the emergent buoyancy Reynolds and Péclet numbers, $Re_b = \alpha ^2 Re$ and $Pe_b = Pr Re_b$, respectively, where $\alpha$ is the aspect ratio of the large-scale turbulent flow structures, and $Re$ is the outer-scale Reynolds number. Scaling relationships relating the aspect ratio, the characteristic vertical velocity and the strength of the stratification (measured by the Froude number $Fr$) naturally emerge from the analysis. When $Pe_b \ll \alpha$, the dynamics at all scales is dominated by buoyancy diffusion, and our results recover the scaling laws empirically obtained from direct numerical simulations by Cope et al. (J. Fluid Mech., vol. 903, 2020, A1). For $Pe_b \ge O(1)$, diffusion is negligible (or at least subdominant) at all scales and our results are consistent with those of Chini et al. (J. Fluid Mech., vol. 933, 2022) for strongly stratified geophysical turbulence at $Pr =O(1)$. Finally, we have identified a new regime for $\alpha \ll Pe_b \ll 1$, in which slow, large scales are diffusive while fast, small scales are not. We conclude by presenting a map of parameter space that clearly indicates the transitions between isotropic turbulence, non-diffusive stratified turbulence, diffusive stratified turbulence and viscously dominated flows, and by proposing parameterisations of the buoyancy flux, mixing efficiency and turbulent diffusion coefficient for each regime.

We present a new derivation of the kinetic equation for weak, non-hydrostatic internal gravity wave turbulence. The equation is equivalent to the one obtained by Caillol & Zeitlin (Dyn. Atmos. Oceans, vol. 32, issue 2, 2000, pp. 81–112), but it takes a canonical form. We show that it conserves the energy without involving the resonance condition in frequency, and look for the isotropic part of the steady, scale-invariant solutions. We provide a parametrization of the resonant manifold of non-hydrostatic internal gravity wave triadic interactions. This allows us to simplify the collision integral, and to evaluate the transfer coefficients of all triadic interactions. In the hydrostatic limit, our equation is equivalent to the Hamiltonian description of Lvov & Tabak (Phys. Rev. Lett., vol. 87, issue 16, 2001, 168501).

We analyze a dataset from a numerically simulated, temporally evolving turbulent wake (Zhou, Phys. Rev. Fluids, vol. 7, 2022, 104802) that exhibits spontaneous anisotropic layering under strong stratification, alongside significant spatiotemporal variability within the flow. This analysis focuses on the irreversible flux coefficient, $\varGamma$, defined as the ratio between turbulent potential and kinetic energy dissipation rates. We find that the volume-averaged $\varGamma$ initially rises, reaches a plateau between 0.45 and 0.49 when the layering dynamics become dominant, and then decreases as viscosity plays a larger role. These peak $\varGamma$ values are consistent with those from prior simulations under strongly stratified conditions. Such efficient mixing occurs when the Ozmidov to Thorpe length scale ratio is between 0.37 and 0.52, consistent with numerical and field data reported by Mashayek et al. (J. Fluid Mech., vol. 826, 2017, pp. 522–552). To account for the coexistence of dynamically distinct regions within the flow, we perform conditional sampling of $\varGamma$ against a locally defined gradient Richardson number, ${\textit {Ri}}$. This reveals a flux-gradient relation between $\varGamma$ and ${\textit {Ri}}$ that remains largely consistent over time. This relation features a large, approximately constant value of $\varGamma$ for ${\textit {Ri}}$ values greater than one, echoing the ‘constant-power’ scenario postulated by Balmforth et al. (J. Fluid Mech., vol. 355, 1998, pp. 329–358).

An experimental study has been conducted on the near wake of a 6 : 1 spheroid, in both uniform and stratified backgrounds. The pitch angle, $\theta$, was varied from $0^\circ \text { to }20^\circ$. When $\theta = 0^\circ$, stratification decreases the characteristic wake element spacing so a characteristic Strouhal number ($St$) increases from 0.32 to 0.4. However, a similar measure scaled on wake momentum thickness shows the wake spacing to converge on those measured for other bluff and streamlined bodies. There is an apparent effect of Reynolds number, which changes the location of separation lines and hence the initial wake thickness. When $\theta > 0^\circ$, the wake is a combination of the usual drag wake together with a collection of streamwise vortices that have separated from the body, and this wake geometry can evolve in ways that are measurably different from the zero incidence case. These differences may be limited to the near wake, as the later evolution appears to converge with previous bluff- and streamlined bodies, with normalised wake height, $L_V = 0.5$ and centreline velocity, $\bar {u}_0 = 0.3$ at $Nt = 10$, as the early wake enters the non-equilibrium regime with similar values to previously studied stratified wakes. In the presence of density stratification, the inclined wake itself generates large-scale internal wave undulations with time scale $2{\rm \pi} /N$, even when the background stratification is not strong and a body-based Froude number is $O(10)$. The geometry and strengths of the primary streamwise vortices are not symmetric, mirroring previous results from experiments and computations in the literature.

We investigate the dynamics of a columnar Taylor–Green vortex array under strong stratification, focusing on Froude numbers $0.125\leq Fr \leq 1.0$, with the aim of identifying and understanding the primary instabilities that lead to the vortices’ breakdown. Linear stability analysis reveals that the fastest-growing vertical wavenumber scales with $Fr^{-1}$, while the dimensionless growth rate remains approximately constant. The most unstable eigenmode, identified as the mixed hyperbolic mode by Hattori et al. (J. Fluid Mech., vol. 909, 2021, A4), bears significant similarities to the zigzag instability, first discovered by Billant & Chomaz (J. Fluid Mech., vol. 418, 2000, pp. 167–188). Direct numerical simulations further confirm that the zigzag instability is crucial in amplifying initial random perturbations to finite amplitude, with the flow structure and modal growth rate consistent with the linear stability analysis. In particular, the characteristic vertical length scale of turbulence matches that of the fastest-growing linear mode. These findings underscore the broader relevance of the zigzag instability mechanism beyond its initial discovery in vortex pairs, demonstrating its role in facilitating direct energy transfer from vertically uniform vortical motions to a characteristic vertical length scale proportional to $Fr$ in strongly stratified flows.

This study proposes a new mechanism that can lead to layering or convection from the finite amplitude perturbation acting on the double diffusive convection with uniform background shear. We focus on the double diffusive convection in the diffusive regime with the cold fresh water laying above the warm salty water. We demonstrate that, although the unperturbed system is linearly stable, the finite amplitude perturbation can trigger the initial flow motions which subsequently obtain energy from the gravitational potential energy and from the uniform background shear, and evolve to layering or convection. By using the linear stability analysis for the initial growth stage and the energy analysis for the following transitional stage, the critical Richardson number can be predicted theoretically. Here the Richardson number measures the relative strength of stratification to the background shear. The dominant wavenumbers and the growth rates of the corresponding modes given by linear theory agree well with the two-dimensional direct numerical simulations, and so does the critical Richardson number predicted by the theoretical model. The layering state is dominated by the double diffusion process, while the convection state at smaller Richardson number exhibits stronger influences from shear and generates smaller heat and salinity fluxes. The theoretical model is further applied to the parameter range which is relevant to the real oceanic environments and reveals that for the typical density ratio observed in the staircase regions in the Arctic Ocean, the current mechanism can lead to layering for relatively weak shear.

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.

Recent theoretical progress using multiscale asymptotic analysis has revealed various possible regimes of stratified turbulence. Notably, buoyancy transport can either be dominated by advection or diffusion, depending on the effective Péclet number of the flow. Two types of asymptotic models have been proposed, which yield measurably different predictions for the characteristic vertical velocity and length scale of the turbulent eddies in both diffusive and non-diffusive regimes. The first, termed a ‘single-scale model’, is designed to describe flow structures having large horizontal and small vertical scales, while the second, termed a ‘multiscale model’, additionally incorporates flow features with small horizontal scales, and reduces to the single-scale model in their absence. By comparing predicted vertical velocity scaling laws with direct numerical simulation data, we show that the multiscale model correctly captures the properties of strongly stratified turbulence within regions dominated by small-scale isotropic motions, whose volume fraction decreases as the stratification increases. Meanwhile its single-scale reduction accurately describes the more orderly, layer-like, quiescent flow outside those regions.

To date, a growing body of literature has documented the existence and impacts of coherent structures known as large- and very-large-scale motions within wall-bounded turbulent flows under neutral and unstable thermal stratification. These coherent structures can account for a considerable fraction of the overall turbulent transport and have been found to modulate small-scale turbulent fluctuations near the wall. In the context of stably stratified flows, however, the examination of such coherent structures has garnered relatively little attention. Stable stratification limits vertical transport and turbulent mixing within flows, which makes it unclear the extent to which previous findings on coherent structures under unstable and neutral stratification are applicable to stably stratified flows. In this study, we investigate the existence and characteristics of coherent structures under stable stratification with a wide range of statistical and spectral analyses. Outer peaks in premultiplied spectrograms under weak stability indicate the presence of large-scale motions, but these peaks become weaker and eventually vanish with increasing stability. Quadrant analysis of turbulent transport efficiencies (the ratio of net fluxes to their respective downgradient components) demonstrates dependencies on both stability and height above ground, which is evidence of morphological differences in the coherent structures under increasing stability. Amplitude modulation by large-scale streamwise velocity was found to decrease with increasing gradient Richardson number, whereas modulation by large-scale vertical velocity was approximately zero across all stability ranges. For sufficiently stable stratification, large eddies are suppressed enough to limit any inner–outer scale interactions.

New experiments in highly turbulent, steady, subcritical and uniform water open-channel flows have been carried out to measure the mean turbulent kinetic energy (TKE) budget of sediment-laden boundary layer flows with two sizes (dp = 3 mm and 1 mm) of Plexiglas particles $({\rm relative\ density}\ = 1.192)$. The experiments covered energetic sediment transport conditions (Shields number of $0.35 < \theta < 1.2$) ranging from non-capacity to full-capacity flows in bedload-to-suspension-dominated transport modes (suspension number of $0.5 < {w_s}/{u_\ast } < 1.3$ where ${w_s}$ is the settling velocity and $u_*$ is the friction velocity) and for weakly to highly inertial, finite size turbulence-particle conditions (Stokes number of 0.1 < St < 3.5 and dp/η > 10 where $\eta$ is the Kolmogorov length scale). It was shown that the effects of sediments on the TKE budget are very pronounced in all large particle experiments for which a bedload layer of several grain diameter thickness is developed above the channel bed. When compared with the corresponding reference clear-water flows, the TKE shear-production rate for the 3 mm particle flows is strongly reduced in the wall region corresponding to the bedload layer. This turbulence damping is seen to increase with sediment load until full capacity for flows with constant Shields value, as well as with Shields number value. Inside this damped TKE shear-production zone, a distinct peak of maximal turbulence production appears to coincide with the upper edge of the bedload layer delimited by a sharp gradient in mean sediment concentration. This vertically upshifted peak of TKE production is accompanied by an enhanced net downward oriented TKE flux when compared with the reference clear-water flows. The downward diffused TKE is found to act in the bedload layer as a local energy source in reasonable balance with the sediment transport term. The mechanism behind this downward TKE transport was further analysed on the basis of coherent flow structure dynamics controlled by ejection- and sweep-type events. The agreement between the height of downward directed mean TKE flux and the height below which sweep-type events dominate the Reynolds shear-stress contribution over ejections, revealed the leading role played by sweeps in mean TKE transport. This agreement holds for all reference clear-water flows supporting the well-known wall-roughness-induced dominance of the sweep contribution in turbulent, rough clear-water boundary layer flows. Furthermore, for all 3 mm particle flows, the two referred to transition levels were significantly and similarly upshifted to the upper edge of the bedload layer. Only for these sediment-laden flows, the bedload layer thickness is seen to exceed the wall-roughness sublayer of the reference clear-water flows. This supports a strong analogy between wall-roughness effects in clear-water flows and bedload layer effects in sediment-laden flows, on the mean TKE budget induced by a similarly modified coherent flow structure dynamics. The bedload layer-controlled wall roughness is finally confirmed by the good prediction of the wall-roughness parameter ks of the logarithmic velocity distribution. An empirical formulation fitting the presented measurements is presented, valid over the range of Shields number values covered herein.

Studies of Kelvin–Helmholtz (KH) instability have typically modelled the initial mean flow as an isolated stratified shear layer. However, geophysical flows frequently exhibit multiple layers. As a step towards understanding these flows, we examine the case of two adjacent stratified shear layers using both linear stability analysis and direct numerical simulations. With sufficiently large layer separation, the characteristics of instability and mixing converge towards the familiar KH turbulence, and similarly when the separation is near zero and the layers add to make a single layer, albeit with a reduced Richardson number. Here, our focus is on intermediate separations, which produce new and complex phenomena. As the separation distance $D$ increases from zero to a critical value $D_c$, approximately half the thickness of the shear layer, the growth rate and wavenumber both decrease monotonically. The minimum Richardson number is relatively low, potentially inducing pairing, and shear-aligned convective instability (SCI) is the primary mechanism for transition. Consequently, mixing is relatively strong and efficient. When $D\sim D_c$, billow length is increased but growth is slowed. Despite the modest growth rate, mixing is strong and efficient, engendered primarily by secondary shear instability manifested on the braids, and by SCI occurring on the eyelids. Shear-aligned vortices are driven in part by buoyancy production; however, shear production and vortex stretching are equally important mechanisms. When $D>D_c$, neighbouring billow interactions suppress the growth of both KH instability and SCI. Strength and efficiency of mixing decrease abruptly as $D_c$ is exceeded. As turbulence decays, layers of marginal instability may arise.

We present evidence for layered anisotropic stratified turbulence (LAST) and mixing produced in a freely evolving uniformly stratified shear layer where the direction of shear is orthogonal to gravity. As originally reported by Basak & Sarkar (J. Fluid. Mech., vol. 568, 2006, pp. 19–54), such a flow develops a rich three-dimensional structure in the form of interlocking columnar vortices formed by horizontal shear instability that remain coherent at large scales due to the stabilising vertical stratification. Here, we modify the initial velocity field by introducing additional small-amplitude vertical perturbations designed to be representative of pre-existing horizontal layers often observed in strongly stratified ocean environments. This reveals a novel finite amplitude, non-normal mode growth mechanism through which the vertical shear between layers may be rapidly amplified by its interaction with the horizontal shear layer prior to the growth of shear instability, leading to a rapid turbulent transition instigated by the subsequent interaction of the layers with the emerging columnar vortices. Through a consideration of relevant flow statistics and associated dimensionless parameters, we demonstrate that turbulence can enter the LAST regime, thereby indicating a generic mechanism leading to the transient development of regions of strongly stratified turbulence in the ocean. We discuss the properties of mixing and the parameterisation of mixing efficiency in terms of the relationship between turbulent length scales in the flow, in particular highlighting links to models based on the classical vertical shear instability paradigm typically associated with more weakly stratified flows that produce isolated turbulent ‘patches’.

Relatively strongly stratified turbulent flows tend to self-organise into a ‘layered anisotropic stratified turbulence’ (LAST) regime, characterised by relatively deep and well-mixed density ‘layers’ separated by relatively thin ‘interfaces’ of enhanced density gradient. Understanding the associated mixing dynamics is a central problem in geophysical fluid dynamics. It is challenging to study LAST mixing, as it is associated with Reynolds numbers $Re := UL/\nu \gg 1$ and Froude numbers $Fr :=(2{\rm \pi} U)/(L N) \ll 1$ ($U$ and $L$ being characteristic velocity and length scales, $\nu$ the kinematic viscosity and $N$ the buoyancy frequency). Since a sufficiently large dynamic range (largely) unaffected by stratification and viscosity is required, it is also necessary for the buoyancy Reynolds number $Re_{b} := \epsilon /(\nu N^{2}) \gg 1$, where $\epsilon$ is the (appropriately volume-averaged) turbulent kinetic energy dissipation rate. This requirement is exacerbated for oceanically relevant flows, as the Prandtl number $Pr := \nu /\kappa = {O}(10)$ in thermally stratified water (where $\kappa$ is the thermal diffusivity), thus leading (potentially) to even finer density field structures. We report here on four forced fully resolved direct numerical simulations of stratified turbulence at various Froude ($Fr=0.5, 2$) and Prandtl ($Pr=1, 7$) numbers forced so that $Re_{b}=50$, with resolutions up to $30\,240 \times 30\,240 \times 3780$. We find that, as $Pr$ increases, emergent ‘interfaces’ become finer and their contribution to bulk mixing characteristics decreases at the expense of the small-scale density structures populating the well-mixed ‘layers’. However, extreme mixing events (as quantified by significantly elevated local destruction rates of buoyancy variance $\chi _0$) are always preferentially found in the (statically stable) interfaces, irrespective of the value of $Pr$.

Can the similarity between rotation and stratification provide quantitative predictions for complex turbulence? In this work, we focus on plane Couette turbulence as the background flow, which allows us to eliminate the key differences between Rayleigh–Bénard and Taylor–Couette turbulence, and facilitates a quantitative mapping across many complex turbulent systems involving rotation, stratification and curvature effects. To characterize the separated or coupled effects of rotation and stratification, we introduce an overall Richardson number ${Ri}_\chi$ which is the sum of the Coriolis Richardson number $Ri_\theta$ and the buoyancy Richardson number $Ri_T$. When the Prandtl number $Pr=1$, the heat and inertial-frame momentum transport almost coincide and mainly depend on ${Ri}_\chi$ and the Reynolds number, regardless of the specific ratio between $Ri_\theta$ and $Ri_T$. When $Pr$ varies, the weighted average ${Nu}_\chi$ of the transport coefficients can still remain approximately invariant in most cases. Furthermore, the large-scale structures in purely rotating and purely stratified cases exhibit strikingly similar features. These findings not only indicate a more reasonable analogy for the ultimate Taylor–Couette turbulence but also pave the way for developing new predictive models for natural and industrial processes.

This article concerns the time-evolution, spectral structure and scaling of weak turbulence subject to rotation and stable stratification. The flow is expressed as a combination of particular solutions, referred to as modes, of the linearised governing equations without viscosity or diffusion. Modes are of two types: oscillatory ones which represent inertial-gravity waves and time-independent ones that express a non-propagating (NP) component of the flow. The presence of the NP component, which plays an active role in the dynamics apart from in the case of pure rotation, renders wave-turbulence analysis problematic because the NP mode is non-dispersive. Equations are derived for the time evolution of the modal amplitudes, evolution which is due to nonlinearity and visco-diffusion. Subsequent analysis assumes that one or other (or both) of the Rossby and Froude numbers is small (weak turbulence). Given this assumption, the NP component is found to evolve independently of the wave one and a numerical scheme, similar to, though significantly different from classical direct numerical simulation, is used to determine its time evolution. The treatment of the wave component assumes its amplitude is large compared with the NP one, otherwise there are seemingly intractable difficulties of closure in the analysis. Given this further assumption, the wave component decouples from the NP one. Evolution equations for the wave spectra are derived using wave-turbulence analysis and are integrated numerically. As might be expected, these equations indicate that nonlinear coupling of wave modes is dominated by resonances. Results are given for both the NP and wave components.

The presence of a mean magnetic field aligned with the direction of the acceleration greatly modifies the development of the Rayleigh–Taylor instability (RTI). High resolution direct numerical simulations of the Boussinesq–Navier–Stokes equations under the magnetohydrodynamics approximation reveal that, after an initial damping of the perturbations at the interface between the two miscible fluids, a rapid increase of the mixing layer is observed. Structures are significantly stretched in the vertical direction because magnetic tension prevents small-scale shear instabilities. When the vertical turbulent velocity exceeds the Alfvén velocity, the flow transitions to turbulence, structures break and an enhanced mixing occurs with strong dissipation. Afterwards, the mixing zone slows down and its growth rate is decreased compared to the hydrodynamic case. For larger magnitudes of the mean magnetic field, a strong anisotropy persists, and an increased fraction of potential energy injected into the system is lost into turbulent magnetic energy: as a consequence, the mixing zone growth rate is decreased even more. This phenomenology is embedded in a general buoyancy-drag equation, derived from simplified equations that reflect the large-scale dynamics, in which the drag coefficient is increased by the presence of turbulent magnetic energy.

We introduce a geometric analysis of turbulent mixing in density-stratified flows based on the alignment of the density gradient in two orthogonal bases that are locally constructed from the velocity gradient tensor. The first basis connects diapycnal mixing to rotation and shearing motions, building on the recent ‘rortex–shear decomposition’ in stratified shear layers (Jiang et al., J. Fluid Mech., vol. 947, 2022, A30), while the second basis connects mixing to the principal axes of the viscous dissipation tensor. Applying this framework to datasets taken in the stratified inclined duct laboratory experiment reveals that density gradients in locations of high shear tend to align preferentially (i) along the direction of minimum dissipation and (ii) normal to the plane spanned by the rortex and shear vectors. The analysis of the local alignment across increasingly turbulent flows offers new insights into the intricate relationship between the density gradient and dissipation, and thus diapycnal mixing.