SSC/Stability/n1PolytropeLAWE/Pt2: Difference between revisions
Created page with "__FORCETOC__ =Radial Oscillations of n = 1 Polytropic Spheres (Pt 2)= <table border="1" align="center" width="100%" colspan="8"> <tr> <td align="center" rowspan="1" bgcolor="lightblue" width="33%"><br />Part I: Search for Analytic Solutions<br /> </td> <td align="center" rowspan="1" bgcolor="lightblue" width="33%"><br />Part II: New Ideas<br /> </td> <td align="center" ro..." |
No edit summary |
||
| Line 8: | Line 8: | ||
</tr> | </tr> | ||
</table> | </table> | ||
==New Idea Involving Logarithmic Derivatives== | |||
===Simplistic Layout=== | |||
Let's begin, again, with the relevant LAWE, as [[#Attempt_at_Deriving_an_Analytic_Eigenvector_Solution|provided above]]. After dividing through by <math>~x</math>, we have, | |||
<div align="center"> | |||
<math> | |||
(\sin\xi )\frac{\xi^2}{x} \cdot \frac{d^2x}{d\xi^2} + 2 \biggl[ \sin\xi + \xi \cos \xi \biggr] \frac{\xi}{x} \cdot \frac{dx}{d\xi} + | |||
\biggl[ \sigma^2 \xi^3 - 2\alpha ( \sin\xi - \xi \cos \xi ) \biggr] = 0 \, , | |||
</math><br /> | |||
</div> | |||
<br /> | |||
where, | |||
<div align="center"> | |||
<table border="0" cellpadding="5"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\sigma^2</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math> | |||
~\frac{\omega^2}{2\pi G\rho_c \gamma_g} \, , | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\alpha</math> | |||
</td> | |||
<td align="center"> | |||
<math>~\equiv</math> | |||
</td> | |||
<td align="left"> | |||
<math> | |||
~3-\frac{4}{\gamma_g} | |||
\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Now, in addition to recognizing that, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{\xi}{x} \cdot \frac{dx}{d\xi} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{d\ln x}{d\ln \xi} \, ,</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
in a [[SSC/Stability/BiPolytrope00Details#Idea_Involving_Logarithmic_Derivatives|separate context]], we showed that, quite generally, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~ \frac{\xi^2}{x} \cdot \frac{d^2x}{d\xi^2} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d}{d\ln\xi} \biggl[ \frac{d\ln x}{d\ln \xi} \biggr] | |||
- \biggl[ 1 - \frac{d\ln x}{d\ln \xi} \biggr]\cdot \frac{d\ln x}{d\ln \xi} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Hence, if we ''assume'' that the eigenfunction is a power-law of <math>~\xi</math>, that is, ''assume'' that, | |||
<div align="center"> | |||
<math>~x = a_0 \xi^{c_0} \, ,</math> | |||
</div> | |||
then the logarithmic derivative of <math>~x</math> is a constant, namely, | |||
<div align="center"> | |||
<math>~\frac{d\ln x}{d\ln\xi} = c_0 \, ,</math> | |||
</div> | |||
and the two key derivative terms will be, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{\xi}{x} \cdot \frac{dx}{d\xi} = c_0 \, ,</math> | |||
</td> | |||
<td align="center"> | |||
and | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{\xi^2}{x} \cdot \frac{d^2x}{d\xi^2} = c_0(c_0-1) \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
In this case, the LAWE is no longer a differential equation but, instead, takes the form, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~-\sigma^2 \xi^3</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
c_0(c_0-1) \sin\xi + 2c_0 [ \sin\xi + \xi \cos \xi ] - 2\alpha ( \sin\xi - \xi \cos \xi ) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\sin\xi [c_0(c_0-1) +2c_0 -2\alpha ] + \xi \cos \xi [2(c_0+\alpha) ] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\sin\xi [c_0^2 + c_0 -2\alpha ] + \xi \cos \xi [2(c_0+\alpha) ] \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Now, the cosine term will go to zero if <math>~c_0 = -\alpha</math>; and the sine term will go to zero if, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\alpha</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~3 </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~ \gamma_g</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\infty \, . </math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
If these two — rather strange — conditions are met, then we have a marginally unstable configuration because, <math>~\sigma^2 = 0</math>. This, in and of itself, is not very physically interesting. However, it may give us a clue regarding how to more generally search for a physically reasonable radial eigenfunction. | |||
===More general Assumption=== | |||
Try, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\xi^{c_0} \biggl[a_0 + b_0\sin\xi + d_0 \xi\cos\xi \biggr] </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~\frac{dx}{d\xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\xi^{c_0} \frac{d}{d\xi}\biggl[a_0 + b_0\sin\xi + d_0 \xi\cos\xi \biggr] + c_0\xi^{c_0-1} \biggl[a_0 + b_0\sin\xi + d_0 \xi\cos\xi \biggr] </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\xi^{c_0} \biggl[ b_0\cos\xi - d_0 \xi\sin\xi +d_0\cos\xi\biggr] + c_0\xi^{c_0-1} \biggl[a_0 + b_0\sin\xi + d_0 \xi\cos\xi \biggr] </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~\frac{d\ln x}{d\ln \xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\xi \biggl[ b_0\cos\xi - d_0 \xi\sin\xi +d_0\cos\xi\biggr]\biggl[a_0 + b_0\sin\xi + d_0 \xi\cos\xi \biggr]^{-1} | |||
+ c_0 </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\biggl[ (b_0+d_0)\xi\cos\xi - d_0 \xi^2\sin\xi \biggr]\biggl[a_0 + b_0\sin\xi + d_0 \xi\cos\xi \biggr]^{-1} | |||
+ c_0 </math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
===Another Viewpoint=== | |||
====Development==== | |||
Multiplying through the [[#Simplistic_Layout|above LAWE]] by <math>~(x \xi^{-3})</math> gives, | |||
<div align="center"> | |||
<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{\sin\xi }{\xi} \cdot \frac{d^2x}{d\xi^2} + 2 \biggl[\frac{ \sin\xi + \xi \cos \xi }{\xi^2}\biggr] \frac{dx}{d\xi} + | |||
\biggl[ \sigma^2 - 2\alpha \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr]x | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Notice that, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{d}{d\xi}\biggl[\frac{\sin\xi}{\xi}\biggr]</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \frac{\sin\xi}{\xi^2} + \frac{\cos\xi}{\xi} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl[ \frac{\xi\cos\xi - \sin\xi }{\xi^2} \biggr] \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
And, hence, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{d^2}{d\xi^2}\biggl[\frac{\sin\xi}{\xi}\biggr]</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d}{d\xi}\biggl[ \frac{\cos\xi }{\xi} - \frac{\sin\xi }{\xi^2} \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\frac{\cos\xi}{\xi^2} -\frac{\sin\xi}{\xi} + \frac{2\sin\xi}{\xi^3} - \frac{\cos\xi}{\xi^2} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
-\frac{\sin\xi}{\xi} + 2\biggl[ \frac{\sin\xi -\xi\cos\xi}{\xi^3} \biggr] \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
So, we can write, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{d^2}{d\xi^2} \biggl\{ \biggl( \frac{\sin\xi}{\xi}\biggr)x \biggr\}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{d}{d\xi} \biggl\{ | |||
\biggl(\frac{\sin\xi}{\xi}\biggr)\frac{dx}{d\xi} + x\frac{d}{d\xi} \biggl[ \biggl(\frac{\sin\xi}{\xi}\biggr) \biggr] | |||
\biggr\}</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{\sin\xi}{\xi} \cdot \frac{d^2 x}{d\xi^2} | |||
+ 2\frac{dx}{d\xi} \cdot \biggl[\frac{d}{d\xi}\biggr(\frac{\sin\xi}{\xi}\biggr) \biggr] | |||
+ x \cdot \frac{d^2}{d\xi^2} \biggl(\frac{\sin\xi}{\xi}\biggr) | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{\sin\xi}{\xi} \cdot \frac{d^2 x}{d\xi^2} | |||
+ 2\frac{dx}{d\xi} \cdot \biggl[ \frac{\xi\cos\xi - \sin\xi }{\xi^2} \biggr] | |||
+ x \cdot \biggl\{ -\frac{\sin\xi}{\xi} + 2\biggl[ \frac{\sin\xi -\xi\cos\xi}{\xi^3} \biggr] \biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
This means that we can rewrite the LAWE as, | |||
<div align="center"> | |||
<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^2}{d\xi^2} \biggl\{ \biggl( \frac{\sin\xi}{\xi}\biggr)x \biggr\} - 2\frac{dx}{d\xi} \cdot \biggl[ \frac{\xi\cos\xi - \sin\xi }{\xi^2} \biggr] | |||
- x \cdot \biggl\{ -\frac{\sin\xi}{\xi} + 2\biggl[ \frac{\sin\xi -\xi\cos\xi}{\xi^3} \biggr] \biggr\} | |||
+ 2 \biggl[\frac{ \sin\xi + \xi \cos \xi }{\xi^2}\biggr] \frac{dx}{d\xi} + | |||
\biggl[ \sigma^2 - 2\alpha \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr]x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d^2}{d\xi^2} \biggl\{ \biggl( \frac{\sin\xi}{\xi}\biggr)x \biggr\} | |||
+ 4 \biggl[\frac{ \sin\xi }{\xi^2}\biggr] \frac{dx}{d\xi} | |||
+ \biggl\{ \frac{\sin\xi}{\xi} | |||
+ \sigma^2 - 2(1+\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr\} \cdot x \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
We recognize, also, that, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{1}{\xi} \cdot \frac{d}{d\xi}\biggl[ \biggl(\frac{\sin\xi}{\xi} \biggr) x \biggr]</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl[ \frac{\xi\cos\xi - \sin\xi }{\xi^3} \biggr]x + \biggl(\frac{\sin\xi}{\xi^2} \biggr)\frac{dx}{d\xi} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\Rightarrow ~~~ | |||
4\biggl(\frac{\sin\xi}{\xi^2} \biggr)\frac{dx}{d\xi} | |||
</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{4}{\xi} \cdot \frac{d}{d\xi}\biggl[ \biggl(\frac{\sin\xi}{\xi} \biggr) x \biggr] | |||
+ 4\biggl[ \frac{\sin\xi - \xi\cos\xi }{\xi^3} \biggr]x | |||
\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
So the LAWE becomes, | |||
<div align="center"> | |||
<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^2}{d\xi^2} \biggl\{ \biggl( \frac{\sin\xi}{\xi}\biggr)x \biggr\} | |||
+ \frac{4}{\xi} \cdot \frac{d}{d\xi}\biggl[ \biggl(\frac{\sin\xi}{\xi} \biggr) x \biggr] | |||
+ 4\biggl[ \frac{\sin\xi - \xi\cos\xi }{\xi^3} \biggr]x | |||
+ \biggl\{ \frac{\sin\xi}{\xi} | |||
+ \sigma^2 - 2(1+\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr\} \cdot x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d^2}{d\xi^2} \biggl\{ \biggl( \frac{\sin\xi}{\xi}\biggr)x \biggr\} | |||
+ \frac{4}{\xi} \cdot \frac{d}{d\xi}\biggl[ \biggl(\frac{\sin\xi}{\xi} \biggr) x \biggr] | |||
+ \biggl\{ \frac{\sin\xi}{\xi} | |||
+ \sigma^2 + [4- 2(1+\alpha)] \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr\} \cdot x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d^2 \Upsilon}{d\xi^2} | |||
+ \frac{4}{\xi} \cdot \frac{d\Upsilon}{d\xi} | |||
+ \Upsilon | |||
+ \biggl[ \sigma^2 + 2(1-\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr] \cdot x \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
where we have introduced the new, modified eigenfunction, | |||
<div align="center"> | |||
<math>\Upsilon \equiv \biggl( \frac{\sin\xi}{\xi} \biggr) x \, .</math> | |||
</div> | |||
Alternatively, the LAWE may be written as, | |||
<div align="center"> | |||
<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^2 \Upsilon}{d\xi^2} | |||
+ \frac{4}{\xi} \cdot \frac{d\Upsilon}{d\xi} | |||
+ \biggl[ \sigma^2 + 2(1-\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) + \frac{\sin\xi}{\xi} \biggr] \cdot x \, ; | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
or, | |||
<div align="center"> | |||
<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{\xi^2}{\Upsilon} \cdot \frac{d^2 \Upsilon}{d\xi^2} | |||
+ \frac{4\xi}{\Upsilon} \cdot \frac{d\Upsilon}{d\xi} | |||
+ \biggl[ \sigma^2 + 2(1-\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) + \frac{\sin\xi}{\xi} \biggr] \cdot \frac{\xi^3}{\sin\xi} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{\xi^2}{\Upsilon} \cdot \frac{d^2 \Upsilon}{d\xi^2} | |||
+ \frac{4\xi}{\Upsilon} \cdot \frac{d\Upsilon}{d\xi} | |||
+ \biggl[ \sigma^2 \biggl(\frac{\xi^3}{\sin\xi} \biggr) + 2(1-\alpha) \biggl( 1 - \xi \cot \xi \biggr) + \xi^2 \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Now, if we adopt the homentropic convention that arises from setting, <math>~\gamma = (n+1)/n</math>, then for our <math>~n=1</math> polytropic configuration, we should set, <math>~\gamma = 2</math> and, hence, <math>~\alpha = 1</math>. This will mean that the lat term in this LAWE naturally goes to zero. Hence, we have, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~- \sigma^2 x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d^2 \Upsilon}{d\xi^2} | |||
+ \frac{4}{\xi} \cdot \frac{d\Upsilon}{d\xi} + \Upsilon | |||
\, ; | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
or, | |||
<div align="center"> | |||
<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^2 \Upsilon}{d\xi^2} | |||
+ \frac{4}{\xi} \cdot \frac{d\Upsilon}{d\xi} + \biggl[1 + \sigma^2 \biggl(\frac{\xi}{\sin\xi}\biggr) \biggr] \Upsilon \, ; | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
or, | |||
<div align="center"> | |||
<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{\xi^2}{\Upsilon} \cdot \frac{d^2 \Upsilon}{d\xi^2} | |||
+ \frac{4\xi}{\Upsilon} \cdot \frac{d\Upsilon}{d\xi} + \biggl[\xi^2 + \sigma^2 \biggl(\frac{\xi^3}{\sin\xi}\biggr) \biggr] \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Does this help? | |||
====Check for Mistakes==== | |||
Given the definition of <math>~\Upsilon</math>, its first derivative is, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{d\Upsilon}{d\xi} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{dx}{d\xi} +x\biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
and its second derivative is, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{d^2\Upsilon}{d\xi^2} </math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{d}{d\xi} \biggl\{ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{dx}{d\xi} +x\biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{d^2x}{d\xi^2} + 2 \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \cdot \frac{dx}{d\xi} | |||
+ x \cdot \frac{d}{d\xi} \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{d^2x}{d\xi^2} + 2 \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \cdot \frac{dx}{d\xi} | |||
+ x \biggl[ -\frac{\sin\xi}{\xi} - \frac{2\cos\xi}{\xi^2} + \frac{2\sin\xi}{\xi^3} \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Hence, the "upsilon" LAWE becomes, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~-\sigma^2 x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d^2 \Upsilon}{d\xi^2} | |||
+ \frac{4}{\xi} \cdot \frac{d\Upsilon}{d\xi} | |||
+ \Upsilon | |||
+ \biggl[ 2(1-\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr] \cdot x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{d^2x}{d\xi^2} + 2 \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \cdot \frac{dx}{d\xi} | |||
+ x \biggl[ -\frac{\sin\xi}{\xi} - \frac{2\cos\xi}{\xi^2} + \frac{2\sin\xi}{\xi^3} \biggr] | |||
+ \frac{4}{\xi} \cdot \biggl\{ \biggl( \frac{\sin\xi}{\xi} \biggr) \frac{dx}{d\xi} +x\biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \biggr\} | |||
+ \biggl[\frac{\sin\xi}{\xi} + 2(1-\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr] \cdot x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{d^2x}{d\xi^2} + \biggl\{\biggl( \frac{4\sin\xi}{\xi^2} \biggr) + 2 \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \biggr\}\cdot \frac{dx}{d\xi} | |||
+ \biggl[ -\frac{\sin\xi}{\xi} - \frac{2\cos\xi}{\xi^2} + \frac{2\sin\xi}{\xi^3} | |||
+ \frac{4\cos\xi}{\xi^2} - \frac{4\sin\xi}{\xi^3} | |||
+ \frac{\sin\xi}{\xi} + 2(1-\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr] \cdot x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{d^2x}{d\xi^2} + \biggl[ \frac{2\cos\xi}{\xi} + \frac{2\sin\xi}{\xi^2} \biggr] \cdot \frac{dx}{d\xi} | |||
+ \biggl[- 2\biggl( \frac{\sin\xi -\xi\cos\xi}{\xi^3} \biggr) | |||
+ (2-2\alpha) \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr] \cdot x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl( \frac{\sin\xi}{\xi} \biggr) \frac{d^2x}{d\xi^2} + 2\biggl[ \frac{\sin\xi}{\xi^2} + \frac{\cos\xi}{\xi} \biggr] \cdot \frac{dx}{d\xi} | |||
+ \biggl[-2\alpha \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr] \cdot x \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
This should be compared with the first expression, [[#Development|above]], namely, | |||
<div align="center"> | |||
<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{\sin\xi }{\xi} \cdot \frac{d^2x}{d\xi^2} + 2 \biggl[\frac{ \sin\xi + \xi \cos \xi }{\xi^2}\biggr] \frac{dx}{d\xi} + | |||
\biggl[ \sigma^2 - 2\alpha \biggl( \frac{\sin\xi - \xi \cos \xi}{\xi^3} \biggr) \biggr]x \, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
and it matches! Q.E.D. | |||
==Motivated by Yabushita's Discovery== | |||
===Initial Exploration=== | |||
This subsection is being developed following our realization — see the [[SSC/Stability/InstabilityOnsetOverview#Overview:_Marginally_Unstable_Pressure-Truncated_Configurations|accompanying overview]] — that the eigenfunction ''is'' known analytically for marginally unstable, pressure-truncated configurations having <math>~3 \le n \le \infty</math>. Specifically, from the work of [http://adsabs.harvard.edu/abs/1975MNRAS.172..441Y Yabushita (1975)] we have the following, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="center" colspan="3"><font color="maroon"><b>Exact Solution to the Isothermal LAWE</b></font></td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\sigma_c^2 = 0</math> | |||
</td> | |||
<td align="center"> | |||
and | |||
</td> | |||
<td align="left"> | |||
<math>~x = 1 - \biggl( \frac{1}{\xi e^{-\psi}}\biggr) \frac{d\psi}{d\xi} \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
And from our own recent work, we have discovered the following, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="center" colspan="3"><font color="maroon"><b>Precise Solution to the Polytropic LAWE</b></font></td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
<math>~\sigma_c^2 = 0</math> | |||
</td> | |||
<td align="center"> | |||
and | |||
</td> | |||
<td align="left"> | |||
<math>~x_P \equiv \frac{3(n-1)}{2n}\biggl[1 + \biggl(\frac{n-3}{n-1}\biggr) \biggl( \frac{1}{\xi \theta^{n}}\biggr) \frac{d\theta}{d\xi}\biggr] </math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
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>. Using this polytropic displacement function as a guide, let's try for the case of <math>~n=1</math>, an expression of the form, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~A - B\biggl[ \biggl( \frac{1}{\xi \theta}\biggr) \frac{d\theta}{d\xi} \biggr]</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~A - B \biggl[ \biggl( \frac{1}{\sin\xi}\biggr) \frac{d}{d\xi} \biggl( \frac{\sin\xi}{\xi} \biggr)\biggr]</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~A - B \biggl( \frac{1}{\sin\xi}\biggr) \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~A + \frac{B}{\xi^2} \biggl( 1-\frac{\xi \cos\xi}{\sin\xi} \biggr) \, ,</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
in which case, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{dx}{d\xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~- B \biggl\{ | |||
\biggl( \frac{- \cos\xi}{\sin^2\xi}\biggr) \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] | |||
+\biggl( \frac{1}{\sin\xi}\biggr) \biggl[ -\frac{\sin\xi}{\xi} - \frac{\cos\xi}{\xi^2} - \frac{\cos\xi}{\xi^2} + \frac{2\sin\xi}{\xi^3} \biggr] | |||
\biggr\}</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~- \frac{B}{\xi^3} \biggl\{ | |||
\biggl( \frac{\cos\xi}{\sin^2\xi}\biggr) \biggl[ - \xi^2 \cos\xi + \xi \sin\xi \biggr] | |||
+\biggl[ 2 -\xi^2 - \frac{2\xi\cos\xi}{\sin\xi} \biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \frac{B}{\xi^3} \biggl\{ | |||
2 -\xi^2 - \frac{\xi\cos\xi}{\sin\xi} | |||
- \frac{\xi^2 \cos^2\xi}{\sin^2\xi} | |||
\biggr\} | |||
\, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
<table border="1" cellpadding="8" align="center" width="85%"><tr><td align="left"> | |||
What if, instead, we try the more generalized form, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~A + \frac{B}{(\lambda \xi)^2} \biggl[ 1-\frac{\lambda \xi \cos(\lambda \xi)}{\sin(\lambda \xi)} \biggr] \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Then we have, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{1}{\lambda B} \cdot \frac{dx}{d\xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \frac{1}{ (\lambda \xi)^3} \biggl\{ | |||
2 - (\lambda \xi)^2 - \frac{\lambda \xi\cos(\lambda \xi)}{\sin(\lambda \xi)} | |||
- \frac{(\lambda \xi)^2 \cos^2(\lambda\xi)}{\sin^2(\lambda\xi)} | |||
\biggr\} | |||
\, , | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
Probably this also means, | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{1}{\lambda^2 B} \cdot \frac{dx}{d\xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{1}{(\lambda\xi)^4} \biggl\{ 6 - 2 (\lambda \xi)^2 | |||
- \frac{2 \lambda \xi\cos(\lambda\xi)}{\sin(\lambda\xi)} | |||
- \frac{2(\lambda\xi)^2 \cos^2(\lambda\xi)}{\sin^2(\lambda\xi)} | |||
- \frac{2(\lambda\xi)^3 \cos(\lambda\xi)}{\sin(\lambda\xi)} | |||
- \frac{2(\lambda\xi)^3 \cos^3(\lambda\xi)}{\sin^3(\lambda\xi)} \biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</td></tr></table> | |||
Let's check against the [[SSC/Stability/Isothermal#Derivation_of_Polytropic_Displacement_Function|more general derivation]], which gives after recognizing that, <math>~B \leftrightarrow (3-n)/(n-1)</math>, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{dx}{d\xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl(\frac{3-n}{n-1}\biggr) \biggl\{ \frac{1}{\xi} | |||
+ \frac{n(\theta^')^2 }{\xi \theta^{n+1}} + \frac{3\theta^' }{\xi^2 \theta^{n}} \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{B}{\xi^3} \biggl\{ \xi^2 | |||
+ \xi^2 \biggl( \frac{\xi}{\sin\xi}\biggr)^2 \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr]^2 + \frac{3\xi^2}{\sin\xi} \biggl[ \frac{\cos\xi}{\xi} - \frac{\sin\xi}{\xi^2} \biggr] \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{B}{\xi^3} \biggl\{ \xi^2 | |||
+ 3\biggl[ \frac{\xi \cos\xi}{\sin\xi} - 1 \biggr] + \biggl[ \frac{\xi \cos\xi}{\sin\xi} - 1 \biggr]^2 \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{B}{\xi^3} \biggl\{ \xi^2 | |||
+ 3\biggl[ \frac{\xi \cos\xi}{\sin\xi} - 1 \biggr] | |||
+ \biggl[ \biggl(\frac{\xi \cos\xi}{\sin\xi}\biggr)^2 - 2\biggl(\frac{\xi \cos\xi}{\sin\xi}\biggr) + 1 \biggr] \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{B}{\xi^3} \biggl\{ \xi^2 | |||
+ \biggl[ \frac{\xi \cos\xi}{\sin\xi} - 2 \biggr] | |||
+ \biggl(\frac{\xi \cos\xi}{\sin\xi}\biggr)^2 \biggr\} | |||
\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
This matches the preceding, direct derivation. | |||
Also, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~\frac{d^2x}{d\xi^2}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{3B}{\xi^4} \biggl\{ | |||
\biggl( \frac{\cos\xi}{\sin^2\xi}\biggr) \biggl[ - \xi^2 \cos\xi + \xi \sin\xi \biggr] | |||
+\biggl[ 2 -\xi^2 - \frac{2\xi\cos\xi}{\sin\xi} \biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~- \frac{B}{\xi^3} \biggl\{ | |||
\biggl[- \frac{1}{\sin\xi} - \frac{2\cos^2\xi}{\sin^3\xi} \biggr] \biggl[ - \xi^2 \cos\xi + \xi \sin\xi \biggr] | |||
+ \biggl( \frac{\cos\xi}{\sin^2\xi}\biggr) \biggl[ - 2\xi \cos\xi + \sin\xi + \xi^2 \sin\xi + \xi \cos\xi \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+\biggl[ -2\xi - \frac{2\cos\xi}{\sin\xi} + \frac{2\xi\sin\xi}{\sin\xi} + \frac{2\xi\cos^2\xi}{\sin^2\xi}\biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{B}{\xi^4} \biggl\{ | |||
\biggl( \frac{3\cos\xi}{\sin^2\xi}\biggr) \biggl[ - \xi^2 \cos\xi + \xi \sin\xi \biggr] | |||
+\biggl[ 6 - 3\xi^2 - \frac{6\xi\cos\xi}{\sin\xi} \biggr] | |||
+ \biggl[\frac{1}{\sin\xi} + \frac{2\cos^2\xi}{\sin^3\xi} \biggr] \biggl[ - \xi^3 \cos\xi + \xi^2 \sin\xi \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ \biggl( \frac{\cos\xi}{\sin^2\xi}\biggr) \biggl[ 2\xi^2 \cos\xi - \xi \sin\xi - \xi^3 \sin\xi - \xi^2 \cos\xi \biggr] | |||
+\biggl[ 2\xi^2 + \frac{2\xi \cos\xi}{\sin\xi} - \frac{2\xi^2\sin\xi}{\sin\xi} - \frac{2\xi^2\cos^2\xi}{\sin^2\xi}\biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{B}{\xi^4} \biggl\{ | |||
\biggl[ - \frac{3\xi^2 \cos^2\xi}{\sin^2\xi} + \frac{3\xi \cos\xi}{\sin\xi} \biggr] | |||
+\biggl[ 6 - 3\xi^2 - \frac{6\xi\cos\xi}{\sin\xi} \biggr] | |||
+ \biggl[ - \frac{\xi^3 \cos\xi}{\sin\xi} + \xi^2 \biggr] | |||
+ \biggl[ - \frac{2\xi^3 \cos^3\xi}{\sin^3\xi} + \frac{2\xi^2 \cos^2\xi}{\sin^2\xi} \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ \biggl[ \frac{2\xi^2 \cos^2\xi}{\sin^2\xi} - \frac{\xi \cos\xi}{\sin\xi} - \frac{\xi^3 \cos\xi}{\sin\xi} - \frac{\xi^2 \cos^2\xi}{\sin^2\xi} \biggr] | |||
+\biggl[ 2\xi^2 + \frac{2\xi \cos\xi}{\sin\xi} - \frac{2\xi^2\sin\xi}{\sin\xi} - \frac{2\xi^2\cos^2\xi}{\sin^2\xi}\biggr] | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{B}{\xi^4} \biggl\{ 6 - 2\xi^2 | |||
- \frac{2\xi\cos\xi}{\sin\xi} | |||
- \frac{2\xi^2 \cos^2\xi}{\sin^2\xi} | |||
- \frac{2\xi^3 \cos\xi}{\sin\xi} | |||
- \frac{2\xi^3 \cos^3\xi}{\sin^3\xi} \biggr\} \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Let's also check this against the [[SSC/Stability/Isothermal#Derivation_of_Polytropic_Displacement_Function|more general derivation]], which gives after again recognizing that, <math>~B \leftrightarrow (3-n)/(n-1)</math>, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>\frac{d^2 x}{d\xi^2}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl(\frac{n-3}{n-1}\biggr) \biggl\{ \frac{4}{\xi^2} + \frac{2n(\theta^')}{\xi \theta} | |||
+ \frac{12\theta^' }{\xi^3 \theta^{n}}+ \frac{8n(\theta^')^2}{\xi^2 \theta^{n+1}} | |||
+ (n+1) \frac{n(\theta^')^3 }{\xi \theta^{n+2}} \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-B | |||
\biggl\{ \frac{4}{\xi^2} + \frac{2}{\xi \theta} \biggl[ \frac{\sin\xi}{\xi^2}\biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)\biggr] | |||
+ \frac{12}{\xi^3 \theta}\biggl[ \frac{\sin\xi}{\xi^2}\biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)\biggr] | |||
+ \frac{8 }{\xi^2 \theta^{2}} \biggl[ \frac{\sin\xi}{\xi^2}\biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)\biggr]^2 | |||
+ \frac{2 }{\xi \theta^{3}} \biggl[ \frac{\sin\xi}{\xi^2}\biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)\biggr]^3\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-\frac{B}{\xi^4} | |||
\biggl\{ 4\xi^2 + 2\xi^2 \biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr) | |||
+ 12\biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr) | |||
+ 8 \biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)^2 | |||
+ 2\biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)^3\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-\frac{2B}{\xi^4} | |||
\biggl\{ 2\xi^2 + \frac{\xi^3\cos\xi}{\sin\xi} - \xi^2+ \frac{6\xi\cos\xi}{\sin\xi} - 6 | |||
+ \frac{4\xi^2\cos^2\xi}{\sin^2\xi} - \frac{8\xi\cos\xi}{\sin\xi} + 4 | |||
+ \biggl(\frac{\xi^2\cos^2\xi}{\sin^2\xi} - \frac{2\xi\cos\xi}{\sin\xi} + 1\biggr) \biggl(\frac{\xi\cos\xi}{\sin\xi} - 1\biggr)\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-\frac{2B}{\xi^4} | |||
\biggl\{-2 + \xi^2 - \frac{2\xi\cos\xi}{\sin\xi} + \frac{\xi^3\cos\xi}{\sin\xi} | |||
+ \frac{4\xi^2\cos^2\xi}{\sin^2\xi} | |||
- \biggl(\frac{\xi^2\cos^2\xi}{\sin^2\xi} - \frac{2\xi\cos\xi}{\sin\xi} + 1\biggr) | |||
+ \frac{\xi^3\cos^3\xi}{\sin^3\xi} - \frac{2\xi^2\cos^2\xi}{\sin^2\xi} + \frac{\xi\cos\xi}{\sin\xi} | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~-\frac{2B}{\xi^4} | |||
\biggl\{-3 + \xi^2 + \frac{\xi\cos\xi}{\sin\xi} + \frac{\xi^3\cos\xi}{\sin\xi} | |||
+ \frac{\xi^2\cos^2\xi}{\sin^2\xi} | |||
+ \frac{\xi^3\cos^3\xi}{\sin^3\xi} | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{B}{\xi^4} | |||
\biggl\{6 -2 \xi^2 - \frac{2\xi\cos\xi}{\sin\xi} | |||
- \frac{2\xi^2\cos^2\xi}{\sin^2\xi} | |||
- \frac{2\xi^3\cos\xi}{\sin\xi} | |||
- \frac{2\xi^3\cos^3\xi}{\sin^3\xi} | |||
\biggr\} | |||
\, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
A cross-check with the first attempt to derive this second derivative expression initially unveiled a couple of coefficient errors. These have now been corrected and both expressions agree. | |||
===Succinct Demonstration=== | |||
Given that, for <math>~n=1</math>, we should set <math>~\gamma_\mathrm{g} = (n+1)/n = 2 \Rightarrow \alpha = (3-4/\gamma_\mathrm{g}) = +1</math>, and, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~Q \equiv - \frac{d\ln\theta}{d\ln\xi}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \frac{\xi^2}{\sin\xi} \cdot \frac{d}{d\xi}\biggl[ \frac{\sin\xi}{\xi}\biggr] | |||
= | |||
1 - \xi \cot\xi \, . | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
If we then employ the displacement function, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~x</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~A + \frac{B}{\xi^2} \biggl[ 1 - \xi \cot\xi \biggr] \, ,</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
the LAWE becomes, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
LAWE | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\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> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{d^2x}{d\xi^2} + \biggl[4 - 2Q \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} + | |||
\biggl[ \biggl(\frac{\sigma_c^2}{6 } \biggr) \frac{\xi^3}{\sin\xi} | |||
- 2Q\biggr] \frac{ x}{\xi^2} </math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{d^2x}{d\xi^2} | |||
+ \biggl[2 + \frac{2\xi\cos\xi}{\sin\xi} \biggr] \frac{1}{\xi} \cdot \frac{dx}{d\xi} | |||
+ \biggl[- 2 + \frac{2\xi\cos\xi}{\sin\xi} \biggr] \frac{ x}{\xi^2} | |||
+ \biggl[ \biggl(\frac{\sigma_c^2}{6 } \biggr) \frac{\xi}{\sin\xi} \biggr] x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{2B}{\xi^4} | |||
\biggl\{3 - \xi^2 - \frac{\xi\cos\xi}{\sin\xi} | |||
- \biggl(\frac{\xi\cos\xi}{\sin\xi} \biggr)^2 | |||
- \frac{\xi^3\cos\xi}{\sin\xi} | |||
- \biggl( \frac{\xi\cos\xi}{\sin\xi} \biggr)^3 | |||
\biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ \frac{2B}{\xi^4}\biggl[1 + \frac{\xi\cos\xi}{\sin\xi} \biggr] \biggl\{ \xi^2 | |||
- 2 + \frac{\xi \cos\xi}{\sin\xi} | |||
+ \biggl(\frac{\xi \cos\xi}{\sin\xi}\biggr)^2 \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ \biggl[- 2 + \frac{2\xi\cos\xi}{\sin\xi} \biggr] \biggl[ \frac{A}{\xi^2} + \frac{B}{\xi^4} \biggl( 1-\frac{\xi \cos\xi}{\sin\xi} \biggr)\biggr] | |||
+ \biggl[ \biggl(\frac{\sigma_c^2}{6 } \biggr) \frac{\xi}{\sin\xi} \biggr] x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{2B}{\xi^4} | |||
\biggl\{3 - \xi^2 - \frac{\xi\cos\xi}{\sin\xi} | |||
- \biggl(\frac{\xi\cos\xi}{\sin\xi} \biggr)^2 | |||
- \frac{\xi^3\cos\xi}{\sin\xi} | |||
- \biggl( \frac{\xi\cos\xi}{\sin\xi} \biggr)^3 | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
+ \xi^2\biggl( \frac{\xi\cos\xi}{\sin\xi} \biggr) - 2\biggl( \frac{\xi\cos\xi}{\sin\xi} \biggr) + \biggl(\frac{\xi \cos\xi}{\sin\xi} \biggr)^2 | |||
+ \biggl(\frac{\xi \cos\xi}{\sin\xi}\biggr)^3 | |||
+ \xi^2 - 2 + \frac{\xi \cos\xi}{\sin\xi} | |||
+ \biggl(\frac{\xi \cos\xi}{\sin\xi}\biggr)^2 \biggr\} | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
| |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
- \frac{2B}{\xi^4} \biggl[ 1 - \frac{2\xi\cos\xi}{\sin\xi} + \biggl(\frac{\xi \cos\xi}{\sin\xi} \biggr)^2 \biggr] | |||
+ \frac{2A}{\xi^2}\biggl[\frac{\xi\cos\xi}{\sin\xi} -1\biggr] | |||
+ \biggl[ \biggl(\frac{\sigma_c^2}{6 } \biggr) \frac{\xi}{\sin\xi} \biggr] x | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{2A}{\xi^2}\biggl[\frac{\xi\cos\xi}{\sin\xi} -1\biggr] | |||
+ \biggl[ \biggl(\frac{\sigma_c^2}{6 } \biggr) \frac{\xi}{\sin\xi} \biggr] x | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
Pretty amazing degree of cancelation! So the above-hypothesized displacement function ''does'' satisfy the <math>~n=1</math>, polytropic LAWE — for any value of the coefficient, <math>~B</math> — if we set <math>~A = 0</math> and <math>~\sigma_c^2=0</math>. If we set <math>~B = 3</math>, the function will be normalized such that it goes to unity at the center. In summary, then, we have, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
<math>~x_P\biggr|_{n=1}</math> | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~\frac{3}{\xi^2} \biggl[ 1 - \xi \cot\xi \biggr] \, .</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
<!-- | |||
Let's play with this a bit more to see if we can uncover a displacement function that works for nonzero values of <math>~\omega_c^2</math>. Leaving both <math>~A</math> and <math>~B</math> unspecified for the time being, we have, | |||
<div align="center"> | |||
<table border="0" cellpadding="5" align="center"> | |||
<tr> | |||
<td align="right"> | |||
LAWE | |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\frac{2A}{\xi^2}\biggl[\frac{\xi\cos\xi}{\sin\xi} -1\biggr] | |||
+ \biggl(\frac{\sigma_c^2}{6 } \biggr) \frac{\xi}{\sin\xi} \biggl[ A + \frac{B}{\xi^2} \biggl( 1-\frac{\xi \cos\xi}{\sin\xi} \biggr) \biggr] | |||
</math> | |||
</td> | |||
</tr> | |||
<tr> | |||
<td align="right"> | |||
| |||
</td> | |||
<td align="center"> | |||
<math>~=</math> | |||
</td> | |||
<td align="left"> | |||
<math>~ | |||
\biggl(\frac{A\sigma_c^2}{6 } \biggr) \frac{\xi}{\sin\xi} | |||
+ \biggl[ \biggl(\frac{B\sigma_c^2}{6 } \biggr) \frac{1}{\xi\sin\xi} | |||
- \frac{2A}{\xi^2} \biggr] \biggl( 1 - \frac{\xi\cos\xi}{\sin\xi} \biggr) | |||
</math> | |||
</td> | |||
</tr> | |||
</table> | |||
</div> | |||
--> | |||
Revision as of 22:08, 20 January 2024
Radial Oscillations of n = 1 Polytropic Spheres (Pt 2)
Part I: Search for Analytic Solutions |
Part II: New Ideas |
III: What About Bipolytropes? |
New Idea Involving Logarithmic Derivatives
Simplistic Layout
Let's begin, again, with the relevant LAWE, as provided above. After dividing through by , we have,
where,
|
|
|
|
|
|
|
|
Now, in addition to recognizing that,
|
|
|
|
in a separate context, we showed that, quite generally,
|
|
|
|
Hence, if we assume that the eigenfunction is a power-law of , that is, assume that,
then the logarithmic derivative of is a constant, namely,
and the two key derivative terms will be,
|
|
and |
|
In this case, the LAWE is no longer a differential equation but, instead, takes the form,
|
|
|
|
|
|
|
|
|
|
|
|
Now, the cosine term will go to zero if ; and the sine term will go to zero if,
|
|
|
|
|
|
|
|
If these two — rather strange — conditions are met, then we have a marginally unstable configuration because, . This, in and of itself, is not very physically interesting. However, it may give us a clue regarding how to more generally search for a physically reasonable radial eigenfunction.
More general Assumption
Try,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Another Viewpoint
Development
Multiplying through the above LAWE by gives,
|
|
|
|
Notice that,
|
|
|
|
|
|
|
|
And, hence,
|
|
|
|
|
|
|
|
|
|
|
|
So, we can write,
|
|
|
|
|
|
|
|
|
|
|
|
This means that we can rewrite the LAWE as,
|
|
|
|
|
|
|
|
We recognize, also, that,
|
|
|
|
|
|
|
|
So the LAWE becomes,
|
|
|
|
|
|
|
|
|
|
|
|
where we have introduced the new, modified eigenfunction,
Alternatively, the LAWE may be written as,
|
|
|
|
or,
|
|
|
|
|
|
|
|
Now, if we adopt the homentropic convention that arises from setting, , then for our polytropic configuration, we should set, and, hence, . This will mean that the lat term in this LAWE naturally goes to zero. Hence, we have,
|
|
|
|
or,
|
|
|
|
or,
|
|
|
|
Does this help?
Check for Mistakes
Given the definition of , its first derivative is,
|
|
|
|
and its second derivative is,
|
|
|
|
|
|
|
|
|
|
|
|
Hence, the "upsilon" LAWE becomes,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
This should be compared with the first expression, above, namely,
|
|
|
|
and it matches! Q.E.D.
Motivated by Yabushita's Discovery
Initial Exploration
This subsection is being developed following our realization — see the accompanying overview — that the eigenfunction is known analytically for marginally unstable, pressure-truncated configurations having . Specifically, from the work of Yabushita (1975) we have the following,
| Exact Solution to the Isothermal LAWE | ||
|
|
and |
|
And from our own recent work, we have discovered the following,
| Precise Solution to the Polytropic LAWE | ||
|
|
and |
|
if the adiabatic exponent is assigned the value, , in which case the parameter, . Using this polytropic displacement function as a guide, let's try for the case of , an expression of the form,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
in which case,
|
|
|
|
|
|
|
|
|
|
|
|
|
What if, instead, we try the more generalized form,
Then we have,
Probably this also means,
|
Let's check against the more general derivation, which gives after recognizing that, ,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
This matches the preceding, direct derivation.
Also,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Let's also check this against the more general derivation, which gives after again recognizing that, ,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A cross-check with the first attempt to derive this second derivative expression initially unveiled a couple of coefficient errors. These have now been corrected and both expressions agree.
Succinct Demonstration
Given that, for , we should set , and,
|
|
|
|
If we then employ the displacement function,
|
|
|
|
the LAWE becomes,
|
LAWE |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pretty amazing degree of cancelation! So the above-hypothesized displacement function does satisfy the , polytropic LAWE — for any value of the coefficient, — if we set and . If we set , the function will be normalized such that it goes to unity at the center. In summary, then, we have,
|
|
|
|
See Also
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |