1 !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC 2 ! C 3 ! Module name: physprop C 4 ! Purpose: Common block containing physical property data C 5 ! C 6 ! Author: M. Syamlal Date: dd-mmm-yy C 7 ! Reviewer: Date: dd-mmm-yy C 8 ! C 9 ! Literature/Document References: C 10 ! C 11 ! C 12 !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C 13 MODULE physprop 14 15 !----------------------------------------------- 16 ! Modules 17 !----------------------------------------------- 18 Use param, only: dim_m, dim_n, dim_n_g, dim_n_s 19 !----------------------------------------------- 20 21 22 ! Scale factor for gas turbulence length scale 23 DOUBLE PRECISION :: K_scale 24 25 ! Number of solids phases 26 INTEGER :: MMAX 27 28 ! Real number of solids phases for GHD theory 29 INTEGER :: SMAX 30 31 ! Particle diameters 32 DOUBLE PRECISION :: D_p0(DIM_M) 33 34 ! Constant solids phase densities. 35 DOUBLE PRECISION :: RO_s0(DIM_M) 36 37 ! Baseline/Unreacted solids phase phase density. This value is only used 38 ! for variable solids density simulations in EOSS. 39 DOUBLE PRECISION :: BASE_ROs(DIM_M) 40 41 ! Constant solids phase species mass fractions. These values delinate 42 ! the baseline/initial solids phase composition for ariable density 43 DOUBLE PRECISION :: X_S0(DIM_M, DIM_N_s) 44 45 ! Density of solid species (constant) 46 DOUBLE PRECISION :: RO_Xs0(DIM_M, DIM_N_s) 47 48 ! The index of an inert solids phase species. This is needed for 49 ! calculating the variable solids phase density. 50 INTEGER :: INERT_SPECIES(DIM_M) 51 52 ! Particle shape factor 53 DOUBLE PRECISION :: SHAPE_FACTOR(DIM_M) 54 55 ! Specified constant solids viscosity 56 DOUBLE PRECISION MU_s0 57 58 ! Flag indicates whether the phase becomes close-packed at ep_star 59 LOGICAL :: CLOSE_PACKED (DIM_M) 60 61 ! Specified constant gas density 62 DOUBLE PRECISION RO_g0 63 64 ! Specified constant gas viscosity 65 DOUBLE PRECISION MU_g0 66 67 ! Virtual (added) mass coefficient Cv 68 DOUBLE PRECISION Cv 69 70 ! Gas viscosity 71 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: MU_g 72 73 ! Average molecular weight of gas 74 DOUBLE PRECISION MW_AVG 75 76 ! Constant constant-pressure specific heat of gas 77 DOUBLE PRECISION C_pg0 78 79 ! Reference temperature for enthalpy calculations (K) 80 DOUBLE PRECISION, PARAMETER :: T_ref = 298 81 82 ! Constant pressure specific heat of gas 83 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: C_pg 84 85 ! Constant constant-pressure specific heat of solids 86 DOUBLE PRECISION C_ps0(DIM_M) 87 88 ! Constant pressure specific heat of solids 89 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: C_ps 90 91 ! Specified constant gas conductivity 92 DOUBLE PRECISION :: K_g0 93 94 ! Conductivity of gas 95 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: K_g 96 97 ! Specified constant solids conductivity 98 DOUBLE PRECISION K_s0(DIM_M) 99 100 ! Conductivity of solids 101 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: K_s 102 103 ! Granular Temperature Conductivity (associated with temperature grad) 104 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: Kth_s 105 106 ! Granular Temperature Conductivity (associated with volume fraction grad) 107 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: Kphi_s 108 109 ! Specified constant gas diffusivity 110 DOUBLE PRECISION DIF_g0 111 112 ! Diffusivity of gas species N 113 DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE :: DIF_g 114 115 ! Specified constant solids diffusivity 116 DOUBLE PRECISION DIF_s0 117 118 ! Diffusivity of solids species N 119 DOUBLE PRECISION, DIMENSION(:, :, :), ALLOCATABLE :: DIF_s 120 121 ! Total number of gas or solids species 122 INTEGER :: NMAX(0:DIM_M) ! Runtime (all phases) 123 INTEGER :: NMAX_g ! Number of gas phase species 124 INTEGER :: NMAX_s(DIM_M) ! Number of solids phase species 125 126 ! Molecular weight of gas species 127 DOUBLE PRECISION :: MW_g (DIM_N_g) 128 129 ! Molecular weight of solids species 130 DOUBLE PRECISION :: MW_s (DIM_M, DIM_N_s) 131 132 ! Molecular weight of gas mixture 133 DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE :: MW_MIX_g 134 135 ! Logical for reading thermochemical database. 136 LOGICAL :: DATABASE_READ = .TRUE. 137 138 ! Polynomical coefficients for calculating specific heat. 139 DOUBLE PRECISION Alow (7,0:DIM_M, DIM_N) ! Tlow --> Tcom 140 DOUBLE PRECISION Ahigh(7,0:DIM_M, DIM_N) ! Tcom --> Thigh 141 142 ! Range where the polynomials are valid. 143 DOUBLE PRECISION Thigh(0:DIM_M, DIM_N) ! Upper bound 144 DOUBLE PRECISION Tlow (0:DIM_M, DIM_N) ! Lower bound 145 DOUBLE PRECISION Tcom (0:DIM_M, DIM_N) ! Switch from low to high 146 147 ! Heat of formation at Tref divided by the gas constant. 148 DOUBLE PRECISION HfrefoR(0:DIM_M, DIM_N) 149 150 ! Reference values. 151 DOUBLE PRECISION ICpoR_l(0:DIM_M, DIM_N) 152 DOUBLE PRECISION ICpoR_h(0:DIM_M, DIM_N) 153 154 END MODULE physprop 155