Classical Electromagnetism, Fall 2013 Stress-Tensor Notes

Total Page:16

File Type:pdf, Size:1020Kb

Classical Electromagnetism, Fall 2013 Stress-Tensor Notes Physics 112: Classical Electromagnetism, Fall 2013 Stress-Tensor Notes Please read through these notes. There are a few bulleted questions in the discussion which you should include in your solutions for next seminar. There are also a few questions at the end of these notes which you should solve. 1 Introduction to the stress-tensor Press your hands together in the y − z plane and feel the force that one hand exerts on the other across a tiny area A{ say, one square millimeter of your hands' palms. That force, of course, is a vector F~ and has three independent components: it has a normal component (along the x- direction); it also has a tangential component: if you try to slide your hands past each other, you feel a component of force along the surface{ a \shear" force in the y and z directions. The normal and tangential (i.e. shear) directions are the natural directions into which we decompose F~ . Not only is the force F~ a vector, but we can associate a vector with the surface itself: A~ = Ax^. Where the direction of the surface element points perpendicularly to the surface itself. In this case we have chosen the negative side of the surface to be the −x side and the positive side to be +x. With this choice, the force F~ is that which the negative hand, on the −x side, exerts on the positive hand. Now, I assert that the force F~ is a linear function of the surface A~. This is another way of saying that if I double the magnitude of A~ then the force F~ also doubles. This is clearly true if the force per unit area is a constant, and we assume that the area over which the force acts is small enough we can take the force per unit area as constant (what is another word for `force per unit area'?). It should seem obvious that two vectors which are linearly related to each other fall under the purview of linear algebra. In particular, in general there exists a matrix which relates them. I am going to work in index notation in these notes using the Einstein summation notation. In this notation the force Fi and the surface Ai are related by a matrix Tij: Fi = TijAj: (1) The matrix Tij is also known as a tensor{ and in particular is the stress tensor. It is the set of coefficients which takes a small surface and reports the force that is exerted across that surface. 1 Let us write out in words what Tij is: i−component of force per unit area Tij = (2) across a surface perpendicular toe ^j 0 i−component of momentum that crosses a unit area 1 = @ which is perpendicular toe ^j per unit time; A (3) with the crossing pointing from − e^j to +e ^j Notice the fact that both of these descriptions are identical since force is just the rate of change of momentum. However, we will see that the second description will better describe the situation we have with the electromagnetic fields. The stresses inside a table with a heavy weight on it are described by the stress tensor Tij, as are the stress in a flowing fluid or plasma, in the electromagnetic field, and in any other physical medium. Accordingly, we shall use the stress tensor as an important (really central) mathematical tool in our study of force balance in kinetic theory, elasticity theory, fluid mechanics, plasma physics, general relativity, electromagnetism, or any other field theory. Figure 1: Two shear forces acting on a cube of length L. It is not obvious from its definition, but the stress tensor Tij is always symmetric in its two indices. To see this, consider small cube with side L in any medium (or field). The medium outside the cube exerts forces, and therefore also exerts torques, on the cube's faces. The z-component of the torque is produced by the shear forces on the front and back faces and on the left and right{ see Fig. 1. Let us use the stress-tensor in order to calculate the torque about the z-axis. In particular, the force acting on the negative y surface (with a normal pointing in the −y^-direction) is 2 Fi = −TiyL : (4) The torque caused by this force about thez ^-axis is then: L L3 N = − (^y × F~ ) ! N = T : (5) z 2 z z 2 zyk ky 2 We can use the properties of the Levi-Civita symbol to show that the only component of the stress-tensor which contributes to this component of the torque is Txy: L3 N = T : (6) z 2 zyx xy Finally, we know that xyz = +1 and that each subsequent permutation flips the sign{ therefore zyx = −1{ so we have L3 N = − T : (7) z 2 xy The shear forces on the other three faces of the cube also contribute to the torque aboutz ^. • Work out the torque on the remaining three faces and show that the total torque about thez ^ axis is 3 Nz = L (Tyx − Txy): (8) Next note that the moment of inertia of a uniform cube about its center is I = (1=12)ρL5 and we see that the cube has an angular acceleration about thez ^-direction: 1 dΩ L3(T − T ) = ρL5 z : (9) yx xy 12 dt Now imagine taking the size of the cube to zero{ L ! 0{ if Tyx − Txy remains non-zero, then what happens to the angular acceleration? Using this result, answer the question: • Why do we require Tyx = Txy? Convince yourself that the same argument works for the other faces of the cube so that we conclude the stress-tensor must be symmetric. 2 Examples of stress-tensors An example will make the stress-tensor more concrete. Inside a static perfect fluid there is an isotropic pressure P , so Txx = Tyy = Tzz = P (the normal forces per unit area across surfaces in the y − z plane, the z − x plane and the x − y plane are all equal to P . Using index notation the stress-tensor takes the form Tij = P δij: (10) To check this result, consider a surface A~ = An^ with area A oriented perpendicular to some arbitrary unit vectorn ^. The force that the fluid exerts across A~, in index notation, is Fi = TikAk = P δikAnk = P Anj; (11) so we can see, as expected, it is a normal force with a magnitude equal to the fluid pressure P times the surface area A. This is what it should be! 3 Momentum conservation and the stress-tensor The stress tensor plays a central role in the Newtonian law of momentum conservation. Recall the physical interpretation of Tij as the i-component of momentum that crosses a unit area perpendicular 3 to the j-direction per unit time. Apply this definition to the little cube in Fig. 1{ and now imagine that the cube is not a solid object, but is a volume into which and out of which momentum can flow. 2 The momentum that flows into the cube in unit time across the front face (at y = 0) is Tjy(0)L , 2 and across the back face (at y = L) is −Tjy(L)L ; and their sum is @T T (0)L2 − T (L)L2 ' jy L3: (12) jy jy @y • Argue a similar relationship is true for the other faces and that the rate of change of momentum- density inside of the cube is given by @(momentum density) @p i + r T = i + r T = 0; (13) @t j ij @t j ij where ~p is the momentum density. This has the standard form for any local conservation law: the time derivative of the density of some quantity (here momentum), plus the divergence of the flux of that quantity (here the momentum flux is the stress-tensor), is zero. 4 So, what is the stress-tensor? As we've seen in these notes, the stress tensor can be ‘defined’ in at least three, equivalent, ways (though we did not prove the equivalence explicitly!): I. The stress-tensor is a matrix which tells us the linear relationship between force acting on a ~ ~ surface, A, and the surface A itself: Fi = TijAj = TijAn^j. II. The stress-tensor tells us the momentum flux that crosses a unit area, A~. III. The stress-tensor is a tensor whose negative divergence is equal to the time rate of change of momentum-density. It is important to realize that all three of these ‘definitions’ are equivalent and each may be useful in different circumstances. 5 Problems You don't have to hand in solutions, but you should sketch out the solutions for yourself! 1. Equations of motion for a perfect fluid (a) Consider a perfect fluid that moves with a velocity ~v that varies in space and time{ ~v(x; y; z; t). Explain why the fluid’s momentum density is ρ~v, and explain why its momentum flux (i.e., stress- tensor) is Tij = P δij + ρvivj: (14) (b) Explain why the law of mass conservation for this fluid is @ρ + r~ · (ρ~v) = 0: (15) @t 4 (c) Explain why the total-time derivative d @ ≡ + ~v · r~ ; (16) dt @t describes the rate of change as measured by an observer who moves locally with the fluid, i.e. with velocity ~v. This is called the fluid's advective time derivative. (d) Use the advective time derivative to rewrite mass conservation as 1 dρ = −r~ · ~v: (17) ρ dt In words this says that the divergence of the fluid’s velocity field is minus the fractional rate of change of its density, as measured in the fluid’s local rest frame.
Recommended publications
  • Chapter 3 Dynamics of the Electromagnetic Fields
    Chapter 3 Dynamics of the Electromagnetic Fields 3.1 Maxwell Displacement Current In the early 1860s (during the American civil war!) electricity including induction was well established experimentally. A big row was going on about theory. The warring camps were divided into the • Action-at-a-distance advocates and the • Field-theory advocates. James Clerk Maxwell was firmly in the field-theory camp. He invented mechanical analogies for the behavior of the fields locally in space and how the electric and magnetic influences were carried through space by invisible circulating cogs. Being a consumate mathematician he also formulated differential equations to describe the fields. In modern notation, they would (in 1860) have read: ρ �.E = Coulomb’s Law �0 ∂B � ∧ E = − Faraday’s Law (3.1) ∂t �.B = 0 � ∧ B = µ0j Ampere’s Law. (Quasi-static) Maxwell’s stroke of genius was to realize that this set of equations is inconsistent with charge conservation. In particular it is the quasi-static form of Ampere’s law that has a problem. Taking its divergence µ0�.j = �. (� ∧ B) = 0 (3.2) (because divergence of a curl is zero). This is fine for a static situation, but can’t work for a time-varying one. Conservation of charge in time-dependent case is ∂ρ �.j = − not zero. (3.3) ∂t 55 The problem can be fixed by adding an extra term to Ampere’s law because � � ∂ρ ∂ ∂E �.j + = �.j + �0�.E = �. j + �0 (3.4) ∂t ∂t ∂t Therefore Ampere’s law is consistent with charge conservation only if it is really to be written with the quantity (j + �0∂E/∂t) replacing j.
    [Show full text]
  • This Chapter Deals with Conservation of Energy, Momentum and Angular Momentum in Electromagnetic Systems
    This chapter deals with conservation of energy, momentum and angular momentum in electromagnetic systems. The basic idea is to use Maxwell’s Eqn. to write the charge and currents entirely in terms of the E and B-fields. For example, the current density can be written in terms of the curl of B and the Maxwell Displacement current or the rate of change of the E-field. We could then write the power density which is E dot J entirely in terms of fields and their time derivatives. We begin with a discussion of Poynting’s Theorem which describes the flow of power out of an electromagnetic system using this approach. We turn next to a discussion of the Maxwell stress tensor which is an elegant way of computing electromagnetic forces. For example, we write the charge density which is part of the electrostatic force density (rho times E) in terms of the divergence of the E-field. The magnetic forces involve current densities which can be written as the fields as just described to complete the electromagnetic force description. Since the force is the rate of change of momentum, the Maxwell stress tensor naturally leads to a discussion of electromagnetic momentum density which is similar in spirit to our previous discussion of electromagnetic energy density. In particular, we find that electromagnetic fields contain an angular momentum which accounts for the angular momentum achieved by charge distributions due to the EMF from collapsing magnetic fields according to Faraday’s law. This clears up a mystery from Physics 435. We will frequently re-visit this chapter since it develops many of our crucial tools we need in electrodynamics.
    [Show full text]
  • Meaning of the Maxwell Tensor
    Meaning of the Maxwell Tensor We define the Maxwell Stress Tensor as 1 1 1 T = ε E E − δ E E + B B − δ B B ij 0 i j ij k k µ i j ij k k 2 0 2 While it looks scary, the Maxwell tensor is just a different way of writing the fact that electric fields exert forces on charges, and that magnetic fields exert forces on cur- rents. Dimensionally, Momentum T = (Area)(Time) This is the dimensionality of a pressure: Force = dP dt = Momentum Area Area (Area)(Time) Meaning of the Maxwell Tensor (2) If we have a uniform electric field in the x direction, what are the components of the tensor? E = (E,0,0) ; B = (0,0,0) 2 1 0 0 1 0 0 ε 1 0 0 2 1 1 0 E = ε 0 0 0 − 0 1 0 + ( ) = 0 −1 0 Tij 0 E 0 µ − 0 0 0 2 0 0 1 0 2 0 0 1 This says that there is negative x momentum flowing in the x direction, and positive y momentum flowing in the y direction, and positive z momentum flowing in the z di- rection. But it’s a constant, with no divergence, so if we integrate the flux of any component out of (not through!) a closed volume, it’s zero Meaning of the Maxwell Tensor (3) But what if the Maxwell tensor is not uniform? Consider a parallel-plate capacitor, and do the surface integral with part of the surface between the plates, and part out- side of the plates Now the x-flux integral will be finite on one end face, and zero on the other.
    [Show full text]
  • Section 2: Maxwell Equations
    Section 2: Maxwell’s equations Electromotive force We start the discussion of time-dependent magnetic and electric fields by introducing the concept of the electromotive force . Consider a typical electric circuit. There are two forces involved in driving current around a circuit: the source, fs , which is ordinarily confined to one portion of the loop (a battery, say), and the electrostatic force, E, which serves to smooth out the flow and communicate the influence of the source to distant parts of the circuit. Therefore, the total force per unit charge is a circuit is = + f fs E . (2.1) The physical agency responsible for fs , can be any one of many different things: in a battery it’s a chemical force; in a piezoelectric crystal mechanical pressure is converted into an electrical impulse; in a thermocouple it’s a temperature gradient that does the job; in a photoelectric cell it’s light. Whatever the mechanism, its net effect is determined by the line integral of f around the circuit: E =fl ⋅=d f ⋅ d l ∫ ∫ s . (2.2) The latter equality is because ∫ E⋅d l = 0 for electrostatic fields, and it doesn’t matter whether you use f E or fs . The quantity is called the electromotive force , or emf , of the circuit. It’s a lousy term, since this is not a force at all – it’s the integral of a force per unit charge. Within an ideal source of emf (a resistanceless battery, for instance), the net force on the charges is zero, so E = – fs. The potential difference between the terminals ( a and b) is therefore b b ∆Φ=− ⋅ = ⋅ = ⋅ = E ∫Eld ∫ fs d l ∫ f s d l .
    [Show full text]
  • Application Limits of the Airgap Maxwell Tensor
    Application Limits of the Airgap Maxwell Tensor Raphaël Pile, Guillaume Parent, Emile Devillers, Thomas Henneron, Yvonnick Le Menach, Jean Le Besnerais, Jean-Philippe Lecointe To cite this version: Raphaël Pile, Guillaume Parent, Emile Devillers, Thomas Henneron, Yvonnick Le Menach, et al.. Application Limits of the Airgap Maxwell Tensor. 2019. hal-02169268 HAL Id: hal-02169268 https://hal.archives-ouvertes.fr/hal-02169268 Preprint submitted on 1 Jul 2019 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Application Limits of the Airgap Maxwell Tensor Raphael¨ Pile1,2, Guillaume Parent3, Emile Devillers1,2, Thomas Henneron2, Yvonnick Le Menach2, Jean Le Besnerais1, and Jean-Philippe Lecointe3 1EOMYS ENGINEERING, Lille-Hellemmes 59260, France 2Univ. Lille, Arts et Metiers ParisTech, Centrale Lille, HEI, EA 2697 - L2EP - Laboratoire d’Electrotechnique et d’Electronique de Puissance, F-59000 Lille, France 3Univ. Artois, EA 4025, Laboratoire Systemes` Electrotechniques´ et Environnement (LSEE), F-62400 Bethune,´ France Airgap Maxwell tensor is widely used in numerical simulations to accurately compute global magnetic forces and torque but also to estimate magnetic surface force waves, for instance when evaluating magnetic stress harmonics responsible for electromagnetic vibrations and acoustic noise in electrical machines.
    [Show full text]
  • Chapter 6. Maxwell Equations, Macroscopic Electromagnetism, Conservation Laws
    Chapter 6. Maxwell Equations, Macroscopic Electromagnetism, Conservation Laws 6.1 Displacement Current The magnetic field due to a current distribution satisfies Ampere’s law, (5.45) Using Stokes’s theorem we can transform this equation into an integral form: ∮ ∫ (5.47) We shall examine this law, show that it sometimes fails, and find a generalization that is always valid. Magnetic field near a charging capacitor loop C Fig 6.1 parallel plate capacitor Consider the circuit shown in Fig. 6.1, which consists of a parallel plate capacitor being charged by a constant current . Ampere’s law applied to the loop C and the surface leads to ∫ ∫ (6.1) If, on the other hand, Ampere’s law is applied to the loop C and the surface , we find (6.2) ∫ ∫ Equations 6.1 and 6.2 contradict each other and thus cannot both be correct. The origin of this contradiction is that Ampere’s law is a faulty equation, not consistent with charge conservation: (6.3) ∫ ∫ ∫ ∫ 1 Generalization of Ampere’s law by Maxwell Ampere’s law was derived for steady-state current phenomena with . Using the continuity equation for charge and current (6.4) and Coulombs law (6.5) we obtain ( ) (6.6) Ampere’s law is generalized by the replacement (6.7) i.e., (6.8) Maxwell equations Maxwell equations consist of Eq. 6.8 plus three equations with which we are already familiar: (6.9) 6.2 Vector and Scalar Potentials Definition of and Since , we can still define in terms of a vector potential: (6.10) Then Faradays law can be written ( ) The vanishing curl means that we can define a scalar potential satisfying (6.11) or (6.12) 2 Maxwell equations in terms of vector and scalar potentials and defined as Eqs.
    [Show full text]
  • PHYS 110B - HW #4 Fall 2005, Solutions by David Pace Equations Referenced As ”EQ
    PHYS 110B - HW #4 Fall 2005, Solutions by David Pace Equations referenced as ”EQ. #” are from Griffiths Problem statements are paraphrased [1.] Problem 8.2 from Griffiths Reference problem 7.31 (figure 7.43). (a) Let the charge on the ends of the wire be zero at t = 0. Find the electric and magnetic fields in the gap, E~ (s, t) and B~ (s, t). (b) Find the energy density and Poynting vector in the gap. Verify equation 8.14. (c) Solve for the total energy in the gap; it will be time-dependent. Find the total power flowing into the gap by integrating the Poynting vector over the relevant surface. Verify that the input power is equivalent to the rate of increasing energy in the gap. (Griffiths Hint: This verification amounts to proving the validity of equation 8.9 in the case where W = 0.) Solution (a) The electric field between the plates of a parallel plate capacitor is known to be (see example 2.5 in Griffiths), σ E~ = zˆ (1) o where I define zˆ as the direction in which the current is flowing. We may assume that the charge is always spread uniformly over the surfaces of the wire. The resultant charge density on each “plate” is then time-dependent because the flowing current causes charge to pile up. At time zero there is no charge on the plates, so we know that the charge density increases linearly as time progresses. It σ(t) = (2) πa2 where a is the radius of the wire and It is the total charge on the plates at any instant (current is in units of Coulombs/second, so the total charge on the end plate is the current multiplied by the length of time over which the current has been flowing).
    [Show full text]
  • PHYS 332 Homework #1 Maxwell's Equations, Poynting Vector, Stress
    PHYS 332 Homework #1 Maxwell's Equations, Poynting Vector, Stress Energy Tensor, Electromagnetic Momentum, and Angular Momentum 1 : Wangsness 21-1. A parallel plate capacitor consists of two circular plates of area S (an effectively infinite area) with a vacuum between them. It is connected to a battery of constant emf . The plates are then slowly oscillated so that the separation d between ~ them is described by d = d0 + d1 sin !t. Find the magnetic field H between the plates produced by the displacement current. Similarly, find H~ if the capacitor is first disconnected from the battery and then the plates are oscillated in the same manner. 2 : Wangsness 21-7. Find the form of Maxwell's equations in terms of E~ and B~ for a linear isotropic but non-homogeneous medium. Do not assume that Ohm's Law is valid. 3 : Wangsness 21-9. Figure 21-3 shows a charging parallel plate capacitor. The plates are circular with radius a (effectively infinite). Find the Poynting vector on the bounding surface of region 1 of the figure. Find the total rate at which energy is entering region 1 and then show that it equals the rate at which the energy of the capacitor is increasing. 4a: Extra-2. You have a parallel plate capacitor (assumed infinite) with charge densities σ and −σ and plate separation d. The capacitor moves with velocity v in a direction perpendicular to the area of the plates (NOT perpendicular to the area vector for the plates). Determine S~. b: This time the capacitor moves with velocity v parallel to the area of the plates.
    [Show full text]
  • Energy, Momentum, and Symmetries 1 Fields 2 Complex Notation
    Energy, Momentum, and Symmetries 1 Fields The interaction of charges was described through the concept of a field. A field connects charge to a geometry of space-time. Thus charge modifies surrounding space such that this space affects another charge. While the abstraction of an interaction to include an intermedi- ate step seems irrelevant to those approaching the subject from a Newtonian viewpoint, it is a way to include relativity in the interaction. The Newtonian concept of action-at-a-distance is not consistent with standard relativity theory. Most think in terms of “force”, and in particular, force acting directly between objects. You should now think in terms of fields, i.e. a modification to space-time geometry which then affects a particle positioned at that geometric point. In fact, it is not the fields themselves, but the potentials (which are more closely aligned with energy) which will be fundamental to the dynamics of interactions. 2 Complex notation Note that a description of the EM interaction involves the time dependent Maxwell equa- tions. The time dependence can be removed by essentially employing a Fourier transform. Remember in an earlier lecture, the Fourier transformation was introduced. ∞ (ω) = 1 dt F (t) e−iωt F √2π −∞R ∞ F (t) = 1 dt (ω) eiωt √2π −∞R F Now assume that all functions have the form; F (~x, t) F (~x)eiωt → These solutions are valid for a particular frequency, ω. One can superimpose solutions of different frequencies through a weighted, inverse Fourier transform to get the time depen- dence of any function. Also note that this assumption results in the introduction of complex functions to describe measurable quantities which in fact must be real.
    [Show full text]
  • The Maxwell Stress Tensor and Electromagnetic Momentum
    Progress In Electromagnetics Research Letters, Vol. 94, 151–156, 2020 The Maxwell Stress Tensor and Electromagnetic Momentum Artice Davis1, * and Vladimir Onoochin2 Abstract—In this paper, we discuss two well-known definitions of electromagnetic momentum, ρA and 0[E × B]. We show that the former is preferable to the latter for several reasons which we will discuss. Primarily, we show in detail—and by example—that the usual manipulations used in deriving the expression 0[E×B] have a serious mathematical flaw. We follow this by presenting a succinct derivation for the former expression. We feel that the fundamental definition of electromagnetic momentum should rely upon the interaction of a single particle with the electromagnetic field. Thus, it contrasts with the definition of momentum as 0[E × B] which depends upon a (defective) integral over an entire region, usually all space. 1. INTRODUCTION The correct form of the electromagnetic energy-momentum tensor has been debated for almost a century. A large number of papers have been written [1], some of which concentrate upon a single particle interaction with the electromagnetic field in free space [2, 3]; others discuss the interaction between a field and a dielectric material upon which it impinges [4]. In works which discuss the Minkowski-Abraham controversy, the system considered is usually on the macroscopic level. On this level, singularities caused by point charges cannot be treated. We show that neglecting the effect due the properties of classical charges creates difficulties which cannot be resolved within the framework typically discussed by textbooks on classical electrodynamics [5]. Our present concept of electromagnetic field momentum density has two interpretations: (i) As 0[E × B], the scaled Poynting vector.
    [Show full text]
  • Chapter 5 Energy and Momentum
    Chapter 5 Energy and Momentum The equations established so far describe the behavior of electric and magnetic fields. They are a direct consequence of Maxwell’s equations and the properties of matter. Although the electric and magnetic fields were initially postulated to explain the forces in Coulomb’s and Ampere’s laws, Maxwell’s equations do not provide any information about the energy content of an electromagnetic field. As we shall see, Poynting’s theorem provides a plausible relationship between the electromagnetic field and its energy content. 5.1 Poynting’s Theorem If the scalar product of the field E and Eq. (1.33) is subtracted from the scalar product of the field H and Eq. (1.32) the following equation is obtained ∂B ∂D H ( E) E ( H) = H E j E . (5.1) · ∇ × − · ∇ × − · ∂t − · ∂t − 0 · Noting that the expression on the left is identical to (E H), integrating both ∇· × sides over space and applying Gauss’ theorem the equation above becomes ∂B ∂D (E H) n da = H + E + j0 E dV (5.2) ∂V × · − V · ∂t · ∂t · Z Z h i Although this equation already forms the basis of Poynting’s theorem, more insight is provided when B and D are substituted by the generally valid equations (1.19). 59 60 CHAPTER 5. ENERGY AND MOMENTUM Equation (5.2) then reads as ∂ 1 (E H) n da + D E + B H dV = j0 E dV (5.3) ∂V × · ∂t V 2 · · − V · Z Z h i Z 1 ∂P ∂E µ ∂M ∂H E P dV 0 H M dV.
    [Show full text]
  • Maxwell's Stress Tensor and Electromagnetic Momentum
    Maxwell's Stress Tensor and Electromagnetic Momentum Dr. Priti Mishra Momentum conservation is rescued in electrodynamics by the realization that the fields themselves carry momentum. This is not so surprising when you consider that we have already attributed energy to the fields. In the case of the two point charges in Fig. 1, whatever momentum is lost to the particles is gained by the fields. Only when the field momentum is added to the mechanical momentum of the charges is momentum conservation restored. Let us see how this works out quantitatively in the following sections. Figure 1: A point charge q1 moving along x-axis encounters an identical one, q2 proceeding in at the same speed along the y-axis. 1 Maxwell's Stress Tensor Let us calculate the total electromagnetic force on the charges in volume V: As we know the force acting on a point charge q moving with velocity v in an electric field E and magnetic field B is given by Lorentz force F = q(E + v × B) (1) 1 Similarly, the force per unit volume acting on a charge distribution ρ in a volume V : f = ρ(E + v × B) (2) As we know ρv = J where J is the current density. ) f = ρE + J × B (3) Let us write the above equation in terms of fields alone, eliminating ρ and J by using Maxwell's equations (i) and (iv): ρ r:E = 0 ) ρ = 0r:E (4) @E r × B = µ J + µ 0 0 0 @t 1 @E ) J = r × B − 0 (5) µ0 @t Substituting ρ and J from equations (4) and (5) into equation (3) we get 1 @E f = 0(r:E)E + r × B − 0 × B (6) µ0 @t Now @ @E @B (E × B) = × B + E × (7) @t @t @t or @E @ @B × B = (E × B) − E × (8) @t @t @t and from Maxwell's third equation (Faraday's law) @B = −r × E (9) @t Substituting above in equation (8) we get @E @ × B = (E × B) + E × (r × E) (10) @t @t Substituting equation (10) in equation (6) we obtain 1 @ f = 0(r:E)E+ (r × B) × B −0 (E × B) + E × (r × E) (11) µ0 @t or 1 @ f = 0 [(r:E)E − E × (r × E)] + (r × B) × B − 0 (E × B) µ0 @t (12) 2 Just to make things look more symmetrical, let's add a term (r:B)B; since r:B = 0, this costs us nothing.
    [Show full text]