Difference between revisions of "User:Tohline/SSC/Stability/n3PolytropeLAWE"

From VistrailsWiki
Jump to navigation Jump to search
(Created page with '__FORCETOC__ <!-- __NOTOC__ will force TOC off --> =Radial Oscillations of n = 3 Polytropic Spheres= {{LSU_HBook_header}} ==Groundwork== In an [[User:Tohline/SSC/Perturbations…')
 
Line 66: Line 66:
==Numerical Integration==
==Numerical Integration==


It can be shown straightforwardly that this matches the LAWE used by [[User:Tohline/SSC/Perturbations#Schwarzschild_.281941.29|Schwarzschild (1941)]], if <math>~n</math> is set to 3.  Here we use the [[User:Tohline/SSC/Stability/Polytropes#Numerical_Integration_from_the_Center.2C_Outward|finite-difference algorithm described separately]] to integrate the discretized LAWE from the center of the polytropic configuration, outward to its surface, which in this case is located at the polytropic-coordinate location,
It can be shown straightforwardly that this matches the LAWE used by [[User:Tohline/SSC/Perturbations#Schwarzschild_.281941.29|Schwarzschild (1941)]], if <math>~n</math> is set to 3.  Here we use the [[User:Tohline/SSC/Stability/Polytropes#Numerical_Integration_from_the_Center.2C_Outward|finite-difference algorithm described separately]] to integrate the discretized LAWE from the center of the polytropic configuration, outward to its surface, which in this case &#8212; see, for example, p. 77 of Horedt &#8212; is located at the polytropic-coordinate location,
<div align="center">
<div align="center">
<math>~\xi_\mathrm{max} = 6.89684862 \, .</math>
<math>~\xi_\mathrm{max} = 6.89684862 \, .</math>
</div>
</div>


The algorithm is as follows:
<ul>
<ul>
   <li>Establish an equally spaced radial-coordinate grid containing <math>~N</math> grid zones (and, accordingly, <math>~N+1</math> grid ''lines''), in which case the grid-spacing parameter, <math>~\Delta_\xi \equiv \xi_\mathrm{max}/N</math>. </li>
   <li>Establish an equally spaced radial-coordinate grid containing <math>~N</math> grid zones (and, accordingly, <math>~N+1</math> grid ''lines''), in which case the grid-spacing parameter, <math>~\Delta_\xi \equiv \xi_\mathrm{max}/N</math>. </li>
Line 82: Line 83:
   <li>Determine the value of the eigenfunction at the first grid ''line'' &#8212; having coordinate location, <math>~\xi_1 = \Delta_\xi </math> &#8212; via the expression,</li>
   <li>Determine the value of the eigenfunction at the first grid ''line'' &#8212; having coordinate location, <math>~\xi_1 = \Delta_\xi </math> &#8212; via the expression,</li>
</ul>
</ul>


=Related Discussions=
=Related Discussions=

Revision as of 18:07, 17 February 2017

Radial Oscillations of n = 3 Polytropic Spheres

Whitworth's (1981) Isothermal Free-Energy Surface
|   Tiled Menu   |   Tables of Content   |  Banner Video   |  Tohline Home Page   |

Groundwork

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

Adiabatic Wave (or Radial Pulsation) Equation

LSU Key.png

<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 = 0 </math>

whose solution gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations. Because this widely used form of the radial pulsation equation is not dimensionless but, rather, has units of inverse length-squared, we have found it useful to also recast it in the following dimensionless form:

<math> \frac{d^2x}{d\chi_0^2} + \biggl[\frac{4}{\chi_0} - \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} \biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \biggr] \frac{dx}{d\chi_0} + \biggl(\frac{\rho_0}{\rho_c}\biggr) \biggl(\frac{P_0}{P_c}\biggr)^{-1} \biggl(\frac{1}{\gamma_\mathrm{g}} \biggr)\biggl[\tau_\mathrm{SSC}^2 \omega^2 + (4 - 3\gamma_\mathrm{g})\biggl(\frac{g_0}{g_\mathrm{SSC}}\biggr) \frac{1}{\chi_0} \biggr] x = 0 , </math>

where,

<math>~g_\mathrm{SSC} \equiv \frac{P_c}{R\rho_c} \, ,</math>       and       <math>~\tau_\mathrm{SSC} \equiv \biggl[\frac{R^2 \rho_c}{P_c}\biggr]^{1/2} \, .</math>

In a separate discussion, we showed that specifically for isolated, polytropic configurations, this linear adiabatic wave equation (LAWE) can be rewritten as,

<math>~\frac{d^2x}{d\xi^2} + \biggl[\frac{4 - (n+1)V(\xi)}{\xi} \biggr] \frac{dx}{d\xi} + \biggl[\omega^2 \biggl(\frac{a_n^2 \rho_c }{\gamma_g P_c} \biggr) \frac{\theta_c}{\theta} - \biggl(3-\frac{4}{\gamma_g}\biggr) \cdot \frac{(n+1)V(x)}{\xi^2} \biggr] x </math>

<math>~=</math>

<math>0 \, ,</math>

where we have adopted the function notation,

<math>~V(\xi)</math>

<math>~\equiv</math>

<math>~- \frac{\xi}{\theta} \frac{d \theta}{d\xi} \, .</math>

Here we perform a numerical integration of the governing LAWE for <math>~n=3</math> polytropes. We can directly compare our results with Schwarzschild's (1941) published work on "Overtone Pulsations for the Standard [Stellar] Model."

Numerical Integration

It can be shown straightforwardly that this matches the LAWE used by Schwarzschild (1941), if <math>~n</math> is set to 3. Here we use the finite-difference algorithm described separately to integrate the discretized LAWE from the center of the polytropic configuration, outward to its surface, which in this case — see, for example, p. 77 of Horedt — is located at the polytropic-coordinate location,

<math>~\xi_\mathrm{max} = 6.89684862 \, .</math>

The algorithm is as follows:

  • Establish an equally spaced radial-coordinate grid containing <math>~N</math> grid zones (and, accordingly, <math>~N+1</math> grid lines), in which case the grid-spacing parameter, <math>~\Delta_\xi \equiv \xi_\mathrm{max}/N</math>.
  • Specify a value of the adiabatic exponent, <math>~\gamma</math>, which, in turn, determines the value of the parameter, <math>~\alpha \equiv (3-4/\gamma) \, .</math>
  • Choose a value for the dimensionless oscillation frequency parameter,

    <math>~\mathfrak{F} = \frac{\sigma_c^2}{\gamma} - 2\alpha \, .</math>

  • Set the central value of the eigenfunction to unity, that is, set <math>~x_0 = 1</math>.
  • Determine the value of the eigenfunction at the first grid line — having coordinate location, <math>~\xi_1 = \Delta_\xi </math> — via the expression,

Related Discussions