Magnetic diffusion

Last updated April 11, 2024

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 (also referred to as the induction equation ) is

{\frac {\partial {\vec {B}}}{\partial t}}=\nabla \times \left[{\vec {v}}\times {\vec {B}}\right]+{\frac {1}{\mu _{0}\sigma }}\nabla ^{2}{\vec {B}}

where $\mu _{0}$ is the permeability of free space and $\sigma$ is the electrical conductivity of the material, which is assumed to be constant. ${\vec {v}}$ 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, $R_{m}$ .

In the case of a non-uniform conductivity the magnetic diffusion equation is

{\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]

however, it becomes significantly harder to solve.

Derivation

Starting from the generalized Ohm's law:^[1]^[2]

{\vec {J}}=\sigma \left({\vec {E}}+{\vec {v}}\times {\vec {B}}\right)

and the curl equations for small displacement currents (i.e. low frequencies)

\nabla \times {\vec {B}}=\mu _{0}{\vec {J}}+\epsilon _{0}\mu _{0}{\frac {\partial {\vec {E}}}{\partial t}}\approx \mu _{0}{\vec {J}}

\nabla \times {\vec {E}}=-{\frac {\partial {\vec {B}}}{\partial t}}

substitute ${\vec {J}}$ into the Ampere-Maxwell law to get

{\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}}.

Taking the curl of the above equation and substituting into Faraday's law,

\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}}.

This expression can be simplified further by writing it in terms of the i-th component of ${\vec {B}}$ and the Levi-Cevita tensor $\varepsilon _{ijk}$ :

{\begin{aligned}-{\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{aligned}}

Using the identity^[3] $\varepsilon _{kij}\varepsilon _{klm}=\delta _{il}\delta _{jm}-\delta _{im}\delta _{jl}$ and recalling $\partial _{j}B_{j}=0$ , the cross products can be eliminated:

{\begin{aligned}-{\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{aligned}}

Written in vector form, the final expression is

{\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}}

where ${\frac {D}{Dt}}={\frac {\partial }{\partial t}}+{\vec {v}}\cdot \nabla$ is the material derivative. This can be rearranged into a more useful form using vector calculus identities and $\nabla \cdot {\vec {B}}=0$ :

{\frac {\partial {\vec {B}}}{\partial t}}=\nabla \times [{\vec {v}}\times {\vec {B}}]+{\frac {1}{\mu _{0}\sigma }}\nabla ^{2}{\vec {B}}

In the case ${\vec {v}}=0$ , this becomes a diffusion equation for the magnetic field,

{\frac {\partial {\vec {B}}}{\partial t}}={\frac {1}{\mu _{0}\sigma }}\nabla ^{2}{\vec {B}}=\eta \nabla ^{2}{\vec {B}}

where $\eta ={\frac {1}{\mu _{0}\sigma }}$ 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

R_{m}={\frac {vL}{\eta }}

where $\eta$ is the diffusivity, $v$ is the magnitude of the plasma's velocity and $L$ is a characteristic length of the plasma.

$(R_{m})$	Physical Condition	Dominating Term	Magnetic Diffusion Equation	Examples
$\gg 1$	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.	${\frac {\partial {\vec {B}}}{\partial t}}\approx \nabla \times [{\vec {v}}\times {\vec {B}}]$	The Sun $(R_{m}\approx 10^{6})$ or the core of the earth $(R_{m}\approx 10^{3})$
$\ll 1$	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.	${\frac {\partial {\vec {B}}}{\partial t}}\approx {\frac {1}{\mu _{0}\sigma }}\nabla ^{2}{\vec {B}}$	Solar flares or created in laboratories using mercury or other liquid metals.

Relation to Skin Effect

At low frequencies, the skin depth $\delta$ for the penetration of an AC electromagnetic field into a conductor is:

\delta ={\sqrt {\frac {2}{\mu \sigma \omega }}}

Comparing with the formula for $\eta$ , the skin depth is the diffusion length of the field over one period of oscillation:

\delta ={\sqrt {\frac {2\eta }{\omega }}}={\sqrt {\frac {\eta T}{\pi }}}

Examples and Visualization

Example of magnetic field frozen into fluid flow. Magnetic shear flow.gif — Example of magnetic field frozen into fluid flow.

For the limit $R_{m}\gg 1$ , 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

{\vec {v}}=v_{0}\sin(ky){\hat {x}}

with a uniform initial magnetic field ${\vec {B}}\left({\vec {r}},0\right)=B_{0}{\hat {y}}$ . The equation for this limit, ${\frac {\partial {\vec {B}}}{\partial t}}=\nabla \times [{\vec {v}}\times {\vec {B}}]$ , has the solution^[4]

{\vec {B}}\left({\vec {r}},t\right)=B_{0}kv_{0}t\cos(ky){\hat {x}}+B_{0}{\hat {y}}

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 $R_{m}\ll 1$ , the magnetic diffusion equation ${\frac {\partial {\vec {B}}}{\partial t}}={\frac {1}{\mu _{0}\sigma }}\nabla ^{2}{\vec {B}}$ 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 $(R_{m}=0)$ with simple geometries a time constant called magnetic diffusion time can be derived.^[5] 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.

Related Research Articles

Maxwell's equations, or Maxwell–Heaviside equations, are a set of coupled partial differential equations that, together with the Lorentz force law, form the foundation of classical electromagnetism, classical optics, electric and magnetic circuits. The equations provide a mathematical model for electric, optical, and radio technologies, such as power generation, electric motors, wireless communication, lenses, radar, etc. They describe how electric and magnetic fields are generated by charges, currents, and changes of the fields. The equations are named after the physicist and mathematician James Clerk Maxwell, who, in 1861 and 1862, published an early form of the equations that included the Lorentz force law. Maxwell first used the equations to propose that light is an electromagnetic phenomenon. The modern form of the equations in their most common formulation is credited to Oliver Heaviside.

Magnetohydrodynamics is a model of electrically conducting fluids that treats all interpenetrating particle species together as a single continuous medium. It is primarily concerned with the low-frequency, large-scale, magnetic behavior in plasmas and liquid metals and has applications in numerous fields including geophysics, astrophysics, and engineering.

Noether's theorem states that every continuous symmetry of the action of a physical system with conservative forces has a corresponding conservation law. This is the first of two theorems proven by mathematician Emmy Noether in 1915 and published in 1918. The action of a physical system is the integral over time of a Lagrangian function, from which the system's behavior can be determined by the principle of least action. This theorem only applies to continuous and smooth symmetries of physical space.

In physics, a Langevin equation is a stochastic differential equation describing how a system evolves when subjected to a combination of deterministic and fluctuating ("random") forces. The dependent variables in a Langevin equation typically are collective (macroscopic) variables changing only slowly in comparison to the other (microscopic) variables of the system. The fast (microscopic) variables are responsible for the stochastic nature of the Langevin equation. One application is to Brownian motion, which models the fluctuating motion of a small particle in a fluid.

<span class="mw-page-title-main">Fokker–Planck equation</span> Partial differential equation

In statistical mechanics and information theory, the Fokker–Planck equation is a partial differential equation that describes the time evolution of the probability density function of the velocity of a particle under the influence of drag forces and random forces, as in Brownian motion. The equation can be generalized to other observables as well. The Fokker-Planck equation has multiple applications in information theory, graph theory, data science, finance, economics etc.

Linear elasticity is a mathematical model of how solid objects deform and become internally stressed due to prescribed loading conditions. It is a simplification of the more general nonlinear theory of elasticity and a branch of continuum mechanics.

In electromagnetism, displacement current density is the quantity $\partial D /\partial t$ appearing in Maxwell's equations that is defined in terms of the rate of change of $D$ , the electric displacement field. Displacement current density has the same units as electric current density, and it is a source of the magnetic field just as actual current is. However it is not an electric current of moving charges, but a time-varying electric field. In physical materials, there is also a contribution from the slight motion of charges bound in atoms, called dielectric polarization.

In differential geometry, the four-gradient $is the four-vector analogue of the gradient from vector calculus.$

In electromagnetism, the Lorenz gauge condition or Lorenz gauge is a partial gauge fixing of the electromagnetic vector potential by requiring $The name is frequently confused with Hendrik Lorentz, who has given his name to many concepts in this field. The condition is Lorentz invariant. The Lorenz gauge condition does not completely determine the gauge: one can still make a gauge transformation where is the four-gradient and is any harmonic scalar function: that is, a scalar function obeying the equation of a massless scalar field.$

In electromagnetism, the electromagnetic tensor or electromagnetic field tensor is a mathematical object that describes the electromagnetic field in spacetime. The field tensor was first used after the four-dimensional tensor formulation of special relativity was introduced by Hermann Minkowski. The tensor allows related physical laws to be written very concisely, and allows for the quantization of the electromagnetic field by Lagrangian formulation described below.

The Grad–Shafranov equation is the equilibrium equation in ideal magnetohydrodynamics (MHD) for a two dimensional plasma, for example the axisymmetric toroidal plasma in a tokamak. This equation takes the same form as the Hicks equation from fluid dynamics. This equation is a two-dimensional, nonlinear, elliptic partial differential equation obtained from the reduction of the ideal MHD equations to two dimensions, often for the case of toroidal axisymmetry. Taking $as the cylindrical coordinates, the flux function is governed by the equation,$

In physics, the Majorana equation is a relativistic wave equation. It is named after the Italian physicist Ettore Majorana, who proposed it in 1937 as a means of describing fermions that are their own antiparticle. Particles corresponding to this equation are termed Majorana particles, although that term now has a more expansive meaning, referring to any fermionic particle that is its own anti-particle.

The electromagnetic wave equation is a second-order partial differential equation that describes the propagation of electromagnetic waves through a medium or in a vacuum. It is a three-dimensional form of the wave equation. The homogeneous form of the equation, written in terms of either the electric field $E$ or the magnetic field $B$ , takes the form:

The Maxwell stress tensor is a symmetric second-order tensor used in classical electromagnetism to represent the interaction between electromagnetic forces and mechanical momentum. In simple situations, such as a point charge moving freely in a homogeneous magnetic field, it is easy to calculate the forces on the charge from the Lorentz force law. When the situation becomes more complicated, this ordinary procedure can become impractically difficult, with equations spanning multiple lines. It is therefore convenient to collect many of these terms in the Maxwell stress tensor, and to use tensor arithmetic to find the answer to the problem at hand.

The covariant formulation of classical electromagnetism refers to ways of writing the laws of classical electromagnetism in a form that is manifestly invariant under Lorentz transformations, in the formalism of special relativity using rectilinear inertial coordinate systems. These expressions both make it simple to prove that the laws of classical electromagnetism take the same form in any inertial coordinate system, and also provide a way to translate the fields and forces from one frame to another. However, this is not as general as Maxwell's equations in curved spacetime or non-rectilinear coordinate systems.

There are various mathematical descriptions of the electromagnetic field that are used in the study of electromagnetism, one of the four fundamental interactions of nature. In this article, several approaches are discussed, although the equations are in terms of electric and magnetic fields, potentials, and charges with currents, generally speaking.

In mathematical physics, spacetime algebra (STA) is the application of Clifford algebra Cl_1,3(R), or equivalently the geometric algebra G(M⁴) to physics. Spacetime algebra provides a "unified, coordinate-free formulation for all of relativistic physics, including the Dirac equation, Maxwell equation and General Relativity" and "reduces the mathematical divide between classical, quantum and relativistic physics."

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.

In magnetohydrodynamics, the induction equation is a partial differential equation that relates the magnetic field and velocity of an electrically conductive fluid such as a plasma. It can be derived from Maxwell's equations and Ohm's law, and plays a major role in plasma physics and astrophysics, especially in dynamo theory.

In mathematical physics, the Gordon decomposition of the Dirac current is a splitting of the charge or particle-number current into a part that arises from the motion of the center of mass of the particles and a part that arises from gradients of the spin density. It makes explicit use of the Dirac equation and so it applies only to "on-shell" solutions of the Dirac equation.

References

↑ Holt, E. H.; Haskell, R. E. (1965). Foundations of Plasma Dynamics . New York: Macmillan. pp. 429-431.
↑ Chen, Francis F. (2016). Introduction to Plasma Physics and Controlled Fusion (3rd ed.). Heidelberg: Springer. pp. 192–194. ISBN 978-3-319-22308-7.
↑ Landau, L. D.; Lifshitz, E. M. (2013). The Classical Theory of Fields (4th revised ed.). New York: Elsevier. ISBN 9781483293288.
↑ Longcope, Dana (2002). "Notes on Magnetohydrodynamics" (PDF). Montana State University - Department of Physics. Retrieved 30 April 2019.
↑ Brauer, J. R. (2014). Magnetic Actuators and Sensors (2nd ed.). Hoboken NJ: Wiley IEEE Press. ISBN 978-1-118-50525-0.

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Holt, E. H.; Haskell, R. E. (1965). Foundations of Plasma Dynamics . New York: Macmillan. pp. 429-431.

[2] Chen, Francis F. (2016). Introduction to Plasma Physics and Controlled Fusion (3rd ed.). Heidelberg: Springer. pp. 192–194. ISBN 978-3-319-22308-7.

[3] Landau, L. D.; Lifshitz, E. M. (2013). The Classical Theory of Fields (4th revised ed.). New York: Elsevier. ISBN 9781483293288.

[4] Longcope, Dana (2002). "Notes on Magnetohydrodynamics" (PDF). Montana State University - Department of Physics. Retrieved 30 April 2019.

[5] Brauer, J. R. (2014). Magnetic Actuators and Sensors (2nd ed.). Hoboken NJ: Wiley IEEE Press. ISBN 978-1-118-50525-0.

[1]

[2]

[3]

[4]

[5]