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() 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 use kintheory, only: CALC_IA_ENERGY_DISSIPATION_SS, CALC_GD_99_ENERGY_DISSIPATION_SS, CALC_GTSH_ENERGY_DISSIPATION_SS 40 41 implicit none 42 43 ! Local variables 44 !----------------------------------------------------------------------- 45 ! Loop counter 46 INTEGER :: M ! Solids phase 47 48 49 IF (VISC(0)) CALL CALC_MU_G() ! Fluid viscosity 50 IF (COND(0)) CALL CALC_K_G() ! Fluid conductivity 51 IF (DIFF(0)) CALL CALC_DIF_G() ! Fluid diffusivity 52 53 DO M = 1, MMAX 54 55 ! Particle-Particle Energy Dissipation 56 ! for gtsh theory this call needs to be done before calc_mu_s so that 57 ! the cooling rate is available for mu_s 58 IF (GRAN_DISS(M)) THEN 59 SELECT CASE (KT_TYPE_ENUM) 60 CASE (IA_2005) 61 CALL CALC_IA_ENERGY_DISSIPATION_SS(M) 62 CASE(GD_1999) 63 CALL CALC_GD_99_ENERGY_DISSIPATION_SS(M) 64 CASE(GTSH_2012) 65 CALL CALC_GTSH_ENERGY_DISSIPATION_SS(M) 66 END SELECT 67 ENDIF 68 ! these were moved after gran_diss since some quantities above are 69 ! needed in the subsequent gtsh calculations 70 IF (COND(M)) CALL CALC_K_S (M) ! Solids conductivity 71 IF (VISC(M)) CALL CALC_MU_S (M) ! Solids viscosity 72 IF (DIFF(M)) CALL CALC_DIF_S (M) ! Solids diffusivity 73 ENDDO 74 75 RETURN 76 END SUBROUTINE TRANSPORT_PROP 77