Electromagnetic Induction
Faraday's law, mutual and self inductance, magnetic energy storage, and the physics of time-varying magnetic fields.
3.1 Faraday's Law of Induction
Faraday's law states that a changing magnetic flux through a circuit induces an electromotive force (EMF):
In differential form: $\nabla\times\mathbf{E} = -\partial\mathbf{B}/\partial t$. The minus sign is Lenz's law.
The EMF can arise from three distinct mechanisms: (1) a time-varying magnetic field with a stationary loop, (2) a moving loop in a static field, or (3) both simultaneously. The universal flux rule encompasses all cases.
Derivation 1: Faraday's Law from the Lorentz Force
For a moving loop in a static field, the EMF arises from the magnetic force on the charge carriers.
Step 1: EMF as work per unit charge
The EMF around a loop is $\mathcal{E} = \oint(\mathbf{f}/q)\cdot d\boldsymbol{\ell}$, where $\mathbf{f}$ is the force on a charge. For a conductor moving with velocity $\mathbf{v}$ in field $\mathbf{B}$:
$$\mathcal{E} = \oint(\mathbf{v}\times\mathbf{B})\cdot d\boldsymbol{\ell}$$
Step 2: Relate to flux change
The change in flux through the loop during time $dt$ arises from the area swept by each element $d\boldsymbol{\ell}$ moving with velocity $\mathbf{v}$: $d(\Phi_B) = \oint\mathbf{B}\cdot(\mathbf{v}\times d\boldsymbol{\ell})\,dt$.
Step 3: Use the scalar triple product identity
$\mathbf{B}\cdot(\mathbf{v}\times d\boldsymbol{\ell}) = (\mathbf{v}\times\mathbf{B})\cdot d\boldsymbol{\ell} = -d\boldsymbol{\ell}\cdot(\mathbf{B}\times\mathbf{v})$. Therefore:
$$\frac{d\Phi_B}{dt} = -\oint(\mathbf{v}\times\mathbf{B})\cdot d\boldsymbol{\ell} = -\mathcal{E}$$
Step 4: Recover Faraday's law
$$\boxed{\mathcal{E} = -\frac{d\Phi_B}{dt}}$$
Step 5: Generalization
When $\mathbf{B}$ itself is time-dependent, the electric field acquires a non-conservative part from $\nabla\times\mathbf{E} = -\partial\mathbf{B}/\partial t$. Both contributions combine into the universal flux rule.
3.2 Mutual and Self Inductance
Mutual inductance: The flux through loop 2 due to current $I_1$ in loop 1 is $\Phi_{21} = M_{21}I_1$. The Neumann formula gives:
Self inductance: A single loop with current $I$ links its own flux: $\Phi = LI$. The self-inductance $L$ depends only on the geometry of the loop. The back-EMF is $\mathcal{E} = -L\,dI/dt$.
Derivation 2: Neumann Formula for Mutual Inductance
We derive the symmetric Neumann formula starting from the vector potential.
Step 1: Vector potential of loop 1
The vector potential due to current $I_1$ in loop $\mathcal{C}_1$ is:
$$\mathbf{A}_1(\mathbf{r}) = \frac{\mu_0 I_1}{4\pi}\oint_{\mathcal{C}_1}\frac{d\boldsymbol{\ell}_1}{|\mathbf{r}-\mathbf{r}_1|}$$
Step 2: Flux through loop 2
By Stokes' theorem, $\Phi_{21} = \int_{\mathcal{S}_2}\mathbf{B}_1\cdot d\mathbf{a}_2 = \oint_{\mathcal{C}_2}\mathbf{A}_1\cdot d\boldsymbol{\ell}_2$.
Step 3: Substitute the vector potential
$$\Phi_{21} = \frac{\mu_0 I_1}{4\pi}\oint_{\mathcal{C}_2}\oint_{\mathcal{C}_1}\frac{d\boldsymbol{\ell}_1\cdot d\boldsymbol{\ell}_2}{|\mathbf{r}_2-\mathbf{r}_1|}$$
Step 4: Extract the mutual inductance
Since $\Phi_{21} = M_{21}I_1$, we read off: $$M_{21} = \frac{\mu_0}{4\pi}\oint_{\mathcal{C}_1}\oint_{\mathcal{C}_2}\frac{d\boldsymbol{\ell}_1\cdot d\boldsymbol{\ell}_2}{|\mathbf{r}_1-\mathbf{r}_2|}$$
Step 5: Symmetry
The double integral is symmetric under $1\leftrightarrow 2$, proving $M_{12} = M_{21} \equiv M$. This is the Neumann formula, establishing that mutual inductance is a purely geometric quantity independent of which loop carries the current.
3.3 Energy Stored in Magnetic Fields
The energy stored in an inductance carrying current $I$ is $W = \frac{1}{2}LI^2$. More generally, the magnetic field energy is:
Derivation 3: Magnetic Energy from Work Against Back-EMF
We derive the energy stored in an inductor by computing the work done against the back-EMF.
Step 1: Power delivered to an inductor
The back-EMF is $\mathcal{E} = -L\,dI/dt$. The power required to drive current against this EMF is $P = -\mathcal{E}\,I = LI\,dI/dt$.
Step 2: Integrate to find total energy
$$W = \int_0^T P\,dt = L\int_0^{I_f}I\,dI = \frac{1}{2}LI_f^2$$
Step 3: Express in terms of B-field
For a solenoid: $L = \mu_0 n^2 V$ and $B = \mu_0 nI$, so $W = \frac{1}{2}\mu_0 n^2 V I^2 = \frac{B^2}{2\mu_0}V$.
Step 4: Generalize to arbitrary fields
$$W = \frac{1}{2\mu_0}\int|\mathbf{B}|^2\,d^3r$$
The energy density is $u_B = B^2/(2\mu_0)$ J/m$^3$.
Step 5: Alternative form using A and J
Using $\mathbf{B} = \nabla\times\mathbf{A}$ and the vector identity $\mathbf{B}\cdot(\nabla\times\mathbf{A}) = \nabla\cdot(\mathbf{A}\times\mathbf{B}) + \mathbf{A}\cdot(\nabla\times\mathbf{B})$, with the surface term vanishing at infinity:
$$\boxed{W = \frac{1}{2}\int\mathbf{A}\cdot\mathbf{J}\,d^3r}$$
3.4 Lenz's Law and Eddy Currents
Lenz's law states that the induced current opposes the change in flux that produced it. This is a consequence of energy conservation and manifests as a minus sign in Faraday's law. In conducting media, time-varying fields induce eddy currents that resist flux penetration.
Derivation 4: Skin Depth in a Conductor
We derive the penetration depth of an oscillating field into a good conductor.
Step 1: Start from Maxwell's equations in a conductor
In a conductor with conductivity $\sigma$, Ohm's law gives $\mathbf{J} = \sigma\mathbf{E}$. For a good conductor, the displacement current is negligible compared to the conduction current, so $\nabla\times\mathbf{B} \approx \mu_0\sigma\mathbf{E}$.
Step 2: Derive the diffusion equation
Taking the curl of Faraday's law and substituting: $$\nabla^2\mathbf{B} = \mu_0\sigma\frac{\partial\mathbf{B}}{\partial t}$$
This is a diffusion equation with diffusivity $D = 1/(\mu_0\sigma)$.
Step 3: Assume harmonic time dependence
For $\mathbf{B} \propto e^{-i\omega t}$: $\nabla^2\mathbf{B} = -i\omega\mu_0\sigma\mathbf{B}$. For a planar geometry with $B = B_0 e^{ikz}$:
$$k^2 = i\omega\mu_0\sigma, \quad k = \frac{1+i}{\delta}$$
Step 4: Identify the skin depth
$$\boxed{\delta = \sqrt{\frac{2}{\omega\mu_0\sigma}}}$$
Step 5: Physical interpretation
The field decays as $B_0 e^{-z/\delta}\cos(z/\delta - \omega t)$. For copper at 60 Hz, $\delta \approx 8.5$ mm. At 1 GHz, $\delta \approx 2$ $\mu$m. This is why RF currents flow in a thin surface layer.
3.5 Coupled Circuits and Transformers
For two magnetically coupled circuits, the flux linkages are:
The coupling coefficient $k = M/\sqrt{L_1 L_2}$ satisfies $0 \leq k \leq 1$. A perfect transformer ($k = 1$) has $V_2/V_1 = N_2/N_1$ and $I_2/I_1 = N_1/N_2$.
Derivation 5: Energy in Coupled Inductors
We derive the total energy stored in two coupled inductors and the constraint $M^2 \leq L_1 L_2$.
Step 1: Build up current in loop 1 first
With $I_2 = 0$, bring $I_1$ to its final value. The work is $W_1 = \frac{1}{2}L_1 I_1^2$.
Step 2: Then build up current in loop 2
Holding $I_1$ fixed, bring $I_2$ to its final value. The work against loop 2's self-inductance is $\frac{1}{2}L_2 I_2^2$. Additionally, the changing $I_2$ induces an EMF in loop 1, requiring extra work $MI_1 I_2$.
Step 3: Total energy
$$\boxed{W = \frac{1}{2}L_1 I_1^2 + MI_1 I_2 + \frac{1}{2}L_2 I_2^2}$$
Step 4: Energy must be non-negative
Since $W \geq 0$ for all $I_1, I_2$, the quadratic form must be positive semi-definite. This requires the discriminant to be non-positive:
$$M^2 \leq L_1 L_2$$
Step 5: Coupling coefficient bound
Defining $k = M/\sqrt{L_1 L_2}$, the constraint becomes $k^2 \leq 1$. Equality holds only when all flux from one loop links the other (perfect coupling). In practice, $k \approx 0.95$-$0.99$ for iron-core transformers.
Historical Notes
Faraday (1831): Michael Faraday discovered electromagnetic induction through meticulous experiments with coils and magnets. His concept of "lines of force" and flux provided the geometric intuition that Maxwell later formalized mathematically.
Henry (1831): Joseph Henry independently discovered self-induction in America, finding that a long coiled wire produced a larger spark on breaking the circuit. The unit of inductance is named in his honor.
Lenz (1834): Heinrich Lenz formulated his law governing the direction of induced currents. The minus sign in Faraday's law encodes this principle, which is fundamentally a consequence of energy conservation.
Neumann (1845): Franz Neumann gave the first mathematical formulation of electromagnetic induction and derived the mutual inductance formula, placing Faraday's empirical discoveries on rigorous mathematical footing.
Applications
Electric Generators
Rotating coils in magnetic fields generate AC power via Faraday's law. Modern generators produce GW of power with efficiencies exceeding 98%, using the $\mathcal{E} = NBAw\sin(\omega t)$ principle.
Induction Heating
Eddy currents induced by RF coils heat conducting workpieces. The skin depth determines the heating profile. Used in metal hardening, cooking, and semiconductor processing.
SQUID Magnetometers
Superconducting loops exploit flux quantization ($\Phi = n\Phi_0$) and Josephson junctions to measure magnetic fields with sensitivity below $10^{-15}$ T, used in brain imaging and geology.
Wireless Charging
Mutual inductance between transmitter and receiver coils transfers power wirelessly. The coupling coefficient $k$ and quality factor $Q$ determine efficiency at the resonant frequency.
Simulation: Electromagnetic Induction
This simulation demonstrates Faraday's law with flux and EMF waveforms, mutual inductance dependence on coil separation, magnetic energy storage in a solenoid, and RL circuit transient response.
Electromagnetic Induction: Faraday's Law, Inductance & Energy
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server