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 ! Gas volumetric flow rate through the boundary 159 DOUBLE PRECISION BC_VOLFLOW_g (DIMENSION_BC) 160 ! 161 ! Solids volumetric flow rate through the boundary 162 DOUBLE PRECISION BC_VOLFLOW_s (DIMENSION_BC, DIM_M) 163 ! 164 ! Gas mass flow rate through the boundary 165 DOUBLE PRECISION BC_MASSFLOW_g (DIMENSION_BC) 166 ! 167 ! Solids mass flow rate through the boundary 168 DOUBLE PRECISION BC_MASSFLOW_s (DIMENSION_BC, DIM_M) 169 ! 170 ! Logical variable to determine whether a bc is defined 171 LOGICAL BC_DEFINED (DIMENSION_BC) 172 !start kapil&anuj 01/19/98 173 ! 174 ! Logical variable to determine whether the partial 175 ! slip bc of Johnson and Jackson is used 176 INTEGER BC_JJ_PS (DIMENSION_BC) 177 !end kapil&anuj 01/19/98 178 ! 179 ! Character variable with values W, E, S, N, B, and T 180 ! to determine the flow plane of a flow cell 181 CHARACTER BC_PLANE (DIMENSION_BC) 182 ! 183 ! The interval at the beginning when normal vel. is equal to 184 ! BC_Jet_g0 185 DOUBLE PRECISION BC_DT_0 (DIMENSION_BC) 186 ! 187 ! Stored value of normal velocity 188 DOUBLE PRECISION BC_Jet_g (DIMENSION_BC) 189 ! 190 ! Value of normal vel. during the initial interval BC_DT_0 191 DOUBLE PRECISION BC_Jet_g0 (DIMENSION_BC) 192 ! 193 ! The interval when normal vel. is equal to BC_Jet_gh 194 DOUBLE PRECISION BC_DT_h (DIMENSION_BC) 195 ! 196 ! Value of normal vel. during the initial interval BC_DT_h 197 DOUBLE PRECISION BC_Jet_gh (DIMENSION_BC) 198 ! 199 ! The interval when normal vel. is equal to BC_Jet_gl 200 DOUBLE PRECISION BC_DT_l (DIMENSION_BC) 201 ! 202 ! Value of normal vel. during the initial interval BC_DT_l 203 DOUBLE PRECISION BC_Jet_gl (DIMENSION_BC) 204 ! 205 ! Time to update a transient boundary condition 206 DOUBLE PRECISION BC_TIME (DIMENSION_BC) 207 ! 208 ! Area of boundary surfaces 209 DOUBLE PRECISION BC_AREA (DIMENSION_BC) 210 ! 211 ! Volume of boundary cells 212 DOUBLE PRECISION BC_VOL (DIMENSION_BC) 213 ! 214 ! Gas species mass fractions in a boundary region 215 DOUBLE PRECISION BC_X_g (DIMENSION_BC, DIM_N_g) 216 ! 217 ! Solids species mass fractions in a boundary region 218 DOUBLE PRECISION BC_X_s (DIMENSION_BC, DIM_M, DIM_N_s) 219 ! 220 ! Accumulated or average mass outflow rate of gas 221 DOUBLE PRECISION BC_MOUT_g(DIMENSION_BC) 222 ! 223 ! Accumulated or average mass outflow rate of solids 224 DOUBLE PRECISION BC_MOUT_s(DIMENSION_BC, DIM_M) 225 ! 226 ! Accumulated or average volumetric outflow rate of gas 227 DOUBLE PRECISION BC_VOUT_g(DIMENSION_BC) 228 ! 229 ! Accumulated or average volumetric outflow rate of solids 230 DOUBLE PRECISION BC_VOUT_s(DIMENSION_BC, DIM_M) 231 ! 232 ! Number of outflow rate values accumulated 233 INTEGER BC_OUT_N (DIMENSION_BC) 234 ! 235 ! Pressure drop specified for cyclic b.c. in X 236 DOUBLE PRECISION DELP_X 237 ! 238 ! Pressure drop specified for cyclic b.c. in Y 239 DOUBLE PRECISION DELP_Y 240 ! 241 ! Pressure drop specified for cyclic b.c. in Z 242 DOUBLE PRECISION DELP_Z 243 ! 244 ! Specified mass flux (e.g., g/cm^2.s) in the cyclic 245 ! direction with specified pressure drop (only one 246 ! direction is allowed). 247 DOUBLE PRECISION Flux_g 248 ! 249 ! Average gas velocity in X direction (for cyclic bc) 250 DOUBLE PRECISION U_g0 251 ! 252 ! Average gas velocity in Y direction (for cyclic bc) 253 DOUBLE PRECISION V_g0 254 ! 255 ! Average gas velocity in Z direction (for cyclic bc) 256 DOUBLE PRECISION W_g0 257 ! 258 ! Average solids velocity in X direction (for cyclic bc) 259 DOUBLE PRECISION U_s0 (DIM_M) 260 ! 261 ! Average solids velocity in Y direction (for cyclic bc) 262 DOUBLE PRECISION V_s0 (DIM_M) 263 ! 264 ! Average solids velocity in Z direction (for cyclic bc) 265 DOUBLE PRECISION W_s0 (DIM_M) 266 ! 267 ! IJK location where P_g is fixed for cyclic b.c's 268 INTEGER IJK_P_g 269 ! 270 ! Coefficient in partial slip condition -- gas 271 DOUBLE PRECISION BC_hw_g (DIMENSION_BC) 272 ! 273 ! Coefficient in partial slip condition -- solids 274 DOUBLE PRECISION BC_hw_s (DIMENSION_BC, DIM_M) 275 ! 276 ! Wall velocity for partial slip condition -- gas 277 DOUBLE PRECISION BC_Uw_g (DIMENSION_BC) 278 ! 279 ! Wall velocity for partial slip condition -- gas 280 DOUBLE PRECISION BC_Vw_g (DIMENSION_BC) 281 ! 282 ! Wall velocity for partial slip condition -- gas 283 DOUBLE PRECISION BC_Ww_g (DIMENSION_BC) 284 ! 285 ! Wall velocity for partial slip condition -- solids 286 DOUBLE PRECISION BC_Uw_s (DIMENSION_BC, DIM_M) 287 ! 288 ! Wall velocity for partial slip condition -- solids 289 DOUBLE PRECISION BC_Vw_s (DIMENSION_BC, DIM_M) 290 ! 291 ! Wall velocity for partial slip condition -- solids 292 DOUBLE PRECISION BC_Ww_s (DIMENSION_BC, DIM_M) 293 294 ! 295 ! Coefficient in heat transfer boundary condition -- gas 296 DOUBLE PRECISION BC_hw_T_g (DIMENSION_BC) 297 ! 298 ! Coefficient in heat transfer boundary condition -- solids 299 DOUBLE PRECISION BC_hw_T_s (DIMENSION_BC, DIM_M) 300 ! 301 ! Wall temperature in heat transfer boundary condition -- gas 302 DOUBLE PRECISION BC_Tw_g (DIMENSION_BC) 303 ! 304 ! Wall temperature in heat transfer boundary condition -- solids 305 DOUBLE PRECISION BC_Tw_s (DIMENSION_BC, DIM_M) 306 ! 307 ! Coefficient in heat transfer boundary condition -- gas 308 DOUBLE PRECISION BC_C_T_g (DIMENSION_BC) 309 ! 310 ! Coefficient in heat transfer boundary condition -- solids 311 DOUBLE PRECISION BC_C_T_s (DIMENSION_BC, DIM_M) 312 313 314 ! 315 ! Coefficient in granular temp boundary condition -- solids 316 DOUBLE PRECISION BC_hw_Theta_m (DIMENSION_BC, DIM_M) 317 ! 318 ! Wall granular temperature -- solids 319 DOUBLE PRECISION BC_Thetaw_m (DIMENSION_BC, DIM_M) 320 ! 321 ! Coefficient in granualr tempearure boundary condition -- solids 322 DOUBLE PRECISION BC_C_Theta_m (DIMENSION_BC, DIM_M) 323 324 ! 325 ! Coefficient in mass transfer boundary condition -- gas 326 DOUBLE PRECISION BC_hw_X_g (DIMENSION_BC, DIM_N_g) 327 ! 328 ! Coefficient in mass transfer boundary condition -- solids 329 DOUBLE PRECISION BC_hw_X_s (DIMENSION_BC, DIM_M, DIM_N_s) 330 ! 331 ! Wall value in mass transfer boundary condition -- gas 332 DOUBLE PRECISION BC_Xw_g (DIMENSION_BC, DIM_N_g) 333 ! 334 ! Wall value in mass transfer boundary condition -- solids 335 DOUBLE PRECISION BC_Xw_s (DIMENSION_BC, DIM_M, DIM_N_s) 336 ! 337 ! Coefficient in mass transfer boundary condition -- gas 338 DOUBLE PRECISION BC_C_X_g (DIMENSION_BC, DIM_N_g) 339 ! 340 ! Coefficient in mass transfer boundary condition -- solids 341 DOUBLE PRECISION BC_C_X_s (DIMENSION_BC, DIM_M, DIM_N_s) 342 ! 343 ! 344 ! User-defined Scalars 345 ! 346 ! 347 ! scalar value in a boundary region 348 DOUBLE PRECISION BC_Scalar (DIMENSION_BC, DIM_scalar) 349 350 ! 351 ! Coefficient in boundary condition 352 DOUBLE PRECISION BC_hw_Scalar (DIMENSION_BC, DIM_scalar) 353 354 ! 355 ! Coefficient in boundary condition 356 DOUBLE PRECISION BC_C_Scalar (DIMENSION_BC, DIM_scalar) 357 358 ! 359 ! Wall value in boundary condition 360 DOUBLE PRECISION BC_ScalarW (DIMENSION_BC, DIM_scalar) 361 ! 362 ! 363 ! K and Epsilon for Gas turbulence 364 ! 365 ! 366 ! K & Epsilon values in a boundary region 367 DOUBLE PRECISION BC_K_Turb_G (DIMENSION_BC) 368 DOUBLE PRECISION BC_E_Turb_G (DIMENSION_BC) 369 370 LOGICAL:: CG_MI_CONVERTED_TO_PS(DIMENSION_BC) 371 372 373 ! Flag to specify the constant number of particles per cell 374 ! for the PIC solids 375 ! Statistical weight of parcels will be calculated by the code 376 INTEGER :: BC_PIC_MI_CONST_NPC(DIMENSION_BC, DIM_M) 377 378 ! Flag to specify the constant statistical weight. 379 ! for the PIC solids 380 ! Number of computational particles/parcels will be calculated by the code 381 DOUBLE PRECISION :: BC_PIC_MI_CONST_STATWT(DIMENSION_BC, DIM_M) 382 383 END MODULE bc 384