SSC/Structure/BiPolytropes/51RenormaizePart3: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Created page with "__FORCETOC__ <!-- __NOTOC__ will force TOC off --> =BiPolytrope with n<sub>c</sub> = 5 and n<sub>e</sub> = 1= {| class="PGEclass" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black" |- ! style="height: 125px; width: 125px; background-color:white;" |<font size="-1"><b>Eggleton, Faulkner<br />& Cannon (1998)<br /><br />Analytic</b><br />(n<sub>c</sub>, n<sub>e</sub>) = (5, 1)</font> |} F..."
 
Line 3: Line 3:


=BiPolytrope with n<sub>c</sub> = 5 and n<sub>e</sub> = 1=
=BiPolytrope with n<sub>c</sub> = 5 and n<sub>e</sub> = 1=
{| class="PGEclass" style="float:left; margin-right: 20px; border-style: solid; border-width: 3px border-color: black"
|-
! style="height: 125px; width: 125px; background-color:white;" |<font size="-1">[[H_BookTiledMenu#MoreModels|<b>Eggleton, Faulkner<br />&amp; Cannon (1998)<br /><br />Analytic</b>]]<br />(n<sub>c</sub>, n<sub>e</sub>) = (5, 1)</font>
|}
[[File:CommentButton02.png|right|100px|Comment by J. E. Tohline on 30 March 2013:  As far as I have been able to determine, this analytic structural model has not previously been published in a refereed, archival journal.      Subsequent comment by J. E. Tohline on 23 June 2013:  Last night I stumbled upon an article by Eagleton, Faulkner, and Cannon (1998) in which this identical analytically definable bipolytrope has been presented.  Insight drawn from this article is presented in an additional subsection, below.]]
Here we construct a [[SSC/Structure/BiPolytropes#BiPolytropes|bipolytrope]] in which the core has an <math>~n_c=5</math> polytropic index and the envelope has an <math>~n_e=1</math> polytropic index.  This system is particularly interesting because the entire structure can be described by closed-form, analytic expressions.  In deriving the properties of this model, we will follow the [[SSC/Structure/BiPolytropes#Solution_Steps|general solution steps for constructing a bipolytrope]] that we have outlined elsewhere. 
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />
&nbsp;<br />


==Steps 2 &amp; 3==
After studying {{ Yabushita75full }} in depth, here we renormalize [[SSC/Structure/BiPolytropes/Analytic51|our original construction]] of bipolytropic models with <math>(n_c, n_e) = (5, 1)</math> such that both entropy values, <math>(K_c, K_e)</math>, are held fixed along each model sequence.


==Original Derivation==
===Throughout the Core===
Drawing from our original derivation, [[SSC/Structure/BiPolytropes/Analytic51#Step_4:_Throughout_the_core_(0_≤_ξ_≤_ξi)|throughout the core]] &hellip;
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="center" colspan="3">
Specify:  <math>K_c</math> and <math>\rho_0 ~\Rightarrow</math>
  </td>
  <td colspan="2">
&nbsp;
  </td>
</tr>
<tr>
  <td align="right">
<math>\rho</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \theta^{n_c}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-5/2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>P</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{1+1/n_c} \theta^{n_c + 1}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{6/5} \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{1/2} \rho_0^{(1-n_c)/(2n_c)} \xi</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c}{G\rho_0^{4/5}} \biggr]^{1/2} \biggl(\frac{3}{2\pi}\biggr)^{1/2} \xi</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>M_r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>4\pi \biggl[ \frac{(n_c + 1)K_c}{4\pi G} \biggr]^{3/2} \rho_0^{(3-n_c)/(2n_c)} \biggl(-\xi^2 \frac{d\theta}{d\xi} \biggr)</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5} } \biggr]^{1/2} \biggl( \frac{2\cdot 3}{\pi } \biggr)^{1/2} \biggl[ \xi^3 \biggl( 1 + \frac{1}{3}\xi^2 \biggr)^{-3/2} \biggr]</math>
  </td>
</tr>
</table>
</div>
===Throughout the Envelope===
And [[SSC/Structure/BiPolytropes/Analytic51#Step_8:_Throughout_the_envelope_(ηi_≤_η_≤_ηs)|throughout the envelope]],
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="center" colspan="5">
&nbsp;
  </td>
  <td align="left" colspan="2">
Knowing:  <math>K_e/K_c</math> and <math>\rho_e/\rho_0</math> from Step 5 &nbsp; <math>\Rightarrow</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\rho</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_e \phi^{n_e}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \biggl(\frac{\rho_e}{\rho_0}\biggr) \phi</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0 \biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{5}_i \phi</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>P</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_e \rho_e^{1+1/n_e} \phi^{n_e + 1}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{6/5} \biggl(\frac{K_e \rho_0^{4/5}}{K_c}\biggr) \biggl(\frac{\rho_e}{\rho_0}\biggr)^{2} \phi^{2}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>K_c \rho_0^{6/5}  \theta^{6}_i \phi^{2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{1/2} \rho_e^{(1-n_e)/(2n_e)} \eta</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c}{G \rho_0^{4/5}} \biggr]^{1/2} \biggl( \frac{K_e \rho_0^{4/5}}{K_c} \biggr)^{1/2} (2\pi)^{-1/2}\eta</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c}{G \rho_0^{4/5}} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-1} \theta^{-2}_i (2\pi)^{-1/2}\eta</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>M_r</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>4\pi \biggl[ \frac{(n_e + 1)K_e}{4\pi G} \biggr]^{3/2} \rho_e^{(3-n_e)/(2n_e)} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2} \biggl( \frac{K_e \rho_0^{4/5}}{K_c} \biggr)^{3/2} \biggl(\frac{\rho_e}{\rho_0}\biggr) \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{K_c^3}{G^3 \rho_0^{2/5}} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-1}_i \biggl( \frac{2}{\pi} \biggr)^{1/2} \biggl(-\eta^2 \frac{d\phi}{d\eta} \biggr)</math>
  </td>
</tr>
</table>
</div>
===Interface Conditions===
And [[SSC/Structure/BiPolytropes/Analytic51#Step_5:_Interface_Conditions|at the interface]] &hellip;
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td colspan="3">
&nbsp;
  </td>
  <td align="left" colspan="2">
Setting <math>n_c=5</math>, <math>n_e=1</math>, and <math>\phi_i = 1 ~~~~\Rightarrow</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\frac{\rho_e}{\rho_0}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{n_c}_i \phi_i^{-n_e}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl( \frac{\mu_e}{\mu_c} \biggr) \theta^{5}_i </math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\biggl( \frac{K_e}{K_c} \biggr) </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0^{1/n_c - 1/n_e}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-(1+1/n_e)} \theta^{1 - n_c/n_e}_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0^{-4/5}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-4}_i</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\frac{\eta_i}{\xi_i}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{n_c + 1}{n_e+1} \biggr]^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \theta_i^{(n_c-1)/2} \phi_i^{(1-n_e)/2}</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>3^{1/2} \biggl( \frac{\mu_e}{\mu_c}\biggr) \theta_i^{2}</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\biggl( \frac{d\phi}{d\eta} \biggr)_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\biggl[ \frac{n_c + 1}{n_e + 1} \biggr]^{1/2} \theta_i^{- (n_c + 1)/2} \phi_i^{(n_e+1)/2} \biggl( \frac{d\theta}{d\xi} \biggr)_i</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>3^{1/2} \theta_i^{- 3} \biggl( \frac{d\theta}{d\xi} \biggr)_i</math>
  </td>
</tr>
</table>
</div>
==New Normalization==
From one of the interface conditions, we see that,
<table border="0" align="center" cellpadding="8">
<tr>
  <td align="right">
<math>\biggl( \frac{K_e}{K_c} \biggr) </math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>\rho_0^{-4/5}\biggl( \frac{\mu_e}{\mu_c} \biggr)^{-2} \theta^{-4}_i</math>
  </td>
</tr>
<tr>
  <td align="right">
<math>\Rightarrow ~~~  \rho_0</math>
  </td>
  <td align="center">
<math>=</math>
  </td>
  <td align="left">
<math>
\biggl[ \biggl( \frac{K_e}{K_c} \biggr)^{-1} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{- 2} \theta^{-4}_i \biggr]^{5 / 4}
=
\biggl[ \biggl( \frac{K_e}{K_c} \biggr)^{-5 / 4} \biggl( \frac{\mu_e}{\mu_c} \biggr)^{-5 / 2} \theta^{-5}_i \biggr]
\, .
</math>
  </td>
</tr>
</table>
Hence, throughout the core, we have,


=See Also=
=See Also=

Revision as of 20:14, 10 November 2023


BiPolytrope with nc = 5 and ne = 1

After studying 📚 S. Yabushita (1975, MNRAS, Vol. 172, pp. 441 - 453) in depth, here we renormalize our original construction of bipolytropic models with (nc,ne)=(5,1) such that both entropy values, (Kc,Ke), are held fixed along each model sequence.

Original Derivation

Throughout the Core

Drawing from our original derivation, throughout the core

Specify: Kc and ρ0

 

ρ

=

ρ0θnc

=

ρ0(1+13ξ2)5/2

P

=

Kcρ01+1/ncθnc+1

=

Kcρ06/5(1+13ξ2)3

r

=

[(nc+1)Kc4πG]1/2ρ0(1nc)/(2nc)ξ

=

[KcGρ04/5]1/2(32π)1/2ξ

Mr

=

4π[(nc+1)Kc4πG]3/2ρ0(3nc)/(2nc)(ξ2dθdξ)

=

[Kc3G3ρ02/5]1/2(23π)1/2[ξ3(1+13ξ2)3/2]

Throughout the Envelope

And throughout the envelope,

 

Knowing: Ke/Kc and ρe/ρ0 from Step 5  

ρ

=

ρeϕne

=

ρ0(ρeρ0)ϕ

=

ρ0(μeμc)θi5ϕ

P

=

Keρe1+1/neϕne+1

=

Kcρ06/5(Keρ04/5Kc)(ρeρ0)2ϕ2

=

Kcρ06/5θi6ϕ2

r

=

[(ne+1)Ke4πG]1/2ρe(1ne)/(2ne)η

=

[KcGρ04/5]1/2(Keρ04/5Kc)1/2(2π)1/2η

=

[KcGρ04/5]1/2(μeμc)1θi2(2π)1/2η

Mr

=

4π[(ne+1)Ke4πG]3/2ρe(3ne)/(2ne)(η2dϕdη)

=

[Kc3G3ρ02/5]1/2(Keρ04/5Kc)3/2(ρeρ0)(2π)1/2(η2dϕdη)

=

[Kc3G3ρ02/5]1/2(μeμc)2θi1(2π)1/2(η2dϕdη)

Interface Conditions

And at the interface

 

Setting nc=5, ne=1, and ϕi=1

ρeρ0

=

(μeμc)θincϕine

=

(μeμc)θi5

(KeKc)

=

ρ01/nc1/ne(μeμc)(1+1/ne)θi1nc/ne

=

ρ04/5(μeμc)2θi4

ηiξi

=

[nc+1ne+1]1/2(μeμc)θi(nc1)/2ϕi(1ne)/2

=

31/2(μeμc)θi2

(dϕdη)i

=

[nc+1ne+1]1/2θi(nc+1)/2ϕi(ne+1)/2(dθdξ)i

=

31/2θi3(dθdξ)i

New Normalization

From one of the interface conditions, we see that,

(KeKc)

=

ρ04/5(μeμc)2θi4

ρ0

=

[(KeKc)1(μeμc)2θi4]5/4=[(KeKc)5/4(μeμc)5/2θi5].

Hence, throughout the core, we have,

See Also


Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |