-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy pathcoldstart_ug.par
130 lines (92 loc) · 3.26 KB
/
coldstart_ug.par
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
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
# Runtime parameters for the Sedov explosion problem.
# Parameters for initial model
# Ambient pressure and density, explosion energy, initial
# radius of perturbation, and location of explosion center
sim_pAmbient = 1.E-5
sim_rhoAmbient = 1.
sim_expEnergy = 1.
sim_rInit = 0.013671875
sim_xctr = 0.5
sim_yctr = 0.5
sim_zctr = 0.5
sim_profFileName = "sedovSolutionModLong.dat"
# Gas ratio of specific heats
gamma = 1.4
# Computational volume parameters
# Grid dimensionality and geometry
geometry = "cartesian"
# Size of computational volume
xmin = 0.
xmax = 1.
ymin = 0.
ymax = 1.
zmin = 0.
zmax = 1.
# Boundary conditions
xl_boundary_type = "outflow"
xr_boundary_type = "outflow"
yl_boundary_type = "outflow"
yr_boundary_type = "outflow"
zl_boundary_type = "outflow"
zr_boundary_type = "outflow"
# Simulation time, I/O parameters
cfl = 0.8
basenm = "sedov_"
restart = .false.
# checkpoint file output parameters
checkpointFileIntervalTime = 0.01
checkpointFileIntervalStep = 0
checkpointFileNumber = 0
# plotfile output parameters
plotfileIntervalTime = 0.
plotfileIntervalStep = 0
plotfileNumber = 0
nend = 10000
tmax = 0.03
run_comment = "Sedov explosion"
log_file = "sedov.log"
eintSwitch = 1.e-4
plot_var_1 = "dens"
plot_var_2 = "pres"
plot_var_3 = "temp"
# Adaptive Grid refinement parameters
lrefine_max = 6
refine_var_1 = "dens"
refine_var_2 = "pres"
# Uniform Grid specific parameters
# see note below for more explanation
#iGridSize = 8 #global number of gridpoints along x, excluding gcells
#jGridSize = 8 #global number of gridpoints along y, excluding gcells
#kGridSize = 1
iProcs = 2 #num procs in i direction
jProcs = 2 #num procs in j direction
kProcs = 1
# When using UG, iProcs, jProcs and kProcs must be specified.
# These are the processors along each of the dimensions
#FIXEDBLOCKSIZE mode ::
# When using fixed blocksize, iGridSize etc are redundant in
# runtime parameters. These quantities are calculated as
# iGridSize = NXB*iprocs
# jGridSize = NYB*jprocs
# kGridSize = NZB*kprocs
#NONFIXEDBLOCKSIZE mode ::
# iGridSize etc must be specified. They constitute the global
# number of grid points in the physical domain without taking
# the guard cell into account. The local blocksize is calculated
# as iGridSize/iprocs etc.
## -------------------------------------------------------------##
## SWITCHES SPECIFIC TO THE UNSPLIT HYDRO SOLVER ##
# I. INTERPOLATION SCHEME:
order = 3 # Interpolation order (First/Second/Third order)
slopeLimiter = "mc" # Slope limiters (minmod, mc, vanLeer, hybrid, limited)
LimitedSlopeBeta= 1. # Slope parameter for the "limited" slope by Toro
charLimiting = .true. # Characteristic limiting vs. Primitive limiting
use_steepening = .true. # Contact Steepening for Third order PPM
# II. RIEMANN SOLVERS:
RiemannSolver = "hllc" # Roe, HLL, HLLC, LF
shockInstabilityFix = .false. # Carbuncle instability fix for the Roe solver
EOSforRiemann = .false. # Call EoS for Riemann solves
entropy = .false.
# III. STRONG SHOCK HANDELING SCHEME:
shockDetect = .false. # Shock Detect for numerical stability
## -------------------------------------------------------------##