-
Notifications
You must be signed in to change notification settings - Fork 2
/
Copy pathbc_mod.f
executable file
·350 lines (333 loc) · 13.9 KB
/
bc_mod.f
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
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
! C
! Module name: bc.inc C
! Purpose: Common block containing boundary conditions data C
! C
! Author: M. Syamlal Date: dd-mmm-yy C
! Reviewer: Date: dd-mmm-yy C
! C
! Revision Number: C
! Purpose: C
! Author: Date: dd-mmm-yy C
! Reviewer: Date: dd-mmm-yy C
! C
! Literature/Document References: None C
! C
! Variables referenced: None C
! Variables modified: None C
! C
! Local variables: None C
! C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
MODULE bc
Use param
Use param1
!
! x coordinate of the west face of a region where
! boundary conditions are specified
DOUBLE PRECISION BC_X_w (DIMENSION_BC)
!
! x coordinate of the east face of a region where
! boundary conditions are specified
DOUBLE PRECISION BC_X_e (DIMENSION_BC)
!
! y coordinate of the south face of a region where
! boundary conditions are specified
DOUBLE PRECISION BC_Y_s (DIMENSION_BC)
!
! y coordinate of the north face of a region where
! boundary conditions are specified
DOUBLE PRECISION BC_Y_n (DIMENSION_BC)
!
! z coordinate of the bottom face of a region where
! boundary conditions are specified
DOUBLE PRECISION BC_Z_b (DIMENSION_BC)
!
! z coordinate of the top face of a region where
! boundary conditions are specified
DOUBLE PRECISION BC_Z_t (DIMENSION_BC)
!
! i index of the west face of a region where
! boundary conditions are specified
INTEGER BC_I_w (DIMENSION_BC)
!
! i index of the east face of a region where
! boundary conditions are specified
INTEGER BC_I_e (DIMENSION_BC)
!
! j index of the south face of a region where
! boundary conditions are specified
INTEGER BC_J_s (DIMENSION_BC)
!
! j index of the north face of a region where
! boundary conditions are specified
INTEGER BC_J_n (DIMENSION_BC)
!
! k index of the bottom face of a region where
! boundary conditions are specified
INTEGER BC_K_b (DIMENSION_BC)
!
! k index of the top face of a region where
! boundary conditions are specified
INTEGER BC_K_t (DIMENSION_BC)
!
! Void fraction in a specified boundary region
DOUBLE PRECISION BC_EP_g (DIMENSION_BC)
!
! Gas pressure in a specified boundary region
DOUBLE PRECISION BC_P_g (DIMENSION_BC)
!
! Microscopic density of gas in a specified
! boundary region
DOUBLE PRECISION BC_RO_g (DIMENSION_BC)
!
! Macroscopic density of gas in a specified
! boundary region
DOUBLE PRECISION BC_ROP_g (DIMENSION_BC)
!
! Macroscopic density of solids phases in a
! specified boundary region
DOUBLE PRECISION BC_ROP_s (DIMENSION_BC, DIM_M)
!
! Gas phase temperature in a specified boundary
! region
DOUBLE PRECISION BC_T_g (DIMENSION_BC)
!
! Solids phase temperature in a specified
! boundary region
DOUBLE PRECISION BC_T_s (DIMENSION_BC, DIM_M)
!
! Solids phase granular temperature in a specified
! boundary region
DOUBLE PRECISION BC_Theta_m (DIMENSION_BC, DIM_M)
!
! x-component of gas velocity in a specified
! boundary region
DOUBLE PRECISION BC_U_g (DIMENSION_BC)
!
! x-component of solids phase velocity in a
! specified boundary region
DOUBLE PRECISION BC_U_s (DIMENSION_BC, DIM_M)
!
! y-component of gas velocity in a specified
! boundary region
DOUBLE PRECISION BC_V_g (DIMENSION_BC)
!
! y-component of solids phase velocity in a
! specified boundary region
DOUBLE PRECISION BC_V_s (DIMENSION_BC, DIM_M)
!
! z-component of gas velocity in a specified
! boundary region
DOUBLE PRECISION BC_W_g (DIMENSION_BC)
!
! z-component of solids phase velocity in a
! specified boundary region
DOUBLE PRECISION BC_W_s (DIMENSION_BC, DIM_M)
!
! Type of boundary: MASS_INFLOW, MASS_OUTFLOW,
! P_INFLOW, P_OUTFLOW, FREE_SLIP_WALL, NO_SLIP_WALL
CHARACTER*16 BC_TYPE (DIMENSION_BC)
!
! Gas volumetric flow rate through the boundary
DOUBLE PRECISION BC_VOLFLOW_g (DIMENSION_BC)
!
! Solids volumetric flow rate through the boundary
DOUBLE PRECISION BC_VOLFLOW_s (DIMENSION_BC, DIM_M)
!
! Gas mass flow rate through the boundary
DOUBLE PRECISION BC_MASSFLOW_g (DIMENSION_BC)
!
! Solids mass flow rate through the boundary
DOUBLE PRECISION BC_MASSFLOW_s (DIMENSION_BC, DIM_M)
!
! Logical variable to determine whether a bc is defined
LOGICAL BC_DEFINED (DIMENSION_BC)
!start kapil&anuj 01/19/98
!
! Logical variable to determine whether the partial
! slip bc of Johnson and Jackson is used
INTEGER BC_JJ_PS (DIMENSION_BC)
!end kapil&anuj 01/19/98
!
! Character variable with values W, E, S, N, B, and T
! to determine the flow plane of a flow cell
CHARACTER BC_PLANE (DIMENSION_BC)
!
! The interval at the beginning when normal vel. is equal to
! BC_Jet_g0
DOUBLE PRECISION BC_DT_0 (DIMENSION_BC)
!
! Stored value of normal velocity
DOUBLE PRECISION BC_Jet_g (DIMENSION_BC)
!
! Value of normal vel. during the initial interval BC_DT_0
DOUBLE PRECISION BC_Jet_g0 (DIMENSION_BC)
!
! The interval when normal vel. is equal to BC_Jet_gh
DOUBLE PRECISION BC_DT_h (DIMENSION_BC)
!
! Value of normal vel. during the initial interval BC_DT_h
DOUBLE PRECISION BC_Jet_gh (DIMENSION_BC)
!
! The interval when normal vel. is equal to BC_Jet_gl
DOUBLE PRECISION BC_DT_l (DIMENSION_BC)
!
! Value of normal vel. during the initial interval BC_DT_l
DOUBLE PRECISION BC_Jet_gl (DIMENSION_BC)
!
! Time to update a transient boundary condition
DOUBLE PRECISION BC_TIME (DIMENSION_BC)
!
! Area of boundary surfaces
DOUBLE PRECISION BC_AREA (DIMENSION_BC)
!
! Gas species mass fractions in a boundary region
DOUBLE PRECISION BC_X_g (DIMENSION_BC, DIM_N_g)
!
! Solids species mass fractions in a boundary region
DOUBLE PRECISION BC_X_s (DIMENSION_BC, DIM_M, DIM_N_s)
!
! Accumulated or average mass outflow rate of gas
DOUBLE PRECISION BC_MOUT_g(DIMENSION_BC)
!
! Accumulated or average mass outflow rate of solids
DOUBLE PRECISION BC_MOUT_s(DIMENSION_BC, DIM_M)
!
! Accumulated or average volumetric outflow rate of gas
DOUBLE PRECISION BC_VOUT_g(DIMENSION_BC)
!
! Accumulated or average volumetric outflow rate of solids
DOUBLE PRECISION BC_VOUT_s(DIMENSION_BC, DIM_M)
!
! Number of outflow rate values accumulated
INTEGER BC_OUT_N (DIMENSION_BC)
!
! Pressure drop specified for cyclic b.c. in X
DOUBLE PRECISION DELP_X
!
! Pressure drop specified for cyclic b.c. in Y
DOUBLE PRECISION DELP_Y
!
! Pressure drop specified for cyclic b.c. in Z
DOUBLE PRECISION DELP_Z
!
! Specified mass flux (e.g., g/cm^2.s) in the cyclic
! direction with specified pressure drop (only one
! direction is allowed).
DOUBLE PRECISION Flux_g
!
! Average gas velocity in X direction (for cyclic bc)
DOUBLE PRECISION U_g0
!
! Average gas velocity in Y direction (for cyclic bc)
DOUBLE PRECISION V_g0
!
! Average gas velocity in Z direction (for cyclic bc)
DOUBLE PRECISION W_g0
!
! Average solids velocity in X direction (for cyclic bc)
DOUBLE PRECISION U_s0 (DIM_M)
!
! Average solids velocity in Y direction (for cyclic bc)
DOUBLE PRECISION V_s0 (DIM_M)
!
! Average solids velocity in Z direction (for cyclic bc)
DOUBLE PRECISION W_s0 (DIM_M)
!
! IJK location where P_g is fixed for cyclic b.c's
INTEGER IJK_P_g
!
! Coefficient in partial slip condition -- gas
DOUBLE PRECISION BC_hw_g (DIMENSION_BC)
!
! Coefficient in partial slip condition -- solids
DOUBLE PRECISION BC_hw_s (DIMENSION_BC, DIM_M)
!
! Wall velocity for partial slip condition -- gas
DOUBLE PRECISION BC_Uw_g (DIMENSION_BC)
!
! Wall velocity for partial slip condition -- gas
DOUBLE PRECISION BC_Vw_g (DIMENSION_BC)
!
! Wall velocity for partial slip condition -- gas
DOUBLE PRECISION BC_Ww_g (DIMENSION_BC)
!
! Wall velocity for partial slip condition -- solids
DOUBLE PRECISION BC_Uw_s (DIMENSION_BC, DIM_M)
!
! Wall velocity for partial slip condition -- solids
DOUBLE PRECISION BC_Vw_s (DIMENSION_BC, DIM_M)
!
! Wall velocity for partial slip condition -- solids
DOUBLE PRECISION BC_Ww_s (DIMENSION_BC, DIM_M)
!
! Coefficient in heat transfer boundary condition -- gas
DOUBLE PRECISION BC_hw_T_g (DIMENSION_BC)
!
! Coefficient in heat transfer boundary condition -- solids
DOUBLE PRECISION BC_hw_T_s (DIMENSION_BC, DIM_M)
!
! Wall temperature in heat transfer boundary condition -- gas
DOUBLE PRECISION BC_Tw_g (DIMENSION_BC)
!
! Wall temperature in heat transfer boundary condition -- solids
DOUBLE PRECISION BC_Tw_s (DIMENSION_BC, DIM_M)
!
! Coefficient in heat transfer boundary condition -- gas
DOUBLE PRECISION BC_C_T_g (DIMENSION_BC)
!
! Coefficient in heat transfer boundary condition -- solids
DOUBLE PRECISION BC_C_T_s (DIMENSION_BC, DIM_M)
!
! Coefficient in granular temp boundary condition -- solids
DOUBLE PRECISION BC_hw_Theta_m (DIMENSION_BC, DIM_M)
!
! Wall granular temperature -- solids
DOUBLE PRECISION BC_Thetaw_m (DIMENSION_BC, DIM_M)
!
! Coefficient in granualr tempearure boundary condition -- solids
DOUBLE PRECISION BC_C_Theta_m (DIMENSION_BC, DIM_M)
!
! Coefficient in mass transfer boundary condition -- gas
DOUBLE PRECISION BC_hw_X_g (DIMENSION_BC, DIM_N_g)
!
! Coefficient in mass transfer boundary condition -- solids
DOUBLE PRECISION BC_hw_X_s (DIMENSION_BC, DIM_M, DIM_N_s)
!
! Wall value in mass transfer boundary condition -- gas
DOUBLE PRECISION BC_Xw_g (DIMENSION_BC, DIM_N_g)
!
! Wall value in mass transfer boundary condition -- solids
DOUBLE PRECISION BC_Xw_s (DIMENSION_BC, DIM_M, DIM_N_s)
!
! Coefficient in mass transfer boundary condition -- gas
DOUBLE PRECISION BC_C_X_g (DIMENSION_BC, DIM_N_g)
!
! Coefficient in mass transfer boundary condition -- solids
DOUBLE PRECISION BC_C_X_s (DIMENSION_BC, DIM_M, DIM_N_s)
!
!
! User-defined Scalars
!
!
! scalar value in a boundary region
DOUBLE PRECISION BC_Scalar (DIMENSION_BC, DIM_scalar)
!
! Coefficient in boundary condition
DOUBLE PRECISION BC_hw_Scalar (DIMENSION_BC, DIM_scalar)
!
! Coefficient in boundary condition
DOUBLE PRECISION BC_C_Scalar (DIMENSION_BC, DIM_scalar)
!
! Wall value in boundary condition
DOUBLE PRECISION BC_ScalarW (DIMENSION_BC, DIM_scalar)
!
!
! K and Epsilon for Gas turbulence
!
!
! K & Epsilon values in a boundary region
DOUBLE PRECISION BC_K_Turb_G (DIMENSION_BC)
DOUBLE PRECISION BC_E_Turb_G (DIMENSION_BC)
END MODULE bc