• 沒有找到結果。

Problem Set 5. Originally appeared at: http://sites.google.com/site/peeterjoot/math2011/relElectroDynProblemSet5.pdf Peeter Joot — peeter.joot@gmail.com

N/A
N/A
Protected

Academic year: 2022

Share "Problem Set 5. Originally appeared at: http://sites.google.com/site/peeterjoot/math2011/relElectroDynProblemSet5.pdf Peeter Joot — peeter.joot@gmail.com"

Copied!
15
0
0

加載中.... (立即查看全文)

全文

(1)

Problem Set 5.

Originally appeared at:

http://sites.google.com/site/peeterjoot/math2011/relElectroDynProblemSet5.pdf Peeter Joot — peeter.joot@gmail.com

Mar 14, 2011 relElectroDynProblemSet5.tex

Contents

1 Problem 1. Sinusoidal current density on an infinite flat conducting sheet. 1

1.1 Statement. . . 1

1.2 1-2. Determining the electromagnetic potentials. . . 2

1.3 3. Find the electric and magnetic fields outside the plane. . . 3

1.4 4. Give a physical interpretation of the results of the previous section. . . 4

1.5 5. Find the direction and magnitude of the energy flux outside the plane. . . 4

1.6 6. Sketch the intensity of the electromagnetic field far from the plane. . . 4

1.7 7. Continuity across the plane? . . . 5

2 Problem 2. Fields generated by an arbitrarily moving charge. 5 2.1 Statement. . . 5

2.2 0. Solution. Gradient and time derivatives of the retarded time function. . . 6

2.3 1. Solution. Computing the EM fields from the Lienard-Wiechert potentials. . . 8

2.4 2. Solution. EM fields from a uniformly moving source. . . 11

2.5 2. Solution. EM fields from a uniformly moving source along x axis.. . . 13

3 Problem 3. 15

4 Grading notes. 15

1. Problem 1. Sinusoidal current density on an infinite flat conducting sheet.

1.1. Statement

An infinitely thin flat conducting surface lying in the x−z plane carries a surface current density:

κ=e3θ(t)κ0sin ωt (1)

Here e3is a unit vector in the z direction, κ0is the peak value of the current density, and θ(t)is the theta function: θ(t <0) −0, θ(t >0) =1.

1. Write down the equations determining the electromagnetic potentials. Specify which gauge you choose to work in.

2. Find the electromagnetic potentials outside the plane.

3. Find the electric and magnetic fields outside the plane.

(2)

4. Give a physical interpretation of the results of the previous section. Do they agree with your qualitative expectations?

5. Find the direction and magnitude of the energy flux outside the plane.

6. Consider a point at some distance from the plane. Sketch the intensity of the electromagnetic field near this point as a function of time. Explain physically.

7. Consider now a point near the plane. Are the electric and magnetic fields you found contin- uous across the conducting plane? Explain.

1.2. 1-2. Determining the electromagnetic potentials.

Augmenting the surface current density with a delta function we can form the current density for the system

J=δ(y)κ=e3θ(t)δ(y)κ0sin ωt. (2) With only a current distribution specified use of the Coulomb gauge allows for setting the scalar potential on the surface equal to zero, so that we have

A= 4πJ

c (3)

E= −1 c

∂A

∂t (4)

B= ∇ ×B (5)

Utilizing our Green’s function

G(x, t) = δ(t− |x|/c)

|x| =δ3(x)δ(t), (6) we can invert our vector potential equation, solving for A

A(x, t) =

Z

d3x0dt0x0,t0G(xx0, t−t0)A(x0, t0)

=

Z

d3x0dt0G(xx0, t−t0)4πJ(x0, t0) c

=

Z

d3x0dt0δ(t−t0− |xx0|/c) |xx0|

4πJ(x0, t0) c

=

Z

d3x0J(x0, t− |xx0|/c c|xx0|

= 1 c

Z

dx0dy0dz0e3θ(t− xx0

/c)δ(y)κ0sin(ω(t− xx0

/c)) 1

|xx0|

= e3κ0 c

Z

dx0dz0θ(t− x− (x0, 0, z0) /c)sin(ω(t− x− (x0, 0, z0) /c)) 1

|x− (x0, 0, z0)|

= e3κ0 c

Z

dx0dz0θ

 t− 1

c q

(x−x0)2+y2+ (z−z0)2

sin ω



t−1cp(x−x0)2+y2+ (z−z0)2

p(x−x0)2+y2+ (z−z0)2

(3)

Now a switch to polar coordinates makes sense. Let’s use

x0−x=r cos α (7)

z0−z=r sin α (8)

This gives us

A(x, t) = e3κ0 c

Z

r=0

Z

α=0rdrdαθ

 t−1

c q

r2+y2

sin ω



t−1cpr2+y2

pr2+y2

= 2πe3κ0 c

Z

r=0rdrθ

 t−1

c q

r2+y2

sin

 ω



t−1cpr2+y2

pr2+y2

Since the theta function imposes a t− 1

c q

r2+y2 >0 (9)

constraint, equivalent to

c2t2>r2+y2, (10)

we can reduce the upper range of the integral and drop the theta function explicitly

A(x, t) = 2πe3κ0 c

Z

c2t2y2

r=0 rdr

sin ω



t−1cpr2+y2

pr2+y2 (11)

Here I got lazy and used Mathematica to help evaluate this integral, for an end result of A(x, t) = 2πκ0ω

c2 e3(1−cos(ω(t− |y|/c))). (12) 1.3. 3. Find the electric and magnetic fields outside the plane.

Our electric field can be calculated by inspection E= −1

c

∂A

∂t

= −2πκ0ω

2

c3 e3sin(ω(t− |y|/c)). (13) For the magnetic field we have

B= ∇ ×A

= −2πκ0ω

c2 e3× ∇(1−cos(ω(t− |y|/c))))

= 2πκ0ω

c2 (−sin(ω(t− |y|/c))e3× ∇ω(t− |y|/c)

= 2πκ0ω

2

c3 sin(ω(t− |y|/c))e3× ∇|y|

= 2πκ0ω

2

c3 sin(ω(t− |y|/c))e3×e2,

(4)

which gives us

B= −2πκ0ω

2

c3 e1sin(ω(t− |y|/c) (14) 1.4. 4. Give a physical interpretation of the results of the previous section.

It was expected that the lack of boundary on the conducting sheet would make the potential away from the plane only depend on the y components of the spatial distance, and this is precisely what we find performing the grunt work of the integration.

Given that we had a sinusoidal forcing function for our wave equation, it seems logical that we also find our non-homogeneous solution to the wave equation has sinusoidal dependence. We find that the sinusoidal current results in sinusoidal potentials and fields very much like one has in the electric circuits problem that we solve with phasors in engineering applications.

We find that the electric and magnetic fields are oriented parallel to the plane containing the surface current density, with the electric field in the direction of the current, and the magnetic field perpendicular to that, but having energy propagate outwards from the plane.

It’s kind of curious that despite introducing a step function in time for the current, that the end result appears to have no constraints on time t > 0, and that we have fields at all points in space and time from this current distribution, even the t < 0 case (ie: we have a non-causal system). I’d have expected a transient response to switching on the current. Perhaps this is because instantaneously inducing a current on all points in an infinite sheet is not physically realizable?

1.5. 5. Find the direction and magnitude of the energy flux outside the plane.

Our energy flux, the Poynting vector, is

S= c

 2πκ0ω2 c3

2

sin2(ω(t− |y|/c)e3×e1. (15) This is

S= πκ

20ω4

c5 sin2(ω(t− |y|/c)e2 = πκ

20ω4

2c5 (1−cos((t− |y|/c)))e2. (16) This energy flux is directed outwards along the y axis, with magnitude oscillating around an average value of

|hSi| = πκ

20ω4

2c5 . (17)

1.6. 6. Sketch the intensity of the electromagnetic field far from the plane.

I’m assuming here that this question does not refer to the flux intensity hSi, since that is con- stant, and boring to sketch.

The time varying portion of either the electric or magnetic field is proportional to

sin(ωtω|y|/c) (18)

We have a sinusoid as a function of time, of period T = 2π/ω where the phase is adjusted at each position by the factor ω|y|/c. Every increase of∆y=2πc/ω shifts the waveform back.

A sketch is attached.

(5)

1.7. 7. Continuity across the plane?

It is sufficient to consider either the electric or magnetic field for the continuity question since the continuity is dictated by the sinusoidal term for both fields.

The point in time only changes the phase, so let’s consider the electric field at t = 0, and an infinitesimal distance y= ±ec/ω. At either point we have

E(0,±ec/ω, 0, 0) = 2πκ0ω

2

c3 e3e (19)

In the limit as e → 0 the field strength matches on either side of the plane (and happens to equal zero for this t=0 case).

We have a discontinuity in the spatial derivative of either field near the plate, but not for the fields themselves. A plot illustrates this nicely

Figure 1: sin(t− |y|)

2. Problem 2. Fields generated by an arbitrarily moving charge.

2.1. Statement

Show that for a particle moving on a worldline parametrized by(ct, xc(t)), the retarded time trwith respect to an arbitrary space time point(ct, x), defined in class as:

|xxc(tr)| =c(t−tr) (20) obeys

∇tr= − xxc(tr)

c|xxc(tr)| =c(t−tr) −vc(tr)cot(xxc(tr)) (21)

(6)

and

∂tr

∂t = c|xxc(tr)|

c|xxc(tr)| =c(t−tr) −vc(tr)cot(xxc(tr)) (22) 1. Then, use these to derive the expressions for E and B given in the book (and in the class

notes).

2. Finally, re-derive the already familiar expressions for the EM fields of a particle moving with uniform velocity.

2.2. 0. Solution. Gradient and time derivatives of the retarded time function.

Let’s use notation something like our text [1], where the solution to this problem is outlined in

§63, and write

R(tr) =xxc(tr) (23)

R=|R| (24)

where

∂R

∂tr

= −vc. (25)

From R2=R·Rwe also have

2R∂R

∂tr

=2R·∂R

∂tr, (26)

so if we write

Rˆ = R

R, (27)

we have

R0(tr) = −Rˆ ·vc. (28) Proceeding in the manner of the text, we have

∂R

∂t

= ∂R

∂tr

∂tr

∂t

= −Rˆ ·vc∂tr

∂t. (29)

From20we also have

R=|xxc(tr)| =c(t−tr), (30) so

∂R

∂t =c

 1− ∂tr

∂t



. (31)

This and29gives us

(7)

∂tr

∂t = 1

1−Rˆ ·vcc (32)

For the gradient we operate on the implicit equation30again. This gives us

∇R= ∇(ct−ctr) = −c∇tr. (33) However, we can also use the spatial definition of R = |xxc(t0)|. Note that this distance R= R(tr)is a function of space and time, since tr = tr(x, t)is implicitly a function of the spatial and time positions at which the retarded time is to be measured.

∇R= ∇ q

(xxc(tr))2

= 1

2R∇(xxc(tr))2

= 1

R(xβ−xβc)eαα(xβ−xcβ(tr))

= 1

R(R)βeα(δαβαxcβ(tr))

We have only this bit ∂αxβc(tr)to expand, but that’s just going to require a chain rule expansion.

This is easier to see in a more generic form

∂ f(g)

∂xα = ∂ f

∂g

∂g

∂xα, (34)

so we have

∂xcβ(tr)

∂xα = ∂x

β c(tr)

∂tr

∂tr

∂xα, (35)

which gets us close to where we want to be

∇R= 1

R R− (R)β∂x

βc(tr)

∂tr eα

∂tr

∂xα

!

= 1

R RR·∂x

βc(tr)

∂tr

∇tr

!

Putting the pieces together we have only minor algebra left since we can now equate the two expansions of∇R

−c∇tr= RˆRˆ ·vc(tr)∇tr. (36) This is given in the text, but these in between steps are left for us and for our homework assignments! From this point we can rearrange to find the desired result

∇tr = −1 c

Rˆ

1−Rˆ ·vcc = −Rˆ c

∂tr

∂t (37)

(8)

2.3. 1. Solution. Computing the EM fields from the Lienard-Wiechert potentials.

Now we are ready to derive the values of E and B that arise from the Lienard-Wiechert poten- tials. We have for the electric field.

We’ll evaluate

E= −1 c

∂A

∂t − ∇φ B= ∇ ×B

For the electric field we’ll use the chain rule on the vector potential

∂A

∂t = ∂tr

∂t

∂A

∂tr. (38)

Similarly for the gradient of the scalar potential we have

φ=eα

∂φ

∂xα

=eα

∂φ

∂tr

∂tr

∂xα

= ∂φ

∂tr∇tr

= −Rˆ c

∂tr

∂t. Our electric field is thus

E= −∂tr

∂t

 1 c

∂A

∂tr

Rˆ c

∂φ

∂tr



(39) For the magnetic field we have

∇ ×A=eα× ∂A

∂xα

=eα× ∂A

∂tr

∂tr

∂xα. The magnetic field will therefore be found by evaluating

B= (∇tr) × ∂A

∂tr = −∂tr

∂t Rˆ

c × ∂A

∂tr (40)

Let’s compare this to ˆR×E

Rˆ ×E=Rˆ ×



∂tr

∂t

 1 c

∂A

∂trRˆ c

∂φ

∂tr



=Rˆ ×



∂tr

∂t 1 c

∂A

∂tr



(9)

This equals40, verifying that we have

B=Rˆ ×E, (41)

something that we can determine even without fully evaluating E.

We are now left to evaluate the retarded time derivatives found in39. Our potentials are

φ(x, t) = e R(tr)

∂tr

∂t (42)

A(x, t) = evc(tr) cR(tr)

∂tr

∂t (43)

It’s clear that the quantity ∂tr/∂t is going to show up all over the place, so let’s label it γtr. This is justified by comparing to a particle’s boosted rest frame worldline

ct0 x0



=γ

 1 −β

β 1

 ct 0



=

 γct

γβct



, (44)

where we have ∂t0/∂t=γ, so for the remainder of this part of this problem we’ll write γtr∂tr

∂t

= 1

1−Rˆ ·vcc. (45)

Using primes to denote partial derivatives with respect to the retarded time trwe have

φ0 =e



R

0

R2γtr +γ

t0r

R



(46) A0 =evc

c



R

0

R2γtr +γ

t0r

R

 +eac

c γtr

R , (47)

so the electric field is

E= −γtr 1

cA0Rˆ cφ0



= −tr c

 vc c



R

0

R2γtr+ γ

0tr

R

 + ac

c γtr

R −Rˆ



R

0

R2γtr + γ

0tr

R



= −tr c

 vc c

 c

R2γtr+ γ

0tr

R

 + ac

c γtr

R −Rˆ

 c

R2γtr +γ

t0r

R



= −tr cR



γtr ac

c +vc R −Rcˆ

R

 +γ0tr

vc

c −Rˆ

 . Here’s where things get slightly messy.

γ0tr =

∂tr 1 1−vcc ·Rˆ

= −γ2tr

∂tr

 1− vc

c ·Rˆ

=γ2tr

ac

c ·Rˆ +vc c ·Rˆ0

,

(10)

and messier

Rˆ0 =

∂tr R R

= R

0

R − RR

0

R2

= −vc

R − Rˆ(−c) R

= 1

R −vc+c ˆR , then a bit unmessier

γ0tr =γ2tr

ac

c ·Rˆ +vc c ·Rˆ0

=γ2tr

ac

c ·Rˆ + vc

cR · (−vc+c ˆR)

=γ2tr

 Rˆ ·ac

c + vc R

− v2c cR

 .

Now we are set to plug this back into our electric field expression and start grouping terms

E= −

2tr

cR

 ac

c + vc R − Rcˆ

R +γtr

 Rˆ ·ac

c + vc R

− v

2c

cR

vc

c −Rˆ

= −

3tr

cR

 ac c +vc

R −Rcˆ R



1−Rˆ ·vc c

 +

 Rˆ ·ac

c +vc R

− v

2c

cR

vc

c −Rˆ

= −

3tr

c2R

 ac



1−Rˆ ·vc c



+Rˆ ·ac

vc

c −Rˆ

3tr

cR

 vc

R −Rcˆ R



1−Rˆ · vc c

 +



Rˆ ·vc R

− v

2c

cR

vc

c −Rˆ

Using

a× (b×c) =b(a·c) −c(a·b) (48) We can verify that

ac

1−Rˆ ·vc c



+Rˆ ·acvc

c −Rˆ

= −ac+a ˆR·v

c −Rˆ ·acvc

c +Rˆ ·acRˆ

=Rˆ × ˆRvc c

ac, which gets us closer to the desired end result

E=

t3r

c2RRˆ × ˆRvc c

ac

−

3tr

cR2



vcRcˆ 

1−Rˆ ·vc c

 +



Rˆ ·vcv

2c

c

vc

c −Rˆ . (49)

(11)

It is also easy to show that the remaining bit reduces nicely, since all the dot product terms conveniently cancel



vcRcˆ 

1−Rˆ ·vc c

 +



Rˆ ·vcv

2c

c

vc

c −Rˆ

=c

 1− v

2c

c2

 ˆRv c



(50) This completes the exercise, leaving us with

E=

3tr

c2RRˆ × ˆRvc c

ac+

3tr

R2

 1−v2c

c2

 ˆRvc c

 B=Rˆ ×E.

(51)

Looking back to45where γtr was defined, we see that this compares to (63.8-9) in the text.

2.4. 2. Solution. EM fields from a uniformly moving source.

For a uniform source moving in space at constant velocity

xc(t) =vt, (52)

our retarded time measured from the spacetime point(ct, x)is defined implicitly by

R=|xxc(tr)| =c(t−tr). (53) Squaring this we have

x2+v2t2r−2trx·v= c2t2+c2t2r −2cttr, (54) or

(c2v2)t2r+2tr(−ct+x·v) =x2−c2t2. (55) Rearranging to complete the square we have



pc2v2trtc2x·v

√c2v2

2

=x2−c2t2+ (tc2x·v)2 c2v2

= (x2−c2t2)(c2v2) + (tc2x·v)2 c2v2

= x

2c2x2v2c4t2+c2t2v2+t2c4+ (x·v)22tc2(x·v) c2v2

= c

2(x2+t2v2−2t(x·v)) −x2v2+ (x·v)2 c2v2

= c

2(xvt)2− (x×v)2 c2v2

Taking roots (and keeping the negative so that we have tr = t− |x|/c for the v = 0 case, we have

(12)

r 1−v2

c2ctr = q 1 1−vc22

ct−x·v c −

r

(xvt)2x×v c

2!

, (56)

or with β= v/c, this is ctr = 1

1−β2



ct−x·β− q

(xvt)2− (x×β)2



. (57)

What’s our retarded distance R =ct−ctr? We get

R= β· (xvt) +p(xvt)2− (x×β)2 1−β2

. (58)

For the vector distance we get (with β· (xβ) = (β·x)β2) R= xvt+β· (xβ) +βp

(xvt)2− (x×β)2 1−β2

. (59)

For the unit vector ˆR=R/R we have

Rˆ = xvt+β· (xβ) +βp

(xvt)2− (x×β)2

β· (xvt) +p(xvt)2− (x×β)2 . (60) The acceleration term in the electric field is zero, so we are left with just

E=

3tr

R2

 1− v

2c

c2

 ˆRvc c



. (61)

Leading to γtr, we have

Rˆ ·β= β· (xvt+Rβ)

β· (xvt) +R , (62) where, following §38 of the text we write

R = q

(xvt)2− (x×β)2 (63) This gives us

γtr = β· (xvt) +R

R(1−β2) . (64)

Observe that this equals one when β=0 as expected.

We can also compute

Rˆ −β= x+β· (xβ) −vt+βp

(xvt)2− (x×β)2 β· (xvt) +p(xvt)2− (x×β)2β

= (xvt)(1−β2)

β· (xvt) +p(xvt)2− (x×β)2.

(13)

Our long and messy expression for the field is therefore

E=3tr 1

R2(1−β2)(Rˆβ)

=e β· (xvt) +R R(1−β2)

!3

(1−β2)2

(β· (xvt) +R)2(1−β2)(xvt)(1−β2) β· (xvt) +R

This gives us our final result

E=e 1

(R)3(1−β2)(xvt) (65) As a small test we observe that we get the expected result

E=e x

|x|3 (66)

for the β=0 case.

When v = Ve1this also recovers equation (38.6) from the text as desired, and if we switch to primed coordinates

x0 =γ(x−vt) (67)

y0 =y (68)

z0 =z (69)

(1−β2)r02= (x−vt)2+ (y2+z2)(1−β2), (70) we recover the field equation derived twice before in previous problem sets

E= e

(r0)3(x0, γy0, γz0) (71) 2.5. 2. Solution. EM fields from a uniformly moving source along x axis.

Initially I had errors in the vector treatment above, so tried with the simpler case using uniform velocity v along the x axis instead. Comparison of the two showed where my errors were in the vector algebra, and that’s now also fixed up.

Performing all the algebra to solve for trin

|x−vtre1| =c(t−tr), (72) I get

ctr = ctp(x−vt)2+ (y2+z2)(1β2)

1−β2 = −γ(βx0+r0) (73)

This matches the vector expression from 57 with the special case of v = ve1 so we at least started off on the right foot.

(14)

For the retarded distance R =ct−ctrwe get

R= β(x−vt) +p(x−vt)2+ (y2+z2)(1−β2)

1−β2 = γ(βx0+r0) (74)

This also matches 58, so things still seem okay with the vector approach. What’s our vector retarded distance

R=xβctre1

= (x−βctr, y, z)

= x−vt+βp

(x−vt)2+ (y2+z2)(1−β2)

1−β2 , y, z

!

= γ(x0+βr0), y0, z0 So

Rˆ = 1

γ(βx0+r0) γ(x0+βr0), y0, z0

= 1

βx0+r0



x0+βr0,y0 γ,z0

γ



Rˆ −β= 1

γ(βx0+r0) γ(x0+βr0), y0, z0− (β, 0, 0)

= 1

βx0+r0



x0(1−β2),y0 γ,z0

γ



= 1

γ(βx0+r0)(x−vt, y, z)

For ∂tr/∂t, using ˆRcalculated above, or from73calculating directly I get

∂tr

∂t = r

0+βx0

r0(1−β2) = γ(r0+βx0)

R , (75)

where, as in §38 of the text, we write R =

q

(x−vt)2+ (y2+z2)(1β2). (76) Putting all the pieces together I get

E=e(1−β2)(x−vt, y, z)

 γ(βx0+r0)



 γ(r0+βx0)

R

3

1

 γ2(βx0+r0)2 so we have

E= e1−β2

(R)3 (x−vt, y, z) (77)

This matches equation (38.6) in the text.

(15)

3. Problem 3.

FIXME: TODO.

4. Grading notes.

Only the first question was graded (I lost 1.5 marks). I got my units wrong when I integrated 11, and my ω/c factor in the result is circled. Should have also done a dimensional analysis check.

There was also a remark that the integral is zero if t <y/c, which introduces a θ function. I think that I incorrectly dropped my θ function, and should have retained.

FIXME: Revisit both of these and make sure to understand where exactly I went wrong.

References

[1] L.D. Landau and E.M. Lifshitz. The classical theory of fields. Butterworth-Heinemann, 1980. 2.2

參考文獻

相關文件

了⼀一個方案,用以尋找滿足 Calabi 方程的空 間,這些空間現在通稱為 Calabi-Yau 空間。.

(12%) Among all planes that are tangent to the surface x 2 yz = 1, are there the ones that are nearest or farthest from the origin?. Find such tangent planes if

After zooming in, the surface and the tangent plane become almost indistinguishable, as shown in the second graph.. (Here, the tangent plane is above the surface.) If we zoom

[r]

[r]

With the help of the pictures and the words below, write a journal entry about what happened.. Write at least

The molal-freezing-point-depression constant (Kf) for ethanol is 1.99 °C/m. The density of the resulting solution is 0.974 g/mL.. 21) Which one of the following graphs shows the

Please create a timeline showing significant political, education, legal and social milestones for women of your favorite country.. Use the timeline template to record key dates