Alfvén's theorem

In ideal magnetohydrodynamics, Alfvén's theorem, or the frozen-in flux theorem, states that electrically conducting fluids and embedded magnetic fields are constrained to move together in the limit of large magnetic Reynolds numbers. It is named after Hannes Alfvén, who put the idea forward in 1943.

Alfvén's theorem implies that the magnetic topology of a fluid in the limit of a large magnetic Reynolds number cannot change. This approximation breaks down in current sheets, where magnetic reconnection can occur.

History
The concept of magnetic fields being frozen into fluids with infinite electrical conductivity was first proposed by Hannes Alfvén in a 1943 paper titled "On the Existence of Electromagnetic-Hydrodynamic Waves", published in the journal Arkiv för matematik, astronomi och fysik. He wrote:

"On the Existence of Electromagnetic-Hydrodynamic Waves" interpreted the results of Alfvén's earlier paper "Existence of Electromagnetic-Hydrodynamic Waves", published in the journal Nature in 1942.

Later in life, Alfvén advised against the use of his own theorem.

Overview
Informally, Alfvén's theorem refers to the fundamental result in ideal magnetohydrodynamic theory that electrically conducting fluids and the magnetic fields within are constrained to move together in the limit of large magnetic Reynolds numbers (Rm)—such as when the fluid is a perfect conductor or when velocity and length scales are infinitely large. Motions of the two are constrained in that all bulk fluid motions perpendicular to the magnetic field result in matching perpendicular motion of the field at the same velocity and vice versa.

Formally, the connection between the movement of the fluid and the movement of the magnetic field is detailed in two primary results, often referred to as magnetic flux conservation and magnetic field line conservation. Magnetic flux conservation implies that the magnetic flux through a surface moving with the bulk fluid velocity is constant, and magnetic field line conservation implies that, if two fluid elements are connected by a magnetic field line, they will always be.

Flux tubes and field lines
Alfvén's theorem is frequently expressed in terms of magnetic flux tubes and magnetic field lines.

A magnetic flux tube is a tube- or cylinder-like region of space containing a magnetic field such that its sides are everywhere parallel to the field. Consequently, the magnetic flux through these sides is zero, and the cross-sections along the tube's length have constant, equal magnetic flux. In the limit of a large magnetic Reynolds number, Alfvén's theorem requires that these surfaces of constant flux move with the fluid that they are embedded in. As such, magnetic flux tubes are frozen into the fluid.

The intersection of the sides of two magnetic flux tubes form a magnetic field line, a curve that is everywhere parallel to the magnetic field. In fluids where flux tubes are frozen-in, it then follows that magnetic field lines must also be frozen-in. However, the conditions for frozen-in field lines are weaker than the conditions for frozen-in flux tubes, or, equivalently, for conservation of flux.

Mathematical statement
In mathematical terms, Alfvén's theorem states that, in an electrically conducting fluid in the limit of a large magnetic Reynolds number, the magnetic flux $S_{1}$ through an orientable, open material surface advected by a macroscopic, space- and time-dependent velocity field $S_{2}$ is constant, or
 * $$\frac{D\Phi_B}{Dt} = 0 ,$$

where $S_{1}$ is the advective derivative.

Flux conservation
In ideal magnetohydrodynamics, magnetic induction dominates over magnetic diffusion at the velocity and length scales being studied. The diffusion term in the governing induction equation is then assumed to be small relative to the induction term and is neglected. The induction equation then reduces to its ideal form:
 * $$\frac{\partial\mathbf{B}}{\partial t} = \nabla \times \left(\mathbf{v}\times\mathbf{B}\right).$$

The conservation of magnetic flux through material surfaces embedded in the fluid follows directly from the ideal induction equation and the assumption of no magnetic monopoles through Gauss's law for magnetism.

In an electrically conducting fluid with a space- and time-dependent magnetic field $S_{2}$ and velocity field $&Phi;_{B}$, an arbitrary, orientable, open surface $v$ at time $t$ is advected by $v$ in a small time $&delta;t$ to the surface $D/Dt = &part;/&part;t + (v &sdot; &nabla;)$. The rate of change of the magnetic flux through the surface as it is advected from $S_{1}$ to $S_{2}$ is then
 * $$\frac{D\Phi_B}{Dt} = \lim_{\delta t \to 0} \frac{\iint_{S_2} \mathbf{B}(t+\delta t) \cdot d\mathbf{S}_2 - \iint_{S_1} \mathbf{B}(t) \cdot d\mathbf{S}_1}{\delta t}.$$

The surface integral over $S_{3}$ can be re expressed by applying Gauss's law for magnetism to assume that the magnetic flux through a closed surface formed by $B$, $v$, and the surface $S_{1}$ that connects the boundaries of $v$ and $S_{2}$ is zero. At time $S_{1}$, this relationship can be expressed as
 * $$0 = -\iint_{S_1} \mathbf{B}(t+\delta t)\cdot d\mathbf{S}_1 + \iint_{S_2} \mathbf{B}(t+\delta t)\cdot d\mathbf{S}_2 + \iint_{S_3} \mathbf{B}(t+\delta t)\cdot d\mathbf{S}_3 ,$$

where the sense of $S_{2}$ was reversed so that $S_{2}$ points outwards from the enclosed volume. In the surface integral over $S_{1}$, the differential surface element $S_{2}$ where $S_{3}$ is the line element around the boundary $S_{1}$ of the surface $S_{2}$. Solving for the surface integral over $t + &delta;t$ then gives
 * $$\iint_{S_2} \mathbf{B}(t+\delta t)\cdot d\mathbf{S}_2 = \iint_{S_1} \mathbf{B}(t+\delta t)\cdot d\mathbf{S}_1 - \oint_{\partial S_1} \left(\mathbf{v}\ \delta t \times \mathbf{B}(t)\right) \cdot d\mathbf{l},$$

where the final term was rewritten using the properties of scalar triple products and a first-order approximation was taken. Substituting this into the expression for $S_{1}$ and simplifying results in
 * $$\begin{align}

\frac{D\Phi_B}{Dt} = \lim_{\delta t \to 0} \iint_{S_1} \frac{\mathbf{B}(t+\delta t) - \mathbf{B}(t)}{\delta t} \cdot d\mathbf{S}_1 - \oint_{\partial S_1} \left(\mathbf{v} \times \mathbf{B}(t)\right) \cdot d\mathbf{l}. \end{align}$$ Applying the definition of a partial derivative to the integrand of the first term, applying Stokes' theorem to the second term, and combining the resultant surface integrals gives
 * $$\frac{D\Phi_B}{Dt} = \iint_{S_1} \left( \frac{\partial\mathbf{B}}{\partial t} - \nabla\times\left(\mathbf{v} \times \mathbf{B}\right) \right) \cdot d\mathbf{S}_1.$$

Using the ideal induction equation, the integrand vanishes, and
 * $$\frac{D\Phi_B}{Dt} = 0.$$

Field line conservation
Field line conservation can also be derived mathematically using the ideal induction equation, Gauss's law for magnetism, and the mass continuity equation.

The ideal induction equation can be rewritten using a vector identity and Gauss's law for magnetism as
 * $$\frac{\partial\mathbf{B}}{\partial t} = (\mathbf{B}\cdot\nabla)\mathbf{v} - (\mathbf{v}\cdot\nabla)\mathbf{B} - \mathbf{B}(\nabla\cdot\mathbf{v}).$$

Using the mass continuity equation,
 * $$\frac{\partial\rho}{\partial t} + (\mathbf{v}\cdot\nabla)\rho = -\rho\nabla\cdot\mathbf{v},$$

the ideal induction equation can be further rearranged to give
 * $$\frac{D}{Dt} \left(\frac{\mathbf{B}}{\rho}\right) = \left(\frac{\mathbf{B}}{\rho}\cdot\nabla\right)\mathbf{v}.$$

Similarly, for a line segment $dS_{1}$ where $S_{3}$ is the bulk plasma velocity at one end and $dS_{3} = dl &times; v &delta;t$ is the velocity at the other end, the differential velocity between the two ends is $dl$ and
 * $$\frac{D\delta\mathbf{l}}{Dt} = (\delta\mathbf{l}\cdot\nabla)\mathbf{v}$$,

which has the same form as the equation obtained previously for $&part;S_{1}$. Therefore, if $S_{1}$ and $S_{2}$ are initially parallel, they will remain parallel.

While flux conservation implies field line conservation (see ), the conditions for the latter are weaker than the conditions for the former. Unlike the conditions for flux conservation, the conditions for field line conservation can be satisfied when an additional, source term parallel to the magnetic field is present in the ideal induction equation.

Mathematically, for field lines to be frozen-in, the fluid must satisfy
 * $$\left( \frac{\partial\mathbf{B}}{\partial t} - \nabla \times \left(\mathbf{v}\times\mathbf{B}\right)\right) \times \mathbf{B} = 0,$$

whereas, for flux to be conserved, the fluid must satisfy the stronger condition imposed by the ideal induction equation.

Kelvin's circulation theorem
Kelvin's circulation theorem states that vortex tubes moving with an ideal fluid are frozen to the fluid, analogous to how magnetic flux tubes moving with a perfectly conducting ideal-MHD fluid are frozen to the fluid. The ideal induction equation takes the same form as the equation for vorticity $D&Phi;_{B}/Dt$ in an ideal fluid where $&delta;l$ is the velocity field:
 * $$\frac{\partial \boldsymbol{\omega}}{\partial t} = \nabla \times (\mathbf{v}\times\boldsymbol{\omega}).$$

However, the induction equation is linear, whereas there is a nonlinear relationship between $v$ and $v + &delta;v$ in the vorticity equation.

Implications
Alfvén's theorem indicates that the magnetic field topology cannot change in a perfectly conducting fluid. However, in the case of complicated or turbulent flows, this would lead to highly tangled magnetic fields with very complicated topologies that should impede the fluid motions. Astrophysical plasmas with high electrical conductivities do not generally show such complicated tangled fields. Magnetic reconnection seems to occur in these plasmas unlike what would be expected from the flux freezing conditions. This has important implications for magnetic dynamos. In fact, a very high electrical conductivity translates into high magnetic Reynolds numbers, which indicates that the plasma will be turbulent.

Resistive fluids
Even for the non-ideal case, in which the electric conductivity is not infinite, a similar result can be obtained by defining the magnetic flux transporting velocity by writing:

\nabla \times (\bf{w}\times \bf{B})=\eta \nabla^2 \bf{B} + \nabla \times (\bf{v} \times \bf{B}), $$ in which, instead of fluid velocity $&delta;v = (&delta;l &sdot; &nabla;)v$, the flux velocity $B/&rho;$ has been used. Although, in some cases, this velocity field can be found using magnetohydrodynamic equations, the existence and uniqueness of this vector field depends on the underlying conditions.

Stochastic flux freezing
The views on flux freezing in highly conducting plasmas are inconsistent with the phenomenon of spontaneous stochasticity. In some textbooks, however, magnetic flux freezing should hold increasingly better as magnetic diffusivity tends to zero (non-dissipative regime). But the subtlety is that very large magnetic Reynolds numbers (i.e., small electric resistivity or high electrical conductivities) are usually associated with high kinetic Reynolds numbers (i.e., very small viscosities). If kinematic viscosity tends to zero simultaneously with the resistivity, and if the plasma becomes turbulent (associated with high Reynolds numbers), then Lagrangian trajectories will no longer be unique. The flux freezing argument, discussed above, does not apply in general, and stochastic flux freezing must be employed.

The stochastic flux-freezing theorem for resistive magnetohydrodynamics generalizes ordinary flux-freezing discussed above. This generalized theorem states that magnetic field lines of the fine-grained magnetic field $B$ are "frozen-in" to the stochastic trajectories solving the following stochastic differential equation, known as the Langevin equation:
 * $$ d{\bf{x}}={\bf{u}}({\bf{x}},t)dt+\sqrt{2\eta} d{\bf{W}}(t)$$

in which $&eta;$ is magnetic diffusivity and $W$ is the three-dimensional Gaussian white noise (see also Wiener process.) The many "virtual" field-vectors that arrive at the same final point must be averaged to obtain the physical magnetic field at that point.