Difference between revisions of "User:Tohline/Appendix/Ramblings/SphericalWaveEquation"
(→Assembling the Key Relations: Work toward establishing eigenvalue problem in terms of W) |
|||
Line 103: | Line 103: | ||
</td> | </td> | ||
<td align="left"> | <td align="left"> | ||
<math>~\frac{P_1}{\rho_0} = \biggl(\frac{P_0}{\rho_0}\biggr)</math> | <math>~\frac{P_1}{\rho_0} = \biggl(\frac{P_0}{\rho_0}\biggr) p \, .</math> | ||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
The second expression then becomes, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~x(4g_0 + \omega^2 r_0)</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{P_0}{\rho_0} \frac{d}{dr_0}\biggl(\frac{W\rho_0}{P_0}\biggr) - \biggl(\frac{g_0 \rho_0}{P_0}\biggr)W</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{dW}{dr_0} | |||
+ \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} | |||
- \frac{W }{P_0} \frac{dP_0}{dr_0} | |||
- \biggl(\frac{g_0 \rho_0}{P_0}\biggr)W</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{dW}{dr_0}+ \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Taking the derivative of this expression with respect to <math>~r_0</math> gives, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{dx}{dr_0}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{d}{dr_0}\biggl\{ | |||
(4g_0 + \omega^2 r_0)^{-1}\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
(4g_0 + \omega^2 r_0)^{-1}\frac{d}{dr_0} | |||
\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] | |||
+\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr]\frac{d}{dr_0} | |||
(4g_0 + \omega^2 r_0)^{-1} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
(4g_0 + \omega^2 r_0)^{-1} \biggl\{ | |||
\frac{d^2W}{dr^2_0} + \frac{d}{dr_0}\biggl[\frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} | |||
-(4g_0 + \omega^2 r_0)^{-2}\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] | |||
\biggl\{ 4\frac{dg_0}{dr_0} + \omega^2 | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow~~~~ | |||
(4g_0 + \omega^2 r_0)^{2} \biggl[ \frac{dx}{dr_0} \biggr] | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
(4g_0 + \omega^2 r_0)\biggl\{ | |||
\frac{d^2W}{dr^2_0} + \frac{d}{dr_0}\biggl[\frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} | |||
-\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] | |||
\biggl\{ 4\frac{dg_0}{dr_0} + \omega^2 | |||
\biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Hence, the linearized equation of continuity becomes, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~- (4g_0 + \omega^2 r_0)^{2}\biggl(\frac{W\rho_0}{\gamma_g r_0P_0}\biggr) </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~(4g_0 + \omega^2 r_0)^{2} \biggl[ \frac{dx}{dr_0} \biggr] +\frac{3 (4g_0 + \omega^2 r_0)}{r_0} \biggl[ (4g_0 + \omega^2 r_0)x \biggr] </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
(4g_0 + \omega^2 r_0)\biggl\{ | |||
\frac{d^2W}{dr^2_0} + \frac{d}{dr_0}\biggl[\frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} | |||
-\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] | |||
\biggl\{ 4\frac{dg_0}{dr_0} + \omega^2 | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+\frac{3 (4g_0 + \omega^2 r_0)}{r_0} \biggl[ \frac{dW}{dr_0}+ \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] | |||
</math> | |||
</td> | </td> | ||
</tr> | </tr> |
Revision as of 01:21, 15 May 2016
Playing With Spherical Wave Equation
| Tiled Menu | Tables of Content | Banner Video | Tohline Home Page | |
The traditional presentation of the (spherically symmetric) adiabatic wave equation focuses on fractional radial displacements, <math>~x \equiv \delta r/r_0</math>, of spherical mass shells. After studying in depth various stability analyses of Papaloizou-Pringle tori, I have begun to wonder whether the wave equation for spherical polytropes might look simpler if we focus, instead, on fluctuations in the fluid entropy.
Assembling the Key Relations
In the traditional approach, the following three linearized equations describe the physical relationship between the three dimensionless perturbation amplitudes <math>~p(r_0) \equiv P_1/P_0</math>, <math>~d(r_0) \equiv \rho_1/\rho_0</math> and <math>~x(r_0) \equiv r_1/r_0</math>, for various characteristic eigenfrequencies, <math>~\omega</math>:
Linearized Linearized Linearized |
Let's switch from the perturbation variable, <math>~p</math>, to an enthalpy-related variable,
<math>~W</math> |
<math>~\equiv</math> |
<math>~\frac{P_1}{\rho_0} = \biggl(\frac{P_0}{\rho_0}\biggr) p \, .</math> |
The second expression then becomes,
<math>~x(4g_0 + \omega^2 r_0)</math> |
<math>~=</math> |
<math>~\frac{P_0}{\rho_0} \frac{d}{dr_0}\biggl(\frac{W\rho_0}{P_0}\biggr) - \biggl(\frac{g_0 \rho_0}{P_0}\biggr)W</math> |
|
<math>~=</math> |
<math>~\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} - \frac{W }{P_0} \frac{dP_0}{dr_0} - \biggl(\frac{g_0 \rho_0}{P_0}\biggr)W</math> |
|
<math>~=</math> |
<math>~\frac{dW}{dr_0}+ \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \, . </math> |
Taking the derivative of this expression with respect to <math>~r_0</math> gives,
<math>~\frac{dx}{dr_0}</math> |
<math>~=</math> |
<math>~\frac{d}{dr_0}\biggl\{ (4g_0 + \omega^2 r_0)^{-1}\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} </math> |
|
<math>~=</math> |
<math>~ (4g_0 + \omega^2 r_0)^{-1}\frac{d}{dr_0} \biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] +\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr]\frac{d}{dr_0} (4g_0 + \omega^2 r_0)^{-1} </math> |
|
<math>~=</math> |
<math>~ (4g_0 + \omega^2 r_0)^{-1} \biggl\{ \frac{d^2W}{dr^2_0} + \frac{d}{dr_0}\biggl[\frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} -(4g_0 + \omega^2 r_0)^{-2}\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggl\{ 4\frac{dg_0}{dr_0} + \omega^2 \biggr\} </math> |
<math>~\Rightarrow~~~~ (4g_0 + \omega^2 r_0)^{2} \biggl[ \frac{dx}{dr_0} \biggr] </math> |
<math>~=</math> |
<math>~ (4g_0 + \omega^2 r_0)\biggl\{ \frac{d^2W}{dr^2_0} + \frac{d}{dr_0}\biggl[\frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} -\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggl\{ 4\frac{dg_0}{dr_0} + \omega^2 \biggr\} \, . </math> |
Hence, the linearized equation of continuity becomes,
<math>~- (4g_0 + \omega^2 r_0)^{2}\biggl(\frac{W\rho_0}{\gamma_g r_0P_0}\biggr) </math> |
<math>~=</math> |
<math>~(4g_0 + \omega^2 r_0)^{2} \biggl[ \frac{dx}{dr_0} \biggr] +\frac{3 (4g_0 + \omega^2 r_0)}{r_0} \biggl[ (4g_0 + \omega^2 r_0)x \biggr] </math> |
|
<math>~=</math> |
<math>~ (4g_0 + \omega^2 r_0)\biggl\{ \frac{d^2W}{dr^2_0} + \frac{d}{dr_0}\biggl[\frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggr\} -\biggl[\frac{dW}{dr_0} + \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] \biggl\{ 4\frac{dg_0}{dr_0} + \omega^2 \biggr\} </math> |
|
|
<math>~ +\frac{3 (4g_0 + \omega^2 r_0)}{r_0} \biggl[ \frac{dW}{dr_0}+ \frac{W}{\rho_0} \frac{d\rho_0}{dr_0} \biggr] </math> |
See Also
© 2014 - 2021 by Joel E. Tohline |