SSC/Stability/n1PolytropeLAWE/Pt2

From JETohlineWiki
Jump to navigation Jump to search

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⁡ξ)ξ2x⋅d2xdξ2+2[sin⁡ξ+ξcos⁡ξ]ξx⋅dxdξ+[σ2ξ3−2α(sin⁡ξ−ξcos⁡ξ)]=0,


where,

σ2

≡

ω22πGρcγg,

α

≡

3−4γg.

Now, in addition to recognizing that,

ξx⋅dxdξ

=

dln⁡xdln⁡ξ,

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

ξ2x⋅d2xdξ2

=

ddln⁡ξ[dln⁡xdln⁡ξ]−[1−dln⁡xdln⁡ξ]⋅dln⁡xdln⁡ξ.

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,

dln⁡xdln⁡ξ=c0,

and the two key derivative terms will be,

ξx⋅dxdξ=c0,

      and      

ξ2x⋅d2xdξ2=c0(c0−1).

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

−σ2ξ3

=

c0(c0−1)sin⁡ξ+2c0[sin⁡ξ+ξcos⁡ξ]−2α(sin⁡ξ−ξcos⁡ξ)

 

=

sin⁡ξ[c0(c0−1)+2c0−2α]+ξcos⁡ξ[2(c0+α)]

 

=

sin⁡ξ[c02+c0−2α]+ξ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ξc0−1[a0+b0sin⁡ξ+d0ξcos⁡ξ]

 

=

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

⇒dln⁡xdln⁡ξ

=

ξ[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ξ+[σ2−2α(sin⁡ξ−ξcos⁡ξξ3)]x

Notice that,

ddξ[sin⁡ξξ]

=

−sin⁡ξξ2+cos⁡ξξ

 

=

[ξcos⁡ξ−sin⁡ξξ2].

And, hence,

d2dξ2[sin⁡ξξ]

=

ddξ[cos⁡ξξ−sin⁡ξξ2]

 

=

−cos⁡ξξ2−sin⁡ξξ+2sin⁡ξξ3−cos⁡ξξ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⁡ξξ)]+x⋅d2dξ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ξ+[σ2−2α(sin⁡ξ−ξcos⁡ξξ3)]x

 

=

d2dξ2{(sin⁡ξξ)x}+4[sin⁡ξξ2]dxdξ+{sin⁡ξξ+σ2−2(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⁡ξξ+σ2−2(1+α)(sin⁡ξ−ξcos⁡ξξ3)}⋅x

 

=

d2dξ2{(sin⁡ξξ)x}+4ξ⋅ddξ[(sin⁡ξξ)x]+{sin⁡ξξ+σ2+[4−2(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ξ+x⋅ddξ[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⁡ξξ2−4sin⁡ξξ3+sin⁡ξξ+2(1−α)(sin⁡ξ−ξcos⁡ξξ3)]⋅x

 

=

(sin⁡ξξ)d2xdξ2+[2cos⁡ξξ+2sin⁡ξξ2]⋅dxdξ+[−2(sin⁡ξ−ξcos⁡ξξ3)+(2−2α)(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ξ+[σ2−2α(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 3≤n≤∞. 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      

xP≡3(n−1)2n[1+(n−3n−1)(1ξθn)dθdξ]

if the adiabatic exponent is assigned the value, γg=(n+1)/n, in which case the parameter, α=(3−n)/(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

=

A−B[(1ξθ)dθdξ]

 

=

A−B[(1sin⁡ξ)ddξ(sin⁡ξξ)]

 

=

A−B(1sin⁡ξ)[cos⁡ξξ−sin⁡ξξ2]

 

=

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

in which case,

dxdξ

=

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

 

=

−Bξ3{(cos⁡ξsin2ξ)[−ξ2cos⁡ξ+ξsin⁡ξ]+[2−ξ2−2ξ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λB⋅dxdξ

=

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

Probably this also means,

1λ2B⋅dxdξ

=

1(λξ)4{6−2(λξ)2−2λξ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↔(3−n)/(n−1),

dxdξ

=

(3−nn−1){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⁡ξ)2−2(ξ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−ξ2−2ξ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⁡ξ]+[6−3ξ2−6ξ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⁡ξ]+[6−3ξ2−6ξ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{6−2ξ2−2ξ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↔(3−n)/(n−1),

d2xdξ2

=

(n−3n−1){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+ξ2−2ξ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{6−2ξ2−2ξ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⇒α=(3−4/γg)=+1, and,

Q≡−dln⁡θ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+[4−2Q]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⁡ξ]{ξ2−2+ξ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+ξ2−2+ξcos⁡ξsin⁡ξ+(ξcos⁡ξsin⁡ξ)2}

 

 

−2Bξ4[1−2ξ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 |