Today in Physics 28: the classic conservation laws in electrodynamics Poynting s theorem Energy conservation in electrodynamics The Maxwell stress tensor (which gets rather messy) Momentum conservation in electrodynamics Electromagnetism on the sun, doing work on matter and emitting radiation (TRACE satellite; Stanford U/Lockheed/NASA) 23 January 2004 Physics 28, Spring 2004
Poynting s theorem Suppose a collection of charges and currents lies entirely within a volume If they are released at some point in time, electromagnetic forces will begin to do work on them Consider, for instance, the work done by the forces on charges and currents in an infinitesimal volume dτ during a time dt: Thus, 2 v d Wmech = F d = dq E+ B vdt c = Ev dqdt ie B does no work = E Jdτ dt using dq = ρdτ, J = ρv dwmech = E Jdτ dt 23 January 2004 Physics 28, Spring 2004 2
Poynting s theorem (continued) We ve seen this before, of course; it s just another way to write P = I But we will learn something by elimination of J using Ampère s law: c E J = B : dwmech E = dτ ce ( B) E dt According to Product Rule #6, ( ) ( ) ( ) E B = E B + B E B = ( E B) B c t Use Faraday s law: 23 January 2004 Physics 28, Spring 2004 3
Poynting s theorem (continued) The time-derivative terms in all of these expressions can be written as E 2 E = ( E E) = E, 2 2 B 2 and, similarly, B = B 2 Using the last three results in the expression for dw/dt gives us dwmech B E = dτ c ( E B) B E dt = dτ c + B + E 2 ( ) ( 2 2 E B ) 23 January 2004 Physics 28, Spring 2004 4
Poynting s theorem (continued) is bounded by surface S Apply the divergence theorem to the first term in the integral, converting it to a surface integral over S: dwmech B E = dτ c ( E B) B E dt c = + 8π S ( ) ( 2 2 E B da dτ B E ) c d ( ) ( 2 2 = E B da B + E ) dτ 8π dt S In the last step we have used the fact that time is the only variable that survives the integration 23 January 2004 Physics 28, Spring 2004 5 Poynting s theorem
Energy conservation in electrodynamics Poynting s theorem of course expresses energy conservation We have long known that the energy density stored in electric and magnetic fields that is, the work required to assemble the charges and currents in the configuration they re in is ( 2 2 u ) EB = E + B 8π So the integral of u is evidently the energy within stored in the E and B fields, and the term containing it just gives the rate at which the stored energy changes 23 January 2004 Physics 28, Spring 2004 6
Energy conservation in electrodynamics (continued) The surface integral term contains a vector field that deserves its own name: c S= E B Poynting vector Thus, dwmech d dw = S da uebd d dt dt τ = S a dt S S This means that the rate at which work is done on the charges and currents by the fields is balanced not just by the rate at which the stored energy decreases, but also by a new term which, since it involves a flux integral over the surface bounding, must be the rate at which the fields carry energy out of EB 23 January 2004 Physics 28, Spring 2004 7
Energy conservation in electrodynamics (continued) Put another way: if the stored energy decreases, not all of the resulting work goes into kinetic energy of the charges some is radiated away The Poynting vector tells us the energy per unit area and time that is radiated Rewrite the last result as d dt d u dτ = Sdτ u dτ mech dt u mech We may equate the integrands now, and obtain: 23 January 2004 Physics 28, Spring 2004 8 EB ueb dτ = Sdτ dτ
Energy conservation in electrodynamics (continued) umech ueb = S, or ( umech + ueb ) + S = 0, an expression that invites comparison with the (chargecurrent) continuity equation, 0 t ρ + J = Think of the Poynting vector S therefore as the current density of energy 23 January 2004 Physics 28, Spring 2004 9
The Maxwell stress tensor We won t use tensors very much in this class, and you won t have to cope with this particular tensor after today But the path to expressions for momentum conservation goes straight through the Maxwell stress tensor, and this tensor does turn out to be a valuable tool in more-advanced courses, so it won t hurt to introduce it here Return to the Lorentz force law, and make a new definition: F = q E+ v B F dτ = ρ E+ v B dτ, where c c = ρ E+ v B = ρe+ J B c c is the force per unit volume exerted by E and B 23 January 2004 Physics 28, Spring 2004 0 F ( J = ρv)
The Maxwell stress tensor (continued) Using Gauss s law and Ampère s law, in the forms c E ρ = E, J = B, we can eliminate the source terms, in favor of the fields: E F= ( E) E+ B B c The very last term can be put into a more useful form by noting that E B ( E B) = B+ E c c c E = B E E c ( ) 23 January 2004 Physics 28, Spring 2004 Use Faraday s law:
The Maxwell stress tensor (continued) Thus, F= ( E) E+ ( B) B+ ( E B) E ( E) = E E E E B B E B c ( ) ( ) ( ) ( ) Since B = 0, it changes nothing if we add B Bto the second square brackets to make the whole thing look more symmetrical: ( ) F= ( ) ( ) E E E E + c ( B) B B ( B) ( E B) 23 January 2004 Physics 28, Spring 2004 2
The Maxwell stress tensor (continued) We can simplify the terms in square brackets if we reintroduce tensor notation Consider the three Cartesian directions to be represented by x, y, z =, 2, 3; for example, x r E x E r = y = r2, or E = Ey = E2 z r 3 E E z 3 Then the x-component of the first term in [] can be written as E E2 E3 ( E) E E ( E) = E + + r r2 r3 E2 E E E3 E2 + E3 r r2 r3 r 23 January 2004 Physics 28, Spring 2004 3
The Maxwell stress tensor (continued) It will pay to multiply this out and regroup: E E2 E3 [] = E + E + E r r r 2 3 E E E E + E + E E 2 3 2 2 3 3 r r2 r3 r = + + 2 r r r 2 ( E ) ( E E ) ( E E ) 2 3 2 3 2 ( E ) ( E ) 2 3 r 2 r 2 = + + + + r r r 2 r 2 ( ) ( ) ( ) ( 2 2 2 E ) EE2 EE3 E E2 E3 2 3 23 January 2004 Physics 28, Spring 2004 4 2 E
The Maxwell stress tensor (continued) 2 [ ] ( ) ( ) ( ) ( 2 2 2 = E ) + EE2 + EE3 E + E2 + E3 r r r 2 r 2 3 3 3 ( ) 2 2 = EE j E = EE j Eδj, r 2 j r r 2 j= j= j where we have reintroduced the Kronecker delta, δ This is for component ; for component i, we d get ij, i = j = 0, i j 2 ( E) E E ( E) = EE i i j E δij r 2 23 January 2004 Physics 28, Spring 2004 5 3 j= j
The Maxwell stress tensor (continued) Similarly, 3 2 ( B) B B ( B) = BB i i j B δij r 2 j= j Let us now define a nine-component object, ( 2 2 T = E E B B E B ) + + δ 2 ij i j i j ij with which we can re-write the ith component of both of the ugly terms in []: ( ) ( ) E E E E i 3 + ( ) ( ) = Tij B B B B i r 23 January 2004 Physics 28, Spring 2004 6, Maxwell stress tensor j= j
The Maxwell stress tensor (continued) The best way to envision the sum on the right-hand side is as matrix multiplication: 3 T T2 T3 Tij = T2 T22 T32 r j j r r2 r = 3 T3 T23 T33 The result is represented by a three-component object; that is, a vector We can use a vector-algebra-like symbolism for this operation, by using T to denote the second-rank tensor whose nine components are T ij The sum is represented thus: 3 Tij T, r j= j since it now looks so much like a divergence 23 January 2004 Physics 28, Spring 2004 7
The Maxwell stress tensor (continued) Finally, the force per unit volume becomes, in terms of the Maxwell stress tensor, F= ( ) ( ) E E E E + ( ) ( ) ( ) B B B B E B c = T geb, t where EB c = Momentum g E B S 2 c density 23 January 2004 Physics 28, Spring 2004 8
Momentum conservation in electrodynamics The total force on the charges within volume can be found by integrating F: dpmech = Fdτ = T geb dτ Use divergence dt theorem: d = T da gebdτ, dt S We can now identify peb gebdτ as the total momentum stored in the fields E and B, so that geb = E B is the c momentum density in the fields 23 January 2004 Physics 28, Spring 2004 9
Momentum conservation in electrodynamics (continued) Rearranging slightly, we have d ( p + peb ) = T da dt Interpretation: mech Changes in mechanical momentum (mvs for all the charges) and momentum stored in fields within are caused by the pressure and shear ( T ) distributed over the boundary surface S, and vice versa The terms on the diagonal of T, 2 2 ( 2 2) Ei + Bi E + B, 2 represent pressures in the E and B fields S Momentum conservation in electrodynamics 23 January 2004 Physics 28, Spring 2004 20
Momentum conservation in electrodynamics (continued) ( ) The off-diagonal terms,, 4 EE i j + BB i j π represent shear in the fields the kind of stress that causes strain in a direction different than the stress One can also write the momentum-conservation relation in differential form: dpmech = gmech dτ = T geb dτ, or dt ( gmech + geb ) T = 0, where g mech is the mechanical (mv-type) momentum per unit volume for the charges within 23 January 2004 Physics 28, Spring 2004 2
Summary Electric and magnetic fields can store energy and momentum: ( 2 2 u ) EB = E + B, energy/volume 8π geb E B, momentum/volume c or, for that matter, even angular momentum: angular momentum/ L EB = r geb = r ( E B), c volume and can transport energy (and thus all the other quantities): c S= E B energy/area/time Electromagnetic waves are the normal method of transport 23 January 2004 Physics 28, Spring 2004 22