Reduced-Order Quasilinear Model of Ocean Boundary-Layer Turbulence

The combined effectiveness of model reduction and the quasilinear approximation for the reproduction of the low-order statistics of oceanic surface boundary layer turbulence is investigated. Idealized horizontally homogeneous problems of surface-forced thermal convection and Langmuir turbulence are studied in detail. Model reduction is achieved with a Galerkin projection of the governing equations onto a subset of modes determined by proper orthogonal decomposition (POD). When applied to boundary layers that are horizontally homogeneous, POD and a horizontally averaged quasilinear approximation both assume flow features that are horizontally wavelike, making the pairing very efficient. For less than 0.2% of the modes retained, the reduced quasilinear model is able to reproduce vertical profiles of horizontal mean fields as well as certain energetically important second-order turbulent transport statistics and energies to within 30% error. Reduced-basis quasilinear statistics must approach the full-basis statistics as the basis size approaches completion; however, some quasilinear statistics resemble those found in the fully nonlinear simulations at smaller basis truncations. Thus, model reduction could possibly improve upon the accuracy of quasilinear dynamics.

Wave kinetic equation for inhomogeneous drift-wave turbulence beyond the quasilinear approximation

The formation of zonal flows from inhomogeneous drift-wave (DW) turbulence is often described using statistical theories derived within the quasilinear approximation. However, this approximation neglects wave–wave collisions. Hence, some important effects such as the Batchelor–Kraichnan inverse-energy cascade are not captured within this approach. Here we derive a wave kinetic equation that includes a DW collision operator in the presence of zonal flows. Our derivation makes use of the Weyl calculus, the quasinormal statistical closure and the geometrical-optics approximation. The obtained model conserves both the total enstrophy and energy of the system. The derived DW collision operator breaks down at the Rayleigh–Kuo threshold. This threshold is missed by homogeneous-turbulence theory but expected from a full-wave quasilinear analysis. In the future, this theory might help better understand the interactions between drift waves and zonal flows, including the validity domain of the quasilinear approximation that is commonly used in the literature.

Generalized quasilinear approximation of the interaction of convection and mean flows in a thermal annulus

In this paper, we examine the interaction of convection, rotation and mean flows in a thermal annulus. In this system, mean flows are driven by correlations induced by rotation leading to non-trivial Reynolds stresses. The mean flows act back on the convective turbulence acting as a barrier to transport. For this system, we demonstrate that the generalized quasilinear approximation (Marston et al . 2016 Phys. Rev. Lett. 116 , 214501. ( doi:10.1103/PhysRevLett.116.214501 )) may provide a much better approximation to the complicated full nonlinear dynamics than the widely used quasilinear approximation. This result will enable the construction of more accurate statistical theories for the description of geophysical and astrophysical flows.

On non-local energy transfer via zonal flow in the Dimits shift

The two-dimensional Terry–Horton equation is shown to exhibit the Dimits shift when suitably modified to capture both the nonlinear enhancement of zonal/drift-wave interactions and the existence of residual Rosenbluth–Hinton states. This phenomenon persists through numerous simplifications of the equation, including a quasilinear approximation as well as a four-mode truncation. It is shown that the use of an appropriate adiabatic electron response, for which the electrons are not affected by the flux-averaged potential, results in an $\boldsymbol{E}\times \boldsymbol{B}$ nonlinearity that can efficiently transfer energy non-locally to length scales of the order of the sound radius. The size of the shift for the nonlinear system is heuristically calculated and found to be in excellent agreement with numerical solutions. The existence of the Dimits shift for this system is then understood as an ability of the unstable primary modes to efficiently couple to stable modes at smaller scales, and the shift ends when these stable modes eventually destabilize as the density gradient is increased. This non-local mechanism of energy transfer is argued to be generically important even for more physically complete systems.

Three-dimensional rotating Couette flow via the generalised quasilinear approximation

We examine the effectiveness of the generalised quasilinear (GQL) approximation introduced by Marston et al. (Phys. Rev. Lett., vol. 116 (21), 2016, 214501). This approximation splits the variables into large and small scales in directions where there is a translational symmetry and removes nonlinear interactions involving only small scales. We utilise as a paradigm problem three-dimensional, turbulent, rotating Couette flow. We compare the results obtained from direct numerical solution of the equations with those from quasilinear (QL) and GQL calculations. In this three-dimensional setting, there is a choice of cutoff wavenumber for the GQL approximation both in the streamwise and in the spanwise directions. We demonstrate that the GQL approximation significantly improves the accuracy of mean flows, spectra and two-point correlation functions over models that are quasilinear in any of the translationally invariant directions, even if only a few streamwise and spanwise modes are included. We argue that this provides significant support for a programme of direct statistical simulation utilising the GQL approximation.