Difference between revisions of "User:Tohline/Appendix/Ramblings/BordeauxSequences"

From VistrailsWiki
Jump to navigation Jump to search
Line 124: Line 124:
</table>
</table>


Let's examine the "high-&omega; limit" that has been ''fit'' by their equation (3), noting that <math>~\omega^2</math> and <math>~\rho</math> are uniform throughout each two-component, equilibrium configuration.  We expect each (thin) ring component to be spinning with the Keplerian frequency,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\frac{r^2 \omega^2_\mathrm{K} }{r}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{GM_\mathrm{pt}}{r^2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \omega_K^2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{GM_\mathrm{pt}}{r^3}
=
\biggl(\frac{G}{r^3}\biggr) \frac{4\pi {\bar{a}}^3 \rho}{3}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>~\Rightarrow ~~~\frac{ \omega_K^2}{4\pi G\rho}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl(\frac{{\bar{a}}^3}{3r^3}\biggr) \, .
</math>
  </td>
</tr>
</table>
   
   
=====Central Object=====
=====Central Object=====
Line 380: Line 335:
   <td align="left">
   <td align="left">
<math>~  
<math>~  
2\pi^2 R_\mathrm{eq}^3 \rho ~
\biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} \biggl\{ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} -1
\biggr\}^2
\, .
</math>
  </td>
</tr>
</table>
=====Combined Configuration=====
The total (maximum) mass of the system is,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~M_\mathrm{tot}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~M_c + M_\mathrm{max}</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{4\pi}{3} \rho R_\mathrm{eq}^3 (1 - e^2)^{1 / 2}
+
2\pi^2 R_\mathrm{eq}^3 \rho ~  
2\pi^2 R_\mathrm{eq}^3 \rho ~  
\biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} \biggl\{ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} -1  
\biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} \biggl\{ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} -1  
Line 387: Line 379:
</tr>
</tr>
</table>
</table>




{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Revision as of 04:52, 23 November 2020

Université de Bordeaux (Part 2)

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

Exterior Gravitational Potential of Toroids

J. -M. Huré, B. Basillais, V. Karas, A. Trova, & O. Semerák (2020), MNRAS, 494, 5825-5838 have published a paper titled, The Exterior Gravitational Potential of Toroids. Here we examine how their work relates to the published work by C.-Y. Wong (1973, Annals of Physics, 77, 279), which we have separately discussed in detail.

We discuss this topic in a separate, accompanying chapter.

Spheroid-Ring Systems

Through a research collaboration at the Université de Bordeaux, B. Basillais & J. -M. Huré (2019), MNRAS, 487, 4504-4509 have published a paper titled, Rigidly Rotating, Incompressible Spheroid-Ring Systems: New Bifurcations, Critical Rotations, and Degenerate States.

Key References

Here are some relevant publications:

Especially,

Key Figures

Eriguchi & Sugimoto (1981)

Fig. 1 extracted without modification from p. 1873 of Eriguchi & Sugimoto (1981)

"Another Equilibrium Sequence of Self-Gravitating and Rotating Incompressible Fluid"

Progress of Theoretical Physics,

vol. 65, pp. 1870-1875 © Progress of Theoretical Physics

Figure 3 from Eriguchi & Hachisu (1983)

CAPTION (modified here):   The squared angular velocity is plotted against <math>~j^2</math> for a segment of the Maclaurin sequence (dashed curve), for the Dyson-Wong sequence (dotted curve), and for the new configurations reported in this 1981 paper by Eriguchi & Sugimoto (solid curve). The "×" mark denotes the neutral point on the Maclaurin sequence against the <math>~P_4(\eta)</math> perturbation. The dotted curve is plotted by using the values which are read from the curve of Fig. 6 of Wong (1974), so it may contain errors to some extent.


Eriguchi & Hachisu (1983)

Fig. 3 extracted without modification from p. 1134 of Eriguchi & Hachisu (1983)

"Two Kinds of Axially Symmetric Equilibrium Sequences of Self-Gravitating and Rotating Incompressible Fluids:
Two-Ring Sequence and Core-Ring Sequence
"

Progress of Theoretical Physics,

vol. 69, pp. 1131-1136 © Progress of Theoretical Physics

Figure 3 from Eriguchi & Hachisu (1983)

CAPTION:  The angular momentum-angular velocity relations. Solid curves represent uniformly rotating equilibrium sequences.

  • MS:   Maclaurin spheroid sequence
  • JE:   Jacobi ellipsoid sequence
  • OR:   one-ring sequence

The number and letter R or C attached to a curve denote mass ratio and two-ring or core-ring sequence, respectively. If differential rotation is allowed, the equilibrium sequences may continue to exist as shown by the dashed curves.

AKM (2003)

Fig. 2 extracted without modification from p. 517 of Ansorg, Kleinwächter & Meinel (2003)

"Uniformly rotating axisymmetric fluid configurations bifurcating from highly flattened Maclaurin spheroids"

MNRAS, vol. 339, pp. 515-523 © Royal Astronomical Society

Figure 2 from Ansorg, Kleinwächter & Meinel (2003)

CAPTION:   For the first five axisymmetric sequences, <math>~\omega_0^2</math> is plotted against the dimensionless squared angular momentum, <math>~j^2</math>, using the same normalizations as Eriguchi & Hachisu (1983). Dotted and dashed curves again refer to the Maclaurin sequence and the Dyson approximation respectively. The full circles mark the bifurcation points on the Maclaurin sequence, and the open square the transition configuration of spheroidal to toroidal bodies on the Dyson ring sequence.

Basillais & Huré (2019)

Fig. 4 extracted without modification from p. 4507 of Basillais & Huré (2019)

"Rigidly rotating, incompressible spheroid-ring systems:   new bifurcations, critical rotations, and degenerate states"

MNRAS, vol. 487, pp. 4504-4509 © Royal Astronomical Society

Figure 4 from Basillais & Huré (2019)

CAPTION:   The spheroid-ring solutions (grey dots) populate the <math>~\omega_0^2 - j^2</math> diagram in between the MLS, the high-ω limit, and the high-j limit. The MLS, ORS, Jacobi sequence, Hamburger sequence, and ε2-sequence are also shown (plain lines). Points labelled a to f (cross) correspond to equilibria shown in Figure 3; see also Table 1. There is a band of degeneracy rightward to the ORS (green dashed zone).


Central Object

Assume that the central object is exactly a Maclaurin spheroid. Then from Figure 1 (and Table 1) of our review of equilibrium models along the Maclaurin spheroid sequence, we appreciate that all we have to do is specify the eccentricity, <math>~0 \le e \le 1</math>, and <math>~\Omega^2 \equiv \omega_0^2/(4\pi G\rho)</math> is known. For example, if we choose <math>~e = 0.60</math>, then from that Table 1, <math>~\omega_0^2/(2\pi G\rho) = 0.1007~~\Rightarrow~~\Omega^2 \approx 0.05</math>. Other properties of this "central" spheroid — such as its mass, moment of inertia, and angular momentum — are given the the following expressions:

<math>~M_c</math>

<math>~=</math>

<math>~\frac{4\pi}{3} \rho R_\mathrm{eq}^2 Z = \frac{4\pi}{3} \rho R_\mathrm{eq}^3 (1 - e^2)^{1 / 2} \, ;</math>

<math>~I_c</math>

<math>~=</math>

<math>~\frac{2}{5} M_c R_\mathrm{eq}^2 \, ;</math>

<math>~L_c</math>

<math>~=</math>

<math>~I_c \omega_0 = \frac{2}{5} M_c R_\mathrm{eq}^2 \omega_0 </math>

 

<math>~=</math>

<math>~\frac{2}{5} \biggl[ \frac{4\pi}{3} \rho R_\mathrm{eq}^3 (1 - e^2)^{1 / 2} \biggr] R_\mathrm{eq}^2 \omega_0 </math>

 

<math>~=</math>

<math>~\frac{2^3\pi}{3\cdot 5} (1 - e^2)^{1 / 2} R_\mathrm{eq}^5 \rho \omega_0 </math>

<math>~\Rightarrow ~~~L^2_c</math>

<math>~=</math>

<math>~4\pi G\biggl( \frac{2^3\pi}{3\cdot 5}\biggr)^2 (1 - e^2) R_\mathrm{eq}^{10} \rho^3 \Omega^2 \, . </math>

We note as well that the Keplerian frequency for a massless particle orbiting in the equatorial plane of this central object will be,

<math>~\omega_K^2</math>

<math>~=</math>

<math>~\frac{GM_c}{r^3}</math>

 

<math>~=</math>

<math>~\frac{G}{r^3} \biggl[ \frac{4\pi}{3} \rho R_\mathrm{eq}^3 (1 - e^2)^{1 / 2} \biggr] \, .</math>

So, if we force this orbital frequency to also equal the spin-frequency of the Maclaurin spheroid, the radius of the orbit must be,

<math>~\Omega^2 (4\pi G\rho)</math>

<math>~=</math>

<math>~\frac{G}{r^3} \biggl[ \frac{4\pi}{3} \rho R_\mathrm{eq}^3 (1 - e^2)^{1 / 2} \biggr] </math>

<math>~\Rightarrow~~~\Omega^2 </math>

<math>~=</math>

<math>~\frac{1}{r^3} \biggl[ \frac{1}{3} R_\mathrm{eq}^3 (1 - e^2)^{1 / 2} \biggr] </math>

<math>~\Rightarrow~~~\frac{r}{R_\mathrm{eq}} </math>

<math>~=</math>

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

For example, when <math>~(e, \Omega^2) = (0.6, 0.05)</math>, we have, <math>~r/R_\mathrm{eq} = (5.3333)^{1 / 3} = 1.747</math>.


Surrounding Torus

We'll assume that the surrounding 2nd object is a thin torus (1) with the same density as the central object, (2) with a major axis, <math>~a</math>, which ensures that the torus is spinning with the Keplerian frequency prescribed by the mass of the central object, (3) and with a minor cross-sectional radius, <math>~b</math>. The second of these constraints means that,

<math>~\frac{a}{R_\mathrm{eq}} </math>

<math>~=</math>

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

The mass of the torus is given by the expression,

<math>~M_t</math>

<math>~=</math>

<math>~2\pi a (\pi b^2)\rho = 2\pi^2 R_\mathrm{eq}^3 \rho \biggl[ \frac{a}{R_\mathrm{eq}} \biggl( \frac{b}{R_\mathrm{eq}} \biggr)^2 \biggr] \, .</math>

Given that the major radius of the torus, <math>~a</math>, is known once the eccentricity of the central Maclaurin spheroid has been selected, and the density of the torus must match the density of the central object, the mass of the torus will only depend on the choice of <math>~b</math>. The maximum allowed value of <math>~b</math> is set by ensuring that the torus lies entirely outside of the central spheroid. This means,

<math>~b_\mathrm{max}</math>

<math>~=</math>

<math>~R_\mathrm{eq}\biggl[\frac{a}{R_\mathrm{eq}}-1 \biggr]</math>

So, the maximum torus mass is,

<math>~M_\mathrm{max}</math>

<math>~=</math>

<math>~ 2\pi^2 R_\mathrm{eq}^3 \rho ~ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} \biggl\{ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} -1 \biggr\}^2 \, . </math>

Combined Configuration

The total (maximum) mass of the system is,

<math>~M_\mathrm{tot}</math>

<math>~=</math>

<math>~M_c + M_\mathrm{max}</math>

 

<math>~=</math>

<math>~ \frac{4\pi}{3} \rho R_\mathrm{eq}^3 (1 - e^2)^{1 / 2} + 2\pi^2 R_\mathrm{eq}^3 \rho ~ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} \biggl\{ \biggl[\frac{1}{3\Omega^2} (1 - e^2)^{1 / 2} \biggr]^{1 / 3} -1 \biggr\}^2 \, . </math>


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