We consider the effect of stratification on systematic, large-scale flows generated in anelastic convection. We present results from three-dimensional numerical simulations of convection in a rotating plane layer in which the angle between the axis of rotation and gravity is allowed to vary. This model is representative of different latitudes of a spherical body. We consider two distinct parameter regimes: (i) weakly rotating and (ii) rapidly rotating. In each case, we examine the effect of stratification on the flow structure and heat transport properties focussing on the difference between Boussinesq and anelastic convection. Furthermore, we show that regimes (i) and (ii) generate very different large-scale flows and we investigate the role stratification has in modifying these flows. The stratified flows possess a net helicity not present in the Boussinesq cases which we suggest, when combined with the self-generated shear flows, could be important for dynamo action.

We continue our investigation into the nonlinear evolution of the Goldreich-Schubert-Fricke (GSF) instability in differentially rotating radiation zones. This instability may be a key player in transporting angular momentum in stars and giant planets, but its nonlinear evolution remains mostly unexplored. In a previous paper we considered the equatorial instability, whereas here we simulate the instability at a general latitude for the first time. We adopt a local Cartesian Boussinesq model in a modified shearing box for most of our simulations, but we also perform some simulations with stress-free, impenetrable, radial boundaries. We first revisit the linear instability and derive some new results, before studying its nonlinear evolution. The instability is found to behave very differently compared with its behaviour at the equator. In particular, here we observe the development of strong zonal jets ("layering" in the angular momentum), which can considerably enhance angular momentum transport, particularly in axisymmetric simulations. The jets are, in general, tilted with respect to the local gravity by an angle that corresponds initially with that of the linear modes, but which evolves with time and depends on the strength of the flow. The instability transports angular momentum much more efficiently (by several orders of magnitude) than it does at the equator, and we estimate that the GSF instability could contribute to the missing angular momentum transport required in both red giant and subgiant stars. It could also play a role in the long-term evolution of the solar tachocline and the atmospheric dynamics of hot Jupiters.

Classification of matter through topological phases and topological edge states between distinct materials has been a subject of great interest recently. While lattices have been the main setting for these studies, a relatively unexplored realm for this physics is that of continuum fluids. In the typical case of a fluid model with a point spectrum, nontrivial topology and associated edge modes have been observed previously. However, another possibility is that a continuous spectrum can conexist with the point spectrum. Here we demonstrate that a fluid plasma model can harbor nontrivial topology within its continuous spectrum, and that there are boundary modes at the interface between topologically distinct regions. We consider the ideal magnetohydrodynamics (MHD) model. In the presence of magnetic shear, we find nontrivial topology in the {\Alfven} continuum. For strong shear, the Chern number is $\pm 1$, depending on the sign of the shear. If the magnetic shear changes sign within the plasma, a topological phase transition occurs, and bulk-boundary correspondence then suggests a mode localized to the layer of zero magnetic shear. We confirm the existence of this mode numerically. Moreover, this reversed-shear {\Alfven} eigenmode (RSAE) is well known within the magnetic fusion community as it has been observed in several tokamaks. A potential implication of our findings is that the RSAE may be topologically robust to three-dimensional perturbations such as magnetic islands. Moreover, this work raises the broader possibility that fusion devices could be carefully tailored to produce other topological states and potentially useful behavior. In considering how the MHD model might be regularized at small scales, we also consider the electron limit of Hall MHD. We show that the whistler band, which plays an important role in planetary ionospheres, has nontrivial topology.

Direct Statistical Simulation (DSS) solves the equations of motion for the statistics of turbulent flows in place of the traditional route of accumulating statistics by Direct Numerical Simulation (DNS). That low-order statistics usually evolve slowly compared with instantaneous dynamics is one important advantage of DSS. Depending on the symmetry of the problem and the choice of averaging operation, however, DSS is usually more expensive computationally than DNS because even low order statistics typically have higher dimension than the underlying fields. Here we show that it is possible to go much further by using Proper Orthogonal Decomposition (POD) to address the "curse of dimensionality." We apply POD directly to DSS in the form of expansions in the equal-time cumulants to second order (CE2). We explore two averaging operations (zonal and ensemble) and test the approach on two idealized barotropic models on a rotating sphere (a jet that relaxes deterministically towards an unstable profile, and a stochastically-driven flow that spontaneously organizes into jets). Order-of-magnitude savings in computational cost are obtained in the reduced basis, potentially enabling access to parameter regimes beyond the reach of DNS.

The generation of magnetic field in an electrically conducting fluid generally involves the complicated nonlinear interaction of flow turbulence, rotation and field. This dynamo process is of great importance in geophysics, planetary science and astrophysics, since magnetic fields are known to play a key role in the dynamics of these systems. This paper gives an introduction to dynamo theory for the fluid dynamicist. It proceeds by laying the groundwork, introducing the equations and techniques that are at the heart of dynamo theory, before presenting some simple dynamo solutions. The problems currently exercising dynamo theorists are then introduced, along with the attempts to make progress. The paper concludes with the argument that progress in dynamo theory will be made in the future by utilising and advancing some of the current breakthroughs in neutral fluid turbulence such as those in transition, self-sustaining processes, turbulence/mean-flow interaction, statistical methods and maintenance and loss of balance.

We present an investigation into the nonlinear evolution of the Goldreich-Schubert-Fricke (GSF) instability using axisymmetric and three-dimensional simulations near the equator of a differentially rotating radiation zone. This instability may provide an important contribution to angular momentum transport in stars and planets. We adopt a local Boussinesq Cartesian shearing box model, which represents a small patch of a differentially rotating stellar radiation zone. Complementary simulations are also performed with stress-free, impenetrable boundaries in the local radial direction. The linear and nonlinear evolution of the equatorial axisymmetric instability is formally equivalent to the salt fingering instability. This is no longer the case in 3D, but we find that the instability behaves nonlinearly in a similar way to salt fingering. Axisymmetric simulations -- and those in 3D with short dimensions along the local azimuthal direction -- quickly develop strong jets along the rotation axis, which inhibit the instability and lead to predator-prey-like temporal dynamics. In 3D, the instability initially produces homogeneous turbulence and enhanced momentum transport, though in some cases jets form on a much longer timescale. We propose and validate numerically a simple theory for nonlinear saturation of the GSF instability and its resulting angular momentum transport. This theory is straightforward to implement in stellar evolution codes incorporating rotation. We estimate that the GSF instability could contribute towards explaining the missing angular momentum transport required in red giant stars, and play a role in the long-term evolution of the solar tachocline.

Global magnetohydrodynamic (MHD) instabilities are investigated in a computationally tractable two-dimensional model of the solar tachocline. The model's differential rotation yields stability in the absence of a magnetic field, but if a magnetic field is present, a joint instability is observed. We analyze the nonlinear development of the instability via fully nonlinear direct numerical simulation, the generalized quasilinear approximation (GQL), and direct statistical simulation (DSS) based upon low-order expansion in equal-time cumulants. As the magnetic diffusivity is decreased, the nonlinear development of the instability becomes more complicated until eventually a set of parameters are identified that produce a previously unidentified long-term cycle in which energy is transformed from kinetic energy to magnetic energy and back. We find that the periodic transitions, which mimic some aspects of solar variability -- for example, the quasiperiodic seasonal exchange of energy between toroidal field and waves or eddies -- are unable to be reproduced when eddy-scattering processes are excluded from the model.

We consider rotating, kinematic dynamos at low magnetic Prandtl number $Pm$. We show that the inclusion of rotation leads to an increase in spatio-temporal coherence and a modification of the turbulent spectrum. These effects make the flow more efficient in driving the dynamo, in the sense that the energy injection rate required to reach the critical value of the magnetic Reynolds number $Rm_c$ is reduced in comparison with a non-rotating dynamo (Seshasayanan et al. 2017). For random dynamos it is known that the growth-rate would largely be determined by the spectral index of the flow at the resistive scale. Here, however, we demonstrate that the dynamo growth-rate in rotating flows is increased by the rotationally induced long-lived large scale eddies with a coherence time greater than the local turnover time. These eddies play the major role in determining the dynamo growth-rate.

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 Generalised Quasilinear Approximation (GQL) (Marston et al 2016) may provide a much better approximation to the complicated full nonlinear dynamics than the widely used Quasilinear Approximation (QL). This result will enable the construction of more accurate statistical theories for the description of geophysical and astrophysical flows.

In this paper we examine the role of weak magnetic fields in breaking Kelvin's circulation theorem and in vortex breakup in two-dimensional magnetohydrodynamics for the physically important case of a low magnetic Prandtl number (low $Pm$) fluid. We consider three canonical inviscid solutions for the purely hydrodynamical problem, namely a Gaussian vortex, a circular vortex patch and an elliptical vortex patch. We examine how magnetic fields lead to an initial loss of circulation $\Gamma$ and attempt to derive scaling laws for the loss of circulation as a function of field strength and diffusion as measured by two non-dimensional parameters. We show that for all cases the loss of circulation depends on the integrated effects of the Lorentz force, with the patch cases leading to significantly greater circulation loss. For the case of the elliptical vortex the loss of circulation depends on the total area swept out by the rotating vortex and so this leads to more efficient circulation loss than for a circular vortex.

We examine the role of stratification in determining the scale for turbulent anelastic convection. Motivated by the range of scales observed in convection at the solar photosphere, we perform local numerical simulations of convection for a range of density contrasts in large domains. We analyse both the Eulerian and Lagrangian statistics of the convection and demonstrate that increasing the stratification shifts the scale of the most energetic structures in the flow to smaller scales; furthermore, the relative amplitude of vertical to horizontal flows in the convection decreases with increasing stratification. We discuss the implications of our results to the issue of solar mesogranulation.

Compositional convection is thought to be an important energy source for magnetic field generation within planetary interiors. The Prandtl number, $Pr$, characterizing compositional convection is significantly larger than unity, suggesting that the inertial force may not be important on the small scales of convection as long as the buoyancy force is not too strong. We develop asymptotic dynamo models for the case of small Rossby number and large Prandtl number in which inertia is absent on the convective scale. The relevant diffusivity parameter for this limit is the compositional Roberts number, $q = D/\eta$, which is the ratio of compositional and magnetic diffusivities. Dynamo models are developed for both order one $q$ and the more geophysically relevant low $q$ limit. For both cases the ratio of magnetic to kinetic energy densities, $M$, is asymptotically large and reflects the fact that Alfv\'en waves have been filtered from the dynamics. Along with previous investigations of asymptotic dynamo models for $Pr=O(1)$, our results show that the ratio $M$ is not a useful indicator of dominant force balances in the momentum equation since many different asymptotic limits of $M$ can be obtained without changing the leading order geostrophic balance. Furthermore, the present models show that inertia is not a requirement for driving low $q$, large-scale dynamos.

We consider dynamo action driven by three-dimensional rotating anelastic convection in a spherical shell. Motivated by the behaviour of the solar dynamo, we examine the interaction of hydromagnetic modes with different symmetries and demonstrate how complicated interactions between convection, differential rotation and magnetic fields may lead to modulation of the basic cycle. For some parameters, Type 1 modulation occurs by the transfer of energy between modes of different symmetries with little change in the overall amplitude, for other parameters, the modulation is of Type 2, where the amplitude is significantly affected (leading to grand minima in activity) without significant changes in symmetry. Most importantly, we identify the presence of "supermodulation" in the solutions, where the activity switches chaotically between Type 1 and Type 2 modulation, this is believed to be an important process in solar activity.

We review progress that has been made in utilizing one form of Direct Statistical Simulation (DSS) to describe geophysical and astrophysical flows that are anisotropic and inhomogeneous. We first explain the approach, which is based upon a systematic and conservative expansion of the equations of motion for low-order equal-time cumulants. We place the method into context with other statistical procedures. Truncation at second order in the hierarchy of cumulants is equivalent to retaining the interaction between zonal mean flows and eddies. Eddy-eddy interactions appear at higher orders, but care must be taken to keep the higher-order expansions realizable with non-negative probability distribution functions. The strengths and weaknesses of different levels of approximation are assessed with numerical experiments on the fiducial problem of a stochastically forced jet on a spherical surface. The results give an insight into the mechanisms that may control jet spacing and strength, and indicate interesting avenues for future research.

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).

The onset of dynamo action is investigated within the context of a newly developed low Rossby, low magnetic Prandtl number, convection-driven dynamo model. This multiscale model represents an asymptotically exact form of an $\alpha^2$ mean field dynamo model in which the small-scale convection is represented explicitly by finite amplitude, single mode solutions. Both steady and oscillatory convection are considered for a variety of horizontal planforms. The kinetic helicity is observed to be a monotonically increasing function of the Rayleigh number. As a result, very small magnetic Prandtl number dynamos can be found for sufficiently large Rayleigh numbers. All dynamos are found to be oscillatory with an oscillation frequency that increases as the strength of the convection is increased and the magnetic Prandtl number is reduced. Kinematic dynamo action is strongly controlled by the profile of the helicity; single mode solutions which exhibit boundary layer behavior in the helicity show a decrease in the efficiency of dynamo action due to the enhancement of magnetic diffusion in the boundary layer regions. For a given value of the Rayleigh number, lower magnetic Prandtl number dynamos are excited for the case of oscillatory convection in comparison to steady convection. With regards to planetary dynamos, these results suggest that the low magnetic Prandtl number dynamos typical of liquid metals are more easily driven by thermal convection than by compositional convection.

Quasilinear theory is often utilized to approximate the dynamics of fluids exhibiting significant interactions between mean flows and eddies. In this paper we present a generalization of quasilinear theory to include dynamic mode interactions on the large scales. This generalized quasilinear (GQL) approximation is achieved by separating the state variables into large and small zonal scales via a spectral filter rather than by a decomposition into a formal mean and fluctuations. Nonlinear interactions involving only small zonal scales are then removed. The approximation is conservative and allows for scattering of energy between small-scale modes via the large scale (through non-local spectral interactions). We evaluate GQL for the paradigmatic problems of the driving of large-scale jets on a spherical surface and on the beta-plane and show that it is accurate even for a small number of large-scale modes. As this approximation is formally linear in the small zonal scales it allows for the closure of the system and can be utilized in direct statistical simulation schemes that have proved an attractive alternative to direct numerical simulation for many geophysical and astrophysical problems.

We live in an age in which high-performance computing is transforming the way we do science. Previously intractable problems are now becoming accessible by means of increasingly realistic numerical simulations. One of the most enduring and most challenging of these problems is turbulence. Yet, despite these advances, the extreme parameter regimes encountered in space physics and astrophysics (as in atmospheric and oceanic physics) still preclude direct numerical simulation. Numerical models must take a Large Eddy Simulation (LES) approach, explicitly computing only a fraction of the active dynamical scales. The success of such an approach hinges on how well the model can represent the subgrid-scales (SGS) that are not explicitly resolved. In addition to the parameter regime, heliophysical and astrophysical applications must also face an equally daunting challenge: magnetism. The presence of magnetic fields in a turbulent, electrically conducting fluid flow can dramatically alter the coupling between large and small scales, with potentially profound implications for LES/SGS modeling. In this review article, we summarize the state of the art in LES modeling of turbulent magnetohydrodynamic (MHD) flows. After discussing the nature of MHD turbulence and the small-scale processes that give rise to energy dissipation, plasma heating, and magnetic reconnection, we consider how these processes may best be captured within an LES/SGS framework. We then consider several specific applications in heliophysics and astrophysics, assessing triumphs, challenges, and future directions.

We investigate the intermittency of energy dissipation in magnetohydrodynamic (MHD) turbulence by identifying dissipative structures and measuring their characteristic scales. We find that the probability distribution of energy dissipation rates exhibits a power law tail with index very close to the critical value of -2.0, which indicates that structures of all intensities contribute equally to energy dissipation. We find that energy dissipation is uniformly spread among coherent structures with lengths and widths in the inertial range. At the same time, these structures have thicknesses deep within the dissipative regime. As the Reynolds number is increased, structures become thinner and more numerous, while the energy dissipation continues to occur mainly in large-scale coherent structures. This implies that in the limit of high Reynolds number, energy dissipation occurs in thin, tightly packed current sheets which nevertheless span a continuum of scales up to the system size, exhibiting features of both coherent structures and nanoflares previously conjectured as a coronal heating mechanism.

We consider the generation of magnetic activity --- dynamo waves --- in the astrophysical limit of very large magnetic Reynolds number. We consider kinematic dynamo action for a system consisting of helical flow and large-scale shear. We demonstrate that large-scale dynamo waves persist at high $Rm$ if the helical flow is characterised by a narrow band of spatial scales and the shear is large enough. However for a wide band of scales the dynamo becomes small-scale with a further increase of $Rm$, with dynamo waves re-emerging only if the shear is then increased. We show that at high $Rm$ the key effect of the shear is to suppress small-scale dynamo action, allowing large-scale dynamo action to be observed. We conjecture that this supports a general "suppression principle" --- large-scale dynamo action can only be observed if there is a mechanism that suppresses the small-scale fluctuations.

The predominant force balance in rapidly rotating planetary cores is between Coriolis, pressure, buoyancy and Lorentz forces. This magnetostrophic balance leads to a Taylor state where the spatially averaged azimuthal Lorentz force is compelled to vanish on cylinders aligned with the rotation axis. Any deviation from this state leads to a torsional oscillation, signatures of which have been observed in the Earth's secular variation and are thought to influence length of day variations via angular momentum conservation. In order to investigate the dynamics of torsional oscillations, we perform several three-dimensional dynamo simulations in a spherical shell. We find torsional oscillations, identified by their propagation at the correct Alfv\'{e}n speed, in many of our simulations. We find that the frequency, location and direction of propagation of the waves are influenced by the choice of parameters. Torsional waves are observed within the tangent cylinder and also have the ability to pass through it. Several of our simulations display waves with core travel times of 4 to 6 years. We calculate the driving terms for these waves and find that both the Reynolds force and ageostrophic convection acting through the Lorentz force are important in driving torsional oscillations.

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.

A review of MHD dynamos and turbulence.

In this paper we introduce the concept of Direct Statistical Simulation (DSS) for astrophysical flows. This technique may be appropriate for problems in astrophysical fluids where the instantaneous dynamics of the flows are of secondary importance to their statistical properties. We give examples of such problems including mixing and transport in planets, stars and disks. The method is described for a general set of evolution equations, before we consider the specific case of a spectral method optimised for problems on a spherical surface. The method is illustrated for the simplest non-trivial example of hydrodynamics and MHD on a rotating spherical surface. We then discuss possible extensions of the method both in terms of computational methods and the range of astrophysical problems that are of interest.

The operation of an interface dynamo (as has been suggested for the Sun and other stars with convective envelopes) relies crucially upon the effective transport of magnetic flux between two spatially disjoint generation regions. In the simplest models communication between the two regions is achieved solely by diffusion. Here we incorporate a highly simplified anisotropic transport mechanism in order to model the net effect of flux conveyance by magnetic pumping and by magnetic buoyancy. We investigate the influence of this mechanism on the efficiency of kinematic dynamo action. It is found that the effect of flux transport on the efficiency of the dynamo is dependent upon the spatial profile of the transport. Typically, transport hinders the onset of dynamo action and increases the frequency of the dynamo waves. However, in certain cases there exists a preferred magnitude of transport for which dynamo action is most efficient. Furthermore, we demonstrate the importance of the imposition of boundary conditions in drawing conclusions on the role of transport.

Motivated by the problem of the formation of active regions from a deep-seated solar magnetic field, we consider the nonlinear three-dimensional evolution of magnetic buoyancy instabilities resulting from a smoothly stratified horizontal magnetic field. By exploring the case for which the instability is continuously driven we have identified a new mechanism for the formation of concentrations of magnetic flux.

We discuss the difficulties of predicting the solar cycle using mean-field models. Here we argue that these difficulties arise owing to the significant modulation of the solar activity cycle, and that this modulation arises owing to either stochastic or deterministic processes. We analyse the implications for predictability in both of these situations by considering two separate solar dynamo models. The first model represents a stochastically-perturbed flux transport dynamo. Here even very weak stochastic perturbations can give rise to significant modulation in the activity cycle. This modulation leads to a loss of predictability. In the second model, we neglect stochastic effects and assume that generation of magnetic field in the Sun can be described by a fully deterministic nonlinear mean-field model -- this is a best case scenario for prediction. We designate the output from this deterministic model (with parameters chosen to produce chaotically modulated cycles) as a target timeseries that subsequent deterministic mean-field models are required to predict. Long-term prediction is impossible even if a model that is correct in all details is utilised in the prediction. Furthermore, we show that even short-term prediction is impossible if there is a small discrepancy in the input parameters from the fiducial model. This is the case even if the predicting model has been tuned to reproduce the output of previous cycles. Given the inherent uncertainties in determining the transport coefficients and nonlinear responses for mean-field models, we argue that this makes predicting the solar cycle using the output from such models impossible.