Dynamo theory

 ( Geomagnetism ) 

In physics, the dynamo theory proposes a mechanism by which a celestial body such as Earth or a star generates a magnetic field. The dynamo theory describes the process through which a rotating, convection, and conducting fluid can maintain a magnetic field over astronomical time scales. A dynamo is thought to be the source of the Earth's magnetic field and the magnetic fields of Mercury and the Jovian planets.

---

History of theory When William Gilbert published De Magnete in 1600, he concluded that the Earth is magnetic and proposed the first hypothesis for the origin of this magnetism: permanent magnetism such as that found in lodestone. In 1822, André-Marie Ampère proposed that internal currents are responsible for Earth's magnetism. In 1919, Joseph Larmor proposed that a dynamo might be generating the field. Reprinted in Engineering, vol. 108, pages 461ff (3 October 1919). However, even after he advanced his hypothesis, some prominent scientists advanced alternative explanations. The Nobel Prize winner Patrick Blackett did a series of experiments looking for a fundamental relation between angular momentum and magnetic moment, but found none. claim that in 1905, shortly after composing his special relativity paper, Albert Einstein described the origin of the Earth's magnetic field as being one of the great unsolved problems facing modern . However, they do not provide details on where he made this statement.

Walter M. Elsasser, considered a "father" of the presently accepted dynamo theory as an explanation of the Earth's magnetism, proposed that this magnetic field resulted from electric currents induced in the fluid outer core of the Earth. He revealed the history of the Earth's magnetic field through pioneering the study of the magnetic orientation of minerals in rocks.

In order to maintain the magnetic field against decay (which would occur for the dipole field in 20,000 years), the outer core must be convecting. The convection is likely some combination of thermal and compositional convection. The mantle controls the rate at which heat is extracted from the core. Heat sources include gravitational energy released by the compression of the core, gravitational energy released by the rejection of light elements (probably sulfur, oxygen, or silicon) at the inner core boundary as it grows, latent heat of crystallization at the inner core boundary, and radioactivity of potassium, uranium and thorium.

At the dawn of the 21st century, numerical modeling of the Earth's magnetic field is far from precise. Initial models are focused on field generation by convection in the planet's fluid outer core. It was possible to show the generation of a strong, Earth-like field when the model assumed a uniform core-surface temperature and exceptionally high viscosities for the core fluid. Computations which incorporated more realistic parameter values yielded magnetic fields that were less Earth-like, but indicated that model refinements may ultimately lead to an accurate analytic model. Slight variations in the core-surface temperature, in the range of a few millikelvins, result in significant increases in convective flow and produce more realistic magnetic fields.

---

Formal definition Dynamo theory describes the process through which a rotating, convecting, and electrically conducting fluid acts to maintain a magnetic field. This theory is used to explain the presence of anomalously long-lived magnetic fields in astrophysical bodies. The conductive fluid in the geodynamo is liquid iron in the outer core, and in the solar dynamo is ionized gas at the tachocline. Dynamo theory of astrophysical bodies uses magnetohydrodynamic equations to investigate how the fluid can continuously regenerate the magnetic field.

It was once believed that the dipole, which comprises much of the Earth's magnetic field and is misaligned along the rotation axis by 11.3 degrees, was caused by permanent magnetization of the materials in the earth. This means that dynamo theory was originally used to explain the Sun's magnetic field in its relationship with that of the Earth. However, this hypothesis, which was initially proposed by Joseph Larmor in 1919, has been modified due to extensive studies of magnetic secular variation, paleomagnetism (including polarity reversals), seismology, and the solar system's abundance of elements. Also, the application of the theories of Carl Friedrich Gauss to magnetic observations showed that Earth's magnetic field had an internal, rather than external, origin.

There are three requisites for a dynamo to operate:

• An electrically conductive fluid medium
• Kinetic energy provided by planetary rotation
• An internal energy source to drive convective motions within the fluid.
  E. Pallé (2025). 9781441916839, Springer. . ISBN 9781441916839

In the case of the Earth, the magnetic field is induced and constantly maintained by the convection of liquid iron in the outer core. A requirement for the induction of field is a rotating fluid. Rotation in the outer core is supplied by the Coriolis effect caused by the rotation of the Earth. The Coriolis force tends to organize fluid motions and electric currents into columns (also see ) aligned with the rotation axis. Induction or generation of magnetic field is described by the induction equation: $$\backslash frac\{\backslash partial\; \backslash mathbf\{B\}\}\{\backslash partial\; t\}\; =\; \backslash eta\; \backslash nabla^2\; \backslash mathbf\{B\}\; +\; \backslash nabla\; \backslash times\; (\backslash mathbf\{u\}\; \backslash times\; \backslash mathbf\{B\})$$ where is velocity, is magnetic field, is time, and $\backslash eta\; =\; 1\; /\; (\backslash sigma\backslash mu)$ is the magnetic diffusivity with $\backslash sigma$ electrical conductivity and $\backslash mu$ permeability. The ratio of the second term on the right hand side to the first term gives the magnetic Reynolds number, a dimensionless ratio of advection of magnetic field to diffusion.

---

Tidal heating supporting a dynamo Tidal forces between celestial orbiting bodies cause friction that heats up their interiors. This is known as tidal heating, and it helps keep the interior in a liquid state. A liquid interior that can conduct electricity is required to produce a dynamo. Saturn's Enceladus and Jupiter's Io have enough tidal heating to liquify their inner cores, but they may not create a dynamo because they cannot conduct electricity. Mercury, despite its small size, has a magnetic field, because it has a conductive liquid core created by its iron composition and friction resulting from its highly elliptical orbit. It is theorized that the Moon once had a magnetic field, based on evidence from magnetized lunar rocks, due to its short-lived closer distance to Earth creating tidal heating. An orbit and rotation of a planet helps provide a liquid core, and supplements kinetic energy that supports a dynamo action.

---

Kinematic dynamo theory In kinematic dynamo theory the velocity field is prescribed, instead of being a dynamic variable: The model makes no provision for the flow distorting in response to the magnetic field. This method cannot provide the time variable behaviour of a fully nonlinear chaotic dynamo, but can be used to study how magnetic field strength varies with the flow structure and speed.

Using Maxwell's equations simultaneously with the curl of Ohm's law, one can derive what is basically a linear eigenvalue equation for magnetic fields ( ), which can be done when assuming that the magnetic field is independent from the velocity field. One arrives at a critical magnetic Reynolds number, above which the flow strength is sufficient to amplify the imposed magnetic field, and below which the magnetic field dissipates.

---

Practical measure of possible dynamos The most functional feature of kinematic dynamo theory is that it can be used to test whether a velocity field is or is not capable of dynamo action. By experimentally applying a certain velocity field to a small magnetic field, one can observe whether the magnetic field tends to grow (or not) in response to the applied flow. If the magnetic field does grow, then the system is either capable of dynamo action or is a dynamo, but if the magnetic field does not grow, then it is simply referred to as “not a dynamo”.

An analogous method called the membrane paradigm is a way of looking at that allows for the material near their surfaces to be expressed in the language of dynamo theory.

---

Spontaneous breakdown of a topological supersymmetry Kinematic dynamo can be also viewed as the phenomenon of the spontaneous breakdown of the topological supersymmetry of the associated stochastic differential equation related to the flow of the background matter. Within stochastic supersymmetric theory, this supersymmetry is an intrinsic property of all stochastic differential equations, its interpretation is that the model's phase space preserves continuity via continuous time flows. When the continuity of that flow spontaneously breaks down, the system is in the stochastic state of Chaos theory. In other words, kinematic dynamo arises because of chaotic flow in the underlying background matter.

---

Nonlinear dynamo theory The kinematic approximation becomes invalid when the magnetic field becomes strong enough to affect the fluid motions. In that case the velocity field becomes affected by the Lorentz force, and so the induction equation is no longer linear in the magnetic field. In most cases this leads to a quenching of the amplitude of the dynamo. Such dynamos are sometimes also referred to as hydromagnetic dynamos. Virtually all dynamos in astrophysics and geophysics are hydromagnetic dynamos.

The main idea of the theory is that any small magnetic field existing in the outer core creates currents in the moving fluid there due to Lorentz force. These currents create further magnetic field due to Ampere's law. With the fluid motion, the currents are carried in a way that the magnetic field gets stronger (as long as $\backslash ;\; \backslash mathbf\{u\}\; \backslash cdot\; (\; \backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}\; )\; \backslash ;$ is negative). Thus a "seed" magnetic field can get stronger and stronger until it reaches some value that is related to existing non-magnetic forces.

Numerical models are used to simulate fully nonlinear dynamos. The following equations are used:

• The induction equation, presented above.
• Maxwell's equations for negligible electric field: $$\backslash begin\{align\}$$

&\nabla \cdot \mathbf{B} = 0 \\[1ex]
&\nabla \times \mathbf{B} = \mu_0 \mathbf{J}
     

\end{align}

• The continuity equation for conservation of mass, for which the Boussinesq approximation is often used: $$\backslash nabla\; \backslash cdot\; \backslash mathbf\{u\}\; =\; 0,$$
• The Navier-Stokes equation for conservation of momentum, again in the same approximation, with the magnetic force and gravitation force as the external forces: $$\backslash frac\{D\backslash mathbf\{u\}\}\{Dt\}\; =\; -\backslash frac\{1\}\{\backslash rho\_0\}\backslash nabla\; p\; +\; \backslash nu\; \backslash nabla^2\; \backslash mathbf\{u\}\; +\; \backslash rho\text{'}\; \backslash mathbf\{g\}\; +\; 2\; \backslash boldsymbol\{\backslash Omega\}\; \backslash times\; \backslash mathbf\{u\}\; +\; \backslash boldsymbol\{\backslash Omega\}\; \backslash times\; \backslash boldsymbol\{\backslash Omega\}\; \backslash times\; \backslash mathbf\{R\}\; +\; \backslash frac\{1\}\{\backslash rho\_0\}\backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}\; ~,$$ where $\backslash ,\backslash nu\backslash ,$ is the kinematic viscosity, $\backslash ,\backslash rho\_0\backslash ,$is the mean density and $\backslash rho\text{'}$ is the relative density perturbation that provides buoyancy (for thermal convection $\backslash ;\backslash rho\text{'}\; =\; \backslash alpha\; \backslash Delta\; T\backslash ;$ where $\backslash ,\backslash alpha\backslash ,$ is coefficient of thermal expansion), $\backslash ,\backslash Omega\backslash ,$ is the rotation rate of the Earth, and $\backslash ,\backslash mathbf\{J\}\backslash ,$ is the electric current density.
• A transport equation, usually of heat (sometimes of light element concentration): $$\backslash frac\{\backslash ,\backslash partial\; T\backslash ,\}\{\backslash partial\; t\}\; =\; \backslash kappa\; \backslash nabla^2\; T\; +\; \backslash varepsilon$$ where is temperature, $\backslash ;\backslash kappa\; =\; k\; /\; \backslash rho\; c\_p\; \backslash ;$ is the thermal diffusivity with thermal conductivity, $\backslash ,c\_p\backslash ,$ heat capacity, and $\backslash rho$ density, and $\backslash varepsilon$ is an optional heat source. Often the pressure is the dynamic pressure, with the hydrostatic pressure and centripetal potential removed.

These equations are then non-dimensionalized, introducing the non-dimensional parameters, $$R\_\backslash mathsf\{a\}\; =\; \backslash frac\{\backslash ,\; g\; \backslash alpha\; T\; D^3\backslash ,\}\{\backslash nu\; \backslash kappa\}\; \backslash ;,\backslash quad\; E\; =\; \backslash frac\{\backslash nu\}\{\backslash ,\backslash Omega\; D^2\backslash ,\}\; \backslash ;,\backslash quad\; P\_\backslash mathsf\{r\}\; =\; \backslash frac\{\backslash ,\backslash nu\backslash ,\}\{\backslash kappa\}\; \backslash ;,\backslash quad\; P\_\backslash mathsf\{m\}\; =\; \backslash frac\{\backslash ,\backslash nu\backslash ,\}\{\backslash eta\}$$ where is the Rayleigh number, the Ekman number, and the Prandtl number and magnetic Prandtl number. Magnetic field scaling is often in Elsasser number units $B\; =\; (\backslash rho\; \backslash Omega/\backslash sigma)^\{1/2\}.$

---

Energy conversion between magnetic and kinematic energy The scalar product of the above form of Navier-Stokes equation with $\backslash ;\backslash rho\_0\; \backslash mathbf\{u\}\backslash ;$ gives the rate of increase of kinetic energy density, $\backslash ;\; \backslash tfrac\{1\}\{2\}\; \backslash rho\_0\; u^2\; c\; \backslash ;$, on the left-hand side. The last term on the right-hand side is then $\backslash ;\; \backslash mathbf\{u\}\; \backslash cdot\; (\; \backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}\; )\; \backslash ;$, the local contribution to the kinetic energy due to Lorentz force.

The scalar product of the induction equation with $\backslash tfrac\{1\}\{\backslash mu\_0\}\backslash mathbf\{B\}$ gives the rate of increase of the magnetic energy density, $\backslash ;\backslash tfrac\{1\}\{2\}\; \backslash mu\_0\; B^2\backslash ;$, on the left-hand side. The last term on the right-hand side is then $\backslash tfrac\{1\}\{\backslash mu\_0\}\; \backslash mathbf\{B\}\; \backslash cdot\; \backslash left(\; \backslash nabla\; \backslash times\; \backslash left(\; \backslash mathbf\{u\}\; \backslash times\; \backslash mathbf\{B\}\; \backslash right)\; \backslash right)\; \backslash ;.$ Since the equation is volume-integrated, this term is equivalent up to a boundary term (and with the double use of the scalar triple product identity) to $\backslash ;\; -\backslash mathbf\{u\}\; \backslash cdot\; \backslash left(\; \backslash tfrac\{1\}\{\backslash mu\_0\}\; \backslash left(\; \backslash nabla\; \backslash times\; \backslash mathbf\{B\}\; \backslash right)\; \backslash times\; \backslash mathbf\{B\}\; \backslash right)\; =\; -\backslash mathbf\{u\}\; \backslash cdot\; \backslash left(\; \backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}\; \backslash right)\; ~$ (where one of Maxwell's equations was used). This is the local contribution to the magnetic energy due to fluid motion.

Thus the term $\backslash ;-\backslash mathbf\{u\}\; \backslash cdot\; (\; \backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}\; )\; \backslash ;$ is the rate of transformation of kinetic energy to magnetic energy. This has to be non-negative at least in part of the volume, for the dynamo to produce magnetic field.

From the diagram above, it is not clear why this term should be positive. A simple argument can be based on consideration of net effects. To create the magnetic field, the net electric current must wrap around the axis of rotation of the planet. In that case, for the term to be positive, the net flow of conducting matter must be towards the axis of rotation. The diagram only shows a net flow from the poles to the equator. However mass conservation requires an additional flow from the equator toward the poles. If that flow was along the axis of rotation, that implies the circulation would be completed by a flow from the ones shown towards the axis of rotation, producing the desired effect.

---

Order of magnitude of the magnetic field created by Earth's dynamo The above formula for the rate of conversion of kinetic energy to magnetic energy, is equivalent to a rate of work done by a force of $\backslash ;\backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}\backslash ;$ on the outer core matter, whose velocity is $\backslash mathbf\{u\}$. This work is the result of non-magnetic forces acting on the fluid.

Of those, the gravitational force and the centrifugal force are conservative and therefore have no overall contribution to fluid moving in closed loops. Ekman number (defined above), which is the ratio between the two remaining forces, namely the viscosity and Coriolis force, is very low inside Earth's outer core, because its viscosity is low (1.2–1.5 ×10 pascal-second) due to its liquidity.

Thus the main time-averaged contribution to the work is from Coriolis force, whose size is $\backslash ;-2\backslash rho\backslash ,\backslash mathbf\{\backslash Omega\}\; \backslash times\; \backslash mathbf\{u\}\; \backslash ;,$ though this quantity and $\backslash mathbf\{J\}\; \backslash times\; \backslash mathbf\{B\}$ are related only indirectly and are not in general equal locally (thus they affect each other but not in the same place and time).

The current density is itself the result of the magnetic field according to Ohm's law. Again, due to matter motion and current flow, this is not necessarily the field at the same place and time. However these relations can still be used to deduce orders of magnitude of the quantities in question.

In terms of order of magnitude, $\backslash ;\; J\; \backslash ,\; B\; \backslash sim\; \backslash rho\backslash ,\; \backslash Omega\backslash ,\; u\; \backslash ;$ and $\backslash ;\; J\; \backslash sim\; \backslash sigma\; u\; B\backslash ;$, giving $\backslash ;\backslash sigma\backslash ,u\backslash ,\; B^2\; \backslash sim\; \backslash rho\backslash ,\; \backslash Omega\backslash ,u\; \backslash ;,$ or: $$B\; \backslash sim\; \backslash sqrt\{\backslash frac\{\backslash ,\backslash rho\backslash ,\backslash Omega\backslash ,\}\{\backslash sigma\}\backslash ;\}$$

The exact ratio between both sides is the square root of Elsasser number.

Note that the magnetic field direction cannot be inferred from this approximation (at least not its sign) as it appears squared, and is, indeed, sometimes reversed, though in general it lies on a similar axis to that of $\backslash mathbf\{\backslash Omega\}$.

For earth outer core, is approximately 10⁴ kg/m³,de Wijs, G. A., Kresse, G., Vočadlo, L., Dobson, D., Alfe, D., Gillan, M. J., & Price, G. D. (1998). target="_blank" rel="nofollow"> The viscosity of liquid iron at the physical conditions of the Earth's core. Nature, 392(6678), 805.   = 2/day = 7.3×10⁻⁵/second   and     is approximately   10⁷Ω⁻¹m⁻¹ .Ohta, K., Kuwayama, Y., Hirose, K., Shimizu, K., & Ohishi, Y. (2016). Experimental determination of the electrical resistivity of iron at Earth’s core conditions. Nature, 534(7605), 95. Link to a summary This gives   2.7×10⁻⁴ Tesla.

The magnetic field of a magnetic dipole has an inverse cubic dependence in distance, so its order of magnitude at the earth surface can be approximated by multiplying the above result with giving   2.5×10⁻⁵ Tesla, not far from the measured value of 3×10⁻⁵ Tesla at the equator.

---

Numerical models Broadly, models of the geodynamo attempt to produce magnetic fields consistent with observed data given certain conditions and equations as mentioned in the sections above. Implementing the magnetohydrodynamic equations successfully was of particular significance because they pushed dynamo models to self-consistency. Though geodynamo models are especially prevalent, dynamo models are not necessarily restricted to the geodynamo; solar and general dynamo models are also of interest. Studying dynamo models has utility in the field of geophysics as doing so can identify how various mechanisms form magnetic fields like those produced by astrophysical bodies like Earth and how they cause magnetic fields to exhibit certain features, such as pole reversals.

The equations used in numerical models of dynamo are highly complex. For decades, theorists were confined to two dimensional kinematic dynamo models described above, in which the fluid motion is chosen in advance and the effect on the magnetic field calculated. The progression from linear to nonlinear, three dimensional models of dynamo was largely hindered by the search for solutions to magnetohydrodynamic equations, which eliminate the need for many of the assumptions made in kinematic models and allow self-consistency. The first self-consistent dynamo models, ones that determine both the fluid motions and the magnetic field, were developed by two groups in 1995, one in Japan and one in the United States. The latter was made as a model with regards to the geodynamo and received significant attention because it successfully reproduced some of the characteristics of the Earth's field. Following this breakthrough, there was a large swell in development of reasonable, three dimensional dynamo models.

Though many self-consistent models now exist, there are significant differences among the models, both in the results they produce and the way they were developed. Given the complexity of developing a geodynamo model, there are many places where discrepancies can occur such as when making assumptions involving the mechanisms that provide energy for the dynamo, when choosing values for parameters used in equations, or when normalizing equations. In spite of the many differences that may occur, most models have shared features like clear axial dipoles. In many of these models, phenomena like secular variation and geomagnetic polarity reversals have also been successfully recreated.

---

Observations Many observations can be made from dynamo models. Models can be used to estimate how magnetic fields vary with time and can be compared to observed Paleomagnetism data to find similarities between the model and the Earth. Due to the uncertainty of paleomagnetic observations, however, comparisons may not be entirely valid or useful. Simplified geodynamo models have shown relationships between the dynamo number (determined by variance in rotational rates in the outer core and mirror-asymmetric convection (e.g. when convection favors one direction in the north and the other in the south)) and magnetic pole reversals as well as found similarities between the geodynamo and the Sun's dynamo. In many models, it appears that magnetic fields have somewhat random magnitudes that follow a normal trend that average to zero. In addition to these observations, general observations about the mechanisms powering the geodynamo can be made based on how accurately the model reflects actual data collected from Earth.

---

Modern modelling

The complexity of dynamo modelling is so great that models of the geodynamo are limited by the current power of , particularly because calculating the Ekman number and Rayleigh number number of the outer core is extremely difficult and requires a vast number of computations.

Many improvements have been proposed in dynamo modelling since the self-consistent breakthrough in 1995. One suggestion in studying the complex magnetic field changes is applying to simplify computations. Ultimately, until considerable improvements in computer power are made, the methods for computing realistic dynamo models will have to be made more efficient, so making improvements in methods for computing the model is of high importance for the advancement of numerical dynamo modelling.

---

Notable people

• Stanislav I. Braginsky, research geophysicist

---

See also

• Antidynamo theorem
• Rotating magnetic field
• Secular variation

• (1996). 9780124912465, Academic Press. ISBN 9780124912465

Categories: Geomagnetism, Plasma Theory And Modeling, Magnetohydrodynamics, Structure Of The Earth, Computational Physics, Unsolved Problems In Physics, Magnetism In Astronomy

Post Comment