-
Notifications
You must be signed in to change notification settings - Fork 0
/
inputs.3d.cuda
193 lines (170 loc) · 8.37 KB
/
inputs.3d.cuda
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
#----------------------DOMAIN DEFINITION------------------------
geometry.is_periodic = 0 0 0 # For each dir, 0: non-perio, 1: periodic
geometry.coord_sys = 0 # 0 => cart, 1 => RZ
geometry.prob_lo = -0.00128 -0.00128 0.
geometry.prob_hi = .00128 0.00128 0.01024
# >>>>>>>>>>>>> BC FLAGS <<<<<<<<<<<<<<<<
# Interior, Inflow, Outflow, Symmetry,
# SlipWallAdiab, NoSlipWallAdiab, SlipWallIsotherm, NoSlipWallIsotherm
peleLM.lo_bc = SlipWallAdiab SlipWallAdiab Inflow
peleLM.hi_bc = SlipWallAdiab SlipWallAdiab Outflow
#-------------------------AMR CONTROL----------------------------
amr.n_cell = 64 64 256 # Level 0 number of cells in each direction
amr.v = 0 # AMR verbose
amr.max_level = 2 # maximum level number allowed
amr.ref_ratio = 2 2 2 2 # refinement ratio
amr.regrid_int = 8 # how often to regrid
amr.n_error_buf = 1 1 2 2 # number of buffer cells in error est
amr.grid_eff = 0.95 # what constitutes an efficient grid
amr.blocking_factor = 16 # block factor in grid generation (min box size)
amr.max_grid_size = 64 # max box size
#----------------------TIME STEPING CONTROL----------------------
max_step = 40000 # maximum number of time steps
stop_time = .5e-1 # final physical time
ns.cfl = 0.1 # cfl number for hyperbolic system
ns.init_dt = 5.e-6
ns.init_shrink = 1. # scale back initial timestep
ns.change_max = 1.1 # max timestep size increase
ns.dt_cutoff = 1.e-12 # level 0 timestep below which we halt
#-------------------------IO CONTROL----------------------------
amr.checkpoint_files_output = 1
amr.check_file = chk # root name of checkpoint file
amr.check_int = 100 # number of timesteps between checkpoints
amr.plot_file = plt # root name of plot file
amr.plot_int = 1 # number of timesteps between plot files
#amr.derive_plot_vars = rhoRT mag_vort avg_pressure gradpx gradpy sumRhoYdot mass_fractions molweight rhominsumrhoY cc_transport_coeffs cpmix mixfrac_only enthalpy
#amr.derive_plot_vars = mag_vort sumRhoYdot mass_fractions enthalpy
#amr.grid_log = grdlog # name of grid logging file
#----------------------- PROBLEM PARAMETERS---------------------
prob.P_mean = 6079500.
# Properties of injected particles
prob.part_temp = 363.
# Log-normal distribution of particle diameters
# Mean particle diameter
prob.part_mean_dia = 5E-6
# Standard deviation of the particle diameter
prob.part_stdev_dia = 0.
# Jet properties
prob.jet_dia = 4e-5
prob.spray_angle_deg = 21.
prob.jets_per_dir = 1
prob.jet_start_time = 0.
prob.jet_end_time = .5E-3
prob.jet_vel = 30.
prob.jet_dx_mod = 42.
#prob.mass_flow_rate = 2.0349e-5
prob.mass_flow_rate = 1.0174e-5
#--------------------SPRAY PARTICLE DATA-----------------------
peleLM.do_spray_particles = 1
particles.v = 1
particles.mom_transfer = 1
particles.mass_transfer = 1
particles.heat_transfer = 1
particles.particle_init_function = 1 # Sprays are initialized using SprayParticleInitInsert.cpp
particles.write_spray_ascii_files = 1
particles.parcel_size = 5. # Number of particles per parcel
particles.fuel_species = NC12H26
# Properties for liquid decane
particles.fuel_crit_temp = 658.2 # K
particles.fuel_boil_temp = 489. # K
particles.fuel_ref_temp = 298.15
particles.fuel_latent = 3.0299E5
particles.fuel_cp = 2.207E3 # Cp at 298
#particles.fuel_latent = 3.0294e9
particles.fuel_rho = 669.43
particles.use_splash_model = false
particles.NC12H26_psat = 4.10549 1625.928 -92.839 1.E5
#
#--------------------REFINEMENT CONTROL------------------------
amr.refinement_indicators = gradT magvort
amr.gradT.max_level = 2
amr.gradT.adjacent_difference_greater = .1
amr.gradT.field_name = temp
amr.magvort.max_level = 2
amr.magvort.vorticity_greater = 1e4
#---------------------PHYSICS CONTROL------------------------
ns.fuelName = NC12H26 # Fuel species name
ns.unity_Le = 0 # Use unity Le number transport ?
ns.hack_nochem = 0 # Bypass chemistry ? 0: no, 1: yes
ns.gravity = 0 # body force (gravity in MKS units)
ns.htt_tempmax = 19000.
peleLM.mixtureFraction.format = Cantera
peleLM.mixtureFraction.type = mole
peleLM.mixtureFraction.oxidTank = O2:0.21 N2:0.79
peleLM.mixtureFraction.fuelTank = NC12H26:1.
peleLM.chem_integrator = "ReactorCvode"
#--------------------NUMERICS CONTROL------------------------
ns.init_iter = 3 # Number of initial iterations to def pressure
ns.num_divu_iters = 1 # Number of DivU iteration at initialization
ns.sdc_iterMAX = 2 # Number of SDC iterations
ns.num_mac_sync_iter = 2 # Number of mac_sync iterations
# ------------------ INPUTS TO CVODE CHEMISTRY ----------------------
#peleLM.use_typ_vals_chem = 1 # Use species/temp typical values in CVODE
#peleLM.relative_tol_chem = 1.0e-8 # CVODE relative tolerance
#peleLM.absolute_tol_chem = 1.0e-6 # CVODE absolute tolerance
#cvode.solve_type = dense # CVODE Linear solve type (for Newton direction)
#ode.analytical_jacobian = 0 # Provide analytical jacobian (from Fuego) ?
ode.analytical_jacobian = 0
cvode.solve_type = GMRES
# ------------------ INPUTS TO ACTIVE CONTROL ----------------------
active_control.on = 0 # Use AC ?
active_control.use_temp = 1 # Default in fuel mass, rather use iso-T position ?
active_control.temperature = 800.0 # Value of iso-T ?
active_control.tau = 1.0e-7 # Control tau (should ~ 10 dt)
active_control.height = 1.000 # Where is the flame held ? Default assumes coordinate along Y in 2D or Z in 3D.
active_control.v = 0 # verbose
active_control.velMax = 50. # Optional: limit inlet velocity
active_control.changeMax = 0.2 # Optional: limit inlet velocity changes (absolute)
active_control.flameDir = 2 # Optional: flame main direction. Default: AMREX_SPACEDIM-1
# ---------------------------------------------------------------
# ------------------ ADVANCED USER INPUTS ----------------------
# ---------------------------------------------------------------
#---------------- ADV ALGORITHM CONTROL -------------------
ns.sum_interval = 5 # timesteps between computing mass
ns.do_reflux = 1 # 1 => do refluxing
ns.do_mac_proj = 1 # 1 => do MAC projection
ns.do_sync_proj = 1 # 1 => do Sync Project
ns.divu_relax_factor = 0.0
ns.be_cn_theta = 0.5
ns.do_temp = 1
ns.do_diffuse_sync = 1
ns.do_reflux_visc = 1
ns.zeroBndryVisc = 1
ns.v = 1
#
ns.divu_ceiling = 1
ns.divu_dt_factor = .4
ns.min_rho_divu_ceiling = .01
# ------------------ INPUTS TO DIFFUSION CLASS --------------------
ns.visc_tol = 1.0e-12 # tolerence for viscous solves
ns.visc_abs_tol = 1.0e-12 # tolerence for viscous solves
diffuse.max_order = 4
diffuse.tensor_max_order = 4
diffuse.v = 0
# ------------------ INPUTS TO PROJECTION CLASS -------------------
proj.proj_tol = 1.0e-12 # tolerence for projections
proj.proj_abs_tol = 1.0e-12
proj.sync_tol = 1.0e-12 # tolerence for projections
proj.rho_wgt_vel_proj = 0 # 0 => const den proj, 1 => rho weighted
proj.do_outflow_bcs = 0
proj.divu_minus_s_factor = .5
proj.divu_minus_s_factor = 0.
proj.proj_2 = 1
nodal_proj.verbose = 0
proj.v = 0
# ------------------ INPUTS TO MACPROJ CLASS -------------------
mac.mac_tol = 1.0e-12 # tolerence for mac projections
mac.mac_sync_tol = 1.0e-12 # tolerence for mac SYNC projection
mac.mac_abs_tol = 1.0e-12
mac.do_outflow_bcs = 0
mac.v = 0
#--------------------------OMP TILE INPUTS-----------------------------
#fabarray.mfiter_tile_size = 8 8 8
#--------------------------DEBUG/REGTESTS INPUTS-----------------------------
amrex.regtest_reduction=1
amrex.fpe_trap_invalid = 1
amrex.fpe_trap_zero = 1
amrex.fpe_trap_overflow = 1
#ns.fixed_dt = 1e-7
amr.restart = chk03420
amr.regrid_on_restart = 1