SSC/Stability/n1PolytropeLAWE/Pt2: Difference between revisions

From JETohlineWiki
Jump to navigation Jump to search
Joel2 (talk | contribs)
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..."
 
Joel2 (talk | contribs)
 
(One intermediate revision by the same user not shown)
Line 3: Line 3:
<table border="1" align="center" width="100%" colspan="8">
<table border="1" align="center" width="100%" colspan="8">
<tr>
<tr>
   <td align="center" rowspan="1" bgcolor="lightblue" width="33%"><br />[[SSC/Stability/n1PolytropeLAWE|Part I: &nbsp; Search for Analytic Solutions]]<br />&nbsp;</td>
   <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/n1PolytropeLAWE|Part I: &nbsp; Search for Analytic Solutions]]<br />&nbsp;</td>
   <td align="center" rowspan="1" bgcolor="lightblue" width="33%"><br />[[SSC/Stability/n1PolytropeLAWE/Pt2|Part II:&nbsp; New Ideas]]<br />&nbsp;</td>
   <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/n1PolytropeLAWE/Pt2|Part II:&nbsp; New Ideas]]<br />&nbsp;</td>
   <td align="center" rowspan="1" bgcolor="lightblue"><br />[[SSC/Stability/n1PolytropeLAWE/Pt3|III:&nbsp; What About Bipolytropes?]]<br />&nbsp;</td>
   <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/n1PolytropeLAWE/Pt3|Part III:&nbsp; What About Bipolytropes?]]<br />&nbsp;</td>
  <td align="center" rowspan="1" bgcolor="lightblue"><br />[[SSC/Stability/n1PolytropeLAWE/Pt4|Part IV:&nbsp; Most General Structural Solution]]<br />&nbsp;</td>
</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">
&nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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 &#8212; rather strange &#8212; 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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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 &#8212; see the [[SSC/Stability/InstabilityOnsetOverview#Overview:_Marginally_Unstable_Pressure-Truncated_Configurations|accompanying overview]] &#8212; 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">
&nbsp;and &nbsp;
  </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">
&nbsp; &nbsp; &nbsp; and &nbsp; &nbsp; &nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </td>
  <td align="center">
&nbsp;
  </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">
&nbsp;
  </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 &#8212; for any value of the coefficient, <math>~B</math> &#8212; 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">
&nbsp;
  </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>
-->


=See Also=
=See Also=


{{ SGFfooter }}
{{ SGFfooter }}

Latest revision as of 13:52, 15 July 2025

Radial Oscillations of n = 1 Polytropic Spheres (Pt 2)[edit]


Part I:   Search for Analytic Solutions
 

Part II:  New Ideas
 

Part III:  What About Bipolytropes?
 

Part IV:  Most General Structural Solution
 


New Idea Involving Logarithmic Derivatives[edit]

Simplistic Layout[edit]

Let's begin, again, with the relevant LAWE, as provided above. After dividing through by x, we have,

(sinξ)ξ2xd2xdξ2+2[sinξ+ξcosξ]ξxdxdξ+[σ2ξ32α(sinξξcosξ)]=0,


where,

σ2

ω22πGρcγg,

α

34γg.

Now, in addition to recognizing that,

ξxdxdξ

=

dlnxdlnξ,

in a separate context, we showed that, quite generally,

ξ2xd2xdξ2

=

ddlnξ[dlnxdlnξ][1dlnxdlnξ]dlnxdlnξ.

Hence, if we assume that the eigenfunction is a power-law of ξ, that is, assume that,

x=a0ξc0,

then the logarithmic derivative of x is a constant, namely,

dlnxdlnξ=c0,

and the two key derivative terms will be,

ξxdxdξ=c0,

      and      

ξ2xd2xdξ2=c0(c01).

In this case, the LAWE is no longer a differential equation but, instead, takes the form,

σ2ξ3

=

c0(c01)sinξ+2c0[sinξ+ξcosξ]2α(sinξξcosξ)

 

=

sinξ[c0(c01)+2c02α]+ξcosξ[2(c0+α)]

 

=

sinξ[c02+c02α]+ξcosξ[2(c0+α)].

Now, the cosine term will go to zero if c0=α; and the sine term will go to zero if,

α

=

3

γg

=

.

If these two — rather strange — conditions are met, then we have a marginally unstable configuration because, σ2=0. 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[edit]

Try,

x

=

ξc0[a0+b0sinξ+d0ξcosξ]

dxdξ

=

ξc0ddξ[a0+b0sinξ+d0ξcosξ]+c0ξc01[a0+b0sinξ+d0ξcosξ]

 

=

ξc0[b0cosξd0ξsinξ+d0cosξ]+c0ξc01[a0+b0sinξ+d0ξcosξ]

dlnxdlnξ

=

ξ[b0cosξd0ξsinξ+d0cosξ][a0+b0sinξ+d0ξcosξ]1+c0

 

=

[(b0+d0)ξcosξd0ξ2sinξ][a0+b0sinξ+d0ξcosξ]1+c0


Another Viewpoint[edit]

Development[edit]

Multiplying through the above LAWE by (xξ3) gives,

0

=

sinξξd2xdξ2+2[sinξ+ξcosξξ2]dxdξ+[σ22α(sinξξcosξξ3)]x

Notice that,

ddξ[sinξξ]

=

sinξξ2+cosξξ

 

=

[ξcosξsinξξ2].

And, hence,

d2dξ2[sinξξ]

=

ddξ[cosξξsinξξ2]

 

=

cosξξ2sinξξ+2sinξξ3cosξξ2

 

=

sinξξ+2[sinξξcosξξ3].

So, we can write,

d2dξ2{(sinξξ)x}

=

ddξ{(sinξξ)dxdξ+xddξ[(sinξξ)]}

 

=

sinξξd2xdξ2+2dxdξ[ddξ(sinξξ)]+xd2dξ2(sinξξ)

 

=

sinξξd2xdξ2+2dxdξ[ξcosξsinξξ2]+x{sinξξ+2[sinξξcosξξ3]}.

This means that we can rewrite the LAWE as,

0

=

d2dξ2{(sinξξ)x}2dxdξ[ξcosξsinξξ2]x{sinξξ+2[sinξξcosξξ3]}+2[sinξ+ξcosξξ2]dxdξ+[σ22α(sinξξcosξξ3)]x

 

=

d2dξ2{(sinξξ)x}+4[sinξξ2]dxdξ+{sinξξ+σ22(1+α)(sinξξcosξξ3)}x.

We recognize, also, that,

1ξddξ[(sinξξ)x]

=

[ξcosξsinξξ3]x+(sinξξ2)dxdξ.

4(sinξξ2)dxdξ

=

4ξddξ[(sinξξ)x]+4[sinξξcosξξ3]x.

So the LAWE becomes,

0

=

d2dξ2{(sinξξ)x}+4ξddξ[(sinξξ)x]+4[sinξξcosξξ3]x+{sinξξ+σ22(1+α)(sinξξcosξξ3)}x

 

=

d2dξ2{(sinξξ)x}+4ξddξ[(sinξξ)x]+{sinξξ+σ2+[42(1+α)](sinξξcosξξ3)}x

 

=

d2Υdξ2+4ξdΥdξ+Υ+[σ2+2(1α)(sinξξcosξξ3)]x,

where we have introduced the new, modified eigenfunction,

Υ(sinξξ)x.

Alternatively, the LAWE may be written as,

0

=

d2Υdξ2+4ξdΥdξ+[σ2+2(1α)(sinξξcosξξ3)+sinξξ]x;

or,

0

=

ξ2Υd2Υdξ2+4ξΥdΥdξ+[σ2+2(1α)(sinξξcosξξ3)+sinξξ]ξ3sinξ

 

=

ξ2Υd2Υdξ2+4ξΥdΥdξ+[σ2(ξ3sinξ)+2(1α)(1ξcotξ)+ξ2]


Now, if we adopt the homentropic convention that arises from setting, γ=(n+1)/n, then for our n=1 polytropic configuration, we should set, γ=2 and, hence, α=1. This will mean that the lat term in this LAWE naturally goes to zero. Hence, we have,

σ2x

=

d2Υdξ2+4ξdΥdξ+Υ;

or,

0

=

d2Υdξ2+4ξdΥdξ+[1+σ2(ξsinξ)]Υ;

or,

0

=

ξ2Υd2Υdξ2+4ξΥdΥdξ+[ξ2+σ2(ξ3sinξ)].

Does this help?

Check for Mistakes[edit]

Given the definition of Υ, its first derivative is,

dΥdξ

=

(sinξξ)dxdξ+x[cosξξsinξξ2],

and its second derivative is,

d2Υdξ2

=

ddξ{(sinξξ)dxdξ+x[cosξξsinξξ2]}

 

=

(sinξξ)d2xdξ2+2[cosξξsinξξ2]dxdξ+xddξ[cosξξsinξξ2]

 

=

(sinξξ)d2xdξ2+2[cosξξsinξξ2]dxdξ+x[sinξξ2cosξξ2+2sinξξ3]

Hence, the "upsilon" LAWE becomes,

σ2x

=

d2Υdξ2+4ξdΥdξ+Υ+[2(1α)(sinξξcosξξ3)]x

 

=

(sinξξ)d2xdξ2+2[cosξξsinξξ2]dxdξ+x[sinξξ2cosξξ2+2sinξξ3]+4ξ{(sinξξ)dxdξ+x[cosξξsinξξ2]}+[sinξξ+2(1α)(sinξξcosξξ3)]x

 

=

(sinξξ)d2xdξ2+{(4sinξξ2)+2[cosξξsinξξ2]}dxdξ+[sinξξ2cosξξ2+2sinξξ3+4cosξξ24sinξξ3+sinξξ+2(1α)(sinξξcosξξ3)]x

 

=

(sinξξ)d2xdξ2+[2cosξξ+2sinξξ2]dxdξ+[2(sinξξcosξξ3)+(22α)(sinξξcosξξ3)]x

 

=

(sinξξ)d2xdξ2+2[sinξξ2+cosξξ]dxdξ+[2α(sinξξcosξξ3)]x.

This should be compared with the first expression, above, namely,

0

=

sinξξd2xdξ2+2[sinξ+ξcosξξ2]dxdξ+[σ22α(sinξξcosξξ3)]x,

and it matches! Q.E.D.

Motivated by Yabushita's Discovery[edit]

Initial Exploration[edit]

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 3n. Specifically, from the work of Yabushita (1975) we have the following,

Exact Solution to the Isothermal LAWE

σc2=0

 and  

x=1(1ξeψ)dψdξ.

And from our own recent work, we have discovered the following,

Precise Solution to the Polytropic LAWE

σc2=0

      and      

xP3(n1)2n[1+(n3n1)(1ξθn)dθdξ]

if the adiabatic exponent is assigned the value, γg=(n+1)/n, in which case the parameter, α=(3n)/(n+1). Using this polytropic displacement function as a guide, let's try for the case of n=1, an expression of the form,

x

=

AB[(1ξθ)dθdξ]

 

=

AB[(1sinξ)ddξ(sinξξ)]

 

=

AB(1sinξ)[cosξξsinξξ2]

 

=

A+Bξ2(1ξcosξsinξ),

in which case,

dxdξ

=

B{(cosξsin2ξ)[cosξξsinξξ2]+(1sinξ)[sinξξcosξξ2cosξξ2+2sinξξ3]}

 

=

Bξ3{(cosξsin2ξ)[ξ2cosξ+ξsinξ]+[2ξ22ξcosξsinξ]}

 

=

Bξ3{2ξ2ξcosξsinξξ2cos2ξsin2ξ},


What if, instead, we try the more generalized form,

x

=

A+B(λξ)2[1λξcos(λξ)sin(λξ)].

Then we have,

1λBdxdξ

=

1(λξ)3{2(λξ)2λξcos(λξ)sin(λξ)(λξ)2cos2(λξ)sin2(λξ)},

Probably this also means,

1λ2Bdxdξ

=

1(λξ)4{62(λξ)22λξcos(λξ)sin(λξ)2(λξ)2cos2(λξ)sin2(λξ)2(λξ)3cos(λξ)sin(λξ)2(λξ)3cos3(λξ)sin3(λξ)}.



Let's check against the more general derivation, which gives after recognizing that, B(3n)/(n1),

dxdξ

=

(3nn1){1ξ+n(θ')2ξθn+1+3θ'ξ2θn}

 

=

Bξ3{ξ2+ξ2(ξsinξ)2[cosξξsinξξ2]2+3ξ2sinξ[cosξξsinξξ2]}

 

=

Bξ3{ξ2+3[ξcosξsinξ1]+[ξcosξsinξ1]2}

 

=

Bξ3{ξ2+3[ξcosξsinξ1]+[(ξcosξsinξ)22(ξcosξsinξ)+1]}

 

=

Bξ3{ξ2+[ξcosξsinξ2]+(ξcosξsinξ)2}.

This matches the preceding, direct derivation.

Also,

d2xdξ2

=

3Bξ4{(cosξsin2ξ)[ξ2cosξ+ξsinξ]+[2ξ22ξcosξsinξ]}

 

 

Bξ3{[1sinξ2cos2ξsin3ξ][ξ2cosξ+ξsinξ]+(cosξsin2ξ)[2ξcosξ+sinξ+ξ2sinξ+ξcosξ]

 

 

+[2ξ2cosξsinξ+2ξsinξsinξ+2ξcos2ξsin2ξ]}

 

=

Bξ4{(3cosξsin2ξ)[ξ2cosξ+ξsinξ]+[63ξ26ξcosξsinξ]+[1sinξ+2cos2ξsin3ξ][ξ3cosξ+ξ2sinξ]

 

 

+(cosξsin2ξ)[2ξ2cosξξsinξξ3sinξξ2cosξ]+[2ξ2+2ξcosξsinξ2ξ2sinξsinξ2ξ2cos2ξsin2ξ]}

 

=

Bξ4{[3ξ2cos2ξsin2ξ+3ξcosξsinξ]+[63ξ26ξcosξsinξ]+[ξ3cosξsinξ+ξ2]+[2ξ3cos3ξsin3ξ+2ξ2cos2ξsin2ξ]

 

 

+[2ξ2cos2ξsin2ξξcosξsinξξ3cosξsinξξ2cos2ξsin2ξ]+[2ξ2+2ξcosξsinξ2ξ2sinξsinξ2ξ2cos2ξsin2ξ]}

 

=

Bξ4{62ξ22ξcosξsinξ2ξ2cos2ξsin2ξ2ξ3cosξsinξ2ξ3cos3ξsin3ξ}.

Let's also check this against the more general derivation, which gives after again recognizing that, B(3n)/(n1),

d2xdξ2

=

(n3n1){4ξ2+2n(θ')ξθ+12θ'ξ3θn+8n(θ')2ξ2θn+1+(n+1)n(θ')3ξθn+2}

 

=

B{4ξ2+2ξθ[sinξξ2(ξcosξsinξ1)]+12ξ3θ[sinξξ2(ξcosξsinξ1)]+8ξ2θ2[sinξξ2(ξcosξsinξ1)]2+2ξθ3[sinξξ2(ξcosξsinξ1)]3}

 

=

Bξ4{4ξ2+2ξ2(ξcosξsinξ1)+12(ξcosξsinξ1)+8(ξcosξsinξ1)2+2(ξcosξsinξ1)3}

 

=

2Bξ4{2ξ2+ξ3cosξsinξξ2+6ξcosξsinξ6+4ξ2cos2ξsin2ξ8ξcosξsinξ+4+(ξ2cos2ξsin2ξ2ξcosξsinξ+1)(ξcosξsinξ1)}

 

=

2Bξ4{2+ξ22ξcosξsinξ+ξ3cosξsinξ+4ξ2cos2ξsin2ξ(ξ2cos2ξsin2ξ2ξcosξsinξ+1)+ξ3cos3ξsin3ξ2ξ2cos2ξsin2ξ+ξcosξsinξ}

 

=

2Bξ4{3+ξ2+ξcosξsinξ+ξ3cosξsinξ+ξ2cos2ξsin2ξ+ξ3cos3ξsin3ξ}

 

=

Bξ4{62ξ22ξcosξsinξ2ξ2cos2ξsin2ξ2ξ3cosξsinξ2ξ3cos3ξsin3ξ}.

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[edit]

Given that, for n=1, we should set γg=(n+1)/n=2α=(34/γg)=+1, and,

Qdlnθdlnξ

=

ξ2sinξddξ[sinξξ]=1ξcotξ.

If we then employ the displacement function,

x

=

A+Bξ2[1ξcotξ],

the LAWE becomes,

LAWE

=

d2xdξ2+[4(n+1)Q]1ξdxdξ+(n+1)[(σc26γg)ξ2θαQ]xξ2

 

=

d2xdξ2+[42Q]1ξdxdξ+[(σc26)ξ3sinξ2Q]xξ2

 

=

d2xdξ2+[2+2ξcosξsinξ]1ξdxdξ+[2+2ξcosξsinξ]xξ2+[(σc26)ξsinξ]x

 

=

2Bξ4{3ξ2ξcosξsinξ(ξcosξsinξ)2ξ3cosξsinξ(ξcosξsinξ)3}

 

 

+2Bξ4[1+ξcosξsinξ]{ξ22+ξcosξsinξ+(ξcosξsinξ)2}

 

 

+[2+2ξcosξsinξ][Aξ2+Bξ4(1ξcosξsinξ)]+[(σc26)ξsinξ]x

 

=

2Bξ4{3ξ2ξcosξsinξ(ξcosξsinξ)2ξ3cosξsinξ(ξcosξsinξ)3

 

 

+ξ2(ξcosξsinξ)2(ξcosξsinξ)+(ξcosξsinξ)2+(ξcosξsinξ)3+ξ22+ξcosξsinξ+(ξcosξsinξ)2}

 

 

2Bξ4[12ξcosξsinξ+(ξcosξsinξ)2]+2Aξ2[ξcosξsinξ1]+[(σc26)ξsinξ]x

 

=

2Aξ2[ξcosξsinξ1]+[(σc26)ξsinξ]x

Pretty amazing degree of cancelation! So the above-hypothesized displacement function does satisfy the n=1, polytropic LAWE — for any value of the coefficient, B — if we set A=0 and σc2=0. If we set B=3, the function will be normalized such that it goes to unity at the center. In summary, then, we have,

xP|n=1

=

3ξ2[1ξcotξ].


See Also[edit]

Tiled Menu

Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS |