[mesa-users] Entropy Profiles

[Josiah Schwab]

Hi Kenny,

> Has anyone had any success with freezing the entropy profile of a star? I
> have tried applying an eps_grav_factor = 0 as I thought this would result
> in dS/dt = 0, but this doesn't appear to be the case. Speaking to Bill it
> sounds as though mixing is the culprit for the change in entropy, so I was
> wondering if anyone has had any success in forcing the entropy profile to
> remain the same.

If you take a star that's in thermal equilibrium and turn off composition
changes due to burning

    dxdt_nuc_factor = 0

composition changes due to mixing

    mix_factor = 0

and neutrino losses

    non_nuc_neu_factor = 0

Then I would expect that your MESA model should be happy to sit forever with

  dL/dM = eps_nuc

and hence Tds/dt = 0.

Josiah