SSC/Structure/BiiPolytropes/FreeEnergy51

From JETohlineWiki
Revision as of 14:58, 18 January 2024 by Joel2 (talk | contribs) (The Envelope)
Jump to navigation Jump to search

Free Energy of BiPolytrope with (nc, ne) = (5, 1)


Part I:  Mass Profile

 


Part II:  Gravitational Potential Energy

 


Part III:  Thermal Energy Reservoir

 


Free Energy
of
Bipolytropes

(nc, ne) = (5, 1)

Here we present a specific example of the equilibrium structure of a bipolytrope as determined from a free-energy analysis. The example is a bipolytrope whose core has a polytropic index, nc=5, and whose envelope has a polytropic index, ne=1. The details presented here build upon an overview of the free energy of bipolytropes that has been presented elsewhere.  
 
 
 
 

Preliminaries

Mass Profile

The Core

The core has nc=5γc=1+1/nc=6/5. Referring to the general relation as established in our accompanying overview, and using ρ0 to represent the central density, we can write,

(For0xq)       Mr

=

Mtot(νq3)(ρ0ρ¯core)eq0x3[ρ(x)ρ0]corex2dx.

Drawing on the derivation of detailed force-balance models of (nc,ne)=(5,1) bipolytropes, the density profile throughout the core is,

[ρ(ξ)ρ0]core

=

(1+13ξ2)5/2,

where the dimensionless radial coordinate is,

ξ

=

[Gρ04/5Kc]1/2(2π3)1/2r.

Switching to the normalizations that have been adopted in the broad context of our discussion of configurations in virial equilibrium and inserting the adiabatic index of the core (γc=6/5) into all normalization parameters, we have,

Rnorm=[(GKc)Mtot2γ]1/(43γ)

Rnorm=(G5Mtot4Kc5)1/2,

ρnorm=34π[Kc3G3Mtot2]1/(43γ)

ρnorm=34π(Kc3G3Mtot2)5/2.

Hence, we can rewrite,

ξ

=

(rRnorm)(ρ0ρnorm)2/5[GKc]1/2(2π3)1/2Rnormρnorm2/5

 

=

r*(ρ0*)2/5[GKc]1/2(2π3)1/2(G5Mtot4Kc5)1/2(34π)2/5(Kc3G3Mtot2)

 

=

r*(ρ0*)2/5[(2π3)5(34π)4]1/10=r*(ρ0*)2/5[π233]1/10.

Now, following the same approach as was used in our introductory discussion and appreciating that our aim here is to redefine the coordinate, ξ, in terms of normalized parameters evaluated in the equilibrium configuration, we will set,

r*

xχeq;

ρ0*

[ρ0ρ¯]core(ρ¯coreρnorm)=[ρ0ρ¯]coreνMtot/(q3Redge3)eqMtot/Rnorm3=νq3[ρ0ρ¯]coreχeq3.

Then we can set,

ξ

=

(3aξ)1/2x,

in which case,

[ρ(x)ρ0]core

=

(1+aξx2)5/2,

where the coefficient,

(3aξ)1/2

χeq[νq3(ρ0ρ¯)coreχeq3]2/5(π233)1/10=χeq1/5[νq3(ρ0ρ¯)core]eq2/5(π233)1/10

aξ

13{χeq1/5[νq3(ρ0ρ¯)core]eq2/5(π233)1/10}2=χeq2/5[νq3(ρ0ρ¯)core]eq4/5(π2336)1/5.

We therefore have,

Mr|core

=

Mtot[νq3(ρ0ρ¯)core]eq0x3(1+aξx2)5/2x2dx

 

=

Mtot[νq3(ρ0ρ¯)core]eq[x3(1+aξx2)3/2].

Note that, when xq, Mr|coreMcore=νMtot. Hence, this last expression gives,

νMtot

=

Mtot[νq3(ρ0ρ¯)core]eq[q3(1+aξq2)3/2]

[(ρ0ρ¯)core]eq

=

(1+aξq2)3/2.

Hence, finally,

Mr|core

=

νMtot(x3q3)[1+aξx21+aξq2]3/2;

and, after the equilibrium radius, χeq, has been determined from the free-energy analysis, the coefficient, aξ, can be determined via the relation,

χeq2

=

(π2336)(νq3)4(1+aξq2)6aξ5.

MORE USEFUL:

Letting, ξ/3,

aξ

=

(iq)2,

𝔣~Mcore

=

(ρ¯ρc)core=(1+i2)3/2.

The Envelope

The envelope has ne=1γe=1+1/ne=2. Again, referring to the general relation as established in our accompanying overview, and continuing to use ρ0 to represent the central density, we can write,

(Forqx1)       Mr

=

Mtot{ν+(1ν1q3)xix3[ρ(x)ρ¯]envx2dx}.

Drawing on the derivation of detailed force-balance models of (nc,ne)=(5,1) bipolytropes, the density profile throughout the envelope is,

[ρ(η)ρ0]env

=

A(μeμc)θi5[sin(ηB)η],

where definitions of the constants A and B are given in an accompanying table of parameter values, and the dimensionless radial coordinate is,

η

=

[Gρ04/5Kc]1/2(μeμc)θi2(2π)1/2r.

Using the same radial and mass-density normalizations as defined, above, for the core, we can write,

η

=

r*(ρ0*)2/5[GKc]1/2(μeμc)θi2(2π)1/2Rnormρnorm2/5

 

=

r*(ρ0*)2/5(34π)2/5(μeμc)θi2(2π)1/2.

Next, we set,

r*

xχeq;

ρ0*

[ρ0ρ¯]env(ρ¯envρnorm)=[ρ0ρ¯]env(1ν)Mtot/[(1q3)Redge3]eqMtot/Rnorm3=1ν1q3(ρ0ρ¯)envχeq3.

Hence, we can write,

η=bηx,

where,

bη

χeq1/5[1ν1q3(ρ0ρ¯)env]2/5(34π)2/5(μeμc)θi2(2π)1/2.

In which case,

[ρ(x)ρ0]env

=

A(μeμc)θi5[sin(bηxB)bηx],

so,

Mr|env

=

Mtot{ν+[(1ν1q3)(ρ0ρ¯)env]eqqx3A(μeμc)θi5[sin(bηxB)bηx]x2dx}

 

=

νMtot+Mtot[(1ν1q3)(ρ0ρ¯)env]eq(μeμc)3Aθi5bηqxsin(bηxB)xdx

 

=

νMtot+Mtot[(1ν1q3)(ρ0ρ¯)env]eq(μeμc)3Aθi5bη[sin(Bbηx)+bηxcos(Bbηx)bη2]qx

 

=

νMtot+Mtot[(1ν1q3)(ρ0ρ¯)env]eq(μeμc)3Aθi5bη3[C1sin(Bbηx)xbηcos(Bbηx)],

where, C1 is a constant obtained by evaluating the integral at the interface (x=xi=q), specifically,

C1sin(Bbηq)+bηqcos(Bbηq).

Now, this expression can be significantly simplified by drawing on earlier results of this section as well as on attributes of the corresponding detailed force-balanced model. First, independent of the specific density profiles that define the structure of a bipolytrope, the ratio of the mean densities of the two structural regions is,

ρ¯eρ¯c

=

q3(1ν)ν(1q3).

Hence the bracketed pre-factor of the second term of the expression for Mr|env may be rewritten as,

[(1ν1q3)(ρ0ρ¯)env]eq

=

[νq3(ρ0ρ¯)core]eq.

But, from the above derivation of the mass profile in the core, we know that,

[(ρ0ρ¯)core]eq

=

(1+aξq2)3/2=(1+13ξi2)3/2=θi3,

where the final step comes from knowledge of the expression for θi drawn from the detailed force-balanced model (see, for example, the associated Parameter Values table). Hence, we can write,

Mr|env

=

νMtot+Mtot(μeμc)3νAθi2(bηq)3[C1sin(Bbηx)xbηcos(Bbηx)],

and note that the expression for the coefficient, bη, becomes simpler as well, specifically,

bη

=

χeq1/5(34π)2/5(μeμc)θi2(2π)1/2(νq3θi3)2/5

χeq

=

bη5(34π23)1/2(μeμc)5ν2θi4q6.

Next — and, again, drawing from knowledge of the internal structure of the detailed force-balanced model, in particular, realizing that,

bηq=ηi=31/2(μeμc)θi2ξi,

— note that the constant, C1, can be rewritten as,

C1

=

bηqcos(bηqB)sin(bηqB)

 

=

ηicos(ηiB)sin(ηiB)

 

=

ηi2A(dϕdη)i=1A(μeμc)231/2θi4ξi3

 

=

1A(μeμc)231/2θi4[31/2(μeμc)1θi2bηq]3

 

=

13Aθi2(μeμc)1(bηq)3,

which means,

Mr|env

=

νMtotνMtot{11C1[sin(Bbηx)+xbηcos(Bbηx)]}

 

=

νMtotC1[sin(Bbηx)+xbηcos(Bbηx)].

MORE USEFUL:

Letting, ξ/3,

bη=ηs

       and       

bηq=ηi=3(μeμc)i(1+i2)1,

𝔣~Menv(ρ¯ρc)env

=

q3(1ν)ν(1q3)𝔣~Mcore

See Also

Tiled Menu

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