1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: bc.inc C 4 ! Purpose: Common block containing boundary conditions data C 5 ! C 6 ! Author: M. Syamlal Date: dd-mmm-yy C 7 ! Reviewer: Date: dd-mmm-yy C 8 ! C 9 ! Revision Number: C 10 ! Purpose: C 11 ! Author: Date: dd-mmm-yy C 12 ! Reviewer: Date: dd-mmm-yy C 13 ! C 14 ! Literature/Document References: None C 15 ! C 16 ! Variables referenced: None C 17 ! Variables modified: None C 18 ! C 19 ! Local variables: None C 20 ! C 21 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 22 23 24 MODULE bc 25 26 27 Use param 28 Use param1 29 30 31 ! 32 ! x coordinate of the west face of a region where 33 ! boundary conditions are specified 34 DOUBLE PRECISION BC_X_w (DIMENSION_BC) 35 ! 36 ! x coordinate of the east face of a region where 37 ! boundary conditions are specified 38 DOUBLE PRECISION BC_X_e (DIMENSION_BC) 39 ! 40 ! y coordinate of the south face of a region where 41 ! boundary conditions are specified 42 DOUBLE PRECISION BC_Y_s (DIMENSION_BC) 43 ! 44 ! y coordinate of the north face of a region where 45 ! boundary conditions are specified 46 DOUBLE PRECISION BC_Y_n (DIMENSION_BC) 47 ! 48 ! z coordinate of the bottom face of a region where 49 ! boundary conditions are specified 50 DOUBLE PRECISION BC_Z_b (DIMENSION_BC) 51 ! 52 ! z coordinate of the top face of a region where 53 ! boundary conditions are specified 54 DOUBLE PRECISION BC_Z_t (DIMENSION_BC) 55 ! 56 ! i index of the west face of a region where 57 ! boundary conditions are specified 58 INTEGER BC_I_w (DIMENSION_BC) 59 ! 60 ! i index of the east face of a region where 61 ! boundary conditions are specified 62 INTEGER BC_I_e (DIMENSION_BC) 63 ! 64 ! j index of the south face of a region where 65 ! boundary conditions are specified 66 INTEGER BC_J_s (DIMENSION_BC) 67 ! 68 ! j index of the north face of a region where 69 ! boundary conditions are specified 70 INTEGER BC_J_n (DIMENSION_BC) 71 ! 72 ! k index of the bottom face of a region where 73 ! boundary conditions are specified 74 INTEGER BC_K_b (DIMENSION_BC) 75 ! 76 ! k index of the top face of a region where 77 ! boundary conditions are specified 78 INTEGER BC_K_t (DIMENSION_BC) 79 ! 80 ! Void fraction in a specified boundary region 81 DOUBLE PRECISION BC_EP_g (DIMENSION_BC) 82 ! 83 ! Gas pressure in a specified boundary region 84 DOUBLE PRECISION BC_P_g (DIMENSION_BC) 85 ! 86 ! Microscopic density of gas in a specified 87 ! boundary region 88 DOUBLE PRECISION BC_RO_g (DIMENSION_BC) 89 ! 90 ! Macroscopic density of gas in a specified 91 ! boundary region 92 DOUBLE PRECISION BC_ROP_g (DIMENSION_BC) 93 ! 94 ! Macroscopic density of solids phases in a 95 ! specified boundary region 96 DOUBLE PRECISION BC_ROP_s (DIMENSION_BC, DIM_M) 97 DOUBLE PRECISION BC_EP_s (DIMENSION_BC, DIM_M) 98 ! 99 ! Gas phase temperature in a specified boundary 100 ! region 101 DOUBLE PRECISION BC_T_g (DIMENSION_BC) 102 ! 103 ! Solids phase temperature in a specified 104 ! boundary region 105 DOUBLE PRECISION BC_T_s (DIMENSION_BC, DIM_M) 106 ! 107 ! Solids phase granular temperature in a specified 108 ! boundary region 109 DOUBLE PRECISION BC_Theta_m (DIMENSION_BC, DIM_M) 110 ! 111 ! x-component of gas velocity in a specified 112 ! boundary region 113 DOUBLE PRECISION BC_U_g (DIMENSION_BC) 114 ! 115 ! x-component of solids phase velocity in a 116 ! specified boundary region 117 DOUBLE PRECISION BC_U_s (DIMENSION_BC, DIM_M) 118 ! 119 ! y-component of gas velocity in a specified 120 ! boundary region 121 DOUBLE PRECISION BC_V_g (DIMENSION_BC) 122 ! 123 ! y-component of solids phase velocity in a 124 ! specified boundary region 125 DOUBLE PRECISION BC_V_s (DIMENSION_BC, DIM_M) 126 ! 127 ! z-component of gas velocity in a specified 128 ! boundary region 129 DOUBLE PRECISION BC_W_g (DIMENSION_BC) 130 ! 131 ! z-component of solids phase velocity in a 132 ! specified boundary region 133 DOUBLE PRECISION BC_W_s (DIMENSION_BC, DIM_M) 134 ! 135 ! JFD: For cut cells, define the magnitude of velocity that will be enforced 136 ! perpendicular to the cut face, for CG_MI boundary condition 137 ! 138 ! magnitude of gas velocity in a specified 139 ! boundary region 140 DOUBLE PRECISION BC_VELMAG_g (DIMENSION_BC) 141 ! 142 ! magnitude of solids phase velocity in a 143 ! specified boundary region 144 DOUBLE PRECISION BC_VELMAG_s (DIMENSION_BC, DIM_M) 145 146 ! 147 ! Type of boundary: MASS_INFLOW, MASS_OUTFLOW, 148 ! P_INFLOW, P_OUTFLOW, FREE_SLIP_WALL, NO_SLIP_WALL 149 CHARACTER(LEN=16) BC_TYPE (DIMENSION_BC) 150 151 152 ! FLAG to specify if this PO BC applies to solid phase 153 ! in discrete implementation or not. For example, setting 154 ! Pressure outflow only for gas-phase. 155 156 LOGICAL BC_PO_APPLY_TO_DES (DIMENSION_BC) 157 158 ! FLAG to specify if this MI BC plane will be 159 ! seen as a wall by discrete phase (default behavior) 160 ! or allow the particles to go through the inflow plane 161 162 LOGICAL BC_MI_AS_WALL_FOR_DES (DIMENSION_BC) 163 164 ! 165 ! Gas volumetric flow rate through the boundary 166 DOUBLE PRECISION BC_VOLFLOW_g (DIMENSION_BC) 167 ! 168 ! Solids volumetric flow rate through the boundary 169 DOUBLE PRECISION BC_VOLFLOW_s (DIMENSION_BC, DIM_M) 170 ! 171 ! Gas mass flow rate through the boundary 172 DOUBLE PRECISION BC_MASSFLOW_g (DIMENSION_BC) 173 ! 174 ! Solids mass flow rate through the boundary 175 DOUBLE PRECISION BC_MASSFLOW_s (DIMENSION_BC, DIM_M) 176 ! 177 ! Logical variable to determine whether a bc is defined 178 LOGICAL BC_DEFINED (DIMENSION_BC) 179 !start kapil&anuj 01/19/98 180 ! 181 ! Logical variable to determine whether the partial 182 ! slip bc of Johnson and Jackson is used 183 INTEGER BC_JJ_PS (DIMENSION_BC) 184 !end kapil&anuj 01/19/98 185 ! 186 ! Character variable with values W, E, S, N, B, and T 187 ! to determine the flow plane of a flow cell 188 CHARACTER BC_PLANE (DIMENSION_BC) 189 ! 190 ! The interval at the beginning when normal vel. is equal to 191 ! BC_Jet_g0 192 DOUBLE PRECISION BC_DT_0 (DIMENSION_BC) 193 ! 194 ! Stored value of normal velocity 195 DOUBLE PRECISION BC_Jet_g (DIMENSION_BC) 196 ! 197 ! Value of normal vel. during the initial interval BC_DT_0 198 DOUBLE PRECISION BC_Jet_g0 (DIMENSION_BC) 199 ! 200 ! The interval when normal vel. is equal to BC_Jet_gh 201 DOUBLE PRECISION BC_DT_h (DIMENSION_BC) 202 ! 203 ! Value of normal vel. during the initial interval BC_DT_h 204 DOUBLE PRECISION BC_Jet_gh (DIMENSION_BC) 205 ! 206 ! The interval when normal vel. is equal to BC_Jet_gl 207 DOUBLE PRECISION BC_DT_l (DIMENSION_BC) 208 ! 209 ! Value of normal vel. during the initial interval BC_DT_l 210 DOUBLE PRECISION BC_Jet_gl (DIMENSION_BC) 211 ! 212 ! Time to update a transient boundary condition 213 DOUBLE PRECISION BC_TIME (DIMENSION_BC) 214 ! 215 ! Area of boundary surfaces 216 DOUBLE PRECISION BC_AREA (DIMENSION_BC) 217 ! 218 ! Volume of boundary cells 219 DOUBLE PRECISION BC_VOL (DIMENSION_BC) 220 ! 221 ! Gas species mass fractions in a boundary region 222 DOUBLE PRECISION BC_X_g (DIMENSION_BC, DIM_N_g) 223 ! 224 ! Solids species mass fractions in a boundary region 225 DOUBLE PRECISION BC_X_s (DIMENSION_BC, DIM_M, DIM_N_s) 226 ! 227 ! Accumulated or average mass outflow rate of gas 228 DOUBLE PRECISION BC_MOUT_g(DIMENSION_BC) 229 ! 230 ! Accumulated or average mass outflow rate of solids 231 DOUBLE PRECISION BC_MOUT_s(DIMENSION_BC, DIM_M) 232 ! 233 ! Accumulated or average volumetric outflow rate of gas 234 DOUBLE PRECISION BC_VOUT_g(DIMENSION_BC) 235 ! 236 ! Accumulated or average volumetric outflow rate of solids 237 DOUBLE PRECISION BC_VOUT_s(DIMENSION_BC, DIM_M) 238 ! 239 ! Number of outflow rate values accumulated 240 INTEGER BC_OUT_N (DIMENSION_BC) 241 ! 242 ! Pressure drop specified for cyclic b.c. in X 243 DOUBLE PRECISION DELP_X 244 ! 245 ! Pressure drop specified for cyclic b.c. in Y 246 DOUBLE PRECISION DELP_Y 247 ! 248 ! Pressure drop specified for cyclic b.c. in Z 249 DOUBLE PRECISION DELP_Z 250 ! 251 ! Specified mass flux (e.g., g/cm^2.s) in the cyclic 252 ! direction with specified pressure drop (only one 253 ! direction is allowed). 254 DOUBLE PRECISION Flux_g 255 ! 256 ! Average gas velocity in X direction (for cyclic bc) 257 DOUBLE PRECISION U_g0 258 ! 259 ! Average gas velocity in Y direction (for cyclic bc) 260 DOUBLE PRECISION V_g0 261 ! 262 ! Average gas velocity in Z direction (for cyclic bc) 263 DOUBLE PRECISION W_g0 264 ! 265 ! Average solids velocity in X direction (for cyclic bc) 266 DOUBLE PRECISION U_s0 (DIM_M) 267 ! 268 ! Average solids velocity in Y direction (for cyclic bc) 269 DOUBLE PRECISION V_s0 (DIM_M) 270 ! 271 ! Average solids velocity in Z direction (for cyclic bc) 272 DOUBLE PRECISION W_s0 (DIM_M) 273 ! 274 ! IJK location where P_g is fixed for cyclic b.c's 275 INTEGER IJK_P_g 276 ! 277 ! Coefficient in partial slip condition -- gas 278 DOUBLE PRECISION BC_hw_g (DIMENSION_BC) 279 ! 280 ! Coefficient in partial slip condition -- solids 281 DOUBLE PRECISION BC_hw_s (DIMENSION_BC, DIM_M) 282 ! 283 ! Wall velocity for partial slip condition -- gas 284 DOUBLE PRECISION BC_Uw_g (DIMENSION_BC) 285 ! 286 ! Wall velocity for partial slip condition -- gas 287 DOUBLE PRECISION BC_Vw_g (DIMENSION_BC) 288 ! 289 ! Wall velocity for partial slip condition -- gas 290 DOUBLE PRECISION BC_Ww_g (DIMENSION_BC) 291 ! 292 ! Wall velocity for partial slip condition -- solids 293 DOUBLE PRECISION BC_Uw_s (DIMENSION_BC, DIM_M) 294 ! 295 ! Wall velocity for partial slip condition -- solids 296 DOUBLE PRECISION BC_Vw_s (DIMENSION_BC, DIM_M) 297 ! 298 ! Wall velocity for partial slip condition -- solids 299 DOUBLE PRECISION BC_Ww_s (DIMENSION_BC, DIM_M) 300 301 ! 302 ! Coefficient in heat transfer boundary condition -- gas 303 DOUBLE PRECISION BC_hw_T_g (DIMENSION_BC) 304 ! 305 ! Coefficient in heat transfer boundary condition -- solids 306 DOUBLE PRECISION BC_hw_T_s (DIMENSION_BC, DIM_M) 307 ! 308 ! Wall temperature in heat transfer boundary condition -- gas 309 DOUBLE PRECISION BC_Tw_g (DIMENSION_BC) 310 ! 311 ! Wall temperature in heat transfer boundary condition -- solids 312 DOUBLE PRECISION BC_Tw_s (DIMENSION_BC, DIM_M) 313 ! 314 ! Coefficient in heat transfer boundary condition -- gas 315 DOUBLE PRECISION BC_C_T_g (DIMENSION_BC) 316 ! 317 ! Coefficient in heat transfer boundary condition -- solids 318 DOUBLE PRECISION BC_C_T_s (DIMENSION_BC, DIM_M) 319 320 321 ! 322 ! Coefficient in granular temp boundary condition -- solids 323 DOUBLE PRECISION BC_hw_Theta_m (DIMENSION_BC, DIM_M) 324 ! 325 ! Wall granular temperature -- solids 326 DOUBLE PRECISION BC_Thetaw_m (DIMENSION_BC, DIM_M) 327 ! 328 ! Coefficient in granualr tempearure boundary condition -- solids 329 DOUBLE PRECISION BC_C_Theta_m (DIMENSION_BC, DIM_M) 330 331 ! 332 ! Coefficient in mass transfer boundary condition -- gas 333 DOUBLE PRECISION BC_hw_X_g (DIMENSION_BC, DIM_N_g) 334 ! 335 ! Coefficient in mass transfer boundary condition -- solids 336 DOUBLE PRECISION BC_hw_X_s (DIMENSION_BC, DIM_M, DIM_N_s) 337 ! 338 ! Wall value in mass transfer boundary condition -- gas 339 DOUBLE PRECISION BC_Xw_g (DIMENSION_BC, DIM_N_g) 340 ! 341 ! Wall value in mass transfer boundary condition -- solids 342 DOUBLE PRECISION BC_Xw_s (DIMENSION_BC, DIM_M, DIM_N_s) 343 ! 344 ! Coefficient in mass transfer boundary condition -- gas 345 DOUBLE PRECISION BC_C_X_g (DIMENSION_BC, DIM_N_g) 346 ! 347 ! Coefficient in mass transfer boundary condition -- solids 348 DOUBLE PRECISION BC_C_X_s (DIMENSION_BC, DIM_M, DIM_N_s) 349 ! 350 ! 351 ! User-defined Scalars 352 ! 353 ! 354 ! scalar value in a boundary region 355 DOUBLE PRECISION BC_Scalar (DIMENSION_BC, DIM_scalar) 356 357 ! 358 ! Coefficient in boundary condition 359 DOUBLE PRECISION BC_hw_Scalar (DIMENSION_BC, DIM_scalar) 360 361 ! 362 ! Coefficient in boundary condition 363 DOUBLE PRECISION BC_C_Scalar (DIMENSION_BC, DIM_scalar) 364 365 ! 366 ! Wall value in boundary condition 367 DOUBLE PRECISION BC_ScalarW (DIMENSION_BC, DIM_scalar) 368 ! 369 ! 370 ! K and Epsilon for Gas turbulence 371 ! 372 ! 373 ! K & Epsilon values in a boundary region 374 DOUBLE PRECISION BC_K_Turb_G (DIMENSION_BC) 375 DOUBLE PRECISION BC_E_Turb_G (DIMENSION_BC) 376 377 LOGICAL:: CG_MI_CONVERTED_TO_PS(DIMENSION_BC) 378 379 380 ! Flag to specify the constant number of particles per cell 381 ! for the PIC solids 382 ! Statistical weight of parcels will be calculated by the code 383 INTEGER :: BC_PIC_MI_CONST_NPC(DIMENSION_BC, DIM_M) 384 385 ! Flag to specify the constant statistical weight. 386 ! for the PIC solids 387 ! Number of computational particles/parcels will be calculated by the code 388 DOUBLE PRECISION :: BC_PIC_MI_CONST_STATWT(DIMENSION_BC, DIM_M) 389 390 END MODULE bc 391