Difference between revisions of "User:Tohline/SSC/Structure/BiPolytropes/Analytic0 0"

From VistrailsWiki
Jump to navigation Jump to search
(→‎Thermal Energy Content: Finished integrations to obtain both S_env and S_core)
(→‎Thermal Energy Content: Simplified expressions and introduced capital pi)
Line 507: Line 507:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>6\pi R^3 \int_0^q  \biggl\{ P_i + \biggl( \frac{\lambda}{q^2} \biggr) (q^2 -  x^2) \biggr\} x^2 dx </math>
<math>6\pi R^3 \int_0^q  \biggl\{ P_i + \Pi (q^2 -  x^2) \biggr\} x^2 dx </math>
   </td>
   </td>
</tr>
</tr>
Line 519: Line 519:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>6\pi R^3 \biggl\{ \biggl[ \biggl(P_i + \lambda\biggr)\frac{x^3}{3} \biggr]_0^q  
<math>6\pi R^3 \biggl\{ \biggl[ \biggl(P_i + q^2 \Pi \biggr)\frac{x^3}{3} \biggr]_0^q  
- \biggl[\biggl( \frac{\lambda}{q^2} \biggr) \frac{x^5}{5} \biggr]_0^q \biggr\}  </math>
- \biggl[ \Pi \frac{x^5}{5} \biggr]_0^q \biggr\}  </math>
   </td>
   </td>
</tr>
</tr>
Line 532: Line 532:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>\biggl( \frac{2\pi}{5} \biggr) R^3 q^3 (5P_i + 2\lambda)  \, ,</math>
<math>\biggl( \frac{4\pi}{5} \biggr) R^3 q^5 \biggl (\frac{5P_i}{2q^2} + \Pi \biggr)  \, ,</math>
   </td>
   </td>
</tr>
</tr>
Line 544: Line 544:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math> 6\pi R^3 \int_q^1 \biggl\{ P_i - \biggl( \frac{\lambda}{q^2} \biggr) \frac{\rho_e}{\rho_0} \biggl( a_e + \frac{b_e}{x} + c_e x^2\biggr) \biggr\} x^2 dx </math>
<math> 6\pi R^3 \int_q^1 \biggl\{ P_i - \Pi \frac{\rho_e}{\rho_0} \biggl( a_e + \frac{b_e}{x} + c_e x^2\biggr) \biggr\} x^2 dx </math>
   </td>
   </td>
</tr>
</tr>
Line 557: Line 557:
   <td align="left">
   <td align="left">
<math> 6\pi R^3 \biggl\{ \biggl[ P_i  
<math> 6\pi R^3 \biggl\{ \biggl[ P_i  
- a_e \biggl( \frac{\lambda}{q^2} \biggr) \biggl( \frac{\rho_e}{\rho_0} \biggr)\biggr] \frac{x^3}{3} \biggr|_q^1  
- a_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr)\biggr] \frac{x^3}{3} \biggr|_q^1  
- \biggl[ b_e \biggl( \frac{\lambda}{q^2} \biggr) \biggl( \frac{\rho_e}{\rho_0} \biggr) \frac{x^2}{2} \biggr]_q^1
- \biggl[ b_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \frac{x^2}{2} \biggr]_q^1
- \biggl[ c_e \biggl( \frac{\lambda}{q^2} \biggr) \biggl( \frac{\rho_e}{\rho_0} \biggr) \frac{x^5}{5} \biggr]_q^1 \biggr\}  </math>
- \biggl[ c_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \frac{x^5}{5} \biggr]_q^1 \biggr\}  </math>
   </td>
   </td>
</tr>
</tr>
Line 572: Line 572:
   <td align="left">
   <td align="left">
<math> \frac{\pi R^3}{5} \biggl\{ 10 \biggl[ P_i  
<math> \frac{\pi R^3}{5} \biggl\{ 10 \biggl[ P_i  
- a_e \biggl( \frac{\lambda}{q^2} \biggr) \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggr] (1-q^3)  
- a_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggr] (1-q^3)  
- 15 \biggl[ b_e \biggl( \frac{\lambda}{q^2} \biggr) \biggl( \frac{\rho_e}{\rho_0} \biggr) (1-q^2) \biggr]
- 15 \biggl[ b_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) (1-q^2) \biggr]
- 6\biggl[ c_e \biggl( \frac{\lambda}{q^2} \biggr) \biggl( \frac{\rho_e}{\rho_0} \biggr) (1-q^5) \biggr] \biggr\}  </math>
- 6\biggl[ c_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) (1-q^5) \biggr] \biggr\}  </math>
   </td>
   </td>
</tr>
</tr>


<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=~</math>
  </td>
  <td align="left">
<math> \frac{\pi R^3}{5} \biggl\{ 10 P_i  (1-q^3) - \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggl[ 10 a_e (1-q^3)
+ 15 b_e (1-q^2) +  6 c_e  (1-q^5) \biggr] \biggr\}  </math>
  </td>
</tr>
</table>
</table>
</div>
</div>
Line 583: Line 595:
<div align="center">
<div align="center">
<math>
<math>
\lambda \equiv \frac{3}{2^3 \pi}  
\Pi \equiv \frac{3}{2^3 \pi}  
\biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^4} \biggr) \, ,
\biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) \, ,
</math>
</math>
</div>
</div>
Line 598: Line 610:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>\biggl(2 - 3\frac{\rho_e}{\rho_0} \biggr) q^2 \, ,</math>
<math>3 \biggl(\frac{2}{3} - \frac{\rho_e}{\rho_0} \biggr) q^2 \, ,</math>
   </td>
   </td>
</tr>
</tr>

Revision as of 00:03, 11 February 2014

BiPolytrope with <math>n_c = 0</math> and <math>n_e=0</math>

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

Here we construct a bipolytrope in which both the core and the envelope have uniform densities, that is, the structure of both the core and the envelope will be modeled using an <math>n = 0</math> polytropic index. It should be possible for the entire structure to be described by closed-form, analytic expressions. Generally, we will follow the general solution steps for constructing a bipolytrope that we have outlined elsewhere. [On 1 February 2014, J. E. Tohline wrote: This particular system became of interest to me during discussions with Kundan Kadam about the relative stability of bipolytropes.]

Step 4: Throughout the core (<math>0 \le \chi \le \chi_i</math>)

Specify: <math>~P_0</math> and <math>\rho_0 ~\Rightarrow</math>

 

<math>~\rho</math>

  <math>~=</math> 

<math>~\rho_0</math>

 

 

<math>~P</math>

  <math>~=</math> 

<math>P_0 - \frac{2}{3} \pi G \rho_0^2 r^2</math>

  <math>~=</math> 

<math>P_0 \biggl( 1 - \frac{2\pi}{3}\chi^2 \biggr)</math>

<math>~r</math>

  <math>~=</math> 

<math>\biggl[ \frac{P_0}{G \rho_0^2} \biggr]^{1/2} \chi</math>

  <math>~=</math> 

<math>\biggl[ \frac{P_0}{G \rho_0^2} \biggr]^{1/2} \chi</math>

<math>~M_r</math>

  <math>~=</math> 

<math>\frac{4\pi}{3} \rho_0 r^3</math>

  <math>~=</math> 

<math>\frac{4\pi}{3} \rho_0 \biggl[ \frac{P_0}{G \rho_0^2} \biggr]^{3/2} \chi^3 = \frac{4\pi}{3} \biggl[ \frac{P_0^3}{G^3 \rho_0^4} \biggr]^{1/2} \chi^3</math>

Step 5: Interface Conditions

Specify: <math>~\chi_i</math> and <math>~\rho_e/\rho_0</math>, and demand …

 

<math>~P_{ei}</math>

  <math>~=</math> 

<math>~P_{ci}</math>

  <math>~=</math> 

<math>P_0 \biggl( 1 - \frac{2\pi}{3}\chi_i^2 \biggr)</math>

Step 6: Envelope Solution (<math>~\chi > \chi_i</math>)

<math>~\rho</math>

  <math>~=</math> 

<math>~\rho_e</math>

<math>~P</math>

  <math>~=</math> 

<math>P_{ei} + \biggl(\frac{2}{3} \pi G \rho_e\biggr) \biggl[ 2(\rho_0 - \rho_e) r_i^3\biggl( \frac{1}{r} - \frac{1}{r_i}\biggr) - \rho_e(r^2 - r_i^2) \biggr]</math>

 

  <math>~=</math> 

<math>P_{ei} + \frac{2\pi}{3} \biggl(\frac{\rho_e}{\rho_0}\biggr) P_0 \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \chi_i^3\biggl( \frac{1}{\chi} - \frac{1}{\chi_i}\biggr) - \frac{\rho_e}{\rho_0} (\chi^2 - \chi_i^2) \biggr]</math>

<math>~\frac{P}{P_0}</math>

  <math>~=</math> 

<math>1 - \frac{2\pi}{3}\chi_i^2 + \frac{2\pi}{3} \biggl(\frac{\rho_e}{\rho_0}\biggr) \chi_i^2 \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( \frac{1}{\xi} - 1\biggr) - \frac{\rho_e}{\rho_0} (\xi^2 - 1) \biggr]</math>

<math>~M_r</math>

  <math>~=</math> 

<math>\frac{4\pi}{3} \biggl[ \rho_0 r_i^3 + \rho_e(r^3 - r_i^3) \biggr]</math>

 

  <math>~=</math> 

<math>\frac{4\pi}{3} \biggl[ \frac{P_0^3}{G^3 \rho_0^4} \biggr]^{1/2} \biggl[\chi_i^3 +\frac{\rho_e}{\rho_0} \biggl( \chi^3 - \chi_i^3 \biggr) \biggr]</math>

 

  <math>~=</math> 

<math>\frac{4\pi}{3} \biggl[ \frac{P_0^3}{G^3 \rho_0^4} \biggr]^{1/2} \chi_i^3\biggl[1 +\frac{\rho_e}{\rho_0} \biggl( \xi^3 - 1\biggr) \biggr]</math>


Step 7: Surface Boundary Condition

At the surface (that is, at <math>r = R</math> and <math>M_r = M_\mathrm{tot}</math>), <math>P/P_0 = 0</math> and <math>\xi = \xi_s = R/r_i = 1/q</math>. Also, we can write,

<math> \chi_i = q\biggl[ \frac{G\rho_0^2 R^2}{P_0} \biggr]^{1/2} \, ; </math>

and, from earlier derivations,

<math> R^3 = \frac{3M_\mathrm{tot}}{4\pi \bar\rho} = \frac{3M_\mathrm{tot}}{4\pi \rho_0} \biggl( \frac{\nu}{q^3} \biggr) \, ; </math>

<math> \frac{\rho_e}{\rho_0} = \frac{q^3}{\nu} \biggl( \frac{1-\nu}{1-q^3}\biggr) \, . </math>

Therefore, setting the pressure to zero at the surface means,

<math>~\frac{3}{2\pi}\chi_i^{-2}</math>

  <math>~=</math> 

<math>1 - \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( q - 1\biggr) - \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr]</math>

<math>\Rightarrow ~~~~~\biggl( \frac{3}{2\pi} \biggr) q^{-2} \biggl[ \frac{P_0}{G\rho_0^2 R^2} \biggr]</math>

  <math>~=</math> 

<math>1 + \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1- q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] </math>

<math>\Rightarrow ~~~~~\biggl( \frac{3}{2\pi} \biggr) \biggl( \frac{4\pi}{3} \biggr)^{2/3} \nu^{-2/3} \biggl[ \frac{P_0^3}{G^3\rho_0^4 M_\mathrm{tot}^2} \biggr]^{1/3}</math>

  <math>~=</math> 

<math>1 + \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1- q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] </math>

<math>\Rightarrow ~~~~~\biggl\{ \biggl( \frac{6}{\pi} \biggr) \frac{1}{\nu^{2} } \biggl[ \frac{P_0^3}{G^3\rho_0^4 M_\mathrm{tot}^2} \biggr] \biggr\}^{1/3}</math>

  <math>~=</math> 

<math>1 + \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1- q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] </math>

It therefore seems prudent to define a function,

<math> g(\nu,q) \equiv \biggl\{ 1 + \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl[ 2 \biggl(1 - \frac{\rho_e}{\rho_0} \biggr) \biggl( 1-q \biggr) + \frac{\rho_e}{\rho_0} \biggl(\frac{1}{q^2} - 1\biggr) \biggr] \biggr\}^{1/2}\, , </math>

in which case the expressions for the equilibrium radius and equilibrium total mass are, respectively,

<math> \biggl[ \frac{G\rho_0^2}{P_0} \biggr]^{1/2} R</math>

  <math>~=</math> 

<math>\biggl( \frac{3}{2\pi} \biggr)^{1/2} \frac{1}{q g} \, ;</math>

<math> \biggl[ \frac{G^3\rho_0^4}{P_0^3} \biggr]^{1/2} M_\mathrm{tot}</math>

  <math>~=</math> 

<math>\biggl( \frac{6}{\pi} \biggr)^{1/2} \frac{1}{\nu g^3} \, .</math>

We can also combine these two expressions and eliminate direct reference to the central density, <math>\rho_0</math>, obtaining,

<math> \biggl[ \frac{R^4}{GM_\mathrm{tot}^2} \biggr] P_0</math>

  <math>~=</math> 

<math>\biggl( \frac{3}{2^3\pi} \biggr) \frac{\nu^2 g^2}{q^4} \, .</math>

Thermal Energy Content

Knowing the pressure distribution throughout this bipolytropic configuration allows us to obtain an analytic expression for the configuration's total thermal content. Specifically, the differential contribution to the total thermal energy that is made by each spherical shell is,

<math>dS = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) dm = \frac{3}{2} \biggl( \frac{P}{\rho} \biggr) 4\pi \rho r^2 dr = 6\pi R^3 P(x) x^2 dx \, ,</math>

where,

<math>x \equiv \frac{r}{R} \, .</math>

(We are switching to a new normalization of the radial coordinate — using <math>~x</math> in preference to <math>~\chi</math> — because, at the interface between the core and the envelope, <math>~x=q</math>, while <math>~x=1 </math> at the surface of the bipolytropic configuration.) Hence, the thermal content of the core and of the envelope will be given by performing the following integrals, respectively:

<math>~S_\mathrm{core}</math>

<math>~=~</math>

<math>6\pi R^3 \int_0^q P_\mathrm{core}(x) x^2 dx \, ;</math>

<math>~S_\mathrm{env}</math>

<math>~=~</math>

<math> 6\pi R^3 \int_q^1 P_\mathrm{env}(x) x^2 dx \, .</math>

Drawing from the expressions obtained in step #4 and step #6, above, the relevant functions <math>~P(x)</math> are,

<math>~P_\mathrm{core}</math>

  <math>~=</math> 

<math>P_0 - \frac{2}{3} \pi G \rho_0^2 r^2 = P_i + \frac{2}{3} \pi G \rho_0^2(r_i^2 - r^2)</math>

 

  <math>~=</math> 

<math>P_i + \frac{2}{3} \pi G \rho_0^2 R^2 (q^2 - x^2) = P_i + \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) (q^2 - x^2) \, ;</math>

<math>~P_\mathrm{env}</math>

  <math>~=</math> 

<math>P_{i} - \biggl(\frac{2}{3} \pi G \rho_e\biggr) \biggl[ 2(\rho_0 - \rho_e) r_i^3\biggl( \frac{1}{r_i} - \frac{1}{r} \biggr) + \rho_e(r^2 - r_i^2) \biggr]</math>

 

  <math>~=</math> 

<math>P_{i} - \biggl(\frac{2}{3} \pi G \rho_0^2 R^2\biggr) \frac{\rho_e}{\rho_0} \biggl[ 2\biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\biggl( \frac{1}{q} - \frac{1}{x} \biggr) + \frac{\rho_e}{\rho_0} (x^2 - q^2) \biggr]</math>

 

  <math>~=</math> 

<math>P_{i} - \biggl(\frac{2}{3} \pi G \rho_0^2 R^2\biggr) \frac{\rho_e}{\rho_0} \biggl[ 2\biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\biggl( \frac{1}{q} - \frac{1}{x} \biggr) + \frac{\rho_e}{\rho_0} (x^2 - q^2) \biggr]</math>

 

  <math>~=</math> 

<math>P_{i} - \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) \frac{\rho_e}{\rho_0} \biggl[ 2\biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\biggl( \frac{1}{q} - \frac{1}{x} \biggr) + \frac{\rho_e}{\rho_0} (x^2 - q^2) \biggr] \, .</math>


Reminder: For a given choice of the parameter set <math>~(M_\mathrm{tot}, R, \nu, q)</math>, — and remembering that choosing the parameter pair <math>~(\nu,q)</math> sets the density ratio via the expression, <math>~\rho_e/\rho_0 = [q^3(1-\nu)]/[\nu(1-q^3)]</math> — the value of the pressure at the interface, <math>~P_i</math>, is determined by setting boundary conditions at the surface of the configuration, namely, by setting <math>~x=1</math> and <math>~P_\mathrm{env} = 0</math> in the second expression. Then, the central pressure <math>~(P_\mathrm{core} = P_0)</math> is determined by setting <math>~x=0</math> and inserting the determined value for <math>~P_i</math> in the first expression.

Hence,

<math>~S_\mathrm{core}</math>

<math>~=~</math>

<math>6\pi R^3 \int_0^q \biggl\{ P_i + \Pi (q^2 - x^2) \biggr\} x^2 dx </math>

 

<math>~=~</math>

<math>6\pi R^3 \biggl\{ \biggl[ \biggl(P_i + q^2 \Pi \biggr)\frac{x^3}{3} \biggr]_0^q - \biggl[ \Pi \frac{x^5}{5} \biggr]_0^q \biggr\} </math>

 

<math>~=~</math>

<math>\biggl( \frac{4\pi}{5} \biggr) R^3 q^5 \biggl (\frac{5P_i}{2q^2} + \Pi \biggr) \, ,</math>

<math>~S_\mathrm{env}</math>

<math>~=~</math>

<math> 6\pi R^3 \int_q^1 \biggl\{ P_i - \Pi \frac{\rho_e}{\rho_0} \biggl( a_e + \frac{b_e}{x} + c_e x^2\biggr) \biggr\} x^2 dx </math>

 

<math>~=~</math>

<math> 6\pi R^3 \biggl\{ \biggl[ P_i - a_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr)\biggr] \frac{x^3}{3} \biggr|_q^1 - \biggl[ b_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \frac{x^2}{2} \biggr]_q^1 - \biggl[ c_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \frac{x^5}{5} \biggr]_q^1 \biggr\} </math>

 

<math>~=~</math>

<math> \frac{\pi R^3}{5} \biggl\{ 10 \biggl[ P_i - a_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggr] (1-q^3) - 15 \biggl[ b_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) (1-q^2) \biggr] - 6\biggl[ c_e \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) (1-q^5) \biggr] \biggr\} </math>

 

<math>~=~</math>

<math> \frac{\pi R^3}{5} \biggl\{ 10 P_i (1-q^3) - \Pi \biggl( \frac{\rho_e}{\rho_0} \biggr) \biggl[ 10 a_e (1-q^3) + 15 b_e (1-q^2) + 6 c_e (1-q^5) \biggr] \biggr\} </math>

where,

<math> \Pi \equiv \frac{3}{2^3 \pi} \biggl( \frac{GM_\mathrm{tot}^2}{R^4} \biggr) \biggl( \frac{\nu^2}{q^6} \biggr) \, , </math>

and, in the expression for <math>~S_\mathrm{env}</math>, we have used the shorthand notation,

<math>~a_e</math>

<math>~\equiv~</math>

<math>3 \biggl(\frac{2}{3} - \frac{\rho_e}{\rho_0} \biggr) q^2 \, ,</math>

<math>~b_e</math>

<math>~\equiv~</math>

<math>- 2\biggl(1 - \frac{\rho_e}{\rho_0} \biggr) q^3\ \, ,</math>

<math>~c_e</math>

<math>~\equiv~</math>

<math>\frac{\rho_e}{\rho_0} \, .</math>

Related Discussions

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