[mesa-users] 15Msun MESA explosion

Fri Jun 12 16:12:45 EDT 2015

I am about to play with this now, it's a lot of Fun! Thanks a lot for 
the help and info on getting this started, Bill. Will let you know
how it goes.

Manos

On 6/12/15 2:28 PM, Bill Paxton wrote:
> Hi,
>
> A quick check shows a problem in my somewhat-out-of-date inlist_adjust 
> that converts from 204 isos to 22. But I do have an old saved model 
> from when that inlist was working, so lets just skip ahead to 
> inlist_remove_core using this starting model.
>
>
>
>
> The remove-core inlist runs very quickly.
>
>
>
> The energy deposition takes longer, but with the small net it isn't 
> too long -- for me, it stops in a few minutes after 364 steps
>
>
>
> The explosion run is of course the most time consuming.  How long it 
> takes depends on how far you want to follow the shock.  Keep in mind 
> that we are not yet prepared to deal with the complexities of surface 
> shock break-out, so you'll need to stop before then.
>
>
>
> Cheers,
> Bill
>
>
>
>
>
>
> On Jun 12, 2015, at 12:01 PM, Bill Paxton wrote:
>
>> Hi Manos & Tug,
>>
>> Please note that the post core collapse runs in section 6.2 all start 
>> with a saved onset-of-core-collapse model saved when infall speed 
>> reaches 1000 km/sec.  Then there are a series of preparation inlists 
>> that remove the core and deposit energy to launch the explosive 
>> shock.  See the inlists for section 6.2 figures.
>>
>> >>>>> GAK!!!! I just discovered that the inlists for section 6.2 are 
>> not yet posted!!!!  Sorry for that oversight.
>>
>> Until that is fixed (soon!), here are some temporary replacements for 
>> you to play with if you are interested in getting started now 
>> (understanding that these may not be exactly the best way to do things).
>>
>> I have a test directory ccsn_15M in a working directory just below 
>> mesa/star.  It has these inlists and a saved starting model that uses 
>> the full 204 iso net.
>>
>> <cc_15m_204i.mod>
>>
>> <inlist>
>> run the following inlists in order by making the appropriate edits to 
>> the main inlist
>>
>> inlist_adjust -- loads the 204 iso starting model and changes it to 
>> use approx21_cr60_plus_co56.net.
>> <inlist_adjust>
>>
>> inlist_remove_core -- remove iron core at 1.5Msun
>> <inlist_remove_core>
>>
>> inlist_prepare_and_edep -- inject energy to start shock
>> enough to give total energy of 1.25e51, injected in 0.05 seconds 
>> uniformly in 0.02 Msun at base
>> <inlist_prepare_and_edep>
>>
>> inlist_explosion -- load the model after energy injection and let it run
>> <inlist_explosion>
>>
>> Let me know if you decide to try this.
>>
>> Cheers,
>> Bill
>>
>>
>>
>> On Jun 12, 2015, at 11:33 AM, Tuguldur Sukhbold wrote:
>>
>>> Manos,
>>>
>>> I am about to play with it as well, but haven't started yet. Did you 
>>> try removing the limiter? i.e.:
>>>
>>> fe_core_infall_limit = 1d99 ! cm/sec
>>>
>>> when I look at section 4 inlists, this seems to be switching to the 
>>> "hydro mode".
>>>
>>> Tuguldur
>>>
>>> On Fri, Jun 12, 2015 at 9:57 AM, Manos Chatzopoulos 
>>> <manolis at flash.uchicago.edu <mailto:manolis at flash.uchicago.edu>> wrote:
>>>
>>>     Dear MESA community,
>>>
>>>     I downloaded the inlists corresponding to Section 6 of the
>>>     instrument
>>>     paper (on the 15Msun explosion) and run them with the latest
>>>     version of MESA (7624). The only thing I changed in the inlist
>>>     was the
>>>     nuclear reaction network to approx21_cr60_plus_co56.net
>>>     <http://approx21_cr60_plus_co56.net/> in order
>>>     to get a faster run and turned the Ledoux criterion to "false"
>>>     (experimentally). The calculation run up to step number 3848
>>>     and then terminated with the message "fe_core_infall_limit" , "stop
>>>     because fe_core_infall > fe_core_infall_limit". Shouldn't with
>>>     the new
>>>     MESA Hydro mode be advanced past that phase and toward SN shock
>>>     propagation after collapse? Am I missing something? Again the
>>>     inlists
>>>     are exactly the same minus the reaction network and the convection
>>>     criterion (that I think for this experimental run shouldn't
>>>     matter in
>>>     terms of switching to hydro).
>>>
>>>     Thanks!
>>>
>>>     Manos
>>>
>>>     --
>>>     Dr. Emmanouil Chatzopoulos
>>>     Enrico Fermi Postdoctoral Fellow
>>>     Department of Astronomy & Astrophysics
>>>     FLASH Center for Computational Science
>>>     University of Chicago
>>>     https://www.linkedin.com/pub/manos-chatzopoulos/a1/688/a03
>>>     https://astro.uchicago.edu/people/manos-chatzopoulos.php
>>>
>>>
>>>     ------------------------------------------------------------------------------
>>>     _______________________________________________
>>>     mesa-users mailing list
>>>     mesa-users at lists.sourceforge.net
>>>     <mailto:mesa-users at lists.sourceforge.net>
>>>     https://lists.sourceforge.net/lists/listinfo/mesa-users
>>>
>>>     !DSPAM:10228,557b0fc5649164031471!
>>>
>>>
>>> ------------------------------------------------------------------------------
>>> _______________________________________________
>>> mesa-users mailing list
>>> mesa-users at lists.sourceforge.net 
>>> <mailto:mesa-users at lists.sourceforge.net>
>>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>>
>> ------------------------------------------------------------------------------
>> _______________________________________________
>> mesa-users mailing list
>> mesa-users at lists.sourceforge.net 
>> <mailto:mesa-users at lists.sourceforge.net>
>> https://lists.sourceforge.net/lists/listinfo/mesa-users
>

-- 
Dr. Emmanouil Chatzopoulos
Enrico Fermi Postdoctoral Fellow
Department of Astronomy & Astrophysics
FLASH Center for Computational Science
University of Chicago
https://www.linkedin.com/pub/manos-chatzopoulos/a1/688/a03
https://astro.uchicago.edu/people/manos-chatzopoulos.php

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20150612/0b27ce43/attachment.html>