User:Tohline/Appendix/Ramblings/PatrickMotl

Discussing Patrick Motl's 2019 Simulations

May 1 (from Patrick)

Hi Joel, I hope things are well for you and yours. I did finally have a couple research students this semester that were able to slog their way into unix, programming, etc. far enough to do some useful things. They ran two spherical bipolytropes with my old cylindrical code. These were n = 5 cores with n = 1 envelopes. No density discontinuity. One model is a little below the stability boundary in xi, the other is a little above the stability boundary. What I can see from the evolutions, especially now that I made plots of the entropy, is that they are convectively unstable and that is just a mess with the core convecting into the envelope. Did you happen to have any thoughts on what might be a better case to run? When I get done grading finals I am going to work through the equations with kappa set explicitly in the envelope so the entropy profile is flat and just live with whatever density discontinuity that gives me. cheers, Patrick

May 4 (from Joel)

Patrick, My initial response to your 1 May email follows. In an effort to make sure we are on the same page when referencing the structural properties and the stability properties of various configurations, my response includes links to various chapters/subsections of my online wiki-based book. Thanks for pursuing this problem. I am still very interested in Its solution. All the best, Joel

Chosen Initial Models

I presume that when you constructed your pair of initial (spherical) models, you used the analytically prescribed properties found in my chapter titled, "BiPolytrope with <math>n_c = 5</math> and <math>n_e=1</math>" and that when you reference the "stability boundary in <math>~\xi</math>" you are drawing from the summary (Table 3 and Figure 3) found near the end of this chapter in the subsection titled, "Stability Condition." In particular, based on the information contained in this chapter, I presume that you are assuming that the marginally unstable model along the <math>~\mu_e/\mu_c = 1</math> sequence is the configuration having the properties …

<math>~(\xi_i, q, \nu) = (2.416, 0.5952, 0.6830) \, ,</math>

and that your pair of initial models have values of <math>~\xi_i</math> that are somewhat greater than and somewhat smaller than 2.416. If this does not properly describe your choice of initial models, please clarify.

Instability Toward Convection

In my chapter titled, "Axisymmetric Instabilities to Avoid," I discuss — among other things — the Schwarzschild Criterion. Toward the very end of that chapter, in a subsection titled, Modeling Implications and Advice, I specifically describe how convection might arise when modeling spherical polytropic configurations. Would you please read this short subsection and let me know whether you agree with my explanation or not?

In the context of this discussion of the Schwarzschild criterion (i.e., convective instabilities) in polytropes, my question to you is: What value of the adiabatic index, <math>~\gamma</math>, did you use when you evolved the bipolytrope? Did you use the same index for the entire configuration, or did you use a value in the envelope that is different from the value used in the core?

If you used <math>~\gamma = 6/5</math> for the core, then the core should have uniform specific entropy and therefore would be (only) marginally stable against convection. Likewise, if you used <math>~\gamma = 2</math> for the envelope, it should be (only) marginally stable against convection. Incidentally, these are the two separate values of the adiabatic index that I used when I used a free-energy analysis to examine stability.

NOTE: In yet another chapter of my book, I discuss the bipolytropic stability analysis that was performed by Murphy & Fiedler (1985b). They examined the stability of bipolytropes having the core/envelope polytropic indexes swapped, that is, their equilibrium models had, <math>~(n_c, n_e) = (1, 5)</math>. Their stability analysis was performed while assuming that the adiabatic index was <math>~\gamma = 5/3</math> throughout the entire configuration. In a short subsection near the beginning of my chapter titled, "Aside Regarding Convectively Unstable Core," I have pointed out that the cores of these models should have all been convectively unstable, according to the Schwarzschild criterion. Please read this short subsection and let me know if you agree with my stated expectation for those models.

Other Approaches to Stability Analysis

After you send me feedback on the comments & questions I have presented, above, I have a good deal more to tell you about my analysis of the stability of our <math>~(n_c, n_e) = (5, 1)</math> bipolytropic configurations. I am still a bit confused, but I'm pretty sure that the transition from stable to unstable configurations along the <math>~\mu_e/\mu_c = 1</math> occurs at a value of <math>~\xi_i \approx 1.67</math>, rather than the value of <math>~\xi_i = 2.416</math> that I obtained via the free-energy analysis. But, as I've said, there is a lot to tell you about this.

Other Suggestions

You asked what alternative models might be examined. It might be smarter to first evolve some simpler models, such as pressure-truncated polytropes — not bipolytropes. But we should discuss this after I read your feedback on the above.

May 4 (from Patrick)

Hi Joel, Replying to your web page The models have xi_interface values of 2.39184 and 2.44016 For the evolutions, I used passive scalars to identify fluids as core or envelope and use effective polytropic indices and exponents based on the proportion of the fluid in a given cell that is core or envelope. There is slightly more to it than that as the “vacuum” is n = 1.5 and a third fluid type. When I do that, the stars are definitely convectively unstable. s.ps is a quick plot of entropy vs radial cell number. The fluid in the envelope is lower in entropy than fluid in the core. The entropy plotted is c_p ln( tau / rho) = R (gamma/(gamma-1)) ln (tau / rho) where R is the gas constant. The movie shows the meridional cross section of density through the star, the core envelope boundary corrugates pretty quickly in the run and then everything goes off the rails. There may be more than one problem but the convective instability is definitely one problem. cheers, Patrick

May 5 (following a phone conversation with Patrick)

Tying Expressions into H_Book Context

In our wiki-based chapter titled, "First Law of Thermodynamics," we have introduced the concept of an entropy tracer, <math>~\tau</math>. In the subsubsection of this chapter that is titled, "Substantiation," we show that an expression for the specific entropy of a fluid element is,

<math>~s = c_P \ln\biggl( \frac{\tau}{\rho} \biggr) + \mathrm{constant} \, .</math>

In addition, from our wiki-based chapter titled, "Ideal Gas Equation of State," we find the relations,

<math>~c_P - c_V </math>

<math>~=</math>

<math>~\frac{\Re}{\bar\mu} </math>

and,

<math>~\gamma_g </math>

<math>~=</math>

<math>~\frac{c_P}{c_V} </math>

<math>~\Rightarrow</math>

<math>~c_P </math>

<math>~=</math>

<math>~\frac{\gamma_g}{(\gamma_g-1)} \biggl( \frac{\Re}{\bar\mu} \biggr) \, .</math>

Hence this expression for the entropy may be rewritten as,

<math>~s = \frac{\gamma_g}{\gamma_g-1} \biggl( \frac{\Re}{\bar\mu} \biggr) \ln\biggl( \frac{\tau}{\rho} \biggr) + \mathrm{constant} \, .</math>

Aside from the factor of <math>~({\bar\mu})^{-1}</math> that appears on the RHS — more on this below — these are the expressions that Patrick has used to generate the s.ps plot, where the (unlabeled) ordinate is the normalized specific entropy, <math>~s/\Re</math>.

At the end of another subsubsection titled, "Initial Recognition," we also find a relevant expression, namely,

<math>~\tau \equiv (\rho\epsilon)^{1/\gamma_g} = \biggl[ \frac{P}{(\gamma_g - 1)} \biggr]^{1/\gamma_g} \, .</math>

Hence, ignoring the additive constant, in general we may write,

<math>~\frac{s}{\Re/\bar{\mu}}</math>	<math>~=</math>	<math>~ \frac{1}{(\gamma_g-1)}\ln \biggl(\frac{\tau}{\rho}\biggr)^{\gamma_g} </math>
	<math>~=</math>	<math>~ \frac{1}{(\gamma_g-1)}\ln \biggl[ \frac{P}{(\gamma_g-1)\rho^{\gamma_g}} \biggr] \, . </math>

Understanding the Step Function at the Core-Envelope Interface

Now, turning to the accompanying tabular summary of the Radial Profile of Various Physical Variables, we are able to determine how the specific entropy behaves throughout the core and, separately, throughout the envelope in <math>~(n_c, n_e) = (5, 1)</math> bipolytropes.

CORE:   Throughout the core, we see that, <math>~P^*</math> <math>~=</math> <math>~\biggl(1+ \frac{\xi^2}{3}\biggr)^{-3}</math>     and     <math>~\rho^*</math> <math>~=</math> <math>~\biggl(1+ \frac{\xi^2}{3}\biggr)^{-5/2}</math>     and     <math>~\gamma_g</math> <math>~=</math> <math>~\frac{6}{5} \, .</math> Hence, independent of the radial location, <math>~\xi</math>, throughout the core, <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{core}</math> <math>~=</math> <math>~ 5\ln (5) \, . </math> ENVELOPE:   Throughout the envelope, we see that, <math>~P^*</math> <math>~=</math> <math>~\theta_i^6 [\phi(\eta)]^2</math>     and     <math>~\rho^*</math> <math>~=</math> <math>~\biggl(\frac{\mu_e}{\mu_c}\biggr)\theta_i^5 [\phi(\eta)]</math>     and     <math>~\gamma_g</math> <math>~=</math> <math>~2 \, .</math> Hence, independent of the radial location, <math>~\eta</math>, throughout the envelope, <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{env}</math> <math>~=</math> <math>~ \ln \biggl[ \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-2} \theta_i^{-4} \biggr] </math> <math>~=</math> <math>~ \ln \biggl[ \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-2} \biggl( 1 + \frac{\xi_i^2}{3}\biggr)^2 \biggr] \, . </math>

It is therefore clear that the core is a uniform specific-entropy sphere and the envelope is a uniform specific-entropy spherical shell, but in general the specific entropy of material in the envelope is different from the specific entropy of material in the core.

According to our free-energy based evaluation of the stability of bipolytropes having <math>~(n_c, n_e) = (5, 1)</math>, the marginally unstable model along the <math>~\mu_e/\mu_c = 1</math> sequence has the following properties: <math>~(\xi_i, q, \nu, \rho_c/\bar\rho) = (2.41610822, 0.59520261, 0.68306067, 16.3788)</math>; the red arrow in the following diagram points to this model's position in the <math>~q-\nu</math> parameter plane. In an effort to test whether or not this model does identify the transition from stable to unstable configurations along the <math>~\mu_e/\mu_c = 1</math> sequence, Patrick picked a pair of models (highlighted in green in the following table) that straddle the location of the marginally unstable model, and followed the dynamical evolution of both models using a fully 3-D hydrodynamics code. In particular, for the pair of models that Patrick evolved, we find:

<math>~\gamma_c = 6/5</math>     and     <math>~\gamma_e = 2</math>   Initial Model Parameters for Patrick's Pair of Simulations (green background) <math>~\frac{\mu_e}{\mu_c}</math> <math>~\xi_i</math> <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{core}</math> <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{env}</math> 1 2.39184 8.04719 2.13422 1 2.41610822 8.04719 2.16080 1 2.44016 8.04719 2.18706 For more details, see the accompanying discussion titled, Free-Energy Stability Evaluation

These tabulated values of the normalized specific entropy in the core and, separately, in the envelope — also see the plot shown here on the right — appear to be consistent with Patrick's s.ps plot of specific entropy. In particular, this confirms that a step function should appear at the core-envelope interface and that the specific entropy of the envelope material should be lower than the specific entropy of the core material. Therefore, the Schwarzschild criterion is violated at the interface and we should not have been surprised to see convective motions develop — initially, only at the interface.

Can the Step Function be Flipped or Erased

Assume a Mean-Molecular-Weight Ratio of Unity

If we continue to examine equilibrium models that have <math>~\mu_e/\mu_c = 1</math>, is there a value of the interface radius, <math>~\xi_i</math>, for which the entropy step-function disappears and above which the step function flips? The answer appears to be, "Yes." It occurs along the equilibrium sequence where the normalized specific entropy has the same value in the core and as in the envelope. That is, it occurs when,

<math>~\gamma_c = 6/5</math>     and     <math>~\gamma_e = 2</math>   <math>~ \ln \biggl[ \biggl( 1 + \frac{\xi_i^2}{3}\biggr)^2 \biggr]_\mathrm{smooth} </math> <math>~=</math> <math>~5\ln(5)</math> <math>~ \Rightarrow ~~~ \biggl( 1 + \frac{\xi_i^2}{3}\biggr)_\mathrm{smooth} </math> <math>~=</math> <math>~5^{5 / 2}</math> <math>~ \Rightarrow ~~~ [\xi_i ]_\mathrm{smooth} </math> <math>~=</math> <math>~\biggl[3 \biggl(5^{5 / 2} - 1 \biggr)\biggr]^{1 / 2} \approx 12.83375</math> …

Prediction:   Any initial model with <math>~\mu_e/\mu_c = 1</math> and <math>~\xi_i > [\xi_i]_\mathrm{smooth}</math> will be stable against convection, but will be globally dynamically unstable.

Pick a Different Molecular-Weight Ratio

Caution:   This subsection contains a derivation that is based on an interpretation of the specific entropy normalization at the interface that may not be physically correct. It needs more thought.

In principle, we can determine in a similar fashion the values of <math>~[\xi_i]_\mathrm{smooth}</math> that are relevant to equilibrium model sequences having <math>~\mu_e/\mu_c < 1</math>. But in doing this, we must take into account that in most of our above derivations the mean-molecular-weight appears in the denominator of the (LHS) expression for the normalized specific entropy. More generally, the prescription for <math>~[\xi_i]_\mathrm{smooth}</math> should come from the demand that,

<math>~0</math>	<math>~=</math>	<math>~\frac{s_\mathrm{env}}{\Re} - \frac{s_\mathrm{core}}{\Re}</math>
<math>~\Rightarrow ~~~ 0</math>	<math>~=</math>	<math>~\frac{1}{\mu_c} \biggl\{\biggl( \frac{\mu_c}{\mu_e} \biggr) \ln \biggl[ \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-2} \biggl( 1 + \frac{\xi_i^2}{3}\biggr)^2 \biggr]_\mathrm{smooth} - 5 \ln(5) \biggr\}</math>
<math>~\Rightarrow ~~~ \ln \biggl[ \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-2} \biggl( 1 + \frac{\xi_i^2}{3}\biggr)^2 \biggr]_\mathrm{smooth} </math>	<math>~=</math>	<math>~ \ln\biggl[ 5 \biggr]^{5(\mu_e/\mu_c)} </math>
<math>~\Rightarrow ~~~ \biggl( 1 + \frac{\xi_i^2}{3}\biggr)_\mathrm{smooth} </math>	<math>~=</math>	<math>~ 5^{2.5(\mu_e/\mu_c)} \biggl(\frac{\mu_e}{\mu_c}\biggr) </math>
<math>~\Rightarrow ~~~ [\xi_i]_\mathrm{smooth} </math>	<math>~=</math>	<math>~ 3^{1 / 2}\biggl[ 5^{2.5(\mu_e/\mu_c)} \biggl(\frac{\mu_e}{\mu_c}\biggr) -1\biggr]^{1 / 2} \, .</math>

Here are a few examples:

<math>~\gamma_c = 6/5</math>     and     <math>~\gamma_e = 2</math>   <math>~\frac{\mu_e}{\mu_c}</math> <math>~[\xi_i]_\mathrm{smooth}</math> 1 12.83375 1/2 2.86620 1/3 0.90754

Use the Same Ratio of Specific Heats Throughout

Let's examine the initial model's entropy profile under the assumption that the system is evolved with <math>~\gamma_g = 5/3</math> throughout the bipolytrope. From the above analysis, in this case the relevant general expression for the specific entropy profile should be,

<math>~\frac{s}{\Re/\bar{\mu}}</math>

<math>~=</math>

<math>~ \frac{3}{2} \ln \biggl[ \frac{3P}{2\rho^{5/3}} \biggr] </math>

<math>~=</math>

<math>~ \ln \biggl[ \biggl(\frac{3P}{2}\biggr)^{3 / 2} \rho^{-5/2} \biggr] \, . </math>

CORE:   Given that, throughout the core, <math>~P^*</math> <math>~=</math> <math>~\biggl(1+ \frac{\xi^2}{3}\biggr)^{-3}</math>     and     <math>~\rho^*</math> <math>~=</math> <math>~\biggl(1+ \frac{\xi^2}{3}\biggr)^{-5/2} \, ,</math> we have, <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{core}</math> <math>~=</math> <math>~ \ln \biggl\{ \biggl(\frac{3}{2}\biggr)^{3 / 2} \biggl[ \biggl(1+ \frac{\xi^2}{3}\biggr)^{-9/2} \biggr] \biggl[ \biggl(1+ \frac{\xi^2}{3}\biggr)^{25/4} \biggr] \biggr\} = \frac{1}{4}\cdot\ln \biggl[ \biggl(\frac{3}{2}\biggr)^{6} \biggl(1+ \frac{\xi^2}{3}\biggr)^{7} \biggr] \, . </math> ENVELOPE:   Given that, throughout the envelope, <math>~P^*</math> <math>~=</math> <math>~\theta_i^6 [\phi(\eta)]^2</math>     and     <math>~\rho^*</math> <math>~=</math> <math>~\biggl(\frac{\mu_e}{\mu_c}\biggr)\theta_i^5 [\phi(\eta)] \, ,</math> we have, <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{env}</math> <math>~=</math> <math>~ \ln \biggl\{ \biggl(\frac{3}{2}\biggr)^{3 / 2} \biggl[ \theta_i^6 [\phi(\eta)]^2 \biggr]^{3/2} \biggl[ \biggl(\frac{\mu_e}{\mu_c}\biggr)\theta_i^5 [\phi(\eta)] \biggr]^{-5/2} \biggr\} </math>   <math>~=</math> <math>~ \frac{1}{2} \cdot \ln \biggl\{ \biggl(\frac{3}{2}\biggr)^{3 } \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-5} \theta_i^{-7} [\phi(\eta)] \biggr\} </math>   <math>~=</math> <math>~ \frac{1}{4} \cdot \ln \biggl\{ \biggl(\frac{3}{2}\biggr)^{6 } \biggl(1+ \frac{\xi_i^2}{3}\biggr)^{7} \biggl(\frac{\mu_e}{\mu_c}\biggr)^{-10} [\phi(\eta)]^2 \biggr\} \, . </math>

Eigenvectors of Marginally Unstable Models

In preparation for our examination of the relative stability of bipolytropic structures having <math>~(n_c, n_e) = (5, 1)</math> — via numerical integration of the Linear-Adiabatic Wave Equation (LAWE) — we have demonstrated that we understand technically how to solve this type of eigenvalue problem by quantitatively reproducing related analyses that have been previously published by other groups.

Good Comparisons With Previously Published Studies

• Schwarzschild (1941)

  Taff & Van Horn (1974)

  Murphy & Fiedler (1985b)

  Schwarzschild (1941) numerically integrated the LAWE for isolated <math>~n=3</math> polytropic spheres to find eigenvectors (i.e., the spatially discrete eigenfunction and corresponding eigenfrequency) for five separate oscillation modes (the fundamental mode, plus the 1^st, 2^nd, 3^rd, and 4^th overtones) for models having four different adopted adiabatic indexes <math>~\gamma_g = \tfrac{4}{3}, \tfrac{10}{7}, \tfrac{20}{13}, \tfrac{5}{3})</math>. In an accompanying chapter of this H_Book, we demonstrate that we have been able to reproduce in detail Schwarzschild's results for the specific case of <math>~\gamma_g = \tfrac{20}{13}</math>.

• Taff & Van Horn (1974) examined radial oscillations in pressure-truncated isothermal spheres, assuming that the configurations remain isothermal during the oscillations. For models having nine different truncation radii — chosen to straddle the position along the equilibrium sequence where the marginally unstable model was expected to arise — they determined and published the fundamental-mode eigenvalues. For three of these models they also determined and published eigenvalues for the first harmonic mode of oscillation; the radial eigenfunctions associated with both the fundamental mode and the first harmonic mode of these three models also has been displayed in their Figure 1. In a separate accompanying discussion, we demonstrate that we have been able to reproduce in detail the subset of eigenfunctions and associated eigenvalues that have been previously published by Taff & Van Horn (1974).
• In their published study of bipolytropes having <math>~(n_c, n_e) = (1,5)</math> with <math>~\mu_e/\mu_c = 1</math>, Murphy & Fiedler (1985b) integrated a coupled pair of LAWEs — one for the core and another for the envelope — to determine the eigenfunctions and corresponding eigenvalues of various radial modes of oscillation in more than a dozen different equilibrium models, assuming that during the oscillations, <math>~\gamma_g = 5/3</math> throughout both the core and the envelope. In an accompanying chapter of this H_Book titled, Review of the BiPolytrope Stability Analysis by Murphy & Fiedler (1985b), we show that we have been able to duplicate in quantitative detail the eigenvectors associated with their equilibrium Models 10 and 17.

Building upon this set of successful comparisons with stability analyses published by other groups, we have carried out numerical integrations of the relevant LAWE to identify the eigenvectors associated with the fundamental-mode of radial oscillation in pressure-truncated, <math>~n = 5</math> polytropic configurations. Details of this analysis are provided in yet another chapter of this H_Book. The following animation sequence illustrates the results of this analysis. As far as we have been able to determine, an analysis of this type has not previously been conducted for pressure-truncated, <math>~n = 5</math> polytropes.

In a subsection of this separate chapter, we have also shown that, at the maximum-mass turning point along the pressure-truncated <math>~n=5</math> equilibrium sequence — identified by the green circular marker in the left-hand panel of this animation — the fundamental-mode eigenfrequency is precisely zero and the associated eigenfunction is described exactly by the formula for a parabola.

Our Numerical Analysis of Bipolytropes Having (n_c, n_e) = (5, 1)

In an accompanying discussion we have shown that we can integrate the linear adiabatic wave equation (LAWE) — that is, we effectively have been able to solve the eigenvalue problem — to obtain the eigenvector associated with marginally unstable models along equilibrium sequences for bipolytropes having <math>~(n_c, n_e) = (5, 1)</math>. The marginally unstable models that have been identified via this more rigorous approach — see the orange triangles in the following figure — fall at different points along each equilibrium sequence than the points that were identified via our free-energy analysis — see the red-dashed demarcation curve. Assuming that the algorithm that we developed to integrate the LAWE was basically error-free, we trust the model identifications generated via this more rigorous technique.

Our analysis indicates that, along the <math>~\mu_e/\mu_c = 1</math> equilibrium sequence, the core-envelope interface of the marginally unstable model is located at <math>~\xi_i = 1.6686460157</math>; in the following figure, the red arrow points to this location along that sequence. In the brief table that accompanies the figure, we have listed values for the dimensionless specific entropy in the core and, separately, in the envelope, along with the interface location, <math>~\xi_i</math>.

<math>~\gamma_c = 6/5</math>     and     <math>~\gamma_e = 2</math>   Initial Model Parameters for LAWE-Determined Marginally Unstable Model <math>~\frac{\mu_e}{\mu_c}</math> <math>~\xi_i</math> <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{core}</math> <math>~\frac{s}{\Re/\bar{\mu}}\biggr|_\mathrm{env}</math> 1 1.6686460157 8.04719 1.31310 For more details, see the accompanying discussion titled, Eigenvectors from Solution of LAWE

May 11 (from Joel)

Patrick, Following our phone conversation earlier this week, I have added a number of paragraphs, derivations, and tables/figures to my VisTrails wiki page — see the above subsection titled, May 5. Here is a quick summary, assuming evolutions are conducted with (gamma_c, gamma_e) = (6/5, 2): I completely agree with your expression for the specific entropy, except I believe that the gas constant needs to be divided by the mean-molecular weight (R/mu). I concur with your s.ps plot. In particular, for the bipolytropic model having mu_e/mu_c = 1 and xi_i = 2.4161 (this is the model that my free-energy-based analysis identifies as the marginally unstable model), I get s_core = 8.04719 and s_env = 2.16080. The size of entropy jump **does** depend on the location of the interface. If we stay on the mu_e/mu_c = 1 sequence, the entropy of the envelope will equal the entropy of the core if xi_i = 12.83375. (I label this critical value, [x_i]_smooth.) At larger values of xi_i, the envelope's entropy will be larger than the core's entropy. If you choose a sequence for which mu_e/mu_c = 0.5, then my analysis says that [x_i]_smooth = 2.86620. If you set mu_e/mu_c = 1 but assume evolutions are conducted with gamma_c = gamma_e = 5/3: The entropy discontinuity disappears at the interface, but the entropy is no longer constant throughout the core or the envelope. Moving radially outward from the center, the entropy steadily increases through the core (convectively stable) then it steadily drops through the envelope (convectively unstable). Finally, we should consider the following before embarking on more 3D simulations: Along any given equilibrium sequence, a free-energy analysis only provides an **approximate** value of xi_i at which a dynamical instability sets in. The **precise** location can be identified by solving the relevant eigenvalue problem and searching for the model at which the fundamental-mode oscillation frequency goes to zero. I have taught myself how to solve this type of eigenvalue problem; in other chapters of my wiki-book, I show that I can match quantitatively related, previously published analyses. I have tentative eigenvalue-analysis results for (n_c, n_e) = (5, 1) bipolytropes; all marginally unstable models arise **sooner** along each sequence -- that is at smaller values of xi_i -- than was predicted by the free-energy analysis. I am very interested in discussing all of this with you, if you are interested and when you have the time. Just let me know. Cheers, Joel

© 2014 - 2021 by Joel E. Tohline |   H_Book Home   |   YouTube   | Appendices: | Equations | Variables | References | Ramblings | Images | myphys.lsu | ADS | Recommended citation:   Tohline, Joel E. (2021), The Structure, Stability, & Dynamics of Self-Gravitating Fluids, a (MediaWiki-based) Vistrails.org publication, https://www.vistrails.org/index.php/User:Tohline/citation