# Plasmas/Magnetohydrodynamics

(Redirected from Magnetohydrodynamics)
The Mitsubishi experimental boat Yamato 1 with magnetohydrodynamic drive is now on display in front of the Maritime Museum in Kobe, Japan. Credit: Geofrog.

“The word magnetohydrodynamics (MHD) is derived from magneto- meaning magnetic field, and hydro- meaning liquid, and -dynamics meaning movement.”[1]

## Theoretical magnetohydrodynamics

Def. "the study of the interaction of electrically conducting fluids with magnetic fields"[2] is called magnetohydrodynamics (MHD).

## Plasma object entities

Def. "the fraction of a given part [or the whole] of [a stellar surface] covered by magnetic field"[3] is called the magnetic filling factor.

## Strong forces

"The incompressible MHD equations are

${\displaystyle {\begin{array}{lcl}{\frac {\partial \mathbf {u} }{\partial t}}+\mathbf {u} \cdot \nabla \mathbf {u} &=&-\nabla p+\mathbf {B} \cdot \nabla \mathbf {B} +\nu \nabla ^{2}\mathbf {u} \\{\frac {\partial \mathbf {B} }{\partial t}}+\mathbf {u} \cdot \nabla \mathbf {B} &=&\mathbf {B} \cdot \nabla \mathbf {u} +\eta \nabla ^{2}\mathbf {B} \\\nabla \cdot \mathbf {u} &=&0\\\nabla \cdot \mathbf {B} &=&0.\end{array}}}$

where u, B, p represent the velocity, magnetic, and total pressure (thermal+magnetic) fields, ${\displaystyle \nu }$ and ${\displaystyle \eta }$ represent kinematic viscosity and magnetic diffusivity. The third equation is the incompressibility condition. In the above equation, the magnetic field is in Alfvén units (same as velocity units)."[4]

"The total magnetic field can be split into two parts: ${\displaystyle \mathbf {B} =\mathbf {B_{0}} +\mathbf {b} }$ (mean + fluctuations)."[4]

"The above equations in terms of Elsässer variables (${\displaystyle \mathbf {z} ^{\pm }=\mathbf {u} \pm \mathbf {b} }$) are"[4]

${\displaystyle {\frac {\partial {\mathbf {z} ^{\pm }}}{\partial t}}\mp \left(\mathbf {B} _{0}\cdot {\mathbf {\nabla } }\right){\mathbf {z} ^{\pm }}+\left({\mathbf {z} ^{\mp }}\cdot {\mathbf {\nabla } }\right){\mathbf {z} ^{\pm }}=-{\mathbf {\nabla } }p+\nu _{+}\nabla ^{2}\mathbf {z} ^{\pm }+\nu _{-}\nabla ^{2}\mathbf {z} ^{\mp }}$

where ${\displaystyle \nu _{\pm }=\nu \pm \eta }$. Nonlinear interactions occur between the Alfvénic fluctuations ${\displaystyle z^{\mp }}$.

"The important nondimensional parameters for MHD are"[4]

${\displaystyle {\begin{array}{lcl}{\text{Reynold number }}Re&=&UL/\nu \\{\text{Magnetic Reynold number }}Re_{M}&=&UL/\eta \\\end{array}}}$

"The Magnetic Reynolds number is a dimensionless group that gives an estimate of the effects of magnetic advection to magnetic diffusion"[5].

${\displaystyle {\begin{array}{lcl}{\text{Magnetic Prandtl number }}P_{M}&=&\nu /\eta .\end{array}}}$

"The magnetic Prandtl number is an important property of the fluid. Liquid metals have small magnetic Prandtl numbers, for example, liquid sodium's ${\displaystyle P_{M}}$ is around ${\displaystyle 10^{-5}}$. But plasmas have large ${\displaystyle P_{M}}$."[4]

"The Reynolds number is the ratio of the nonlinear term ${\displaystyle \mathbf {u} \cdot \nabla \mathbf {u} }$ of the Navier-Stokes equation to the viscous term. While the magnetic Reynolds number is the ratio of the nonlinear term and the diffusive term of the induction equation."[4]

"In many practical situations, the Reynolds number ${\displaystyle Re}$ of the flow is quite large. For such flows typically the velocity and the magnetic fields are random. Such flows are called to exhibit MHD turbulence. Note that ${\displaystyle Re_{M}}$ need not be large for MHD turbulence. ${\displaystyle Re_{M}}$ plays an important role in dynamo (magnetic field generation) problem."[4]

"The mean magnetic field plays an important role in MHD turbulence, for example it can make the turbulence anisotropic; suppress the turbulence by decreasing energy cascade etc. The earlier MHD turbulence models assumed isotropy of turbulence, while the later models have studied anisotropic aspects. In the following discussions will summarize these models. More discussions on MHD turbulence can be found in Biskamp[6] and Verma[7]."[4]

## Isotropic turbulence

"Iroshnikov[8] and Kraichnan[9] formulated the first phenomenological theory of MHD turbulence. They argued that in the presence of a strong mean magnetic field, ${\displaystyle z^{+}}$ and ${\displaystyle z^{-}}$ wavepackets travel in opposite directions with the phase velocity of ${\displaystyle B_{0}}$, and interact weakly. The relevant time scale is Alfven time ${\displaystyle (B_{0}k)^{-1}}$. As a results the energy spectra is

${\displaystyle E^{u}(k)\approx E^{b}(k)\approx A(\Pi V_{A})^{1/2}k^{-3/2}.}$

where ${\displaystyle \Pi }$ is the energy cascade rate."[4]

"Later Dobrowolny et al. [10]derived the following generalized formulas for the cascade rates of ${\displaystyle z^{\pm }}$ variables:

${\displaystyle \Pi ^{+}\approx \Pi ^{-}\approx \tau _{k}^{\pm }E^{+}(k)E^{-}(k)k^{4}\approx E^{+}(k)E^{-}(k)k^{3}/B_{0}}$

where ${\displaystyle \tau ^{\pm }}$ are the interaction time scales of ${\displaystyle z^{\pm }}$ variables."[4]

"Iroshnikov and Kraichnan's phenomenology follows once we choose ${\displaystyle \tau ^{\pm }\approx 1/(kV_{A})}$."[4]

"Marsch[11] chose the nonlinear time scale ${\displaystyle T_{NL}^{\pm }\approx (kz_{k}^{\mp })^{-1}}$ as the interaction time scale for the eddies and derived Kolmogorov-like energy spectrum for the Elsasser variables:

${\displaystyle E^{\pm }(k)=K^{\pm }(\Pi ^{\pm })^{4/3}(\Pi ^{\mp })^{-2/3}k^{-5/3}}$

where ${\displaystyle \Pi ^{+}}$ and ${\displaystyle \Pi ^{-}}$ are the energy cascade rates of ${\displaystyle z^{+}}$ and ${\displaystyle z^{-}}$ respectively, and ${\displaystyle K^{\pm }}$ are constants."[4]

"Matthaeus and Zhou[12] attempted to combine the above two time scales by postulating the interaction time to be the harmonic mean of Alfven time and nonlinear time.[4]

"The main difference between the two competing phenomenologies (-3/2 and -5/3) is the chosen time scales for the interaction time. The main underlying assumption in that Iroshnikov and Kraichnan's phenomenology should work for strong mean magnetic field, whereas Marsh's phenomenology should work when the fluctuations dominate the mean magnetic field (strong turbulence)."[4]

## Anisotropic turbulence

"Mean magnetic field makes turbulence anisotropic. This aspect has been studied in last two decades. In the limit ${\displaystyle \delta z^{\pm }\ll B_{0}}$, Galtier et al.[13] showed using kinetic equations that

${\displaystyle E(k)\sim (\Pi B_{0})^{1/2}k_{||}^{1/2}k_{\perp }^{-2}}$

where ${\displaystyle k_{||}}$ and ${\displaystyle k_{\perp }}$ are components of the wavenumber parallel and perpendicular to mean magnetic field. The above limit is called the weak turbulence limit."[4]

"Under the strong turbulence limit, ${\displaystyle \delta z^{\pm }\sim B_{0}}$, Goldereich and Sridhar[14] argue that ${\displaystyle k_{\perp }z_{k_{\perp }}\sim k_{||}B_{0}}$ ("critical balanced state") which implies that"[4]

${\displaystyle {\begin{array}{lcl}E(k)&\propto &k_{\perp }^{-5/3};\\k_{||}&\propto &k_{\perp }^{2/3}\end{array}}}$

"The above anisotropic turbulence phenomenology has been extended for large cross helicity MHD."[4]

## Electromagnetics

"The Stuart number (N, also known as magnetic interaction parameter) is a dimensionless number of fluids, i.e. gases or liquids."[15]

"It is defined as the ratio of electromagnetic to inertial forces, which gives an estimate of the relative importance of a magnetic field on a flow. The Stuart number is relevant for flows of conducting fluids, e.g. in fusion reactors, steel casters or plasmas.[16]"[15]

${\displaystyle {\it {N}}={\frac {Ha^{2}}{Re}}={\frac {B^{2}\cdot L_{c}\cdot \sigma }{\rho U}}}$
• ${\displaystyle Ha}$ – Hartmann number
• ${\displaystyle Re}$ – Reynolds number
• ${\displaystyle B}$ – magnetic flux density
• ${\displaystyle L_{c}}$ – characteristic length
• ${\displaystyle \sigma }$ – electric conductivity
• ${\displaystyle \mu }$ – dynamic viscosity, sometimes denoted as ${\displaystyle \eta }$

## Electrohydrodynamics

Electrohydrodynamics (EHD), also known as electro-fluid-dynamics (EFD) or electrokinetics, is the study of the dynamics of electrically charged fluids.[17] It is the study of the motions of ionised particles or molecules and their interactions with electric fields and the surrounding fluid. The term may be considered to be synonymous with the rather elaborate electrostrictive hydrodynamics. EHD covers the following types of particle and fluid transport mechanisms: Electrophoresis, electrokinesis, dielectrophoresis, electro-osmosis, and electrorotation. In general, the phenomena relate to the direct conversion of electrical energy into kinetic energy, and vice versa."[18]

"In the first instance, shaped electrostatic fields create hydrostatic pressure (or motion) in dielectric media. When such media are fluids, a flow is produced. If the dielectric is a vacuum or a solid, no flow is produced. Such flow can be directed against the electrodes, generally to move the electrodes. In such case, the moving structure acts as an electric motor.”[18]

## Hartmann numbers

A Hartmann number is the ratio of electromagnetic force to the viscous force first introduced by Hartmann.[19] It is defined by:

${\displaystyle H_{a}=BL{\sqrt {\frac {\sigma }{\mu }}}}$

where

• ${\displaystyle B}$ is the magnetic field
• ${\displaystyle L}$ is the characteristic length scale
• ${\displaystyle \sigma }$ is the electrical conductivity
• ${\displaystyle \mu }$ is the viscosity[20]

## Computational magnetohydrodynamics

"Computational magnetohydrodynamics (CMHD) is a rapidly developing branch of magnetohydrodynamics that uses numerical methods and algorithms to solve and analyze problems that involve electrically conducting fluids. Most of the methods used in CMHD are borrowed from the well established techniques employed in Computational fluid dynamics. The complexity mainly arises due to the presence of a magnetic field and its coupling with the fluid. One of the important issues is to numerically maintain the ${\displaystyle \nabla \cdot {\mathbf {B} }=0}$ (conservation of magnetic flux) condition, from Maxwell's equations, to avoid any unphysical effects."[21]

"Numerical resistivity is a problem in computer simulations of ideal magnetohydrodynamics (MHD). It is a form of numerical diffusion. In near-ideal MHD systems, the magnetic field can diffuse only very slowly through the plasma or fluid of the system; it is rate-limited by the resistivity of the fluid. In Eulerian simulations where the field is arbitrarily aligned compared to the simulation grid, the numerical diffusion rate takes the form similar to an additional resistivity, causing non-physical and sometimes bursty magnetic reconnection in the simulation. Numerical resistivity is a function of resolution, alignment of the magnetic field with the grid, and numerical method. In general, numerical resistivity will not behave isotropically, and there can be different effective numerical resistivities in different parts of the computational domain. For current (2005) simulations of the solar corona and inner heliosphere, this numerical effect can be several orders of magnitude larger than the physical resistivity of the plasma."[22]

## Weak forces

"[M]otions resulting from [a linear magnetohydrodynamic] instability act as a dynamo to sustain the magnetic field."[23] "Supersonic flows are initially generated by the Balbus-Hawley magnetic shear instability."[23]

"A plasma with local magnetohydrodynamic instabilities creates mechanical turbulence, motion, or shear (a dynamo) which in turn generates or sustains the local magnetic field."[24]

"Magnetic field generation by dynamo action is often studied within the theoretical framework of magnetohydrodynamics (MHD)."[25]

"For small values of the Hall parameter, the small-scale dynamo is more efficient, displaying faster growth and saturating at larger amplitudes of the magnetic field. For larger values of the Hall parameter, saturation of the magnetic field is reached at smaller amplitudes than in the MHD case."[25]

When this magnetohydrodynamic dynamo occurs between or within radiative layers, a radiative dynamo is operating.

"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."[26]

"Numerical models are used to simulate fully nonlinear dynamos. A minimum of 5 equations are needed. They are as follows. The induction equation, see above. Maxwell's equation:"[26]

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

"The (sometimes) Boussinesq conservation of mass:"[26]

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

"The (sometimes) Boussinesq conservation of momentum, also known as the Navier-Stokes equation:

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

where ${\displaystyle \nu }$ is the kinematic viscosity, ${\displaystyle \rho ^{'}}$ is the density perturbation that provides buoyancy (for thermal convection ${\displaystyle \rho ^{'}=\alpha \Delta T}$, ${\displaystyle \Omega }$ is the rotation rate of the Earth, and ${\displaystyle \mathbf {J} }$ is the electrical current density."[26]

"Finally, 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."[26]

"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 }$."[26]

## Antidynamo theorem

Def. “one of several results that restrict the type of magnetic fields that may be produced by dynamo action”[27] is called an antidynamo theorem.

Def. “no axisymmetric magnetic field can be maintained through a self-sustaining dynamo action by an axially symmetric current”[27][28] is called Cowling's antidynamo theorem.

“The impact of the known anti-dynamo theorems is that successful dynamos do not possess a high degree of symmetry.”[27]

## Plasma objects

This is a photo of the kink instability in action - the 3 by 25 cm pyrex tube at Aldermaston. Credit: Alan Sykes, UK Atomic Energy Authority.

"Plasma instabilities can be divided into two general groups:

1. hydrodynamic instabilities
2. kinetic instabilities."[29]

"Plasma instabilities are also categorised into different modes:[30]"[29]

 Mode (azimuthal wave number) Note Description Radial modes Description m=0 Sausage instability: displays harmonic variations of beam radius with distance along the beam axis n=0 Axial hollowing n=1 Standard sausaging n=2 Axial bunching m=1 Sinuous, kink or hose instability: represents transverse displacements of the beam cross-section without change in the form or in a beam characteristics other than the position of its center of mass m=2 Filamentation modes: growth leads towards the breakup of the beam into separate filaments. Gives an elliptic cross-section m=3 Gives a pyriform (pear-shaped) cross-section

"Ideal MHD instabilities driven by current or pressure gradients".[29]

"A kink instability, also oscillation or mode, is a class of magnetohydrodynamic instabilities which sometimes develop in a thin plasma column carrying a strong axial current. If a "kink" begins to develop in a column the magnetic forces on the inside of the kink become larger than those on the outside, which leads to growth of the perturbation.[31] As it develops at fixed areas in the plasma, kinks belong to the class of "absolute plasma instabilities", as opposed to convective processes."[32]

## Space plasma characteristics

"Space plasma pioneers Hannes Alfvén and Carl-Gunne Fälthammar divided the plasmas in the solar system into three different categories:"[33]

Classification of Magnetic Cosmic Plasmas

 Characteristic Space plasma density categories (Note that density does not refer to only particle density) Ideal comparison High density Medium Density Low Density Criterion λ << ρ λ << ρ << lc lc << λ lc << λD Examples Stellar interior Solar photosphere Solar chromosphere/corona Interstellar/intergalactic space Ionosphere above 70 km Magnetosphere during magnetic disturbance. Interplanetary space Single charges in a high vacuum Diffusion Isotropic Anisotropic Anisotropic and small No diffusion Conductivity Isotropic Anisotropic Not defined Not defined Electric field parallel to B in completely ionized gas Small Small Any value Any value Particle motion in plane perpendicular to B Almost straight path between collisions Circle between collisions Circle Circle Path of guiding centre parallel to B Straight path between collisions Straight path between collisions Oscillations (e.g. between mirror points) Oscillations (e.g. between mirror points) Debye Distance λD λD << lc λD << lc λD << lc λD >> lc Magnetohydrodynamics suitability Yes Approximately No No

λ=Mean free path. ρ= Larmor radius (gyroradius) of electron. λD=Debye length. lc=Characteristic length
Adapted From Cosmical Electrodynamics (2nd Ed. 1952) Alfvén and Fälthammar

## Gaseous objects

"Supersonic flows are initially generated by the Balbus-Hawley magnetic shear instability."[23]

“The magnetorotational instability or MRI is a fluid instability that arises when the angular velocity of a magnetized fluid decreases as the distance from the rotation center increases. It can also be known as the Velikhov-Chandrasekhar instability or Balbus-Hawley instability in the literature; not to be confused with the Velikhov instability which is the electrothermal instability. The MRI is of particular relevance in astrophysics where it is an important part of the dynamics in accretion discs. Balbus and Hawley were the first to realize the astrophysical importance of this instability and explain its physical mechanism.”[34]

“A rotating hydrodynamic fluid disc (e.g. a nonmagnetic accretion disc) will remain in a laminar flow state as long as the angular momentum (per unit mass) increases outwards. This is also known as the Rayleigh stability criterion:

${\displaystyle {\frac {\partial (r^{2}\Omega )}{\partial r}}>0}$

where ${\displaystyle \Omega }$ is the angular velocity of a fluid element and ${\displaystyle r}$ is its distance to the rotation center. Instabilities leading ultimately to turbulence occur if a magnetic field is present and the angular velocity decreases with radius:

${\displaystyle {\frac {\partial \Omega }{\partial r}}<0}$[34].

“The MRI was first noticed in a non-astrophysical context by Evgeny Velikhov in 1959 when considering the stability of Couette flow of an ideal hydromagnetic fluid.[35] His result was later generalized by S. Chandrasekhar in 1960.[36] This mechanism was proposed by Acheson & Hide (1973) to perhaps play a role in the context of the Earth's geodynamo problem.[37] The applicability of this mechanism to the problem of accretion disks was appreciated when S. A. Balbus and J. F. Hawley established that weak magnetic fields can substantially alter the stability character of accretion disks.[38][34]

“The observed accretion rates in astrophysical objects cannot be explained by a molecular viscosity, as the outward angular momentum transport in that case would not be enough to account for the inward flow of mass. The MRI provides a mechanism to account for the additional outward angular momentum transport. It is intrinsically a magnetohydrodynamic (MHD) phenomenon, with no hydrodynamic analog. Dynamo action is usually invoked to justify the existence of the magnetic field required for the instability to set in.[39][34]

## Liquid objects

“The resistive skin time is a characteristic time of typical magnetohydrodynamic (MHD) phenomena.

The resistive skin time is defined as: ${\displaystyle \tau _{R}={\frac {\mu _{0}a^{2}}{\eta }}}$ where ${\displaystyle \eta }$ is the resistivity, ${\displaystyle a}$ is a typical radius of the device and ${\displaystyle \mu _{0}}$ is the magnetic permeability.”[40]

## Waves

"There are several distinct kinds of MHD modes which have quite different dispersive, polarisation, and propagation properties:

• Kink (or transverse) modes, which are oblique fast magnetoacoustic (also known as magnetosonic waves) guided by the plasma structure; the mode causes the displacement of the axis of the plasma structure. These modes are weakly compressible, but could nevertheless be observed with imaging instruments as periodic standing or propagating displacements of coronal structures, e.g. coronal loops. The frequency of transverse or "kink" modes is given by the following expression:
${\displaystyle \omega _{K}={\sqrt {\frac {2k_{z}B^{2}}{\mu (\rho _{i}+\rho _{e})}}}}$

For kink modes the parameter ${\displaystyle m}$ is equal to 1.

• Sausage modes, which are also oblique fast magnetoacoustic waves guided by the plasma structure; the mode causes expansions and contractions of the plasma structure, but does not displace its axis. These modes are compressible and cause significant variation of the absolute value of the magnetic field in the oscillating structure. The frequency of sausage modes is given by the following expression:
${\displaystyle \omega _{S}={\sqrt {\frac {k_{z}^{2}B^{2}}{\mu \rho _{e}}}}}$

For sausage modes the parameter ${\displaystyle m}$ is equal to 0.

• Longitudinal (or slow, or acoustic) modes, which are slow magnetoacoustic waves propagating mainly along the magnetic field in the plasma structure; these mode are essentially compressible. The magnetic field perturbation in these modes is negligible. The frequency of slow modes is given by the following expression:
${\displaystyle \omega _{L}={\sqrt {k_{z}^{2}\left({\frac {C_{s}^{2}C_{A}^{2}}{C_{s}^{2}+C_{A}^{2}}}\right)}}}$

Where we define ${\displaystyle C_{s}}$ as the sound speed and ${\displaystyle C_{A}}$ as the Alfvèn speed.

• Torsional (Alfvén or twist) modes are incompressible transverse perturbations of the magnetic field along certain individual magnetic surfaces. In contrast with kink modes, torsional modes cannot be observed with imaging instruments, as they do not cause the displacement of either the structure axis or its boundary.
${\displaystyle \omega _{A}={\sqrt {\frac {k_{z}^{2}B^{2}}{\mu \rho _{i}}}}}$"[41]

"Waves in plasmas are an interconnected set of particles and fields which propagates in a periodically repeating fashion. A plasma is a quasineutral, electrically conductive fluid. In the simplest case, it is composed of electrons and a single species of positive ions, but it may also contain multiple ion species including negative ions as well as neutral particles. Due to its electrical conductivity, a plasma couples to electric and magnetic fields. This complex of particles and fields supports a wide variety of waves."[42]

"Waves in plasmas can be classified as electromagnetic or electrostatic according to whether or not there is an oscillating magnetic field. Applying Faraday's law of induction to plane waves, we find ${\displaystyle \mathbf {k} \times {\tilde {\mathbf {E} }}=\omega {\tilde {\mathbf {B} }}}$, implying that an electrostatic wave must be purely longitudinal. An electromagnetic wave, in contrast, must have a transverse component, but may also be partially longitudinal."[42]

"Waves can be further classified by the oscillating species. In most plasmas of interest, the electron temperature is comparable to or larger than the ion temperature. This fact, coupled with the much smaller mass of the electron, implies that the electrons are much faster than the ions. An electron mode depends on the mass of the electrons, but the ions may be assumed to be infinitely massive, i.e. stationary. An ion mode depends on the ion mass, but the electrons are assumed to be massless and to redistribute themselves instantaneously according to the Boltzmann relation. Only rarely, e.g. in the lower hybrid oscillation, will a mode depend on both the electron and the ion mass."[42]

"The various modes can also be classified according to whether they propagate in an unmagnetized plasma or parallel, perpendicular, or oblique to the stationary magnetic field. Finally, for perpendicular electromagnetic electron waves, the perturbed electric field can be parallel or perpendicular to the stationary magnetic field."[42]

Summary of elementary plasma waves[42]
EM character oscillating species conditions dispersion relation name
electrostatic electrons ${\displaystyle {\vec {B}}_{0}=0\ {\rm {or}}\ {\vec {k}}\|{\vec {B}}_{0}}$ ${\displaystyle \omega ^{2}=\omega _{p}^{2}+3k^{2}v_{th}^{2}}$ plasma oscillation (or Langmuir wave)
${\displaystyle {\vec {k}}\perp {\vec {B}}_{0}}$ ${\displaystyle \omega ^{2}=\omega _{p}^{2}+\omega _{c}^{2}=\omega _{h}^{2}}$ upper hybrid oscillation
ions ${\displaystyle {\vec {B}}_{0}=0\ {\rm {or}}\ {\vec {k}}\|{\vec {B}}_{0}}$ ${\displaystyle \omega ^{2}=k^{2}v_{s}^{2}=k^{2}{\frac {\gamma _{e}KT_{e}+\gamma _{i}KT_{i}}{M}}}$ ion acoustic wave
${\displaystyle {\vec {k}}\perp {\vec {B}}_{0}}$ (nearly) ${\displaystyle \omega ^{2}=\Omega _{c}^{2}+k^{2}v_{s}^{2}}$ electrostatic ion cyclotron wave
${\displaystyle {\vec {k}}\perp {\vec {B}}_{0}}$ (exactly) ${\displaystyle \omega ^{2}=[(\Omega _{c}\omega _{c})^{-1}+\omega _{i}^{-2}]^{-1}}$ lower hybrid oscillation
electromagnetic electrons ${\displaystyle {\vec {B}}_{0}=0}$ ${\displaystyle \omega ^{2}=\omega _{p}^{2}+k^{2}c^{2}}$ light wave
${\displaystyle {\vec {k}}\perp {\vec {B}}_{0},\ {\vec {E}}_{1}\|{\vec {B}}_{0}}$ ${\displaystyle {\frac {c^{2}k^{2}}{\omega ^{2}}}=1-{\frac {\omega _{p}^{2}}{\omega ^{2}}}}$ O wave
${\displaystyle {\vec {k}}\perp {\vec {B}}_{0},\ {\vec {E}}_{1}\perp {\vec {B}}_{0}}$ ${\displaystyle {\frac {c^{2}k^{2}}{\omega ^{2}}}=1-{\frac {\omega _{p}^{2}}{\omega ^{2}}}\,{\frac {\omega ^{2}-\omega _{p}^{2}}{\omega ^{2}-\omega _{h}^{2}}}}$ X wave
${\displaystyle {\vec {k}}\|{\vec {B}}_{0}}$ (right circ. pol.) ${\displaystyle {\frac {c^{2}k^{2}}{\omega ^{2}}}=1-{\frac {\omega _{p}^{2}/\omega ^{2}}{1-(\omega _{c}/\omega )}}}$ R wave (whistler mode)
${\displaystyle {\vec {k}}\|{\vec {B}}_{0}}$ (left circ. pol.) ${\displaystyle {\frac {c^{2}k^{2}}{\omega ^{2}}}=1-{\frac {\omega _{p}^{2}/\omega ^{2}}{1+(\omega _{c}/\omega )}}}$ L wave
ions ${\displaystyle {\vec {B}}_{0}=0}$   none
${\displaystyle {\vec {k}}\|{\vec {B}}_{0}}$ ${\displaystyle \omega ^{2}=k^{2}v_{A}^{2}}$ Alfvén wave
${\displaystyle {\vec {k}}\perp {\vec {B}}_{0}}$ ${\displaystyle {\frac {\omega ^{2}}{k^{2}}}=c^{2}\,{\frac {v_{s}^{2}+v_{A}^{2}}{c^{2}+v_{A}^{2}}}}$ magnetosonic wave

${\displaystyle \omega }$ - wave frequency, ${\displaystyle k}$ - wave number, ${\displaystyle c}$ - speed of light, ${\displaystyle \omega _{p}}$ - plasma frequency, ${\displaystyle \omega _{i}}$ - ion plasma frequency, ${\displaystyle \omega _{c}}$ - electron gyrofrequency, ${\displaystyle \Omega _{c}}$ - proton gyrofrequency, ${\displaystyle \omega _{h}}$ - upper hybrid frequency, ${\displaystyle v_{s}}$ - plasma "sound" speed, ${\displaystyle v_{A}}$ - plasma Alfven speed

## Rocky objects

"The jump conditions across an MHD shock or discontinuity are referred as the Rankine-Hugoniot equations for MHD. In the frame moving with the shock/discontinuity, those jump conditions can be written:

${\displaystyle \rho _{1}v_{n1}=\rho _{2}v_{n2},}$
${\displaystyle B_{n1}=B_{n2},}$
${\displaystyle \rho _{1}v_{n1}^{2}+p_{1}+{\frac {B_{t1}^{2}}{2\mu _{0}}}=\rho _{2}v_{n2}^{2}+p_{2}+{\frac {B_{t2}^{2}}{2\mu _{0}}},}$
${\displaystyle \rho _{1}v_{n1}\mathbf {v_{t1}} -{\frac {\mathbf {B_{t1}} B_{n1}}{\mu _{0}}}=\rho _{2}v_{n2}\mathbf {v_{t2}} -{\frac {\mathbf {B_{t2}} B_{n2}}{\mu _{0}}},}$
${\displaystyle \left({\frac {\gamma }{\gamma -1}}{\frac {p_{1}}{\rho _{1}}}+{\frac {v_{1}^{2}}{2}}\right)\rho _{1}v_{n1}+{\frac {v_{n1}B_{t1}^{2}}{\mu _{0}}}-{\frac {B_{n1}(\mathbf {B_{t1}} \cdot \mathbf {v_{t1}} )}{\mu _{0}}}=\left({\frac {\gamma }{\gamma -1}}{\frac {p_{2}}{\rho _{2}}}+{\frac {v_{2}^{2}}{2}}\right)\rho _{2}v_{n2}+{\frac {v_{n2}B_{t2}^{2}}{\mu _{0}}}-{\frac {B_{n2}(\mathbf {B_{t2}} \cdot \mathbf {v_{t2}} )}{\mu _{0}}},}$
${\displaystyle (\mathbf {v} \times \mathbf {B} )_{t1}=(\mathbf {v} \times \mathbf {B} )_{t2},}$

where ${\displaystyle \rho }$, v, p, B are the plasma density, velocity, (thermal) pressure and magnetic field respectively. The subscripts t and n refer to the tangential and normal components of a vector (with respect to the shock/discontinuity front). The subscripts 1 and 2 refer to the two states of the plasma on each side of the shock/discontinuity".[43]

## Shocks

"Shocks are transition layers across which there is a transport of particles. There are three types of shocks in MHD: slow-mode, intermediate and fast-mode shocks."[43]

"Intermediate shocks are non-compressive (meaning that the plasma density does not change across the shock). A special case of the intermediate shock is referred to as a rotational discontinuity. They are isentropic. All thermodynamic quantities are continuous across the shock, but the tangential component of the magnetic field can "rotate". Intermediate shocks in general however, unlike rotational discontinuities, can have a discontinuity in the pressure."[43]

"Slow-mode and fast-mode shocks are compressive and are associated with an increase in entropy. Across slow-mode shock, the tangential component of the magnetic field decreases. Across fast-mode shock it increases."[43]

"The type of shocks depend on the relative magnitude of the upstream velocity in the frame moving with the shock with respect to some characteristic speed. Those characteristic speeds, the slow and fast magnetosonic speeds, are related to the Alfvén speed, ${\displaystyle v_{A}}$ and the sonic speed, ${\displaystyle c_{s}}$ as follows:

${\displaystyle a_{\mathrm {slow} }^{2}={\frac {1}{2}}\left[\left(c_{s}^{2}+V_{A}^{2}\right)-{\sqrt {\left(c_{s}^{2}+V_{A}^{2}\right)^{2}-4c_{s}^{2}V_{A}^{2}\cos ^{2}\theta _{Bn}}}\,\right],}$
${\displaystyle a_{\mathrm {fast} }^{2}={\frac {1}{2}}\left[\left(c_{s}^{2}+V_{A}^{2}\right)+{\sqrt {\left(c_{s}^{2}+V_{A}^{2}\right)^{2}-4c_{s}^{2}V_{A}^{2}\cos ^{2}\theta _{Bn}}}\,\right],}$

where ${\displaystyle V_{A}}$ is the Alfvén speed and ${\displaystyle \theta _{Bn}}$ is the angle between the incoming magnetic field and the shock normal vector."[43]

"The normal component of the slow shock propagates with velocity ${\displaystyle a_{\mathrm {slow} }}$ in the frame moving with the upstream plasma, that of the intermediate shock with velocity ${\displaystyle V_{An}}$ and that of the fast shock with velocity ${\displaystyle a_{\mathrm {fast} }}$. The fast mode waves have higher phase velocities than the slow mode waves because the density and magnetic field are in phase, whereas the slow mode wave components are out of phase."[43]

## Magnetic tension

"The magnetic tension force is a restoring force (SI unit: Pa·m−1) that acts to straighten bent magnetic field lines. It equals:"[44]

${\displaystyle {\frac {\left(\mathbf {B} \cdot \nabla \right)\mathbf {B} }{\mu _{0}}}}$

"It is analogous to rubber bands and their restoring force. The force is directed antiradially. Although magnetic tension is referred to as a force, it is actually a pressure gradient (Pa m−1) which is also a force density (N m−3)."[44]

"The magnetic pressure is the energy density of the magnetic field and it increases as magnetic field lines convene with each other. In contrast, magnetic tension force is determined by how much the magnetic pressure changes with distance. Magnetic tension forces also rely on vector current densities ${\displaystyle \mathbf {J} }$ and their interaction with the magnetic field ${\displaystyle \mathbf {B} }$. Plotting magnetic tension along adjacent field lines can give a picture as to their divergence and convergence with respect to each other as well as current densities ${\displaystyle \mathbf {J} }$."[44]

"Magnetic tension is particularly important in plasma physics and magnetohydrodynamics, where it controls dynamics of some systems and the shape of magnetized structures. In magnetohydrodynamics, the magnetic tension force can be derived from the momentum equation

${\displaystyle \rho \left({\frac {\partial }{\partial t}}+\mathbf {V} \cdot \nabla \right)\mathbf {V} =\mathbf {J} \times \mathbf {B} -\nabla p}$

using the relation ${\displaystyle \mu _{0}\mathbf {J} =\nabla \times \mathbf {B} }$. The first term on the right hand side of the above equation represents electromagnetic forces and the second term represents pressure gradient forces."[44]

"[A] Taylor state is the minimum energy state of a plasma satisfying the constraint of conserving magnetic helicity.[45]"[46]

"Consider a closed, simply-connected, flux-conserving, perfectly conducting surface ${\displaystyle S}$ surrounding a plasma with negligible thermal energy (${\displaystyle \beta \rightarrow 0}$)."[46]

"Since ${\displaystyle {\vec {B}}.{\vec {ds}}=0}$ on ${\displaystyle S}$. This implies that ${\displaystyle {\vec {A}}_{||}=0}$."[46]

"As discussed above, the plasma would relax towards a minimum energy state while conserving its magnetic helicity. Since the boundary is perfectly conducting, there cannot be any change in the associated flux. This implies ${\displaystyle \delta {\vec {B}}.{\vec {ds}}=0}$ and ${\displaystyle \delta {\vec {A}}_{||}=0}$ on ${\displaystyle S}$."[46]

"We formulate a variational problem of minimizing the plasma energy ${\displaystyle W=\int d^{3}rB^{2}/2\mu _{\circ }}$ while conserving magnetic helicity ${\displaystyle K=\int d^{3}r{\vec {A}}.{\vec {B}}}$."[46]

"The variational problem is ${\displaystyle \delta W-\lambda \delta K=0}$.

After some algebra this leads to the following constraint for the minimum energy state ${\displaystyle \nabla \times {\vec {B}}=\lambda {\vec {B}}}$."[46]

## Sun (star)

Main sources: Stars/Sun and Sun (star)

"The solar dynamo is the physical process that generates the Sun's magnetic field. The Sun is permeated by an overall dipole magnetic field, as are many other celestial bodies such as the Earth. The dipole field is produced by a circular electric current flowing deep within the star, following Ampère's law. The current is produced by shear (stretching of material) between different parts of the Sun that rotate at different rates, and the fact that the Sun itself is a very good electrical conductor (and therefore governed by the laws of magnetohydrodynamics)."[47]

"Any electrically conducting fluid can form a dynamo simply by shear within the fluid itself, because of a consequence of Lenz's law of induction: moving the fluid through a pre-existing magnetic field will induce electrical currents in the fluid that distort the pre-existing magnetic field. The direction of the distortion is such that the existing field lines tend to be dragged along with the fluid, like threads of dye embedded in taffy or syrup. If the flow has a strong shear component then the individual field lines are stretched by the flow, amplifying the existing magnetic field. Such systems are called MHD dynamos."[47]

"Depending on the structure of the flow, the dynamo may be self-exciting and stable, self-exciting and chaotic, or decaying. The Sun's dynamo is self-exciting: the direction of the field reverses itself about every 11 years, causing the sunspot cycle as ropes of magnetic field lines rise to the surface of the Sun and manifest as sunspots on the surface. The solar dynamo is thought to be located at the tachocline, a region inside the Sun that exhibits a very large shear profile."[47]

## Coronal clouds

In a coronal cloud are magnetohydrodynamic plasma flux tubes along magnetic field lines.[48]

## Interplanetary magnetic field

Notation: let the symbol IMF stand for interplanetary magnetic field.

"[T]he IMF switched from northward to southward while Mariner [10] was in the Mercurian magnetosphere."[49]

## Magnetospheres

“A magnetosphere is formed when a stream of charged particles, such as the solar wind, interacts with and is deflected by the magnetic field of a planet or similar body.”[50]

“Planets which generate magnetic fields in their interiors ... are surrounded by invisible magnetospheres. ... [I]n many respects, the magnetosphere of Venus is a scaled-down version of Earth’s. ... Earth’s magnetosphere is 10 times larger [than that of Venus]”[51]

## Mercury

Main source: Mercury
The graph shows the relative strength of Mercury's magnetic field. Credit: .

"Despite its small size and slow 59-day-long rotation, Mercury has a significant, and apparently global, magnetic field. According to measurements taken by Mariner 10, it is about 1.1% as strong as the Earth's. The magnetic field strength at the Mercurian equator is about 300 nT.[52][53] Like that of Earth, Mercury's magnetic field is dipolar.[54] Unlike Earth, Mercury's poles are nearly aligned with the planet's spin axis.[55] Measurements from both the Mariner 10 and MESSENGER space probes have indicated that the strength and shape of the magnetic field are stable.[55]"[56]

"During its second flyby of the planet on October 6, 2008, MESSENGER discovered that Mercury's magnetic field can be extremely "leaky." The spacecraft encountered magnetic "tornadoes" – twisted bundles of magnetic fields connecting the planetary magnetic field to interplanetary space – that were up to 800 km wide or a third of the radius of the planet. These 'tornadoes' form when magnetic fields carried by the solar wind connect to Mercury's magnetic field. As the solar wind blows past Mercury's field, these joined magnetic fields are carried with it and twist up into vortex-like structures. These twisted magnetic flux tubes, technically known as flux transfer events, form open windows in the planet's magnetic shield through which the solar wind may enter and directly impact Mercury's surface.[57]"[56]

"The process of linking interplanetary and planetary magnetic fields, called magnetic reconnection, is common throughout the cosmos. It occurs in Earth's magnetic field, where it generates magnetic tornadoes as well. The MESSENGER observations show the reconnection rate is ten times higher at Mercury. Mercury's proximity to the Sun only accounts for about a third of the reconnection rate observed by MESSENGER.[57]"[56]

There are "observed MHD fluctuations in the magnetosphere"[49] of Mercury.

## Venus

Main source: Venus

"In 1967, Venera-4 found the Venusian magnetic field is much weaker than that of Earth. This magnetic field is induced by an interaction between the ionosphere and the solar wind,[58][59] ... Venus's small induced magnetosphere provides negligible protection to the atmosphere against cosmic radiation. This radiation may result in cloud-to-cloud lightning discharges.[60]"[61]

"The weak magnetosphere around Venus means the solar wind is interacting directly with the outer atmosphere of the planet. Here, ions of hydrogen and oxygen are being created by the dissociation of neutral molecules from ultraviolet radiation. The solar wind then supplies energy that gives some of these ions sufficient velocity to escape the planet's gravity field. This erosion process results in a steady loss of low-mass hydrogen, helium, and oxygen ions, while higher-mass molecules, such as carbon dioxide, are more likely to be retained."[61]

## Earth

Main source: Earth
This is a schematic diagram of Earth's interior. The outer core is the source of the geomagnetic field. Credit: NASA.
Reversals of Earth's magnetic field are diagrammed. Credit: NASA.
This is a schematic of the Birkeland currents and their connection to Earth's ionospheric current systems. Credit: Le, G., J. A. Slavin, and R. J. Strangeway.
Auroral-like Birkeland currents were created by scientist Kristian Birkeland in his terrella, featuring a magnetised anode globe in an evacuated chamber. Credit: Kristian Birkeland.
The complex self-constricting magnetic field lines and current paths in a Birkeland current that may develop in a plasma (Figure 15.3.2, Alfvén and Arrhenius, 1976)[62] Credit: .

"Beneath the Earth's mantle, lies the core which is made up of two parts - the solid inner core and liquid outer core - both have significant quantities of iron. The liquid outer core moves in the presence of the magnetic field and eddies are set up into the same due to the Coriolis effect. These eddies develop a magnetic field which boosts Earth's original magnetic field - a process which is self-sustaining and is called as the geomagnetic dynamo.[63]"[1]

"Based on the MHD equations, Glatzmaier and Paul Roberts have made a supercomputer model of the Earth's interior. After running the simulations for thousands of years in virtual time, the changes in Earth's magnetic field can be studied. The simulation results are in good agreement with the observations as the simulations have correctly predicted that the Earth's magnetic field flips every few thousands of years. During the flips, the magnetic field doesn't vanish altogether - it just gets more complicated.[64]"[1]

"The Atmospheric dynamo is a pattern of electrical currents that are set up in the Earth's ionosphere by multiple effects, mostly the Sun's solar wind, but also the tides of the Moon and Sun.[65][66] The currents flow in circuits between the poles and the equator, but they are not well understood.[67]"[68]

"Magnetite typically carries the dominant magnetic signature in rocks, and so it has been a critical tool in paleomagnetism, a science important in discovering and understanding plate tectonics and as historic data for magnetohydrodynamics and other scientific fields. The relationships between magnetite and other iron-rich oxide minerals such as ilmenite, hematite, and ulvospinel have been much studied, as the complicated [metamorphic] reactions between these minerals and oxygen influence how and when magnetite preserves records of the Earth's magnetic field."[69]

The schematic at right describes the Birkeland or field-aligned currents and the ionospheric current systems they connect to.[70]

"A Birkeland current is a set of currents which flow along geomagnetic field lines connecting the Earth’s magnetosphere to the Earth's high latitude ionosphere. Lately, the term Birkeland currents has been expanded by some authors to include magnetic field aligned currents in space plasmas in general. In the Earth’s magnetosphere, the currents are driven by the solar wind and interplanetary magnetic field and by bulk motions of plasma through the magnetosphere (convection which is indirectly driven by the interplanetary environment). The strength of the Birkeland currents changes with activity in the magnetosphere (e.g. during substorms). Small scale variations in the upward current sheets (downward flowing electrons) accelerate magnetospheric electrons which, when they reach the upper atmosphere, create the Aurora Borealis and Australis. In the high latitude ionosphere (or auroral zones), the Birkeland currents close through the region of the auroral electrojet, which flows perpendicular to the local magnetic field in the ionosphere. The Birkeland currents occur in two pairs of field-aligned current sheets. One pair extends from noon through the dusk sector to the midnight sector. The other pair extends from noon through the dawn sector to the midnight sector. The sheet on the high latitude side of the auroral zone is referred to as the Region 1 current sheet and the sheet on the low latitude side is referred to as the Region 2 current sheet."[71]

"A reason why Birkeland currents are particularly interesting is that, in the plasma forced to carry them, they cause a number of plasma physical processes to occur (waves, instabilities, fine structure formation). These in turn lead to consequences such as acceleration of charged particles, both positive and negative, and element separation (such as preferential ejection of oxygen ions). Both of these classes of phenomena should have a general astrophysical interest far beyond that of understanding the space environment of our own Earth."[72]

"Auroral Birkeland currents carry about 100,000 amperes during quiet times[73] and more than 1 million amperes during geomagnetically disturbed times.[74] Birkeland had estimated currents "at heights of several hundred kilometres, and strengths of up to a million amperes" in 1908.[75] The ionospheric currents which connect the field-aligned currents heat up the upper atmosphere due to the finite conductivity of the ionosphere. The heat (also known as Joule heat) is transferred from the ionospheric plasma to the gas of the upper atmosphere which rises and increases drag on low-altitude satellites."[71]

"Birkeland currents can also be created in the laboratory with multi-terawatt pulsed power generators. The resulting cross-section pattern indicates a hollow beam of electrons in the form of a circle of vortices, a formation called the diocotron instability[76] (similar to, but different from, the Kelvin-Helmholtz instability), that subsequently leads to filamentation. Such vortices can be seen in aurora as "auroral curls".[77]"[71]

"Birkeland currents are also one of a class of plasma phenomena called a z-pinch, so named because the azimuthal magnetic fields produced by the current pinches the current into a filamentary cable. This can also twist, producing a helical pinch that spirals like a twisted or braided rope, and this most closely corresponds to a Birkeland current. Pairs of parallel Birkeland currents will also interact due to Ampère's force law: parallel Birkeland currents moving in the same direction will attract each other with an electromagnetic force inversely proportional to their distance apart whilst parallel Birkeland currents moving in opposite directions will repel each other. There is also a short-range circular component to the force between two Birkeland currents that is opposite to the longer-range parallel forces.[78]"[71]

"Electrons moving along a Birkeland current may be accelerated by a plasma double layer. If the resulting electrons approach relativistic velocities (i.e. if they approach the speed of light) they may subsequently produce a Bennett pinch, which in a magnetic field causes the electrons to spiral and emit synchrotron radiation that may include radio, optical (i.e. visible light), x-rays, and gamma rays."[71]

"[T]he currents there are imagined as having come into existence mainly as a secondary effect of the electric corpuscles from the sun drawn in out of space, and thus far come under the second of the possibilities mentioned above". And p. 105, "Fig. 50a represents those in which the current-directions at the storm-centre are directed westwards, and 50b those in which the currents move eastwards"[75].

"[C]urrents there [in the aurora] are imagined as having come into existence mainly as a secondary effect of the electric corpuscles from the sun drawn in out of space,"[75]

## Mars

Main source: Mars

Mars has a magnetosphere.[79]

## 1963-38C

"A U.S. Navy satellite 1963-38C, launched in 1963 and carrying a magnetometer above the ionosphere.[80] ... In 1966 Alfred Zmuda, J.H. Martin, and F.T.Heuring[81] analysed the satellite magnetometer results and reported their findings of magnetic disturbance in the aurora"[71]

## Research

Main source: Research

Hypothesis:

1. Magnetohydrodynamics works even when the electromagnetism occurs first.

## Control groups

Main source: Control groups
This is an image of a Lewis rat. Credit: Charles River Laboratories.

The findings demonstrate a statistically systematic change from the status quo or the control group.

“In the design of experiments, treatments [or special properties or characteristics] are applied to [or observed in] experimental units in the treatment group(s).[82] In comparative experiments, members of the complementary group, the control group, receive either no treatment or a standard treatment.[83]"[84]

## Proof of concept

Main source: Proof of concept

Def. a “short and/or incomplete realization of a certain method or idea to demonstrate its feasibility"[85] is called a proof of concept.

Def. evidence that demonstrates that a concept is possible is called proof of concept.

The proof-of-concept structure consists of

1. background,
2. procedures,
3. findings, and
4. interpretation.[86]

## References

1. "Magnetohydrodynamics, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. October 1, 2012. Retrieved 2012-11-14.
2. "magnetohydrodynamics, In: Wiktionary". San Francisco, California: Wikimedia Foundation, Inc. November 10, 2012. Retrieved 2012-11-14.
3. Y. C. Unruh, S. K. Solanki, & M. Fligge (May 1999). "The spectral dependence of facular contrast and solar irradiance variations". Astronomy and Astrophysics 345 (5): 635-42.
4. "Magnetohydrodynamic turbulence, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. December 25, 2011. Retrieved 2012-11-15.
5. "Magnetic Reynolds number, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. June 22, 2012. Retrieved 2012-11-15.
6. D. Biskamp (2003), Magnetohydrodynamical Turbulence, (Cambridge University Press, Cambridge.)
7. M. K. Verma (2004), Statistical theory of magnetohydrodynamic turbulence, Phys. Rep., 401, 229.
8. P. S. Iroshnikov (1964), Turbulence of a Conducting Fluid in a Strong Magnetic Field, Soviet Astronomy, 7, 566.
9. R. Kraichnan(1965), Inertial-Range Spectrum of Hydromagnetic Turbulence, Physics of Fluids, 8, 1385.
10. M. Dobrowlny, A. Mangeney, P. Veltri (1980), Fully developed anisotropic hydromagnetic turbulence in interplanetary plasma, Phys. Rev. Lett., 45, 144.
11. E. Marsch (1990), Turbulence in the solar wind, in: G. Klare (Ed.), Reviews in Modern Astronomy, Springer, Berlin, p. 43.
12. W. H. Matthaeus, Y. Zhou (1989), Extended inertial range phenomenology of magnetohydrodynamic turbulence, Phys. Fluids B, 1, 1929.
13. S. Galtier, S. V. Nazarenko, A. C. Newell, A. Pouquet (2000), A weak turbulence theory for incompressible magnetohydrodynamics, Journal of Plasma Physics, 63, 447
14. Goldreich, P. & Sridhar, S. (1995), Toward a theory of interstellar turbulence. 2: Strong Alfvénic turbulence, Astrophysical Journal, 438, 763
15. "Stuart number, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 10, 2012. Retrieved 2012-11-15.
16. Journal of Fluid Mechanics 439: 367–94. 2001. doi:10.1017/S0022112001004621.
17. Castellanos, A. (1998). Electrohydrodynamics.
18. "Electrohydrodynamics, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. May 11, 2012. Retrieved 2012-11-14.
19. Kenneth R.Cramer, Shih-l Pai (1973). Magnetofluid Dynamics for Engineers and Applied Physicists. Scripta Publishing Company. ISBN 0-07-013425-1.
20. "Hartmann number, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. June 26, 2012. Retrieved 2012-11-14.
21. "Computational magnetohydrodynamics, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. May 9, 2011. Retrieved 2012-11-14.
22. "Numerical resistivity, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. August 22, 2011. Retrieved 2012-11-14.
23. Axel Brandenburg, Åke Nordlund, Robert F. Stein, and Ulf Torkelsson (June 1995). "Dynamo-generated Turbulence and Large-Scale Magnetic Fields in a Keplerian Shear Flow". The Astrophysical Journal 446 (6): 741-54. doi:10.1086/175831. Retrieved 2012-01-18.
24. Marshallsumter (January 18, 2012). "Radiative dynamo, In: Wikiversity". Retrieved 2012-11-14.
25. DO Gómez, PD Mininni, P Dmitruk (September 27, 2010). "Hall-magnetohydrodynamic small-scale dynamos". Physical Review E 82 (3): 10. doi:10.1103/PhysRevE.82.036406. Retrieved 2012-11-14.
26. "Dynamo theory, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. October 12, 2012. Retrieved 2012-11-17.
27. "Antidynamo theorem, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. October 17, 2012. Retrieved 2012-11-14.
28. "The Magnetic Field of Sunspots" by T.G. Cowling. Monthly Notices of the Royal Astronomical Society, 94: 39-48 (1934)
29. "Plasma stability, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. June 27, 2012. Retrieved 2012-11-15.
30. Andre Gsponer, "Physics of high-intensity high-energy particle beam propagation in open air and outer-space plasmas" (2004)Physics of high-intensity high-energy particle beam propagation in open air and outer-space plasmas
31. Plasma Dictionary
32. "Kink instability, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. February 14, 2011. Retrieved 2012-11-16.
33. "Astrophysical plasma, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 13, 2012. Retrieved 2012-11-16.
34. "Magnetorotational instability, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 9, 2011. Retrieved 2012-11-14.
35. Velikhov, E. P. (1959). "Stability of an Ideally Conducting Liquid Flowing Between Cylinders Rotating in a Magnetic Field". J. Exptl. Theoret. Phys. 36: 1398–1404.
36. Chandrasekhar, S. (1960). "The stability of non-dissipative Couette flow in hydromagnetics". Proc. Natl. Acad. Sci. 46 (2): 253–257. doi:10.1073/pnas.46.2.253.
37. Acheson, D. J.; Hide, R. (1973). "Hydromagnetics of Rotating Fluids". Reports on Progress in Physics 36 (2): 159–221. doi:10.1088/0034-4885/36/2/002.
38. Balbus, Steven A.; Hawley, John F. (1991). "A powerful local shear instability in weakly magnetized disks. I - Linear analysis. II - Nonlinear evolution". Astrophysical Journal 376: 214–233. doi:10.1086/170270.
39. Rüdiger, Günther; Hollerbach, Rainer (2004). The Magnetic Universe: Geophysical and Astrophysical Dynamo Theory. Wiley-VCH. ISBN 3-527-40409-0.
40. "Resistive skin time, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. August 27, 2012. Retrieved 2012-11-14.
41. "Coronal seismology, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. September 13, 2012. Retrieved 2012-11-15.
42. "Waves in plasmas, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 8, 2012. Retrieved 2012-11-18.
43. "Shocks and discontinuities (magnetohydrodynamics), In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 24, 2011. Retrieved 2012-11-16.
44. "Magnetic tension force, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. September 14, 2012. Retrieved 2012-11-16.
45. Paul M. Bellan (2000). Spheromaks: A Practical Application of Magnetohydrodynamic dynamos and plasma self-organization. pp. 71–79. ISBN 1-86094-141-9.
46. "Taylor state, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. June 1, 2012. Retrieved 2012-11-18.
47. "Solar dynamo, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. January 10, 2012. Retrieved 2012-11-15.
48. Markus J. Aschwanden (2007). Erdelyi R. ed. "Fundamental Physical Processes in Coronae: Waves, Turbulence, Reconnection, and Particle Acceleration In: Waves & Oscillations in the Solar Atmosphere: Heating and Magneto-Seismology". Proceedings IAU Symposium 3 (S247): 257–68. doi:10.1017/S1743921308014956.
49. C. T. Russell, D. N. Baker and J. A. Slavin (January 1, 1988). Faith Vilas, Clark R. Chapman, Mildred Shapley Matthews. ed. The Magnetosphere of Mercury, In: Mercury. Tucson, Arizona, United States of America: University of Arizona Press. pp. 514-61. ISBN 0816510857. Bibcode: 1988merc.book..514R. Retrieved 2012-08-23.
50. "Magnetosphere, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 11, 2012. Retrieved 2012-11-14.
51. Tielong Zhang and Håkan Svedhem (April 05, 2012). "A magnetic surprise for Venus Express". The Netherlands: European Space Agency. Retrieved 2012-04-07.
52. Michael A. Seeds (2004). Astronomy: The Solar System and Beyond (4th ed.). Brooks Cole. ISBN 0-534-42111-3.
53. David R. Williams (January 6, 2005). "Planetary Fact Sheets". NASA National Space Science Data Center. Retrieved 2006-08-10.
54. J. Kelly Beatty, Carolyn Collins Petersen, Andrew Chaikin (1999). The New Solar System. Cambridge University Press. ISBN 0-521-64587-5.
55. Staff (January 30, 2008). "Mercury's Internal Magnetic Field". NASA. Retrieved 2008-04-07.
56. "Mercury (planet), In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. April 3, 2013. Retrieved 2013-04-03.
57. Bill Steigerwald (June 2, 2009). "Magnetic Tornadoes Could Liberate Mercury's Tenuous Atmosphere". NASA Goddard Space Flight Center. Retrieved 2009-07-18.
58. Dolginov, Nature of the Magnetic Field in the Neighborhood of Venus, COsmic Research, 1969
59. Kivelson G. M., Russell, C. T. (1995). Introduction to Space Physics. Cambridge University Press. ISBN 0-521-45714-9.
60. Upadhyay, H. O.; Singh, R. N. (April 1995). "Cosmic ray Ionization of Lower Venus Atmosphere". Advances in Space Research 15 (4): 99–108. doi:10.1016/0273-1177(94)00070-H.
61. "Venus, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. September 27, 2012. Retrieved 2012-10-04.
62. Hannes Alfvén (1976). Evolution of the Solar System. Washington. D.C., USA: Scientific and Technical Information Office, National Aeronautics and Space Administration.
63. NOVA | Magnetic Storm | What Drives Earth's Magnetic Field? | PBS
64. Earth's Inconstant Magnetic Field - NASA Science
65. An introduction to the ionosphere ... - Google Books
66. On The Atmospheric Dynamo Theory
67. NASA - Sounding Rockets Study How Winds In Space Drive Currents in the Upper Atmosphere
68. "Atmospheric dynamo, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. January 11, 2012. Retrieved 2012-11-15.
69. "Magnetite, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. November 16, 2012. Retrieved 2012-11-16.
70. G. Le, J. A. Slavin, and R. J. Strangeway (2010). "Space Technology 5 observations of the imbalance of regions 1 and 2 field-aligned currents and its implication to the cross-polar cap Pedersen currents". Journal of Geophys. Research 115 (A07202). doi:10.1029/2009JA014979.
71. "Birkeland current, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. October 15, 2012. Retrieved 2012-11-17.
72. Carl-Gunne Fälthammar (December 1986). "Magnetosphere-Ionosphere Interactions. Near Earth Manifestations of the Plasma Universe". IEEE Transactions on Plasma Science PS-14 (6): 616–28. doi:10.1109/TPS.1986.4316613.
73. Akira Suzuki, Naoshi Fukushima (1998). "Space current around the earth obtained with Ampère’s law applied to the MAGSAT orbit and data". Earth Planets Space 50 (1): 43–56.
74. B. J. Anderson, J. B. Gary, T. A. Potemra, R. A. Frahm, J. R. Sharber, and J. D. Winningham (1998). "UARS observations of Birkeland currents and Joule heating rates for the November 4, 1993, storm". J. Geophys. Research 103 (A11): 26323–35. doi:10.1029/98JA01236.
75. Birkeland, Kristian (1908 (section 1), 1913 (section 2)). The Norwegian Aurora Polaris Expedition 1902-1903. New York and Christiania (now Oslo): H. Aschehoug & Co.  out-of-print, full text online
76. Plasma phenomena - instabilities
77. Pseudo-color, white-light images of curl formations in auroral arcs
78. Electromagnetic Forces
79. Dolginov, S. S.; Yeroshenko, Y. G.; Zhuzgov, L. N.; Sharova, V. A.; Gringauz, K. I.; Bezrukikh, V. V.; Breus, T. K.; Verigin, M. I.; Remizov, A. P. (1976). Magnetic field and plasma inside and outside of the Martian magnetosphere, In: Solar-Wind Interaction with the Planets Mercury, Venus, and Mars. Greenbelt, Maryland USA: NASA Goddard Space Flight Center. pp. 1-20. Retrieved 2012-04-07.
80. Peratt, A. L., Peter, W., & Snell, C. M. (June 19–23, 1989). 3-dimensional particle-in-cell simulations of spiral galaxies. Heidelberg, Federal Republic of Germany: Kluwer Academic Publishers. pp. 143-150. Bibcode: 1990IAUS..140..143P. Retrieved 16 May 2012.
81. Alfred Zmuda, J.H. Martin and F.T.Heuring (1966). "Transverse Magnetic Disturbances at 1100 Kilometers in the Auroral Region". Journal of Geophysical Research 71 (21): 5033–45. doi:10.1029/JZ071i021p05033.
82. Klaus Hinkelmann, Oscar Kempthorne (2008). Design and Analysis of Experiments, Volume I: Introduction to Experimental Design (2nd ed.). Wiley. ISBN 978-0-471-72756-9.
83. R. A. Bailey (2008). Design of comparative experiments. Cambridge University Press. ISBN 978-0-521-68357-9.
84. "Treatment and control groups, In: Wikipedia". San Francisco, California: Wikimedia Foundation, Inc. May 18, 2012. Retrieved 2012-05-31.
85. "proof of concept, In: Wiktionary". San Francisco, California: Wikimedia Foundation, Inc. November 10, 2012. Retrieved 2013-01-13.
86. Ginger Lehrman and Ian B Hogue, Sarah Palmer, Cheryl Jennings, Celsa A Spina, Ann Wiegand, Alan L Landay, Robert W Coombs, Douglas D Richman, John W Mellors, John M Coffin, Ronald J Bosch, David M Margolis (August 13, 2005). "Depletion of latent HIV-1 infection in vivo: a proof-of-concept study". Lancet 366 (9485): 549-55. doi:10.1016/S0140-6736(05)67098-5. Retrieved 2012-05-09.