SSC/Stability/n1PolytropeLAWE/Pt2
Radial Oscillations of n = 1 Polytropic Spheres (Pt 2)
Part I: Search for Analytic Solutions |
Part II: New Ideas |
III: What About Bipolytropes? |
New Idea Involving Logarithmic Derivatives
Simplistic Layout
Let's begin, again, with the relevant LAWE, as provided above. After dividing through by , we have,
where,
|
|
|
|
|
|
|
|
Now, in addition to recognizing that,
|
|
|
|
in a separate context, we showed that, quite generally,
|
|
|
|
Hence, if we assume that the eigenfunction is a power-law of , that is, assume that,
then the logarithmic derivative of is a constant, namely,
and the two key derivative terms will be,
|
|
and |
|
In this case, the LAWE is no longer a differential equation but, instead, takes the form,
|
|
|
|
|
|
|
|
|
|
|
|
Now, the cosine term will go to zero if ; and the sine term will go to zero if,
|
|
|
|
|
|
|
|
If these two — rather strange — conditions are met, then we have a marginally unstable configuration because, . This, in and of itself, is not very physically interesting. However, it may give us a clue regarding how to more generally search for a physically reasonable radial eigenfunction.
More general Assumption
Try,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Another Viewpoint
Development
Multiplying through the above LAWE by gives,
|
|
|
|
Notice that,
|
|
|
|
|
|
|
|
And, hence,
|
|
|
|
|
|
|
|
|
|
|
|
So, we can write,
|
|
|
|
|
|
|
|
|
|
|
|
This means that we can rewrite the LAWE as,
|
|
|
|
|
|
|
|
We recognize, also, that,
|
|
|
|
|
|
|
|
So the LAWE becomes,
|
|
|
|
|
|
|
|
|
|
|
|
where we have introduced the new, modified eigenfunction,
Alternatively, the LAWE may be written as,
|
|
|
|
or,
|
|
|
|
|
|
|
|
Now, if we adopt the homentropic convention that arises from setting, , then for our polytropic configuration, we should set, and, hence, . This will mean that the lat term in this LAWE naturally goes to zero. Hence, we have,
|
|
|
|
or,
|
|
|
|
or,
|
|
|
|
Does this help?
Check for Mistakes
Given the definition of , its first derivative is,
|
|
|
|
and its second derivative is,
|
|
|
|
|
|
|
|
|
|
|
|
Hence, the "upsilon" LAWE becomes,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
This should be compared with the first expression, above, namely,
|
|
|
|
and it matches! Q.E.D.
Motivated by Yabushita's Discovery
Initial Exploration
This subsection is being developed following our realization — see the accompanying overview — that the eigenfunction is known analytically for marginally unstable, pressure-truncated configurations having . Specifically, from the work of Yabushita (1975) we have the following,
| Exact Solution to the Isothermal LAWE | ||
|
|
and |
|
And from our own recent work, we have discovered the following,
| Precise Solution to the Polytropic LAWE | ||
|
|
and |
|
if the adiabatic exponent is assigned the value, , in which case the parameter, . Using this polytropic displacement function as a guide, let's try for the case of , an expression of the form,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
in which case,
|
|
|
|
|
|
|
|
|
|
|
|
|
What if, instead, we try the more generalized form,
Then we have,
Probably this also means,
|
Let's check against the more general derivation, which gives after recognizing that, ,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
This matches the preceding, direct derivation.
Also,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Let's also check this against the more general derivation, which gives after again recognizing that, ,
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
A cross-check with the first attempt to derive this second derivative expression initially unveiled a couple of coefficient errors. These have now been corrected and both expressions agree.
Succinct Demonstration
Given that, for , we should set , and,
|
|
|
|
If we then employ the displacement function,
|
|
|
|
the LAWE becomes,
|
LAWE |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Pretty amazing degree of cancelation! So the above-hypothesized displacement function does satisfy the , polytropic LAWE — for any value of the coefficient, — if we set and . If we set , the function will be normalized such that it goes to unity at the center. In summary, then, we have,
|
|
|
|
See Also
|
Appendices: | VisTrailsEquations | VisTrailsVariables | References | Ramblings | VisTrailsImages | myphys.lsu | ADS | |