Editing
SSC/Stability/n1PolytropeLAWE/Pt2
Jump to navigation
Jump to search
Warning:
You are not logged in. Your IP address will be publicly visible if you make any edits. If you
log in
or
create an account
, your edits will be attributed to your username, along with other benefits.
Anti-spam check. Do
not
fill this in!
__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="25%"><br />[[SSC/Stability/n1PolytropeLAWE|Part I: Search for Analytic Solutions]]<br /> </td> <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/n1PolytropeLAWE/Pt2|Part II: New Ideas]]<br /> </td> <td align="center" rowspan="1" bgcolor="lightblue" width="25%"><br />[[SSC/Stability/n1PolytropeLAWE/Pt3|Part III: What About Bipolytropes?]]<br /> </td> <td align="center" rowspan="1" bgcolor="lightblue"><br />[[SSC/Stability/n1PolytropeLAWE/Pt4|Part IV: Most General Structural Solution]]<br /> </td> </tr> </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> --> =See Also= {{ SGFfooter }}
Summary:
Please note that all contributions to JETohlineWiki may be edited, altered, or removed by other contributors. If you do not want your writing to be edited mercilessly, then do not submit it here.
You are also promising us that you wrote this yourself, or copied it from a public domain or similar free resource (see
JETohlineWiki:Copyrights
for details).
Do not submit copyrighted work without permission!
Cancel
Editing help
(opens in new window)
Template used on this page:
Template:SGFfooter
(
edit
)
Navigation menu
Personal tools
Not logged in
Talk
Contributions
Log in
Namespaces
Page
Discussion
English
Views
Read
Edit
View history
More
Search
Navigation
Main page
Tiled Menu
Table of Contents
Old (VisTrails) Cover
Appendices
Variables & Parameters
Key Equations
Special Functions
Permissions
Formats
References
lsuPhys
Ramblings
Uploaded Images
Originals
Recent changes
Random page
Help about MediaWiki
Tools
What links here
Related changes
Special pages
Page information