Dynamo theory

src: i.ytimg.com

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, convecting, and electrically 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, as well as the magnetic fields of other planets.

Video Dynamo theory

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 1919, Joseph Larmor proposed that a dynamo might be generating the field. However, even after he advanced his hypothesis, some prominent scientists advanced alternative explanations. Einstein believed that there might be an asymmetry between the charges of the electron and proton so that the Earth's magnetic field would be produced by the entire Earth. 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.

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 ohmic 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 has not been successfully demonstrated, but appears to be in reach. 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 also point the way to model refinements which 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.

Maps Dynamo theory

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.

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 Taylor columns) aligned with the rotation axis. Induction or creation of magnetic field is described by the induction equation:

${\displaystyle {\frac {\partial \mathbf {B} }{\partial t}}=\eta \nabla ^{2}\mathbf {B} +\nabla \times (\mathbf {u} \times \mathbf {B} )}$

where u is velocity, B is magnetic field, t is time, and ${\displaystyle \eta =1/(\sigma \mu )}$ is the magnetic diffusivity with ${\displaystyle \sigma }$ electrical conductivity and ${\displaystyle \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 liquid. 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.

src: cdn.vox-cdn.com

Kinematic dynamo theory

In kinematic dynamo theory the velocity field is prescribed, instead of being a dynamic variable. This method cannot provide the time variable behavior of a fully nonlinear chaotic dynamo but is useful in studying 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 the linear eigenvalue equation for magnetic fields (B) 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 it decays.

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 applying a certain velocity field to a small magnetic field, it can be determined through observation whether the magnetic field tends to grow or not in reaction 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 non-dynamo.

The membrane paradigm is a way of looking at black holes that allows for the material near their surfaces to be expressed in the language of dynamo theory.

As a spontaneous breakdown of 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 supersymmetric theory of stochastics, this supersymmetry is an intrinsic property of all stochastic differential equations, its meaning is the preservation of the continuity of the phase space of the model by continuous time flows, and its spontaneous breakdown is the stochastic generalization of the concept of deterministic chaos. In other words, kinematic dynamo is a manifestation of the chaoticity of the underlying flow of the background matter.

src: i.ytimg.com

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 Lorenz 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 ${\displaystyle \mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$ 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:

${\displaystyle \nabla \cdot \mathbf {B} =0}$

${\displaystyle \nabla \times \mathbf {B} =\mu _{0}\mathbf {J} }$

• The continuity equation for conservation of mass, for which the Boussinesq approximation is often used:

${\displaystyle \nabla \cdot \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:

${\displaystyle {\frac {D\mathbf {u} }{Dt}}=-{\frac {1}{\rho _{0}}}\nabla p+\nu \nabla ^{2}\mathbf {u} +\rho '\mathbf {g} +2\mathbf {\Omega } \times \mathbf {u} +\mathbf {\Omega } \times \mathbf {\Omega } \times \mathbf {R} +{\frac {1}{\rho _{0}}}\mathbf {J} \times \mathbf {B} ,}$

where ${\displaystyle \nu }$ is the kinematic viscosity, ${\displaystyle \rho _{0}}$is the mean density and ${\displaystyle \rho '}$ is the relative density perturbation that provides buoyancy (for thermal convection ${\displaystyle \rho '=\alpha \Delta T}$ where ${\displaystyle \alpha }$ is coefficient of thermal expansion), ${\displaystyle \Omega }$ is the rotation rate of the Earth, and ${\displaystyle \mathbf {J} }$ is the electric current density.

• A transport equation, usually of heat (sometimes of light element concentration):

${\displaystyle {\frac {\partial T}{\partial t}}=\kappa \nabla ^{2}T+\epsilon }$

where T is temperature, ${\displaystyle \kappa =k/\rho c_{p}}$ is the thermal diffusivity with k thermal conductivity, ${\displaystyle c_{p}}$ heat capacity, and ${\displaystyle \rho }$ density, and ${\displaystyle \epsilon }$ 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,

${\displaystyle Ra={\frac {g\alpha TD^{3}}{\nu \kappa }},E={\frac {\nu }{\Omega D^{2}}},Pr={\frac {\nu }{\kappa }},Pm={\frac {\nu }{\eta }}}$

where Ra is the Rayleigh number, E the Ekman number, Pr and Pm the Prandtl and magnetic Prandtl number. Magnetic field scaling is often in Elsasser number units ${\displaystyle B=(\rho \Omega /\sigma )^{1/2}}$.

Energy conversion between magnetic and kinematic energy

The scalar product of the above form of Navier-Stokes equation with ${\displaystyle \rho _{0}\mathbf {u} }$ gives the rate of increase of kinetic energy density, ${\displaystyle (1/2)\rho _{0}u^{2}}$, on the left-hand side. The last term on the right-hand side is then ${\displaystyle \mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$, the local contribution to the kinetic energy due to Lorentz force.

The scalar product of the induction equation with ${\displaystyle (1/\mu _{0})\mathbf {B} }$ gives the rate of increase of the magnetic energy density, ${\displaystyle (1/2\mu _{0})B^{2}}$, on the left-hand side. The last term on the right-hand side is then ${\displaystyle (1/\mu _{0})\mathbf {B} \cdot \left(\nabla \times (\mathbf {u} \times \mathbf {B} )\right)}$. 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 ${\displaystyle -\mathbf {u} \cdot \left((1/\mu _{0})(\nabla \times \mathbf {B} )\times \mathbf {B} )\right)=-\mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$ (where one of Maxwell's equations was used). This is the local contribution to the magnetic energy due to fluid motion.

Thus the term ${\displaystyle -\mathbf {u} \cdot (\mathbf {J} \times \mathbf {B} )}$ 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.

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 ${\displaystyle \mathbf {J} \times \mathbf {B} }$ on the outer core matter, whose velocity is ${\displaystyle \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 x10^-2 pascal-second ) due to its liquidity.

Thus the main time-averaged contribution to the work is from Coriolis force, whose size is ${\displaystyle -2\rho \,\mathbf {\Omega } \times \mathbf {u} }$, though this quantity and ${\displaystyle \mathbf {J} \times \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 J 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, ${\displaystyle J\,B\sim \rho \,\Omega \,u}$ and ${\displaystyle J\sim \sigma uB}$, giving ${\displaystyle \sigma \,u\,B^{2}\sim \rho \,\Omega \,u}$, or:

${\displaystyle B\sim {\sqrt {\frac {\rho \,\Omega }{\sigma }}}}$

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 ${\displaystyle \mathbf {\Omega } }$.

For earth outer core, ? is approximately 10⁴ kg/m³, ?=2?/day = 7.3x10^-5 seconds and ? is approximately 10⁷?^-1m^-1. This gives 2.7x10^-4 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 (R_{outer core}/R_Earth)³ = (2890/6370)³ = 0.093, giving 2.5x10^-5 Tesla, not far from to the measured value of 3x10^-5 Tesla at the equator.

src: cdn.vox-cdn.com

Numerical models

The equations for the geodynamo are enormously difficult to solve, and the realism of the solutions is limited mainly by computer power. For decades, theorists were confined to kinematic dynamo models described above, in which the fluid motion is chosen in advance and the effect on the magnetic field calculated. Kinematic dynamo theory was mainly a matter of trying different flow geometries and seeing whether they could sustain a dynamo.

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 received significant attention because it successfully reproduced some of the characteristics of the Earth's field, including geomagnetic reversals.

src: helios.augustana.edu

See also

• Antidynamo theorem
• Rotating magnetic field
• Secular variation

src: cdn.vox-cdn.com

References

Source of the article : Wikipedia