[Mesa-users] Supernova post-shock break-out evolution with MESA

[Bill Paxton]

Hi Manos,

Good to hear from you again.  I just checked this in test case example_ccsn_IIp in the current version and after a couple of inlist tweaks it seems to work.

Replace inlist_shock_part6_header by this new version:


Then

cp inlist_shock_part6_header inlist
./rn1

The run should continue from where part5 ended.  It will need subsecond timesteps as the shock approaches the surface, then post breakout the photosphere velocity should rise quickly to about 8000 kms and the timesteps will settle in at the max allowed 1000 seconds per step.

If you want to include part6 as a standard part in the runs, edit the rn script.

The test case runs to this ending at 15 days.

        runtime (minutes), retries, backups, steps       10.87        91         2      8596
__________________________________________________________________________________________________________________________________________________

       step    lg_Tcntr    Teff     lg_LH       lg_Lnuc     Mass       H_rich     H_cntr     N_cntr     Y_surf     X_avg     eta_cntr  zones retry
  lg_dt_sec    lg_Dcntr    lg_R     lg_L3a      lg_Lneu     lg_Mdot    He_core    He_cntr    O_cntr     Z_surf     Y_avg     gam_cntr  iters bckup
   age_days    lg_Pcntr    lg_L     lg_gsurf    lg_Psurf    lg_Dsurf   C_core     C_cntr     Ne_cntr    Si_cntr    Z_avg     v_div_cs     dt_limit
__________________________________________________________________________________________________________________________________________________

       8596   4.639027   7053.675 -99.000000 -99.000000  14.455983   5.575878   0.064167   0.000721   0.362797   0.429304 -32.181741   1489     91
   2.912958 -13.276783   4.191200 -99.000000 -99.000000 -99.000000   8.880104   0.162113   0.545453   0.020353   0.382525   0.000669      3      2
 1.5000E+01   3.957852   8.729544 -99.000000  -2.010296 -13.797006   0.000000   0.019193   0.008655   0.103680  1.882E-01  0.496E+02  max increase
                               rel_E_err    9.7813287035892916D-14
                       log_rel_run_E_err       -5.5576749914832924
                           photosphere_r    1.2881719408513769D+04
                           photosphere_m       14.4392497506957405
                      photosphere_v_km_s     6639.6591976505696948
                      photosphere_cell_k            20
                                m_center        1.5780295776351911
                          total_ni_co_56    9.3329629956430135D-16
                            total_energy    7.9498235985958987D+50

Cheers,
Bill



> On Oct 9, 2019, at 10:28 AM, Manos Chatzopoulos via Mesa-users <mesa-users at lists.mesastar.org> wrote:
> 
> Dear MESA community,
> 
> I am trying to reproduce some of the results presented in the MESA IV Instrument paper that include Type II supernovae and handoffs to the STELLA code for calculation of light curves. More specifically, I am trying to reproduce Figure 31 in the instrument paper, i.e. examine the effect of handing off to STELLA at different times, including times after shock break-out.
> 
> For this reason I have downloaded all the publicly available inlists used for these test problems from mesastar.org (https://zenodo.org/record/2596697#.XZ4YaS-ZNuU). I have a successfully reproduced some of the explosions and their lightcurves, as shown in the instrument paper, but I have not determined what parameter should be changed to allow the shock to evolve past break-out in MESA, before handing off to STELLA. Any insights on that would be appreciated.
> 
> With best regards,
> 
> Manos Chatzopoulos
> 
> -- 
> Dr. Manos Chatzopoulos
> Assistant Professor
> Department of Physics & Astronomy
> Louisiana State University
> Nicholson #275
> Office # 225-578-2907
> Baton Rouge, LA 70803
> United States of America
> http://www.lsu.edu/physics/people/faculty/chatzopoulos.php
> 
> _______________________________________________
> mesa-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20191009/71822cae/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: inlist_shock_part6_header
Type: application/octet-stream
Size: 1068 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20191009/71822cae/attachment.obj>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20191009/71822cae/attachment-0001.html>