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 |
---|
7 | ye_unburned = 0.5 |
---|
8 | sumyi_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 |
---|
17 | ye_burned = 0.5 |
---|
18 | sumyi_burned = 0.017857 |
---|
19 | flame_deltae = 7.8587e17 |
---|
20 | |
---|
21 | use_flame_speed = .true. |
---|
22 | flame_speed = 6.0e5 |
---|
23 | |
---|
24 | # flame parameters (set width and sharpness) |
---|
25 | fl_epsilon_0 = 1.e-3 |
---|
26 | fl_epsilon_1 = 1.e-3 |
---|
27 | fl_kpp_fact = 1.309e0 |
---|
28 | fl_b = 3.2e0 |
---|
29 | |
---|
30 | # initial stuff |
---|
31 | flame_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 |
---|
39 | spert_ampl1 = 2e5 |
---|
40 | spert_wl1 = 15e5 |
---|
41 | spert_phase1 = 0.0 |
---|
42 | # one the size of lambda_fp = 6*pi*s^2/A/g |
---|
43 | # out of phase with above |
---|
44 | spert_ampl2 = 2e5 |
---|
45 | spert_wl2 = 0.11628e5 |
---|
46 | spert_phase2 = 0.235243 |
---|
47 | |
---|
48 | dens_unburned = 1.e8 |
---|
49 | temp_unburned = 1.e8 |
---|
50 | gconst = -1.9e9 |
---|
51 | |
---|
52 | ######################## |
---|
53 | ## runtimes and dumptimes |
---|
54 | basenm = "rtflame_512_s6_" |
---|
55 | |
---|
56 | tmax = 1.0 |
---|
57 | dtinit = 1.e-9 |
---|
58 | dtmin = 1.e-9 |
---|
59 | |
---|
60 | cfl = 0.8 |
---|
61 | nend = 999999 |
---|
62 | |
---|
63 | |
---|
64 | restart = .true. |
---|
65 | checkpointFileNumber = 129 |
---|
66 | plotFileNumber = 250 |
---|
67 | particleFileNumber = 250 |
---|
68 | |
---|
69 | checkpointFileIntervalTime = 0.05 |
---|
70 | wall_clock_checkpoint = 10800 |
---|
71 | |
---|
72 | wall_clock_time_limit = 19800 |
---|
73 | |
---|
74 | plot_var_1 = "dens" |
---|
75 | plot_var_2 = "temp" |
---|
76 | plot_var_3 = "flam" |
---|
77 | plot_var_4 = "velx" |
---|
78 | plot_var_5 = "vely" |
---|
79 | plot_var_6 = "velz" |
---|
80 | plotFileIntervalTime = 0.001 |
---|
81 | |
---|
82 | particleFileIntervalTime = 0.001 |
---|
83 | |
---|
84 | ########## |
---|
85 | ## Grid |
---|
86 | geometry = "cartesian" |
---|
87 | |
---|
88 | nblockx = 54 |
---|
89 | xmin = 0.e0 |
---|
90 | xmax = 810e5 |
---|
91 | nblocky = 1 |
---|
92 | ymin = -7.5e5 |
---|
93 | ymax = 7.5e5 |
---|
94 | nblockz = 1 |
---|
95 | zmin = -7.5e5 |
---|
96 | zmax = 7.5e5 |
---|
97 | |
---|
98 | xl_boundary_type = "reflect" |
---|
99 | xr_boundary_type = "diode" |
---|
100 | yl_boundary_type = "periodic" |
---|
101 | yr_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 |
---|
109 | refine_uniform_region = .true. |
---|
110 | refine_region_size = 60e5 |
---|
111 | refine_region_stepdown_size = 45e5 |
---|
112 | refine_lead = 2e5 |
---|
113 | refine_buf = 1e5 |
---|
114 | |
---|
115 | # this leaves enough room to run for about 1 second |
---|
116 | sim_ParticleRefineRegion = .true. |
---|
117 | sim_ParticleRefineRegionLevel = 4 |
---|
118 | sim_ParticleRefineRegionBottom = 50e5 |
---|
119 | sim_ParticleRefineRegionTop = 280e5 |
---|
120 | |
---|
121 | lrefine_min = 1 |
---|
122 | lrefine_max = 6 |
---|
123 | |
---|
124 | ########## |
---|
125 | # Particles |
---|
126 | useParticles = .true. |
---|
127 | pt_maxPerProc = 100000 |
---|
128 | |
---|
129 | # this is enough particles for about 1 second |
---|
130 | pt_initialXMin = 80e5 |
---|
131 | pt_initialXMax = 230e5 |
---|
132 | pt_initialYMin = -7.5e5 |
---|
133 | pt_initialYMax = 7.5e5 |
---|
134 | pt_initialZMin = -7.5e5 |
---|
135 | pt_initialZMax = 7.5e5 |
---|
136 | pt_numX = 1280 |
---|
137 | pt_numY = 128 |
---|
138 | pt_numZ = 128 |
---|
139 | |
---|
140 | |
---|
141 | |
---|
142 | |
---|
143 | ##### |
---|
144 | # hydro |
---|
145 | use_steepening = .false. |
---|
146 | |
---|
147 | eachProcWritesSummary = .false. |
---|