Difference between revisions of "User:Tohline/SSC/Structure/IsothermalSphere"

From VistrailsWiki
Jump to navigation Jump to search
(→‎Emden Numerical Solution: Rearrange layout of Table 14 and Figure 1)
 
(37 intermediate revisions by the same user not shown)
Line 1: Line 1:
__FORCETOC__
__FORCETOC__
=Isothermal Sphere=
{{LSU_HBook_header}}
{{LSU_HBook_header}}
=Isothermal Sphere (structure)=
Here we supplement the [[User:Tohline/SphericallySymmetricConfigurations/PGE|simplified set of principal governing equations]] with an isothermal equation of state, that is, {{User:Tohline/Math/VAR_Pressure01}} is related to {{User:Tohline/Math/VAR_Density01}} through the relation,  
Here we supplement the [[User:Tohline/SphericallySymmetricConfigurations/PGE|simplified set of principal governing equations]] with an isothermal equation of state, that is, {{User:Tohline/Math/VAR_Pressure01}} is related to {{User:Tohline/Math/VAR_Density01}} through the relation,  
<div align="center">
<div align="center">
<math>P = c_s^2 \rho \, ,</math>
<math>P = c_s^2 \rho \, ,</math>
</div>
</div>
where, <math>c_s</math> is the isothermal sound speed.  Comparing this {{User:Tohline/Math/VAR_Pressure01}}-{{User:Tohline/Math/VAR_Density01}} relationship to
where, <math>~c_s</math> is the isothermal sound speed.  Comparing this {{User:Tohline/Math/VAR_Pressure01}}-{{User:Tohline/Math/VAR_Density01}} relationship to


<div align="center">
<div align="center">
Line 19: Line 19:
<math>c_s^2 = \frac{\Re T}{\bar{\mu}} = \frac{k T}{m_u \bar{\mu}} \, ,</math>
<math>c_s^2 = \frac{\Re T}{\bar{\mu}} = \frac{k T}{m_u \bar{\mu}} \, ,</math>
</div>
</div>
where, {{User:Tohline/Math/C_GasConstant}}, {{User:Tohline/Math/C_BoltzmannConstant}}, {{User:Tohline/Math/C_AtomicMassUnit}}, and {{User:Tohline/Math/MP_MeanMolecularWeight}} are all defined in the accompanying [[User:Tohline/Appendix/Variables_templates|variables appendix]].  It will be useful to note that, for an isothermal gas, {{User:Tohline/Math/VAR_Enthalpy01}} is related to {{User:Tohline/Math/VAR_Density01}} via the expression,
where, {{User:Tohline/Math/C_GasConstant}}, {{User:Tohline/Math/C_BoltzmannConstant}}, {{User:Tohline/Math/C_AtomicMassUnit}}, and {{User:Tohline/Math/MP_MeanMolecularWeight}} are all defined in the accompanying [[User:Tohline/Appendix/Variables_templates|variables appendix]].  It will be useful to note that, for an isothermal gas, the enthalpy, {{User:Tohline/Math/VAR_Enthalpy01}}, is related to {{User:Tohline/Math/VAR_Density01}} via the expression,
 
<div align="center">
<div align="center">
<math>
<math>
Line 25: Line 26:
</math>
</math>
</div>
</div>


==Governing Relations==
==Governing Relations==


Adopting [[User:Tohline/SphericallySymmetricConfigurations/SolutionStrategies#Technique_2|solution technique #2]], we need to solve the following second-order ODE relating the two unknown functions, {{User:Tohline/Math/VAR_Density01}} and {{User:Tohline/Math/VAR_Enthalpy01}}:
Adopting [[User:Tohline/SphericallySymmetricConfigurations/SolutionStrategies#Technique_2|solution technique #2]], we need to solve the following second-order ODE relating the two unknown functions, {{User:Tohline/Math/VAR_Density01}} and {{User:Tohline/Math/VAR_Enthalpy01}}:
<div align="center">
<div align="center">
<math>\frac{1}{r^2} \frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =-  4\pi G \rho</math> .
<math>\frac{1}{r^2} \frac{d}{dr}\biggl( r^2 \frac{dH}{dr} \biggr) =-  4\pi G \rho</math> .
</div>
</div>
Using the {{User:Tohline/Math/VAR_Enthalpy01}}-{{User:Tohline/Math/VAR_Density01}} relationship for an isothermal gas presented above, this can be rewritten entirely in terms of the density as,
Using the {{User:Tohline/Math/VAR_Enthalpy01}}-{{User:Tohline/Math/VAR_Density01}} relationship for an isothermal gas presented above, this can be rewritten entirely in terms of the density as,
<div align="center">
<div align="center">
<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>
<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>
</div>
</div>
or, equivalently,
 
<span id="keyExpression">or, equivalently,</span>
<div align="center">
<div align="center">
<math>
<math>
Line 42: Line 48:
</math>  
</math>  
</div>
</div>
where,
where,
<div align="center">
<div align="center">
<math>
<math>
Line 48: Line 56:
</math>
</math>
</div>
</div>
This matches the governing ODE whose derivation was published on p. 131 of [http://books.google.com/books?id=MiDQAAAAMAAJ&printsec=frontcover#v=onepage&q&f=true Robert Emden's (1907) book titled, ''Gaskugeln''].  
This matches the governing ODE whose derivation was published on p. 131 of [http://books.google.com/books?id=MiDQAAAAMAAJ&printsec=frontcover#v=onepage&q&f=true Robert Emden's (1907) book titled, ''Gaskugeln''].  


Line 63: Line 72:
   <td align="center">
   <td align="center">
[[File:EmdenIsothermalDerivation.jpg|500px|center|Emden (1907)]]
[[File:EmdenIsothermalDerivation.jpg|500px|center|Emden (1907)]]
  </td>
</tr>
<tr>
  <td align="center" colspan="2">
<font size="-1">
Note that, in Emden's derivation, <math>H</math> is not enthalpy but, rather, the effective gas constant, <math>H = c_s^2/T</math>.
</font>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
By adopting the following dimensionless variables,
By adopting the following dimensionless variables,
<div align="center">
<div align="center">
<math>
<math>
Line 73: Line 91:
</math>
</math>
</div>
</div>
where <math>\rho_c</math> is the configuration's central density, the governing ODE can be rewritten in dimensionless form as,
 
where <math>~\rho_c</math> is the configuration's central density, the governing ODE can be rewritten in dimensionless form as,
 
<div align="center">
<div align="center">
<math>
<math>
\frac{d^2v_1}{d\mathfrak{r}_1^2} +\frac{2}{\mathfrak{r}_1} \frac{dv_1}{dr} + e^{v_1} = 0 \, ,
\frac{d^2v_1}{d\mathfrak{r}_1^2} +\frac{2}{\mathfrak{r}_1} \frac{dv_1}{d\mathfrak{r}_1} + e^{v_1} = 0 \, ,
</math>  
</math>  
</div>
</div>
which is exactly the equation numbered (II"a) that can be found on p. 133 of [http://books.google.com/books?id=MiDQAAAAMAAJ&printsec=frontcover#v=onepage&q&f=true Emden (1907)].
which is exactly the equation numbered (II"a) that can be found on p. 133 of [http://books.google.com/books?id=MiDQAAAAMAAJ&printsec=frontcover#v=onepage&q&f=true 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 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.  This table has been reproduced [[#Emden.27s_Numerical_Solution|immediately below]], primarily for historical purposes.
 
Note that a somewhat more extensive tabulation of the structural properties of isothermal spheres is provided by [http://adsabs.harvard.edu/abs/1949ApJ...109..551C Chandrasekhar &amp; Wares (1949, ApJ, 109, 551)].  In this published work as well as in &sect;22 of Chapter IV in [<b>[[User:Tohline/Appendix/References#C67|<font color="red">C67</font>]]</b>], Chandrasekhar has written the governing ODE in a form that we will refer to as the,
<div align="center" id="Chandrasekhar">
<font color="maroon"><b>Isothermal Lane-Emden Equation</b></font><br />
{{ User:Tohline/Math/EQ_SSLaneEmden02 }}
</div>
It is straightforward to show that this is identical to Emden's governing expression after making the variable substitutions:
<div align="center">
<math>~\mathfrak{r}_1 \rightarrow \xi</math>&nbsp; &nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp; &nbsp; <math>~v_1 \rightarrow -\psi </math>.
</div>
Across the astrophysics community, Chadrasekhar's notation has been widely &#8212; although not universally &#8212; adopted as the standard.


==Emden Numerical Solution==
==Emden's Numerical Solution==
<div align="center">
<div align="center">
<table border="1" cellpadding="1">
<table border="1" cellpadding="1">
Line 93: Line 125:
   <td align="center">
   <td align="center">
<font size="-1">
<font size="-1">
Note:  The <math>3^\mathrm{rd}</math> column table entry highlighted in blue must be a typesetting error.
Note:  The entry highlighted in blue in the <math>3^\mathrm{rd}</math> column must be a typesetting error.
</font>
  </td>
</tr>
<tr>
  <td align="left">
<font size="-1">
A more recent and more extensive tabulation of the structural properties of isothermal spheres is provided by:
* [http://adsabs.harvard.edu/abs/1949ApJ...109..551C Chandrasekhar &amp; Wares (1949, ApJ, 109, 551)]: &nbsp; ''The Isothermal Function''
* [http://adsabs.harvard.edu/abs/1986Ap%26SS.126..357H G. P. Horedt (1986, Astrophys. &amp; Space Science, 126, 357-408)]:&nbsp; ''Seven-Digit Tables of Lane-Emden Functions'' &#8212; See, in particular, pp. 405-406 (Sphere of polytropic index <math>~n = \pm \infty</math>).
 
An analytic &#8212; but ''approximate'' &#8212; solution to the isothermal Lane-Emden equation can be found:
* [http://adsabs.harvard.edu/abs/1996MNRAS.281.1197L F. K. Liu (1996, MNRAS, 281, 1197-1205)]: &nbsp; ''Polytropic Gas Spheres:  An Approximate Analytic Solution of the Lane-Emden Equation''
* [http://adsabs.harvard.edu/abs/1997MNRAS.286..268N Priyamvada Natarajan &amp; Donald Lynden-Bell (1997, MNRAS, 286, 268-270)]: &nbsp; '' An Analytic Approximation to the Isothermal Sphere''
* [http://adsabs.harvard.edu/abs/2013RMxAA..49...63R A. C. Raga, J. C. Rodr&iacute;guez-Ram&iacute;rez, M. Villasante, A. Rodr&iacute;guez-Gonz&aacute;lez, &amp; V. Lora (2013, Revista Mexicana de Astronom&iacute;a y Astrof&iacute;sica, 49, 63-69)]: &nbsp; ''A New Analytic Approximation to the Isothermal, Self-Gravitating Sphere''
 
</font>
</font>
   </td>
   </td>
Line 99: Line 146:
</table>
</table>
</div>
</div>
A plot of <math>~v_1</math> versus <math>~\ln\mathfrak{r}_1</math>, as shown below in Figure 1a, translates into a log-log plot of the equilibrium configuration's <math>~\rho(r)</math> density profile.  Notice that this isolated isothermal configuration extends to infinity and that, at large radii, the density profile displays a simple power-law behavior &#8212; specifically, <math>~ \rho \propto r^{-2}</math>.  This is consistent with our general discussion, [[User:Tohline/SSC/Structure/PowerLawDensity#Isothermal_Equation_of_State|presented elsewhere]], of power-law density distributions as solutions of the Lane-Emden equation. 


<div align="center">
<div align="center">
<table border="0" cellpadding="5" width="85%">
<table border="0" cellpadding="5" width="85%">
<tr>
<tr>
   <td align="center" colspan="2">
   <td align="center" colspan="2">
&nbsp;
'''Figure 1:  Emden's Numerical Solution'''
  </td>
</tr>
<tr>
  <td align="center" colspan="2">
'''Figure 1:  Emden's Numerical Solutions'''
   </td>
   </td>
</tr>
</tr>
Line 117: Line 164:
   </td>
   </td>
   <td align="center" valign="top">
   <td align="center" valign="top">
[[File:IsothermalDensityPlot.jpg|350px|center|Plotted from Emden's (1907) tabulated data]]
[[File:EmdenMassProfile.jpg|350px|center|Plotted from Emden's (1907) tabulated data]]
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td valign="top">
   <td valign="top">
The set of <math>(x,y)</math> data points plotted in blue are drawn directly from column 1 &#8212; actually, <math>x = \ln(\mathfrak{r}_1)</math> &#8212; and column 3 &#8212; specifically, <math>y = v_1</math> &#8212; 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.
(a) The <math>~(x,y)</math> locations of the data points plotted in blue are drawn directly from column 1 and column 3 of Emden's Table 14 &#8212; specifically, <math>~x = \ln(\mathfrak{r}_1)</math> and <math>~y = v_1</math>.  The dashed red line has a slope of <math>~-2</math> and serves to illustrate that, at large radii, the [[User:Tohline/SSC/Structure/PowerLawDensity#Isothermal_Equation_of_State|isothermal density profile tends toward a <math>~\rho \propto r^{-2}</math> distribution]].
   </td>
   </td>
   <td valign="top">
   <td valign="top">
The set of <math>(x,y)</math> data points plotted in blue are drawn directly from column 1 &#8212; actually, <math>x = \ln(\mathfrak{r}_1)</math> &#8212; and column 3 &#8212; specifically, <math>y = v_1</math> &#8212; 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.
(b) The <math>~(x,y)</math> locations of the data points plotted in purple are drawn directly from column 1 and column 7 of Emden's Table 14 &#8212; specifically, <math>~x = \ln(\mathfrak{r}_1)</math> and <math>~y = \mathfrak{r}_1^2 v_1'</math>.  The dashed green line has a slope of <math>~+1</math> and serves to illustrate that, at large radii, the isothermal <math>~M(r)</math> distribution tends toward a <math>~M_r \propto r</math> distribution.
  </td>
</tr>
</table>
</div>
 
==Mass Profile==
The mass enclosed within a given radius, <math>~M_r</math>, can be determined by performing an appropriate volume-weighted integral over the density distribution.  Specifically, based on the key expression for,
 
<div align="center">
<span id="HydrostaticBalance"><font color="#770000">'''Mass Conservation'''</font></span><br />
 
{{User:Tohline/Math/EQ_SSmassConservation01}}
</div>
in spherically symmetric configurations, the relevant integral is,
<div align="center">
<math>
~M_r = \int_0^r 4\pi r^2 \rho(r) dr \, .
</math>
</div>
But <math>~M_r</math> also can be determined from the information provided in column 7 of Emden's Table 14 &#8212; that is, from knowledge of the first derivative of <math>~v_1</math>.  The appropriate expression can be obtained from the mathematical prescription for
 
<div align="center">
<span id="HydrostaticBalance"><font color="#770000">'''Hydrostatic Balance'''</font></span><br />
 
{{User:Tohline/Math/EQ_SShydrostaticBalance01}}
</div>
in a spherically symmetric configuration.  Since, for an isothermal equation of state (see above),
 
<div align="center">
<math>
~\frac{dP}{\rho} = c_s^2 {d\ln\rho} \, ,
</math>
</div>
 
the statement of hydrostatic balance can be rewritten as,
 
<div align="center">
<math>
~M_r = \frac{c_s^2}{G} \biggl[ - r^2 \frac{d\ln\rho}{dr} \biggr] = \frac{c_s^2}{G \rho_c^{1/2} \beta} \biggl[ - \mathfrak{r}_1^2 \frac{dv_1}{d\mathfrak{r}_1} \biggr]
= \biggl( \frac{c_s^6}{4\pi G^3 \rho_c} \biggr)^{1/2}  \biggl[ - \mathfrak{r}_1^2 v_1' \biggr] \, .
</math>
</div>
 
The quantity tabulated in column 7 of Emden's Table 14 is precisely the dimensionless term inside the square brackets of this last expression; having units of mass, the coefficient out front sets the mass scale of the equilibrium configuration and depends only on the choice of central density and isothermal sound speed.  Hence, a plot of <math>~\ln(\mathfrak{r}_1^2 v_1')</math> versus <math>~\ln\mathfrak{r}_1</math>, as shown above in Figure 1b, translates into a log-log plot of the equilibrium configuration's <math>~M_r</math> mass profile.  Notice that, along with the radius, the mass of this isolated isothermal configuration extends to infinity and that, at large radii, the mass profile displays a simple power-law behavior &#8212; specifically, <math> ~M_r \propto r^{+1}</math>.
 
As was realized independently by [http://adsabs.harvard.edu/abs/1955ZA.....37..217E Ebert (1955)] and [http://adsabs.harvard.edu/abs/1956MNRAS.116..351B Bonnor (1956)], a spherically symmetric isothermal equilibrium configuration of finite radius and finite mass can be constructed if the system is embedded in a pressure-confining external medium.  We discuss their findings [[User:Tohline/SSC/Structure/BonnorEbert|elsewhere]].
 
<font color="darkblue">
==Summary==
</font>
 
Based on the above derivations, the internal structural properties of an equilibrium isothermal sphere can be described in terms of the tabulated quantities provided in Emden's Table 14 as follows:
* <font color="red">Radial Coordinate Position</font>:
: Given the isothermal sound speed, <math>c_s</math>, and the central density, <math>\rho_c</math>, the radial coordinate is,
<div align="center">
<math>r = ( \rho_c \beta^2 )^{-1/2} \mathfrak{r}_1 = \biggl( \frac{c_s^2}{4\pi G \rho_c} \biggr)^{1/2} \mathfrak{r}_1 </math> .
</div>
 
* <font color="red">Density &amp; Pressure</font>:
: As a function of the radial coordinate, <math>r(\mathfrak{r}_1)</math>, the density profile is,
<div align="center">
<math>\rho(r(\mathfrak{r}_1))= \rho_c e^{v_1(\mathfrak{r}_1)}</math>;
</div>
: and the pressure profile is,
<div align="center">
<math>P(r(\mathfrak{r}_1))= (c_s^2 \rho_c) e^{v_1(\mathfrak{r}_1)}</math>.
</div>
: As has been explicitly pointed out in the above discussion associated with Figure 1a, the density profile &#8212; and, hence, also the pressure profile &#8212; extends to infinity and, at large radii, behaves as a power law; specifically, <math>\rho \propto r^{-2}</math>.
 
* <font color="red">Mass</font>: 
: Given <math>c_s</math> and <math>\rho_c</math>, the natural mass scale is,
<div align="center">
<math>M_0 \equiv \biggl( \frac{c_s^6}{4\pi G^3 \rho_c} \biggr)^{1/2}  </math> ;
</div>
: and, expressed in terms of <math>M_0</math>, the mass that lies interior to radius <math>r</math> is,
<div align="center">
<math>
M_r = M_0 [ - \mathfrak{r}_1^2 v_1' ] \, .
</math>
</div>
: As discussed above in the context of Figure 1b, at large radii, the mass increases linearly with <math>r</math>.  Because the density and pressure profiles extend to infinity, this means that the mass of an isolated isothermal sphere is infinite.
 
* <font color="red">Enthalpy &amp; Gravitational Potential</font>:
: To within an additive constant, the enthalpy distribution is,
<div align="center">
<math>H(r(\mathfrak{r}_1))= c_s^2  [- v_1(\mathfrak{r}_1)]</math>;
</div>
: and the gravitational potential is,
<div align="center">
<math>\Phi(r(\mathfrak{r}_1)) = - H(r(\mathfrak{r}_1))= c_s^2  v_1(\mathfrak{r}_1)</math>.
</div>
 
* <font color="red">Mean-to-Local Density Ratio</font>:
: The ratio of the configuration's mean density, inside a given radius, to its local density at that radius is,
<div align="center">
<math>\frac{\bar{\rho}}{\rho} = \frac{3M_r}{4\pi r^3 \rho} = 3\biggl[- \frac{v_1'}{\mathfrak{r}_1 e^{v_1}} \biggr] </math> .
</div>
: As Figure 2 shows, at large <math>r</math> this density ratio goes to the value of 3, which means that the term inside the square brackets goes to unity at large <math>r</math>.  This behavior is consistent with the limiting power-law behavior of both <math>M_r</math> and <math>\rho</math>, discussed above.
 
<div align="center">
<table border="0" cellpadding="5" width="360">
<tr>
  <td align="center">
'''Figure 2:  From Emden's Tabulated Data'''
  </td>
</tr>
<tr>
  <td align="center">
[[File:PlotMeanToLocalDensity.jpg|350px|center|Plot based on data from Emden's (1907) Table 14]]
  </td>
</tr>
<tr>
  <td align="left">
The blue curve displays an evaluation of the density ratio, <math>[- 3v_1'/(\mathfrak{r}_1 e^{v_1}) ]</math>, as a function of <math>\ln (\mathfrak{r}_1)</math>, as determined from the data presented in Emden's Table 14, shown above.
  </td>
</tr>
</table>
</div>
 
=Our Numerical Integration=
The [[#keyExpression|above governing relation]] &#8212; see especially [[#Chandrasekhar|Chandrasekhar's notation]] &#8212; may be rewritten as (see also, for example, &sect;19.8, eq. 19.35 of [<b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]]</b>]),
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{d^2w}{dr^2} +\frac{2}{r} \frac{d w}{dr}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~e^{-w} \, ,</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
where we appreciate that,
<div align="center">
<math>~w \equiv \ln\biggl(\frac{\rho}{\rho_c}\biggr) \, .</math>
</div>
</div>


==Bonnor and Ebert==
We'll adopt the following finite-difference approximations for the first and second derivatives on a grid of radial spacing, <math>~\Delta_r</math>:
In their studies of pressure-bounded isothermal spheres, [http://adsabs.harvard.edu/abs/1955ZA.....37..217E Ebert] (1955, ZA, 37, 217) and [http://adsabs.harvard.edu/abs/1956MNRAS.116..351B 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.
<div align="center">
<table border="0" cellpadding="5" align="center">


<tr>
  <td align="right">
<math>~w_i'</math>
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~\frac{w_+ - w_-}{2\Delta_r}</math>
  </td>
</tr>
</table>
</div>
and,
<div align="center">
<div align="center">
<table border="1" cellpadding="4">
<table border="0" cellpadding="5" align="center">
 
<tr>
<tr>
   <td align="center">Derivation by [http://adsabs.harvard.edu/abs/1956MNRAS.116..351B Bonnor] (edited)</td>
   <td align="right">
   <td align="center">''translation''</td>
<math>~w_i''</math>
   <td align="center">Derivation by [http://adsabs.harvard.edu/abs/1955ZA.....37..217E Ebert] (edited)</td>
  </td>
   <td align="center">
<math>~\approx</math>
  </td>
   <td align="left">
<math>~\frac{w_+ - 2w_i +w_-}{\Delta_r^2} \, .</math>
  </td>
</tr>
</tr>
</table>
</div>
Our finite-difference approximation of the governing equation is, then,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="center" rowspan="5">
   <td align="right">
[[File:BonnorDerivation01.jpg|300px|center|Bonnor (1956, MNRAS, 116, 351)]]
<math>~r_i \biggl[ \frac{w_+ - 2w_i +w_-}{\Delta_r^2} \biggr] + 2\biggl[  \frac{w_+ - w_-}{2\Delta_r} \biggr]
</math>
  </td>
  <td align="center">
<math>~=</math>
   </td>
   </td>
   <td align="center"><math>G \Leftrightarrow \gamma</math></td>
   <td align="left">
  <td align="center" rowspan="5">
<math>~r_i e^{-w_i} </math>
[[File:EbertDerivation01.jpg|300px|center|Ebert (1955, ZA, 37, 217)]]
   </td>
   </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center"><math>\rho_c \Leftrightarrow \rho_0</math></td>
   <td align="right">
<math>~\Rightarrow ~~~ r_i [ w_+ - 2w_i +w_- ] + \Delta_r [ w_+ - w_- ]
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\Delta_r^2 r_i e^{-w_i} </math>
  </td>
</tr>
</tr>
<tr>
<tr>
   <td align="center"><math>\frac{kT}{m} \Leftarrow c_s^2 \Rightarrow \frac{\Re T_0}{\mu}</math></td>
  <td align="right">
<math>~\Rightarrow ~~~w_+
</math>
  </td>
   <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\Delta_r^2 r_i e^{-w_i} +  2r_i w_i  + w_- (\Delta_r - r_i)}{( \Delta_r + r_i)} \, .</math>
  </td>
</tr>
</tr>
</table>
</div>
Now, for the first two steps away from the center &#8212; where, <math>~w_i = w_0 = 0</math> and <math>~r_i = r_0 = 0</math> &#8212;  we will use the following [[User:Tohline/Appendix/Ramblings/PowerSeriesExpressions#IsothermalLaneEmden|power-series expansion]] (see, for example, eq. 377 from &sect;22 in Chapter IV of [<b>[[User:Tohline/Appendix/References#C67|<font color="red">C67</font>]]</b>]) to determine the value of <math>~w_i</math>:
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="center"><math>\beta^{1/2}\lambda^{-1/2} \Leftrightarrow l_0</math></td>
  <td align="right">
<math>~w_1
</math>
  </td>
   <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\Delta_r^2}{6} - \frac{\Delta_r^4}{120} + \frac{\Delta_r^6}{1890} \, ,</math>
  </td>
</tr>
</tr>
</table>
</div>
and,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
<tr>
   <td align="center"><math>e^{-\psi} \Leftrightarrow \eta</math></td>
  <td align="right">
<math>~w_2
</math>
  </td>
   <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{(2\Delta_r)^2}{6} - \frac{(2\Delta_r)^4}{120} + \frac{(2\Delta_r)^6}{1890} \, .</math>
  </td>
</tr>
</tr>
</table>
</table>
</div>
</div>


=Related Wikipedia Discussions=
=Related Discussions=
* [http://en.wikipedia.org/wiki/Isothermal_Sphere Isothermal Sphere]
==Journal Articles==
* [http://adsabs.harvard.edu/abs/1870AmJS...50...57L J. H. Lane (1870)], '''American Journal of Science''', ''On the Theoretical Temperature of the Sun''
* [https://archive.org/details/mobot31753002152772/page/56 J. H. Lane (1870)], '''The American Journal of Science and Arts''', Vol. 50, pp. 57 - 74: ''On the Theoretical Temperature of the Sun under the Hypothesis of a Gaseous Mass Maintaining Its Volume by Its Internal Heat and Depending on the Laws of Gases Known to Terrestrial Experiment''
 
==Wikipedia==
* [https://en.wikipedia.org/wiki/Robert_Emden Robert Emden]
* [https://en.wikipedia.org/wiki/Emden–Chandrasekhar_equation Emden-Chandrasekhar equation]
* [https://en.wikipedia.org/wiki/Jonathan_Homer_Lane Jonathan Home Lane]
* [https://en.wikipedia.org/wiki/Lane–Emden_equation Lane-Emden equation]




{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 23:02, 8 March 2019


Isothermal Sphere

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

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,

LSU Key.png

<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, the enthalpy, <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)

Emden (1907)
Emden (1907)

Note that, in Emden's derivation, <math>H</math> is not enthalpy but, rather, the effective gas constant, <math>H = c_s^2/T</math>.

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}{d\mathfrak{r}_1} + 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. This table has been reproduced immediately below, primarily for historical purposes.

Note that a somewhat more extensive tabulation of the structural properties of isothermal spheres is provided by Chandrasekhar & Wares (1949, ApJ, 109, 551). In this published work as well as in §22 of Chapter IV in [C67], Chandrasekhar has written the governing ODE in a form that we will refer to as the,

Isothermal Lane-Emden Equation

LSU Key.png

<math>~\frac{1}{\xi^2} \frac{d}{d\xi}\biggl( \xi^2 \frac{d\psi}{d\xi} \biggr) = e^{-\psi}</math>

It is straightforward to show that this is identical to Emden's governing expression after making the variable substitutions:

<math>~\mathfrak{r}_1 \rightarrow \xi</math>        and         <math>~v_1 \rightarrow -\psi </math>.

Across the astrophysics community, Chadrasekhar's notation has been widely — although not universally — adopted as the standard.

Emden's Numerical Solution

Emden's (1907) Table 14

Note: The entry highlighted in blue in the <math>3^\mathrm{rd}</math> column must be a typesetting error.

A more recent and more extensive tabulation of the structural properties of isothermal spheres is provided by:

An analytic — but approximate — solution to the isothermal Lane-Emden equation can be found:


A plot of <math>~v_1</math> versus <math>~\ln\mathfrak{r}_1</math>, as shown below in Figure 1a, translates into a log-log plot of the equilibrium configuration's <math>~\rho(r)</math> density profile. Notice that this isolated isothermal configuration extends to infinity and that, at large radii, the density profile displays a simple power-law behavior — specifically, <math>~ \rho \propto r^{-2}</math>. This is consistent with our general discussion, presented elsewhere, of power-law density distributions as solutions of the Lane-Emden equation.


Figure 1: Emden's Numerical Solution

Plotted from Emden's (1907) tabulated data
Plotted from Emden's (1907) tabulated data

(a) The <math>~(x,y)</math> locations of the data points plotted in blue are drawn directly from column 1 and column 3 of Emden's Table 14 — specifically, <math>~x = \ln(\mathfrak{r}_1)</math> and <math>~y = v_1</math>. 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.

(b) The <math>~(x,y)</math> locations of the data points plotted in purple are drawn directly from column 1 and column 7 of Emden's Table 14 — specifically, <math>~x = \ln(\mathfrak{r}_1)</math> and <math>~y = \mathfrak{r}_1^2 v_1'</math>. The dashed green line has a slope of <math>~+1</math> and serves to illustrate that, at large radii, the isothermal <math>~M(r)</math> distribution tends toward a <math>~M_r \propto r</math> distribution.

Mass Profile

The mass enclosed within a given radius, <math>~M_r</math>, can be determined by performing an appropriate volume-weighted integral over the density distribution. Specifically, based on the key expression for,

Mass Conservation

LSU Key.png

<math>~\frac{dM_r}{dr} = 4\pi r^2 \rho</math>

in spherically symmetric configurations, the relevant integral is,

<math> ~M_r = \int_0^r 4\pi r^2 \rho(r) dr \, . </math>

But <math>~M_r</math> also can be determined from the information provided in column 7 of Emden's Table 14 — that is, from knowledge of the first derivative of <math>~v_1</math>. The appropriate expression can be obtained from the mathematical prescription for

Hydrostatic Balance

LSU Key.png

<math>~\frac{dP}{dr} = - \frac{GM_r \rho}{r^2}</math>

in a spherically symmetric configuration. Since, for an isothermal equation of state (see above),

<math> ~\frac{dP}{\rho} = c_s^2 {d\ln\rho} \, , </math>

the statement of hydrostatic balance can be rewritten as,

<math> ~M_r = \frac{c_s^2}{G} \biggl[ - r^2 \frac{d\ln\rho}{dr} \biggr] = \frac{c_s^2}{G \rho_c^{1/2} \beta} \biggl[ - \mathfrak{r}_1^2 \frac{dv_1}{d\mathfrak{r}_1} \biggr] = \biggl( \frac{c_s^6}{4\pi G^3 \rho_c} \biggr)^{1/2} \biggl[ - \mathfrak{r}_1^2 v_1' \biggr] \, . </math>

The quantity tabulated in column 7 of Emden's Table 14 is precisely the dimensionless term inside the square brackets of this last expression; having units of mass, the coefficient out front sets the mass scale of the equilibrium configuration and depends only on the choice of central density and isothermal sound speed. Hence, a plot of <math>~\ln(\mathfrak{r}_1^2 v_1')</math> versus <math>~\ln\mathfrak{r}_1</math>, as shown above in Figure 1b, translates into a log-log plot of the equilibrium configuration's <math>~M_r</math> mass profile. Notice that, along with the radius, the mass of this isolated isothermal configuration extends to infinity and that, at large radii, the mass profile displays a simple power-law behavior — specifically, <math> ~M_r \propto r^{+1}</math>.

As was realized independently by Ebert (1955) and Bonnor (1956), a spherically symmetric isothermal equilibrium configuration of finite radius and finite mass can be constructed if the system is embedded in a pressure-confining external medium. We discuss their findings elsewhere.

Summary

Based on the above derivations, the internal structural properties of an equilibrium isothermal sphere can be described in terms of the tabulated quantities provided in Emden's Table 14 as follows:

  • Radial Coordinate Position:
Given the isothermal sound speed, <math>c_s</math>, and the central density, <math>\rho_c</math>, the radial coordinate is,

<math>r = ( \rho_c \beta^2 )^{-1/2} \mathfrak{r}_1 = \biggl( \frac{c_s^2}{4\pi G \rho_c} \biggr)^{1/2} \mathfrak{r}_1 </math> .

  • Density & Pressure:
As a function of the radial coordinate, <math>r(\mathfrak{r}_1)</math>, the density profile is,

<math>\rho(r(\mathfrak{r}_1))= \rho_c e^{v_1(\mathfrak{r}_1)}</math>;

and the pressure profile is,

<math>P(r(\mathfrak{r}_1))= (c_s^2 \rho_c) e^{v_1(\mathfrak{r}_1)}</math>.

As has been explicitly pointed out in the above discussion associated with Figure 1a, the density profile — and, hence, also the pressure profile — extends to infinity and, at large radii, behaves as a power law; specifically, <math>\rho \propto r^{-2}</math>.
  • Mass:
Given <math>c_s</math> and <math>\rho_c</math>, the natural mass scale is,

<math>M_0 \equiv \biggl( \frac{c_s^6}{4\pi G^3 \rho_c} \biggr)^{1/2} </math> ;

and, expressed in terms of <math>M_0</math>, the mass that lies interior to radius <math>r</math> is,

<math> M_r = M_0 [ - \mathfrak{r}_1^2 v_1' ] \, . </math>

As discussed above in the context of Figure 1b, at large radii, the mass increases linearly with <math>r</math>. Because the density and pressure profiles extend to infinity, this means that the mass of an isolated isothermal sphere is infinite.
  • Enthalpy & Gravitational Potential:
To within an additive constant, the enthalpy distribution is,

<math>H(r(\mathfrak{r}_1))= c_s^2 [- v_1(\mathfrak{r}_1)]</math>;

and the gravitational potential is,

<math>\Phi(r(\mathfrak{r}_1)) = - H(r(\mathfrak{r}_1))= c_s^2 v_1(\mathfrak{r}_1)</math>.

  • Mean-to-Local Density Ratio:
The ratio of the configuration's mean density, inside a given radius, to its local density at that radius is,

<math>\frac{\bar{\rho}}{\rho} = \frac{3M_r}{4\pi r^3 \rho} = 3\biggl[- \frac{v_1'}{\mathfrak{r}_1 e^{v_1}} \biggr] </math> .

As Figure 2 shows, at large <math>r</math> this density ratio goes to the value of 3, which means that the term inside the square brackets goes to unity at large <math>r</math>. This behavior is consistent with the limiting power-law behavior of both <math>M_r</math> and <math>\rho</math>, discussed above.

Figure 2: From Emden's Tabulated Data

Plot based on data from Emden's (1907) Table 14

The blue curve displays an evaluation of the density ratio, <math>[- 3v_1'/(\mathfrak{r}_1 e^{v_1}) ]</math>, as a function of <math>\ln (\mathfrak{r}_1)</math>, as determined from the data presented in Emden's Table 14, shown above.

Our Numerical Integration

The above governing relation — see especially Chandrasekhar's notation — may be rewritten as (see also, for example, §19.8, eq. 19.35 of [KW94]),

<math>~\frac{d^2w}{dr^2} +\frac{2}{r} \frac{d w}{dr} </math>

<math>~=</math>

<math>~e^{-w} \, ,</math>

where we appreciate that,

<math>~w \equiv \ln\biggl(\frac{\rho}{\rho_c}\biggr) \, .</math>

We'll adopt the following finite-difference approximations for the first and second derivatives on a grid of radial spacing, <math>~\Delta_r</math>:

<math>~w_i'</math>

<math>~\approx</math>

<math>~\frac{w_+ - w_-}{2\Delta_r}</math>

and,

<math>~w_i</math>

<math>~\approx</math>

<math>~\frac{w_+ - 2w_i +w_-}{\Delta_r^2} \, .</math>

Our finite-difference approximation of the governing equation is, then,

<math>~r_i \biggl[ \frac{w_+ - 2w_i +w_-}{\Delta_r^2} \biggr] + 2\biggl[ \frac{w_+ - w_-}{2\Delta_r} \biggr] </math>

<math>~=</math>

<math>~r_i e^{-w_i} </math>

<math>~\Rightarrow ~~~ r_i [ w_+ - 2w_i +w_- ] + \Delta_r [ w_+ - w_- ] </math>

<math>~=</math>

<math>~\Delta_r^2 r_i e^{-w_i} </math>

<math>~\Rightarrow ~~~w_+ </math>

<math>~=</math>

<math>~\frac{\Delta_r^2 r_i e^{-w_i} + 2r_i w_i + w_- (\Delta_r - r_i)}{( \Delta_r + r_i)} \, .</math>

Now, for the first two steps away from the center — where, <math>~w_i = w_0 = 0</math> and <math>~r_i = r_0 = 0</math> — we will use the following power-series expansion (see, for example, eq. 377 from §22 in Chapter IV of [C67]) to determine the value of <math>~w_i</math>:

<math>~w_1 </math>

<math>~=</math>

<math>~\frac{\Delta_r^2}{6} - \frac{\Delta_r^4}{120} + \frac{\Delta_r^6}{1890} \, ,</math>

and,

<math>~w_2 </math>

<math>~=</math>

<math>~\frac{(2\Delta_r)^2}{6} - \frac{(2\Delta_r)^4}{120} + \frac{(2\Delta_r)^6}{1890} \, .</math>

Related Discussions

Journal Articles

  • J. H. Lane (1870), American Journal of Science, On the Theoretical Temperature of the Sun
  • J. H. Lane (1870), The American Journal of Science and Arts, Vol. 50, pp. 57 - 74: On the Theoretical Temperature of the Sun under the Hypothesis of a Gaseous Mass Maintaining Its Volume by Its Internal Heat and Depending on the Laws of Gases Known to Terrestrial Experiment

Wikipedia


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