Difference between revisions of "User:Tohline/SSC/Structure/IsothermalSphere"
(→Isothermal Sphere (structure): Add tag to Emden's Table 14) |
(→Emden Numerical Solution: Insert plot to illustrate Emden's solution) |
||
Line 84: | Line 84: | ||
==Emden Numerical Solution== | ==Emden Numerical Solution== | ||
<div align="center"> | <div align="center"> | ||
<table border="0" cellpadding=" | <table border="0" cellpadding="1"> | ||
<tr> | |||
<td> | |||
| |||
</td> | |||
<td align="center" valign="top" rowspan="3"> | |||
[[File:EmdenTable14.jpg|500px|center|Emden's (1907) Table 14]] | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="center" valign="top"> | |||
[[File:IsothermalDensityPlot.jpg|350px|center|Plotted from Emden's (1907) tabulated data]] | |||
</td> | |||
</tr> | |||
<tr> | |||
<td valign="top"> | |||
The set of <math>(x,y)</math> data points plotted in blue are drawn directly from column 1 — actually, <math>x = \ln(\mathfrak{r}_1)</math> — and column 3 — specifically, <math>y = v_1</math> — of Emden's Table 14. The dashed red line has a slope of <math>-2</math> and serves to illustrate that, at large radii, the isothermal density profile tends toward a <math>\rho \propto r^{-2}</math> distribution. | |||
</td> | |||
</tr> | |||
<tr> | <tr> | ||
<td align="center"> | <td align="center"> | ||
| |||
</td> | |||
<td align="center"> | |||
<font size="-1"> | |||
Note: The <math>3^\mathrm{rd}</math> column table entry highlighted in blue must be a typesetting error. | |||
</font> | |||
</td> | </td> | ||
</tr> | </tr> |
Revision as of 22:30, 29 October 2012
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Isothermal Sphere (structure)
Here we supplement the simplified set of principal governing equations with an isothermal equation of state, that is, <math>~P</math> is related to <math>~\rho</math> through the relation,
<math>P = c_s^2 \rho \, ,</math>
where, <math>c_s</math> is the isothermal sound speed. Comparing this <math>~P</math>-<math>~\rho</math> relationship to
Form A
of the Ideal Gas Equation of State,
<math>~P_\mathrm{gas} = \frac{\Re}{\bar{\mu}} \rho T</math> |
we see that,
<math>c_s^2 = \frac{\Re T}{\bar{\mu}} = \frac{k T}{m_u \bar{\mu}} \, ,</math>
where, <math>~\Re</math>, <math>~k</math>, <math>~m_u</math>, and <math>~\bar{\mu}</math> are all defined in the accompanying variables appendix. It will be useful to note that, for an isothermal gas, <math>~H</math> is related to <math>~\rho</math> via the expression,
<math> dH = \frac{dP}{\rho} = c_s^2 d\ln\rho \, . </math>
Governing Relations
Adopting solution technique #2, we need to solve the following second-order ODE relating the two unknown functions, <math>~\rho</math> and <math>~H</math>:
<math>\frac{1}{r^2} \frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =- 4\pi G \rho</math> .
Using the <math>~H</math>-<math>~\rho</math> relationship for an isothermal gas presented above, this can be rewritten entirely in terms of the density as,
<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>
or, equivalently,
<math> \frac{d^2\ln\rho}{dr^2} +\frac{2}{r} \frac{d\ln\rho}{dr} + \beta^2 \rho = 0 \, , </math>
where,
<math> \beta^2 \equiv \frac{4\pi G}{c_s^2} \, . </math>
This matches the governing ODE whose derivation was published on p. 131 of Robert Emden's (1907) book titled, Gaskugeln.
Derivation by Emden (edited) |
|
By adopting the following dimensionless variables,
<math> \mathfrak{r}_1 \equiv \rho_c^{1/2} \beta r \, , ~~~~\mathrm{and}~~~~v_1 \equiv \ln(\rho/\rho_c) \, , </math>
where <math>\rho_c</math> is the configuration's central density, the governing ODE can be rewritten in dimensionless form as,
<math> \frac{d^2v_1}{d\mathfrak{r}_1^2} +\frac{2}{\mathfrak{r}_1} \frac{dv_1}{dr} + e^{v_1} = 0 \, , </math>
which is exactly the equation numbered (II"a) that can be found on p. 133 of Emden (1907). Emden numerically determined the behavior of the function <math>v_1(\mathfrak{r}_1)</math>, its first derivative with respect to <math>\mathfrak{r}_1</math>, <math>v_1'</math>, along with <math>e^{v_1}</math> and several other useful products, and published his results as Table 14, on p. 135 of his book. The table is reproduced here, primarily for historical purposes.
Emden Numerical Solution
|
|
The set of <math>(x,y)</math> data points plotted in blue are drawn directly from column 1 — actually, <math>x = \ln(\mathfrak{r}_1)</math> — and column 3 — specifically, <math>y = v_1</math> — of Emden's Table 14. The dashed red line has a slope of <math>-2</math> and serves to illustrate that, at large radii, the isothermal density profile tends toward a <math>\rho \propto r^{-2}</math> distribution. |
|
|
Note: The <math>3^\mathrm{rd}</math> column table entry highlighted in blue must be a typesetting error. |
Bonnor and Ebert
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) |
<math>G \Leftrightarrow \gamma</math> | ||
<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> |
Related Wikipedia Discussions
© 2014 - 2021 by Joel E. Tohline |