Eddy currents and magnetic diffusion
Derive field diffusion in an Ohmic conductor, identify its assumptions, and connect spatial scale to decay time and skin depth.
On this page
The question for today
What you’ll understand
- Derive magnetic diffusion after explicitly neglecting displacement current.
- Compute modal decay times and sinusoidal skin depth.
- Connect eddy-current heating with magnetic energy loss and geometric scale.
A changing magnetic field drives electric currents inside a bulk conductor even when no external wire closes a visible loop. These eddy currents modify the field that induced them. Combining Faraday's law, Ohm's law, and magnetoquasistatic Ampère's law turns this feedback into a diffusion equation. Its spatial scales determine how long a field pattern survives and how deeply a periodic drive penetrates, rather than describing unrestricted electromagnetic-wave propagation.
Which Maxwell term was discarded?
Assume a stationary homogeneous conductor with constant positive conductivity and a linear isotropic, time-independent permeability. Ohm's law gives current density equal to conductivity times E. Neglect displacement current relative to conduction current, so curl B equals permeability times current density. For a sinusoidal process a basic material check is ; the geometry and frequencies must also justify a magnetoquasistatic description rather than significant radiation or propagation delay.
Taking the curl of Ohm's law combined with Ampère's law, then substituting Faraday's law, yields field diffusion. Zero magnetic divergence turns the curl of a curl into minus the vector Laplacian. Conductivity or permeability variations cannot simply be pulled through these derivatives. Motion would add an induction term, and nonlinear hysteretic magnetization would require a different constitutive model. Listing these exclusions identifies exactly which physical feedback the simple equation captures.
Short spatial scales decay faster
For an axial field depending on one transverse coordinate, a sine spatial mode is consistent with zero divergence. Two spatial derivatives return minus the squared wavenumber times that mode. Its amplitude therefore decreases exponentially with a time constant inversely proportional to wavenumber squared. Doubling the wavenumber shortens the decay time by a factor of four. The diffusion time over a length scale is proportional to permeability times conductivity times length squared, with numerical factors fixed by the actual boundary conditions.
A slab whose transient axial field vanishes at both faces has sine modes with wavenumbers equal to integer multiples of pi divided by slab thickness. A driven boundary field must first be separated from this homogeneous-boundary transient; otherwise the chosen sine expansion will violate the boundary data. Diffusion does not generally erase an imposed steady field. It erases departures from the steady solution while currents dissipate energy as heat.
Skin effect and the destination of field energy
A sinusoidal boundary drive produces a field with amplitude decaying over the skin depth and phase changing with depth. Raising frequency reduces penetration; raising conductivity also reduces penetration because induced currents oppose changing interior flux more effectively. This does not make a normal conductor a universal shield against static magnetic fields. In the zero-frequency limit its ideal skin depth grows without bound, and transient shielding eventually decays at finite conductivity.
In this linear model the magnetic energy density is B squared divided by twice permeability. Poynting's balance relates its change to boundary energy flux and positive Joule heating, with density current squared divided by conductivity. Thin insulated laminations interrupt large transverse current loops and reduce eddy-current loss in a suitable core geometry. This geometrical suppression differs from reducing hysteresis: an electrically resistive material can still dissipate energy through irreversible magnetization processes.
Work it through
A slab's slowest transient
A uniform stationary conductor has μσ=100 s/m² and thickness d=0.020 m. Its transient axial field is initially T and vanishes at x=0 and x=d. Find the decay time and peak field one decay time later.
The diffusivity is . The prescribed initial field is already a single fundamental mode, so no higher coefficients are needed to satisfy either boundary.
Its wavenumber is . Therefore s. The boundary-dependent pi-squared factor distinguishes this result from a rough diffusion-scale estimate.
At t=τ the peak amplitude is multiplied by , giving T. The nodes remain fixed at both slab faces; the profile amplitude decreases without a translating wavefront in this reduced model.
The fundamental decay time is 4.05 ms and the peak falls from 0.10 T to 0.0368 T after that interval. The result assumes the neglected displacement and material effects remain negligible.
The living notebook
Explore the model
Interactive lab · heat
The same modal diffusion mechanism
Reset and set both controls to Spatial mode n=2 and Diffusivity κ=0.4. Predict amplitude exp(−0.8), approximately 0.449, after one Step to t=0.5. For a fresh comparison, Reset and set both controls to n=4 and κ=0.4; one Step gives exp(−3.2), approximately 0.0408. The exact periodic scalar heat mode shares the decay mathematics of a permitted axial magnetic component.
- Current amplitude
- 1
- Spatial L² norm
- 1.772
- Decay rate · κn²
- 3.6 /s
The exact periodic heat-equation mode is u(x,t) = exp(−κn²t) sin(nx), on [0,2π]. Playback advances physical time. The initial sine is shown separately, and the unnormalized L² norm is √π exp(−κn²t). This analytic instrument does not run a finite-difference solver or test discretization stability. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.
Look forDoubling spatial mode multiplies the exponential decay rate by four. The lab solves one analytic periodic scalar mode; it contains no conducting boundary, skin-depth drive, Joule-heating readout, or Maxwell wave solver.
Pause & reason
Check your understanding
Your turn
Try it yourself
For a semi-infinite conductor with the same μσ=100 s/m², find the sinusoidal skin depth at 50 Hz and at 200 Hz. Assume constant material parameters and a valid magnetoquasistatic regime.
Take a moment to make a prediction or work it out on paper. Open a hint when you need one.
Give me a hint
Convert ordinary frequency to angular frequency using ω=2πf before inserting it into the skin-depth expression.
Compare with the solution
At 50 Hz, m, or 7.98 mm. Increasing the frequency fourfold to 200 Hz reduces the square-root expression by a factor of two, giving 3.99 mm. These are amplitude e-folding distances in sinusoidal steady state. They are neither the slab's thickness nor its transient decay time, and the lower penetration at higher frequency does not imply permanent exclusion of a static field.
Keep the idea
Take it with you
- Magnetic diffusion follows only after specified constitutive assumptions and neglect of displacement current.
- Modal decay time grows with spatial length squared, including boundary-dependent factors.
- Skin effect describes periodic penetration and differs from both static shielding and hysteresis loss.
Read further
Original sources and deeper explanations to follow at your own pace.
Optional progress, saved in this browser when available. Revisit any lesson at any time.
Lesson not marked complete.