Phoebus-LANL-CCSN-GPU-code/inputs/blast_wave.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
# © 2021-2025. Triad National Security, LLC. All rights reserved.  This
# program was produced under U.S. Government contract 89233218CNA000001
# for Los Alamos National Laboratory (LANL), which is operated by Triad
# National Security, LLC for the U.S.  Department of Energy/National
# Nuclear Security Administration. All rights in the program are
# reserved by Triad National Security, LLC, and the U.S. Department of
# Energy/National Nuclear Security Administration. The Government is
# granted for itself and others acting on its behalf a nonexclusive,
# paid-up, irrevocable worldwide license in this material to reproduce,
# prepare derivative works, distribute copies to the public, perform
# publicly and display publicly, and to permit others to do so.

# Non-relativistic Sedov blast wave
# As descriged in in the Castro test suite
# https://amrex-astro.github.io/Castro/docs/Verification.html
# To do a non-relativistic blast, set rescale = small number
# To do a relativistic one, set rescale = 1

<phoebus>
problem = sedov

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

<parthenon/output1>
variables = p.density,  &
      	    c.density,  &
            p.velocity, &
            c.momentum, &
            p.energy,   &
            c.energy,   &
	          pressure,	&
      	    cs, &
      	    p.ye, &
      	    g.c.coord, &
      	    g.n.coord

file_type   = hdf5      # Tabular data dump
dt          = 0.005     # time increment between outputs

<parthenon/time>
nlim        = -1        # cycle limit
tlim        = 0.5       # 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
#refinement = adaptive
#numlevel = 3

nx1         = 512         # Number of zones in X1-direction
x1min       = 0.0         # minimum value of X1
x1max       = 0.3         # maximum value of X1
ix1_bc      = reflect  # Inner-X1 boundary condition flag
ox1_bc      = gr_outflow  # Outer-X1 boundary condition flag

nx2         = 1             # Number of zones in X2-direction
x2min       = 0             # minimum value of X2
x2max       = 3.14159265359 # maximum value of X2. Pi
ix2_bc      = reflect       # Inner-X2 boundary condition flag
ox2_bc      = reflect       # Outer-X2 boundary condition flag

nx3         = 1             # Number of zones in X3-direction
x3min       = 0             # minimum value of X3
x3max       = 6.28318530718 # maximum value of X3. 2*pi
ix3_bc      = periodic      # Inner-X3 boundary condition flag
ox3_bc      = periodic      # Outer-X3 boundary condition flfgag

num_threads = 1         # maximum number of OMP threads

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

<parthenon/refinement0>
field = c.c.bulk.rho
method = derivative_order_1
max_level = 3

<eos>
type = IdealGas
Gamma = 1.6666666666666666667
Cv = 1.0

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

<fluid>
xorder = 2
cfl = 0.1
riemann = hll
recon = mp5
c2p_max_iter = 1000
Ye = false

<sedov>
rho_ambient = 1
p_ambient = 1e-3
explosion_energy = 1
rinner = 0.05
#spherical_coords = false