Search

Numerical simulations of quasi-static magnetoconvection with a vertical magnetic field are carried out up to a Chandrasekhar number of $Q=10^8$ over a broad range of Rayleigh numbers $Ra$. Three magnetoconvection regimes are identified: two of the regimes are magnetically-constrained in the sense that a leading-order balance exists between the Lorentz and buoyancy forces, whereas the third regime is characterized by unbalanced dynamics that is similar to non-magnetic convection. Each regime is distinguished by flow morphology, momentum and heat equation balances, and heat transport behavior. One of the magnetically-constrained regimes appears to represent an `ultimate' magnetoconvection regime in the dual limit of asymptotically-large buoyancy forcing and magnetic field strength; this regime is characterized by an interconnected network of anisotropic, spatially-localized fluid columns aligned with the direction of the imposed magnetic field that remain quasi-laminar despite having large flow speeds. As for non-magnetic convection, heat transport is controlled primarily by the thermal boundary layer. Empirically, the scaling of the heat transport and flow speeds with $Ra$ appear to be independent of the thermal Prandtl number within the magnetically-constrained, high-$Q$ regimes., Comment: 24 pages, 11 figures

We present results from an asymptotic magnetohydrodynamic model that is suited for studying the rapidly rotating, low viscosity regime typical of the electrically conducting fluid interiors of planets and stars. We show that the presence of sufficiently strong magnetic fields prevents the formation of large-scale vortices and saturates the inverse cascade at a finite length-scale. This saturation corresponds to an equilibrated state in which the energetics of the depth-averaged flows are characterized by a balance of convective power input and ohmic dissipation. A quantitative criteria delineating the transition between finite-size flows and domain-filling (large-scale) vortices in electrically conducting fluids is found. By making use of the inferred and observed properties of planetary interiors, our results suggest that convection-driven large-scale vortices do not form in the electrically conducting regions of many bodies., Comment: Supplementary material included as ancillary file

The present numerical study aims at shedding light on the mechanism underlying the precessional instability in a sphere. Precessional instabilities in the form of parametric resonance due to topographic coupling have been reported in a spheroidal geometry both analytically and numerically. We show that such parametric resonances can also develop in spherical geometry due to the conical shear layers driven by the Ekman pumping singularities at the critical latitudes. Scaling considerations lead to a stability criterion of the form, $|P_o|>O(E^{4/5})$, where $P_o$ represents the Poincar\'e number and $E$ the Ekman number. The predicted threshold is consistent with our numerical simulations as well as previous experimental results. When the precessional forcing is supercriticial, our simulations show evidence of an inverse cascade, i.e. small scale flows merging into large scale cyclones with a retrograde drift. Finally, it is shown that this instability mechanism may be relevant to precessing celestial bodies such as the Earth and Earth's moon., Comment: published on PoF 2015

The dependence of the heat transfer, as measured by the nondimensional Nusselt number $Nu$, on Ekman pumping for rapidly rotating Rayleigh-B\'enard convection in an infinite plane layer is examined for fluids with Prandtl number $Pr = 1$. A joint effort utilizing simulations from the Composite Non-hydrostatic Quasi-Geostrophic model (CNH-QGM) and direct numerical simulations (DNS) of the incompressible fluid equations has mapped a wide range of the Rayleigh number $Ra$ - Ekman number $E$ parameter space within the geostrophic regime of rotating convection. Corroboration of the $Nu$-$Ra$ relation at $E = 10^{-7}$ from both methods along with higher $E$ covered by DNS and lower $E$ by the asymptotic model allows for this range of the heat transfer results. For stress-free boundaries, the relation $Nu-1 \propto (Ra E^{4/3} )^{\alpha} $ has the dissipation-free scaling of $\alpha = 3/2$ for all $E \leq 10^{-7}$. This is directly related to a geostrophic turbulent interior that throttles the heat transport supplied to the thermal boundary layers. For no-slip boundaries, the existence of ageostrophic viscous boundary layers and their associated Ekman pumping yields a more complex 2D surface in $Nu(E,Ra)$ parameter space. For $E<10^{-7}$ results suggest that the surface can be expressed as $Nu-1 \propto (1+ P(E)) (Ra E^{4/3} )^{3/2}$ indicating the dissipation-free scaling law is enhanced by Ekman pumping by the multiplicative prefactor $(1+ P(E))$ where $P(E) \approx 5.97 E^{1/8}$. It follows for $E<10^{-7}$ that the geostrophic turbulent interior remains the flux bottleneck in rapidly rotating Rayleigh-B\'enard convection. For $E\sim10^{-7}$, where DNS and asymptotic simulations agree quantitatively, it is found that the effects of Ekman pumping are sufficiently strong to influence the heat transport with diminished exponent $\alpha \approx 1.2$ and $Nu-1 \propto (Ra E^{4/3} )^{1.2}$., Comment: 9 pages, 14 figures

Numerical simulations of 3D, rapidly rotating Rayleigh-Benard convection are performed using an asymptotic quasi-geostrophic model that incorporates the effects of no-slip boundaries through (i) parameterized Ekman pumping boundary conditions, and (ii) a thermal wind boundary layer that regularizes the enhanced thermal fluctuations induced by pumping. The fidelity of the model, obtained by an asymptotic reduction of the Navier-Stokes equations that implicitly enforces a pointwise geostrophic balance, is explored for the first time by comparisons of simulations against the findings of direct numerical simulations and laboratory experiments. Results from these methods have established Ekman pumping as the mechanism responsible for significantly enhancing the vertical heat transport. This asymptotic model demonstrates excellent agreement over a range of thermal forcing for Pr ~1 when compared with results from experiments and DNS at maximal values of their attainable rotation rates, as measured by the Ekman number (E ~ 10^{-7}); good qualitative agreement is achieved for Pr > 1. Similar to studies with stress-free boundaries, four spatially distinct flow morphologies exists. Despite the presence of frictional drag at the upper and/or lower boundaries, a strong non-local inverse cascade of barotropic (i.e., depth-independent) kinetic energy persists in the final regime of geostrophic turbulence and is dominant at large scales. For mixed no-slip/stress-free and no-slip/no-slip boundaries, Ekman friction is found to attenuate the efficiency of the upscale energy transport and, unlike the case of stress-free boundaries, rapidly saturates the barotropic kinetic energy. For no-slip/no-slip boundaries, Ekman friction is strong enough to prevent the development of a coherent dipole vortex condensate. Instead vortex pairs are found to be intermittent, varying in both time and strength., Comment: 20 pages, 10 figures

Precession has been proposed as an alternative power source for planetary dynamos. Previous hydrodynamic simulations suggested that precession can generate very complex flows in planetary liquid cores [Y. Lin, P. Marti, and J. Noir, "Shear-driven parametric instability in a precessing sphere," Physics of Fluids 27, 046601 (2015)]. In the present study, we numerically investigate the magnetohydrodynamics of a precessing sphere. We demonstrate precession driven dynamos in different flow regimes, from laminar to turbulent flows. In particular, we highlight the magnetic field generation by large scale cyclonic vortices, which has not been explored previously. In this regime, dynamos can be sustained at relatively low Ekman numbers and magnetic Prandtl numbers, which paves the way for planetary applications., Comment: 22 pages, 12 figures

It is a well established result of linear theory that the influence of differing mechanical boundary conditions, i.e., stress-free or no-slip, on the primary instability in rotating convection becomes asymptotically small in the limit of rapid rotation. This is accounted for by the diminishing impact of the viscous stresses exerted within Ekman boundary layers and the associated vertical momentum transport by Ekman pumping. By contrast, in the nonlinear regime recent experiments and supporting simulations are now providing evidence that the efficiency of heat transport remains strongly influenced by Ekman pumping in the rapidly rotating limit. In this paper, a reduced model is developed for the case of low Rossby number convection in a plane layer geometry with no-slip upper and lower boundaries held at fixed temperatures. A complete description of the dynamics requires the existence of three distinct regions within the fluid layer: a geostrophically balanced interior where fluid motions are predominately aligned with the axis of rotation, Ekman boundary layers immediately adjacent to the bounding plates, and thermal wind layers driven by Ekman pumping in between. The reduced model uses a classical Ekman pumping parameterization to alleviate the need for spatially resolving the Ekman boundary layers. Results are presented for both linear stability theory and a special class of nonlinear solutions described by a single horizontal spatial wavenumber. It is shown that Ekman pumping allows for significant enhancement in the heat transport relative to that observed in simulations with stress-free boundaries. Without the intermediate thermal wind layer the nonlinear feedback from Ekman pumping would be able to generate a heat transport that diverges to infinity. This layer arrests this blowup resulting in finite heat transport at a significantly enhanced value., Comment: 38 pages, 14 figures

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., Comment: 13 pages, 12 figures, submitted to Physical Review E

The influence of fixed temperature and fixed heat flux thermal boundary conditions on rapidly rotating convection in the plane layer geometry is investigated for the case of stress-free mechanical boundary conditions. It is shown that whereas the leading order system satisfies fixed temperature boundary conditions implicitly, a double boundary layer structure is necessary to satisfy the fixed heat flux thermal boundary conditions. The boundary layers consist of a classical Ekman layer adjacent to the solid boundaries that adjust viscous stresses to zero, and a layer in thermal wind balance just outside the Ekman layers adjusts the temperature such that the fixed heat flux thermal boundary conditions are satisfied. The influence of these boundary layers on the interior geostrophically balanced convection is shown to be asymptotically weak, however. Upon defining a simple rescaling of the thermal variables, the leading order reduced system of governing equations are therefore equivalent for both boundary conditions. These results imply that any horizontal thermal variation along the boundaries that varies on the scale of the convection has no leading order influence on the interior convection.

The linear theory for rotating compressible convection in a plane layer geometry is presented for the astrophysically-relevant case of low Prandtl number gases. When the rotation rate of the system is large, the flow remains geostrophically balanced for all stratification levels investigated and the classical (i.e., incompressible) asymptotic scaling laws for the critical parameters are recovered. For sufficiently small Prandtl numbers, increasing stratification tends to further destabilise the fluid layer, decrease the critical wavenumber and increase the oscillation frequency of the convective instability. In combination, these effects increase the relative magnitude of the time derivative of the density perturbation contained in the conservation of mass equation to non-negligible levels; the resulting convective instabilities occur in the form of compressional quasi-geostrophic oscillations. We find that the anelastic equations, which neglect this term, cannot capture these instabilities and possess spuriously-growing eigenmodes in the rapidly rotating, low Prandtl number regime. It is shown that the Mach number for rapidly rotating compressible convection is intrinsically small for all background states, regardless of the departure from adiabaticity., Comment: 16 pages, 6 figures, 1 table. Changes made: added new Figure 1, added discussion of Mach number, added references, added Table of linear stability data