SSC/Stability/BiPolytropes/51Models: Difference between revisions

From jetwiki
Jump to navigation Jump to search
Line 287: Line 287:
==Stability==
==Stability==


Here we solve the LAWE numerically using a 2<sup>nd</sup>-order accurate,  [[Appendix/Ramblings/51BiPolytropeStability/BetterInterfacePt2#Convert_to_Implicit_Approach|implicit integration scheme]].
Here we solve the LAWE numerically (on a uniformly zoned mesh &#8212; different <math>\Delta\tilde{r}</math> for the separate core/envelope regions) using a 2<sup>nd</sup>-order accurate,  [[Appendix/Ramblings/51BiPolytropeStability/BetterInterfacePt2#Convert_to_Implicit_Approach|implicit integration scheme]] in which the LAWE is broken into a pair of 1<sup>st</sup>-order ODEs.  These results should be compared against a separate [[SSC/Stability/BiPolytropes/SuccinctDiscussion#Stability|succinct discussion]] of our analysis obtained from integrating the LAWE in its standard 2<sup>nd</sup>-order ODE form.


=See Also=
=See Also=


{{ SGFfooter }}
{{ SGFfooter }}

Revision as of 15:47, 27 September 2023

BiPolytrope with nc = 5 and ne = 1

Here we construct and analyze the relative stability of a bipolytrope in which the core has an nc=5 polytropic index and the envelope has an ne=1 polytropic index.

Structure

  1. Individual model profiles, taken from:
  2. (q,ν) sequences of fixed μe/μc, taken from:
  3. νmax model, taken from:
    • SSC/Structure/BiPolytropes/Analytic51#Limiting_Mass
       

      Maximum Fractional Core Mass, ν=Mcore/Mtot (solid green circular markers)
      for Equilibrium Sequences having Various Values of μe/μc

      μeμc

      ξi

      θi

      ηi

      Λi

      A

      ηs

      LHS

      RHS

      qrcoreR

      νMcoreMtot

      Extrema along Various Equilibrium Sequences

      13

      --- --- --- --- --- --- --- 0.0 2π

      0.33

      24.00496 0.0719668 0.0710624 0.2128753 0.0726547 1.8516032 -223.8157 -223.8159 0.038378833 0.52024552

      0.316943

      10.744571 0.1591479 0.1493938 0.4903393 0.1663869 2.1760793 -31.55254 -31.55254 0.068652714 0.382383875

      0.31

      9.014959766 --- --- 0.59835053 --- --- --- --- 0.0755022550 0.3372170064

      0.3090

      8.8301772 0.1924833 0.1750954 0.6130669 0.2053811 2.2958639 -18.47809 -18.47808 0.076265588 0.331475715

      14

      4.9379256 0.3309933 0.2342522 1.4179907 0.4064595 2.761622 -2.601255 -2.601257 0.084824137 0.139370157

      Recall that,

      iξi3;       and       m33(μeμc).


    •  
    • SSC/Structure/BiPolytropes/Analytic51Renormalize#Model_Pairings
       
      file = Dropbox/WorkFolder/Wiki edits/Bipolytrope/Stability/qAndNuMax.xlsx --- worksheet = B-KB74 thru MinuPreparation
      file = Dropbox/WorkFolder/Wiki edits/Bipolytrope/Stability/qAndNuMax.xlsx --- worksheet = B-KB74 thru MinuPreparation
      Bipolytrope with (nc,ne)=(5,1)
      Selected Pairings along the μe/μc=0.31 Sequence
      Pairing ξi Λi ν q
      A 9.014959766 0.59835053 0.3372170064 0.0755022550
      B1 9.12744 0.60069262 0.3372001445 0.0746451491
      B2 8.90394 0.59610192 0.33720014467 0.0763642133


      Bipolytropic (5, 1) Equilibrium Sequences
      Bipolytropic (5, 1) Equilibrium Sequences
      Bipolytropic (5, 1) Equilibrium Sequences
      Bipolytropic (5, 1) Equilibrium Sequences

Stability

Here we solve the LAWE numerically (on a uniformly zoned mesh — different Δr~ for the separate core/envelope regions) using a 2nd-order accurate, implicit integration scheme in which the LAWE is broken into a pair of 1st-order ODEs. These results should be compared against a separate succinct discussion of our analysis obtained from integrating the LAWE in its standard 2nd-order ODE form.

See Also

Tiled Menu

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