User:Tohline/Apps/DysonPotential
From VisTrailsWiki
Contents

Dyson (1893)
 Tiled Menu  Tables of Content  Banner Video  Tohline Home Page  
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 DysonWong tori, beginning from a bifurcation point on the Maclaurin spheroid sequence.
External Potential
His Derived Expression
(See an accompanying Ramblings Chapter for additional derivation details.) On p. 62, in §8 of Dyson (1893a), we find the following approximate expression for the potential at point "P", anywhere exterior to an anchor ring:

Equation image extracted without modification from p. 62 of Dyson (1893a)
The Potential of an Anchor Ring, Phil. Trans. Royal Soc. London. A., Vol. 184 
In Dyson's expression, the leading factor of is the complete elliptic integral of the first kind, namely,



where, . Similarly, is the complete elliptic integral of the second kind.
Comparison With Thin Ring Approximation
In the limit of , Dyson's expression gives,



where we have acknowledged that, in the twentyfirst century, the complete elliptic integral of the first kind is more customarily represented by the letter, . In a separate discussion, we have shown that the gravitational potential of an infinitesimally thin ring is given precisely by the expression,



where, . Is Dyson's expression identical to this one when ?
Proof
Taking a queue from our accompanying discussion of toroidal coordinates, if we adopt the variable notation,
then we can write,



which implies that,



This is the definition of the parameter, , in the expression for . Now, if we employ the Descending Landen Transformation for the complete elliptic integral of the first kind, we can make the substitution,



where, 



But notice that, , in which case,







which is the definition of the parameter, , in the expression for . Hence, we can write,









Aside from the adopted sign convention, this is indeed precisely the expression given by when .
Evaluation
Dyson's Figures
In his effort to illustrate the behavior of equipotential contours in the space exterior to various anchor rings, Dyson evaluated his expression for the potential up through ; that is, he evaluated the function,



Figures 1  6 from Dyson (1893a) — replicated immediately below — show his resulting set of contours for six cases: Tori (anchor rings) having aspect ratios of . Click on an image to view the contour plot at higher resolution. In what follows we present results from our own evaluation of this "V_{2}" function for the single case of an anchor ring having .

Our Attempt to Replicate
First, let's test the accuracy of Dyson's (1893a) "series expansion" expression for the elliptic integrals, and ; in the following table, the highprecision evaluations labeled "Numerical Recipes" have been drawn from the tabulated data that is provided in our accompanying discussion of incomplete elliptic integrals. According to, for example, Wikipedia, the relevant seriesexpansion expressions are:






These expressions — up through — can be found in the middle of p. 58 of Dyson (1893a). We strongly suspect that, in constructing the equipotential contours shown in his figures 16, Dyson used expressions for and that were more accurate than this. For example, we found it necessary to include terms up through in order to match to three digits accuracy the potential contour values and coordinate locations reported by Dyson.
Numerical Recipes  Series expansion up through  Series expansion up through  
^{†}  
0.34202014  1.62002589  1.52379921  1.6198  1.5239  1.6200263  1.5237989 
0.57357644  1.73124518  1.43229097  1.7239  1.4336  1.73124518  1.43230 
0.76604444  1.93558110  1.30553909  1.8773  1.3150  1.93558109  1.3061 
0.90630779  2.30878680  1.16382796  2.042  1.199  2.308784  1.1700 
0.98480775  3.15338525  1.04011440  2.16  1.12  3.150  1.069 
^{†}We actually used the "descending Landen transformation" to evaluate through . 
For and a specification of the ratio, , take the following steps to map out an equipotential curve that has :
 Choose a value of
 Guess a value of and,
 Set
 Evaluate the function,
 If to the desired accuracy, loop back up and guess another value of
 If to the desired accuracy, save the coordinate location, , and loop back up to pick another value of
Tabulated Data
As the data in the following table documents, we have been able to construct equipotential contours that agree with Dyson, not only qualitatively, but quantitatively. For example:
 The dark green contour has been designed to touch the surface of the torus precisely where its outermost edge cuts through the equatorial plane . This means that and . (These four coordinate values are highlighted in pink in the second major column of the table.) When we plugged these values of and into Dyson's expression for , we determined that the value of the potential at this point on the torus surface is 0.8551 — see the yellowhighlighted heading of the second major table column. Compare this to the value of 0.855 that Dyson has printed just below the Figure 3 xaxis where a fiducial identifies the coordinate, . As has been catalogued at the bottom of table column #2, we have found that this darkgreen contour touches the vertical axis at the coordinate location, , for which, .
 By design — see the coordinate values highlighted in pink in table column #1 — our outermost (pink) contour touches the equatorial plane at . When we plugged these values of and into Dyson's expression for , we determined that the value of the potential at this point outside the torus is 0.7737 — see the yellowhighlighted heading of table column #1. Compare this to the value of 0.777 that Dyson has printed just below the Figure 3 xaxis where a fiducial identifies the coordinate, . As has been catalogued at the bottom of table column #1, we have found that this pink contour touches the vertical axis at the coordinate location, , for which, .
 Similarly, we have constructed contours that intersect the equatorial plane at the fiducials marking (red curve & table column #5), (lightgreen curve & table column #6), and (lightblue curve & table column #7). According to our calculations, they correspond, respectively, to values of the potential, (Dyson's corresponding fiducial label is 0.980), (Dyson's corresponding fiducial label is 0.990), and (Dyson's corresponding fiducial label is 1.021).
 Finally, we constructed two contours (blue and orange) by initially specifying the value of the potential, rather than specifying the coordinate values . We used the values of the potential that Dyson associated with the fiducials along the vertical axis at and at : Respectively, — blue contour detailed in our table column #3 — and — orange contour detailed in our table column #4. We determined that these two contour curves intersected the vertical axis at, respectively, and , that is, at coordinate locations that were nearly identical to the locations labeled by Dyson.
Coordinates of Points that Trace Seven Different Equipotential Contours External to the Anchor Ring With 


Column #1  Column #2  Column #3  Column #4  Column #5  Column #6  Column #7  







Intermediate Step
Objective
As has been reprinted above, on p. 62 of Dyson's Part I we find his powerseries expression for the external potential, namely,






where — as in the context of toroidal coordinates — we occasionally will make the substitution, , and therefore,



In order to facilitate matching boundary conditions at the surface of the torus, between the exterior and interior expressions for the gravitational potential, Dyson rewrites this Part I expression for the external potential and — explicitly evaluating it on the torus surface — sets, . Specifically, on p. 1049 of Dyson's Part II we find equation (6), which reads,









In order to obtain this alternate powerseries expression, Dyson …
 Expresses angular variations in terms of the angle, , instead of the angle, ; these two angles are identified in the above schematic.
 Employs powerseries expansions of both elliptic integral functions, and .
 Uses the binomial theorem to develop a number of other powerseries expressions.
In what follows we will attempt to demonstrate that this second (Part II, equation 6) expression is identical to the first.
The Ratio R_{1}/c
Note that, via the law of cosines,






At the surface of the torus, where , we therefore have,



Low Order
Employing the binomial theorem, we can write,





















Higher Order
Adopting the shorthand notation,
and
and employing the binomial theorem, we can write,
























Also, we have,



























Hence,



And, adopting the shorthand notation,
we have,













































Relationship Between Angles
Drawing on the Law of Cosines, as above, we can state that on the torus surface,



Alternatively, applying the Law of Cosines to the angle, , we have,






Therefore, anywhere along the surface of the torus, we can switch from one of these angles to the other via the relation,






Cosine ψ Expansion
Employing the binomial theorem, we therefore can write,
























CosineSquared Expansion
Letting,
via the binomial theorem we have,

























































CosineCubed Expansion
Again, letting,
via the binomial theorem we have,




























































Coefficients of Elliptic Integrals
Rewriting the external potential, as provided in the abovestated objective, and evaluating it at the torus surface,






where,



and,



Given our derived powerseries expressions for various trigonometric functions, these coefficients can be rewritten as,

































and,



























Now, inserting to the appropriate order the above expression for the ratio, — namely,



we have,


















and,


















Alternate "Small" Argument of Elliptic Integrals
Defining the "small parameter,"


















At the surface of the torus, where , we therefore have,





















Alternatively, if — as above — we adopt the shorthand notation, , we can write,
Summary



















where,


















Now, if we adopt the shorthand notation,
we also have,




































And,





















And,


















Elliptic Integral Expressions
Hence, drawing upon the series expansion for the complete elliptic integral of the first kind that can be found, for example, as eq. (8.113.3) in the Fourth Edition of Gradshteyn & Ryzhik (1965), we can write,









Now, we recognize that,









Also, drawing upon the series expansion for the complete elliptic integral of the second kind that can be found, for example, as eq. (8.114.3) in the Fourth Edition of Gradshteyn & Ryzhik (1965), we can write,









External Potential at Torus Surface
Initial Low Resolution
Hence,












Hence,












To order , this gives,


















We are trying to match equation (6) in Dyson's (1893b) "Part II", that is,









High Resolution
Summary



























































































That is,

























































Insert Expressions for K and E
Summary

Remember that (see above),



And (see above),



We are trying to match equation (6) in Dyson's (1893b) "Part II", that is,









To First Order












Hence,



























To Second Order
























where,



Hence,


































































In an effort to compare this expression with equation (6) from Dyson's (1893b) "Part II", we should make the substitutions,
and
This means,















This expression differs from the 2^{nd}order term in Dyson's equation (6) by the amount,


















Interior Potential
In equation (9) on p. 1050 of his "Part II", Dyson (1893b) presents the following powerseries expression for the gravitational potential at points inside the torus:






where,



Note that, for the example illustrated above, and, hence, . Therefore, at any point on the surface of this example torus,












Red Contour
Let's construct an equipotential contour that extends the red contour into the interior region. Let's begin by evaluating Dyson's interior potential expression at the coordinate location where the red contour touches the surface of the torus. According to Column #5, this point on the surface has coordinates, ; or, equivalently, and,






Hence, for this specific point on the torus surface, we find,



Contour (column #) 
(external)  Surface  (radians) 

blue (3)  0.9120  1.117  1.8676 
orange (4)  0.961  0.929  1.392 
red (5)  0.9800  0.867  1.2318 
light green (6)  0.9800  0.838  1.1538 
light blue (7)  1.0212  0.749  0.8925 
See Also
 T. Fukushima (2016, AJ, 152, article id. 35, 31 pp.) — Zonal Toroidal Harmonic Expansions of External Gravitational Fields for Ringlike Objects
 W.T. Kim & S. Moon (2016, ApJ, 829, article id. 45, 22 pp.) — Equilibrium Sequences and Gravitational Instability of Rotating Isothermal Rings
 E. Y. Bannikova, V. G. Vakulik & V. M. Shulga (2011, MNRAS, 411, 557  564) — Gravitational Potential of a Homogeneous Circular Torus: a New Approach
 D. Petroff & S. Horatschek (2008, MNRAS, 389,156  172) — Uniformly Rotating Homogeneous and Polytropic Rings in Newtonian Gravity
The following quotes have been taken from Petroff & Horatschek (2008):
§1: "The problem of the selfgravitating ring captured the interest of such renowned scientists as Kowalewsky (1885), Poincaré (1885a,b,c) and Dyson (1892, 1893). Each of them tackled the problem of an axially symmetric, homogeneous ring in equilibrium by expanding it about the thin ring limit. In particular, Dyson provided a solution to fourth order in the parameter , where provides a measure for the radius of the crosssection of the ring and the distance of the crosssection's centre of mass from the axis of rotation."
§7: "In their work on homogeneous rings, Poincaré and Kowalewsky, whose results disagreed to first order, both had made mistakes as Dyson has shown. His result to fourth order is also erroneous as we point out in Appendix B." 
 P. H. Chavanis (2006, International Journal of Modern Physics B, 20, 3113  3198) — Phase Transitions in SelfGravitating Systems
 M. Ansorg, A. Kleinwächter & R. Meinel (2003, MNRAS, 339, 515) — Uniformly Rotating Axisymmetric Fluid Configurations Bifurcating from Highly Flattened Maclaurin Spheroids
 M. Lombardi & G. Bertin (2001, Astronomy & Astrophysics, 375, 1091  1099) — Boyle's Law and Gravitational Instability
 W. Kley (1996, MNRAS, 282, 234) — Maclurin Discs and Bifurcations to Rings
 J. W. Woodward, J. E. Tohline, & I. Hachisu (1994, ApJ, 420, 247  267) — The Stability of Thick, SelfGravitating Disks in Protostellar Systems
 I. Bonnell & P. Bastien (1991, ApJ, 374, 610  622) — The Collapse of Cylindrical Isothermal and Polytropic Clouds with Rotation
 J. E. Tohline & I. Hachisu (1990, ApJ, 361, 394  407) — The Breakup of SelfGravitating Rings, Tori, and Thick Accretion Disks
 F. Schmitz (1988, Astronomy & Astrophysics, 200, 127  134) — Equilibrium Structures of Differentially Rotating SelfGravitating Gases
 P. Veugelen (1985, Astrophysics & Space Science, 109, 45  55) — Equilibrium Models of Differentially Rotating Polytropic Cylinders
 M. A. Abramowicz, A. Curir, A. SchwarzenbergCzerny, & R. E. Wilson (1984, MNRAS, 208, 279  291) — SelfGravity and the Global Structure of Accretion Discs
 P. Bastien (1983, Astronomy & Astrophysics, 119, 109  116) — Gravitational Collapse and Fragmentation of Isothermal, NonRotating, Cylindrical Clouds
 Y. Eriguchi & D. Sugimoto (1981, Progress of Theoretical Physics, 65, 1870  1875) — Another Equilibrium Sequence of SelfGravitating and Rotating Incompressible Fluid
 J. E. Tohline (1980, ApJ, 236, 160  171) — Ring Formation in Rotating Protostellar Clouds
 T. Fukushima, Y. Eriguchi, D. Sugimoto, & G. S. BisnovatyiKogan (1980, Progress of Theoretical Physics, 63, 1957  1970) — Concave Hamburger Equilibrium of Rotating Bodies
 J. Katz & D. LyndenBell (1978, MNRAS, 184, 709  712) — The Gravothermal Instability in Two Dimensions
 P. S. Marcus, W. H. Press, & S. A. Teukolsky (1977, ApJ, 214, 584 597) — Stablest Shapes for an Axisymmetric Body of Gravitating, Incompressible Fluid (includes torus with nonuniform rotation)
 Shortly after their equation (3.2), Marcus, Press & Teukolsky make the following statement: "… we know that an equilibrium incompressible configuration must rotate uniformly on cylinders (the famous "PoincaréWavre" theorem, cf. Tassoul 1977, &Sect;4.3) …"
 C. J. Hansen, M. L. Aizenman, & R. L. Ross (1976, ApJ, 207, 736  744) — The Equilibrium and Stability of Uniformly Rotating, Isothermal Gas Cylinders
 C.Y. Wong (1974, ApJ, 190, 675  694) — Toroidal Figures of Equilibrium
 C.Y. Wong (1973, Annals of Physics, 77, 279  353) — Toroidal and Spherical Bubble Nuclei
 J. Ostriker (1964, ApJ, 140, 1056) — The Equilibrium of Polytropic and Isothermal Cylinders
 J. Ostriker (1964, ApJ, 140, 1067) — The Equilibrium of SelfGravitating Rings
 J. Ostriker (1964, ApJ, 140, 1529) — On the Oscillations and the Stability of a Homogeneous Compressible Cylinder
 J. Ostriker (1965, ApJ Supplements, 11, 167) — Cylindrical Emden and Associated Functions
 Gunnar Randers (1942, ApJ, 95, 88) — The Equilibrium and Stability of RingShaped 'barred SPIRALS'.
 William Duncan MacMillan (1958), The Theory of the Potential, New York: Dover Publications
 Oliver Dimon Kellogg (1929), Foundations of Potential Theory, Berlin: Verlag Von Julius Springer
 Lord Rayleigh (1917, Proc. Royal Society of London. Series A, 93, 148154) — On the Dynamics of Revolving Fluids
 F. W. Dyson (1893, Philosophical Transaction of the Royal Society London. A., 184, 1041  1106) — The Potential of an Anchor Ring. Part II.
 In this paper, Dyson derives the gravitational potential inside the ring mass distribution
 F. W. Dyson (1893, Philosophical Transaction of the Royal Society London. A., 184, 43  95) — The Potential of an Anchor Ring. Part I.
 In this paper, Dyson derives the gravitational potential exterior to the ring mass distribution
 S. Kowalewsky (1885, Astronomische Nachrichten, 111, 37) — Zusätze und Bemerkungen zu Laplace's Untersuchung über die Gestalt der Saturnsringe
 Poincaré (1885a, C. R. Acad. Sci., 100, 346), (1885b, Bull. Astr., 2, 109), (1885c, Bull. Astr. 2, 405). — references copied from paper by Wong (1974)
© 2014  2019 by Joel E. Tohline 