Difference between revisions of "User:Tohline/SR"

From VistrailsWiki
Jump to navigation Jump to search
(→‎Nonrelativistic ZTF Gas: Begin working on series expansions)
 
(28 intermediate revisions by the same user not shown)
Line 1: Line 1:
__FORCETOC__
__FORCETOC__
=Supplemental Relations=
=Supplemental Relations=
<table border="0" cellpadding="3" align="center" width="60%">
<tr><td align="left">
<font color="darkgreen">
"The equilibrium configuration of a stellar system depends crucially on the equation of state.  Similarly, the time evolution of the perturbation will depend on the relationship that is assumed between the perturbed values of pressure, temperature, and density.  The simplest assumption regarding these perturbations will be to involve an adiabatic relationship between perturbed pressure and density."</font>
</td></tr>
<tr><td align="right">
&#8212; Drawn from &sect;7, p. 172 of  [<b>[[User:Tohline/Appendix/References#P00|<font color="red">P00</font>]]</b>]
</td></tr></table>
{{LSU_HBook_header}}
{{LSU_HBook_header}}
Examining the set of [http://www.vistrails.org/index.php/User:Tohline/PGE principal governing equations] that describe ''adiabatic'' flows, we see that &#8212; apart from the independent variables {{User:Tohline/Math/VAR_Time01}} and {{User:Tohline/Math/VAR_PositionVector01}} &#8212; the equations involve the vector velocity {{User:Tohline/Math/VAR_VelocityVector01}},  and the four scalar variables, {{User:Tohline/Math/VAR_NewtonianPotential01}}, {{User:Tohline/Math/VAR_Pressure01}}, {{User:Tohline/Math/VAR_Density01}}, and {{User:Tohline/Math/VAR_SpecificInternalEnergy01}}. Because the variables outnumber the equations by one, a supplemental relationship between the physical variables must be specified in order to close the set of equations.  (If nonadiabatic flows are considered, additional supplemental relations must be specified because the scalar variables {{User:Tohline/Math/VAR_Temperature01}} and {{User:Tohline/Math/VAR_SpecificEntropy01}} enter the discussion as well.)  Also, in order to complete the unique specification of a particular physical problem, either a steady-state flow field or initial conditions must be specified, depending on whether one is studying a time-independent (structure) or time-dependent (stability or dynamics) problem, respectively.  Throughout this H_Book, the following strategy will be adopted in order to complete the physical specification of each examined system:  (A) For time-independent problems, we will adopt a structural relationship between {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}}, and specify a steady-state flow-field.  (B) For time-dependent problems, we will adopt an equation of state, and specify initial conditions.
Examining the set of [[User:Tohline/PGE#Principal_Governing_Equations|principal governing equations]] that describe ''adiabatic'' flows, we see that &#8212; apart from the independent variables {{User:Tohline/Math/VAR_Time01}} and {{User:Tohline/Math/VAR_PositionVector01}} &#8212; the equations involve the vector velocity {{User:Tohline/Math/VAR_VelocityVector01}},  and the four scalar variables, {{User:Tohline/Math/VAR_NewtonianPotential01}}, {{User:Tohline/Math/VAR_Pressure01}}, {{User:Tohline/Math/VAR_Density01}}, and {{User:Tohline/Math/VAR_SpecificInternalEnergy01}}. Because the variables outnumber the equations by one, a supplemental relationship between the physical variables must be specified in order to close the set of equations.  (If nonadiabatic flows are considered, additional supplemental relations must be specified because the scalar variables {{User:Tohline/Math/VAR_Temperature01}} and {{User:Tohline/Math/VAR_SpecificEntropy01}} enter the discussion as well.)  Also, in order to complete the unique specification of a particular physical problem, either a steady-state flow field or initial conditions must be specified, depending on whether one is studying a time-independent (structure) or time-dependent (stability or dynamics) problem, respectively.  Throughout this H_Book, the following strategy will be adopted in order to complete the physical specification of each examined system:  (A) For time-independent problems, we will adopt a structural relationship between {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}}, and specify a steady-state flow-field.  (B) For time-dependent problems, we will adopt an equation of state, and specify initial conditions.




Line 8: Line 18:
===Barotropic Structure===
===Barotropic Structure===


For ''time-independent'' problems, a structural relationship between {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}} is required to close the system of [http://www.vistrails.org/index.php/User:Tohline/PGE principal governing equations].  [Tassoul (1978) refers to this as a "geometrical" rather than a "structural" relationship; see the discussion associated with his Chapter 4, Eq. 14.] Generally throughout this H_Book, we will assume that all time-independent configurations can be described as [http://en.wikipedia.org/wiki/Barotropic barotropic] structures; that is, we will assume that {{User:Tohline/Math/VAR_Pressure01}} is only a function of {{User:Tohline/Math/VAR_Density01}} throughout such structures.  More specifically, we generally will adopt one of the two ''analytically'' prescribable {{User:Tohline/Math/VAR_Pressure01}}({{User:Tohline/Math/VAR_Density01}}) relationships displayed in the first row of the following Table.
For ''time-independent'' problems, a structural relationship between {{User:Tohline/Math/VAR_Pressure01}} and {{User:Tohline/Math/VAR_Density01}} is required to close the system of [[User:Tohline/PGE#Principal_Governing_Equations|principal governing equations]].  [<b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]]</b>] refers to this as a "geometrical" rather than a "structural" relationship; see the discussion associated with his Chapter 4, Eq. 14. Generally throughout this H_Book, we will assume that all time-independent configurations can be described as barotropic structures; that is, we will assume that {{User:Tohline/Math/VAR_Pressure01}} is only a function of {{User:Tohline/Math/VAR_Density01}} throughout such structures. (The [[User:Tohline/2DStructure/AxisymmetricInstabilities#Poincar.C3.A9-Wavre_Theorem|Poincar&eacute;-Wavre theorem]] provides additional support for, as well as additional implications of, this assumption.) More specifically, we generally will adopt one of the two ''analytically'' prescribable <math>~P(\rho)</math> relationships displayed in the first row of the following Table.


<table width="80%" align="center" border=1 cellpadding=5>
<table width="80%" align="center" border=1 cellpadding=5>
Line 26: Line 36:
<td align="center">
<td align="center">
{{User:Tohline/Math/EQ_ZTFG01}}
{{User:Tohline/Math/EQ_ZTFG01}}
<br />
<font size="-1">
<font size="-1">
Reference (original): [http://adsabs.harvard.edu/abs/1935MNRAS..95..207C Chandraskehar, S. (1935)]
Reference (original): [http://adsabs.harvard.edu/abs/1935MNRAS..95..207C S. Chandraskehar (1935)]
</font>
</font>


Line 35: Line 44:
<tr>
<tr>
<td align="center">
<td align="center">
{{User:Tohline/Math/EQ_Polytrope02}}
<!-- {{User:Tohline/Math/EQ_Polytrope02}}-->
<math>~H = (n+1)K_n \rho^{1/n} = \frac{(n+1)P}{\rho}</math>
</td>
</td>
<td align="center">
<td align="center">
Line 58: Line 68:




In the polytropic relation, the "polytropic index" {{User:Tohline/Math/MP_PolytropicIndex}} and the "polytropic constant" {{User:Tohline/Math/MP_PolytropicConstant}} are assumed to be independent of both {{User:Tohline/Math/VAR_PositionVector01}} and {{User:Tohline/Math/VAR_Time01}}.  In the zero-temperature Fermi gas relation, the two  constants {{User:Tohline/Math/C_FermiPressure}} and {{User:Tohline/Math/C_FermiDensity}} are expressible in terms of various fundamental physical constants, as detailed in the accompanying [http://www.vistrails.org/index.php/User:Tohline/Appendix/Variables_templates variables appendix].  This table also displays (<math>2^{nd}</math> row) the enthalpy as a function of mass density, <math>H(\rho)</math>, and (<math>3^{rd}</math> row) the inverted relation <math>\rho(H)</math> for both barotropic relations, where  
In the polytropic relation, the "polytropic index" {{User:Tohline/Math/MP_PolytropicIndex}} and the "polytropic constant" {{User:Tohline/Math/MP_PolytropicConstant}} are assumed to be independent of both {{User:Tohline/Math/VAR_PositionVector01}} and {{User:Tohline/Math/VAR_Time01}}.  In the zero-temperature Fermi gas relation, the two  constants {{User:Tohline/Math/C_FermiPressure}} and {{User:Tohline/Math/C_FermiDensity}} are expressible in terms of various fundamental physical constants, as detailed in the accompanying [[User:Tohline/Appendix/Variables_templates|variables appendix]].  This table also displays (2<sup>nd</sup> row) the enthalpy as a function of mass density, <math>~H(\rho)</math>, and (3<sup>rd</sup> row) the inverted relation <math>~\rho(H)</math> for both barotropic relations, where  
<div align="center">
<div align="center">
<math>H = \int\frac{dP}{\rho}</math> .
<math>H = \int\frac{dP}{\rho}</math> .
Line 72: Line 82:
</div>
</div>


the zero-temperature Fermi (ZTF) equation of state describes the pressure-density behavior of a nonrelativistic (NR) degenerate electron gas.  We can determine the expression for <math>~P_\mathrm{deg}</math> in this limit by writing the function, <math>~F(\chi)</math>, in terms of two relevant series expansions then keeping only the highest order terms.   
the zero-temperature Fermi (ZTF) equation of state describes the pressure-density behavior of a nonrelativistic (NR) degenerate electron gas.  We can determine the expression for <math>~P_\mathrm{deg}</math> in this limit by writing the function, <math>~F(\chi)</math>, in terms of two relevant series expansions &#8212; for the inverse hyperbolic sine function, see for example, the [http://en.wikipedia.org/wiki/Inverse_hyperbolic_function#Series_expansions Wikipedia presentation] &#8212; then keeping only the highest order terms.   
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">
Line 84: Line 94:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~- 3\chi \biggl(1 - \frac{2}{3}\chi^2 \biggr) \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}+ \cdots  \biggr)
<math>~- 3\chi \biggl(1 - \frac{2}{3}\chi^2 \biggr) \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7} \chi^8+ \cdots  \biggr)
+3 \displaystyle\sum_{m=0}^{\infty} \biggl[ \frac{(-1)^m (2m)!}{2^{2m}(m!)^2} \cdot \frac{\chi^{2m+1}}{(2m + 1)}\biggr]
+3 \displaystyle\sum_{m=0}^{\infty} \biggl[ \frac{(-1)^m (2m)!}{2^{2m}(m!)^2} \cdot \frac{\chi^{2m+1}}{(2m + 1)}\biggr]
</math>
</math>
Line 98: Line 108:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~- 3\chi \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}+ \cdots  \biggr) +
<math>~- 3\chi \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}\chi^8 + \cdots  \biggr) +
2\chi^3 \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}+ \cdots  \biggr)  
2\chi^3 \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}\chi^8 + \cdots  \biggr)
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </td>
  <td align="left">
<math>~
+3\biggl\{\chi - \frac{1}{6}\chi^3 + \frac{3}{2^3 \cdot 5} \chi^5 - \frac{5}{2^4\cdot 7} \chi^7 + \cdots \biggr\}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left"><math>~-3\chi + 3\chi + \chi^3\biggl(-\frac{3}{2} + 2 - \frac{1}{2}  \biggr) + \chi^5 \biggl(\frac{3}{8} +1 +\frac{3^2}{2^3\cdot 5}  \biggr)
+\chi^7 \biggl( -\frac{3}{2^4} -\frac{1}{2^2} - \frac{3\cdot 5}{2^4\cdot 7} \biggr) + \cdots </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left"><math>~\chi^5 \biggl(\frac{15 + 40 + 9}{2^3\cdot 5}  \biggr)
-\chi^7 \biggl( \frac{21 + 28 + 15}{2^4\cdot 7} \biggr) + \cancelto{0}{\cdots} </math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~\approx</math>
  </td>
  <td align="left"><math>~\chi^5 \biggl(\frac{2^3}{5}  \biggr)
-\chi^7 \biggl( \frac{2^2}{7} \biggr) \, .</math>
  </td>
</tr>
</table>
</div>
 
This agrees with, for example, the asymptotic form presented as equation (24) in &sect;X.1 (p. 361) of [<b>[[User:Tohline/Appendix/References#C67|<font color="red">C67</font>]]</b>].  Keeping only the leading term leads to the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~P_\mathrm{deg}\biggr|_\mathrm{NR}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2^3}{5} A_F \biggl( \frac{\rho}{B_F}\biggr)^{5/3}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{2^3}{5} \biggl( \frac{\pi m_e^4 c^5}{3h^3} \biggr) \biggl[ \frac{8\pi m_p}{3} \biggl( \frac{m_e c}{h} \biggr)^3 \mu_e \biggr]^{-5/3} \rho^{5/3}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\mu_e^{-5/3} \biggl[ \frac{2^9 \cdot 3^5 \pi^3}{2^{15}\cdot 3^3\cdot  5^3 \pi^5} \biggr]^{1/3} \biggl( \frac{m_e^4 c^5}{h^3} \biggr)
\biggl( \frac{h}{m_e c} \biggr)^5 \biggl( \frac{\rho}{m_p}\biggr)^{5/3}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2^2 \cdot 5}\biggl( \frac{3}{\pi} \biggr)^{2/3} \biggl( \frac{h^2}{m_e} \biggr) \biggl( \frac{\rho}{m_p \mu_e}\biggr)^{5/3} \, .</math>
  </td>
</tr>
</table>
</div>
 
This matches equation (27) in &sect;X.1 (p. 362) of [<b>[[User:Tohline/Appendix/References#C67|<font color="red">C67</font>]]</b>]; see also, equation (2-32) of [http://adsabs.harvard.edu/abs/1968psen.book.....C D. D. Clayton (1968)], equation (11.42) of [<b>[[User:Tohline/Appendix/References#H87|<font color="red">H87</font>]]</b>], equation (15.23) of [<b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]]</b>], the upper expression labeled as equation (5.163) in [<b>[[User:Tohline/Appendix/References#P00|<font color="red">P00</font>]]</b>], and equations (5.9)-(5.10) of [http://adsabs.harvard.edu/abs/2010asph.book.....C A. R. Choudhuri (2010)].
 
====ZTF Gas in Relativistic Limit====
 
At sufficiently high densities, specifically, when
 
<div align="center">
<math>~\chi \gg 1 \, ,</math>
</div>
 
the zero-temperature Fermi (ZTF) equation of state describes the pressure-density behavior of a degenerate electron gas in the (special) relativistic limit (RL).  We can determine the expression for <math>~P_\mathrm{deg}</math> in this limit by writing the function, <math>~F(\chi)</math>, in terms of two relevant series expansions &#8212; for the inverse hyperbolic sine function, see for example, [http://dlmf.nist.gov/4.38 NIST's ''Digital Library of Mathematical Functions''] &#8212; then keeping only the highest order terms. 
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~F(\chi)</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\biggl(2\chi^4 - 3\chi^2 \biggr) \biggl( 1 + \frac{1}{2\chi^2} - \frac{1}{2^3\chi^4} +\frac{1}{2^4\chi^6} - \frac{5}{2^7\chi^8} + \cdots  \biggr)
+3 \biggl[ \ln(2\chi) + \frac{1}{2^2 \chi^2} - \frac{3}{2^5\chi^4} + \frac{5}{2^5\cdot 3\chi^6} - \cdots \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\biggl( 2\chi^4 + \chi^2 - \frac{1}{2^2} +\frac{1}{2^3\chi^2} - \frac{5}{2^6\chi^4} + \cdots  \biggr)
- \biggl( 3\chi^2 + \frac{3}{2} - \frac{3}{2^3\chi^2} +\frac{3}{2^4\chi^4} + \cdots  \biggr)
</math>
</math>
   </td>
   </td>
Line 112: Line 272:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~+3 \displaystyle\sum_{m=0}^{\infty} \biggl[ \frac{(-1)^m (2m)!}{2^{2m}(m!)^2} \cdot \frac{\chi^{2m+1}}{(2m + 1)}\biggr]
<math>~
+3\biggl\{\chi - \frac{1}{6}\chi^3 +   \biggr\}
+ \biggl[ 3\ln(2\chi) + \frac{3}{2^2 \chi^2} - \frac{3^2}{2^5\chi^4} + \cdots \biggr]
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
2\chi^4 - 2 \chi^2 + 3\ln(2\chi) - \frac{7}{4} +\frac{5}{4\chi^2} - \frac{35}{2^6\chi^4} + \cdots 
</math>
</math>
   </td>
   </td>
Line 119: Line 293:
</table>
</table>
</div>
</div>
where, the series expansion for the inverse hyperbolic sine has been obtained from [http://en.wikipedia.org/wiki/Inverse_hyperbolic_function#Series_expansions Wikipedia's presentation]
 
This agrees with, for example, the asymptotic form presented as equation (25) in &sect;X.1 (p. 361) of [<b>[[User:Tohline/Appendix/References#C67|<font color="red">C67</font>]]</b>].  Keeping only the leading term leads to the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~P_\mathrm{deg}\biggr|_\mathrm{RL}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2 A_F \biggl( \frac{\rho}{B_F}\biggr)^{4/3}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~2 \biggl[ \frac{\pi m_e^4 c^5}{3h^3} \biggr] \biggl[ \frac{8\pi m_p \mu_e}{3} \biggl( \frac{m_e c}{h}\biggr)^3\biggr]^{-4/3} \rho^{4/3}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{1}{2^3}\biggl(\frac{3}{\pi}\biggr)^{1/3} (hc) \biggl(\frac{\rho}{m_p \mu_e}\biggr)^{4/3} \, .</math>
  </td>
</tr>
</table>
</div>
 
This matches equation (28) in &sect;X.1 (p. 362) of [<b>[[User:Tohline/Appendix/References#C67|<font color="red">C67</font>]]</b>]; see also, equation (11.43) of [<b>[[User:Tohline/Appendix/References#H87|<font color="red">H87</font>]]</b>], equation (15.26) of [<b>[[User:Tohline/Appendix/References#KW94|<font color="red">KW94</font>]]</b>], the lower expression labeled as equation (5.163) in [<b>[[User:Tohline/Appendix/References#P00|<font color="red">P00</font>]]</b>], and equations (5.11)-(5.12) of  [http://adsabs.harvard.edu/abs/2010asph.book.....C A. R. Choudhuri (2010)].


===Steady-state Flow-Field Specification===
===Steady-state Flow-Field Specification===


A steady-state velocity flow-field must be specified for time-independent problems.  The specification can be as simple as stating that {{User:Tohline/Math/VAR_VelocityVector01}} = 0 everywhere in space, or that the system has uniform ("solid body") rotation.  Throughout the literature, efforts to generate equilibrium, axisymmetric configurations have adopted a variety of different [[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#SRPtable|''simple rotation profiles'']].  In this H_Book, the flow-field specification will generally vary from chapter to chapter.
A steady-state velocity flow-field must be specified for time-independent problems.  The specification can be as simple as stating that <math>~\vec{v}=0</math> everywhere in space, or that the system has uniform ("solid body") rotation.  Throughout the literature, efforts to generate equilibrium, axisymmetric configurations have adopted a variety of different [[User:Tohline/AxisymmetricConfigurations/SolutionStrategies#SRPtable|''simple rotation profiles'']].  In this H_Book, the flow-field specification will generally vary from chapter to chapter.


==Time-Dependent Problems==
==Time-Dependent Problems==
===Equation of State===
===Equation of State===


For ''time-dependent'' problems we usually will supplement the set of [http://www.vistrails.org/index.php/User:Tohline/PGE principal governing equations] by adopting an equation of state &#8212; ''i.e.,'' a relationship between the state variables {{User:Tohline/Math/VAR_Pressure01}}, {{User:Tohline/Math/VAR_Density01}}, and {{User:Tohline/Math/VAR_Temperature01}} &#8212; that is given by one of the expressions in the following Table, or by [http://www.vistrails.org/index.php/User:Tohline/SR/PressureCombinations some combination of these expressions].  (For example, we could write <math>P_\mathrm{total} = P_\mathrm{gas} + P_\mathrm{deg} + P_\mathrm{rad}</math>.)   
====Components of the Total Pressure====
For ''time-dependent'' problems we usually will supplement the set of [[User:Tohline/PGE#Principal_Governing_Equations|principal governing equations]] by adopting a relationship between the state variables {{User:Tohline/Math/VAR_Pressure01}}, {{User:Tohline/Math/VAR_Density01}}, and {{User:Tohline/Math/VAR_Temperature01}} that is given by one of the expressions in the following Table, or by [[User:Tohline/SR/PressureCombinations#Total_Pressure|some combination of these expressions]].  (For example, we could write <math>~P_\mathrm{total} = P_\mathrm{gas} + P_\mathrm{deg} + P_\mathrm{rad}</math>.)   


<table width="95%" align="center" border=1 cellpadding=5>
<table width="95%" align="center" border=1 cellpadding=5>
<tr>
<tr>
<th colspan=3 align="center">
<th colspan=2 align="center">
<font color="darkblue">Analytic Equations of State</font>
<font color="darkblue">Analytic Equations of State</font>
</th>
</th>
<td align="center" rowspan="2" width="25%"><font color="darkblue">Radiation Pressure</font></td>
</tr>
</tr>
<tr>
<tr>
<td align="center" width="25%"><font color="darkblue">Ideal Gas</font></td>
<td align="center" width="25%"><font color="darkblue">Ideal Gas</font></td>
<td align="center"><font color="darkblue">Degenerate Electron Gas</font></td>
<td align="center"><font color="darkblue">Degenerate Electron Gas</font></td>
<td align="center" width="25%"><font color="darkblue">Radiation</font></td>
</tr>
</tr>
<tr>
<tr>
<td align="center">
<td align="center">
{{User:Tohline/Math/EQ_EOSideal0A}}
{{User:Tohline/Math/EQ_EOSideal0A}}
[http://www.vistrails.org/index.php/User:Tohline/SR/IdealGas [[Image:LSU_OtherFormsButton.jpg|150px]] ]
[[Image:LSU_OtherFormsButton.jpg|150px|link=http://www.vistrails.org/index.php/User:Tohline/SR/IdealGas#Ideal_Gas_Equation_of_State]]
</td>
</td>
<td align="center">
<td align="center">
Line 161: Line 379:
<td colspan=3 align="center">
<td colspan=3 align="center">
{{User:Tohline/Math/EQ_PressureTotal01}}
{{User:Tohline/Math/EQ_PressureTotal01}}
[http://www.vistrails.org/index.php/User:Tohline/SR/PressureCombinations [[Image:LSU_OriginButton.jpg|150px]] ]
[[Image:LSU_OriginButton.jpg|120px|link=http://www.vistrails.org/index.php/User:Tohline/SR/PressureCombinations#Total_Pressure]]
</td>
</td>
</tr>
</tr>
Line 169: Line 387:
In the so-called ''ideal gas'' equation of state, {{User:Tohline/Math/C_GasConstant}} is the gas constant and {{User:Tohline/Math/MP_MeanMolecularWeight}} is the mean molecular weight of the gas. In the equation that gives the electron degeneracy pressure, {{User:Tohline/Math/C_FermiPressure}} is the characteristic Fermi pressure and {{User:Tohline/Math/C_FermiDensity}} is the characteristic Fermi density.  And in the expression for the photon radiation pressure, {{User:Tohline/Math/C_RadiationConstant}} is the radiation constant.  The value of each of these identified physical constants can be found by simply scrolling the computer mouse over the symbol for the constant found in the text of this paragraph, and a definition of each constant can be found in the [[User:Tohline/Appendix/Variables_templates|Variables Appendix]] of this H_Book.
In the so-called ''ideal gas'' equation of state, {{User:Tohline/Math/C_GasConstant}} is the gas constant and {{User:Tohline/Math/MP_MeanMolecularWeight}} is the mean molecular weight of the gas. In the equation that gives the electron degeneracy pressure, {{User:Tohline/Math/C_FermiPressure}} is the characteristic Fermi pressure and {{User:Tohline/Math/C_FermiDensity}} is the characteristic Fermi density.  And in the expression for the photon radiation pressure, {{User:Tohline/Math/C_RadiationConstant}} is the radiation constant.  The value of each of these identified physical constants can be found by simply scrolling the computer mouse over the symbol for the constant found in the text of this paragraph, and a definition of each constant can be found in the [[User:Tohline/Appendix/Variables_templates|Variables Appendix]] of this H_Book.


All three of these equations of state are among the set of [[User:Tohline/Appendix/Equation_templates|key physical equations]] that provide the foundation for our discussion of the '''structure, stability, and dynamics of self-gravitating systems.'''  A discussion of the physical principles that underpin each of these relations can be found in any of a number of different published texts &#8212; see, for example, the set of ''parallel references'' identified in the [[User:Tohline/Appendix/Equation_templates|Equations Appendix]] of this H_Book &#8212; or in the Wiki pages that can be accessed by clicking the linked "other forms" buttons in the above Table.  See also Tassoul (1978) &#8212; specifically the discussion associated with his Chapter 4, Eq. 13 &#8212; for a more general statement related to the proper specification of a supplemental, equation of state relationship.
All three of these equations are among the set of [[User:Tohline/Appendix/Equation_templates#Equations_of_State|key physical equations]] that provide the foundation for our discussion of the '''structure, stability, and dynamics of self-gravitating fluids.'''  A discussion of the physical principles that underpin each of these relations can be found in any of a number of different published texts &#8212; see, for example, the set of ''parallel references'' identified in the [[User:Tohline/Appendix/Equation_templates#Equations_of_State|Equations Appendix]] of this H_Book &#8212; or in the Wiki pages that can be accessed by clicking the linked "other forms" buttons in the above Table.  See also [<b>[[User:Tohline/Appendix/References#T78|<font color="red">T78</font>]]</b>] &#8212; specifically the discussion associated with his Chapter 4, Eq. 13 &#8212; for a more general statement related to the proper specification of a supplemental, equation of state relationship.
 
====The Parameter, &beta;====
It should be pointed out that, in the astrophysics community, the dimensionless quantity <math>~(1-\beta)</math> is sometimes used to denote the relative importance of radiation pressure in a gaseous configuration; specifically,
<div align="center">
<math>~1-\beta = \frac{P_\mathrm{rad}}{P} \, .</math>
</div>
Hence, in the context of our present discussions, the parameter, <math>~\beta</math>, itself is,
<div align="center">
<math>~\beta = \frac{P_\mathrm{gas} + P_\mathrm{deg}}{P} \, .</math>
</div>
Examples include our discussion of [[User:Tohline/SSC/Structure/BiPolytropes/Analytic1.5_3#BiPolytrope_with__and_ne_.3D_3|bipolytropic configurations with]] <math>~(n_c, n_e) = (\tfrac{3}{2}, 3)</math>, as introduced by  [http://adsabs.harvard.edu/abs/1930MNRAS..91....4M E. A. Milne (1930)]; the study of [[User:Tohline/Apps/SMS#Equation_of_State|rotating, supermassive stars]], especially as introduced by [http://adsabs.harvard.edu/abs/1984ApJ...280..825B J. R. Bond, W. D. Arnett, &amp; B. J. Carr (1984)]; and our reference to a [[User:Tohline/SSC/Perturbations#Ledoux_and_Pekeris_.281941.29|derivation of the linear adiabatic wave equation]] by [http://adsabs.harvard.edu/abs/1941ApJ....94..124L P. Ledoux &amp; C. L. Pekeris (1941)].
 
====Adiabatic Exponent====
 
Consider, first, the evolution of a system that is composed entirely of an ideal gas; that is,  <math>~P_\mathrm{rad} = P_\mathrm{deg} = 0</math> and <math>~\beta = 1</math>.  It is widely appreciated that if the entropy of such a system remains constant during a phase of expansion/contraction, then the variation of pressure with density can be properly described by the expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{d\ln P}{d\ln \rho}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\gamma_\mathrm{g} \, ,</math>
  </td>
</tr>
</table>
</div>
where the value of the adiabatic exponent, <math>~\gamma_\mathrm{g}</math>, is given by the ratio of specific heats of the (ideal) gas.  Now, according to the ideal-gas equation of state, changes in the three state variables must always be related to one another via the differential expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~dP\biggr|_\mathrm{gas}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\Re}{\bar{\mu}} \biggl[ \rho dT + T d\rho \biggr]_\mathrm{gas}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ d\ln P</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~d\ln \rho +  d\ln T \, .
</math>
  </td>
</tr>
</table>
</div>
 
We therefore also deduce that,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{d\ln P}{d\ln \rho}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~1 +  \frac{d\ln T}{d\ln\rho}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Rightarrow ~~~ \frac{d\ln T}{d\ln\rho} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\gamma_\mathrm{g}-1  \, ;
</math>
  </td>
</tr>
</table>
</div>
and,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{d\ln P}{d\ln T}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{d\ln \rho}{d\ln T} +1 = \biggl( \frac{1}{\gamma_\mathrm{g} -1}\biggr) + 1
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\gamma_\mathrm{g}}{\gamma_\mathrm{g}-1 } \, .
</math>
  </td>
</tr>
</table>
</div>
 
Following the lead of [[User:Tohline/Appendix/References#C67|[<b><font color="red">C67</font></b>]]], the astrophysics community has found that, when radiation pressure is included in the mix &#8212; that is, when we consider situations in which,
<div align="center">
<math>~P = P_\mathrm{gas} + P_\mathrm{rad}</math>
</div>
&#8212; it can be useful to characterize the adiabatic compression/expansion of fluid elements in terms of three ''separate'' adiabatic exponents, <math>~\Gamma_1, \Gamma_2, \Gamma_3</math>, that are defined via similar differential expressions, namely,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{d\ln P}{d\ln \rho}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\Gamma_1 \, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{d\ln P}{d\ln T}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{\Gamma_2}{\Gamma_2-1} \, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\frac{d\ln T}{d\ln \rho}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\Gamma_3-1 \, .</math>
  </td>
</tr>
</table>
</div>
In this case, though, each of the three adiabatic exponents is a function of <math>~\beta</math> as well as <math>~\gamma_\mathrm{g}</math>; specifically,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\Gamma_1</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\beta + \frac{(4-3\beta)^2 (\gamma_\mathrm{g}-1)}{\beta + 12(\gamma_\mathrm{g}-1)(1-\beta)}
\, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Gamma_2</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1 + \frac{(4-3\beta)(\gamma_\mathrm{g} - 1)}{\beta^2 + 3(\gamma_\mathrm{g} - 1)(1-\beta)(4+\beta)}
\, ;</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\Gamma_3</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
1 + \frac{(4-3\beta)(\gamma_\mathrm{g}-1)}{\beta+12(\gamma_\mathrm{g} - 1)(1-\beta)}
\, .</math>
  </td>
</tr>
</table>
</div>
 
===Entropy Tracer===
We begin with the basic equation of state,
<div align="center">
<math>~P = (\gamma_g - 1)\epsilon\rho \, ,</math>
</div>
and the 1<sup>st</sup> Law of Thermodynamics,
<div align="center">
{{ User:Tohline/Math/EQ_FirstLaw01 }}
</div>
Adopting the concept of an ''entropy tracer'',
<div align="center">
<math>~\tau \equiv (\epsilon\rho)^{1/\gamma_g} \, ,</math>
</div>
the 1<sup>st</sup> becomes,
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\frac{\tau}{c_p} \frac{ds}{dt}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{\partial\tau}{\partial t} + \nabla\cdot (\tau \vec{v}) \, .</math>
  </td>
</tr>
</table>
See equation (33) of Motl et al. (2002)
 


===Initial Conditions===
===Initial Conditions===
For ''time-dependent'' problems, the [http://www.vistrails.org/index.php/User:Tohline/PGE principal governing equations] must be supplemented further through the specification of initial conditions.  Frequently throughout this H_Book, we will select as initial conditions a specification of {{User:Tohline/Math/VAR_Density01}}( {{User:Tohline/Math/VAR_PositionVector01}}, {{User:Tohline/Math/VAR_Time01}} = 0), {{User:Tohline/Math/VAR_Pressure01}}( {{User:Tohline/Math/VAR_PositionVector01}}, {{User:Tohline/Math/VAR_Time01}} = 0), and {{User:Tohline/Math/VAR_VelocityVector01}}( {{User:Tohline/Math/VAR_PositionVector01}}, {{User:Tohline/Math/VAR_Time01}} = 0) that, as a group themselves, define a static or steady-state equilibrium '''structure'''.  Perturbation or computational fluid dynamic (CFD) techniques can be used to test the '''stability''' or nonlinear '''dynamical behavior''' of such structures.
For ''time-dependent'' problems, the [[User:Tohline/PGE#Principal_Governing_Equations|principal governing equations]] must be supplemented further through the specification of initial conditions.  Frequently throughout this H_Book, we will select as initial conditions a specification of <math>~\rho(\vec{x}, t=0)</math>, <math>~P(\vec{x}, t=0)</math>, and <math>~\vec{v}(\vec{x}, t=0)</math> that, as a group themselves, define a static or steady-state equilibrium '''structure'''.  Perturbation or computational fluid dynamic (CFD) techniques can be used to test the '''stability''' or nonlinear '''dynamical behavior''' of such structures.
 
=See Also=
For derivations and additional discussions of the triplet of adiabatic exponents, <math>~\Gamma_1, \Gamma_2, \Gamma_3</math>, see:
* [[User:Tohline/Appendix/References#C67|[<b><font color="red">C67</font></b>]]] &#8212; especially pp. 56-58; note that <math>~d\ln V = - d\ln \rho</math>;
* [http://adsabs.harvard.edu/abs/1968psen.book.....C Donald D. Clayton (1968)]; especially p. 118; note that <math>~d\ln V = - d\ln \rho</math>;
* [[User:Tohline/Appendix/References#KW94|[<b><font color="red">KW94</font></b>]]] &#8212; especially p. 106;
* [http://adsabs.harvard.edu/abs/1994sipp.book.....H Carl J. Hansen &amp; Steven D. Kawaler (1994)] &#8212; especially pp. 136-138;
* Our own brief presentation of [[User:Tohline/Apps/SMS#Expression_for_Adiabatic_Exponent|several equivalent expressions for <math>~\Gamma_1</math>]].
 


{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 15:26, 9 August 2019

Supplemental Relations

"The equilibrium configuration of a stellar system depends crucially on the equation of state. Similarly, the time evolution of the perturbation will depend on the relationship that is assumed between the perturbed values of pressure, temperature, and density. The simplest assumption regarding these perturbations will be to involve an adiabatic relationship between perturbed pressure and density."

— Drawn from §7, p. 172 of [P00]

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

Examining the set of principal governing equations that describe adiabatic flows, we see that — apart from the independent variables <math>~t</math> and <math>~\vec{x}</math> — the equations involve the vector velocity <math>~\vec{v}</math>, and the four scalar variables, <math>~\Phi</math>, <math>~P</math>, <math>~\rho</math>, and <math>~\epsilon</math>. Because the variables outnumber the equations by one, a supplemental relationship between the physical variables must be specified in order to close the set of equations. (If nonadiabatic flows are considered, additional supplemental relations must be specified because the scalar variables <math>~T</math> and <math>~s</math> enter the discussion as well.) Also, in order to complete the unique specification of a particular physical problem, either a steady-state flow field or initial conditions must be specified, depending on whether one is studying a time-independent (structure) or time-dependent (stability or dynamics) problem, respectively. Throughout this H_Book, the following strategy will be adopted in order to complete the physical specification of each examined system: (A) For time-independent problems, we will adopt a structural relationship between <math>~P</math> and <math>~\rho</math>, and specify a steady-state flow-field. (B) For time-dependent problems, we will adopt an equation of state, and specify initial conditions.


Time-Independent Problems

Barotropic Structure

For time-independent problems, a structural relationship between <math>~P</math> and <math>~\rho</math> is required to close the system of principal governing equations. [T78] refers to this as a "geometrical" rather than a "structural" relationship; see the discussion associated with his Chapter 4, Eq. 14. Generally throughout this H_Book, we will assume that all time-independent configurations can be described as barotropic structures; that is, we will assume that <math>~P</math> is only a function of <math>~\rho</math> throughout such structures. (The Poincaré-Wavre theorem provides additional support for, as well as additional implications of, this assumption.) More specifically, we generally will adopt one of the two analytically prescribable <math>~P(\rho)</math> relationships displayed in the first row of the following Table.

Barotropic Relations

Polytropic Zero-temperature Fermi (degenerate electron) Gas

<math>~P = K_\mathrm{n} \rho^{1+1/n}</math>

LSU Key.png

<math>~P_\mathrm{deg} = A_\mathrm{F} F(\chi) </math>

where:  <math>F(\chi) \equiv \chi(2\chi^2 - 3)(\chi^2 + 1)^{1/2} + 3\sinh^{-1}\chi</math>

and:   

<math>\chi \equiv (\rho/B_\mathrm{F})^{1/3}</math>

Reference (original): S. Chandraskehar (1935)

<math>~H = (n+1)K_n \rho^{1/n} = \frac{(n+1)P}{\rho}</math>

<math> H = \frac{8A_\mathrm{F}}{B_\mathrm{F}} \biggl[(\chi^2 + 1)^{1/2} - 1 \biggr] </math>

<math>~\rho = \biggl[ \frac{H}{(n+1)K_\mathrm{n}} \biggr]^n </math>

<math> \rho = B_\mathrm{F} \biggl[\biggl(\frac{HB_\mathrm{F}}{8A_\mathrm{F}} + 1 \biggr)^{2} - 1 \biggr]^{3/2} </math>


In the polytropic relation, the "polytropic index" <math>~n</math> and the "polytropic constant" <math>~K_\mathrm{n}</math> are assumed to be independent of both <math>~\vec{x}</math> and <math>~t</math>. In the zero-temperature Fermi gas relation, the two constants <math>~A_\mathrm{F}</math> and <math>~B_\mathrm{F}</math> are expressible in terms of various fundamental physical constants, as detailed in the accompanying variables appendix. This table also displays (2nd row) the enthalpy as a function of mass density, <math>~H(\rho)</math>, and (3rd row) the inverted relation <math>~\rho(H)</math> for both barotropic relations, where

<math>H = \int\frac{dP}{\rho}</math> .

In both cases, we have chosen an integration constant such that <math>~H</math> is zero when <math>~\rho</math> is zero.

Nonrelativistic ZTF Gas

At sufficiently low densities, specifically, when

<math>~\chi \ll 1 \, ,</math>

the zero-temperature Fermi (ZTF) equation of state describes the pressure-density behavior of a nonrelativistic (NR) degenerate electron gas. We can determine the expression for <math>~P_\mathrm{deg}</math> in this limit by writing the function, <math>~F(\chi)</math>, in terms of two relevant series expansions — for the inverse hyperbolic sine function, see for example, the Wikipedia presentation — then keeping only the highest order terms.

<math>~F(\chi)</math>

<math>~=</math>

<math>~- 3\chi \biggl(1 - \frac{2}{3}\chi^2 \biggr) \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7} \chi^8+ \cdots \biggr) +3 \displaystyle\sum_{m=0}^{\infty} \biggl[ \frac{(-1)^m (2m)!}{2^{2m}(m!)^2} \cdot \frac{\chi^{2m+1}}{(2m + 1)}\biggr] </math>

 

<math>~=</math>

<math>~- 3\chi \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}\chi^8 + \cdots \biggr) + 2\chi^3 \biggl( 1 + \frac{1}{2}\chi^2 - \frac{1}{2^3}\chi^4 +\frac{1}{2^4}\chi^6 - \frac{5}{2^7}\chi^8 + \cdots \biggr) </math>

 

 

<math>~ +3\biggl\{\chi - \frac{1}{6}\chi^3 + \frac{3}{2^3 \cdot 5} \chi^5 - \frac{5}{2^4\cdot 7} \chi^7 + \cdots \biggr\} </math>

 

<math>~=</math>

<math>~-3\chi + 3\chi + \chi^3\biggl(-\frac{3}{2} + 2 - \frac{1}{2} \biggr) + \chi^5 \biggl(\frac{3}{8} +1 +\frac{3^2}{2^3\cdot 5} \biggr)

+\chi^7 \biggl( -\frac{3}{2^4} -\frac{1}{2^2} - \frac{3\cdot 5}{2^4\cdot 7} \biggr) + \cdots </math>

 

<math>~=</math>

<math>~\chi^5 \biggl(\frac{15 + 40 + 9}{2^3\cdot 5} \biggr)

-\chi^7 \biggl( \frac{21 + 28 + 15}{2^4\cdot 7} \biggr) + \cancelto{0}{\cdots} </math>

 

<math>~\approx</math>

<math>~\chi^5 \biggl(\frac{2^3}{5} \biggr)

-\chi^7 \biggl( \frac{2^2}{7} \biggr) \, .</math>

This agrees with, for example, the asymptotic form presented as equation (24) in §X.1 (p. 361) of [C67]. Keeping only the leading term leads to the expression,

<math>~P_\mathrm{deg}\biggr|_\mathrm{NR}</math>

<math>~=</math>

<math>~\frac{2^3}{5} A_F \biggl( \frac{\rho}{B_F}\biggr)^{5/3}</math>

 

<math>~=</math>

<math>~\frac{2^3}{5} \biggl( \frac{\pi m_e^4 c^5}{3h^3} \biggr) \biggl[ \frac{8\pi m_p}{3} \biggl( \frac{m_e c}{h} \biggr)^3 \mu_e \biggr]^{-5/3} \rho^{5/3}</math>

 

<math>~=</math>

<math>~\mu_e^{-5/3} \biggl[ \frac{2^9 \cdot 3^5 \pi^3}{2^{15}\cdot 3^3\cdot 5^3 \pi^5} \biggr]^{1/3} \biggl( \frac{m_e^4 c^5}{h^3} \biggr) \biggl( \frac{h}{m_e c} \biggr)^5 \biggl( \frac{\rho}{m_p}\biggr)^{5/3}</math>

 

<math>~=</math>

<math>~\frac{1}{2^2 \cdot 5}\biggl( \frac{3}{\pi} \biggr)^{2/3} \biggl( \frac{h^2}{m_e} \biggr) \biggl( \frac{\rho}{m_p \mu_e}\biggr)^{5/3} \, .</math>

This matches equation (27) in §X.1 (p. 362) of [C67]; see also, equation (2-32) of D. D. Clayton (1968), equation (11.42) of [H87], equation (15.23) of [KW94], the upper expression labeled as equation (5.163) in [P00], and equations (5.9)-(5.10) of A. R. Choudhuri (2010).

ZTF Gas in Relativistic Limit

At sufficiently high densities, specifically, when

<math>~\chi \gg 1 \, ,</math>

the zero-temperature Fermi (ZTF) equation of state describes the pressure-density behavior of a degenerate electron gas in the (special) relativistic limit (RL). We can determine the expression for <math>~P_\mathrm{deg}</math> in this limit by writing the function, <math>~F(\chi)</math>, in terms of two relevant series expansions — for the inverse hyperbolic sine function, see for example, NIST's Digital Library of Mathematical Functions — then keeping only the highest order terms.

<math>~F(\chi)</math>

<math>~=</math>

<math>~\biggl(2\chi^4 - 3\chi^2 \biggr) \biggl( 1 + \frac{1}{2\chi^2} - \frac{1}{2^3\chi^4} +\frac{1}{2^4\chi^6} - \frac{5}{2^7\chi^8} + \cdots \biggr) +3 \biggl[ \ln(2\chi) + \frac{1}{2^2 \chi^2} - \frac{3}{2^5\chi^4} + \frac{5}{2^5\cdot 3\chi^6} - \cdots \biggr] </math>

 

<math>~=</math>

<math>~ \biggl( 2\chi^4 + \chi^2 - \frac{1}{2^2} +\frac{1}{2^3\chi^2} - \frac{5}{2^6\chi^4} + \cdots \biggr) - \biggl( 3\chi^2 + \frac{3}{2} - \frac{3}{2^3\chi^2} +\frac{3}{2^4\chi^4} + \cdots \biggr) </math>

 

 

<math>~ + \biggl[ 3\ln(2\chi) + \frac{3}{2^2 \chi^2} - \frac{3^2}{2^5\chi^4} + \cdots \biggr] </math>

 

<math>~=</math>

<math>~ 2\chi^4 - 2 \chi^2 + 3\ln(2\chi) - \frac{7}{4} +\frac{5}{4\chi^2} - \frac{35}{2^6\chi^4} + \cdots </math>

This agrees with, for example, the asymptotic form presented as equation (25) in §X.1 (p. 361) of [C67]. Keeping only the leading term leads to the expression,

<math>~P_\mathrm{deg}\biggr|_\mathrm{RL}</math>

<math>~=</math>

<math>~2 A_F \biggl( \frac{\rho}{B_F}\biggr)^{4/3}</math>

 

<math>~=</math>

<math>~2 \biggl[ \frac{\pi m_e^4 c^5}{3h^3} \biggr] \biggl[ \frac{8\pi m_p \mu_e}{3} \biggl( \frac{m_e c}{h}\biggr)^3\biggr]^{-4/3} \rho^{4/3}</math>

 

<math>~=</math>

<math>~\frac{1}{2^3}\biggl(\frac{3}{\pi}\biggr)^{1/3} (hc) \biggl(\frac{\rho}{m_p \mu_e}\biggr)^{4/3} \, .</math>

This matches equation (28) in §X.1 (p. 362) of [C67]; see also, equation (11.43) of [H87], equation (15.26) of [KW94], the lower expression labeled as equation (5.163) in [P00], and equations (5.11)-(5.12) of A. R. Choudhuri (2010).

Steady-state Flow-Field Specification

A steady-state velocity flow-field must be specified for time-independent problems. The specification can be as simple as stating that <math>~\vec{v}=0</math> everywhere in space, or that the system has uniform ("solid body") rotation. Throughout the literature, efforts to generate equilibrium, axisymmetric configurations have adopted a variety of different simple rotation profiles. In this H_Book, the flow-field specification will generally vary from chapter to chapter.

Time-Dependent Problems

Equation of State

Components of the Total Pressure

For time-dependent problems we usually will supplement the set of principal governing equations by adopting a relationship between the state variables <math>~P</math>, <math>~\rho</math>, and <math>~T</math> that is given by one of the expressions in the following Table, or by some combination of these expressions. (For example, we could write <math>~P_\mathrm{total} = P_\mathrm{gas} + P_\mathrm{deg} + P_\mathrm{rad}</math>.)

Analytic Equations of State

Radiation Pressure
Ideal Gas Degenerate Electron Gas

LSU Key.png

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

LSU OtherFormsButton.jpg

LSU Key.png

<math>~P_\mathrm{deg} = A_\mathrm{F} F(\chi) </math>

where:  <math>F(\chi) \equiv \chi(2\chi^2 - 3)(\chi^2 + 1)^{1/2} + 3\sinh^{-1}\chi</math>

and:   

<math>\chi \equiv (\rho/B_\mathrm{F})^{1/3}</math>

LSU Key.png

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

Normalized Total Pressure:

LSU Key.png

<math>~p_\mathrm{total} = \biggl(\frac{\mu_e m_p}{\bar{\mu} m_u} \biggr) 8 \chi^3 \frac{T}{T_e} + F(\chi) + \frac{8\pi^4}{15} \biggl( \frac{T}{T_e} \biggr)^4</math>

LSU OriginButton.jpg

In the so-called ideal gas equation of state, <math>~\Re</math> is the gas constant and <math>~\bar{\mu}</math> is the mean molecular weight of the gas. In the equation that gives the electron degeneracy pressure, <math>~A_\mathrm{F}</math> is the characteristic Fermi pressure and <math>~B_\mathrm{F}</math> is the characteristic Fermi density. And in the expression for the photon radiation pressure, <math>~a_\mathrm{rad}</math> is the radiation constant. The value of each of these identified physical constants can be found by simply scrolling the computer mouse over the symbol for the constant found in the text of this paragraph, and a definition of each constant can be found in the Variables Appendix of this H_Book.

All three of these equations are among the set of key physical equations that provide the foundation for our discussion of the structure, stability, and dynamics of self-gravitating fluids. A discussion of the physical principles that underpin each of these relations can be found in any of a number of different published texts — see, for example, the set of parallel references identified in the Equations Appendix of this H_Book — or in the Wiki pages that can be accessed by clicking the linked "other forms" buttons in the above Table. See also [T78] — specifically the discussion associated with his Chapter 4, Eq. 13 — for a more general statement related to the proper specification of a supplemental, equation of state relationship.

The Parameter, β

It should be pointed out that, in the astrophysics community, the dimensionless quantity <math>~(1-\beta)</math> is sometimes used to denote the relative importance of radiation pressure in a gaseous configuration; specifically,

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

Hence, in the context of our present discussions, the parameter, <math>~\beta</math>, itself is,

<math>~\beta = \frac{P_\mathrm{gas} + P_\mathrm{deg}}{P} \, .</math>

Examples include our discussion of bipolytropic configurations with <math>~(n_c, n_e) = (\tfrac{3}{2}, 3)</math>, as introduced by E. A. Milne (1930); the study of rotating, supermassive stars, especially as introduced by J. R. Bond, W. D. Arnett, & B. J. Carr (1984); and our reference to a derivation of the linear adiabatic wave equation by P. Ledoux & C. L. Pekeris (1941).

Adiabatic Exponent

Consider, first, the evolution of a system that is composed entirely of an ideal gas; that is, <math>~P_\mathrm{rad} = P_\mathrm{deg} = 0</math> and <math>~\beta = 1</math>. It is widely appreciated that if the entropy of such a system remains constant during a phase of expansion/contraction, then the variation of pressure with density can be properly described by the expression,

<math>~\frac{d\ln P}{d\ln \rho}</math>

<math>~=</math>

<math>~\gamma_\mathrm{g} \, ,</math>

where the value of the adiabatic exponent, <math>~\gamma_\mathrm{g}</math>, is given by the ratio of specific heats of the (ideal) gas. Now, according to the ideal-gas equation of state, changes in the three state variables must always be related to one another via the differential expression,

<math>~dP\biggr|_\mathrm{gas}</math>

<math>~=</math>

<math>~\frac{\Re}{\bar{\mu}} \biggl[ \rho dT + T d\rho \biggr]_\mathrm{gas} </math>

<math>~\Rightarrow ~~~ d\ln P</math>

<math>~=</math>

<math>~d\ln \rho + d\ln T \, . </math>

We therefore also deduce that,

<math>~\frac{d\ln P}{d\ln \rho}</math>

<math>~=</math>

<math>~1 + \frac{d\ln T}{d\ln\rho} </math>

<math>~\Rightarrow ~~~ \frac{d\ln T}{d\ln\rho} </math>

<math>~=</math>

<math>~\gamma_\mathrm{g}-1 \, ; </math>

and,

<math>~\frac{d\ln P}{d\ln T}</math>

<math>~=</math>

<math>~\frac{d\ln \rho}{d\ln T} +1 = \biggl( \frac{1}{\gamma_\mathrm{g} -1}\biggr) + 1 </math>

 

<math>~=</math>

<math>~\frac{\gamma_\mathrm{g}}{\gamma_\mathrm{g}-1 } \, . </math>

Following the lead of [C67], the astrophysics community has found that, when radiation pressure is included in the mix — that is, when we consider situations in which,

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

— it can be useful to characterize the adiabatic compression/expansion of fluid elements in terms of three separate adiabatic exponents, <math>~\Gamma_1, \Gamma_2, \Gamma_3</math>, that are defined via similar differential expressions, namely,

<math>~\frac{d\ln P}{d\ln \rho}</math>

<math>~=</math>

<math>~\Gamma_1 \, ;</math>

<math>~\frac{d\ln P}{d\ln T}</math>

<math>~=</math>

<math>~\frac{\Gamma_2}{\Gamma_2-1} \, ;</math>

<math>~\frac{d\ln T}{d\ln \rho}</math>

<math>~=</math>

<math>~\Gamma_3-1 \, .</math>

In this case, though, each of the three adiabatic exponents is a function of <math>~\beta</math> as well as <math>~\gamma_\mathrm{g}</math>; specifically,

<math>~\Gamma_1</math>

<math>~=</math>

<math>~ \beta + \frac{(4-3\beta)^2 (\gamma_\mathrm{g}-1)}{\beta + 12(\gamma_\mathrm{g}-1)(1-\beta)} \, ;</math>

<math>~\Gamma_2</math>

<math>~=</math>

<math>~ 1 + \frac{(4-3\beta)(\gamma_\mathrm{g} - 1)}{\beta^2 + 3(\gamma_\mathrm{g} - 1)(1-\beta)(4+\beta)} \, ;</math>

<math>~\Gamma_3</math>

<math>~=</math>

<math>~ 1 + \frac{(4-3\beta)(\gamma_\mathrm{g}-1)}{\beta+12(\gamma_\mathrm{g} - 1)(1-\beta)} \, .</math>

Entropy Tracer

We begin with the basic equation of state,

<math>~P = (\gamma_g - 1)\epsilon\rho \, ,</math>

and the 1st Law of Thermodynamics,

LSU Key.png

<math>T \frac{ds}{dt} = \frac{d\epsilon}{dt} + P \frac{d}{dt} \biggl(\frac{1}{\rho}\biggr)</math>

Adopting the concept of an entropy tracer,

<math>~\tau \equiv (\epsilon\rho)^{1/\gamma_g} \, ,</math>

the 1st becomes,

<math>~\frac{\tau}{c_p} \frac{ds}{dt}</math>

<math>~=</math>

<math>~ \frac{\partial\tau}{\partial t} + \nabla\cdot (\tau \vec{v}) \, .</math>

See equation (33) of Motl et al. (2002)


Initial Conditions

For time-dependent problems, the principal governing equations must be supplemented further through the specification of initial conditions. Frequently throughout this H_Book, we will select as initial conditions a specification of <math>~\rho(\vec{x}, t=0)</math>, <math>~P(\vec{x}, t=0)</math>, and <math>~\vec{v}(\vec{x}, t=0)</math> that, as a group themselves, define a static or steady-state equilibrium structure. Perturbation or computational fluid dynamic (CFD) techniques can be used to test the stability or nonlinear dynamical behavior of such structures.

See Also

For derivations and additional discussions of the triplet of adiabatic exponents, <math>~\Gamma_1, \Gamma_2, \Gamma_3</math>, see:


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