Setup Parameters

CCSN Code

Parameters in the project/1dmlmix/setup

setup

<Input File>
Data
<Output File>
DataOut
<Input PyTorch Model>
None
<Input Grid Size to interpolate for PyTorch Model>
200
<Dump # to read>
0
<Dump time intervals [s]>
5.d-4
<Max time [s]>
0.7
<Artifiicqial Viscocity #1>
1.0
<Artificial Viscocity #2>
1.0
<External Force (not default)>
1
<EOS option>
5
<EOS Table Path>
Hempel_SFHoEOS_rho222_temp180_ye60_version_1.3_20190605.h5
<If > 1, include core mass>
2.d0
<delp (not default)>
1.d-5
<Number of steps per printout>
100
<Damping Term>
0.01
<Damping zones below this number>
1
<Flux limiter option>
1
<Capture rate option>
1
<Changing the nuclear potential energy (do not touch!)>
1.0
<yefact (not used)>
1.0
<Constant Pturb>
0.0

Extra options

Input PyTorch Model

Name	Description
None	sets turbulent pressure to 0 and doesn't use any ML subroutines
ModelName.pt	uses ML subroutines (model name can be anything)
Constant	uses subroutine turbpress_contant to apply <Constant Pturb> defined to setup. Two mode ['mach', 'rhov2'] are available that need to be specified in subroutine turbpress_constant. Recompilation is required.

EOS options

EOS	Subroutine	Description
1	eospg	perfect gas equation of state
2	eospgr	equation of state that includes gas and radiation pressure
3	eos3	Ocean eos assuming NSE for low density and Swesty-Lattimer eos for high density above rhoswe. Uses energy as the primary variable.
4	eos3	Ocean eos assuming NSE for low density and Swesty-Lattimer eos for high density above rhoswe. Uses entropy as the primary variable.
5	eos5	SFHo EOS Tables with entropy as a primary variable.

Data Preparation

Edit prep_data/setup_prep for initial parameters, then run make data from root directory.

setup_prep

<Input File>
sukhbold2016/s19.0_presn

<Output File>
Data

<Goal Size of the PNS>
600

<Goal Size of Convective Grid>
7000

<Goal Total Resolution>
8000

<Enclosed Mass Cutoff for the PNS (Msol)>
1.1

<Enclosed Mass Cutoff for Convective Region (Msol)>
1.7

<Initial Cell Mass>
4e-4

<Cell Mass Growth Rate past enclosed mass cutoff>
1e-3

<Maximum Radius of the Grid (cm)>
1500000000.0

<EOS option>
5

<EOS Table Path>
../project/1dmlmix/Hempel_SFHoEOS_rho222_temp180_ye60_version_1.3_20190605.h5

The grid in COLLAPSO1D is static and non-uniform, set by mass in each cell (deltam). It is split into 3 different regions:

Region	Resolution	Behavior
PNS	Medium	linearly decreasing deltam
Convection	High	constant
Outer	Low	exponentially growing

For the PNS and Outer regions, growth rates are determined automatically within the code, based on the target grid cells. Convection region is defined as enclosed mass cutoff for (Convective Region - PNS)

The parameters of most importance are as follows:

Parameter
Goal Size of the PNS
Goal Size of Convective Grid
Goal Total Resolution
Enclosed Mass Cutoff for the PNS (Msol)
Enclosed Mass Cutoff for Convective Region (Msol)
Maximum Radius of the Grid (cm)

Convective Region

Mass in the Convective Region is defined as Convective Region cutoff minus PNS cutoff, i.e., lines 5 - 4 from the table above. Even though it is quite thin at any given moment, a lot of mass will be infalling. Since the grid is static Lagrangian, the high-resolution convective region has to account for all of the mass that will eventually pass through it, e.g. 0.15-0.3 \(M_{\odot}\). This way, the region and the shock will be spatially resolved throughout the simulation runtime.

Outdated

After setting the target grid size and the convection region resolution, the last thing we need to actually vary is the enclosed mass of the convection region, since this is a Lagrangian code. We wanted target a convective region resolution of ~1000 to cover the region between 20 and 200 km, so here are the ecnlosed masses:

conv_region=1000

Progenitor Mass	Enclosed Mass
9	1.
10	1.
11	1.31
12	1.31
13	1.33
14	1.35
15	1.35
16	1.37
17	1.37
18	1.35
19	1.37
20	1.39
21	1.
22	1.
23	1.
24	1.
25	1.

conv_region=2000, pns=800, grid=4000

Progenitor Mass	Enclosed Mass
11	1.4
12	1.42
13	1.46
14	1.46
15	1.46
16	1.46
17	1.47
18	1.48
19	1.49
20	1.5

pns_cutoff = Enclosed Mass-0.16 for the 4k grid above. Maxrad was 5e9 cm for those runs.

Note: all of the above checkout ~200km radius ~10ms after bounce

Parameters for near ~4k resolution based on Sukhbold 2016 data. Maximum radius is 1e10 cm, i.e. 1e10 in code units.

OUTDATED: maxrad=1e1

Progenitor Mass	Initial Cell Mass	Result Grid Size
9	1.9e-5	3985
10	2.5e-5	3960
11	3.1e-5	3955
12	4e-5	3975
13	4.75e-5	3958
14	5.7e-5	3982
15	6.9e-5	3980
16	7.5e-5	3962
17	8e-5	3982
18	9.3e-5	3958
19	9.2e-5	3966
20	1.05e-4	3955
21	1.15e-4	4000
22	1.2e-4	3955
23	1.31e-4	3984
24	1.43e-4	3991
25	1.6e-4	3943

OUTDATED: variable maxrad

Progenitor Mass	Initial Cell Mass	maxrad	Result Grid Size
9	1.9e-5	1e1	3985
10	0.000063	3e4	3989
12	0.000085	6.5e4	3991
13	0.00009	6.5e4	3965
15	0.00012	6.5e4	3976
17	0.00014	6.5e4	3921
19	0.00015	6.5e4	3956

Binary to Readable

Parameters in the project/1dmlmix/setup_read

setup_readout

<Data File Name>
DataOut

<Output Basename>
DataOut_readable

<Number of dumps>
10000