From Maxwell to circuit theory
Why power systems can be analyzed as lumped circuits
Foundation
Definition
Power systems are ultimately governed by Maxwell’s equations: the four partial differential equations that describe how electric and magnetic fields evolve in space and time. Yet the working language of power engineering is not field theory, but circuit theory — voltages, currents, lumped resistors, inductors, capacitors, and Kirchhoff’s laws. This page develops the bridge between these descriptions, showing how the quasi-stationary approximation allows electromagnetic fields to be represented by lumped circuit elements at 50/60 Hz, and clarifying what is gained, and what is necessarily lost, in that simplification.
Mathematical notation
The equations below use the standard operators of vector calculus. For readers needing a brief refresher, here is what each symbol means.
| Symbol | Name | What it does |
|---|---|---|
| $\nabla f$ | Gradient of a scalar field $f(x,y,z)$ | A vector pointing in the direction of steepest increase of $f$, with magnitude equal to the rate of change. |
| $\nabla \cdot \mathbf{F}$ | Divergence of a vector field $\mathbf{F}$ | A scalar measuring how much $\mathbf{F}$ spreads outward from a point — positive at sources, negative at sinks, zero where field lines are conserved. |
| $\nabla \times \mathbf{F}$ | Curl of a vector field $\mathbf{F}$ | A vector measuring how much $\mathbf{F}$ rotates around a point. Its direction follows the right-hand rule. |
| $\mathbf{a} \cdot \mathbf{b}$ | Scalar (dot) product | A number: the projection of one vector onto another, $|\mathbf{a}||\mathbf{b}|\cos\theta$. |
| $\mathbf{a} \times \mathbf{b}$ | Vector (cross) product | A vector perpendicular to both $\mathbf{a}$ and $\mathbf{b}$, with magnitude $|\mathbf{a}||\mathbf{b}|\sin\theta$. |
| $\partial / \partial t$ | Partial time derivative | Rate of change of a field with respect to time, holding spatial position fixed. |
Readers seeking deeper treatment may consult Griffiths’ Introduction to Electrodynamics (chapter 1) or Jackson’s Classical Electrodynamics (chapter 1, appendix). The discussion below assumes only operational familiarity with these symbols, not derivations of their identities.
Maxwell’s equations
In their modern (Heaviside) form, the four equations of classical electromagnetism in a medium with free charge density $\rho_f$ and free current density $\mathbf{J}_f$ are:
$$\nabla \cdot \mathbf{D} = \rho_f, \qquad \nabla \cdot \mathbf{B} = 0,$$$$\nabla \times \mathbf{E} = -\frac{\partial \mathbf{B}}{\partial t}, \qquad \nabla \times \mathbf{H} = \mathbf{J}_f + \frac{\partial \mathbf{D}}{\partial t}.$$- $\mathbf{E}$ is the electric field, the force per unit charge experienced by a test charge.
- $\mathbf{D} = \varepsilon\mathbf{E}$ is the electric displacement field, accounting for the response of polarisable media (permittivity $\varepsilon$).
- $\mathbf{B}$ is the magnetic flux density, the quantity whose flux through a surface defines magnetic flux.
- $\mathbf{H} = \mathbf{B}/\mu$ is the magnetic field intensity, related to $\mathbf{B}$ through the permeability $\mu$ of the medium.
- $\mathbf{J}_f = \sigma\mathbf{E}$ is the conduction current density, with $\sigma$ the electrical conductivity (Ohm’s law in field form).
These four equations, together with the continuity equation $\nabla \cdot \mathbf{J} = -\partial \rho / \partial t$, describe every classical electromagnetic phenomenon — from radio waves to lightning, from antennas to transformers. Solving them in their full generality requires PDE techniques over three-dimensional geometries, which is impractical for systems of thousands of buses and tens of thousands of branches. The next two sections show how the geometry of a power system at $50$–$60$ Hz simplifies these PDEs into the algebraic relations of lumped-element circuit theory.
The following four subsections address each equation in turn, with historical context, physical intuition, a worked example, and its specific relevance to power-system analysis.
Gauss’s law for electricity
What it says
The divergence of the electric displacement field at a point equals the free charge density at that point. Equivalently: the electric flux through any closed surface equals the total free charge enclosed, divided by permittivity.
Physical intuition
Imagine electric field lines as streams flowing outward from positive charges and into negative ones. The divergence $\nabla \cdot \mathbf{D}$ at a point measures the net “outflow” of those streams. Where charge is present, the streams have to start or end. Where there is no charge, field lines pass through without diverging.
Relevance to power systems
Gauss’s law is the physical foundation of capacitance. Gauss’s law dictates how charge couples to the electric field in the dielectric between two conductors, yielding the lumped relation $Q = CV$, $\;i = C\,dv/dt$. Every shunt admittance term $b_{ij}$ in the AC power-flow equations ultimately traces back to a Gauss’s-law calculation on the geometry of a transmission line or transformer.
Gauss’s law for magnetism
What it says
The divergence of the magnetic flux density is zero everywhere. Equivalently: the total magnetic flux through any closed surface is zero — every magnetic field line that enters a closed region also leaves it.
Physical intuition
Magnetic field lines, unlike electric ones, never begin or end — they always close on themselves into loops. Around a bar magnet, field lines emerge from the north pole, sweep through space, and return to the south pole through the magnet’s interior. Around a current-carrying wire, field lines form closed circles. This topological constraint is the geometric content of Gauss’s law for magnetism.
Relevance to power systems
This law is the physical reason magnetic flux conservation holds in transformers and induction machines. The ideal transformer relation $V_1/V_2 = N_1/N_2$ is a direct consequence: the same magnetic flux $\Phi$ links both windings, so the induced voltage in each is proportional to its turn count. Without $\nabla \cdot \mathbf{B} = 0$, transformer design would be a far more complicated affair.
Faraday’s law
What it says
The curl of the electric field at a point equals the negative time derivative of the magnetic flux density at the same point. Equivalently: a time-varying magnetic flux through any loop induces an EMF around that loop, proportional to the rate of change of the flux. The minus sign (Lenz’s law) enforces the direction: the induced EMF opposes the change in flux that caused it. Without this sign, energy conservation would be violated.
Physical intuition
A changing magnetic field acts as a source of swirling electric field. If you wave a bar magnet through a coil of wire, the magnetic flux through the loop changes in time, and a circulating electric field is generated around the loop — which, in the presence of free charges, drives a current. The faster you wave, the larger the induced EMF.
Relevance to power systems
Faraday’s law underwrites every voltage source in the grid: hydro generators, thermal generators, wind turbines, and synchronous condensers all rely on a rotating field cutting through windings. It is also the physical basis of inductance: a coil with $N$ turns carrying current $i$ produces a flux $\Phi = Li/N$; if $i$ changes, the self-flux induces a back-EMF $v = L\,di/dt$. Every series reactance $x_{ij}$ in the AC power-flow equations is, ultimately, a Faraday’s-law calculation. Transformers are nothing but two coils linked by a shared magnetic flux: the secondary voltage is the EMF induced by the primary’s time-varying flux, exactly as Faraday demonstrated in 1831.
Ampère–Maxwell law
What it says
The curl of the magnetic field intensity at a point equals the sum of two contributions: the free current density flowing through that point, and the time derivative of the electric displacement field at that point. The second term is the celebrated displacement current — the single most consequential addition to physics in the nineteenth century.
Physical intuition
The Ampère-Maxwell law says that magnetic field loops are created in two ways: by currents (Ampère) and by changing electric fields (Maxwell). The second mechanism is what allows electromagnetic waves to exist far from any source — the wave carries its own field-generation mechanism through space.
Relevance to power systems
In the quasi-stationary regime at $50$–$60$ Hz, the displacement current is negligible everywhere except inside capacitive elements. In a transmission line it accounts for the capacitive coupling between conductors and ground (line charging); in a transformer, for inter-winding capacitance. In OPF, these effects collapse into the shunt susceptances $b^{\text{sh}}$ and admittances $g_{fr}, b_{fr}$ in the $\pi$ model of each branch.
The conduction current $\mathbf{J}_f = \sigma\mathbf{E}$ — the workhorse of Ohm’s law — is what flows in the conductors of transmission lines and through the windings of every electrical machine. Every resistance $r_{ij}$ in the AC power-flow equations traces back to this term in Ampère’s law, applied to a particular wire geometry.
The quasi-stationary regime
Power systems operate at 50 or 60 Hz. The electromagnetic wavelength at these frequencies is:
A transmission line of 500 km is one-tenth of a wavelength; a substation is a hundred-thousandth. For any device whose physical extent is much smaller than λ, the time it takes a wave to traverse it is negligible compared to the oscillation period, and fields can be treated as if they were quasi-stationary — varying in time, but spatially uniform within each device.
In this regime, three simplifications follow immediately: (1) The displacement current ∂D/∂t outside of explicit capacitive elements is negligible. Maxwell’s fourth equation reduces to ∇ × H ≈ J inside conductors. (2) Spatial fields collapse to lumped quantities. Voltage between two points becomes a single scalar, current through a cross-section becomes a single scalar, charge stored on a capacitor becomes a single scalar. (3) Each component is characterized by a small set of parameters — R for resistors, L for inductors, C for capacitors — which encode the integrated effect of distributed physics over the component’s volume. Cross-component coupling reduces to topological wiring: which terminal connects to which.
Deriving Kirchhoff’s laws
Kirchhoff’s current law (KCL)
The continuity equation, integrated over the volume enclosing a node, gives:
In the quasi-stationary regime, charge does not accumulate at nodes. The volume integral vanishes, and the surface integral becomes a sum of currents entering the node:
This is Kirchhoff’s current law — a direct consequence of charge conservation combined with the lumped-element approximation.
Kirchhoff’s voltage law (KVL)
Faraday’s law in integral form is:
For a closed circuit loop, lumping inductive flux into discrete inductors Lk with vk = Lk dik/dt and treating capacitors and resistors similarly, the equation becomes:
This is Kirchhoff’s voltage law — Faraday’s law specialized to the lumped-circuit setting where all flux is contained inside named inductors.
Constitutive element laws
Each lumped element’s voltage-current relation follows from integrating its constitutive physics over its volume.
Resistor
For a conductor with conductivity σ and uniform cross-section, J = σE integrated over length and area gives:
This is Ohm’s law. The resistance R packages the conductor’s geometry (length ℓ, area A) and material (σ) into a single number.
Inductor
For a coil of N turns enclosing a magnetic flux Φ, Faraday’s law gives v = N dΦ/dt. For linear magnetic media, Φ is proportional to current: Φ = (L/N) i, where L is the self-inductance. Combining:
The inductance L encapsulates the coil geometry and the surrounding magnetic medium’s permeability μ.
Capacitor
For two conductors separated by a dielectric, the displacement current ∂D/∂t integrated over the dielectric’s cross-section gives i = C dv/dt:
The capacitance C packages the conductor areas, their separation, and the dielectric’s permittivity ε into a single parameter.
From time to phasor
In sinusoidal steady state at angular frequency ω, the element laws translate to phasor relations (see Phasors foundation for details):
The differential operator d/dt becomes algebraic multiplication by jω. Impedance Z = R + jX absorbs both resistive and reactive behavior into a single complex number per element. Kirchhoff’s laws apply directly to phasors, and the entire network reduces to the linear complex system I = YV analyzed in the Bus admittance matrix foundation.
When the approximation breaks
The quasi-stationary approximation fails when device dimensions approach the electromagnetic wavelength. For 60 Hz this requires sizes on the order of hundreds of kilometers — a regime entered only by the longest transmission lines, where distributed-parameter modeling (the telegrapher’s equations) replaces lumped π-models. Within typical power-system geographic scales, the lumped-circuit approximation is exact to better than 1% accuracy.
Higher-frequency phenomena — switching transients, lightning surges, harmonics from power electronics — push toward shorter effective wavelengths and may require electromagnetic transient (EMT) simulation rather than steady-state phasor analysis. OPF and power-flow problems live in the safe quasi-stationary regime; transient stability and power-quality studies push at its boundaries.
Further reading
- Jackson, J. D. (1999). Classical Electrodynamics (3rd ed.). Wiley. — The canonical electromagnetic theory reference; Chapter 6 derives quasi-stationarity.
- Griffiths, D. J. (2017). Introduction to Electrodynamics (4th ed.). Cambridge University Press. — Accessible derivation of Maxwell-to-circuit reduction.
- Haus, H. A., & Melcher, J. R. (1989). Electromagnetic Fields and Energy. Prentice Hall. — MIT OCW companion text; explicit treatment of lumped-element validity.
- Bergen, A. R., & Vittal, V. (2000). Power Systems Analysis (2nd ed.). Prentice Hall. — Bridge from field theory to power-system equations.
- Kundur, P. (1994). Power System Stability and Control. McGraw-Hill. — Appendix on electromagnetic vs. electromechanical timescales.