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

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


<div align="center">
<div align="center">
<math>~M_u = M_{u,3} \biggl(\frac{K}{c^2}\biggr)^{(n_p-3)/2} \, ,</math>
<math>~M_u = M_\mathrm{norm} \biggl(\frac{K}{c^2}\biggr)^{(n_p-3)/2} \, ,</math>
</div>
</div>


and to determine, first, an expression for the mass-normalization when <math>~n_p = 3</math>, namely,
and to determine, first, an expression for the mass-normalization when <math>~n_p = 3</math>, namely,
<div align="center">
<div align="center">
<math>~M_{u,3} \equiv \biggl( \frac{K}{G}\biggr)^{3/2} .</math>
<math>~M_\mathrm{norm} \equiv \biggl( \frac{K}{G}\biggr)^{3/2} .</math>
</div>
</div>


Line 159: Line 159:
</table>
</table>
</div>
</div>
In terms of the BAC84 terminology, this means that,
Adopting the BAC84 terminology, this means that,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 165: Line 165:
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~M_{u,3}^{2} = \biggl(\frac{K}{G}\biggr)^3</math>
<math>~M_\mathrm{norm}^{2} = \biggl(\frac{K}{G}\biggr)^3</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
Line 183: Line 183:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl( \frac{k Y_T}{m_B}\biggr)^4 \biggl[\sigma^4(1+\sigma^{-1})^3\biggr] \frac{3}{G^3 a_\mathrm{rad}} \, .</math>
<math>~\biggl( \frac{k Y_T}{m_B}\biggr)^4 \biggl[\sigma^4(1+\sigma^{-1})^3\biggr] \frac{3}{G^3 a_\mathrm{rad}} </math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~~ M_\mathrm{norm} = \biggl(\frac{K}{G}\biggr)^{3/2}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~(1+\sigma^{-1})^{3/2} \biggl( \frac{k Y_T}{m_B}\biggr)^2  \biggl(\frac{3}{G^3 a_\mathrm{rad}}\biggr)^{1/2} \sigma^2 </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left">
<math>~\biggl(1+\frac{3}{2\sigma} \biggr) \biggl( \frac{k Y_T}{m_B}\biggr)^2  \biggl(\frac{3}{G^3 a_\mathrm{rad}}\biggr)^{1/2} \sigma^2 \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
When radiation pressure significantly dominates over gas pressure &#8212; that is, in the limit <math>~\sigma >> 1</math> &#8212; the factor of <math>~(1+\sigma^{-1})^3 \approx 1</math>, and we see that this expression for <math>~M_{u,3}^2</math> exactly matches equation (10) of BAC84.
When radiation pressure significantly dominates over gas pressure &#8212; that is, in the limit <math>~\sigma \gg 1</math> &#8212; the factor of <math>~(1+\sigma^{-1})^3 \approx 1</math>, and we see that this expression for <math>~M_\mathrm{norm}</math> exactly matches the expression for <math>~M_{u,3}</math> given by equation (10) of BAC84.


====Polytropic Index Slightly Less Than 3====
====Polytropic Index Slightly Less Than 3====
More generally, equating the two expressions for the total pressure and drawing (twice) on the expression for <math>~\sigma</math> provided above, we have,
More generally, equating the two expressions for the total pressure and drawing (twice) on the expression for <math>~\sigma</math> [[#Ratio_of_Radiation_Pressure_to_Gas_Pressure|provided above]], we have,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 280: Line 304:
</table>
</table>
</div>
</div>
This means that, for any <math>~\sigma >> 1</math>,
This means that, for any <math>~\sigma \gg 1</math>,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 292: Line 316:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1})\biggl[ \frac{Y_T  k \sigma}{m_B} \biggr]^{4/3} \rho^{-1/6\sigma} \, .</math>
<math>~\biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1})\biggl( \frac{Y_T  k \sigma}{m_B} \biggr)^{4/3} \rho^{-1/6\sigma} \, .</math>
   </td>
   </td>
</tr>
</tr>
</table>
</table>
</div>
</div>
This matches exactly expression (7) in BAC84.  Again from above &#8212; and continuing to assume <math>~\sigma >> 1</math> &#8212; we have,
This matches exactly expression (7) in BAC84.  Again from above &#8212; and continuing to assume <math>~\sigma \gg 1</math> &#8212; we have,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 350: Line 374:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~M_{u,3} \biggl(\frac{K}{c^2}\biggr)^{-3/4\sigma}</math>
<math>~M_\mathrm{norm} \biggl(\frac{K}{c^2}\biggr)^{-3/4\sigma}</math>
   </td>
   </td>
</tr>
</tr>
Line 362: Line 386:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~M_{u,3} \biggl[\frac{1}{c^2} \biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1})
<math>~M_\mathrm{norm} \biggl[\frac{1}{c^2} \biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1})
\biggl( \frac{Y_T  k \sigma}{m_B} \biggr)^{4/3} \rho^{-1/6\sigma} \biggr]^{-3/4\sigma} \, .</math>
\biggl( \frac{Y_T  k \sigma}{m_B} \biggr)^{4/3} \rho^{-1/6\sigma} \biggr]^{-3/4\sigma} \, .</math>
   </td>
   </td>
Line 368: Line 392:
</table>
</table>
</div>
</div>
Now, because,
Drawing again from the definition of <math>~\sigma</math> [[#Ratio_of_Radiation_Pressure_to_Gas_Pressure|provided above]], we have,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 385: Line 409:
</table>
</table>
</div>
</div>
this last relation can be rewritten as,
so this last relation can be rewritten as,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 397: Line 421:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~M_{u,3} \biggl[(1+\sigma^{-1})\rho^{-1/6\sigma} \biggr]^{-3/4\sigma}  
<math>~M_\mathrm{norm}\biggl[(1+\sigma^{-1})\rho^{-1/6\sigma} \biggr]^{-3/4\sigma}  
\biggl[\frac{T \sigma Y_T k}{m_B c^2} (n m_B)^{-1/3}\biggr]^{-3/4\sigma} </math>
\biggl[\frac{T \sigma Y_T k}{m_B c^2} (n m_B)^{-1/3}\biggr]^{-3/4\sigma} </math>
   </td>
   </td>
Line 410: Line 434:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~M_{u,3} \biggl[(1+\sigma^{-1})\rho^{-1/6\sigma} \biggr]^{-3/4\sigma}  
<math>~M_\mathrm{norm} \biggl[(1+\sigma^{-1})\rho^{-1/6\sigma} \biggr]^{-3/4\sigma}  
\rho^{1/4\sigma}  
\rho^{1/4\sigma}  
\biggl[\frac{m_B c^2}{T \sigma Y_T k} \biggr]^{3/4\sigma}  
\biggl[\frac{m_B c^2}{T \sigma Y_T k} \biggr]^{3/4\sigma}  

Revision as of 17:15, 18 December 2015

Rotating, Supermassive Stars

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


Here we draw upon the work of J. R. Bond, W. D. Arnett, & B. J. Carr (1984, ApJ, 380, 825; hereafter BAC84) who were among the first to seriously address the question of the fate of very massive (stellar) objects.

Equation of State

Our discussion of the equation of state (EOS) that was used by BAC84 draws on the terminology that has already been adopted in our introductory discussion of supplemental relations and closely parallels our review of the properties of the envelope that E. A. Milne (1930, MNRAS, 91, 4) used to construct a bipolytropic sphere.

Expression for Total Pressure

Ignoring the component due to a degenerate electron gas, <math>~P_\mathrm{deg}</math>, the total gas pressure can be expressed as the sum of two separate components: a component of ideal gas pressure, and a component of radiation pressure. That is, in BAC84 the total pressure is given by the expression,

<math>~P</math>

<math>~=</math>

<math>~P_\mathrm{gas} + P_\mathrm{rad} \, ,</math>

where,

Ideal Gas Radiation

LSU Key.png

<math>~P_\mathrm{gas} = \frac{\Re}{\bar{\mu}} \rho T</math>

LSU Key.png

<math>~P_\mathrm{rad} = \frac{1}{3} a_\mathrm{rad} T^4</math>

Now, BAC84 define the rest-mass density in terms of the mean baryon mass, <math>~m_B</math>, via the expression, <math>~\rho = m_B n</math>, and write (see their equation 1),

<math>~P</math>

<math>~=</math>

<math>~Y_T n k T + \frac{1}{3}a_\mathrm{rad} T^4 \, .</math>

In converting from our notation to theirs we conclude, therefore, that,

<math>~\frac{\Re}{\bar{\mu}} (m_B n) T</math>

<math>~=</math>

<math>~Y_T n k T </math>

<math>~\Rightarrow ~~~~ Y_T </math>

<math>~=</math>

<math>~\frac{\Re}{k} \cdot \frac{m_B}{\bar{\mu}} \, .</math>

Ratio of Radiation Pressure to Gas Pressure

Following Milne (1930), we have defined the parameter, <math>~\beta</math>, as the ratio of gas pressure to total pressure. That is, in the context of BAC84, we have,

<math>\beta \equiv \frac{P_\mathrm{gas} }{P} \, ,</math>

in which case, also,

<math>\frac{P_\mathrm{rad}}{P} = 1-\beta </math>         and         <math>\frac{P_\mathrm{gas}}{P_\mathrm{rad}} = \frac{\beta}{1-\beta} \, .</math>

Using a different notation, BAC84 (see their equation 5) define <math>~\sigma</math> as the ratio of the radiation pressure to the gas pressure. Therefore, in converting from our notation to theirs we have,

<math>\sigma = \frac{1-\beta}{\beta} ~~~~\Rightarrow ~~~~ \beta = (1 + \sigma)^{-1} \, , </math>

as well as,

<math>\sigma = \frac{P_\mathrm{rad}}{P_\mathrm{gas}} = \frac{a_\mathrm{rad} T^3}{3} \cdot \frac{\bar\mu}{(\Re m_B)n} = \frac{a_\mathrm{rad} T^3}{3Y_T n k} \, ,</math>

which is precisely the definition provided in equation (5) of BAC84.

Mass Normalization

Now, according to BAC84 (see their equation 8), when the total pressure is written in polytropic form — specifically, if we set,

<math>P = K\rho^{(1+1/n_p)} </math>

— the mass-scaling for relativistic configurations will depend on <math>~G</math>, <math>~c</math>, <math>~K</math>, and <math>~n_p</math> via the expression,

<math>~M_u = K^{n_p/2} G^{-3/2} c^{3-n_p} = \biggl( \frac{K}{G}\biggr)^{3/2} \biggl(\frac{K}{c^2}\biggr)^{(n_p-3)/2} \, .</math>

It is convenient to rewrite this expression in the form,

<math>~M_u = M_\mathrm{norm} \biggl(\frac{K}{c^2}\biggr)^{(n_p-3)/2} \, ,</math>

and to determine, first, an expression for the mass-normalization when <math>~n_p = 3</math>, namely,

<math>~M_\mathrm{norm} \equiv \biggl( \frac{K}{G}\biggr)^{3/2} .</math>

Polytropic Index Equals 3

Referencing our separate discussion of Milne's (1930) work, when <math>~n_p = 3</math>, the polytropic constant is related to the relevant set of physical parameters via the relation,

<math>~K</math>

<math>~=</math>

<math>~\biggl[ \biggl( \frac{\Re}{\bar\mu}\biggr)^4 \biggl(\frac{1-\beta}{\beta^4}\biggr) \frac{3}{a_\mathrm{rad}} \biggr]^{1/3} \, .</math>

Adopting the BAC84 terminology, this means that,

<math>~M_\mathrm{norm}^{2} = \biggl(\frac{K}{G}\biggr)^3</math>

<math>~=</math>

<math>~\biggl( \frac{\Re}{\bar\mu}\biggr)^4 \biggl[\frac{1-\beta}{\beta^4}\biggr] \frac{3}{G^3 a_\mathrm{rad}} </math>

 

<math>~=</math>

<math>~\biggl( \frac{k Y_T}{m_B}\biggr)^4 \biggl[\sigma^4(1+\sigma^{-1})^3\biggr] \frac{3}{G^3 a_\mathrm{rad}} </math>

<math>~\Rightarrow ~~~~ M_\mathrm{norm} = \biggl(\frac{K}{G}\biggr)^{3/2}</math>

<math>~=</math>

<math>~(1+\sigma^{-1})^{3/2} \biggl( \frac{k Y_T}{m_B}\biggr)^2 \biggl(\frac{3}{G^3 a_\mathrm{rad}}\biggr)^{1/2} \sigma^2 </math>

 

<math>~\approx</math>

<math>~\biggl(1+\frac{3}{2\sigma} \biggr) \biggl( \frac{k Y_T}{m_B}\biggr)^2 \biggl(\frac{3}{G^3 a_\mathrm{rad}}\biggr)^{1/2} \sigma^2 \, .</math>

When radiation pressure significantly dominates over gas pressure — that is, in the limit <math>~\sigma \gg 1</math> — the factor of <math>~(1+\sigma^{-1})^3 \approx 1</math>, and we see that this expression for <math>~M_\mathrm{norm}</math> exactly matches the expression for <math>~M_{u,3}</math> given by equation (10) of BAC84.

Polytropic Index Slightly Less Than 3

More generally, equating the two expressions for the total pressure and drawing (twice) on the expression for <math>~\sigma</math> provided above, we have,

<math>~K\rho^{(1 + 1/n_p)}</math>

<math>~=</math>

<math>~Y_T n k T + \frac{a_\mathrm{rad}}{3} T^4 </math>

 

<math>~=</math>

<math>~\frac{a_\mathrm{rad}}{3} (1+\sigma^{-1})T^4 </math>

 

<math>~=</math>

<math>~\frac{a_\mathrm{rad}}{3} (1+\sigma^{-1})\biggl[ \frac{3Y_T n k \sigma}{a_\mathrm{rad}} \biggr]^{4/3} </math>

 

<math>~=</math>

<math>~\biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1})\biggl[ Y_T n k \sigma \biggr]^{4/3} \, .</math>

Now, from above we have,

<math>~1 + \frac{1}{n_p} = \Gamma</math>

<math>~\approx</math>

<math>~\frac{4}{3} + \frac{1}{6\sigma} \, ,</math>

so the lefthand-side of this last expression can be written as,

<math>~K\rho^{(1+1/n_p)}</math>

<math>~\approx</math>

<math>~K\rho^{(4/3+1/6\sigma)} = K(m_B n)^{4/3} \rho^{1/6\sigma} \, .</math>

This means that, for any <math>~\sigma \gg 1</math>,

<math>~K </math>

<math>~\approx</math>

<math>~\biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1})\biggl( \frac{Y_T k \sigma}{m_B} \biggr)^{4/3} \rho^{-1/6\sigma} \, .</math>

This matches exactly expression (7) in BAC84. Again from above — and continuing to assume <math>~\sigma \gg 1</math> — we have,

<math>~1 + \frac{1}{n_p} \approx \frac{4}{3} + \frac{1}{6\sigma} </math>

      <math>~~~~\Rightarrow ~~~~</math>     

<math>~\frac{1}{n_p} \approx \frac{1}{3}\biggl(1 + \frac{1}{2\sigma}\biggr) </math>

 

      <math>~~~~\Rightarrow ~~~~</math>     

<math>~n_p \approx 3\biggl(1 + \frac{1}{2\sigma}\biggr)^{-1} \approx 3\biggl(1 - \frac{1}{2\sigma}\biggr)</math>

 

      <math>~~~~\Rightarrow ~~~~</math>     

<math>~\frac{(n_p-3)}{2} \approx - \frac{3}{4\sigma} \, .</math>

Hence, when the polytropic index is slightly less than 3, the mass normalization is,

<math>~M_u</math>

<math>~\approx</math>

<math>~M_\mathrm{norm} \biggl(\frac{K}{c^2}\biggr)^{-3/4\sigma}</math>

 

<math>~\approx</math>

<math>~M_\mathrm{norm} \biggl[\frac{1}{c^2} \biggl( \frac{3}{a_\mathrm{rad}} \biggr)^{1/3}(1+\sigma^{-1}) \biggl( \frac{Y_T k \sigma}{m_B} \biggr)^{4/3} \rho^{-1/6\sigma} \biggr]^{-3/4\sigma} \, .</math>

Drawing again from the definition of <math>~\sigma</math> provided above, we have,

<math>~\biggl(\frac{3}{a_\mathrm{rad}}\biggr)^{1/3}</math>

<math>~=</math>

<math>~T (\sigma Y_T n k)^{-1/3} \, ,</math>

so this last relation can be rewritten as,

<math>~M_u</math>

<math>~\approx</math>

<math>~M_\mathrm{norm}\biggl[(1+\sigma^{-1})\rho^{-1/6\sigma} \biggr]^{-3/4\sigma} \biggl[\frac{T \sigma Y_T k}{m_B c^2} (n m_B)^{-1/3}\biggr]^{-3/4\sigma} </math>

 

<math>~=</math>

<math>~M_\mathrm{norm} \biggl[(1+\sigma^{-1})\rho^{-1/6\sigma} \biggr]^{-3/4\sigma} \rho^{1/4\sigma} \biggl[\frac{m_B c^2}{T \sigma Y_T k} \biggr]^{3/4\sigma} </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