File: /nfs/home/0/users/jenkins/mfix.git/model/GhdTheory/transport_coeff_ghd.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
2     !                                                                      C
3     !  Subroutine name: TRANSPORT_COEFF_GHD                                C
4     !  Purpose: Calculate all transport coefficients for use in GHD theory C
5     !                                                                      C
6     !  Author: S. Benyahia                              Date: 03-MAR-09    C
7     !  Reviewer:                                          Date:            C
8     !                                                                      C
9     !  Literature/Document References:                                     C
10     !     C. Hrenya handnotes and Garzo, Hrenya, Dufty papers (PRE, 2007)  C
11     !                                                                      C
12     !                                                                      C
13     !  Variables modified:  All Transport Coefficients and species temp.   C
14     !                                                                      C
15     !                                                                      C
16     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
17     
18           SUBROUTINE TRANSPORT_COEFF_GHD (M)
19     
20     !-----------------------------------------------
21     ! Modules
22     !-----------------------------------------------
23           USE param
24           USE param1
25           USE geometry
26           USE compar
27           USE fldvar
28           USE indices
29           USE visc_s
30           USE ghdtheory
31           USE physprop
32           USE run
33           USE constant
34           USE toleranc
35           USE functions
36           IMPLICIT NONE
37     !-----------------------------------------------
38     ! Local variables
39     !-----------------------------------------------
40     ! Index
41           INTEGER :: IJK, I, J, K
42     ! Solids phase
43           INTEGER :: M, L
44     ! particles properties and GHD transport coefficients
45           DOUBLE PRECISION :: SIGMAI(smax), Mi(smax), phii(smax)
46           DOUBLE PRECISION :: Ti(smax), tmpDT(smax)
47           DOUBLE PRECISION :: tmpZeta0, tmpZetaU, TMix
48           DOUBLE PRECISION :: tmpP, tmpKappa, tmpEta, tmpLambda
49           DOUBLE PRECISION :: tmpLij(smax,smax), tmpDij(smax,smax), &
50                               tmpDF(smax,smax), tmpDijQ(smax,smax)
51     !-----------------------------------------------
52     ! Function subroutines
53     !-----------------------------------------------
54     
55           DO 200 IJK = ijkstart3, ijkend3
56              I = I_OF(IJK)
57              J = J_OF(IJK)
58              K = K_OF(IJK)
59     
60              IF ( FLUID_AT(IJK) ) THEN
61     
62                 DO M = 1, SMAX
63                    SIGMAI(M) = D_P(IJK,M)
64                    Mi(M)     = (PI/6.d0)*SIGMAI(M)**3 * RO_S(IJK,M)
65                    phii(M)   = ROP_S(IJK,M) / RO_S(IJK,M)
66                    Ti(M)= THETA_M(IJK,M)
67                 ENDDO
68                 TMix = THETA_M(IJK,MMAX)
69     
70                 CALL GHD_MODEL(SMAX, SIGMAI, IJK, r_p(:smax,:smax), Mi, &
71                          phii, TMix, tmpZeta0, tmpZetaU, Ti, tmpP, &
72                          tmpKappa, tmpEta, tmpDT, tmpDF, tmpLambda, &
73                          tmpLij, tmpDij, tmpDijQ)
74     
75     ! zeroth order cooling rate
76                    Zeta0(IJK) =  tmpZeta0
77     
78     ! first order cooling rate transport coefficient
79                    ZetaU(IJK) =  tmpZetaU
80     
81     ! species temperature Ti and GHD-specific transport properties
82                    DO M = 1, SMAX
83                      THETA_M(IJK,M) = Ti(M)
84     ! Thermal diffusivity
85                      DiT(IJK,M) = tmpDT(M)
86                      DO L = 1, SMAX
87     ! mass mobility coefficient
88                        DijF(IJK,M,L) = tmpDF(M,L)
89     ! thermal mobility coefficient
90                        Lij(IJK,M,L) = tmpLij(M,L)
91     ! ordinary diffucsion coefficient
92                        Dij(IJK,M,L) = tmpDij(M,L)
93     ! Dufour coefficient:
94                        DijQ(IJK,M,L) = tmpDijQ(M,L)
95                      ENDDO
96                    ENDDO
97     
98     ! mixture solids pressure
99                    P_s_v(IJK) = tmpP
100     
101     ! mixture bulk viscosity
102                    Mu_b_v(IJK) = tmpKappa
103     
104     ! mixture kinetic-collisional viscosity
105                    Mu_s_v(IJK) = tmpEta
106     
107     ! mixture total bulk viscosity
108                    LAMBDA_S_V(IJK) = Mu_b_v(IJK) - (2.d0/3.d0)*Mu_s_v(IJK)
109     
110     ! granular conductivity in Mth solids phase
111                    Kth_s(IJK,M) = tmpLambda
112     
113               ENDIF     ! Fluid_at
114      200  CONTINUE     ! outer IJK loop
115     
116           RETURN
117           END SUBROUTINE TRANSPORT_COEFF_GHD
118