User:Tohline/SSC/Structure/BiPolytropes

From VistrailsWiki
< User:Tohline
Revision as of 15:59, 24 March 2013 by Tohline (talk | contribs) (Start new page: BiPolytropes)
(diff) ← Older revision | Latest revision (diff) | Newer revision → (diff)
Jump to navigation Jump to search

BiPolytropes

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

As has been derived and discussed elsewhere, an isolated isothermal sphere has a density profile that extends to infinity and, correspondingly, an unbounded total mass. In an astrophysical context, neither of these properties is desirable. A more realistic isothermal configuration can be constructed by embedding the structure in a low density, but hot external medium whose pressure, <math>P_e</math>, confines the isothermal configuration to a finite size. In a mathematical model, this can be accomplished by ripping off an outer layer of the isolated isothermal configuration down to the radius — label it <math>\xi_e</math> — at which the configuration's original (internal) pressure equals <math>P_e</math>; the interior of the configuration that remains — containing mass <math>M_{\xi_e}</math> — should be unaltered and in equilibrium. (This will work only for spherically symmetric configurations, as the gravitational acceleration at any location only depends on the mass contained inside that radius.) Ebert (1955) and Bonnor (1956) are credited with constructing the first such models and, most significantly, discovering that, for any specified sound speed and applied external pressure, there is a mass above which no equilibrium configuration exists. We present, here, the salient elements of these (essentially equivalent) derivations.

Prior to studying this discussion of pressure-bounded isothermal spheres, we recommend studying our related discussion of pressure-bounded <math>~n</math> = 5 polytropes. As with isolated isothermal spheres, isolated <math>~n</math> = 5 polytropes extend to infinity. But, unlike their isothermal counterparts, the structure of <math>~n</math> = 5 polytropes is describable analytically. Hence, an analysis of their structure and its extension to pressure-bounded configurations avoids the clutter introduced by a model — such as the isothermal sphere — that can only be described numerically. As it turns out, the pressure-bounded <math>~n</math> = 5 polytrope exhibits a Bonnor-Ebert type limiting mass that is analytically prescribable. Its derivation is mathematically quite clean and provides a firm foundation for understanding the better known — but only numerically prescribable — Bonnor-Ebert limiting mass.

Governing Relation

The equilibrium structure of an isolated isothermal sphere, as derived by Emden (1907), has been discussed elsewhere. From this separate discussion we appreciate that the governing ODE is,

<math>\frac{1}{r^2} \frac{d}{dr}\biggl( r^2 \frac{d\ln\rho}{dr} \biggr) =- \frac{4\pi G}{c_s^2} \rho \, ,</math>

where,

<math>c_s^2 = \frac{\Re T}{\bar{\mu}} = \frac{k T}{m_u \bar{\mu}} \, ,</math>

is the square of the isothermal sound speed. In their studies of pressure-bounded isothermal spheres, Ebert (1955, ZA, 37, 217) and Bonnor (1956, MNRAS, 116, 351) both started with this governing ODE, but developed its solution in different ways. Here we present both developments while highlighting transformations between the two.

Derivation by Bonnor (edited) translation Derivation by Ebert (edited)
Bonnor (1956, MNRAS, 116, 351)
<math>G \Leftrightarrow \gamma</math>
Ebert (1955, ZA, 37, 217)
<math>\rho_c \Leftrightarrow \rho_0</math>
<math>\frac{kT}{m} \Leftarrow c_s^2 \Rightarrow \frac{\Re T_0}{\mu}</math>
<math>\beta^{1/2}\lambda^{-1/2} \Leftrightarrow l_0</math>
<math>e^{-\psi} \Leftrightarrow \eta</math>

Both of these dimensionless governing ODEs — Bonnor's Eq. (2.8) and Ebert's Eq. (17) — are identical to the dimensionless expression derived by Emden (see the presentation elsewhere), namely,

<math> \frac{d^2v_1}{d\mathfrak{r}_1^2} +\frac{2}{\mathfrak{r}_1} \frac{dv_1}{dr} + e^{v_1} = 0 \, . </math>

The translation from Emden-to-Bonnor-to-Ebert is straightforward:

<math> \mathfrak{r}_1 = \xi|_\mathrm{Bonner} = \xi|_\mathrm{Ebert}~~~~\mathrm{and}~~~~e^{v_1} = e^{-\psi} = \eta \, . </math>

Related Wikipedia Discussions

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