transport_prop.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  SUBROUTINE: TRANSPORT_PROP                                          C
!  Purpose: Calculate the indicated transport properties that vary     C
!           with time if directed to do so by the corresponding flag   C
!                                                                      C
!  Author: M. Syamlal                                 Date: 17-JUL-92  C
!  Reviewer: P. Nicoletti                             Date: 11-DEC-92  C
!                                                                      C
!  Revision Number: 1                                                  C
!  Purpose: Mods for MFIX 2.0 (old name CALC_PHYSPROP)                 C
!  Author: M. Syamlal                                 Date: 23-APR-96  C
!  Reviewer:                                          Date: dd-mmm-yy  C
!                                                                      C
!  Literature/Document References:                                     C
!                                                                      C
!  Variables referenced: None                                          C
!  Variables modified: None                                            C
!  Local variables: None                                               C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
      SUBROUTINE transport_prop()

! Global Variables:
!----------------------------------------------------------------------
! Number of solids phases.
      use physprop, only: mmax
! Flags for calculating viscosity.
      use coeff, only: visc
! Flags for calculating conductivity.
      use coeff, only: cond
! Flags for calculating diffusivity.
      use coeff, only: diff
! Flags for calculating particle-particle energy dissipation.
      use coeff, only: gran_diss
! Kinetic theory model.
      use run, only: kt_type_enum
      use run, only: ia_2005, gd_1999, gtsh_2012
      use kintheory, only: calc_ia_energy_dissipation_ss
      use kintheory, only: calc_gd_99_energy_dissipation_ss
      use kintheory, only: calc_gtsh_energy_dissipation_ss

      implicit none

! Local variables
!-----------------------------------------------------------------------
! Loop counter
      INTEGER :: M ! Solids phase


      IF (visc(0)) CALL calc_mu_g()    ! Fluid viscosity
      IF (cond(0)) CALL calc_k_g()     ! Fluid conductivity
      IF (diff(0)) CALL calc_dif_g()   ! Fluid diffusivity

      DO m = 1, mmax
! Particle-Particle Energy Dissipation
! for gtsh theory this call needs to be done before calc_mu_s so that
! the cooling rate is available for mu_s
         IF (gran_diss(m)) THEN
         SELECT CASE (kt_type_enum)
            CASE (ia_2005)
               CALL calc_ia_energy_dissipation_ss(m)
            CASE(gd_1999)
               CALL calc_gd_99_energy_dissipation_ss(m)
            CASE(gtsh_2012)
               CALL calc_gtsh_energy_dissipation_ss(m)
            END SELECT
         ENDIF
! these were moved after gran_diss since some quantities above are
! needed in the subsequent gtsh calculations performed via calc_mu_s
         IF (cond(m)) CALL calc_k_s (m)   ! Solids conductivity
         IF (visc(m)) CALL calc_mu_s (m)  ! Solids viscosity
         IF (diff(m)) CALL calc_dif_s (m) ! Solids diffusivity
      ENDDO

      RETURN
      END SUBROUTINE transport_prop