SSC/Stability/BiPolytropes/RedGiantToPN/Pt2: Difference between revisions

From JETohlineWiki
Jump to navigation Jump to search
Joel2 (talk | contribs)
Joel2 (talk | contribs)
Line 27: Line 27:
</div>
</div>


whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations.  Assuming that the underlying equilibrium structure is that of a bipolytrope having <math>~(n_c, n_e) = (5, 1)</math>, it makes sense to adopt the normalizations used when defining the equilibrium structure, namely,
whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations.  Multiplying this LAWE through by <math>(K_c/G)\rho_c^{-4 / 5}</math> and recognizing that,
<div align="center">
<table border="0" cellpadding="3">
<tr>
  <td align="right">
<math>~\rho^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{\rho_0}{\rho_c}</math>
  </td>
 
  <td align="center">; &nbsp;&nbsp;&nbsp;</td>
 
  <td align="right">
<math>~r^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{r_0}{[K_c^{1/2}/(G^{1/2}\rho_c^{2/5})]}</math>
  </td>
</tr>
 
<tr>
  <td align="right">
<math>~P^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{P_0}{K_c\rho_c^{6/5}}</math>
  </td>
 
  <td align="center">; &nbsp;&nbsp;&nbsp;</td>
 
  <td align="right">
<math>~M_r^*</math>
  </td>
  <td align="center">
<math>~\equiv</math>
  </td>
  <td align="left">
<math>~\frac{M_r}{[K_c^{3/2}/(G^{3/2}\rho_c^{1/5})]}</math>
  </td>
</tr>
 
</table>
</div>
 
We note as well that,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


Line 92: Line 38:
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~\frac{GM(r_0)}{r_0^2}</math>
<math>~\frac{GM(r_0)}{r_0^2}
  </td>
=
</tr>
\frac{G M_r^*}{(r^*)^2} \biggl[ \rho_c^{3 / 5} \biggl( \frac{K_c}{G}\biggr)^{1 / 2} \biggr]  
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
G \biggl[ M_r^* \rho_c^{-1 / 5} \biggl( \frac{K_c}{G}\biggr)^{3 / 2} \biggr] \biggl[ r^* \rho_c^{-2 / 5}\biggl( \frac{K_c}{G}\biggr)^{1 / 2} \biggr]^{-2}
</math>
  </td>
</tr>
 
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{G M_r^*}{(r^*)^2} \biggl[ \rho_c^{3 / 5} \biggl( \frac{K_c}{G}\biggr)^{1 / 2} \biggr] \, .
</math>
</math>
   </td>
   </td>
Line 125: Line 46:
</table>
</table>


Hence, multiplying the LAWE through by <math>~(K_c/G)\rho_c^{-4 / 5}</math> gives,
we have,
 
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">


Line 138: Line 58:
   <td align="left">
   <td align="left">
<math>~
<math>~
\frac{d^2x}{dr_0^2} + \biggl[\frac{4}{r_0} - \biggl(\frac{g_0 \rho_0}{P_0}\biggr) \biggr] \frac{dx}{dr_0}  
\frac{d^2x}{dr*^2} + \biggl\{ \frac{4}{r^*} -\biggl(\frac{\rho^*}{P^*}\biggr)\frac{ M_r^*}{(r^*)^2}\biggr\} \frac{dx}{dr*}  
+ \biggl(\frac{\rho_0}{\gamma_\mathrm{g} P_0} \biggr)\biggl[\omega^2 + (4 - 3\gamma_\mathrm{g})\frac{g_0}{r_0} \biggr] x  
+ \biggl(\frac{\rho^*}{ P^* } \biggr)\biggl\{ \frac{\omega^2}{\gamma_\mathrm{g} G\rho_c} + \biggl(\frac{4}{\gamma_\mathrm{g}} - 3\biggr)\frac{ M_r^*}{(r^*)^3}\biggr\} x  
</math>
</math>
   </td>
   </td>
Line 146: Line 66:




 
Assuming that the underlying equilibrium structure is that of a bipolytrope having <math>~(n_c, n_e) = (5, 1)</math>, it makes sense to adopt the normalizations used when defining the equilibrium structure, namely,
<table border="0" cellpadding="5" align="center">
<div align="center">
 
<table border="0" cellpadding="3">
<tr>
<tr>
   <td align="right">
   <td align="right">
<math>~0</math>
<math>~\rho^*</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{\rho_0}{\rho_c}</math>
\frac{d^2x}{dr*^2} + \biggl[\frac{4}{r^*} ~-~ \rho_c^{-2 / 5} \biggl( \frac{K_c}{G} \biggr)^{1 / 2}\biggl(\frac{g_0 \rho_0}{P_0}\biggr) \biggr] \frac{dx}{dr*}
~+~ \rho_c^{-4 / 5}\biggl( \frac{K_c}{G} \biggr)\biggl(\frac{\rho_0}{\gamma_\mathrm{g} P_0} \biggr)\biggl[\omega^2 ~+~ (4 - 3\gamma_\mathrm{g})\frac{g_0}{r_0} \biggr]  x
</math>
   </td>
   </td>
</tr>


<tr>
  <td align="center">; &nbsp;&nbsp;&nbsp;</td>
 
   <td align="right">
   <td align="right">
&nbsp;
<math>~r^*</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{r_0}{[K_c^{1/2}/(G^{1/2}\rho_c^{2/5})]}</math>
\frac{d^2x}{dr*^2} ~+~ \biggl\{\frac{4}{r^*} ~-~ \rho_c^{-2 / 5} \biggl( \frac{K_c}{G} \biggr)^{1 / 2} \frac{G M_r^*}{(r^*)^2} \biggl[ \rho_c^{3 / 5} \biggl( \frac{K_c}{G}\biggr)^{1 / 2} \biggr] \biggl[ \frac{\rho_c \rho^*}{P^* K_c \rho_c^{6/5}}\biggr]  \biggr\} \frac{dx}{dr*}
~+~ \rho_c^{-4 / 5}\biggl( \frac{K_c}{G} \biggr)\biggl[ \frac{\rho_c \rho^*}{\gamma_\mathrm{g}P^* K_c \rho_c^{6/5}}\biggr]
\biggl\{\omega^2 ~+~ (4 - 3\gamma_\mathrm{g}) \frac{G M_r^*}{(r^*)^2} \biggl[ \rho_c^{3 / 5} \biggl( \frac{K_c}{G}\biggr)^{1 / 2} \biggr] \frac{\rho_c^{2 / 5}}{r^*}\biggl( \frac{G}{K_c}\biggr)^{1 / 2} \biggr\} x
</math>
   </td>
   </td>
</tr>
</tr>
Line 182: Line 95:
<tr>
<tr>
   <td align="right">
   <td align="right">
&nbsp;
<math>~P^*</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{P_0}{K_c\rho_c^{6/5}}</math>
\frac{d^2x}{dr*^2} ~+~ \biggl\{\frac{4}{r^*} ~-~ \frac{M_r^*}{(r^*)^2} \biggl[ \frac{\rho^*}{P^* }\biggr]  \biggr\} \frac{dx}{dr*}
~+~ \biggl( \frac{1}{G\rho_c} \biggr)\biggl[ \frac{ \rho^*}{\gamma_\mathrm{g}P^* }\biggr]
\biggl\{\omega^2 ~+~ (4 - 3\gamma_\mathrm{g}) \frac{G\rho_c M_r^*}{(r^*)^3} \biggr\} x
</math>
   </td>
   </td>
</tr>


<tr>
  <td align="center">; &nbsp;&nbsp;&nbsp;</td>
 
   <td align="right">
   <td align="right">
&nbsp;
<math>~M_r^*</math>
   </td>
   </td>
   <td align="center">
   <td align="center">
<math>~=</math>
<math>~\equiv</math>
   </td>
   </td>
   <td align="left">
   <td align="left">
<math>~
<math>~\frac{M_r}{[K_c^{3/2}/(G^{3/2}\rho_c^{1/5})]}</math>
\frac{d^2x}{dr*^2} + \biggl\{ \frac{4}{r^*} -\biggl(\frac{\rho^*}{P^*}\biggr)\frac{ M_r^*}{(r^*)^2}\biggr\} \frac{dx}{dr*}
+ \biggl(\frac{\rho^*}{ P^* } \biggr)\biggl\{ \frac{\omega^2}{\gamma_\mathrm{g} G\rho_c} + \biggl(\frac{4}{\gamma_\mathrm{g}} - 3\biggr)\frac{ M_r^*}{(r^*)^3}\biggr\} x
</math>
   </td>
   </td>
</tr>
</tr>


<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~
\frac{d^2x}{dr*^2} + \biggl\{ 4 -\biggl(\frac{\rho^*}{P^*}\biggr)\frac{ M_r^*}{(r^*)}\biggr\}\frac{1}{r^*} \frac{dx}{dr*}
+ \biggl(\frac{\rho^*}{ P^* } \biggr)\biggl\{ \frac{2\pi \sigma_c^2}{3\gamma_\mathrm{g}}  ~-~\frac{\alpha_\mathrm{g} M_r^*}{(r^*)^3}\biggr\}  x \, .
</math>
  </td>
</tr>
</table>
</table>
</div>
We note as well that,


=Related Discussions=
=Related Discussions=


{{ SGFfooter }}
{{ SGFfooter }}

Revision as of 14:29, 25 December 2025

Main Sequence to Red Giant to Planetary Nebula (Part 2)


Part I:  Background & Objective

 


Part II: 

 


Part III: 

 


Part IV: 

 

Foundation

In an accompanying discussion, we derived the so-called,

Adiabatic Wave (or Radial Pulsation) Equation

d2xdr02+[4r0(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(43γg)g0r0]x=0

whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations. Multiplying this LAWE through by (Kc/G)ρc4/5 and recognizing that,

g0

=

GM(r0)r02=GMr*(r*)2[ρc3/5(KcG)1/2]

we have,

0

=

d2xdr*2+{4r*(ρ*P*)Mr*(r*)2}dxdr*+(ρ*P*){ω2γgGρc+(4γg3)Mr*(r*)3}x


Assuming that the underlying equilibrium structure is that of a bipolytrope having (nc,ne)=(5,1), it makes sense to adopt the normalizations used when defining the equilibrium structure, namely,

ρ*

ρ0ρc

;    

r*

r0[Kc1/2/(G1/2ρc2/5)]

P*

P0Kcρc6/5

;    

Mr*

Mr[Kc3/2/(G3/2ρc1/5)]

We note as well that,

Related Discussions

Tiled Menu

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