1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv! 2 ! ! 3 ! Module name: DES_THERMO ! 4 ! ! 5 ! Purpose: Common elements for MFIX-DEM heat transfer. ! 6 ! condition. ! 7 ! ! 8 ! Author: J.Musser Date: 16-Jun-10 ! 9 ! ! 10 ! Comments: ! 11 ! ! 12 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^! 13 MODULE DES_THERMO 14 15 USE param 16 17 ! Heat transfer correlation specified in mfix.dat 18 ! Default [RANZ_1952] 19 CHARACTER(LEN=24) :: DES_CONV_CORR 20 21 INTEGER :: DES_CONV_CORR_ENUM 22 INTEGER, PARAMETER :: RANZ_1952 = 0 23 24 ! Particle properties 25 !-------------------------- 26 ! Particle temperature at current time step (S_TIME) 27 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: DES_T_s_OLD !(PARTICLES) 28 ! Particle temperature at previous time step (S_TIME - DT_SOLID) 29 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: DES_T_s_NEW !(PARTICLES) 30 31 ! DES specific heat of particles by particle 32 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: DES_C_ps 33 34 ! Emissivity of particles 35 DOUBLE PRECISION DES_Em(DIM_M) 36 ! Stefan-Boltzmann Constant 37 DOUBLE PRECISION SB_CONST 38 ! Bulk solids temperature for radiative heat transfer 39 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: avgDES_T_s 40 41 42 ! Rates of heat transfer 43 !-------------------------- 44 ! Internal heat generation resulting from a chemical reaction 45 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: Qint !(PARTICLES) 46 ! Heat transfer source to the particle. 47 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: Q_Source 48 ! Previous time step's rate of heat transfer. Used for Adams-Bashforth 49 ! time integration scheme 50 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: Q_Source0 51 52 ! Thermodynamic Neighborhood 53 !-------------------------- 54 ! Fluid Lens Proportion Constant used to calculate the radius of the 55 ! fluid lens that surrounds the particle for particle-fluid-particle 56 ! conduction. Default [ 0.2 ] 57 DOUBLE PRECISION FLPC 58 59 ! Mininum separation distance between the surface of two contacting 60 ! particles. This value is used to remove the singluarity that the 61 ! particle-fluid-particle conduciton model develops at the contact 62 ! interface. [4.0x10^(-10) meters] 63 DOUBLE PRECISION DES_MIN_COND_DIST 64 65 66 ! Run time flags for calculating the various modes of heat transfer 67 LOGICAL :: CALC_CONV_DES 68 LOGICAL :: CALC_COND_DES(DIM_M) 69 LOGICAL :: CALC_RADT_DES(DIM_M) 70 71 ! Fluid/Particle coupling 72 !---------------------------------------------------------------------// 73 ! Source term for TFM gas phase energy equation. 74 DOUBLE PRECISION, ALLOCATABLE :: DES_ENERGY_SOURCE(:) 75 76 77 END MODULE DES_THERMO 78