Editing
Appendix/Ramblings/PowerSeriesExpressions
(section)
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!
===Displacement Function for Polytropic LAWE=== The [[SSC/Stability/Polytropes#Adiabatic_.28Polytropic.29_Wave_Equation|LAWE for polytropic spheres]] 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^2x}{d\xi^2} + \biggl[\frac{4}{\xi} - \frac{(n+1)}{\theta} \biggl(- \frac{d\theta}{d\xi}\biggr)\biggr] \frac{dx}{d\xi} + \frac{(n+1)}{\theta}\biggl[\frac{\sigma_c^2}{6\gamma } - \frac{\alpha}{\xi} \biggl(- \frac{d\theta}{d\xi}\biggr)\biggr] x </math> </td> </tr> <tr> <td align="right"> </td> <td align="center"> <math>~=</math> </td> <td align="left"> <math>~\theta \frac{d^2x}{d\xi^2} + \biggl[4\theta - (n+1)\xi \biggl(- \frac{d\theta}{d\xi}\biggr)\biggr] \frac{1}{\xi}\frac{dx}{d\xi} + \frac{(n+1)}{6} \biggl[\frac{\sigma_c^2}{\gamma } - \frac{6\alpha}{\xi} \biggl(- \frac{d\theta}{d\xi}\biggr)\biggr] x \, ,</math> </td> </tr> </table> </div> where, <math>~\theta(\xi)</math> is the polytropic Lane-Emden function describing the configuration's unperturbed radial density distribution, and <math>~\gamma</math>, <math>~\sigma_c^2</math>, and <math>~\alpha \equiv (3-4/\gamma)</math> are constants. Here we seek a power-series expression for the displacement function, <math>~x(r)</math>, expanded about the center of the configuration, that approximately satisfies this LAWE. First we note that, near the center, an accurate [[#PolytropicLaneEmden|power-series expression for the polytropic Lane-Emden function]] is, <table border="0" cellpadding="5" align="center"> <tr> <td align="right"> <math>~\theta</math> </td> <td align="center"> <math>~=</math> </td> <td align="left"> <math>~ 1 - \frac{\xi^2}{6} + \frac{n}{120} \xi^4 - \frac{n}{378} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^6 + \cdots </math> </td> </tr> </table> Hence, <div align="center"> <table border="0" cellpadding="5" align="center"> <tr> <td align="right"> <math>~-\frac{d\theta}{d\xi}</math> </td> <td align="center"> <math>~\approx</math> </td> <td align="left"> <math>~ \frac{1}{3} \biggl[ \xi - \frac{n}{10} \xi^3 + \frac{n}{21} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^5 \biggr] \, .</math> </td> </tr> </table> </div> Therefore, near the center of the configuration, the LAWE may be written as, <div align="center"> <table border="0" cellpadding="5" align="center"> <tr> <td align="right"> <math>~6~\theta \frac{d^2x}{d\xi^2} + \biggl\{ 12~\theta - (n+1)\xi \biggl[ \xi - \frac{n}{10} \xi^3 + \frac{n}{21} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^5 \biggr] \biggr\} \frac{2}{\xi}\frac{dx}{d\xi}</math> </td> <td align="center"> <math>~\approx</math> </td> <td align="left"> <math>~ - (n+1) \biggl\{ \frac{\sigma_c^2}{\gamma } - \frac{2\alpha}{\xi} \biggl[ \xi - \frac{n}{10} \xi^3 + \frac{n}{21} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^5 \biggr] \biggr\} x </math> </td> </tr> <tr> <td align="right"> <math>\Rightarrow~~~ ~6\biggl[ 1 - \frac{\xi^2}{6} + \frac{n}{120} \xi^4 \biggr] \frac{d^2x}{d\xi^2} + \biggl\{ 12 \biggl[ 1 - \frac{\xi^2}{6} + \frac{n}{120} \xi^4 \biggr] - (n+1)\biggl[ \xi^2 - \frac{n}{10} \xi^4 \biggr] \biggr\} \frac{2}{\xi}\frac{dx}{d\xi}</math> </td> <td align="center"> <math>~\approx</math> </td> <td align="left"> <math>~ - (n+1) \biggl\{ \mathfrak{F} + 2\alpha \biggl[ \frac{n}{10} \xi^2 - \frac{n}{21} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^4 \biggr] \biggr\} x </math> </td> </tr> <tr> <td align="right"> <math>\Rightarrow~~~ ~\biggl( 6 - \xi^2 + \frac{n}{20} \xi^4 \biggr) \frac{d^2x}{d\xi^2} + \biggl[ 12 - (n+3)\xi^2 + \frac{n(n+2)}{10} \xi^4 \biggr] \frac{2}{\xi}\frac{dx}{d\xi}</math> </td> <td align="center"> <math>~\approx</math> </td> <td align="left"> <math>~ - (n+1) \biggl[ \mathfrak{F} + \frac{n\alpha}{5} \xi^2 - \frac{2n\alpha}{21} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^4 \biggr] x \, ,</math> </td> </tr> </table> </div> where, <math>\mathfrak{F} \equiv (\sigma_c^2/\gamma - 2\alpha)</math> and, for present purposes, we have kept terms in the series no higher than <math>~\xi^4</math>. ====Displacement Finite at Center==== Let's adopt a power-series expression for the displacement function of a form that is finite at the center of the configuration, namely, <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>~ 1 + a\xi + b\xi^2 + c\xi^3 + d\xi^4 + e\xi^5 + f\xi^6\cdots </math> </td> </tr> <tr> <td align="right"> <math>~\Rightarrow ~~~ \frac{1}{\xi}\frac{dx}{d\xi}</math> </td> <td align="center"> <math>~=</math> </td> <td align="left"> <math>~ \frac{a}{\xi} + 2b + 3 c\xi + 4d\xi^2 + 5e\xi^3 + 6f\xi^4 +\cdots </math> </td> </tr> </table> </div> and, <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>~ 2b + 6c\xi + 12d\xi^2 + 20e\xi^3 + 30f\xi^4 + \cdots </math> </td> </tr> </table> </div> Substituting these expressions into the LAWE gives, <div align="center"> <table border="0" cellpadding="5" align="center"> <tr> <td align="right"> <math>~\biggl( 6 - \xi^2 + \frac{n}{20} \xi^4 \biggr) \biggl( 2b + 6c\xi + 12d\xi^2 + 20e\xi^3 + 30f\xi^4 \biggr) + \biggl[ 12 - (n+3)\xi^2 + \frac{n(n+2)}{10} \xi^4 \biggr] \biggl( \frac{2a}{\xi} + 4b + 6 c\xi + 8d\xi^2 + 10e\xi^3 + 12f\xi^4 \biggr)</math> </td> <td align="center"> <math>~\approx</math> </td> <td align="left"> <math>~ - (n+1) \biggl[ \mathfrak{F} + \frac{n\alpha}{5} \xi^2 - \frac{2n\alpha}{21} \biggl( \frac{n}{5} - \frac{1}{8} \biggr) \xi^4 \biggr] \biggl( 1 + a\xi + b\xi^2 + c\xi^3 + d\xi^4 \biggr)</math> </td> </tr> </table> </div> Expressions for the various coefficients can now be determined by equating terms on the LHS and RHS that have like powers of <math>~\xi</math>. <div align="center"> <table border="1" cellpadding="5" align="center"> <tr> <td align="center">Term</td> <td align="center">LHS</td> <td align="center">RHS</td> <td align="center">Implication</td> </tr> <tr> <td align="right"> <math>~\xi^{-1}:</math> </td> <td align="center"> <math>~24a</math> </td> <td align="center"> <math>~0</math> </td> <td align="left"> <math>~\Rightarrow ~~~a=0</math> </td> </tr> <tr> <td align="right"> <math>~\xi^{0}:</math> </td> <td align="center"> <math>~(12b + 48b)</math> </td> <td align="center"> <math>~-(n+1)\mathfrak{F}</math> </td> <td align="left"> <math>~\Rightarrow ~~~b = - \frac{(n+1)\mathfrak{F}}{60}</math> </td> </tr> <tr> <td align="right"> <math>~\xi^{1}:</math> </td> <td align="center"> <math>~[36c+72c-2a(n+3)]</math> </td> <td align="center"> <math>~-a(n+1)\mathfrak{F}</math> </td> <td align="left"> <math>~\Rightarrow ~~~108c = 2a(n+3)-a(n+1)\mathfrak{F} \Rightarrow~~c=0</math> </td> </tr> <tr> <td align="right"> <math>~\xi^{2}:</math> </td> <td align="center"> <math>~[72d-2b+96d-4b(n+3)]</math> </td> <td align="center"> <math>~\biggl[-b(n+1)\mathfrak{F}-\frac{n(n+1)\alpha}{5}\biggr]</math> </td> <td align="left"> <math>~\Rightarrow ~~~d = - (n+1)\biggl\{ \frac{n\alpha +\mathfrak{F}[(4n+14)-(n+1)\mathfrak{F} ]}{10080} \biggr\}</math> </td> </tr> </table> </div> In summary, the desired, approximate power-series expression for the polytropic displacement function is: <div align="center" id="PolytropicDisplacement"> <table border="1" width="80%" cellpadding="8" align="center"><tr><td align="center"> <table border="0" cellpadding="5" align="center"> <tr> <td align="right"> <math>~x(\xi)</math> </td> <td align="center"> <math>~=</math> </td> <td align="left"> <math>~ 1 - \frac{(n+1)\mathfrak{F}}{60} \xi^2- (n+1)\biggl\{ \frac{n\alpha +\mathfrak{F}[(4n+14)-(n+1)\mathfrak{F} ]}{10080} \biggr\} \xi^4 + \cdots </math> </td> </tr> </table> </td></tr></table> </div>
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)
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