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

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
2     !                                                                      !
3     !  Module name: DES_RADIATION                                          !
4     !                                                                      !
5     !  Purpose:                                                            !
6     !                                                                      !
7     !  Author: J.Musser                                   Date: 25-Jun-10  !
8     !                                                                      !
9     !  Commen:                                                             !
10     !                                                                      !
11     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
12           SUBROUTINE DES_RADIATION(I, iM, iIJK, FOCUS)
13     
14           USE constant
15           USE des_thermo
16           USE discretelement
17           USE fldvar
18           USE param1
19           USE physprop, only: SMAX
20           USE toleranc
21     
22           IMPLICIT NONE
23     
24     ! Passed variables
25     !---------------------------------------------------------------------//
26     ! Global index of particle
27           INTEGER, INTENT(IN) :: I
28     ! Solid phase of particle I
29           INTEGER, INTENT(IN) :: iM
30     ! Fluid cell index containing particle I
31           INTEGER, INTENT(IN) :: iIJK
32     ! Logical used for debugging
33           LOGICAL, INTENT(IN) :: FOCUS
34     
35     ! Local variables
36     !---------------------------------------------------------------------//
37     ! Surface area of particle
38           DOUBLE PRECISION :: A_S
39     ! Radiative heat transfer
40           DOUBLE PRECISION :: Qrd
41     ! Environment temperature
42           DOUBLE PRECISION :: Tenv
43     ! Particle Emmisivity
44           DOUBLE PRECISION :: lEm
45     
46     ! Set the environment temperature.
47           IF(COMPARE(EP_g(iIJK),ONE)) THEN
48              Tenv = T_g(iIJK)
49           ELSE
50              Tenv = EP_g(iIJK)*T_g(iIJK) + (ONE-EP_g(iIJK))*avgDES_T_s(iIJK)
51           ENDIF
52     
53     ! Set the particle emmisivity. Phase shift needed for TFM/DEM hybrid.
54           lEM= DES_Em(iM + SMAX)
55     
56     ! Calculate the surface area of the particle
57           A_S = 4.0d0 * Pi * DES_RADIUS(I) * DES_RADIUS(I)
58     ! Calculate the heat source.
59           Qrd = SB_CONST * A_s * lEm * (Tenv**4 - (DES_T_s_NEW(I))**4)
60     ! Update the thermal source term.
61           Q_Source(I) = Q_Source(I) + Qrd
62     
63           IF(FOCUS)THEN
64              WRITE(*,"(//5X,A)")'From: DES_RADIATION -'
65              WRITE(*,"(8X,A,D13.6)")'Tp: ',DES_T_s_NEW(I)
66              WRITE(*,"(8X,A,D13.6)")'Tenv: ',Tenv
67              WRITE(*,"(8X,A,D13.6)")'Qrd: ',Qrd
68              WRITE(*,"(5X,25('-')/)")
69           ENDIF
70     
71           RETURN
72     
73           END SUBROUTINE DES_RADIATION
74