Difference between revisions of "User:Tohline/SSC/VirialEquilibrium/UniformDensity"

From VistrailsWiki
Jump to navigation Jump to search
(Begin discussion of uniform-density virial equilibria)
 
 
(9 intermediate revisions by the same user not shown)
Line 1: Line 1:
<!-- __FORCETOC__ will force the creation of a Table of Contents -->
__FORCETOC__ <!-- will force the creation of a Table of Contents -->
<!-- __NOTOC__ will force TOC off -->
<!-- __NOTOC__ will force TOC off -->
=Isolated Uniform-Density Sphere=
=Uniform-Density Sphere=
{{LSU_HBook_header}}
{{LSU_HBook_header}}
==Review==
==Review==
Line 90: Line 90:
Once the pressure exerted by the external medium (<math>~P_e</math>), and the configuration's mass (<math>~M_\mathrm{tot}</math>), angular momentum (<math>~J</math>), and specific entropy (via <math>~K</math>) &#8212; or, in the isothermal case, sound speed (<math>~c_s</math>) &#8212;  have been specified, the values of all of the coefficients are known and <math>~\chi_\mathrm{eq}</math> can be determined.
Once the pressure exerted by the external medium (<math>~P_e</math>), and the configuration's mass (<math>~M_\mathrm{tot}</math>), angular momentum (<math>~J</math>), and specific entropy (via <math>~K</math>) &#8212; or, in the isothermal case, sound speed (<math>~c_s</math>) &#8212;  have been specified, the values of all of the coefficients are known and <math>~\chi_\mathrm{eq}</math> can be determined.


==Specific Case Being Considered Here==
==Adiabatic Evolution of an Isolated Sphere==
Here we seek to determine the equilibrium radius of a non-rotating configuration (<math>~J = 0</math>) that undergoes adiabatic compression/expansion (<math>\delta_{1\gamma_g} =~0</math>) and that is not confined by an external medium (<math>P_e = 0~</math>). In this case, the statement of virial equilibrium is simplified considerably.  Specifically, <math>~\chi_\mathrm{eq}</math> is given by the root(s) of the equation,
Here we seek to determine the equilibrium radius of a non-rotating configuration (<math>~J = 0</math>) that undergoes adiabatic compression/expansion (<math>\delta_{1\gamma_g} =~0</math>) and that is not confined by an external medium (<math>P_e = 0~</math>).  
===Solution===
In this case, the statement of virial equilibrium is simplified considerably.  Specifically, <math>~\chi_\mathrm{eq}</math> is given by the root(s) of the equation,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~A\chi_\mathrm{eq}^{-1} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~3(\gamma_g-1) B\chi_\mathrm{eq}^{3 -3\gamma_g} </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\Rightarrow ~~~~~\chi_\mathrm{eq}^{3\gamma_g-4} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{3(\gamma_g-1) B}{A}
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
<math>
3(\gamma_g-1) B\chi^{3 -3\gamma_g~-~A\chi^{-1}   = 0 \, .
\biggl[ 3K M_\mathrm{tot} \biggl( \frac{3M_\mathrm{tot} }{4\pi R_0^3} \biggr)^{\gamma_g - 1}  \cdot \mathfrak{f}_A \biggr]
\biggl[ \frac{3}{5} \frac{GM_\mathrm{tot} ^2}{R_0} \cdot \mathfrak{f}_W \biggr] ^{-1}
</math>
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl( \frac{1}{R_0} \biggr)^{3\gamma_g-4} \biggl[ 5\biggl( \frac{3}{4\pi} \biggr)^{\gamma_g-1} \biggr(\frac{K}{G}\biggr)
M^{(\gamma_g-2)} \cdot \frac{\mathfrak{f}_A}{\mathfrak{f}_W} \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>
</div>
In other words,
In other words,
<div align="center">
<div align="center">
<math>
<math>
R_\mathrm{eq} = R_0 \chi_\mathrm{eq} = \biggl[ \frac{3(\gamma_g-1) B}{A} \cdot R_0^{(3\gamma_g-4)} \biggr]^{1/(3\gamma_g-4)} = \biggl[ 5\biggl( \frac{3}{4\pi} \biggr)^{\gamma_g-1} \cdot \frac{KM^{(\gamma_g-2)}}{G} \biggr]^{1/(3\gamma_g-4)} \, .
R_\mathrm{eq} = \biggl[ 5\biggl( \frac{3}{4\pi} \biggr)^{\gamma_g-1} \biggr(\frac{K}{G}\biggr)  
M^{(\gamma_g-2)} \cdot \frac{\mathfrak{f}_A}{\mathfrak{f}_W} \biggr]^{1/(3\gamma_g-4)} \, .
</math>
</math>
</div>
</div>
Accordingly, the equilibrium mass-radius relationship for adiabatic configurations of a given specific entropy is,
 
===Comparison with Detailed Force-Balance Model===
This derived solution will look more familiar if, instead of <math>~K</math>, we express the solution in terms of the central pressure,
<div align="center">
<math>P_c = K\rho_0^{\gamma_g} \, ,</math>
</div>
where, for this uniform-density sphere, <math>~\rho_0 = 3M_\mathrm{tot}/(4\pi R_\mathrm{eq}^3)</math>.  Hence,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~K</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>P_c \biggl( \frac{4\pi R_\mathrm{eq}^3}{3M_\mathrm{tot}} \biggr)^{\gamma_g} \, ,</math>
  </td>
</tr>
</table>
</div>
and the solution takes the form,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>R_\mathrm{eq}^{3\gamma_g - 4}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>5\biggl( \frac{4\pi}{3} \biggr) \biggr(\frac{P_c R_\mathrm{eq}^{3\gamma_g}}{GM^2_\mathrm{tot}}\biggr)
\cdot \frac{\mathfrak{f}_A}{\mathfrak{f}_W} </math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>\Rightarrow ~~~~ R_\mathrm{eq}^{4} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{3}{20\pi} \biggr) \biggr(\frac{GM^2_\mathrm{tot}}{P_c}\biggr)
\cdot \frac{\mathfrak{f}_W}{\mathfrak{f}_A} \, .</math>
  </td>
</tr>
</table>
</div>
Or, solving for the central pressure,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~P_c</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{3}{20\pi} \cdot \frac{\mathfrak{f}_W}{\mathfrak{f}_A} \biggr) \frac{GM^2_\mathrm{tot}}{R_\mathrm{eq}^{4} }
\, .</math>
  </td>
</tr>
</table>
</div>
This should be compared with our [[User:Tohline/SSC/Structure/UniformDensity#Summary|detailed force-balance solution]] of the interior structure of an isolated, nonrotating, uniform-density sphere, which gives the precise expression,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~P_c</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{3}{8\pi}  \biggr) \frac{GM^2_\mathrm{tot}}{R_\mathrm{eq}^{4} }
\, .</math>
  </td>
</tr>
</table>
</div>
The expression for <math>~P_c</math> derived from our identification of an extremum in the free energy is identical to the expression derived from the more precise, detailed force-balance analysis, except that the leading numerical coefficients differ by a factor of <math>~(5\mathfrak{f}_A/2\mathfrak{f}_W)</math>. 
 
From a free-energy analysis alone, the best we can do is assume that both structural form factors, <math>~\mathfrak{f}_W</math> and <math>\mathfrak{f}_A</math>, are of order unity.  But knowing the detailed force-balance solution allows us to evaluate both form factors.  From our [[User:Tohline/SphericallySymmetricConfigurations/Virial#FormFactors|introductory discussion of the free energy function]], their respective definitions are,
<div align="center">
<table border="0" cellpadding="5" align="center">
 
<tr>
  <td align="right">
<math>~\mathfrak{f}_W</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~  3\cdot 5 \int_0^1 \biggl\{ \int_0^x  \biggl[ \frac{\rho(x)}{\rho_c}\biggr] x^2 dx \biggr\}  \biggl[ \frac{\rho(x)}{\rho_c}\biggr] x dx\, ,</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~\mathfrak{f}_A</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~ 3\int_0^1 \biggl[ \frac{P(x)}{P_c}\biggr]  x^2 dx \, .</math>
  </td>
</tr>
</table>
</div>
Now, because the configuration under discussion has a uniform density, we should set <math>~\rho(x)/\rho_c = 1</math> in the definition of <math>~\mathfrak{f}_W</math> which, after evaluation of the nested integrals, gives <math>~\mathfrak{f}_W = 1</math>.  But, instead of being uniform throughout the configuration, in the [[User:Tohline/SSC/Structure/UniformDensity#Summary|detailed force-balance model]], the pressure drops from the center to the surface according to the relation,
<div align="center">
<math>\frac{P(x)}{P_c} = 1 - x^2 \, .</math>
</div>
Integrating over this function, in accordance with the definition of <math>~\mathfrak{f}_A</math>, gives,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\mathfrak{f}_A</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~ 3\int_0^1 (1-x^2) x^2 dx = 3 \biggl[ \frac{x^3}{3} - \frac{x^5}{5} \biggr]_0^1 = \frac{2}{5} \, .</math>
  </td>
</tr>
</table>
</div>
Hence, the ratio,
<div align="center">
<div align="center">
<math>
<math>
M^{(\gamma_g - 2)} \propto R_\mathrm{eq}^{(3\gamma_g -4)} \, .
\frac{5\mathfrak{f}_A}{2\mathfrak{f}_W} = 1 \, ,
</math>
</math>
</div>
</div>
Notice that, for <math>\gamma_g=2</math>, the equilibrium radius depends only on the specific entropy of the gas and is independent of the configuration's mass.  Conversely, notice that, for <math>\gamma_g = 4/3</math>, the mass of the configuration is independent of the radiusFor <math>\gamma_g</math> &#x3E; <math> 2</math> or <math>\gamma_g </math>&#x3C; <math>4/3</math>, configurations with larger mass (but the same specific entropy) have larger equilibrium radiiHowever, for <math>\gamma_g</math> in the range, <math>2</math> &#x3E; <math>\gamma_g </math> &#x3E; <math>4/3</math>, configurations with larger mass have smaller equilibrium radii.  Note that the result obtained for the isothermal configuration could have been obtained by setting <math>\gamma_g = 1</math> in this adiabatic solution, because <math>K = c_s^2</math> when  <math>\gamma_g = 1</math>.
which brings into perfect agreement the two separate determinations of the equilibrium expressions for <math>~R_\mathrm{eq}</math> and <math>~P_c</math> in terms of one another and the total mass
 
This demonstrates that the free-energy approach to determining the equilibrium radius of a spherical configuration is only handicapped by its inability to precisely nail down values of the structural form factorsBut this is not a severe limitation as the (dimensionless) form factors are generally of order unityIn contrast, the free-energy analysis brings with it a capability to readily evaluate the global stability of equilibrium configurations.  
 


It is also instructive to write the coefficient <math>B</math> in terms of the average sound speed as defined above. In this case,
==Adiabatic Evolution of Pressure-truncated Sphere==
Here we seek to determine the equilibrium radius of a non-rotating configuration (<math>~J = 0</math>) that undergoes adiabatic compression/expansion (<math>\delta_{1\gamma_g} =~0</math>) and that is embedded in a hot, tenuous external medium whose confining pressure, <math>~P_e</math>, truncates the configuration.  
===Solution===
In this case, virial equilibrium implies that <math>~\chi_\mathrm{eq}</math> is given by the root(s) of the equation,
<div align="center">
<div align="center">
<math>
<math>
R_\mathrm{eq} = R_0 \biggl[ \frac{GM}{5 \bar{c_s}^2 R_0} \biggr]^{1/(4- 3\gamma_g)} \, ,
3(\gamma_g-1) B\chi^{3 -3\gamma_g} ~ -~A\chi^{-1-~ 3D\chi^3 = 0 \, .
</math>
</math>
</div>
</div>
so the equilibrium radius of an isolated, nonrotating, uniform density, adiabatic sphere is,
Hence,
<div align="center">
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~D</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
<math>
R_\mathrm{eq} = R_0 = \frac{GM}{5 \bar{c_s}^2 } \, .
(\gamma_g-1) B\chi_\mathrm{eq}^{-3\gamma_g} -  \frac{A}{3}\chi_\mathrm{eq}^{-4}
</math>
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\Rightarrow ~~~~\frac{4\pi}{3} R_0^3 P_e</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl[ K M_\mathrm{tot} \biggl( \frac{3M_\mathrm{tot} }{4\pi R_0^3} \biggr)^{\gamma_g - 1}  \cdot \mathfrak{f}_A \biggr] \chi_\mathrm{eq}^{-3\gamma_g}
- \biggl[ \frac{1}{5} \frac{GM_\mathrm{tot} ^2}{R_0} \cdot \mathfrak{f}_W  \biggr] \chi_\mathrm{eq}^{-4}
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\Rightarrow ~~~~ P_e</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl[ K \biggl( \frac{3M_\mathrm{tot} }{4\pi R_\mathrm{eq}^3} \biggr)^{\gamma_g}  \cdot \mathfrak{f}_A \biggr]
- \biggl[ \biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{R_\mathrm{eq}^4} \cdot \mathfrak{f}_W  \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
P_c \cdot \mathfrak{f}_A - \biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{R_\mathrm{eq}^4} \cdot \mathfrak{f}_W  \, ,
</math>
  </td>
</tr>
</table>
</div>
where, in the last step as was [[User:Tohline/SSC/VirialEquilibrium/UniformDensity#Comparison_with_Detailed_Force-Balance_Model|recognized above]], we have set,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~K</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>P_c \biggl( \frac{4\pi R_\mathrm{eq}^3}{3M_\mathrm{tot}} \biggr)^{\gamma_g} \, .</math>
  </td>
</tr>
</table>
</div>
Hence, for any external pressure, <math>~P_e < P_c</math>, the pressure-confined equilibrium radius is,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>
~R_\mathrm{eq}
</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl[ \biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{P_c} \cdot \frac{\mathfrak{f}_W}{\mathfrak{f}_A}
\biggl( 1  -  \frac{P_e}{P_c} \cdot \frac{1}{\mathfrak{f}_A} \biggr)^{-1} \biggr]^{1/4} \, .
</math>
  </td>
</tr>
</table>
</div>
===Comparison with Detailed Force-Balance Model===
It is reasonable to ask how close this virial expression for the equilibrium radius is to the [[User:Tohline/SSC/Structure/UniformDensity#Uniform-Density_Sphere_Embedded_in_an_External_Medium|exact result]].  As before, from a free-energy analysis alone, the best we can do is assume that both structural form factors, <math>~\mathfrak{f}_W</math> and <math>\mathfrak{f}_A</math>, are of order unity.  But we can do better than this.  To begin with, because <math>~\rho</math> is uniform throughout the configuration, <math>~\mathfrak{f}_W = 1</math>, even though the configuration is truncated by the imposed external pressure.  We need to reassess how <math>~\mathfrak{f}_A</math> is evaluated, however, because the pressure does not drop to zero at the surface of the configuration.
Going back to our [[User:Tohline/SphericallySymmetricConfigurations/Virial#Energy_Content_for_a_System_of_a_Given_Size_and_Internal_Structure|original definition of the thermodynamic energy reservoir for spherically symmetric adiabatic systems]],
<div align="center">
<math>\mathfrak{W}_A = \frac{1}{({\gamma_g}-1)}  \int_0^R  4\pi r^2 P dr
\, ,</math>
</div>
we begin by normalizing the radial coordinate to <math>~R_0</math>, the radius of the isolated (''i.e.,'' not truncated) sphere, because we know [[User:Tohline/SSC/Structure/UniformDensity#Summary|from the detailed force-balanced solution]] that, structurally, the pressure varies with <math>~r</math> inside the configuration as,
<div align="center">
<math>\frac{P(x)}{P_c} = 1 - x^2 \, ,</math>
</div>
where, <math>~x \equiv r/R_0</math>.  Integrating only out to the edge of the ''truncated'' sphere, which we will identify as <math>~R_e</math> and, correspondingly,
<div align="center">
<math>~x_e \equiv \frac{R_e}{R_0} = \biggl( 1 - \frac{P_e}{P_c} \biggr)^{1/2} \, ,</math>
</div>
we have,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\mathfrak{W}_A</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>\frac{4\pi P_c R_0^3}{({\gamma_g}-1)}  \int_0^{x_e}  ( 1-x^2 ) x^2 dx</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\frac{4\pi P_c R_0^3}{({\gamma_g}-1)}  \biggl[ \frac{x^3}{3}-\frac{x^5}{5} \biggr]_0^{x_e}
= \frac{P_c }{({\gamma_g}-1)} \biggl( \frac{4\pi R_e^3}{3} \biggr)  \biggl[ 1-\frac{3}{5}x_e^2 \biggr]
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\frac{M_\mathrm{tot} }{({\gamma_g}-1)} \biggl( \frac{P_c}{\rho_c} \biggr)  \biggl[ 1-\frac{3}{5}\biggl(1 - \frac{P_e}{P_c}\biggr) \biggr] \, .
</math>
  </td>
</tr>
</table>
</div>
Hence, in the case of a pressure-truncated, uniform-density sphere, we surmise that the relevant structural form factor is,
<div align="center">
<math>
\mathfrak{f}_A = 1-\frac{3}{5}\biggl(1 - \frac{P_e}{P_c}\biggr) = \frac{2}{5} + \frac{3}{5}\frac{P_e}{P_c} \, .
</math>
</div>
Plugging this expression for <math>~\mathfrak{f}_A</math> along with <math>~\mathfrak{f}_W = 1</math> into the [[User:Tohline/SSC/VirialEquilibrium/UniformDensity#Comparison_with_Detailed_Force-Balance_Model_2|just-derived virial equilibrium solution]] gives,
<div align="center">
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~\biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{R_\mathrm{eq}^4} </math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~P_c \cdot \biggl[ \frac{2}{5} + \frac{3}{5}\frac{P_e}{P_c} \biggr] - P_e
</math>
  </td>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
~\frac{2}{5} P_c \biggl( 1 - \frac{P_e}{P_c} \biggr)
</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\Rightarrow ~~~~ R_\mathrm{eq}</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>
\biggl[ \biggl( \frac{3}{2^3\pi} \biggr) \frac{G M^2}{P_c} \biggl( 1 - \frac{P_e}{P_c} \biggr)^{-1} \biggr]^{1/4} \, .
</math>
  </td>
</tr>
</table>
</div>
</div>
This result exactly matches the solution for the equilibrium radius of a pressure-truncated, uniform-density sphere that has been [[User:Tohline/SSC/Structure/UniformDensity#Uniform-Density_Sphere_Embedded_in_an_External_Medium|derived elsewhere]].
=See Also=
<ul>
<li>[[User:Tohline/SphericallySymmetricConfigurations/IndexFreeEnergy#Index_to_Free-Energy_Analyses|Index to a Variety of Free-Energy and/or Virial Analyses]]</li>
</ul>




{{LSU_HBook_footer}}
{{LSU_HBook_footer}}

Latest revision as of 23:15, 10 July 2016

Uniform-Density Sphere

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

Review

In an introductory discussion of the virial equilibrium structure of spherically symmetric configurations — see especially the section titled, Energy Extrema — we deduced that a system's equilibrium radius, <math>~R_\mathrm{eq}</math>, measured relative to a reference length scale, <math>~R_0</math>, i.e., the dimensionless equilibrium radius,

<math>~\chi_\mathrm{eq} \equiv \frac{R_\mathrm{eq}}{R_0} \, ,</math>

is given by the root(s) of the following equation:

<math> 2C \chi^{-2} + ~ (1-\delta_{1\gamma_g})~3(\gamma_g-1) B\chi^{3 -3\gamma_g} +~ \delta_{1\gamma_g} B_I ~-~A\chi^{-1} -~ 3D\chi^3 = 0 \, , </math>

where the definitions of the various coefficients are,

<math>~A</math>

<math>~\equiv</math>

<math>\frac{3}{5} \frac{GM_\mathrm{tot} ^2}{R_0} \cdot \mathfrak{f}_W \, ,</math>

<math>~B</math>

<math>~\equiv</math>

<math> \frac{K M_\mathrm{tot} }{(\gamma_g-1)} \biggl( \frac{3M_\mathrm{tot} }{4\pi R_0^3} \biggr)^{\gamma_g - 1} \cdot \mathfrak{f}_A = \frac{\bar{c_s}^2 M_\mathrm{tot} }{(\gamma_g - 1)} \cdot \mathfrak{f}_A \, , </math>

<math>~B_I</math>

<math>~\equiv</math>

<math> 3c_s^2 M_\mathrm{tot} \cdot \mathfrak{f}_M \, , </math>

<math>~C</math>

<math>~\equiv</math>

<math> \frac{5J^2}{4M_\mathrm{tot} R_0^2} \cdot \mathfrak{f}_T \, , </math>

<math>~D</math>

<math>~\equiv</math>

<math> \frac{4}{3} \pi R_0^3 P_e \, . </math>

Once the pressure exerted by the external medium (<math>~P_e</math>), and the configuration's mass (<math>~M_\mathrm{tot}</math>), angular momentum (<math>~J</math>), and specific entropy (via <math>~K</math>) — or, in the isothermal case, sound speed (<math>~c_s</math>) — have been specified, the values of all of the coefficients are known and <math>~\chi_\mathrm{eq}</math> can be determined.

Adiabatic Evolution of an Isolated Sphere

Here we seek to determine the equilibrium radius of a non-rotating configuration (<math>~J = 0</math>) that undergoes adiabatic compression/expansion (<math>\delta_{1\gamma_g} =~0</math>) and that is not confined by an external medium (<math>P_e = 0~</math>).

Solution

In this case, the statement of virial equilibrium is simplified considerably. Specifically, <math>~\chi_\mathrm{eq}</math> is given by the root(s) of the equation,

<math>~A\chi_\mathrm{eq}^{-1} </math>

<math>~=</math>

<math>~3(\gamma_g-1) B\chi_\mathrm{eq}^{3 -3\gamma_g} </math>

<math>\Rightarrow ~~~~~\chi_\mathrm{eq}^{3\gamma_g-4} </math>

<math>~=</math>

<math>~\frac{3(\gamma_g-1) B}{A} </math>

 

<math>~=</math>

<math> \biggl[ 3K M_\mathrm{tot} \biggl( \frac{3M_\mathrm{tot} }{4\pi R_0^3} \biggr)^{\gamma_g - 1} \cdot \mathfrak{f}_A \biggr] \biggl[ \frac{3}{5} \frac{GM_\mathrm{tot} ^2}{R_0} \cdot \mathfrak{f}_W \biggr] ^{-1} </math>

 

<math>~=</math>

<math> \biggl( \frac{1}{R_0} \biggr)^{3\gamma_g-4} \biggl[ 5\biggl( \frac{3}{4\pi} \biggr)^{\gamma_g-1} \biggr(\frac{K}{G}\biggr) M^{(\gamma_g-2)} \cdot \frac{\mathfrak{f}_A}{\mathfrak{f}_W} \biggr] \, . </math>

In other words,

<math> R_\mathrm{eq} = \biggl[ 5\biggl( \frac{3}{4\pi} \biggr)^{\gamma_g-1} \biggr(\frac{K}{G}\biggr) M^{(\gamma_g-2)} \cdot \frac{\mathfrak{f}_A}{\mathfrak{f}_W} \biggr]^{1/(3\gamma_g-4)} \, . </math>

Comparison with Detailed Force-Balance Model

This derived solution will look more familiar if, instead of <math>~K</math>, we express the solution in terms of the central pressure,

<math>P_c = K\rho_0^{\gamma_g} \, ,</math>

where, for this uniform-density sphere, <math>~\rho_0 = 3M_\mathrm{tot}/(4\pi R_\mathrm{eq}^3)</math>. Hence,

<math>~K</math>

<math>~=</math>

<math>P_c \biggl( \frac{4\pi R_\mathrm{eq}^3}{3M_\mathrm{tot}} \biggr)^{\gamma_g} \, ,</math>

and the solution takes the form,

<math>R_\mathrm{eq}^{3\gamma_g - 4}</math>

<math>~=</math>

<math>5\biggl( \frac{4\pi}{3} \biggr) \biggr(\frac{P_c R_\mathrm{eq}^{3\gamma_g}}{GM^2_\mathrm{tot}}\biggr) \cdot \frac{\mathfrak{f}_A}{\mathfrak{f}_W} </math>

<math>\Rightarrow ~~~~ R_\mathrm{eq}^{4} </math>

<math>~=</math>

<math>\biggl( \frac{3}{20\pi} \biggr) \biggr(\frac{GM^2_\mathrm{tot}}{P_c}\biggr) \cdot \frac{\mathfrak{f}_W}{\mathfrak{f}_A} \, .</math>

Or, solving for the central pressure,

<math>~P_c</math>

<math>~=</math>

<math>\biggl( \frac{3}{20\pi} \cdot \frac{\mathfrak{f}_W}{\mathfrak{f}_A} \biggr) \frac{GM^2_\mathrm{tot}}{R_\mathrm{eq}^{4} } \, .</math>

This should be compared with our detailed force-balance solution of the interior structure of an isolated, nonrotating, uniform-density sphere, which gives the precise expression,

<math>~P_c</math>

<math>~=</math>

<math>\biggl( \frac{3}{8\pi} \biggr) \frac{GM^2_\mathrm{tot}}{R_\mathrm{eq}^{4} } \, .</math>

The expression for <math>~P_c</math> derived from our identification of an extremum in the free energy is identical to the expression derived from the more precise, detailed force-balance analysis, except that the leading numerical coefficients differ by a factor of <math>~(5\mathfrak{f}_A/2\mathfrak{f}_W)</math>.

From a free-energy analysis alone, the best we can do is assume that both structural form factors, <math>~\mathfrak{f}_W</math> and <math>\mathfrak{f}_A</math>, are of order unity. But knowing the detailed force-balance solution allows us to evaluate both form factors. From our introductory discussion of the free energy function, their respective definitions are,

<math>~\mathfrak{f}_W</math>

<math>~\equiv</math>

<math>~ 3\cdot 5 \int_0^1 \biggl\{ \int_0^x \biggl[ \frac{\rho(x)}{\rho_c}\biggr] x^2 dx \biggr\} \biggl[ \frac{\rho(x)}{\rho_c}\biggr] x dx\, ,</math>

<math>~\mathfrak{f}_A</math>

<math>~\equiv</math>

<math>~ 3\int_0^1 \biggl[ \frac{P(x)}{P_c}\biggr] x^2 dx \, .</math>

Now, because the configuration under discussion has a uniform density, we should set <math>~\rho(x)/\rho_c = 1</math> in the definition of <math>~\mathfrak{f}_W</math> which, after evaluation of the nested integrals, gives <math>~\mathfrak{f}_W = 1</math>. But, instead of being uniform throughout the configuration, in the detailed force-balance model, the pressure drops from the center to the surface according to the relation,

<math>\frac{P(x)}{P_c} = 1 - x^2 \, .</math>

Integrating over this function, in accordance with the definition of <math>~\mathfrak{f}_A</math>, gives,

<math>~\mathfrak{f}_A</math>

<math>~\equiv</math>

<math>~ 3\int_0^1 (1-x^2) x^2 dx = 3 \biggl[ \frac{x^3}{3} - \frac{x^5}{5} \biggr]_0^1 = \frac{2}{5} \, .</math>

Hence, the ratio,

<math> \frac{5\mathfrak{f}_A}{2\mathfrak{f}_W} = 1 \, , </math>

which brings into perfect agreement the two separate determinations of the equilibrium expressions for <math>~R_\mathrm{eq}</math> and <math>~P_c</math> in terms of one another and the total mass.

This demonstrates that the free-energy approach to determining the equilibrium radius of a spherical configuration is only handicapped by its inability to precisely nail down values of the structural form factors. But this is not a severe limitation as the (dimensionless) form factors are generally of order unity. In contrast, the free-energy analysis brings with it a capability to readily evaluate the global stability of equilibrium configurations.


Adiabatic Evolution of Pressure-truncated Sphere

Here we seek to determine the equilibrium radius of a non-rotating configuration (<math>~J = 0</math>) that undergoes adiabatic compression/expansion (<math>\delta_{1\gamma_g} =~0</math>) and that is embedded in a hot, tenuous external medium whose confining pressure, <math>~P_e</math>, truncates the configuration.

Solution

In this case, virial equilibrium implies that <math>~\chi_\mathrm{eq}</math> is given by the root(s) of the equation,

<math> 3(\gamma_g-1) B\chi^{3 -3\gamma_g} ~ -~A\chi^{-1} -~ 3D\chi^3 = 0 \, . </math>

Hence,

<math>~D</math>

<math>~=</math>

<math> (\gamma_g-1) B\chi_\mathrm{eq}^{-3\gamma_g} - \frac{A}{3}\chi_\mathrm{eq}^{-4} </math>

<math>\Rightarrow ~~~~\frac{4\pi}{3} R_0^3 P_e</math>

<math>~=</math>

<math> \biggl[ K M_\mathrm{tot} \biggl( \frac{3M_\mathrm{tot} }{4\pi R_0^3} \biggr)^{\gamma_g - 1} \cdot \mathfrak{f}_A \biggr] \chi_\mathrm{eq}^{-3\gamma_g} - \biggl[ \frac{1}{5} \frac{GM_\mathrm{tot} ^2}{R_0} \cdot \mathfrak{f}_W \biggr] \chi_\mathrm{eq}^{-4} </math>

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

<math>~=</math>

<math> \biggl[ K \biggl( \frac{3M_\mathrm{tot} }{4\pi R_\mathrm{eq}^3} \biggr)^{\gamma_g} \cdot \mathfrak{f}_A \biggr] - \biggl[ \biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{R_\mathrm{eq}^4} \cdot \mathfrak{f}_W \biggr] </math>

 

<math>~=</math>

<math> P_c \cdot \mathfrak{f}_A - \biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{R_\mathrm{eq}^4} \cdot \mathfrak{f}_W \, , </math>

where, in the last step as was recognized above, we have set,

<math>~K</math>

<math>~=</math>

<math>P_c \biggl( \frac{4\pi R_\mathrm{eq}^3}{3M_\mathrm{tot}} \biggr)^{\gamma_g} \, .</math>

Hence, for any external pressure, <math>~P_e < P_c</math>, the pressure-confined equilibrium radius is,

<math>

~R_\mathrm{eq}

</math>

<math>~=</math>

<math> \biggl[ \biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{P_c} \cdot \frac{\mathfrak{f}_W}{\mathfrak{f}_A} \biggl( 1 - \frac{P_e}{P_c} \cdot \frac{1}{\mathfrak{f}_A} \biggr)^{-1} \biggr]^{1/4} \, . </math>

Comparison with Detailed Force-Balance Model

It is reasonable to ask how close this virial expression for the equilibrium radius is to the exact result. As before, from a free-energy analysis alone, the best we can do is assume that both structural form factors, <math>~\mathfrak{f}_W</math> and <math>\mathfrak{f}_A</math>, are of order unity. But we can do better than this. To begin with, because <math>~\rho</math> is uniform throughout the configuration, <math>~\mathfrak{f}_W = 1</math>, even though the configuration is truncated by the imposed external pressure. We need to reassess how <math>~\mathfrak{f}_A</math> is evaluated, however, because the pressure does not drop to zero at the surface of the configuration.

Going back to our original definition of the thermodynamic energy reservoir for spherically symmetric adiabatic systems,

<math>\mathfrak{W}_A = \frac{1}{({\gamma_g}-1)} \int_0^R 4\pi r^2 P dr

\, ,</math>

we begin by normalizing the radial coordinate to <math>~R_0</math>, the radius of the isolated (i.e., not truncated) sphere, because we know from the detailed force-balanced solution that, structurally, the pressure varies with <math>~r</math> inside the configuration as,

<math>\frac{P(x)}{P_c} = 1 - x^2 \, ,</math>

where, <math>~x \equiv r/R_0</math>. Integrating only out to the edge of the truncated sphere, which we will identify as <math>~R_e</math> and, correspondingly,

<math>~x_e \equiv \frac{R_e}{R_0} = \biggl( 1 - \frac{P_e}{P_c} \biggr)^{1/2} \, ,</math>

we have,

<math>~\mathfrak{W}_A</math>

<math>~=</math>

<math>\frac{4\pi P_c R_0^3}{({\gamma_g}-1)} \int_0^{x_e} ( 1-x^2 ) x^2 dx</math>

 

<math>~=</math>

<math> \frac{4\pi P_c R_0^3}{({\gamma_g}-1)} \biggl[ \frac{x^3}{3}-\frac{x^5}{5} \biggr]_0^{x_e} = \frac{P_c }{({\gamma_g}-1)} \biggl( \frac{4\pi R_e^3}{3} \biggr) \biggl[ 1-\frac{3}{5}x_e^2 \biggr] </math>

 

<math>~=</math>

<math> \frac{M_\mathrm{tot} }{({\gamma_g}-1)} \biggl( \frac{P_c}{\rho_c} \biggr) \biggl[ 1-\frac{3}{5}\biggl(1 - \frac{P_e}{P_c}\biggr) \biggr] \, . </math>

Hence, in the case of a pressure-truncated, uniform-density sphere, we surmise that the relevant structural form factor is,

<math> \mathfrak{f}_A = 1-\frac{3}{5}\biggl(1 - \frac{P_e}{P_c}\biggr) = \frac{2}{5} + \frac{3}{5}\frac{P_e}{P_c} \, . </math>

Plugging this expression for <math>~\mathfrak{f}_A</math> along with <math>~\mathfrak{f}_W = 1</math> into the just-derived virial equilibrium solution gives,

<math>~\biggl(\frac{3}{20\pi} \biggr) \frac{GM_\mathrm{tot} ^2}{R_\mathrm{eq}^4} </math>

<math>~=</math>

<math> ~P_c \cdot \biggl[ \frac{2}{5} + \frac{3}{5}\frac{P_e}{P_c} \biggr] - P_e </math>

 

<math>~=</math>

<math> ~\frac{2}{5} P_c \biggl( 1 - \frac{P_e}{P_c} \biggr) </math>

<math>\Rightarrow ~~~~ R_\mathrm{eq}</math>

<math>~=</math>

<math> \biggl[ \biggl( \frac{3}{2^3\pi} \biggr) \frac{G M^2}{P_c} \biggl( 1 - \frac{P_e}{P_c} \biggr)^{-1} \biggr]^{1/4} \, . </math>

This result exactly matches the solution for the equilibrium radius of a pressure-truncated, uniform-density sphere that has been derived elsewhere.


See Also


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