Physics:Magnetic diffusion

Magnetic diffusion refers to the motion of magnetic fields, typically in the presence of a conducting solid or fluid such as a plasma. The motion of magnetic fields is described by the magnetic diffusion equation and is due primarily to induction and diffusion of magnetic fields through the material. The magnetic diffusion equation is a partial differential equation commonly used in physics. Understanding the phenomenon is essential to magnetohydrodynamics and has important consequences in astrophysics, geophysics, and electrical engineering.

Equation

The magnetic diffusion equation is [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} = \nabla \times \left[\vec{v} \times \vec{B}\right] + \frac{1}{\mu_0 \sigma}\nabla^2 \vec{B} }[/math] where [math]\displaystyle{ \mu_0 }[/math] is the permeability of free space and [math]\displaystyle{ \sigma }[/math] is the electrical conductivity of the material, which is assumed to be constant. [math]\displaystyle{ \vec{v} }[/math] denotes the (non-relativistic) velocity of the plasma. The first term on the right hand side accounts for effects from induction of the plasma, while the second accounts for diffusion. The latter acts as a dissipation term, resulting in a loss of magnetic field energy to heat. The relative importance of the two terms is characterized by the magnetic Reynolds number, [math]\displaystyle{ R_m }[/math].

In the case of a non-uniform conductivity the magnetic diffusion equation is [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} = \nabla \times \left[\vec{v} \times \vec{B}\right] - \frac{1}{\mu_0} \nabla \times \left[\frac{1}{\sigma} \nabla \times \vec{B} \right] }[/math] however, it becomes significantly harder to solve.

Derivation

Starting from the generalized Ohm's law: [math]\displaystyle{ \vec{J} = \sigma \left(\vec{E}+\vec{v}\times\vec{B} \right) }[/math] and the curl equations for small displacement currents (i.e. low frequencies) [math]\displaystyle{ \nabla\times\vec{B} = \mu_0 \vec{J} + \epsilon_0 \mu_0 \frac{\partial \vec{E}}{\partial t} \approx \mu_0 \vec{J} }[/math] [math]\displaystyle{ \nabla\times\vec{E} = -\frac{\partial \vec{B}}{\partial t} }[/math] substitute [math]\displaystyle{ \vec{J} }[/math] into the Ampere-Maxwell law to get [math]\displaystyle{ \frac{1}{\mu_0 \sigma} \nabla\times\vec{B} = \vec{E} + \vec{v}\times\vec{B} \quad\Rightarrow\quad \vec{E} = \frac{1}{\mu_0 \sigma}\nabla\times\vec{B}-\vec{v}\times\vec{B}. }[/math] Taking the curl of the above equation and substituting into Faraday's law, [math]\displaystyle{ \nabla\times\vec{E} = \nabla\times\left(\frac{1}{\mu_0 \sigma}\nabla\times\vec{B} - \vec{v}\times\vec{B}\right) = -\frac{\partial \vec{B}}{\partial t}. }[/math] This expression can be simplified further by writing it in terms of the i-th component of [math]\displaystyle{ \vec{B} }[/math] and the Levi-Cevita tensor [math]\displaystyle{ \varepsilon_{ijk} }[/math]: [math]\displaystyle{ \begin{align} -\frac{\partial B_i}{\partial t} & = \varepsilon_{ijk} \partial_j \left( \frac{1}{\mu_0 \sigma}\varepsilon_{klm}\partial_l B_m - \varepsilon_{klm}v_l B_m \right)\\ & = \varepsilon_{kij} \varepsilon_{klm} \left(\frac{1}{\mu_0 \sigma}\partial_j\partial_l B_m - \left(v_l \partial_j B_m + B_m \partial_j v_l \right)\right) \end{align} }[/math] Using the identity [math]\displaystyle{ \varepsilon_{kij} \varepsilon_{klm}= \delta_{il}\delta_{jm}-\delta_{im}\delta_{jl} }[/math] and recalling [math]\displaystyle{ \partial_j B_j = 0 }[/math], the cross products can be eliminated: [math]\displaystyle{ \begin{align} -\frac{\partial B_i}{\partial t} & = \frac{1}{\mu_0 \sigma}\left(\partial_i\partial_j B_j - \partial_j \partial_j B_i\right) - \left(v_i \partial_j B_j - v_j \partial_j B_i\right) - \left(B_j \partial_j v_i - B_i \partial_j v_j\right) \\ & = -\frac{1}{\mu_0 \sigma}\partial_j \partial_j B_i + v_j \partial_j B_i - \left(B_j \partial_j v_i - B_i \partial_j v_j\right) \end{align} }[/math] Written in vector form, the final expression is [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t}+\left(\vec{v}\cdot\nabla\right)\vec{B} = \frac{D\vec{B}}{Dt} = \left(\vec{B}\cdot\nabla\right)\vec{v}-\vec{B}\left(\nabla\cdot\vec{v}\right)+\frac{1}{\mu_0 \sigma}\nabla^2 \vec{B} }[/math] where [math]\displaystyle{ \frac{D}{Dt}=\frac{\partial}{\partial t}+\vec{v}\cdot\nabla }[/math] is the material derivative. This can be rearranged into a more useful form using vector calculus identities and [math]\displaystyle{ \nabla \cdot \vec{B}=0 }[/math]: [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t}= \nabla \times [\vec{v} \times \vec{B}] + \frac{1}{\mu_0 \sigma}\nabla^2 \vec{B} }[/math] In the case [math]\displaystyle{ \vec{v}=0 }[/math], this becomes a diffusion equation for the magnetic field, [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} = \frac{1}{\mu_0 \sigma}\nabla^2 \vec{B} = \eta\nabla^2 \vec{B} }[/math] where [math]\displaystyle{ \eta = \frac{1}{\mu_0 \sigma} }[/math] is the magnetic diffusivity.

Limiting Cases

In some cases it is possible to neglect one of the terms in the magnetic diffusion equation. This is done by estimating the magnetic Reynolds number [math]\displaystyle{ R_m = \frac{v L}{\eta} }[/math] where [math]\displaystyle{ \eta }[/math] is the diffusivity, [math]\displaystyle{ v }[/math] is the magnitude of the plasma's velocity and [math]\displaystyle{ L }[/math] is a characteristic length of the plasma.

[math]\displaystyle{ (R_m) }[/math]	Physical Condition	Dominating Term	Magnetic Diffusion Equation	Examples
[math]\displaystyle{ \gg 1 }[/math]	Large electrical conductivity, large length scales or high plasma velocity.	The inductive term dominates in this case. The motion of magnetic fields is determined by the flow of the plasma. This is the case for most naturally occurring plasmas in the universe.	[math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} \approx \nabla \times [\vec{v} \times \vec{B}] }[/math]	The Sun [math]\displaystyle{ (R_m \approx 10^6) }[/math] or the core of the earth [math]\displaystyle{ (R_m \approx 10^3) }[/math]
[math]\displaystyle{ \ll 1 }[/math]	Small electrical conductivity, small length scales or low plasma velocity.	The diffusive term dominates in this case. The motion of the magnetic field obeys the typical (nonconducting) fluid diffusion equation.	[math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} \approx \frac{1}{\mu_0 \sigma}\nabla^2 \vec{B} }[/math]	Solar flares or created in laboratories using mercury or other liquid metals.

Relation to Skin Effect

At low frequencies, the skin depth [math]\displaystyle{ \delta }[/math] for the penetration of an AC electromagnetic field into a conductor is: [math]\displaystyle{ \delta = \sqrt{\frac{2}{\mu \sigma \omega}} }[/math] Comparing with the formula for [math]\displaystyle{ \eta }[/math], the skin depth is the diffusion length of the field over one period of oscillation: [math]\displaystyle{ \delta = \sqrt{\frac{2\eta}{\omega}} = \sqrt{\frac{\eta T}{\pi}} }[/math]

Examples and Visualization

Example of magnetic field frozen into fluid flow.

For the limit [math]\displaystyle{ R_m \gg 1 }[/math], the magnetic field lines become "frozen in" to the motion of the conducting fluid. A simple example illustrating this behavior has a sinusoidally-varying shear flow [math]\displaystyle{ \vec{v} = v_0\sin(k y)\hat{x} }[/math] with a uniform initial magnetic field [math]\displaystyle{ \vec{B}\left(\vec{r},0\right) = B_0\hat{y} }[/math]. The equation for this limit, [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} = \nabla \times [\vec{v} \times \vec{B}] }[/math], has the solution [math]\displaystyle{ \vec{B}\left(\vec{r},t\right) = B_0 k v_0 t\cos(k y)\hat{x}+B_0\hat{y} }[/math] As can be seen in the figure to the right, the fluid drags the magnetic field lines so that they obtain the sinusoidal character of the flow field.

For the limit [math]\displaystyle{ R_m \ll 1 }[/math], the magnetic diffusion equation [math]\displaystyle{ \frac{\partial \vec{B}}{\partial t} = \frac{1}{\mu_0 \sigma} \nabla^2 \vec{B} }[/math] is just a vector-valued form of the heat equation. For a localized initial magnetic field (e.g. Gaussian distribution) within a conducting material, the maxima and minima will asymptotically decay to a value consistent with Laplace's equation for the given boundary conditions. This behavior is illustrated in the figure below.

Diffusion Times for Stationary Conductors

For stationary conductors [math]\displaystyle{ (R_m=0) }[/math] with simple geometries a time constant called magnetic diffusion time can be derived. Different one-dimensional equations apply for conducting slabs and conducting cylinders with constant magnetic permeability. Also, different diffusion time equations can be derived for nonlinear saturable materials such as steel.