drag_mod.f

Go to the documentation of this file.

! -*- f90 -*-
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Module name: drag                                                   C
!  Purpose: Common block containing drag arrays                        C
!                                                                      C
!  Author: M. Syamlal                                 Date: 20-MAY-92  C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

MODULE drag

! Gas-solids drag
  DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  f_gs

! Solids-solids drag
  DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  f_ss

! Off diagonal friction coefficient in HYS drag relation
  DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE ::  beta_ij

CONTAINS

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Function(s): C_DsxRe                                                C
!  Purpose:                                                            C
!     Calculate single sphere drag correlation multiplied by           C
!     the Reynolds number or                                           C
!     Calculate the single sphere drag correlation                     C
!                                                                      C
!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC

! Dalla Valle (1948)
!----------------------------------------------------------------->>>
  DOUBLE PRECISION FUNCTION c_dsxre_dv(RE)
    USE param
    USE param1
    IMPLICIT NONE
    DOUBLE PRECISION, INTENT(IN) :: RE ! Reynolds number

    c_dsxre_dv = (0.63d0*sqrt(re) + 4.8d0)**2
    RETURN
  END FUNCTION c_dsxre_dv

! Schiller and Naumann (1933)
!----------------------------------------------------------------->>>
  DOUBLE PRECISION FUNCTION c_ds_sn(RE)
    USE param
    USE param1
    IMPLICIT NONE
    DOUBLE PRECISION, INTENT(IN) :: RE ! Reynolds number

    c_ds_sn = 24.d0*(1.d0 + 0.15d0*re**0.687d0)/(re+small_number)
    RETURN
  END FUNCTION c_ds_sn
!-----------------------------------------------------------------<<<

END MODULE drag