Motivated by recent advances in Direct Statistical Simulation (DSS) of astrophysical phenomena such as out of equilibrium jets, we perform a Direct Numerical Simulation (DNS) of the helical magnetorotational instability (HMRI) under the generalised quasilinear approximation (GQL). This approximation generalises the quasilinear approximation (QL) to include the self-consistent interaction of large-scale modes, interpolating between fully nonlinear DNS and QL DNS whilst still remaining formally linear in the small scales. In this paper we address whether GQL can more accurately describe low-order statistics of axisymmetric HMRI when compared with QL by performing DNS under various degrees of GQL approximation. We utilise various diagnostics, such as energy spectra in addition to first and second cumulants, for calculations performed for a range of Reynolds and Hartmann numbers (describing rotation and imposed magnetic field strength respectively). We find that GQL performs significantly better than QL in describing the statistics of the HMRI even when relatively few large-scale modes are kept in the formalism. We conclude that DSS based on GQL (GCE2) will be significantly more accurate than that based on QL (CE2).

We present a Direct Statistical Simulation (DSS) of jet formation on a \beta-plane, solving for the statistics of a fluid flow via an expansion in cumulants. Here we compare an expansion truncated at second order (CE2) to statistics accumulated by direct numerical simulations (DNS). We show that, for jets near equilibrium, CE2 is capable of reproducing the jet structure (although some differences remain in the second cumulant). However as the degree of departure from equilibrium is increased (as measured by the zonostrophy parameter) the jets meander more and CE2 becomes less accurate. We discuss a possible remedy by inclusion of higher cumulants.