-
Notifications
You must be signed in to change notification settings - Fork 2
Expand file tree
/
Copy pathmfix.dat
More file actions
executable file
·301 lines (219 loc) · 10.4 KB
/
Copy pathmfix.dat
File metadata and controls
executable file
·301 lines (219 loc) · 10.4 KB
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
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
!
! Generic simulation base for small eruptions and lab. setup
!
!
!
! Run-control section
RUN_NAME = 'lab'
DESCRIPTION = 'lab run'
RUN_TYPE = 'NEW'
UNITS = 'SI'
TIME = 0.0 !Start time
TSTOP = 300.0 !Stop Time
DT = 1.0E-5 !Initial time step
DT_MAX = 5.0e-2 !Max time step
ENERGY_EQ = .TRUE. !Solve Energy (Thermal)?
DT_MIN = 1.0e-9 !Minimum time step before declaring non-convergence
SPECIES_EQ(0)=.FALSE. !do not solve species eq
SPECIES_EQ(1)=.FALSE.
SPECIES_EQ(2)=.FALSE.
SPECIES_EQ(3)=.FALSE.
!MAX_NIT=50
TOL_RESID=1.0e-2 !Residuals for convergence
MAX_INLET_VEL_FAC=500.
GRAVITY=9.81
!-------------------------------------------------------------------------------------------------------!
! Controls for pulsing flow. Make sure all values are specified, and that the number of PARTICLE_INLET_FRAC
! Equals the number of MMAX. Also, ensure that the values sum to 1.
MASS_INFLUX_TIME=.TRUE. ! Circularize the vent and keep constant ratio of particles.
PULSED=.TRUE. ! False for steady, True for pulsing.
FREQUENCY=1.0
MIN_GAS=0.995
MAX_GAS=0.999
PARTICLE_INLET_FRAC(1) = 0.24
PARTICLE_INLET_FRAC(2) = 0.41
PARTICLE_INLET_FRAC(3)= 0.35
!Atmospheric parameters
ATMOSPHERIC=.TRUE.
TROPOPAUSE=11000.0 ! in meters -- make sure this specification matches temp profile in eosg.f
!-------------------------------------------------------------------------------------------------------!
CALL_USR=.TRUE.
GRANULAR_ENERGY=.FALSE.
FRICTION=.FALSE.
DISCRETIZE(1)=2 !Superbee discretization
DISCRETIZE(2)=2
DISCRETIZE(3)=2
DISCRETIZE(4)=2
DISCRETIZE(5)=2
WRITE_LEVEL=1 !Control for write-out (1 is standard ascii write out
LES=.TRUE. !LES Turbulence model
LEAK=.FALSE. !Boundary Leak model
LAGRAN=.FALSE. !Particle tracking
L_SCALE0=50.0
!LAGRAN Specifications-----------------------------------!
!----------------------------------------------------------------!
LAGRAN_UPDATE=0.1
BOT_X_LAG=23800.0 !Region to introduce Tracers
TOP_X_LAG=24200.0
BOT_Y_LAG=300.0
TOP_Y_LAG=500.0
U_LAG=0.0 !Passive tracers should be initialized with 0 Vel. at beginning.
V_LAG=0.0
TIME_START_LAG=0.0 !Time to start introducing tracers
TIME_STOP_LAG= 0.001 !Time to stop introducting tracers
DT_LAG=5000 !Number of tracers introduced per time step
LAGRANWRITEOUT=0.5 !How often information is given about tracers
!The primary output now can be found in Particles1
!This gives: 1. Time 2. Tracer number 3. Bubble diameter 4. Dissolve Water content 5. Fragmentation condition 6. X position 7. Y Position 8. U velocity 9. V velocity
LAG_RAD=0.10
LAG_DENSE=3000.0
WATER_RUN=.FALSE. !Simulation with water boundary (don't use for this model base)
Mu_gmax=1.0e4 !Maximum Eulerian (gas) phase viscosity
DT_FAC=.8 !Adaptive time stepping
DETECT_STALL=.TRUE.
! Geometry Section
COORDINATES = 'cartesian'
XLENGTH = 8000. !width
IMAX = 160 !cells in i direction
YLENGTH = 10E3 !height
JMAX = 200 !cells in j direction
NO_K = .FALSE. !Don't use K? (2D)
ZLENGTH = 8000.0
KMAX = 160
! Gas-phase Section
NMAX(0)=2 !Two gas species
MW_g = 28.47 18.00
! MU_G0=1e-5 !Constant gas viscosity (change if using energy eq.)
! Solids-phase Section
MMAX=3 !Two solid phases
NMAX(1)= 1
NMAX(2)=1
NMAX(3)=1
MW_s=48.0
MW_s=48.0
MW_s=48.0
CLOSE_PACKED(1)=.TRUE.
CLOSE_PACKED(2)=.TRUE.
CLOSE_PACKED(3)=.TRUE.
UR_FAC(1)=.1 !relaxation factors for stiff problems
UR_FAC(2)=.1
C_f = .57 !Coefficient of Friction
RO_s = 2740.0 2740.0 1950.0 !solids density
D_p0 = 3.e-5 2.5e-4 2e-3 !particle diameter
e = 0.5 !restitution coefficient
Phi = 30.0 !angle of internal friction
EP_star = 0.0 !void fraction at minimum
! fluidization
! Initial Conditions Section
! 1. initial condition 1
IC_X_w(1) = 0.0
IC_X_e(1) = 8E3
IC_Y_s(1) = 0.0
IC_Y_n(1) = 10E3
IC_Z_b(1) = 0.0
IC_Z_t(1) = 8E3
IC_EP_g = 1.0 !gas volume fraction
IC_P_Star(1) = 0.0
IC_ROP_s(1,1) = 0.0
IC_ROP_s(1,2)=0.0
IC_ROP_s(1,3)=0.0
IC_T_g(1)=300.0
IC_T_s(1,1)=300.
IC_T_s(1,2)=300.
IC_T_s(1,3)=300.0
IC_X_g(1,1) = 1.0
IC_X_s(1,1,1) = 1.0
IC_U_g(1) = 0.0 !radial gas velocity
IC_V_g(1) = 0.0 !axial gas velocity
IC_W_g(1) = 0.0
IC_U_s(1,1) = 0.0 !radial solids velocity
IC_V_s(1,1) = 0.0 !axial solids velocity
IC_U_s(1,2) = 0.0 !radial solids velocity
IC_V_s(1,2) = 0.0 !axial solids velocity
IC_W_s(1,1) = 0.0
IC_W_s(1,2) = 0.0
IC_U_s(1,3)=0.0
IC_V_s(1,3)=0.0
IC_W_s(1,3)=0.0
IC_Theta_m(1,1)=30.0
IC_Theta_m(1,2)=30.0
IC_Theta_m(1,3)=30.0
! Boundary Conditions Section
!2. Inflow
SETRADIUS=200.0
BC_X_w(2)=3800.
BC_X_e(2)=4200.
BC_Y_s(2)=0.0
BC_Y_n(2)=0.0
BC_Z_b(2)=3800.0
BC_Z_t(2)=4200.
BC_T_g(2)=1100.0
BC_T_s(2,1)=1100.0
BC_T_s(2,2)=1100.0
BC_T_s(2,3)=1100.0
BC_TYPE(2)= 'MASS_INFLOW' !Specified mass inflow
BC_EP_g(2) = 1.0 ! gas volume fraction
BC_ROP_s(2,1) = 0.0 !PARTICLE_INLET_FRAC(1)*RO_s(1)*(1.-BC_EP_g(2)) ! particle volume fraction * effective density
BC_ROP_s(2,2) = 0.0 !PARTICLE_INLET_FRAC(2)*RO_s(2)*(1.-BC_EP_g(2))
BC_ROP_S(2,3) = 0.0 !PARTICLE_INLET_FRAC(3)*RO_s(3)*(1.-BC_EP_g(2))
BC_X_g(2,1) = 1.0
BC_X_s(2,1,1) = 1.0
BC_X_s(2,2,1) = 1.0
BC_X_s(2,3,1)=1.0
BC_Xw_g(2,1)=1.0
BC_Xw_g(2,2)=0.0
BC_U_g(2) = 0.0 !radial gas velocity
BC_V_g(2) = 0.0 !axial gas velocity
BC_W_g(2) = 0.0
BC_U_s(2,1) = 0.0!50.0 !radial solids velocity
BC_V_s(2,1) = 0.0 !axial solids velocity
BC_W_s(2,1) = 0.0
BC_U_s(2,2) = 0.0 !radial solids velocity
BC_V_s(2,2) = 0.0 !axial solids velocity
BC_W_s(2,2) = 0.0
BC_U_s(2,3) = 0.0 !radial solids velocity
BC_V_s(2,3) = 0.0 !axial solids velocity
BC_W_s(2,3) = 0.0
BC_Theta_m(2,1) = 100 ! Granular temperature
BC_Theta_m(2,2) = 100
BC_Theta_m(2,3)=100.
BC_P_g(2)=1.2e5 !Inlet pressure
! 3. Exit
BC_X_w(3) = 0.0 !top exit
BC_X_e(3) = 8E3
BC_Y_s(3) = 10E3
BC_Y_n(3) = 10E3
BC_Z_b(3) = 0.0
BC_Z_t(3) = 8E3
BC_TYPE(3) = 'P_OUTFLOW' !specified pressure outflow
BC_P_g(3) = 2.64e4
!
! Output Control
!
OUT_DT = 5.0 !write text file
TIME_INTERVAL = 5.0 !same as OUT_DT, used in mfixconst for post-processing -- T. Black, 10 July 2015
!
RES_DT = 200. !write binary restart file
!
NLOG = 1e30 !write logfile every NLOG time steps
FULL_LOG = .TRUE. !display residuals on screen
!SPX_DT values determine how often SPx files are written.
!
! EP_g P_g U_g U_s ROP_s T_g X_g
! P_star V_g V_s T_s X_s Theta Scalar
! W_g W_s
SPX_DT = 100.1 100.1 100.1 100.1 100.1 100.1 100.1 100.1 100.1
!Parallel controls (domain decomposition)
NODESI=4
NODESJ=4
NODESK=4
! Sweep Direction
LEQ_SWEEP(1) = 'ISIS'
LEQ_SWEEP(2) = 'ISIS'
LEQ_SWEEP(3) = 'ISIS'
LEQ_SWEEP(4) = 'ISIS'
LEQ_SWEEP(5) = 'ISIS'
LEQ_SWEEP(6) = 'ISIS'
LEQ_SWEEP(7) = 'ISIS'
LEQ_SWEEP(8) = 'ISIS'
LEQ_SWEEP(9) = 'ISIS'