mesa-users] Massive star runs end before planetary nebula stage
Bill Paxton
paxton at kitp.ucsb.edu
Thu Jul 1 13:36:06 EDT 2010
One more thing --- the particular failure mode for the run you sent is related to the convergence problem.
During an attempt to find a solution, the newton iteration has gotten horribly confused and has created a
cell with logT = 3.0 and logRho = 228.3 -- that's a rather large density! The eos not unreasonably gave
back junk leading to a NaN. Rather than doing a stop (which is what happened), the code was supposed
to reject the candidate solution and try again. This is fixed in the newest version (2483). But as I said
in the previous email, that's just a symptom of a deeper problem that remains unsolved.
-B
On Jul 1, 2010, at 10:20 AM, Eric Blais wrote:
> Hi,
>
> I'm having trouble figuring out why the more high mass stars (5 solar masses and up) end their lives in the thermally-pulsing AGB. From the logs, I don't see why the run ends at this point, and toying with some parameters has not solved this problem. Perhaps someone could help me?
>
> I've attached the inlist file. There's nothing special there, except that I changed the net to agb.net once the star went up the AGB, but the run with the basic net ended in the same way.
>
> I've included a saved model also. The error happens 194 steps later and the terminal output for the last few steps is in log5M.dat. I don't know where these variables can be found in the code, which is why I'm having trouble interpreting this error.
>
> On the computer I use, the star is halfway up the asymptotic giant branch in 15 minutes, reaches the point of maximum luminosity on this graph in about 2 hours and then spends many hours on that last, downward curve you can see on the HR diagram.
>
> It would be great if I could get these stars to go just a little bit further in their evolution (say, once the envelope has been reduced to a tenth of a solar mass).
>
> Cheers,
> Eric Blais
> <log5M.dat><inlist_test><tpagb5Msun.mod><HR_5MsunNoOversh.JPG>------------------------------------------------------------------------------
> This SF.net email is sponsored by Sprint
> What will you do first with EVO, the first 4G phone?
> Visit sprint.com/first -- http://p.sf.net/sfu/sprint-com-first_______________________________________________
> mesa-users mailing list
> mesa-users at lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/mesa-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20100701/fa34a906/attachment.html>
More information about the Mesa-users
mailing list