Differential rotation in giant planets maintained by density-stratified turbulent convection

The zonal winds on the surfaces of giant planets vary with latitude. Jupiter and Saturn, for example, have several bands of alternating eastward (prograde) and westward (retrograde) jets relative to the angular velocity of their global magnetic fields. These surface wind profiles are likely manifestations of the variations in depth and latitude of angular velocity deep within the liquid interiors of these planets. Two decades ago it was proposed that this differential rotation could be maintained by vortex stretching of convective fluid columns that span the interiors of these planets from the northern hemisphere surface to the southern hemisphere surface. This now classic mechanism explains the differential rotation seen in laboratory experiments and in computer simulations of, at best, weakly turbulent convection in rotating constant-density fluid spheres. However, these experiments and simulations are poor approximations for the density-stratified strongly-turbulent interiors of giant planets. The long thin global convective columns predicted by the classic geostrophic theory for these planets would likely not develop. Here we propose a much more robust mechanism for maintaining differential rotation in radius based on the local generation of vorticity as rising plumes expand and sinking plumes contract. Our high-resolution two-dimensional computer simulations demonstrate how this mechanism could maintain either prograde or retrograde surface winds in the equatorial region of a giant planet depending on how the density scale height varies with depth.Comment: Geophysical and Astrophysical Fluid Dynamics, in pres

Performance benchmarks for a next generation numerical dynamo model

Numerical simulations of the geodynamo have successfully represented many observable characteristics of the geomagnetic field, yielding insight into the fundamental processes that generate magnetic fields in the Earth's core. Because of limited spatial resolution, however, the diffusivities in numerical dynamo models are much larger than those in the Earth's core, and consequently, questions remain about how realistic these models are. The typical strategy used to address this issue has been to continue to increase the resolution of these quasi-laminar models with increasing computational resources, thus pushing them toward more realistic parameter regimes. We assess which methods are most promising for the next generation of supercomputers, which will offer access to O(106) processor cores for large problems. Here we report performance and accuracy benchmarks from 15 dynamo codes that employ a range of numerical and parallelization methods. Computational performance is assessed on the basis of weak and strong scaling behavior up to 16,384 processor cores. Extrapolations of our weak-scaling results indicate that dynamo codes that employ two-dimensional or three-dimensional domain decompositions can perform efficiently on up to ∼106 processor cores, paving the way for more realistic simulations in the next model generation

Theory of disk accretion onto supermassive black holes

Accretion onto supermassive black holes produces both the dramatic phenomena associated with active galactic nuclei and the underwhelming displays seen in the Galactic Center and most other nearby galaxies. I review selected aspects of the current theoretical understanding of black hole accretion, emphasizing the role of magnetohydrodynamic turbulence and gravitational instabilities in driving the actual accretion and the importance of the efficacy of cooling in determining the structure and observational appearance of the accretion flow. Ongoing investigations into the dynamics of the plunging region, the origin of variability in the accretion process, and the evolution of warped, twisted, or eccentric disks are summarized.Comment: Mostly introductory review, to appear in "Supermassive black holes in the distant Universe", ed. A.J. Barger, Kluwer Academic Publishers, in pres

Investigation of regional variation in core flow models using spherical Slepian functions

Abstract By assuming that changes in the magnetic field in the Earth’s outer core are advection-dominated on short timescales, models of the core surface flow can be deduced from secular variation. Such models are known to be under-determined and thus require other assumptions to produce feasible flows. There are regions where poor knowledge of the core flow dynamics gives rise to further uncertainty, such as within the tangent cylinder, and assumptions about the nature of the flow may lead to ambiguous patches, such as if it is assumed to be strongly tangentially geostrophic. We use spherical Slepian functions to spatially and spectrally separate core flow models, confining the flow to either inside or outside these regions of interest. In each region we examine the properties of the flow and analyze its contribution to the overall model. We use three forms of flow model: (a) synthetic models from randomly generated coefficients with blue, red and white energy spectra, (b) a snapshot of a numerical geodynamo simulation and (c) a model inverted from satellite magnetic field measurements. We find that the Slepian decomposition generates unwanted spatial leakage which partially obscures flow in the region of interest, particularly along the boundaries. Possible reasons for this include the use of spherical Slepian functions to decompose a scalar quantity that is then differentiated to give the vector function of interest, and the spectral frequency content of the models. These results will guide subsequent investigation of flow within localized regions, including applying vector Slepian decomposition methods

A comparison of no-slip, stress-free and inviscid models of rapidly rotating fluid in a spherical shell

We investigate how the choice of either no-slip or stress-free boundary conditions affects numerical models of rapidly rotating flow in Earth's core by computing solutions of the weakly-viscous magnetostrophic equations within a spherical shell, driven by a prescribed body force. For non-axisymmetric solutions, we show that models with either choice of boundary condition have thin boundary layers of depth E^(1/2), where E is the Ekman number, and a free-stream flow that converges to the formally inviscid solution. At Earth-like values of viscosity, the boundary layer thickness is approximately 1m, for either choice of condition. In contrast, the axisymmetric flows depend crucially on the choice of boundary condition, in both their structure and magnitude (either E^(-1/2) or E^(-1)). These very large zonal flows arise from requiring viscosity to balance residual axisymmetric torques. We demonstrate that switching the mechanical boundary conditions can cause a distinct change of structure of the flow, including a sign-change close to the equator, even at asymptotically low viscosity. Thus implementation of stress-free boundary conditions, compared with no-slip conditions, may yield qualitatively different dynamics in weakly-viscous magnetostrophic models of Earth's core. We further show that convergence of the free-stream flow to its asymptotic structure requires E ≤10^(-5)