1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: drag C 4 ! Purpose: Common block containing drag arrays C 5 ! C 6 ! Author: M. Syamlal Date: 20-MAY-92 C 7 ! C 8 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 9 10 MODULE drag 11 12 ! Gas-solids drag 13 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: F_gs 14 15 ! Solids-solids drag 16 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: F_ss 17 18 ! Off diagonal friction coefficient in HYS drag relation 19 DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: beta_ij 20 21 CONTAINS 22 23 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 24 ! C 25 ! Function(s): C_DsxRe C 26 ! Purpose: C 27 ! Calculate single sphere drag correlation multiplied by C 28 ! the Reynolds number or C 29 ! Calculate the single sphere drag correlation C 30 ! C 31 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 32 33 ! Dalla Valle (1948) 34 !----------------------------------------------------------------->>> 35 DOUBLE PRECISION FUNCTION C_DSXRE_DV(RE) 36 USE param 37 USE param1 38 IMPLICIT NONE 39 DOUBLE PRECISION, INTENT(IN) :: RE ! Reynolds number 40 41 C_DSXRE_DV = (0.63D0*SQRT(RE) + 4.8D0)**2 42 RETURN 43 END FUNCTION C_DSXRE_DV 44 45 ! Schiller and Naumann (1933) 46 !----------------------------------------------------------------->>> 47 DOUBLE PRECISION FUNCTION C_DS_SN(RE) 48 USE param 49 USE param1 50 IMPLICIT NONE 51 DOUBLE PRECISION, INTENT(IN) :: RE ! Reynolds number 52 53 C_DS_SN = 24.D0*(1.D0 + 0.15D0*RE**0.687D0)/(RE+SMALL_NUMBER) 54 RETURN 55 END FUNCTION C_DS_SN 56 !-----------------------------------------------------------------<<< 57 58 END MODULE drag 59