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