Search

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.

We simulate numerically Boussinesq convection in non-rotating spherical shells for a fluid with a Prandtl number of unity and for Rayleigh numbers up to 109. In this geometry, curvature and radial variations of the gravitational acceleration yield asymmetric boundary layers. A systematic parameter study for various radius ratios (from η = ri/ro = 0.2 to η = 0.95) and gravity profiles allows us to explore the dependence of the asymmetry on these parameters. We find that the average plume spacing is comparable between the spherical inner and outer bounding surfaces. An estimate of the average plume separation allows us to accurately predict the boundary layer asymmetry for the various spherical shell configurations explored here. The mean temperature and horizontal velocity profiles are in good agreement with classical Prandtl-Blasius laminar boundary layer profiles, provided the boundary layers are analysed in a dynamical frame that fluctuates with the local and instantaneous boundary layer thicknesses. The scaling properties of the Nusselt and Reynolds numbers are investigated by separating the bulk and boundary layer contributions to the thermal and viscous dissipation rates using numerical models with η = 0.6 and with gravity proportional to 1/r2. We show that our spherical models are consistent with the predictions of Grossmann & Lohse's (J. Fluid Mech., vol. 407, 2000, pp. 27-56) theory and that Nu(Ra) and Re(Ra) scalings are in good agreement with plane layer results.

The surface zonal winds observed in the giant planets form a complex jet pattern with alternating prograde and retrograde direction. While the main equatorial band is prograde on the gas giants, both ice giants have a pronounced retrograde equatorial jet.We use three-dimensional numerical models of compressible convection in rotating spherical shells to explore the properties of zonal flows in different regimes where either rotation or buoyancy dominates the force balance. We conduct a systematic parameter study to quantify the dependence of zonal flows on the background density stratification and the driving of convection.In our numerical models, we find that the direction of the equatorial zonal wind is controlled by the ratio of the global-scale buoyancy force and the Coriolis force. The prograde equatorial band maintained by Reynolds stresses is found in the rotation-dominated regime. In cases where buoyancy dominates Coriolis force, the angular momentum per unit mass is homogenized and the equatorial band is retrograde, reminiscent to those observed in the ice giants. In this regime, the amplitude of the zonal jets depends on the background density contrast with strongly stratified models producing stronger jets than comparable weakly stratified cases. Furthermore, our results can help to explain the transition between solar-like (i.e. prograde at the equator) and the "anti-solar" differential rotations (i.e. retrograde at the equator) found in anelastic models of stellar convection zones.In the strongly stratified cases, we find that the leading order force balance can significantly vary with depth. While the flow in the deep interior is dominated by rotation, buoyancy can indeed become larger than Coriolis force in a thin region close to the surface. This so-called "transitional regime" has a visible signature in the main equatorial jet which shows a pronounced dimple where flow amplitudes notably decay towards the equator. A similar dimple is observed on Jupiter, which suggests that convection in the planet interior could possibly operate in this regime. © 2013 Elsevier Inc.