Difference between revisions of "User:Tohline/Apps/DysonPotential"

From VistrailsWiki
Jump to navigation Jump to search
Line 6: Line 6:


==Overview==
==Overview==
In his pioneering work, [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893, Philosophical Transactions of the Royal Society of London. A., 184, 43 - 95)] and [http://adsabs.harvard.edu/abs/1893RSPTA.184.1041D (1893, 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.  [http://adsabs.harvard.edu/abs/1974ApJ...190..675W 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, [http://adsabs.harvard.edu/abs/1981PThPh..65.1870E Y. Eriguchi & D. Sugimoto (1981, Progress of Theoretical Physics, 65, 1870 - 1875)] and [http://adsabs.harvard.edu/abs/1988ApJS...66..315H 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.
In his pioneering work, [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893a, Philosophical Transactions of the Royal Society of London. A., 184, 43 - 95)] and [http://adsabs.harvard.edu/abs/1893RSPTA.184.1041D (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.  [http://adsabs.harvard.edu/abs/1974ApJ...190..675W 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, [http://adsabs.harvard.edu/abs/1981PThPh..65.1870E Y. Eriguchi & D. Sugimoto (1981, Progress of Theoretical Physics, 65, 1870 - 1875)] and [http://adsabs.harvard.edu/abs/1988ApJS...66..315H 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===
===External Potential===
On p. 62 of [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893, Philosophical Transactions of the Royal Society of London. A., 184, 43 - 95)], Dyson presents the following approximate expression for the potential everywhere exterior to an ''anchor ring'':
On p. 62 of [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893a)], we find the following approximate expression for the potential at point "P", anywhere exterior to an [http://www.mathematicsdictionary.com/english/vmd/full/t/torusanchorring.htm anchor ring]:
<table border="1" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
<tr><td align="center" colspan="1">
<tr><td align="center" colspan="1">
'''Equation image extracted without modification from p. 62 of [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893)]'''<p></p>
'''Equation image extracted without modification from p. 62 of [http://adsabs.harvard.edu/abs/1893RSPTA.184...43D F. W. Dyson (1893)]'''<p></p>
Line 20: Line 20:
</td>
</td>
</tr>
</tr>
</table>
<table border="0" cellpadding="10" align="right" width="40%"><tr><td align="center">
<table border="1" cellpadding="5" align="right" width="100%">
<tr>
<tr>
<td align="center">
<td align="left">
''Anchor ring'' schematic, adapted from figure near the top of &sect;2 (on p. 47) of Dyson (1893a)<br />
[http://www.mathematicsdictionary.com/english/vmd/full/t/torusanchorring.htm Anchor ring] schematic, adapted from figure near the top of &sect;2 (on p. 47) of Dyson (1893a)<br />
[[File:DysonTorusIllustration02.png|450px|center|Anchor Ring Schematic]]
[[File:DysonTorusIllustration02.png|300px|center|Anchor Ring Schematic]]
</td>
</td>
</tr>
</tr>
</table>
</table>
</td></tr></table>
In Dyson's expression, the leading factor of <math>~F</math> is the [https://en.wikipedia.org/wiki/Elliptic_integral#Complete_elliptic_integral_of_the_first_kind complete elliptic integral of the first kind], namely,
In Dyson's expression, the leading factor of <math>~F</math> is the [https://en.wikipedia.org/wiki/Elliptic_integral#Complete_elliptic_integral_of_the_first_kind complete elliptic integral of the first kind], namely,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">

Revision as of 03:06, 28 August 2018

Dyson (1893)

Whitworth's (1981) Isothermal Free-Energy Surface
|   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 Dyson-Wong tori, beginning from a bifurcation point on the Maclaurin spheroid sequence.

External Potential

On p. 62 of F. W. 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 F. W. Dyson (1893)

The Potential of an Anchor Ring, Phil. Trans. Royal Soc. London. A., Vol. 184

The Potential Exterior to an Anchor Ring

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

Anchor Ring Schematic

In Dyson's expression, the leading factor of <math>~F</math> is the complete elliptic integral of the first kind, namely,

<math>~F = F(\mu)</math>

<math>~\equiv</math>

<math>~\int_0^{\pi/2} \frac{d\phi}{\sqrt{1 - \mu^2 \sin^2\phi}} \, ,</math>

where, <math>~\mu \equiv (R_1 - R)/(R_1 + R)</math>. Similarly, <math>~E = E(\mu)</math> is the complete elliptic integral of the second kind. In the limit of <math>~a/c \rightarrow 0</math>, Dyson's expression gives,

<math>~V_\mathrm{Dyson}</math>

<math>~=</math>

<math>~\frac{4K(\mu)}{R+R_1} \, ,</math>

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

<math>~\biggl[ \frac{\pi}{GM}\biggr] \Phi_\mathrm{TR}</math>

<math>~=</math>

<math>~- \frac{2K(k)}{R_1} \, ,</math>

where, <math>~k \equiv [1-(R/R_1)^2]^{1 / 2}</math>. Is Dyson's expression identical to this one when <math>~a/c = 0</math> ?


If we adopt the variable notation,

<math>~e^\eta \equiv \frac{R_1}{R} \, ,</math>

then we can write,

<math>~\cosh\eta = \frac{1}{2}\biggl[e^\eta + e^{-\eta}\biggr]</math>

<math>~=</math>

<math>~\frac{R^2 + R_1^2}{2RR_1} \, ,</math>

which implies that,

<math>~\biggl[ \frac{2}{\coth\eta +1} \biggr]^{1 / 2} = [1 - e^{-2\eta}]^{1 / 2}</math>

<math>~=</math>

<math>~\biggl[ 1 - \biggl(\frac{R}{R_1}\biggr)^2 \biggr]^{1 / 2} \, .</math>

This is the definition of the parameter, <math>~k</math>, in the expression for <math>~\Phi_\mathrm{TR}</math>. Now, if we employ the Descending Landen Transformation for the complete elliptic integral of the first kind, we can make the substitution,

<math>~K(k)</math>

<math>~=</math>

<math>~ (1 + k_1)K(k_1) \, , </math>

      where,      

<math>~k_1</math>

<math>~\equiv</math>

<math>~ \frac{1-\sqrt{1-k^2}}{1+\sqrt{1-k^2}} \, . </math>

But notice that, <math>~\sqrt{1-k^2} = e^{-\eta}</math>, in which case,

<math>~k_1 </math>

<math>~=</math>

<math>~ \frac{1-e^{-\eta}}{1+e^{-\eta}} </math>

<math>~=</math>

<math>~ \frac{1-R/R_1}{1+R/R_1} </math>

<math>~=</math>

<math>~ \frac{R_1-R}{R_1+R} \, , </math>

which is the definition of the parameter, <math>~\mu</math>, in the expression for <math>~V_\mathrm{Dyson}</math>. Hence, we can write,

<math>~\biggl[ \frac{\pi}{GM}\biggr] \Phi_\mathrm{TR}</math>

<math>~=</math>

<math>~- \frac{2}{R_1} \biggl[(1+k_1)K(k_1) \biggr] </math>

 

<math>~=</math>

<math>~- \frac{2K(\mu)}{R_1} \biggl[1+\frac{R_1-R}{R_1+R} \biggr] </math>

 

<math>~=</math>

<math>~- \frac{4K(\mu)}{R_1+R} \, ,</math>

which, aside from the adopted sign convention, is indeed precisely the expression given by <math>~V_\mathrm{Dyson}</math> when <math>~a/c = 0</math> .


Figures 1 - 6 extracted without modification from pp. 63-66 of F. W. Dyson (1893)

The Potential of an Anchor Ring, Phil. Trans. Royal Soc. London. A., Vol. 184

The Potential Exterior to an Anchor Ring; R/d = infinity
The Potential Exterior to an Anchor Ring; R/d = 5
The Potential Exterior to an Anchor Ring; R/d = 2.5
The Potential Exterior to an Anchor Ring; R/d = 1.667
The Potential Exterior to an Anchor Ring; R/d = 1.25
The Potential Exterior to an Anchor Ring; R/d = 1

See Also

The following quotes have been taken from Petroff & Horatschek (2008):

§1:   "The problem of the self-gravitating 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 <math>~\sigma = a/b</math>, where <math>~a = r_t</math> provides a measure for the radius of the cross-section of the ring and <math>~b = \varpi_t</math> the distance of the cross-section'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."

  1. 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) …"


 

Whitworth's (1981) Isothermal Free-Energy Surface

© 2014 - 2021 by Joel E. Tohline
|   H_Book Home   |   YouTube   |
Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS |
Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation