User:Tohline/Appendix/Ramblings/NumericallyDeterminedEigenvectors
Numerically Determined Eigenvectors of a Zero-Zero Bipolytrope
Here we build on the analytic foundation summarized in an accompanying chapter and attempt to numerically construct a variety of eigenvectors that describe radial oscillations of bipolytropes for which, <math>~(n_c, n_e) = (0,0)</math>.
|
|---|
| | Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Setup
We'll begin with the linear-adiabatic wave equations that describe oscillations of the core and envelope, separately. We also will immediately restrict our investigation to configurations for which,
<math>~g^2 = \mathcal{B} </math> <math>~\Rightarrow</math> <math>~g^2 = \frac{1+8q^3}{ (1+2q^3)^2 } \, ,</math> and, <math>~q^3 = \mathcal{D} = \biggl[ \frac{\rho_e/\rho_c}{2(1-\rho_e/\rho_c)} \biggr] </math> <math>~\Rightarrow</math> <math>~\frac{\rho_e}{\rho_c} = \frac{2q^3}{1+2q^3} \, .</math>
For the core we have,
|
<math>~0</math> |
<math>~=</math> |
<math>~ (1 - \eta^2)\frac{d^2x}{d\eta^2} + ( 4 - 6\eta^2 ) \frac{1}{\eta} \cdot \frac{dx}{d\eta} + \mathfrak{F}_\mathrm{core} x \, , </math> |
where,
<math>~\eta \equiv \frac{\xi}{g} \, ,</math> and <math>~\mathfrak{F}_\mathrm{core} \equiv \frac{3\omega_\mathrm{core}^2}{2\pi G\gamma_c \rho_c} - 2\alpha_c\, .</math>
And, for the envelope we have,
|
<math>~0</math> |
<math>~=</math> |
<math>~ ( 1 - q^3 \xi^3 ) \frac{d^2x}{d\xi^2} + ( 3 - 6q^3 \xi^3 ) \frac{1}{\xi} \cdot \frac{dx}{d\xi} + \biggl[ q^3 \mathfrak{F}_\mathrm{env} \xi^3 -\alpha_e \biggr]\frac{x}{\xi^2} \, , </math> |
where,
|
<math>~\mathfrak{F}_\mathrm{env}</math> |
<math>~\equiv</math> |
<math>~\frac{3\omega^2_\mathrm{env}}{2\pi G \gamma_e \rho_e} - 2\alpha_e \, . </math> |
Initial Focus
|
|
|---|
|
© 2014 - 2021 by Joel E. Tohline |