| 1 | #@ total_tasks = 1104 |
|---|
| 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 2007 |
|---|
| 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 -site=seaborg.nersc.gov -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 = 80.0000e5 |
|---|
| 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 = 73.6000e5 |
|---|
| 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_pm3bS5_" |
|---|
| 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 |
|---|
![(please configure the [header_logo] section in trac.ini)](/projects/FLASH-PERI/chrome/common/trac_banner.png)
