SSC/Stability/n1PolytropeLAWE/Pt2

From JETohlineWiki
Revision as of 22:08, 20 January 2024 by Joel2 (talk | contribs)
Jump to navigation Jump to search

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 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

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

Development

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

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

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 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

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

Tiled Menu

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