Energy and momentum flux (formal)

Problem:

For a single charge, the rate of doing work by external fields B and E is qv∙E, in which v is the velocity of the charge.
(a)  Find the corresponding expression for a continuous distribution of charge and current and interpret it physically.
(b)  Use Maxwell's equations to express the result from part (a) in terms of the fields alone.
(c)  From your result in part b, verify Poynting's theorem  (∂u/∂t) + ∇∙S = -E∙j.
Find expressions for the terms u and S on the left-hand side, interpret those terms and ultimately the physical significance of the theorem that  (∂u/∂t) + ∇∙S = -E∙j represents.

Solution:

• Concepts:
Maxwell's equations, energy density and energy flux in the electromagnetic field
• Reasoning:
We are explicitly instructed to use Maxwell's equations to convert the expression obtained in part (a) into one that contains the fields alone.  When work is done on a charge by the field, field energy is converted into other forms of energy.  Energy conservation applies.
• Details of the calculation:
(a)  Let ρdV = dq be the amount of charge in a volume element dV.
dW/dt = dqv∙E = ρv∙EdV = j∙EdV = rate at which work is done by the field on the charges in dV.
(b)  E∙j = E∙(1/μ0)(×B) - ε0(∂E/∂t)∙E
= -(1/μ0)∇∙(E×B) + B∙(1/μ0)(×E) - ε0(∂E/∂t)∙E
= -(1/μ0)∇∙(E×B) - (1/μ0)B∙(1/μ0)(∂B/∂t) - ε0(∂E/∂t)∙E
= -(1/μ0)∇∙(E×B) - (∂/∂t)(1/(2μ0))B2 + (ε0/2)E2
Therefore
(∂/∂t)(1/(2μ0))B2 + (ε0/2)E2) + (1/μ0)∇∙(E×B) = -E∙j.
(c)  Interpretation:
Let u = (1/(2μ0))B2 + (ε0/2)E2 be the energy density in the electromagnetic field.
Let S = (1/μ0)(E×B) be the energy flux.
Then (∂u/∂t) + ∇∙S = -E∙j.
V∇∙S  dV = ∮S·n dA,  -∫V ∂u/∂t dV = ∮S·n dA + ∫V E∙j  dV
This is a statement of energy conservation.  The rate at which the field energy in a volume decreases equals the rate at which it leaves across the boundaries plus the rate at which it gets converted into other forms.

Problem:

Beginning with Newton's second law and the force law for a charged particle in an electromagnetic field, show that the time rate of change for the total momentum of the field plus the charge distribution in a volume V is given by
d/dt(Pfield  + Pcharge)i = ∮AjTij nj dA,
where T is Maxwell's stress tensor given by
Tij = ε0EiEj + (1/μ0)BiBj - ½(ε0E2 + (1/μ0)B2ij.
Give the physical significance of the integral on the right-hand side.

Solution:

• Concepts:
Newton's 2nd law, the Lorentz force, Maxwell's equations
• Reasoning:
If we denote the momentum of all the particles in a volume V by Pcharge then
dPcharge/dt =  ∑qq(E + v×B) = ∫Vd3r(ρE + j×B), with j = ρv.  We can now use Maxwell's equations to replace ρ and j by expressions in terms of the fields.
• Details of the calculation:
ρ  = ε0∇∙Ej = (1/μ0)(×B - (1/c2)∂E/∂t)
ρE + j×B = ε0E(∇∙E) + (1/μ0)(×B - (1/c2)∂E/∂t)×B
= ε0E(∇∙E) - (1/μ0)B××B + (1/(c2μ0))B×(∂E/∂t)
= ε0E(∇∙E) + (1/μ0)B(∇∙B) - (1/μ0)B××B - ε0E/∂t×B)
= ε0E(∇∙E) + (1/μ0)B(∇∙B) - (1/μ0)B××B - ε0∂/∂t(E×B) + ε0E×∂B/∂t
= -ε0∂/∂t(E×B) - (1/μ0)B××B - ε0E××E + ε0E(∇∙E) + (1/μ0)B(∇∙B).
Note:  We can add the term (1/μ0)B(∇∙B) since ∇∙B = 0.
Interpretation:
(1/μ0)(E×B) = energy flux, (1/(cμ0))(E×B) = momentum flux,
(1/(c2μ0))(E×B) = ε0(E×B) = momentum density.
We have
dPfield/dt = d/dt ∫Vd3r(1/(c2μ0))(E×B) = ∫Vd3r∂/∂t[(1/(c2μ0))(E×B)]
in some fixed volume V.
We can therefore write
d/dt(Pfield  + Pcharge) =  ∫Vd3r[ε0(E(∇∙E) - E××E) + (1/μ0)(B(∇∙B) - B××B)].
Write the right hand side in component form.
(E(∇∙E) - E××E)i = Ei(∂Ei/∂xi + ∂Ej/∂xj + ∂Ek/∂xk ) - Ej(×E)k + Ek(×E)j
= Ei(∂Ei/∂xi + ∂Ej/∂xj + ∂Ek/∂xk ) - Ej(∂Ej/∂xi - ∂Ei/∂xj) - Ek(∂Ek/∂xi - ∂Ei/∂xk)
= ½∂Ei2/∂xi + ∂EiEj/∂xj + ∂EiEk/∂xk - ½∂/∂xi(Ej2 + Ek2)
= ∂Ei2/∂xi + ∂EiEj/∂xj + ∂EiEk/∂xk - ½∂/∂xi(Ei2 + Ej2 + Ek2)
= ∑l∂/∂xl[EiEl) - ½ E∙E δil]
Similarly
B(∇∙B) - B××B)i = ∑l∂/∂xl[BiBl) - ½ B∙B δil]
Define:
Tij = ε0EiEj + (1/μ0)BiBj - ½(ε0E2 + (1/μ0)B2ij,
then  d/dt(Pfield  + Pcharge)i =  ∑l V ∂/∂xlTil dV.
The components Tii, Tij, Tik  form a vector.
l ∂/∂xlTil is the divergence of a vector field.
Use Stoke's theorem.
l V ∂/∂xlTil dV =  AlTil nl dA.
d/dt(Pfield  + Pcharge)i = ∮AlTil nl dA.
The integral represents the flow of the i = (x, y, z) component of the momentum across the surface of the volume V.