Fluid motion in the liquid iron core of the Earth generates a magnetic field, which protects the surface environment from the harmful effects of solar wind and cosmic rays. The low viscosity of liquid iron (comparable to that of liquid water) ensures that the fluid motion is highly turbulent and time dependent. Accurate numerical simulations of the fluid motion and field generation are beyond the limit of current computational power, so modelers are forced to use values for the fluid viscosity and other physical parameters that are very far from Earth-like values. As a consequence, the nature of the dynamics is altered and our ability to address important scientific questions is limited. Only modest improvements in the models are expected through brute-force improvements in resolution, so there is a pressing need to develop better models for dealing with turbulence in the core.
We propose to develop a second generation of geodynamo model that advances the current models in two ways. First, the new model will employ a method of solution that is better suited to parallel computing. Second, the new model will include a more sophisticated treatment of small-scale turbulence. To achieve these goals we will use the finite-element method to solve the governing equations and implement an adaptive, scale-similarity model to account for the effects of turbulence. The finite-element method reduces the amount of communication on a parallel computer and allows us to locally refine the resolution of the calculation, but only where it is needed. The scale-similarity method is a proven technique that has been successfully implemented in a plane-layer dynamo model to account for complex interactions between the flow and the magnetic field at small scales. We conservatively estimate that the introduction of the scale-similarity model in the geodynamo model will allow us to lower the fluid viscosity by two orders of magnitude, relative to current capabilities. This may be sufficient to reach the dynamical regime (e.g. a Taylor state) where the Earth's geodynamo is thought to operate. The new model will also allow us to probe the geophysically relevant behavior of dynamos at low magnetic Prandtl number Pm, which represents the ratio of fluid viscosity to magnetic diffusivity. Current state-of-the-art numerical simulations assume Pm > 1, whereas the expected value for the Earth's core is Pm ~ 1.0e-6 . The operation of dynamos at low Pm raises a number of important scientific questions and poses a substantial computational challenge.
The Earth's magnetic field protects the surface from the harmful effects of radiation, which contributes to the habitability of our planet. The origin of the magnetic field is thought to be a result of fluid motion in the liquid metal core of the Earth, although the process is not well understood. Numerical simulations of the process offer valuable insights, but a comprehensive understanding is lacking because we can not adequately account for the influence of small-scale turbulence. It is not currently possible to fully resolved the turbulent flow in numerical simulations, so we are forced to develop models to represent the influence of turbulence in simulations. During this project we developed sophisticated turbulence models that account for the effects of planetary rotation and large-scale magnetic fields on the structure of small scale flow. We tested the turbulence models using state-of-the-art computational facilities and implemented these models in numerical geodynamo simulations. One of the unexpected outcomes was the recognition that small-scale turbulence can substantially enhance the generation of magnetic field. In many cases the scale of the turbulent flow was too small to directly generate magnetic field. However, the small-scale flow was able to transfer energy in to larger scale flows, which are much more effective at generating field. In addition, the transfer of energy from small-scale to large-scales occurred precisely in the region where field generation is most intense. We have begun to use our improved geodynamo simulations to gains new insights into the origin of the Earth's magnetic fields. The methods developed in this study should be broadly applicable to other problems, including models of oceanic circulations and astrophysical problems, where the combination of rotation and/or magnetic fields have an important influence on the dynamics.