| 1 |
# Number of processors: 1536*4 |
|---|
| 2 |
##################################################################### |
|---|
| 3 |
# |
|---|
| 4 |
# FLASH file for a Type Ia supernova in 3-d cartesian geometry |
|---|
| 5 |
# |
|---|
| 6 |
# off-set ignition whole-star model without tracer particles |
|---|
| 7 |
# |
|---|
| 8 |
# This file is for the weak convergence tests of July 2008 |
|---|
| 9 |
# |
|---|
| 10 |
# This file is for a 0.2 km resolution run on Seaborg, with a centered |
|---|
| 11 |
# bubble (with varying radius) |
|---|
| 12 |
# |
|---|
| 13 |
# study evaluates convergence with the sharpened KPP ADR flame model, |
|---|
| 14 |
# also first large set of runs with quiet, realistic energy release |
|---|
| 15 |
# and neutronization |
|---|
| 16 |
# |
|---|
| 17 |
# configuration: |
|---|
| 18 |
# block of 16 |
|---|
| 19 |
# "size 16" model |
|---|
| 20 |
# domain [-2457.6,2457.6]^3 km |
|---|
| 21 |
# refine_min=1 |
|---|
| 22 |
# at min refinement, base of 3 blocks in each dimension (27 total) |
|---|
| 23 |
# |
|---|
| 24 |
# |
|---|
| 25 |
# setup with: |
|---|
| 26 |
# ./setup WD_def -3d +cube16 -maxblocks=150 +noio -auto |
|---|
| 27 |
# |
|---|
| 28 |
#----------------------------- |
|---|
| 29 |
# refine_max resolution km |
|---|
| 30 |
#----------------------------- |
|---|
| 31 |
# 1 48 102.4 |
|---|
| 32 |
# 2 96 51.2 |
|---|
| 33 |
# 3 192 25.6 |
|---|
| 34 |
# 4 384 12.8 |
|---|
| 35 |
# 5 768 6.4 |
|---|
| 36 |
# 6 1536 3.2 |
|---|
| 37 |
# 7 3072 1.6 |
|---|
| 38 |
# 8 6144 0.8 |
|---|
| 39 |
# 9 12288 0.4 |
|---|
| 40 |
# 10 24576 0.2 |
|---|
| 41 |
# 11 49152 0.1 |
|---|
| 42 |
# 12 98304 0.05 |
|---|
| 43 |
#----------------------------- |
|---|
| 44 |
# |
|---|
| 45 |
##################################################################### |
|---|
| 46 |
|
|---|
| 47 |
################## |
|---|
| 48 |
# model parameters |
|---|
| 49 |
################## |
|---|
| 50 |
|
|---|
| 51 |
# amplitude of random velocity perturbations |
|---|
| 52 |
# not used |
|---|
| 53 |
v_pert = 0.e5 |
|---|
| 54 |
n_pert = 0 |
|---|
| 55 |
|
|---|
| 56 |
# central velocity dipole size and speed |
|---|
| 57 |
# not used |
|---|
| 58 |
r_conv = 0.e5 |
|---|
| 59 |
v_conv = 0.e5 |
|---|
| 60 |
|
|---|
| 61 |
# Euler angles of the dipole in degrees |
|---|
| 62 |
# not used |
|---|
| 63 |
phi_conv = 0.e0 |
|---|
| 64 |
theta_conv = 0.e0 |
|---|
| 65 |
|
|---|
| 66 |
# flame capturing and nuclear network |
|---|
| 67 |
useFlame = .true. |
|---|
| 68 |
useBurn = .true. |
|---|
| 69 |
|
|---|
| 70 |
# ignition specific |
|---|
| 71 |
# position and size of spherical region set as burned material |
|---|
| 72 |
# at time zero |
|---|
| 73 |
ignite = .true. |
|---|
| 74 |
x_match = 0.e5 |
|---|
| 75 |
y_match = 0.e5 |
|---|
| 76 |
z_match = 0.e5 |
|---|
| 77 |
r_match = 90.5000e5 |
|---|
| 78 |
|
|---|
| 79 |
# damping |
|---|
| 80 |
# not used |
|---|
| 81 |
damp_method = 4 |
|---|
| 82 |
damp_const1 = 1.e-3 |
|---|
| 83 |
damp_const2 = 1.e-3 |
|---|
| 84 |
damp_const3 = 1.8e8 |
|---|
| 85 |
damp_time = 0.5e0 |
|---|
| 86 |
|
|---|
| 87 |
################# |
|---|
| 88 |
# flame capturing |
|---|
| 89 |
################# |
|---|
| 90 |
|
|---|
| 91 |
#eos_coulombMult = 0.0 |
|---|
| 92 |
|
|---|
| 93 |
|
|---|
| 94 |
# front tracking stuff |
|---|
| 95 |
# saves flame speed |
|---|
| 96 |
fl_debug = .true. |
|---|
| 97 |
fl_epsilon = 1.e-12 |
|---|
| 98 |
fl_epsilon_0 = 1.e-3 |
|---|
| 99 |
fl_epsilon_1 = 1.e-3 |
|---|
| 100 |
fl_kpp_fact = 1.309e0 |
|---|
| 101 |
fl_b = 3.2e0 |
|---|
| 102 |
|
|---|
| 103 |
|
|---|
| 104 |
# use an approximate atwood number for calculating |
|---|
| 105 |
# turbulent flame speed |
|---|
| 106 |
approx_atwood = .true. |
|---|
| 107 |
|
|---|
| 108 |
# turbulent flame speed parameters |
|---|
| 109 |
# turbulent flame speed is |
|---|
| 110 |
# 0.5*sqrt(A g dx* fl_subgrid_s1) |
|---|
| 111 |
fl_subgrid_s1 = 0.06e0 |
|---|
| 112 |
fl_subgrid_s2 = 0.e0 |
|---|
| 113 |
flame_speed_mult = 10.e0 |
|---|
| 114 |
pbFqdotMult = 1.e0 |
|---|
| 115 |
|
|---|
| 116 |
# flame quenching at low densities |
|---|
| 117 |
# TODO do we want to change this to 1e6-1e5? |
|---|
| 118 |
use_quenching = .true. |
|---|
| 119 |
quenching_dens0 = 5.e6 |
|---|
| 120 |
quenching_dens1 = 1.e7 |
|---|
| 121 |
|
|---|
| 122 |
################# |
|---|
| 123 |
# nuclear network |
|---|
| 124 |
################# |
|---|
| 125 |
|
|---|
| 126 |
# not used (iburn=0) |
|---|
| 127 |
# except xrenorm_burning is checked |
|---|
| 128 |
|
|---|
| 129 |
dnucmin = 1.e6 |
|---|
| 130 |
dnucmax = 1.e14 |
|---|
| 131 |
tnucmin = 1.e8 |
|---|
| 132 |
tnucmax = 1.e12 |
|---|
| 133 |
#shock_burning = .false. |
|---|
| 134 |
#xrenorm_burning = .true. |
|---|
| 135 |
#ode_steper = 1 |
|---|
| 136 |
|
|---|
| 137 |
################# |
|---|
| 138 |
# mesh definition |
|---|
| 139 |
################# |
|---|
| 140 |
|
|---|
| 141 |
# max and min levels of refinement |
|---|
| 142 |
lrefine_min = 1 |
|---|
| 143 |
lrefine_max = 10 |
|---|
| 144 |
lrefine_del = 0 |
|---|
| 145 |
|
|---|
| 146 |
# use custom refinement criteria |
|---|
| 147 |
refine_var_1 = "none" |
|---|
| 148 |
refine_var_2 = "none" |
|---|
| 149 |
refine_var_3 = "none" |
|---|
| 150 |
refine_var_4 = "none" |
|---|
| 151 |
|
|---|
| 152 |
refine_ddens = 0.5e0 |
|---|
| 153 |
refine_dtvel = 0.2e0 |
|---|
| 154 |
refine_xtvel = 100.e5 |
|---|
| 155 |
refine_dphi1 = 0.1e0 |
|---|
| 156 |
refine_xphi1 = 1.e-15 |
|---|
| 157 |
refine_sphi1 = 1.e-20 |
|---|
| 158 |
refine_xenuc = 1.e18 |
|---|
| 159 |
|
|---|
| 160 |
# minimum density at which refinement criteria other than |
|---|
| 161 |
# flame front are checked |
|---|
| 162 |
refine_dens_min = 5.e5 |
|---|
| 163 |
# max density for weighting refinement |
|---|
| 164 |
# weighted by log(dens/mindens)/log(maxdens/mindens) |
|---|
| 165 |
refine_dens_max = 1.e10 |
|---|
| 166 |
# above uni dens, don't let resolution go coarser than uni_dx |
|---|
| 167 |
refine_uni_dens = 3.e16 |
|---|
| 168 |
refine_uni_dx = 0.1e5 |
|---|
| 169 |
# within uni_radius, don't let resolution go coarser than uni_dx |
|---|
| 170 |
refine_uni_radius = 84.1000e5 |
|---|
| 171 |
# outside max_radius force resolution to be coarser than 2*uni_dx |
|---|
| 172 |
refine_max_radius = 2500.e5 |
|---|
| 173 |
|
|---|
| 174 |
# enhancement of refinement in inner region |
|---|
| 175 |
# ( help decrease number of particles per processor) |
|---|
| 176 |
refine_inner_dens_min = 2.5e10 |
|---|
| 177 |
refine_inner_dens_dx = 10.e5 |
|---|
| 178 |
|
|---|
| 179 |
refine_ignition_time = 10.e0 |
|---|
| 180 |
#The following is unused!? - KW |
|---|
| 181 |
refine_ignition_radius = 30.e5 |
|---|
| 182 |
lrefine_ignition_del = 0 |
|---|
| 183 |
|
|---|
| 184 |
#gcd refinement values |
|---|
| 185 |
|
|---|
| 186 |
gcd_focus_time = 1.5 |
|---|
| 187 |
gcd_focus_angle = 30.0 |
|---|
| 188 |
gcd_focus_min_radius = 1500.e5 |
|---|
| 189 |
gcd_focus_max_radius = 3000.e5 |
|---|
| 190 |
|
|---|
| 191 |
|
|---|
| 192 |
# physical domain |
|---|
| 193 |
geometry = "cartesian" |
|---|
| 194 |
|
|---|
| 195 |
xmin = -2457.6e5 |
|---|
| 196 |
xmax = 2457.6e5 |
|---|
| 197 |
ymin = -2457.6e5 |
|---|
| 198 |
ymax = 2457.6e5 |
|---|
| 199 |
zmin = -2457.6e5 |
|---|
| 200 |
zmax = 2457.6e5 |
|---|
| 201 |
|
|---|
| 202 |
quadrant = .false. |
|---|
| 203 |
octant = .false. |
|---|
| 204 |
|
|---|
| 205 |
nblockx = 3 |
|---|
| 206 |
nblocky = 3 |
|---|
| 207 |
nblockz = 3 |
|---|
| 208 |
|
|---|
| 209 |
# boundary conditions for the whole star |
|---|
| 210 |
xl_boundary_type = "diode" |
|---|
| 211 |
xr_boundary_type = "diode" |
|---|
| 212 |
|
|---|
| 213 |
yl_boundary_type = "diode" |
|---|
| 214 |
yr_boundary_type = "diode" |
|---|
| 215 |
|
|---|
| 216 |
zl_boundary_type = "diode" |
|---|
| 217 |
zr_boundary_type = "diode" |
|---|
| 218 |
|
|---|
| 219 |
#bc_velocity_type = "diode" |
|---|
| 220 |
|
|---|
| 221 |
grav_boundary_type = "isolated" |
|---|
| 222 |
|
|---|
| 223 |
######### |
|---|
| 224 |
# physics |
|---|
| 225 |
######### |
|---|
| 226 |
|
|---|
| 227 |
# self-gravity |
|---|
| 228 |
##igrav = 1 |
|---|
| 229 |
mpole_lmax = 4 |
|---|
| 230 |
grav_3daxisymmetric = .true. |
|---|
| 231 |
grav_dumpMoments = .false. |
|---|
| 232 |
grav_subsamplepotential = 1 |
|---|
| 233 |
|
|---|
| 234 |
# eos |
|---|
| 235 |
eos_coulombMult = 1.e0 |
|---|
| 236 |
|
|---|
| 237 |
##### |
|---|
| 238 |
# I/O |
|---|
| 239 |
##### |
|---|
| 240 |
|
|---|
| 241 |
# names of files |
|---|
| 242 |
basenm = "conv_bgp_S5h_" |
|---|
| 243 |
|
|---|
| 244 |
# for starting a new run |
|---|
| 245 |
# remember to set enuc_tint below!!! |
|---|
| 246 |
restart = .false. |
|---|
| 247 |
|
|---|
| 248 |
# checkpointFileNumber = 1 |
|---|
| 249 |
# plotFileNumber = 1 |
|---|
| 250 |
|
|---|
| 251 |
# dump checkpoint files every trstrt seconds |
|---|
| 252 |
checkpointFileIntervalTime = 1.0E-01 |
|---|
| 253 |
checkpointFileIntervalStep = 0 |
|---|
| 254 |
|
|---|
| 255 |
# dump plot files every tplot seconds |
|---|
| 256 |
plotFileIntervalTime = 5.0e-3 |
|---|
| 257 |
|
|---|
| 258 |
# go for nend steps or tmax seconds, whichever comes first |
|---|
| 259 |
nend = 10 |
|---|
| 260 |
tmax = 0.1e0 |
|---|
| 261 |
wall_clock_checkpoint = 21600.e0 # 6 hrs |
|---|
| 262 |
#wall_clock_checkpoint = 85700.e0 # 23.81 hrs |
|---|
| 263 |
|
|---|
| 264 |
# variables for plotting |
|---|
| 265 |
# TODO do these work without the spaces in them? |
|---|
| 266 |
plot_var_1 = "dens" |
|---|
| 267 |
plot_var_2 = "temp" |
|---|
| 268 |
plot_var_10 = "rpv1" |
|---|
| 269 |
|
|---|
| 270 |
corners = .false. |
|---|
| 271 |
|
|---|
| 272 |
########## |
|---|
| 273 |
# timestep |
|---|
| 274 |
########## |
|---|
| 275 |
|
|---|
| 276 |
# CFL limit |
|---|
| 277 |
cfl = 0.80e0 |
|---|
| 278 |
cfl_ini = 0.1e0 |
|---|
| 279 |
|
|---|
| 280 |
# initial and minimum |
|---|
| 281 |
dtinit = 1.59e-5 |
|---|
| 282 |
dtmin = 1.e-20 |
|---|
| 283 |
dtmax = 8.e-4 |
|---|
| 284 |
tstep_change_factor = 1.2e0 |
|---|
| 285 |
#enuc_factor = 1.e0 |
|---|
| 286 |
|
|---|
| 287 |
# density, temperature and composition of the fluff |
|---|
| 288 |
dens_fluff = 1.e-3 |
|---|
| 289 |
temp_fluff = 3.e7 |
|---|
| 290 |
xc12_fluff = 0.5e0 |
|---|
| 291 |
xo16_fluff = 0.5e0 |
|---|
| 292 |
|
|---|
| 293 |
################## |
|---|
| 294 |
# tracer particles |
|---|
| 295 |
################## |
|---|
| 296 |
|
|---|
| 297 |
# ParticleFileNumber = 1 |
|---|
| 298 |
|
|---|
| 299 |
##ParticleFileIntervalTime = 5.0e-3 |
|---|
| 300 |
##ParticleFileIntervalStep = 0 |
|---|
| 301 |
|
|---|
| 302 |
useParticles = .false. |
|---|
| 303 |
pt_numParticlesWanted = 0 |
|---|
| 304 |
pt_maxPerProc = 0 |
|---|
| 305 |
|
|---|
| 306 |
############# |
|---|
| 307 |
# miscellaneous |
|---|
| 308 |
############# |
|---|
| 309 |
|
|---|
| 310 |
# floor values |
|---|
| 311 |
|
|---|
| 312 |
small = 1.e-100 |
|---|
| 313 |
smalle = 1.e-100 |
|---|
| 314 |
smallt = 1.e7 |
|---|
| 315 |
smallu = 1.e-10 |
|---|
| 316 |
# set to 1/5 of dens_fluff |
|---|
| 317 |
smlrho = 2.e-4 |
|---|
| 318 |
smallp = 1.e15 |
|---|
| 319 |
smallx = 1.e-100 |
|---|
| 320 |
|
|---|
| 321 |
# misc (many are same as defaults) |
|---|
| 322 |
|
|---|
| 323 |
eintSwitch = 0.e0 |
|---|
| 324 |
cvisc = 0.1e0 |
|---|
| 325 |
iplm = 0 |
|---|
| 326 |
use_steepening = .false. |
|---|
| 327 |
ppm_modifystates = .false. |
|---|
| 328 |
leveque = .false. |
|---|
| 329 |
hybrid_riemann = .true. |
|---|
| 330 |
convertToConsvdForMeshCalls = .false. |
|---|
| 331 |
converttoconsvdinmeshinterp = .true. |
|---|
| 332 |
use_cma_advection = .false. |
|---|
| 333 |
dp_sh_md = 1.e0 |
|---|
| 334 |
dp_sh = 1.e0 |
|---|
| 335 |
|
|---|
| 336 |
enuc_tint = 0.e48 # this is a integral quantities thing. |
|---|
| 337 |
|
|---|
| 338 |
#checkpoint_timers = .false. |
|---|
| 339 |
|
|---|
| 340 |
memory_stat_freq = 1 |
|---|
| 341 |
run_number = 2 |
|---|
