File: N:\mfix\model\des\calc_thermo_des.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
2     !                                                                      !
3     !  Module name: CALC_THERMO_DES                                        !
4     !                                                                      !
5     !  Purpose: This subroutine is called from DES routines. It calls      !
6     !  functions that calculate heat and mass transfer.                    !
7     !                                                                      !
8     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
9           SUBROUTINE CALC_THERMO_DES
10     
11           USE compar
12           USE des_rxns
13           USE des_thermo
14           USE discretelement
15           USE fldvar
16           USE functions
17           USE geometry
18           USE indices
19           USE interpolation
20           USE param1
21           USE run
22     
23           IMPLICIT NONE
24     
25     ! Local variables
26     !---------------------------------------------------------------------//
27           INTEGER :: NP
28     ! Functions
29     !---------------------------------------------------------------------//
30     
31     ! This is a quick work-around to keep the thermo routines from causes
32     ! issues while the "check_data" routines are rewritten. Moving forward
33     ! this routine should be split apart to avoid the particle loops for
34     ! cold-flow, non-reacting cases.
35           IF(.NOT.ENERGY_EQ .AND. .NOT.ANY_SPECIES_EQ) RETURN
36     
37     ! Calculate time dependent physical properties
38           FORALL(NP=1:MAX_PIP, PARTICLE_STATE(NP)==NORMAL_PARTICLE) &
39              DES_C_PS(NP) = CALC_CP_DES(NP)
40     
41           IF(DES_EXPLICITLY_COUPLED) THEN
42     ! Apply the convective heat transfer calculated by the gas phase.
43              IF(CALC_CONV_DES) THEN
44                 WHERE(PARTICLE_STATE == NORMAL_PARTICLE) &
45                    Q_Source = Q_Source + CONV_Qs
46              ENDIF
47              IF(ANY_SPECIES_EQ) THEN
48                 WHERE(PARTICLE_STATE == NORMAL_PARTICLE) &
49                    Q_Source = Q_Source + RXNS_Qs
50              ENDIF
51           ELSE
52              IF(CALC_CONV_DES)CALL CONV_GS_DES1
53              IF(any(CALC_RADT_DES)) CALL DES_RADIATION
54              IF(ANY_SPECIES_EQ) CALL RXNS_GS_DES1
55           ENDIF
56     
57           END SUBROUTINE CALC_THERMO_DES
58