SSC/Stability/NeutralMode: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Line 2,164: Line 2,164:
</tr>
</tr>
</table>
</table>
We conclude, therefore, that the radial displacement function (''i.e.,'' the eigenfunction) for the neutral mode of all polytropic configurations is,
We conclude, therefore, that the radial displacement function (''i.e.,'' the eigenfunction) for the neutral <math>(\sigma_c^2 = 0)</math> mode of all polytropic configurations is,
<table border="0" cellpadding="5" align="center">
<table border="0" cellpadding="5" align="center">



Revision as of 13:36, 8 September 2021

LAWE

Most General Form

In an accompanying discussion, we derived the so-called,

Adiabatic Wave (or Radial Pulsation) Equation

d2xdr02+[4r0(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(43γg)g0r0]x=0

where the gravitational acceleration,

g0

GMrr02=1ρ0dP0dr0g0ρ0r0P0=dlnP0dlnr0.

The solution to this equation gives eigenfunctions that describe various radial modes of oscillation in spherically symmetric, self-gravitating fluid configurations. The boundary condition conventionally used in connection with the adiabatic wave equation is,

r0dlnxdr0

=

1γg(43γg+ω2R3GMtot)        at         r0=R.

Polytropic Configurations

Part 1

If the initial, unperturbed equilibrium configuration is a polytropic sphere whose internal structure is defined by the function, θ(ξ), that provides a solution to the,

Lane-Emden Equation

1ξ2ddξ(ξ2dΘHdξ)=ΘHn

then,

r0

=

anξ,

ρ0

=

ρcθn,

P0

=

Kρ0(n+1)/n=Kρc(n+1)/nθn+1,

g0

=

GM(r0)r02=Gr02[4πan3ρc(ξ2dθdξ)],

where,

an

=

[(n+1)K4πGρc(1n)/n]1/2.

Hence, after multiplying through by an2, the above adiabatic wave equation can be rewritten in the form,

d2xdξ2+[4ξg0an(an2ρ0P0)]dxdξ+(an2ρ0γgP0)[ω2+(43γg)g0anξ]x

=

0.

In addition, given that,

g0an

=

4πGρc(dθdξ),

and,

an2ρ0P0

=

(n+1)(4πGρc)θ=an2ρcPcθcθ,

we can write,

0

=

d2xdξ2+[4(n+1)V(ξ)ξ]dxdξ+[ω2(an2ρcγgPc)θcθ(34γg)(n+1)V(x)ξ2]x

 

=

d2xdξ2+[4(n+1)V(ξ)ξ]dxdξ+(n+1)[ω2(an2ρcγgPc)ξ2θc(n+1)θ(34γg)V(x)]xξ2

where we have adopted the function notation,

V(ξ)

ξθdθdξ.

Part 2

Drawing from an accompanying discussion, we have the following:

Polytropic LAWE (linear adiabatic wave equation)

0=d2xdξ2+[4(n+1)Q]1ξdxdξ+(n+1)[(σc26γg)ξ2θαQ]xξ2

where:    Q(ξ)dlnθdlnξ,    σc23ω22πGρc,     and,     α(34γg)

In order to reconcile with the "Part 1" expression, we note first that V(ξ)Q(ξ). We note as well that since,

(an2ρcPc)θc

=

(n+1)4πGρc,

we have,

ω2(an2ρcγgPc)ξ2θc(n+1)θ

ω2γg[(n+1)4πGρc]ξ2(n+1)θ=16γg[3ω22πGρc]ξ2θ=(σc26γg)ξ2θ.

All physically reasonable solutions are subject to the inner boundary condition,

dxdξ=0         at         ξ=0,

but the relevant outer boundary condition depends on whether the underlying equilibrium configuration is isolated (surface pressure is zero), or whether it is a "pressure-truncated" configuration. As is the case with the pressure-truncated isothermal spheres, discussed above, if the polytropic configuration is truncated by the pressure, Pe, of a hot, tenuous external medium, then the solution to the LAWE is subject to the outer boundary condition,

dlnxdlnξ=3         at         ξ=ξ~.

But, for isolated polytropes, the sought-after solution is subject to the more conventional boundary condition,

dlnxdlnξ=(3nn+1)+nσc26(n+1)[ξθ']         at         ξ=ξsurf.

Radial Pulsation Neutral Mode

Background

The integro-differential version of the statement of hydrostatic balance is

dPdr=GMrρr2

From our separate discussion, we have found that,

Exact Solution to the (3n<) Polytropic LAWE

σc2=0

      and      

xP3(n1)2n[1+(n3n1)(1ξθn)dθdξ].

Let's rewrite the significant functional term in this expressions in terms of basic variables. That is,

(1ξθn)dθdξ

=

(anρcr0ρ0)g04πGρcan

 

=

M(r0)4πr03ρ0.

Trial Eigenfunction & Its Derivatives

Let's adopt the following trial solution:

xt

=

abMr4πr03ρ0=abg04πGr0ρ0.

Then we have,

(1b)dxtdr0

=

ddr0[Mr4πr03ρ0]

 

=

[14πr03ρ0]dMrdr0[Mr4πr03ρ02]dρ0dr0[3Mr4πr04ρ0]

 

=

1r0[3Mr4πr04ρ0][Mr4πr03ρ02]dρ0dr0

(1b)d2xtdr02

=

ddr0{1r0[3Mr4πr04ρ0][Mr4πr03ρ02]dρ0dr0}

 

=

1r02[34πr04ρ0]dMrdr0+[3Mr4πr04ρ02]dρ0dr0+4[3Mr4πr05ρ0]

 

 

[Mr4πr03ρ02]d2ρ0dr02[14πr03ρ02]dρ0dr0dMrdr0+[3Mr4πr04ρ02]dρ0dr0+[2Mr4πr03ρ03](dρ0dr0)2

 

=

4r02+3Mr4πr05ρ0[4+dlnρ0dlnr0]+1r02[3Mr4πr03ρ01]dlnρ0dlnr0[Mr4πr03ρ02]d2ρ0dr02+[2Mr4πr05ρ0](dlnρ0dlnr0)2.

Given that,

ΔMr4πr03ρ0

=

14πG(g0r0ρ0)=[P04πGr02ρ02dlnP0dlnr0],

these expression can be rewritten as,

(r02b)dxtdr0

=

r0{13ΔΔdlnρ0dlnr0},

and,

(r02b)d2xtdr02

=

4+3Δ[4+dlnρ0dlnr0]+[3Δ1]dlnρ0dlnr0Δ(r02ρ0)d2ρ0dr02+2Δ(dlnρ0dlnr0)2.

Plug Trial Eigenfunction Into LAWE

 

LAWE

d2xdr02+[4r0(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(43γg)g0r0]x=0


Plugging our trial radial displacement function, xt, into the LAWE gives,

LAWE

=

(r02b)d2xtdr02(r0b)[4+dlnP0dlnr0]dxtdr0(r02b)(ρ0γgP0)[(43γg)g0r0+σc2(2πGρc3)]xt

 

=

(r02b)d2xtdr02(r0b)[4+dlnP0dlnr0]dxtdr0+(1b)1γg(dlnP0dlnr0)(43γg)xt(r02b)(ρ0γgP0)[σc2(2πGρc3)]xt

 

=

4+3Δ[4+dlnρ0dlnr0]+[3Δ1]dlnρ0dlnr0Δ(r02ρ0)d2ρ0dr02+2Δ(dlnρ0dlnr0)2

 

 

+[4+dlnP0dlnr0]{13ΔΔdlnρ0dlnr0}+(1b)1γg(dlnP0dlnr0)(43γg)(abΔ)(1b)(ρ0r02γgP0)[σc2(2πGρc3)](abΔ).

Now, if we set σc2=0 and dlnP0/dlnr0=γg(dlnρ0/dlnr0), this expression becomes,

LAWE

=

4+3Δ[4+dlnρ0dlnr0]+[3Δ1]dlnρ0dlnr0Δ(r02ρ0)d2ρ0dr02+2Δ(dlnρ0dlnr0)2

 

 

+[4+dlnP0dlnr0]{13ΔΔdlnρ0dlnr0}+(1b)1γg(dlnP0dlnr0)(43γg)(abΔ)(1b)(ρ0r02γgP0)[σc2(2πGρc3)](abΔ).

Notice that the key components of this last term may be rewritten as,

(ρ0r02γgP0)[σc2(2πGρc3)]

=

(4πGρ02r02P0)[σc26γg(ρcρ0)]

 

=

(1Δ)dlnP0dlnr0[σc26γg(ρcρ0)].

So, for our trial eigenfunction, we have …

LAWE

=

[2Δ1]dlnρ0dlnr0Δ(r02ρ0)d2ρ0dr02+2Δ(dlnρ0dlnr0)2+dlnP0dlnr0{[13ΔΔdlnρ0dlnr0]+(43γg)γg[abΔ]+σc26γg(ρcρ0)[abΔ1]}.

Consider Polytropic Structures

Referring back to, for example, a separate review of polytropic structures, we recognize that,

Δ=14πG(g0r0ρ0)

=

1ξ3[(ξ2dθdξ)]θn=1ξ(dθdξ)θn=θ'ξθn,

dlnρ0dlnr0

=

n,

dlnP0dlnr0

=

(n+1).

Also,

(r02ρ0)d2ρ0dr02

=

(ξ2ρcθn)ddξ[nρcθn1θ']

 

=

(nξ2θn)[(n1)θn2(θ')2+θn1θ']

 

=

(nξ2θ2)[(n1)(θ')2+θθ']

 

=

(nξ2θ2)[(n1)(θ')2(θn+1+2θθ'ξ)]

 

=

(nξ2θ2)[(n1)(ξθnΔ)2+θn+1]

 

=

n(n1)(ξn+1θn1Δ)2+nξ2θn1.

Hence,

LAWE

=

n(2Δ1)Δ(r02ρ0)d2ρ0dr02+2n2Δ+(n+1)[13ΔnΔ]+(n+1){(43γg)γg[abΔ]+σc26γg(ρcρ0)[abΔ1]}

 

=

2nΔn+2n2Δ+n3nΔn2Δ+13ΔnΔΔ(r02ρ0)d2ρ0dr02+(n+1){(43γg)γg[abΔ]+σc26γg(ρcρ0)[abΔ1]}

 

=

1+n2Δ(2n+3)ΔΔ(r02ρ0)d2ρ0dr02+(n+1)[4γg3][abΔ]+(n+1){σc26γg(ρcρ0)[abΔ1]}.

If, γg=(n+1)/n, we can further simplify and obtain,

LAWE

=

1+n2Δ(2n+3)ΔΔ(r02ρ0)d2ρ0dr02+[n3][abΔ]+(n+1){σc26γg(ρcρ0)[abΔ1]}

 

=

1+(n3)ab+n2Δ+(3n)Δ(2n+3)ΔΔ(r02ρ0)d2ρ0dr02+(n+1){σc26γg(ρcρ0)[abΔ1]}

Try Again

General Form of Wave Equation

 

LAWE

d2xdr02+[4r0(g0ρ0P0)]dxdr0+(ρ0γgP0)[ω2+(43γg)g0r0]x=0

Employing the substitutions,

σc2

3ω22πGρc,

α

34γg=3nn+1,

g0

GMrr02=1ρ0dP0dr0g0ρ0r0P0=dlnP0dlnr0,

Δ

Mr4πr03ρ0=14πG(g0r0ρ0)=[P04πGr02ρ02dlnP0dlnr0],

we have,

LAWE =

d2xdr02+1r0[4g0ρ0r0P0]dxdr0+[(4γg3)g0ρ0r0P0]xr02+(ρ0P0)[4πGρc(σc26γg)]x

  =

d2xdr02+1r0[4+dlnP0dlnr0]dxdr0+[αdlnP0dlnr0]xr021Δ[dlnP0dlnr0ρcρ0(σc26γg)]xr02.


In the context of polytropic configurations (see more below), we appreciate that,

ρ0ρc =

θn,

dlnP0dlnr0 =

(n+1)dlnθdlnξ=(n+1)Q,      and,

1ΔdlnP0dlnr0 =

(n+1)ξ2θn1.

Inserting these into the LAWE expression and multiplying through by the square of the polytropic length scale, an2, we obtain,

LAWE =

d2xdξ2+1ξ[4(n+1)Q]dxdξ[α(n+1)Q]xξ2+[(n+1)ξ2θ(σc26γg)]xξ2.


This is identical to what has been referred to in a separate discussion, as the

Polytropic LAWE

0=d2xdξ2+[4(n+1)Q]1ξdxdξ+(n+1)[(σc26γg)ξ2θαQ]xξ2

where:    Q(ξ)dlnθdlnξ,    σc23ω22πGρc,     and,     α(34γg)

Derivatives of Δ

Here we evaluate the first derivative of Δ with respect to r0,

dΔdr0 =

ddr0{Mr4πr03ρ0}

  =

14πr03ρ0dMrdr03Mr4πr04ρ0Mr4πr03ρ02dρ0dr0

  =

1r03Mr4πr04ρ0Mr4πr03ρ02dρ0dr0

  =

1r0{1Δ[3+dlnρ0dlnr0]}

r0dΔdr0 =

[13ΔΔdlnρ0dlnr0];

and the second derivative of Δ with respect to r0,

d2Δdr02 =

ddr0{1r0}ddr0{3Mr4πr04ρ0}ddr0{Mr4πr03ρ02dρ0dr0}

  =

1r0234π{1r04ρ0dMrdr04Mrr05ρ0Mrr04ρ02dρ0dr0}

   

14π{1r03ρ02dρ0dr0dMrdr03Mrr04ρ02dρ0dr02Mrr03ρ03[dρ0dr0]2+Mrr03ρ02d2ρ0dr02}

  =

1r023r02{1Δ[4+dlnρ0dlnr0]}1r02{[13Δ]dlnρ0dlnr02Δ[dlnρ0dlnr0]2+Δ(r02ρ0)d2ρ0dr02}

r02d2Δdr02 =

1+33Δ[4+dlnρ0dlnr0]+[13Δ]dlnρ0dlnr02Δ[dlnρ0dlnr0]2+Δ(r02ρ0)d2ρ0dr02

  =

412Δ+[16Δ]dlnρ0dlnr02Δ[dlnρ0dlnr0]2+Δ(r02ρ0)d2ρ0dr02.

Trial Eigenfunction

As above, let's adopt a trial eigenfunction of the form,

xt

=

abMr4πr03ρ0=abΔ.

Then we have,

1b[r02×LAWE]trial =

r02d2Δdr02[4+dlnP0dlnr0]r0dΔdr0+[αdlnP0dlnr0](abΔ)1Δ[dlnP0dlnr0ρcρ0(σc26γg)](abΔ)

  =

{412Δ+[16Δ]dlnρ0dlnr02Δ[dlnρ0dlnr0]2+Δ(r02ρ0)d2ρ0dr02}[4+dlnP0dlnr0][13ΔΔdlnρ0dlnr0]

   

+[αdlnP0dlnr0](abΔ)1Δ[dlnP0dlnr0ρcρ0(σc26γg)](abΔ).

Assume Polytropic Relations

If we assume that the equilibrium models are polytropes, then we know that,

ρ0θn

           

dlnρ0dlnr0=ndlnθdlnξ;

P0θn+1

           

dlnP0dlnr0=(n+1)dlnθdlnξ.

We also deduce that,

(r02ρ0)d2ρ0dr02

=

(ξ2θn)ddξ[dθndξ]

 

=

(ξ2θn)ddξ[nθn1θ']

 

=

(nξ2θn)[(n1)θn2(θ')2+θn1θ']

 

=

[n(n1)ξ2θ2](θ')2(nξ2θ)[θn+2ξθ']

 

=

n(n1)[dlnθdlnξ]2nξ2θn12ndlnθdlnξ

 

=

n(n1)Ξ2nξ2θn12nΞ,

where we have introduced the shorthand notation,

Ξ

dlnθdlnξ.


Drawing from our accompanying discussion, for example, we note as well that,

Δ

=

[P04πGr02ρ02dlnP0dlnr0]

1ΔdlnP0dlnr0

=

4πGr02ρ02P0

 

=

4πG(an2ξ2)(ρcθn)2[K1ρc(n+1)/nθ(n+1)]

 

=

4πG[(n+1)K4πGρc(1n)/n](ρc)2[K1ρc(n+1)/n]ξ2θ(n1)

 

=

(n+1)ξ2θ(n1)

Δξ2θn1

=

1(n+1)dlnP0dlnr0=Ξ.

Hence,

1b[r02×LAWE]trial =

412Δ+[16Δ]nΞ2Δ[nΞ]2+Δ[n(n1)Ξ2nξ2θn12nΞ][4+(n+1)Ξ][13ΔΔnΞ]

   

+[α(n+1)Ξ](abΔ)1Δ[(n+1)Ξρcρ0(σc26γg)](abΔ)

  =

nΞ6ΔnΞ2Δ[nΞ]2+Δ[n(n1)Ξ2nξ2θn1+2nΞ]+2(n+1)Ξ+n(n+1)Ξ2Δ

   

+[α(n+1)Ξ](ab)[α(n+1)Ξ]Δ1Δ[(n+1)Ξρcρ0(σc26γg)](abΔ)

  =

Ξ[1+n+2(n+1)+α(n+1)(ab)]+Δ[[4n+α(n+1)]Ξ+n(n1)Ξ2+n(n+1)Ξ22n2Ξ2]

   

1Δ[(n+1)Ξρcρ0(σc26γg)](abΔ)

  =

Ξ[3(n+1)+α(n+1)(ab)]Δ[4n+α(n+1)]Ξ1Δ[(n+1)Ξρcρ0(σc26γg)](abΔ)

Third Time

General Relations

Various general relations taken from above derivations:

LAWE =

d2xdr02+1r0[4+dlnP0dlnr0]dxdr0+[αdlnP0dlnr0]xr021Δ[dlnP0dlnr0ρcρ0(σc26γg)]xr02,

where,

Δ

Mr4πr03ρ0

1ΔdlnP0dlnr0 =

4πGr02ρ02P0;

r0dΔdr0 =

[13ΔΔdlnρ0dlnr0];

r02d2Δdr02 =

412Δ+[16Δ]dlnρ0dlnr02Δ[dlnρ0dlnr0]2+Δ(r02ρ0)d2ρ0dr02.

Polytropes

If polytropic relations are adopted:

Δ =

Qξ2θn1;

1ΔdlnP0dlnr0 =

(n+1)ξ2θn1;

r02×LAWE =

r02d2xdr02+[4(n+1)Q]r0dxdr0+[(n3)Q]x+1Δ[(n+1)Q1θn(σc26γg)]x;

r0dΔdr0 =

[13Δ+nΔQ]=1+(nQ3)Δ;

r02d2Δdr02 =

412Δn[16Δ]Q2n2ΔQ2+Δ(r02ρ0)d2ρ0dr02

  =

412Δn[16Δ]Q2n2ΔQ2+Δ[n(n1)Q2+2nQnQΔ]

  =

4nQ12Δ+6nQΔ2n2ΔQ2+Δ[n(n1)Q2+2nQ]nQ

  =

42nQ+Δ[6nQ2n2Q2+n(n1)Q2+2nQ12].

  =

42nQ+Δ[8nQn2Q2nQ212].

Eigenfunction Choice

Again, let's try the trial eigenfunction,

xt =

abΔ,

in which case,

1b[r02×LAWE] =

r02d2Δdr02[4(n+1)Q]r0dΔdr0+[(n3)Q](abΔ)+1Δ[(n+1)Q1θn(σc26γg)](abΔ)

  =

42nQ+Δ[8nQn2Q2nQ212][4(n+1)Q][1+(nQ3)Δ]+[(n3)Q](abΔ)

   

+1Δ[(n+1)Q1θn(σc26γg)](abΔ)

  =

42nQ+Δ[8nQn2Q2nQ212]4(4nQ12)Δ+(n+1)Q+(n+1)(nQ23Q)Δ+(n3)Q(ab)(n3)QΔ

   

+1Δ[(n+1)Q1θn(σc26γg)](abΔ)

  =

[(n+1)2n+(n3)(ab)]Q+Δ[8nQn2Q2nQ2124nQ+12+(n+1)(nQ2)3Q(n+1)+(3n)Q]

   

+1Δ[(n+1)Q1θn(σc26γg)](abΔ)

  =

[(1n)+(n3)(ab)]Q+Δ[0]+1Δ[(n+1)Q1θn(σc26γg)](abΔ).

Hence, we are left with only the σc2 term if we set,

0 =

(1n)+(n3)(ab)

(ab) =

n1n3.

We conclude, therefore, that the radial displacement function (i.e., the eigenfunction) for the neutral (σc2=0) mode of all polytropic configurations is,

xneutral =

1b[n1n3Δ]=1b[n1n3Qξ2θn1]

  =

1b(n1n3)[1+(n3n1)1ξ2θn1dlnθdlnξ]

  =

1b(n1n3)[1+(n3n1)1ξθndθdξ].

This last expression exactly matches our earlier result found for polytropic configurations if we choose an overall amplitude coefficient of the form,

b =

2n3(n3).

Hooray!

See Also

Tiled Menu

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