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

From JETohlineWiki
Jump to navigation Jump to search
Joel2 (talk | contribs)
Created page with "__FORCETOC__ =Main Sequence to Red Giant to Planetary Nebula (Part 2)= <table border="1" align="center" width="100%" colspan="8"> <tr> <td align="center" bgcolor="lightblue" width="25%"><br />Part I:  Background & Objective   </td> <td align="center" bgcolor="lightblue" width="25%"><br />Part II:    </td> <td align="center" bgcolor="lightblue" wi..."
 
Joel2 (talk | contribs)
Line 18: Line 18:
</table>
</table>


==Preface==
==Foundation==
In an [[SSC/Perturbations#2ndOrderODE|accompanying discussion]], we derived the so-called,


<div align="center" id="2ndOrderODE">
<font color="#770000">'''Adiabatic Wave''' (or ''Radial Pulsation'') '''Equation'''</font><br />
{{Math/EQ_RadialPulsation01}}
</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,
<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">
<tr>
  <td align="right">
<math>~g_0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<math>~\frac{GM(r_0)}{r_0^2}</math>
  </td>
</tr>
<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>
  </td>
</tr>
</table>
Hence, multiplying the LAWE through by <math>~(K_c/G)\rho_c^{-4 / 5}</math> gives,
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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}
+ \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>
</tr>
</table>
<table border="0" cellpadding="5" align="center">
<tr>
  <td align="right">
<math>~0</math>
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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>
</tr>
<tr>
  <td align="right">
&nbsp;
  </td>
  <td align="center">
<math>~=</math>
  </td>
  <td align="left">
<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>
</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>


=Related Discussions=
=Related Discussions=


{{ SGFfooter }}
{{ SGFfooter }}

Revision as of 14:18, 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. 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,

g0

=

GM(r0)r02

 

=

G[Mr*ρc1/5(KcG)3/2][r*ρc2/5(KcG)1/2]2

 

=

GMr*(r*)2[ρc3/5(KcG)1/2].

Hence, multiplying the LAWE through by (Kc/G)ρc4/5 gives,

0

=

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


0

=

d2xdr*2+[4r*ρc2/5(KcG)1/2(g0ρ0P0)]dxdr*+ρc4/5(KcG)(ρ0γgP0)[ω2+(43γg)g0r0]x

 

=

d2xdr*2+{4r*ρc2/5(KcG)1/2GMr*(r*)2[ρc3/5(KcG)1/2][ρcρ*P*Kcρc6/5]}dxdr*+ρc4/5(KcG)[ρcρ*γgP*Kcρc6/5]{ω2+(43γg)GMr*(r*)2[ρc3/5(KcG)1/2]ρc2/5r*(GKc)1/2}x

 

=

d2xdr*2+{4r*Mr*(r*)2[ρ*P*]}dxdr*+(1Gρc)[ρ*γgP*]{ω2+(43γg)GρcMr*(r*)3}x

 

=

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

 

=

d2xdr*2+{4(ρ*P*)Mr*(r*)}1r*dxdr*+(ρ*P*){2πσc23γgαgMr*(r*)3}x.

Related Discussions

Tiled Menu

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