Force balance in numerical geodynamo simulations: a systematic study

T. Schwaiger,T. Gastine,J. Aubert

DOI: https://doi.org/10.1093/gji/ggz192

2019-05-13

Abstract:Dynamo action in the Earth's outer core is expected to be controlled by a balance between pressure, Coriolis, buoyancy and Lorentz forces, with marginal contributions from inertia and viscous forces. Current numerical simulations of the geodynamo, however, operate at much larger inertia and viscosity because of computational limitations. This casts some doubt on the physical relevance of these models. Our work aims at finding dynamo models in a moderate computational regime which reproduce the leading-order force balance of the Earth. By performing a systematic parameter space survey with Ekman numbers in the range $10^{-6} \leq E \leq 10^{-4}$, we study the variations of the force balance when changing the forcing (Rayleigh number, $Ra$) and the ratio between viscous and magnetic diffusivities (magnetic Prandtl number, $Pm$). For dipole-dominated dynamos, we observe that the force balance is structurally robust throughout the investigated parameter space, exhibiting a quasi-geostrophic (QG) balance (balance between Coriolis and pressure forces) at zeroth order, followed by a first-order MAC balance between the ageostrophic Coriolis, buoyancy and Lorentz forces. At second order this balance is disturbed by contributions from inertia and viscous forces. Dynamos with a different sequence of the forces, where inertia and/or viscosity replace the Lorentz force in the first-order force balance, can only be found close to the onset of dynamo action and in the multipolar regime. Our study illustrates that most classical numerical dynamos are controlled by a QG-MAC balance, while cases where viscosity and inertia play a dominant role are the exception rather than the norm.

Geophysics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is whether the force - balance structure can reflect the physical conditions in the Earth's outer core in the geodynamo simulation of the Earth's outer core. Specifically, the paper focuses on how to find a geodynamo model that can reproduce the main force balance in the Earth's outer core (i.e., the balance among pressure, Coriolis force, buoyancy, and Lorentz force) under moderate computational conditions. Since the current numerically - simulated geodynamos usually operate under large inertia and viscosity, which is far from the actual conditions of the Earth's outer core, the physical relevance of these models is questioned. By systematically studying geodynamo models under different parameters, the paper explores the changes in force balance under different Rayleigh numbers (\( \text{Ra} \)) and magnetic Prandtl numbers (\( \text{Pm} \)), especially the interactions among the Coriolis force, buoyancy, and Lorentz force, and the impact of these changes on the geodynamo model.

The main objectives of the paper are:
1. **Verify the structure of force balance**: Confirm whether the geodynamo model can achieve a force - balance structure similar to that of the Earth's outer core under different parameter settings.
2. **Explore the parameter space**: By systematically changing parameters (such as Ekman number, Rayleigh number, and magnetic Prandtl number), study the variation law of force balance.
3. **Evaluate the physical relevance of the model**: By introducing a parameter to quantify the distance between first - order and second - order forces, evaluate the consistency between the model force balance and the expected force balance of the Earth's outer core.

Through these studies, the paper aims to provide a more accurate numerical model for understanding the geodynamo mechanism in the Earth's outer core and explore the applicability and physical significance of these models under different conditions.