1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! SUBROUTINE: TRANSPORT_PROP C 4 ! Purpose: Calculate the indicated transport properties that vary C 5 ! with time if directed to do so by the corresponding flag C 6 ! C 7 ! Author: M. Syamlal Date: 17-JUL-92 C 8 ! Reviewer: P. Nicoletti Date: 11-DEC-92 C 9 ! C 10 ! Revision Number: 1 C 11 ! Purpose: Mods for MFIX 2.0 (old name CALC_PHYSPROP) C 12 ! Author: M. Syamlal Date: 23-APR-96 C 13 ! Reviewer: Date: dd-mmm-yy C 14 ! C 15 ! Literature/Document References: C 16 ! C 17 ! Variables referenced: None C 18 ! Variables modified: None C 19 ! Local variables: None C 20 ! C 21 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 22 SUBROUTINE TRANSPORT_PROP(IER) 23 24 ! Global Variables: 25 !---------------------------------------------------------------------- 26 ! Number of solids phases. 27 use physprop, only: MMAX 28 ! Flags for calculating viscosity. 29 use coeff, only: VISC 30 ! Flags for calculating conductivity. 31 use coeff, only: COND 32 ! Flags for calculating diffusivity. 33 use coeff, only: DIFF 34 ! Flags for calculating particle-particle energy dissipation. 35 use coeff, only: GRAN_DISS 36 ! Kinetic theory model. 37 use run, only: KT_TYPE_enum 38 use run, only: ia_2005, gd_1999, gtsh_2012 39 40 implicit none 41 42 ! Dummy arguments 43 !----------------------------------------------------------------------- 44 ! Error index 45 INTEGER, intent(inout) :: IER 46 47 ! Local variables 48 !----------------------------------------------------------------------- 49 ! Loop counter 50 INTEGER :: M ! Solids phase 51 52 53 IF (VISC(0)) CALL CALC_MU_G (IER) ! Fluid viscosity 54 IF (COND(0)) CALL CALC_K_G() ! Fluid conductivity 55 IF (DIFF(0)) CALL CALC_DIF_G() ! Fluid diffusivity 56 57 DO M = 1, MMAX 58 59 ! Particle-Particle Energy Dissipation 60 ! for gtsh theory this call needs to be done before calc_mu_s so that 61 ! the cooling rate is available for mu_s 62 IF (GRAN_DISS(M)) THEN 63 SELECT CASE (KT_TYPE_ENUM) 64 CASE (IA_2005) 65 CALL CALC_IA_ENERGY_DISSIPATION_SS(M) 66 CASE(GD_1999) 67 CALL CALC_GD_99_ENERGY_DISSIPATION_SS(M, IER) 68 CASE(GTSH_2012) 69 CALL CALC_GTSH_ENERGY_DISSIPATION_SS(M) 70 END SELECT 71 ENDIF 72 ! these were moved after gran_diss since some quantities above are 73 ! needed in the subsequent gtsh calculations 74 IF (COND(M)) CALL CALC_K_S (M) ! Solids conductivity 75 IF (VISC(M)) CALL CALC_MU_S (M, IER) ! Solids viscosity 76 IF (DIFF(M)) CALL CALC_DIF_S (M) ! Solids diffusivity 77 ENDDO 78 79 RETURN 80 END SUBROUTINE TRANSPORT_PROP 81