This PR introduces new logic to how MESA decides whether or not to use a TP-AGB wind (Blocker 1995). Previously, it used the condition s% have_done_TP which never triggers if convective shell overshooting is 0 (which it is in our inlists).
The new logic checks that the star has a He depleted core and that the He shell mass is < 0.1 Msol. The TP-AGB phase is post-He burning, and occurs via unstable He ignition in a thin He shell above the CO core. This logic is meant to identify physical conditions where this should occur, and turn the TP-AGB wind on when they are satisfied.
Additionally, the maximum of either the Reimers (RGB) or Blocker (TP-AGB) wind is used. The wind transitions from Reimers to Blocker as a low mass star evolves, and Blocker overtakes Reimers during thermal pulses.
Previously there was a condition which said that as long as the model's initial mass was greater than 10 solar masses, the star should experience high mass "Dutch" winds. This would mean that any star that may have lost mass and become "low mass" ( < 10 Msol in this case) would still have the "Dutch" wind applied. This condition has been removed. The logic now checks instead whether a star satisfies the condition
if (T1 > s% hot_wind_full_on_T) then
! assign Dutch wind
...
where s% hot_wind_full_on_T is set in the inlists as T = 12000 K.
The changes exist in run_star_extras.f.
TO DO: test that the max(Reimers, Blocker) transition works properly.
This PR introduces new logic to how MESA decides whether or not to use a TP-AGB wind (Blocker 1995). Previously, it used the condition
s% have_done_TPwhich never triggers if convective shell overshooting is 0 (which it is in our inlists).The new logic checks that the star has a He depleted core and that the He shell mass is < 0.1 Msol. The TP-AGB phase is post-He burning, and occurs via unstable He ignition in a thin He shell above the CO core. This logic is meant to identify physical conditions where this should occur, and turn the TP-AGB wind on when they are satisfied.
Additionally, the maximum of either the Reimers (RGB) or Blocker (TP-AGB) wind is used. The wind transitions from Reimers to Blocker as a low mass star evolves, and Blocker overtakes Reimers during thermal pulses.
Previously there was a condition which said that as long as the model's initial mass was greater than 10 solar masses, the star should experience high mass "Dutch" winds. This would mean that any star that may have lost mass and become "low mass" ( < 10 Msol in this case) would still have the "Dutch" wind applied. This condition has been removed. The logic now checks instead whether a star satisfies the condition
where
s% hot_wind_full_on_Tis set in the inlists as T = 12000 K.The changes exist in
run_star_extras.f.TO DO: test that the max(Reimers, Blocker) transition works properly.