File: N:\mfix\model\ic_mod.f

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvv!
2     !                                                                      !
3     !  Module: ic                                                          !
4     !  Author: M. Syamlal                                 Date: dd-mmm-yy  !
5     !                                                                      !
6     !  Purpose: Global initial conditions variables.                       !
7     !                                                                      !
8     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^!
9           MODULE ic
10     
11     ! Maximum number of IC regions.
12           use param, only: DIMENSION_IC
13     ! Maximum number of solids phases.
14           use param, only: DIM_M
15     ! Maximum number of gas phase species
16           use param, only: DIM_N_g
17     ! Maximum number of solids phase species
18           use param, only: DIM_N_s
19     ! Maximum number of scalar equations
20           use param, only: DIM_Scalar
21     
22     ! IC region West face, X-coordinate
23           DOUBLE PRECISION :: IC_X_w (DIMENSION_IC)
24     
25     ! IC region East face, X-coordinate
26           DOUBLE PRECISION :: IC_X_e (DIMENSION_IC)
27     
28     ! IC region South face, Y-coordinate
29           DOUBLE PRECISION :: IC_Y_s (DIMENSION_IC)
30     
31     ! IC region North face, Y-coordinate
32           DOUBLE PRECISION :: IC_Y_n (DIMENSION_IC)
33     
34     ! IC region Bottom face, Z-coordinate
35           DOUBLE PRECISION :: IC_Z_b (DIMENSION_IC)
36     
37     ! IC region Top face, Z-coordinate
38           DOUBLE PRECISION :: IC_Z_t (DIMENSION_IC)
39     
40     ! IC region, West face, I Index
41           INTEGER :: IC_I_w (DIMENSION_IC)
42     
43     ! IC region, East face, I Index
44           INTEGER :: IC_I_e (DIMENSION_IC)
45     
46     ! IC region, South face, J Index
47           INTEGER :: IC_J_s (DIMENSION_IC)
48     
49     ! IC region, North face, J Index
50           INTEGER :: IC_J_n (DIMENSION_IC)
51     
52     ! IC region, Bottom face, K Index
53           INTEGER :: IC_K_b (DIMENSION_IC)
54     
55     ! IC region, Top face, K Index
56           INTEGER :: IC_K_t (DIMENSION_IC)
57     
58     ! Type of initial condition: PATCH
59           CHARACTER(LEN=16) :: IC_TYPE(DIMENSION_IC)
60     
61     ! Initial gas phase volume fraction
62           DOUBLE PRECISION :: IC_EP_g (DIMENSION_IC)
63     
64     ! Initial gas pressure
65           DOUBLE PRECISION :: IC_P_g (DIMENSION_IC)
66     
67     ! Initial gas pressure
68           DOUBLE PRECISION :: IC_P_star(DIMENSION_IC)
69     
70     ! Initial turbulence length scale
71           DOUBLE PRECISION :: IC_L_scale(DIMENSION_IC)
72     
73     ! Initial macroscopic density of solids phases
74           DOUBLE PRECISION :: IC_ROP_s(DIMENSION_IC, DIM_M)
75     
76     ! Initial solids phase volume fraction
77           DOUBLE PRECISION :: IC_EP_s (DIMENSION_IC, DIM_M)
78     
79     ! Initial gas phase temperature
80           DOUBLE PRECISION :: IC_T_g(DIMENSION_IC)
81     
82     ! Initial solids phase temperature
83           DOUBLE PRECISION :: IC_T_s(DIMENSION_IC, DIM_M)
84     
85     ! Initial granular temperature
86           DOUBLE PRECISION :: IC_Theta_m(DIMENSION_IC, DIM_M)
87     
88     ! Initial x-component of gas velocity
89           DOUBLE PRECISION :: IC_U_g(DIMENSION_IC)
90     
91     ! Initial x-component of solids phase velocity
92           DOUBLE PRECISION :: IC_U_s(DIMENSION_IC, DIM_M)
93     
94     ! Initial y-component of gas velocity
95           DOUBLE PRECISION :: IC_V_g(DIMENSION_IC)
96     
97     ! Initial y-component of solids phase velocity
98           DOUBLE PRECISION :: IC_V_s(DIMENSION_IC, DIM_M)
99     
100     ! Initial z-component of gas velocity
101           DOUBLE PRECISION :: IC_W_g(DIMENSION_IC)
102     
103     ! Initial z-component of solids phase velocity
104           DOUBLE PRECISION :: IC_W_s(DIMENSION_IC, DIM_M)
105     
106     ! Logical variable to determine whether an ic is defined
107           LOGICAL :: IC_DEFINED (DIMENSION_IC)
108     
109     ! Initial gas species mass fractions
110           DOUBLE PRECISION :: IC_X_g(DIMENSION_IC, DIM_N_g)
111     
112     ! Initial solids species mass fractions
113           DOUBLE PRECISION :: IC_X_s(DIMENSION_IC, DIM_M, DIM_N_s)
114     
115     ! Gas phase radiation coefficient
116           DOUBLE PRECISION :: IC_GAMA_Rg (DIMENSION_IC)
117     
118     ! Gas phase radiation temperature
119           DOUBLE PRECISION :: IC_T_Rg(DIMENSION_IC)
120     
121     ! Solids phase-1 radiation coefficient
122           DOUBLE PRECISION :: IC_GAMA_Rs(DIMENSION_IC, DIM_M)
123     
124     ! Solids phase-1 radiation temperature
125           DOUBLE PRECISION :: IC_T_Rs(DIMENSION_IC, DIM_M)
126     
127     ! Initial scalar value in a region
128           DOUBLE PRECISION :: IC_Scalar(DIMENSION_IC, DIM_scalar)
129     
130     ! Initial K & Epsilon values in a region
131           DOUBLE PRECISION :: IC_K_Turb_G(DIMENSION_IC)
132           DOUBLE PRECISION :: IC_E_Turb_G(DIMENSION_IC)
133     
134     ! Initial conditions for DES cases (such as, DEM, MPPIC, hybrid)
135     
136     ! Flag to extend the lattice distribution in a given IC to available area
137           LOGICAL :: IC_DES_FIT_TO_REGION (DIMENSION_IC)
138     
139     ! Flag to specify the initial constant number of particles per cell
140     ! for the MPPIC method initialization.
141     ! Statistical weight of parcels will be calculated by the code
142           INTEGER :: IC_PIC_CONST_NPC(DIMENSION_IC, DIM_M)
143     
144     ! Flag to specify the initial constant statistical weight.
145     ! for the MPPIC method initialization.
146     ! Number of computational particles/parcels will be calculated by the code
147           DOUBLE PRECISION :: IC_PIC_CONST_STATWT(DIMENSION_IC, DIM_M)
148     
149     
150           END MODULE ic
151