Phoebus-LANL-CCSN-GPU-code/inputs/thincooling.pin at main · tomekplewa/Phoebus-LANL-CCSN-GPU-code · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
<phoebus>
problem = thin_cooling

<parthenon/job>
problem_id  = cooling       # problem ID: basename of output filenames

<parthenon/output1>
variables = p.density,  &
	          c.density,  &
            p.velocity, &
            p.energy,   &
            p.ye,       &
            temperature

file_type   = hdf5      # Tabular data dump
dt          = 1.e-3     # time increment between outputs

<parthenon/time>
nlim        = -1       # cycle limit
tlim        = 1.e-1           # time limit
integrator  = rk2      # time integration algorithm
ncycle_out  = 1        # interval for stdout summary info
dt_init_fact = 1.e-6

<parthenon/mesh>
nghost = 4

nx1         = 4         # Number of zones in X1-direction
x1min       = 0         # minimum value of X1
x1max       = 5.e-4     # maximum value of X1 (for monte_carlo)
ix1_bc      = periodic  # Inner-X1 boundary condition flag
ox1_bc      = periodic  # Outer-X1 boundary condition flag

nx2         = 1          # Number of zones in X2-direction
x2min       = -1         # minimum value of X2
x2max       = 1          # maximum value of X2
ix2_bc      = periodic   # Inner-X2 boundary condition flag
ox2_bc      = periodic   # Outer-X2 boundary condition flag

nx3         = 1         # Number of zones in X3-direction
x3min       = -1        # minimum value of X3
x3max       = 1         # maximum value of X3
ix3_bc      = periodic  # Outer-X3 boundary condition flfgag
ox3_bc      = periodic  # Outer-X3 boundary condition flfgag

num_threads = 1         # maximum number of OMP threads

<parthenon/meshblock>
nx1 = 4
nx2 = 1
nx3 = 1

<eos>
type = IdealGas
Gamma = 1.6666666666666666667
Cv = 1.995640e+03

<units>
scale_free = false
geom_length_cm = 1.e14
fluid_mass_g = 1.e0

<physics>
hydro = true
he = false
3t = false
rad = true
mhd = false
tracers = false

<fluid>
xorder = 2
cfl = 0.9
riemann = hll
recon = linear
c2p_max_iter = 100
Ye = true

<radiation>
#method = cooling_function
method = monte_carlo
nu_min = 1.e15
nu_max = 0.9999e17
nu_bins = 100
tune_emission = 1.e-43
#tune_emission = 1.e-41
#tune_emission = 1.e-34
num_particles = 16
opacity_model = tophat
absorption = false
do_nu_electron = true
do_nu_electron_anti = false
do_nu_heavy = false

<opacity>
type = tophat
tophat_C = 1.e0
tophat_numin = 1.e15
tophat_numax = 1.e17

<mean_opacity>
numin = 1.e15
numax = 1.e17
nnu = 10

<monte_carlo>
remove_emitted_particles = true

<thincooling>
species = electron