Appendix/Ramblings/Dyson1893Part1

From jetwiki
Jump to navigation Jump to search

Dyson (1893a) Part I: Some Details

This chapter provides some derivation details relevant to our accompanying discussion of Dyson's analysis of the gravitational potential exterior to an anchor ring.

Overview

In his pioneering work, F. W. Dyson (1893a, Philosophical Transactions of the Royal Society of London. A., 184, 43 - 95) and (1893b, Philosophical Transactions of the Royal Society of London. A., 184, 1041 - 1106) used analytic techniques to determine the approximate equilibrium structure of axisymmetric, uniformly rotating, incompressible tori. C.-Y. Wong (1974, ApJ, 190, 675 - 694) extended Dyson's work, using numerical techniques to obtain more accurate — but still approximate — equilibrium structures for incompressible tori having solid body rotation. Since then, Y. Eriguchi & D. Sugimoto (1981, Progress of Theoretical Physics, 65, 1870 - 1875) and I. Hachisu, J. E. Tohline & Y. Eriguchi (1987, ApJ, 323, 592 - 613) have mapped out the full sequence of Dyson-Wong tori, beginning from a bifurcation point on the Maclaurin spheroid sequence.

External Potential in Terms of Angle ψ

Step 1

On p. 59, at the end of §6 of Dyson (1893a), we find the following expression for the potential at point "P", anywhere exterior to an anchor ring:

πV(r,θ)M

=

(r,θ,c)+a2231cddc[(r,θ,c)]a42631cddc{1cddc[(r,θ,c)]}+

 

 

+(1)n+12a2n2n+2[135(2n3)224262(2n)2](1cddc)n[(r,θ,c)]+

where (see beginning of §8 on p. 61),

Anchor Ring Schematic
Anchor Ring Schematic

Caption: Anchor ring schematic, adapted from figure near the top of §2 (on p. 47) of Dyson (1893a)

(r,θ,c)

0πdϕ[r22crsinθcosϕ+c2]1/2

 

=

20π/2dϕ[R12(R12R2)sin2ϕ]1/2

 

=

2R10π/2dϕ[1(R12R2R12)sin2ϕ]1/2

 

=

2K(k)R1,

and, where furthermore,

K(k)

=

0π/2dϕ[1k2sin2ϕ]1/2

      and      

k

[R12R2R12]1/2.

Step 2

Taking a queue from our accompanying discussion of toroidal coordinates, if we adopt the variable notation,

ηln(R1R),

then we can write,

coshη=12[eη+eη]

=

R2+R122RR1,

which implies that,

[2cothη+1]1/2=[1e2η]1/2

=

[1(RR1)2]1/2=k.

Now, if we employ the Descending Landen Transformation for the complete elliptic integral of the first kind, we can make the substitution,

K(k)

=

(1+μ)K(μ),

      where,      

μ

11k21+1k2.

But notice that, 1k2=eη, in which case,

μ

=

1eη1+eη

=

1R/R11+R/R1

=

R1RR1+R.

Hence, we can write,

(r,θ,c)=2K(k)R1

=

2R1[(1+μ)K(μ)]

 

=

2K(μ)R1[1+R1RR1+R]

 

=

4K(μ)R1+R.

This is the expression for (r,θ,c) that was adopted by Dyson at the beginning of his §8.

Step 3

Subsequently, Dyson was able to obtain analytic expressions for successive derivatives of the function, (r,θ,c), by first demonstrating that

Comment by J. E. Tohline on 17 September 2018: In the middle of p. 61 of Dyson(1893a), there appears to be a typographical error in the expression for the derivative of R1 with respect to c; as we have indicated here, the numerator should be 4cR1 instead of 4cR.
Comment by J. E. Tohline on 17 September 2018: In the middle of p. 61 of Dyson(1893a), there appears to be a typographical error in the expression for the derivative of R1 with respect to c; as we have indicated here, the numerator should be 4cR1 instead of 4cR.

dRdc

=

4c2+R2R124cR,

dR1dc

=

4c2+R12R24cR1,       and,

dμdc

=

μccosψ,

where — as shown above in the Anchor ring schematic — ψ is the angle between R and R1 for which (according to the law of cosines),

cosψ

=

R2+R124c22RR1.

It will be useful for us to note that,

d(cosψ)dc

=

ddc[R2+R124c22RR1]

 

=

12RR1ddc[R2+R124c2]+(R2+R124c2)ddc[12RR1]

 

=

12RR1[2RdRdc+2R1dR1dc8c]+(R2+R124c2)[12R2R1dRdc12RR12dR1dc]

 

=

14cRR1[4c2+R2R12+4c2+R12R224c2](R2+R124c2){[4c2+R2R128cR3R1]+[4c2+R12R28cR13R]}

 

=

2cRR1cosψ4cR2R12[4c2(R12+R2)+2R12R2R14R4]

 

=

cosψ4cR2R12[R14+R44c2(R12+R2)2R12R2]2cRR1

 

=

14cR2R12{cosψ[R14+R42R12R24c2(R12+R2)]8c2RR1}

But,

[R14+R42R12R24c2(R12+R2)]

=

(R12R2)24c2(R12+R2)

 

=

[(R1+R)(R1R)]24c2(R12+R2)

 

=

(R12+2R1R+R2)(R1R)24c2(R12+R2)

 

=

2R1R(R1R)2+(R12+R2)(R1R)24c2(R12+R2)

 

=

2R1R(R1R)2+(R12+R2)[R122R1R+R24c2]

 

=

2R1R[(R1R)2+(R12+R2)(cosψ1)]

 

=

2R1R[(R12+R2)cosψ2R1R]

 

=

(2R1R)2[cos2ψ+(2c2R1R)cosψ1]

Hence,

d(cosψ)dc

=

(2R1R)24cR2R12{cosψ[cos2ψ+(2c2R1R)cosψ1]8c2RR1(2R1R)2}

 

=

1c{cosψ[cos2ψ+(2c2R1R)cosψ1]2c2R1R}

2cd(cosψ)dc

=

2cosψ[cos2ψ+(2c2R1R)cosψ1]4c2R1R

Step 4

Then, drawing upon known expressions for the derivatives of elliptic integrals, as are now tidily catalogued online in NIST's Digital Library of Mathematical Functions, Dyson showed that,


LaTeX mathematical expressions cut-and-pasted directly from
NIST's Digital Library of Mathematical Functions

According to §19.4 of NIST's Digital Library of Mathematical Functions,

dK(k)dk

=

E(k)k2K(k)kk2,

d(E(k)k2K(k))dk

=

kK(k),

dE(k)dk

=

E(k)K(k)k,

d(E(k)K(k))dk

=

kE(k)k2,

d2E(k)dk2

=

1kdK(k)dk=k2K(k)E(k)k2k2,

where,

k1k2.


d(r,θ,c)dc=ddc[4K(μ)R1+R]

=

4K(μ)ddc[R1+R]1+4(R1+R)[E(μ)μ2K(μ)μμ2]dμdc

 

=

4K(μ)(R1+R)2[dR1dc+dRdc]4(R1+R)[K(μ)μ]μccosψ+4(R1+R)[E(μ)μμ2]μccosψ

 

=

4K(μ)(R1+R)2[4c2+R12R24cR1+4c2+R2R124cR][4K(μ)c(R1+R)]cosψ+4E(μ)c(R1+R)[(R1+R)24RR1]cosψ

 

=

[(R1+R)E(μ)cRR1]cosψ[4K(μ)c(R1+R)]cosψ4K(μ)4cR1R(R1+R)2[R(4c2+R12R2)+R1(4c2+R2R12)]

 

=

[(R1+R)E(μ)cRR1]cosψ[4K(μ)c(R1+R)]cosψK(μ)cR1R(R1+R)2[(4c2+R1R)(R1+R)(R13+R3)]

 

=

[(R1+R)E(μ)cRR1]cosψ[4K(μ)c(R1+R)]cosψK(μ)cR1R(R1+R)[(4c2+R1R)(R12+R2R1R)]

 

=

[(R1+R)E(μ)cRR1]cosψ[4K(μ)c(R1+R)]cosψK(μ)cR1R(R1+R)[4c2(R1R)2]

 

=

[(R1+R)E(μ)cRR1]cosψ[4K(μ)c(R1+R)]cosψK(μ)c(R1+R)[2(1cosψ)]

 

=

[(R1+R)E(μ)cRR1]cosψ2K(μ)c(R1+R)(1+cosψ)

1cd(r,θ,c)dc

=

1c2{[E(μ)(R+R1)RR1]cosψ[4K(μ)R+R1]cos2ψ2}.

This expression appears at the top of Dyson's p. 62.

Step 5

Differentiating a second time gives,

ddc[1cd(r,θ,c)dc]

=

ddc{[(R1+R)E(μ)c2RR1]cosψ2K(μ)c2(R1+R)(1+cosψ)}

 

=

2c3{[(R1+R)E(μ)RR1]cosψ2K(μ)(R1+R)(1+cosψ)}

 

 

+1c2ddc{[(R1+R)E(μ)RR1]cosψ2K(μ)(R1+R)(1+cosψ)}

 

=

2c3{[(R1+R)E(μ)RR1]cosψ2K(μ)(R1+R)(1+cosψ)}

 

 

+1c2ddc{[(R1+R)E(μ)RR12K(μ)(R1+R)]cosψ}1c2ddc[2K(μ)(R1+R)]

 

=

2c3{[(R1+R)E(μ)RR1]cosψ2K(μ)(R1+R)(1+cosψ)}+1c2[(R1+R)E(μ)RR12K(μ)(R1+R)]d(cosψ)dc

 

 

+cosψc2ddc[(R1+R)E(μ)RR1](1+cosψ)2c2ddc[4K(μ)(R1+R)]

 

=

1c3[4K(μ)(R1+R)][(1+cosψ)(c2)d(cosψ)dc]2c3[(R1+R)E(μ)RR1][cosψ(c2)d(cosψ)dc]

 

 

+cosψc2ddc[(R1+R)E(μ)RR1](1+cosψ)2c3{[E(μ)(R+R1)RR1]cosψ[4K(μ)R+R1]cos2ψ2}

 

=

1c3[4K(μ)(R1+R)][(1+cosψ)222+(1+cosψ)(c2)d(cosψ)dc]

 

 

2c3[(R1+R)E(μ)RR1][(cosψ+cos2ψ)4+cosψ(c2)d(cosψ)dc]+cosψc2ddc[(R1+R)E(μ)RR1]

 

=

14c3[4K(μ)(R1+R)][5+6cosψ+cos2ψ2cd(cosψ)dc]

 

 

12c3[(R1+R)E(μ)RR1][5cosψ+cos2ψ2cd(cosψ)dc]+cosψc2ddc[(R1+R)E(μ)RR1].

Now,

ddc[(R1+R)E(μ)RR1]

=

[(R1+R)RR1]dE(μ)dμdμdc+E(μ)ddc[(R1+R)RR1]

 

=

[(R1+R)RR1][E(μ)K(μ)c]cosψ+E(μ){1RR1[dR1dc+dRdc](R1+R)[1RR12dR1dc+1R2R1dRdc]}

 

=

[(R1+R)RR1][E(μ)K(μ)c]cosψE(μ){[1R2dRdc]+[1R12dR1dc]}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)[(R1+R)cRR1]cosψE(μ){[4c2+R2R124cR3]+[4c2+R12R24cR13]}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)[(R1+R)cRR1]cosψ+E(μ){[R12+R24c24cR3][R22cR3]+[R12+R24c24cR13][R122cR13]}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)[(R1+R)cRR1]cosψ+E(μ)2cR12R2{R13cosψRR12+R3cosψR1R2}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)[(R1+R)cRR1]cosψ+E(μ)2cR12R2{(R13+R3)cosψRR1(R1+R)}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)[(R1+R)cRR1]cosψ+E(μ)(R1+R)2cR12R2{(R12+R24c2)cosψ+4c2cosψRR1(1+cosψ)}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)[(R1+R)cRR1]cosψ+E(μ)(R1+R)cR1R{cos2ψ+2c2cosψR1R12(1+cosψ)}

 

=

K(μ)[(R1+R)cRR1]cosψ+E(μ)(R1+R)cR1R[cos2ψ+12(cosψ1)+2c2cosψR1R].

We therefore have,

1cddc[1cd(r,θ,c)dc]

=

14c4[4K(μ)(R1+R)]{5+6cosψ+cos2ψ2cosψ[cos2ψ+(2c2R1R)cosψ1]+4c2R1R}

 

 

12c4[(R1+R)E(μ)RR1]{5cosψ+cos2ψ2cosψ[cos2ψ+(2c2R1R)cosψ1]+4c2R1R}

 

 

+cosψc4{E(μ)(R1+R)R1R[cos2ψ+12(cosψ1)+2c2cosψR1R]K(μ)[(R1+R)RR1]cosψ}

 

=

1c4[K(μ)(R1+R)]{5+6cosψ+cos2ψ2cosψ[cos2ψ+(2c2R1R)cosψ1]+4c2R1R}

 

 

12c4[(R1+R)E(μ)RR1]{5cosψ+cos2ψ2cosψ[cos2ψ+(2c2R1R)cosψ1]+4c2R1R}

 

 

12c4[E(μ)(R1+R)R1R][2cos3ψcos2ψ+cosψ4c2cos2ψR1R]

 

 

+1c4[K(μ)(R1+R)][(R1+R)2cos2ψRR1]

 

=

1c4[K(μ)(R1+R)]{5+8cosψ+cos2ψ2cos3ψ+4c2R1R[R12+R24c2R1R]cos2ψ[2R1RR1R]cos2ψ[8c2R1R]cos2ψ}

 

 

12c4[(R1+R)E(μ)RR1]{8cosψ4cos3ψ+(4c2R1R)(12cos2ψ)}

 

=

1c4[K(μ)(R1+R)]{5+8cosψcos2ψ4cos3ψ+4c2R1R(12cos2ψ)}

 

 

1c4[(R1+R)E(μ)RR1]{4cosψ2cos3ψ+(2c2R1R)(12cos2ψ)}

 

=

1c4[K(μ)(R1+R)]{5+8cosψcos2ψ4cos3ψ(4c2R1R)cos2ψ}

 

 

+1c4[(R1+R)E(μ)RR1]{4cosψ+2cos3ψ+(2c2R1R)cos2ψ},

which exactly matches the second equation from the top of p. 62 in Dyson (1893a).

Step 6 (Summary)

In summary, then,

πV(r,θ)M

=

(r,θ,c)+a2231cddc[(r,θ,c)]a42631cddc{1cddc[(r,θ,c)]}+𝒪(a5c5)

 

=

4K(μ)R1+R+123a2c2{[E(μ)(R+R1)RR1]cosψ[4K(μ)R+R1]cos2ψ2}

 

 

1263a4c4{[(R1+R)E(μ)RR1][4cosψ+2cos3ψ+(2c2R1R)cos2ψ]

 

 

+[K(μ)(R1+R)][5+8cosψcos2ψ4cos3ψ(4c2R1R)cos2ψ]}+𝒪(a5c5).

In preparation for a comparison between Dyson's expression and the expression derived by Wong (1973), let's regroup terms according to the various trigonometric functions. Keep in mind that,

cos2ψ2

=

12(1+cosψ);

cos2ψ

=

12(1+cos2ψ);

cos3ψ

=

14(3cosψ+cos3ψ).

Hence, we have,

πV(r,θ)M

=

4K(μ)R1+R+123a2c2{[E(μ)(R+R1)RR1]cosψ[2K(μ)R+R1](1+cosψ)}

 

 

1263a4c4{[(R1+R)E(μ)RR1][4cosψ+12(3cosψ+cos3ψ)+(2c2R1R)cos2ψ]

 

 

+[K(μ)(R1+R)][5+8cosψ12(1+cos2ψ)(3cosψ+cos3ψ)(4c2R1R)cos2ψ]}+𝒪(a5c5)

 

=

4K(μ)R1+R123a2c2[2K(μ)R+R1]+123a2c2[E(μ)(R+R1)RR1]cosψ123a2c2[2K(μ)R+R1]cosψ

 

 

1273a4c4[(R1+R)E(μ)RR1][5cosψ+(4c2R1R)cos2ψ+cos3ψ]

 

 

1273a4c4[K(μ)(R1+R)][9+10cosψcos2ψ(8c2R1R)cos2ψ2cos3ψ]+𝒪(a5c5)

 

=

2K(μ)R1+R[2123a2c2328a4c4+𝒪(a5c5)]

 

 

+{5263a4c4[K(μ)(R1+R)]+123a2c2[E(μ)(R+R1)RR1]123a2c2[2K(μ)R+R1]+5273a4c4[(R1+R)E(μ)RR1]}cosψ

 

 

+{1273a4c4[K(μ)(R1+R)][1+(8c2R1R)]1273a4c4[(R1+R)E(μ)RR1](4c2R1R)}cos2ψ

 

 

+{1263a4c4[K(μ)(R1+R)]1273a4c4[(R1+R)E(μ)RR1]}cos3ψ+𝒪(a5c5)

 

=

2K(μ)R1+R[2123a2c2328a4c4+𝒪(a5c5)]

 

 

[2K(μ)R+R1E(μ)(R+R1)RR1][123a2c2+5273a4c4+𝒪(a5c5)]cosψ

 

 

+{[2K(μ)(R1+R)][12+(4c2R1R)][(R1+R)E(μ)RR1](4c2R1R)}[1273a4c4+𝒪(a5c5)]cos2ψ

 

 

+[2K(μ)(R1+R)(R1+R)E(μ)RR1][1273a4c4+𝒪(a5c5)]cos3ψ

 

=

2K(μ)R1+R[2123a2c2328a4c4+𝒪(a5c5)]

 

 

[2K(μ)R+R1E(μ)(R+R1)RR1][123a2c2+5273a4c4+𝒪(a5c5)]cosψ

 

 

+[2K(μ)(R1+R)(R1+R)E(μ)RR1](4c2R1R)[1273a4c4+𝒪(a5c5)]cos2ψ+[2K(μ)(R1+R)][1283a4c4+𝒪(a5c5)]cos2ψ

 

 

+[2K(μ)(R1+R)(R1+R)E(μ)RR1][1273a4c4+𝒪(a5c5)]cos3ψ

Now, drawing from our accompanying discussion of relationships between complete elliptic integrals and relationships between their parameter arguments and, in particular, associating Dyson's parameter, μ=(R1R)/(R1+R), with k1, we understand that,

2E(μ)1+μ(1μ)K(μ)

=

E(k),

where the alternate parameter,

k=2μ1+μ=[1(RR1)2]1/2=[2cothη+1]1/2=1e2η.

Hence, we appreciate that,

E(k)

=

2E(μ)[1+R1RR1+R]1[1R1RR1+R]K(μ)

 

=

2E(μ)[2R1R1+R]1[2RR1+R]K(μ)

 

=

R[E(μ)(R1+R)RR12K(μ)R1+R]

As a result, we can rewrite Dyson's expression for the external potential as,

πV(r,θ)M

=

2K(μ)R1+R[2123a2c2328a4c4+𝒪(a5c5)]+E(k)R[123a2c2+5273a4c4+𝒪(a5c5)]cosψ

 

 

E(k)R(4c2R1R)[1273a4c4+𝒪(a5c5)]cos2ψ+[2K(μ)(R1+R)][1283a4c4+𝒪(a5c5)]cos2ψ

 

 

E(k)R[1273a4c4+𝒪(a5c5)]cos3ψ

 

=

2K(μ)R1+R{[2123a2c2328a4c4+𝒪(a5c5)]+[1283a4c4+𝒪(a5c5)]cos2ψ}

 

 

+E(k)R{[123a2c2+5273a4c4+𝒪(a5c5)]cosψ(4c2R1R)[1273a4c4+𝒪(a5c5)]cos2ψ[1273a4c4+𝒪(a5c5)]cos3ψ}.

External Potential in Terms of Angle χ

Anchor Ring Schematic
Anchor Ring Schematic

Caption: Anchor ring schematic, adapted from figure near the top of §2 (on p. 47) of Dyson (1893a)

Is it relatively straightforward to develop a similar expression for the external potential that is written in terms of the angle, χ, instead of (as above) in terms of the angle, ψ ? This would make the transition to Dyson's Paper II smoother. Initially I have in mind making the transformation via the law of sines, whereby,

R1sinχ

=

2csinψ

=

Rsin(πψχ).

This means that the following associations may be used as well:

sinψ

=

(2cR1)sinχ;

cos2ψ=1sin2ψ

=

1[(2cR1)sinχ]2.

From an accompanying discussion that builds upon the law of cosines, we also may write,

2(R1c)1

=

1+12(ac)cosχ+123(ac)2[3cos2χ1]+124(ac)3[5cos3χ3cosχ]

 

 

+127(ac)4[330cos2χ+35cos4χ]+𝒪(a5c5).


Compare Notations Used by Dyson and Wong

Dyson Anchor Ring Schematic
Dyson Anchor Ring Schematic

Dyson's Notation (slightly modified)

Wong Anchor Ring Schematic
Wong Anchor Ring Schematic

Wong's Notation (slightly modified)


Key relationship, with understanding that, ϵd/c:

a2

=

c2(1ϵ2)

           

ac

=

(1ϵ2)1/2.

Relationship Between Lengths

Subtraction

In both cases, lengths are referenced to the same "external" point whose meridional-plane coordinates are, (ϖ,z). Hence,

R2

=

z2+(cϖ)2

r22

=

z2+(aϖ)2

R2r22

=

(cϖ)2(aϖ)2

 

=

(c22cϖ+ϖ2)(a22aϖ+ϖ2)

 

=

c2a22(ca)ϖ

ϖ

=

c2a2+r22R22(ca).

and,

R12

=

z2+(c+ϖ)2

r12

=

z2+(a+ϖ)2

R12r12

=

(c+ϖ)2(a+ϖ)2

 

=

(c2+2cϖ+ϖ2)(a2+2aϖ+ϖ2)

 

=

c2a2+2(ca)ϖ

ϖ

=

R12r12(c2a2)2(ca).

Put together, then,

c2a2+r22R22(ca)

=

R12r12(c2a2)2(ca)

(R2+R12)(r12+r22)

=

2(c2a2)

 

=

2c22c2(1ϵ2)

 

=

2c2ϵ2

Multiplication

R2R12

=

[z2+(cϖ)2][z2+(c+ϖ)2]

0

=

z4+z2[(cϖ)2+(c+ϖ)2]+[(cϖ)2(c+ϖ)2R2R12]

2z2

=

±[[(cϖ)2+(c+ϖ)2]24[(cϖ)2(c+ϖ)2R2R12]]1/2[(cϖ)2+(c+ϖ)2]

r12r22

=

[z2+(aϖ)2][z2+(a+ϖ)2]

0

=

z4+z2[(aϖ)2+(a+ϖ)2]+[(aϖ)2(a+ϖ)2r12r22]

2z2

=

±[[(aϖ)2+(a+ϖ)2]24[(aϖ)2(a+ϖ)2r12r22]]1/2[(aϖ)2+(a+ϖ)2]

Note: In order to be physically relevant, the quantity, 2z2, must be positive. Given that the second term on the RHS is always negative, we conclude that only a solution with the superior (plus) sign is physically relevant.

Hence,

[[(cϖ)2+(c+ϖ)2]24[(cϖ)2(c+ϖ)2R2R12]]1/2[(cϖ)2+(c+ϖ)2]

=

[[(aϖ)2+(a+ϖ)2]24[(aϖ)2(a+ϖ)2r12r22]]1/2[(aϖ)2+(a+ϖ)2]

Ratios

Let's work with the following (dimensionless) length ratios:

λ22R2r22

=

z2+(cϖ)2z2+(aϖ)2,

λ12R12r12

=

z2+(c+ϖ)2z2+(a+ϖ)2.

Solving for z2 in the first case, we find,

λ22[z2+(aϖ)2]

=

z2+(cϖ)2

z2(λ221)

=

(cϖ)2λ22(aϖ)2

z2

=

(cϖ)2λ22(aϖ)2(λ221);

and in the second case, we find,

λ12[z2+(a+ϖ)2]

=

z2+(c+ϖ)2

z2(λ121)

=

(c+ϖ)2λ12(a+ϖ)2

z2

=

(c+ϖ)2λ12(a+ϖ)2(λ121).

Combined, this gives,

(cϖ)2λ22(aϖ)2(λ221)

=

(c+ϖ)2λ12(a+ϖ)2(λ121)

(λ121)[(cϖ)2λ22(aϖ)2]

=

(λ221)[(c+ϖ)2λ12(a+ϖ)2]

(λ121)[(c22cϖ+ϖ2)λ22(a22aϖ+ϖ2)]

=

(λ221)[(c2+2cϖ+ϖ2)λ12(a2+2aϖ+ϖ2)]

(λ121)[c22cϖ+ϖ2a2λ22+2aϖλ22ϖ2λ22]

=

(λ221)[c2+2cϖ+ϖ2a2λ122aϖλ12ϖ2λ12]

(λ121)[(c2a2λ22)+(2aλ222c)ϖ+ϖ2(1λ22)]

=

(λ221)[(c2a2λ12)+(2c2aλ12)ϖ+ϖ2(1λ12)]

(λ121)(c2a2λ22)+ϖ2(1λ22)(λ121)+ϖ(2aλ222c)(λ121)

=

(λ221)(c2a2λ12)+ϖ(2c2aλ12)(λ221)+ϖ2(1λ12)(λ221).

ϖ2[(1λ22)(λ121)(1λ12)(λ221)]

=

ϖ[(2c2aλ12)(λ221)(2aλ222c)(λ121)]+[(λ221)(c2a2λ12)(λ121)(c2a2λ22)]

Now, because the coefficient of the ϖ2 term (on the left-hand-side of this last expression) is zero, we find that,

2ϖ

=

(λ121)(c2a2λ22)(λ221)(c2a2λ12)(caλ12)(λ221)(aλ22c)(λ121)

 

=

(λ12λ22)(c2a2)(λ12+λ22)(c+a)2c.

Relationship Between Key Angles

Law of Cosines:

(2c)2

=

R12+R22RR1cosψ

(2a)2

=

r12+r222r1r2cosθ

Try Again

Let's transform from Wong's lengths and angles to the ones used by Dyson. (This will permit us to compare Wong's analytic solution to the approximate one presented by Dyson.) Given the coordinate-pair (ϖ,z) and Dyson's length scale, c, we have,

tanχR

=

zcϖ,

R2

=

z2+(cϖ)2.

From this pair of expressions we can transform from (R,χR) to (ϖ,z) as follows:

cϖ

=

zcotχR

R2

=

z2(1+cotχR2)=[zsinχR]2

z

=

RsinχR

ϖ

=

ccotχRRsinχR=cRcosχR.

Next, given the coordinate-pair (ϖ,z) and Wong's length scale, a, we have,

tanχ2

=

zaϖ,

r22

=

z2+(aϖ)2;

and, combining the two sets gives,

r22

=

R2sin2χR+[RcosχR+(ac)]2

 

=

R2sin2χR+[R2cos2χR+2(ac)RcosχR+(ac)2]

 

=

R2+(ac)22(ca)RcosχR;

tanχ2

=

RsinχRRcosχR+(ac).

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |