[mesa-users] Convergence issues when simulating massive eruptions by modifying .mod file

[David Arnett]

MESA is mostly hydrostatic; exceptions are convection (assumed to be steady
state) and mass loss (winds assumed to be steady state). Your goal is to do
an eruption. This is not steady state but hydrodynamic (hyperbolic in the
language of differential equations).
This means that neglected (time dependent) terms are not small, and MESA
loses its solution.

I have done a version of this problem (125 Msun) by using a symplectic
radiation hydrodynamic module in TYCHO. The original hydrostatic model was
simply continued using the hydrodynamic module. It gave a "supernova
imposter", ejecting 10 Msun and brightening, sort of LBV like. The Fe bump
in opacity caused radiative acceleration which did the eruption. Bill and I
have discussed adding such a module to MESA; it is on the "to do" list.

I think Bill's suggestion is a good intermediate step.

-- 
David Arnett
Regents Professor
Steward Observatory
University of Arizona

Facts are stubborn, but statistics are more pliable.  Mark Twain
Facts do not cease to exist because they are ignored.  Aldous Huxley
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20130418/d3a43e41/attachment.html>