WikiStart: RTFlame_512_s6.par

File RTFlame_512_s6.par, 2.8 KB (added by daley, 10 years ago)

Production RTflame flash.par (16384 processes)

Line 
1# parameter file for runs to resolve the gibson scale
2
3#     ./setup RTFlame -3d -auto -nxb=16 -nyb=16 -nzb=16 +parallelio -unit=Particles/ParticlesMain
4
5# unburned and burned material
6# 50/50 CO
7ye_unburned = 0.5
8sumyi_unburned = 0.041667
9
10# to 50/50 O/Mg24
11#ye_burned = 0.5
12#sumyi_burned  =  0.072917
13#flame_deltae  = 2.8e17
14
15# to Ni56
16# gives Atwood number of 0.32 at density of 10^8
17ye_burned = 0.5
18sumyi_burned = 0.017857
19flame_deltae = 7.8587e17
20
21use_flame_speed = .true.
22flame_speed = 6.0e5
23
24# flame parameters (set width and sharpness)
25fl_epsilon_0            = 1.e-3
26fl_epsilon_1            = 1.e-3
27fl_kpp_fact             = 1.309e0
28fl_b                    = 3.2e0
29
30# initial stuff
31flame_initial_position = 80.0e5
32
33# we aren't using a velocity perturbation
34#vel_pert_amp = 10.e5
35#vel_pert_wavelength1 = 7.e5
36
37# parameters for perturbation of initial flame surface
38# one the size of the box
39spert_ampl1 = 2e5
40spert_wl1   = 15e5
41spert_phase1 = 0.0
42# one the size of lambda_fp = 6*pi*s^2/A/g
43# out of phase with above
44spert_ampl2 = 2e5
45spert_wl2   = 0.11628e5
46spert_phase2 = 0.235243
47
48dens_unburned = 1.e8
49temp_unburned = 1.e8
50gconst = -1.9e9
51
52########################
53## runtimes and dumptimes
54basenm  = "rtflame_512_s6_"
55
56tmax = 1.0
57dtinit = 1.e-9
58dtmin  = 1.e-9
59
60cfl = 0.8
61nend = 999999
62
63
64restart = .true.
65checkpointFileNumber = 129
66plotFileNumber = 250
67particleFileNumber = 250
68
69checkpointFileIntervalTime = 0.05
70wall_clock_checkpoint = 10800
71
72wall_clock_time_limit = 19800
73
74plot_var_1 = "dens"
75plot_var_2 = "temp"
76plot_var_3 = "flam"
77plot_var_4 = "velx"
78plot_var_5 = "vely"
79plot_var_6 = "velz"
80plotFileIntervalTime = 0.001
81
82particleFileIntervalTime = 0.001
83
84##########
85##  Grid
86geometry = "cartesian"
87
88nblockx = 54
89xmin = 0.e0
90xmax = 810e5
91nblocky = 1
92ymin = -7.5e5
93ymax = 7.5e5
94nblockz = 1
95zmin = -7.5e5
96zmax = 7.5e5
97
98xl_boundary_type = "reflect"
99xr_boundary_type = "diode"
100yl_boundary_type = "periodic"
101yr_boundary_type = "periodic"
102
103# here we are using uniform refine region
104# instead of adaptive refining
105#refine_var_1 = "dens"
106#refine_var_2 = "flam"
107
108# uniformly refined region near flame
109refine_uniform_region = .true.
110refine_region_size = 60e5
111refine_region_stepdown_size = 45e5
112refine_lead = 2e5
113refine_buf = 1e5
114
115# this leaves enough room to run for about 1 second
116sim_ParticleRefineRegion = .true.
117sim_ParticleRefineRegionLevel = 4
118sim_ParticleRefineRegionBottom = 50e5
119sim_ParticleRefineRegionTop = 280e5
120
121lrefine_min = 1
122lrefine_max = 6
123
124##########
125#  Particles
126useParticles = .true.
127pt_maxPerProc = 100000
128
129# this is enough particles for about 1 second
130pt_initialXMin = 80e5
131pt_initialXMax = 230e5
132pt_initialYMin = -7.5e5
133pt_initialYMax = 7.5e5
134pt_initialZMin = -7.5e5
135pt_initialZMax = 7.5e5
136pt_numX = 1280
137pt_numY = 128
138pt_numZ = 128
139
140
141
142
143#####
144# hydro
145use_steepening = .false.
146
147eachProcWritesSummary = .false.