[Mesa-users] Turning off hydrogen recombination energy

[Natalie Rees]

Dear mailing list,

I am trying to prevent the hydrogen recombination instability in TP-AGB Stars (Wagenhuber & Weiss 1994, https://ui.adsabs.harvard.edu/abs/1994A%26A...290..807W/abstract). To do so I want to try removing the energy produced by hydrogen recombination in the outer envelope temporarily. Does anyone have any guidance to as to how I could achieve this? I’ve tried playing around with the other_energy hook a bit but even if I inject a small amount of negative heat (small compared to the current energy of the shell) it causes convergence difficulties. Any advice or suggestions would be greatly appreciated.

Many thanks,
Natalie Rees



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230727/5de63e9d/attachment.htm>