Skip to content

Charge Conservation Inside Maxwell's Equations

Irodov 3.356DifficultyDifficulty 2 of 5UniquenessUniqueness 3 of 5

Problem Irodov 3.356

Show that the law of conservation of electric charge, ∇⋅j=−∂ρ/∂t\nabla\cdot\vect j = -\partial\rho/\partial t, is a consequence of Maxwell’s equations.

Given
  • ∇×H=j+∂D/∂t\nabla\times\vect H = \vect j + \partial\vect D/\partial t
  • ∇⋅D=ρ\nabla\cdot\vect D = \rho
Find
show that ∇⋅j=−∂ρ/∂t\nabla\cdot\vect j = -\partial\rho/\partial t

Solution

Step 1The divergence of the circulation law

The circulation law for H\vect H ties the curl of H\vect H to the conduction current and the displacement current:

∇×H=j+∂D∂t(1) \nabla\times\vect H = \vect j + \frac{\partial\vect D}{\partial t} \tag{1}

The divergence of a curl is zero for every smooth vector field, since the mixed second derivatives cancel in pairs:

∇⋅(∇×H)=0\nabla\cdot(\nabla\times\vect H) = 0

Taking the divergence of both sides of (1), and interchanging the space and time derivatives in the last term:

0=∇⋅j+∂∂t(∇⋅D)(2) 0 = \nabla\cdot\vect j + \frac{\partial}{\partial t}\left(\nabla\cdot\vect D\right) \tag{2}
Step 2Gauss’s law closes it

The remaining divergence is fixed by the other source equation, Gauss’s law, which names the density of free charge:

∇⋅D=ρ\nabla\cdot\vect D = \rho

Putting it into (2) and moving the time derivative to the right:

∇⋅j=−∂ρ∂t\nabla\cdot\vect j = -\frac{\partial\rho}{\partial t}
(3)

This is the law of charge conservation in local form, derived from the two source equations alone.

Step 3What it says about a volume

To read the law, integrate (3) over a fixed volume VV bounded by a closed surface SS, and turn the volume integral of the divergence into the flux through SS (Gauss’s theorem):

A fixed closed surface S round a region holding the charge q, with current density j leaving through S.

A fixed closed surface SS round a region holding the charge qq, with current density j\vect j leaving through SS.

∮Sj⋅dS=∫V∇⋅j dV=−∫V∂ρ∂t dV=−dqdt\oint_S \vect j\cdot d\vect S = \int_V \nabla\cdot\vect j\,dV = -\int_V \frac{\partial\rho}{\partial t}\,dV = -\frac{dq}{dt}

The current flowing out through SS equals the rate at which the charge qq inside VV decreases: charge is never created or destroyed, only carried.

Check

In a steady state ∂ρ/∂t=0\partial\rho/\partial t = 0, and (3) reduces to ∇⋅j=0\nabla\cdot\vect j = 0: as much current enters any volume as leaves it, which is Kirchhoff’s junction rule. Units: ∇⋅j\nabla\cdot\vect j is (A/m2)/m=C/(m3⋅s)(\text{A/m}^2)/\text{m} = \text{C}/(\text{m}^3\cdot\text{s}), the unit of ∂ρ/∂t\partial\rho/\partial t.

Remark

Without the displacement current the divergence of (1) would give ∇⋅j=0\nabla\cdot\vect j = 0 everywhere and always, which is false at the plate of a capacitor being charged. The term ∂D/∂t\partial\vect D/\partial t is exactly what makes the circulation law consistent with charge conservation; seen this way, Maxwell’s addition was forced rather than guessed.