Difference between revisions of "User:Tohline/SSC/Stability/InstabilityOnsetOverview"
Line 42: | Line 42: | ||
</div> | </div> | ||
Yabushita ([http://adsabs.harvard.edu/abs/1974MNRAS.167…95Y 1974], [http://adsabs.harvard.edu/abs/1975MNRAS.172..441Y 1975]) showed that | Yabushita ([http://adsabs.harvard.edu/abs/1974MNRAS.167…95Y 1974], [http://adsabs.harvard.edu/abs/1975MNRAS.172..441Y 1975]) showed that the following eigenvector specification provides a | ||
<div align="center"> | <div align="center"> | ||
<table border="0" cellpadding="5" align="center"> | <table border="0" cellpadding="5" align="center"> | ||
Line 56: | Line 56: | ||
</td> | </td> | ||
<td align="left"> | <td align="left"> | ||
<math>~x = 1 - \biggl( \frac{1}{\xi e^{-\psi}}\biggr) \frac{d\psi}{d\xi} \, | <math>~x = 1 - \biggl( \frac{1}{\xi e^{-\psi}}\biggr) \frac{d\psi}{d\xi} \, ,</math> | ||
</td> | </td> | ||
</tr> | </tr> | ||
</table> | </table> | ||
</div> | </div> | ||
if the adiabatic exponent is assigned the value, <math>~\gamma_g = 1</math>, in which case the parameter, <math>~\alpha = -1</math>. | |||
When viewed in concert with the surface boundary condition, | When viewed in concert with the surface boundary condition, | ||
<div align="center"> | <div align="center"> |
Revision as of 18:49, 19 March 2017
Overview: Marginally Unstable Pressure-Truncated Configurations
Additional details may be found here.
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
Isothermal
The internal structure of a detailed force-balance model is provided via the function, <math>~\psi(\xi)</math>, which is a solution to the,
Isothermal Lane-Emden Equation
|
Equilibrium sequence for pressure-truncated configurations is displayed in three ways.
Figure 1: Bonnor's P-V Diagram |
|
This equation — in the following, slightly rewritten form — can be found among our selected set of key equations associated with the study of radial pulsation, and will henceforth be referred to as the,
Isothermal LAWE
<math>~0 = \frac{d^2x}{d\xi^2} + \biggl[4 - \xi \biggl( \frac{d\psi}{d\xi} \biggr) \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} + \biggl[ \biggl( \frac{\sigma_c^2}{6\gamma_\mathrm{g}}\biggr)\xi^2 - \alpha \xi \biggl( \frac{d\psi}{d\xi} \biggr) \biggr] \frac{x}{\xi^2} </math> |
|
where: <math>~\sigma_c^2 \equiv \frac{3\omega^2}{2\pi G\rho_c}</math> and, <math>~\alpha \equiv \biggl(3 - \frac{4}{\gamma_\mathrm{g}}\biggr)</math> |
Yabushita (1974, 1975) showed that the following eigenvector specification provides a
Precise Solution to the Isothermal LAWE | ||
<math>~\sigma_c^2 = 0</math> |
and |
<math>~x = 1 - \biggl( \frac{1}{\xi e^{-\psi}}\biggr) \frac{d\psi}{d\xi} \, ,</math> |
if the adiabatic exponent is assigned the value, <math>~\gamma_g = 1</math>, in which case the parameter, <math>~\alpha = -1</math>. When viewed in concert with the surface boundary condition,
<math>~\frac{d\ln x}{d\ln\xi}</math> |
<math>~=</math> |
<math>~- 3 \, ,</math> |
the relevant configuration is precisely defined by the surface condition, xxx, which is identical to the configuration at the turning point.
Polytropic
Given a value of the polytropic index, <math>~n</math>, the internal structure of a detailed force-balance model is provided via the function, <math>~\theta(\xi)</math>, which is a solution of the,
Polytropic Lane-Emden Equation
|
subject to the boundary conditions, <math>~\Theta_H = 1</math> and <math>~d\Theta_H/d\xi = 0</math> at <math>~\xi = 0</math>.
To identify pure radial oscillation modes, we seek solutions to the,
Polytropic LAWE
<math>~0 = \frac{d^2x}{d\xi^2} + \biggl[ 4 - (n+1) Q \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} + (n+1) \biggl[ \biggl( \frac{\sigma_c^2}{6\gamma_g } \biggr) \frac{\xi^2}{\theta} - \alpha Q\biggr] \frac{x}{\xi^2} </math> |
|
where: <math>~Q(\xi) \equiv - \frac{d\ln\theta}{d\ln\xi} \, ,</math> <math>~\sigma_c^2 \equiv \frac{3\omega^2}{2\pi G\rho_c} \, ,</math> and, <math>~\alpha \equiv \biggl(3 - \frac{4}{\gamma_\mathrm{g}}\biggr)</math> |
We have discovered that, for any value of the polytropic index in the range, <math>~3 \le n < \infty</math>, the following eigenvector specification provides a
|
if the adiabatic exponent is assigned the value, <math>~\gamma_g = (n+1)/n</math>, in which case the parameter, <math>~\alpha = (3-n)/(n+1)</math>.
References
© 2014 - 2021 by Joel E. Tohline |