← Part II: Magnetostatics & Induction
Jackson Chapter 5-6

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):

$$\boxed{\mathcal{E} = -\frac{d\Phi_B}{dt} = -\frac{d}{dt}\int_\mathcal{S}\mathbf{B}\cdot d\mathbf{a}}$$

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:

$$\boxed{M_{12} = 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|}}$$

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:

$$\boxed{W = \frac{1}{2\mu_0}\int |\mathbf{B}|^2\,d^3r = \frac{1}{2}\int \mathbf{A}\cdot\mathbf{J}\,d^3r}$$

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:

$$\Phi_1 = L_1 I_1 + M I_2, \qquad \Phi_2 = M I_1 + L_2 I_2$$

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

Python
script.py114 lines

Click Run to execute the Python code

Code will be executed with Python 3 on the server

Rate this chapter: