transport_coeff_ghd.f

Go to the documentation of this file.

!vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
!                                                                      C
!  Subroutine name: TRANSPORT_COEFF_GHD                                C
!  Purpose: Calculate all transport coefficients for use in GHD theory C
!                                                                      C
!  Author: S. Benyahia                              Date: 03-MAR-09    C
!  Reviewer:                                          Date:            C
!                                                                      C
!  Literature/Document References:                                     C
!     C. Hrenya handnotes and Garzo, Hrenya, Dufty papers (PRE, 2007)  C
!                                                                      C
!                                                                      C
!  Variables modified:  All Transport Coefficients and species temp.   C
!                                                                      C
!                                                                      C
!^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C

      SUBROUTINE transport_coeff_ghd (M)

!-----------------------------------------------
! Modules
!-----------------------------------------------
      USE param
      USE param1
      USE geometry
      USE compar
      USE fldvar
      USE indices
      USE visc_s
      USE ghdtheory
      USE physprop
      USE run
      USE constant
      USE toleranc
      USE functions
      IMPLICIT NONE
!-----------------------------------------------
! Local variables
!-----------------------------------------------
! Index
      INTEGER :: IJK, I, J, K
! Solids phase
      INTEGER :: M, L
! particles properties and GHD transport coefficients
      DOUBLE PRECISION :: SIGMAI(smax), Mi(smax), phii(smax)
      DOUBLE PRECISION :: Ti(smax), tmpDT(smax)
      DOUBLE PRECISION :: tmpZeta0, tmpZetaU, TMix
      DOUBLE PRECISION :: tmpP, tmpKappa, tmpEta, tmpLambda
      DOUBLE PRECISION :: tmpLij(smax,smax), tmpDij(smax,smax), &
                          tmpDF(smax,smax), tmpDijQ(smax,smax)
!-----------------------------------------------
! Function subroutines
!-----------------------------------------------

      DO 200 ijk = ijkstart3, ijkend3
         i = i_of(ijk)
         j = j_of(ijk)
         k = k_of(ijk)

         IF ( fluid_at(ijk) ) THEN

            DO m = 1, smax
               sigmai(m) = d_p(ijk,m)
               mi(m)     = (pi/6.d0)*sigmai(m)**3 * ro_s(ijk,m)
               phii(m)   = rop_s(ijk,m) / ro_s(ijk,m)
               ti(m)= theta_m(ijk,m)
            ENDDO
            tmix = theta_m(ijk,mmax)

            CALL ghd_model(smax, sigmai, ijk, r_p(:smax,:smax), mi, &
                     phii, tmix, tmpzeta0, tmpzetau, ti, tmpp, &
                     tmpkappa, tmpeta, tmpdt, tmpdf, tmplambda, &
                     tmplij, tmpdij, tmpdijq)

! zeroth order cooling rate
               zeta0(ijk) =  tmpzeta0

! first order cooling rate transport coefficient
               zetau(ijk) =  tmpzetau

! species temperature Ti and GHD-specific transport properties
               DO m = 1, smax
                 theta_m(ijk,m) = ti(m)
! Thermal diffusivity
                 dit(ijk,m) = tmpdt(m)
                 DO l = 1, smax
! mass mobility coefficient
                   dijf(ijk,m,l) = tmpdf(m,l)
! thermal mobility coefficient
                   lij(ijk,m,l) = tmplij(m,l)
! ordinary diffucsion coefficient
                   dij(ijk,m,l) = tmpdij(m,l)
! Dufour coefficient:
                   dijq(ijk,m,l) = tmpdijq(m,l)
                 ENDDO
               ENDDO

! mixture solids pressure
               p_s_v(ijk) = tmpp

! mixture bulk viscosity
               mu_b_v(ijk) = tmpkappa

! mixture kinetic-collisional viscosity
               mu_s_v(ijk) = tmpeta

! mixture total bulk viscosity
               lambda_s_v(ijk) = mu_b_v(ijk) - (2.d0/3.d0)*mu_s_v(ijk)

! granular conductivity in Mth solids phase
               kth_s(ijk,m) = tmplambda

          ENDIF     ! Fluid_at
 200  CONTINUE     ! outer IJK loop

      RETURN
      END SUBROUTINE transport_coeff_ghd