skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Asymptotic behaviour of rotating convection-driven dynamos in the plane layer geometry
Dynamos driven by rotating convection in the plane layer geometry are investigated numerically for a range of Ekman number ( $$E$$ ), magnetic Prandtl number ( $Pm$ ) and Rayleigh number ( $Ra$ ). The primary purpose of the investigation is to compare results of the simulations with previously developed asymptotic theory that is applicable in the limit of rapid rotation. We find that all of the simulations are in the quasi-geostrophic regime in which the Coriolis and pressure gradient forces are approximately balanced at leading order, whereas all other forces, including the Lorentz force, act as perturbations. Agreement between simulation output and asymptotic scalings for the energetics, flow speeds, magnetic field amplitude and length scales is found. The transition from large-scale dynamos to small-scale dynamos is well described by the magnetic Reynolds number based on the small convective length scale, $$\widetilde {Rm}$$ , with large-scale dynamos preferred when $$\widetilde {Rm} \lesssim O(1)$$ . The magnitude of the large-scale magnetic field is observed to saturate and become approximately constant with increasing Rayleigh number. Energy spectra show that all length scales present in the flow field and the small-scale magnetic field are consistent with a scaling of $$E^{1/3}$$ , even in the turbulent regime. For a fixed value of $$E$$ , we find that the viscous dissipation length scale is approximately constant over a broad range of $Ra$ ; the ohmic dissipation length scale is approximately constant within the large-scale dynamo regime, but transitions to a $$\widetilde {Rm}^{-1/2}$$ scaling in the small-scale dynamo regime.  more » « less
Award ID(s):
1743852 1945270
PAR ID:
10408628
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Journal of Fluid Mechanics
Volume:
951
ISSN:
0022-1120
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Steady two-dimensional Rayleigh–Bénard convection between stress-free isothermal boundaries is studied via numerical computations. We explore properties of steady convective rolls with aspect ratios $${\rm \pi} /5\leqslant \varGamma \leqslant 4{\rm \pi}$$ , where $$\varGamma$$ is the width-to-height ratio for a pair of counter-rotating rolls, over eight orders of magnitude in the Rayleigh number, $$10^3\leqslant Ra\leqslant 10^{11}$$ , and four orders of magnitude in the Prandtl number, $$10^{-2}\leqslant Pr\leqslant 10^2$$ . At large $Ra$ where steady rolls are dynamically unstable, the computed rolls display $$Ra \rightarrow \infty$$ asymptotic scaling. In this regime, the Nusselt number $Nu$ that measures heat transport scales as $$Ra^{1/3}$$ uniformly in $Pr$ . The prefactor of this scaling depends on $$\varGamma$$ and is largest at $$\varGamma \approx 1.9$$ . The Reynolds number $Re$ for large- $Ra$ rolls scales as $$Pr^{-1} Ra^{2/3}$$ with a prefactor that is largest at $$\varGamma \approx 4.5$$ . All of these large- $Ra$ features agree quantitatively with the semi-analytical asymptotic solutions constructed by Chini & Cox ( Phys. Fluids , vol. 21, 2009, 083603). Convergence of $Nu$ and $Re$ to their asymptotic scalings occurs more slowly when $Pr$ is larger and when $$\varGamma$$ is smaller. 
    more » « less
  2. A numerical investigation of an asymptotically reduced model for quasigeostrophic Rayleigh-Bénard convection is conducted in which the depth-averaged flows are numerically suppressed by modifying the governing equations. At the largest accessible values of the Rayleigh number Ra, the Reynolds number and Nusselt number show evidence of approaching the diffusion-free scalings of Re ∼ RaE/Pr and Nu ∼ Pr−1/2Ra3/2E2, respectively, where E is the Ekman number and Pr is the Prandtl number. For large Ra, the presence of depth-invariant flows, such as large-scale vortices, yield heat and momentum transport scalings that exceed those of the diffusion-free scaling laws. The Taylor microscale does not vary significantly with increasing Ra, whereas the integral length scale grows weakly. The computed length scales remain O(1) with respect to the linearly unstable critical wave number; we therefore conclude that these scales remain viscously controlled. We do not find a point-wise Coriolis-inertia-Archimedean (CIA) force balance in the turbulent regime; interior dynamics are instead dominated by horizontal advection (inertia), vortex stretching (Coriolis) and the vertical pressure gradient. A secondary, subdominant balance between the Archimedean buoyancy force and the viscous force occurs in the interior and the ratio of the root mean square (rms) of these two forces is found to approach unity with increasing Ra. This secondary balance is attributed to the turbulent fluid interior acting as the dominant control on the heat transport. These findings indicate that a pointwise CIA balance does not occur in the high Rayleigh number regime of quasigeostrophic convection in the plane layer geometry. Instead, simulations are characterized by what may be termed a nonlocal CIA balance in which the buoyancy force is dominant within the thermal boundary layers and is spatially separated from the interior Coriolis and inertial forces. 
    more » « less
  3. ABSTRACT Numerical simulations are used to investigate large-scale (mean) magnetic field generation in rotating spherical dynamos. Beyond a certain threshold, we find that the magnitude of the mean magnetic field becomes nearly independent of the system rotation rate and buoyancy forcing. The analysis suggests that this saturation arises from the Malkus-Proctor mechanism in which a Coriolis-Lorentz force balance is achieved in the zonal component of the mean momentum equation. When based on the large-scale magnetic field, the Elsasser number is near unity in the saturated regime. The results show that the large and small magnetic field saturate via distinct mechanisms in rapidly rotating dynamos, and that only the axisymmetric component of the magnetic field appears to follow an Elsasser number scaling. 
    more » « less
  4. ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This leads to a larger characteristic size of the, now advected, magnetic field when compared to the local disc height. We characterize the complete dynamo mechanism with two time-scales: one for the local magnetic field generation, $$t_{\rm gen}$$, and one for the large-scale scale advection, $$t_{\rm adv}$$. Whereas the dynamo we describe is non-linear, we explore the potential of linear mean field models to replicate its core features. Our findings indicate that traditional $$\alpha$$-dynamo models, often computed in stratified shearing box simulations, are inadequate and that the effective large-scale dynamics is better described by the shear current effects or stochastic $$\alpha$$-dynamos. 
    more » « less
  5. Summary The large-scale dynamics of convection-driven dynamos in a spherical shell, as relevant to the geodynamo, is analyzed with numerical simulation data and asymptotic theory. An attempt is made to determine the asymptotic size (with the small parameter being the Ekman number, Ek) of the forces, and the associated velocity and magnetic fields. In agreement with previous work, the leading order mean force balance is shown to be thermal wind (Coriolis, pressure gradient, buoyancy) in the meridional plane and Coriolis-Lorentz in the zonal direction. The Lorentz force is observed to be weaker than the mean buoyancy force across a range of Ek and thermal forcing; the relative difference in these forces appears to be O(Ek1/6) within the parameter space investigated. We find that the thermal wind balance requires that the mean zonal velocity scales as O(Ek−1/3), whereas the meridional circulation is asymptotically smaller by a factor of O(Ek1/6). The mean temperature equation shows a balance between thermal diffusion and the divergence of the convective heat flux, indicating the presence of a mean temperature length scale of size O(Ek1/6). Neither the mean nor the fluctuating magnetic field show a strong dependence on the Ekman number, though the simulation data shows evidence of a mean magnetic field length scale of size O(Ek1/6). A consequence of the asymptotic ordering of the forces is that Taylor’s constraint is satisfied to accuracy O(Ek1/6), despite the absence of a leading-order magnetostrophic balance. Further consequences of the force balance are discussed with respect to the large-scale flows thought to be important for the geodynamo. 
    more » « less