Context Navigation

Back to WikiStart

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
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 = 6pis^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.

Download in other formats:

Original Format