File: RELATIVE:/../../../mfix.git/model/calc_k_g.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
2     !                                                                      C
3     !  Module name: CALC_K_g(IER)                                          C
4     !  Purpose: Calculate the effective conductivity of fluid phase        C
5     !                                                                      C
6     !  Author:M. Syamlal                                  Date: 24-APR-96  C
7     !  Reviewer:                                          Date: dd-mmm-yy  C
8     !                                                                      C
9     !  Revision Number: 1                                                  C
10     !  Purpose: allow SI unit                                              C
11     !  Author: S. Dartevelle                              Date: 01-Jul-02  C
12     !  Reviewer:                                          Date: dd-mmm-yy  C
13     !                                                                      C
14     !  Literature/Document References:                                     C
15     !                                                                      C
16     !  Variables referenced:                                               C
17     !  Variables modified:                                                 C
18     !                                                                      C
19     !  Local variables:                                                    C
20     !                                                                      C
21     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
22     
23           SUBROUTINE CALC_K_G()
24     !...Translated by Pacific-Sierra Research VAST-90 2.06G5  12:17:31  12/09/98
25     !...Switches: -xf
26     
27     !-----------------------------------------------
28     !   M o d u l e s
29     !-----------------------------------------------
30           USE param
31           USE param1
32           USE parallel
33           USE physprop
34           USE fldvar
35           USE geometry
36           USE indices
37           USE constant
38           USE compar
39           USE run
40           USE sendrecv
41           USE functions
42           IMPLICIT NONE
43     !-----------------------------------------------
44     !   G l o b a l   P a r a m e t e r s
45     !-----------------------------------------------
46     !-----------------------------------------------
47     !   D u m m y   A r g u m e n t s
48     !-----------------------------------------------
49     !                      Indices
50           INTEGER          IJK
51     !-----------------------------------------------
52     
53           IF (K_G0 /= UNDEFINED) RETURN
54     
55     !!!!$omp parallel do private(ijk) &
56     !!!!$omp& schedule(dynamic,chunk_size)
57     
58           DO IJK = ijkstart3, ijkend3
59              IF (FLUID_AT(IJK)) THEN
60     ! Gas conductivity (air)
61     ! Bird, Stewart, and Lightfoot (1960) -- Temperature dependence from formula
62     ! 8.3-12 on p. 255 and conductivity value at 300 K from p. 263
63                 K_G(IJK) = 6.02D-5*SQRT(T_G(IJK)/300.D0)           ! cal/(s.cm.K)
64              ELSE
65                 K_G(IJK) = ZERO
66              ENDIF
67     ! 1 cal = 4.183925D0 J
68              IF (UNITS == 'SI') K_G(IJK) = 418.3925D0*K_G(IJK)      !J/s.m.K
69     
70           ENDDO
71     
72           CALL send_recv(K_G, 2)
73     
74           RETURN
75           END SUBROUTINE CALC_K_G
76     
77