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

1     MODULE energy
2     
3       USE fldvar
4       USE param
5       USE param1
6     
7     !   Gas-phase heat of reaction
8       DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  HOR_g
9     
10     !   Solids-phase heat of reaction
11       DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  HOR_s
12     
13     !   Gas-solids heat transfer coefficient
14       DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  GAMA_gs
15     
16     !   Gas-phase radiation coefficient
17       DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  GAMA_Rg
18     
19     !   Solids-phase radiation coefficient
20       DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  GAMA_Rs
21     
22     !   Gas-phase radiation temperature
23       DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  T_Rg
24     
25     !   Solids-phase radiation temperature
26       DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  T_Rs
27     
28     CONTAINS
29     
30       !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
31       !  Source terms for the energy equations.  The source term is linearized as
32       !  S = S_c - S_p * T,  where S_c and S_p must be positive.
33       !
34       !  By default the source terms have been coded for radiation sources.
35     
36       !
37       !  S_p for gas phase at i, j, k
38       DOUBLE PRECISION FUNCTION S_Rpg(IJK)
39         IMPLICIT NONE
40         INTEGER IJK
41         S_Rpg = 4.d0 * GAMA_Rg(IJK) *  T_g(IJK)**3
42       END FUNCTION S_Rpg
43     
44       !  S_c for gas phase at i, j, k
45       DOUBLE PRECISION FUNCTION S_Rcg(IJK)
46         IMPLICIT NONE
47         INTEGER IJK
48         S_Rcg = GAMA_Rg(IJK) * ( T_Rg(IJK)**4 + 3.d0 * T_g(IJK)**4 )
49       END FUNCTION S_Rcg
50     
51       !  S_p for solids phase at i, j, k
52       DOUBLE PRECISION FUNCTION S_Rps(IJK, M)
53         IMPLICIT NONE
54         INTEGER IJK, M
55         S_Rps = 4.d0 * GAMA_Rs(IJK, M) *  T_s(IJK, M)**3
56       END FUNCTION S_Rps
57     
58       !  S_c for solids phase at i, j, k
59       DOUBLE PRECISION FUNCTION S_Rcs(IJK, M)
60         IMPLICIT NONE
61         INTEGER IJK, M
62         S_Rcs = GAMA_Rs(IJK, M) * ( T_Rs(IJK, M)**4 &
63              + 3.d0 * T_s(IJK, M)**4 )
64       END FUNCTION S_Rcs
65     
66     END MODULE energy
67