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

1     !vvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvvC
2     !                                                                      C
3     !  Module name: mflux_mod                                              C
4     !  Purpose: Module for mass fluxes and densities at faces               C
5     !                                                                      C
6     !  Author: M. Syamlal                                 Date: dd-mmm-yy  C
7     !  Reviewer:                                          Date: dd-mmm-yy  C
8     !                                                                      C
9     !  Revision Number:                                                    C
10     !  Purpose:                                                            C
11     !  Author:                                            Date: dd-mmm-yy  C
12     !  Reviewer:                                          Date: dd-mmm-yy  C
13     !                                                                      C
14     !  Literature/Document References: None                                C
15     !                                                                      C
16     !  Variables referenced: None                                          C
17     !  Variables modified: None                                            C
18     !                                                                      C
19     !  Local variables: None                                               C
20     !                                                                      C
21     !^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^C
22     
23     
24           MODULE mflux
25     
26     
27           Use param
28           Use param1
29     
30     
31     !
32     !                      x-component of gas mass flux
33           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_gE
34     !
35     !                      x-component of solids mass flux
36           DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  Flux_sE
37     !
38     !                      y-component of gas mass flux
39           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_gN
40     !
41     !                      y-component of solids mass flux
42           DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  Flux_sN
43     !
44     !                      z-component of gas mass flux
45           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_gT
46     !
47     !                      z-component of solids mass flux
48           DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  Flux_sT
49     !
50     !                      Added mass Flux Components to be used for scalar eq.
51     !                      Note: added mass apply only to one solids phase (M=M_AM)
52           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_gSE, Flux_sSE, Flux_gSN
53           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_sSN, Flux_gST, Flux_sST
54     !
55     !
56     !                      macroscopic gas density at east face
57           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ROP_gE
58     !
59     !                      macroscopic solids density at east face
60           DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  ROP_sE
61     !
62     !                      macroscopic gas density at north face
63           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ROP_gN
64     !
65     !                      macroscopic solids density at north face
66           DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  ROP_sN
67     !
68     !                      macroscopic gas density at top face
69           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  ROP_gT
70     !
71     !                      macroscopic solids density at top face
72           DOUBLE PRECISION, DIMENSION(:, :), ALLOCATABLE ::  ROP_sT
73     
74     
75     ! for GHD Theory
76     !                      x-component of solids total number density flux
77           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_nE
78     !                      y-component of solids total number density flux
79           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_nN
80     !                      z-component of solids total number density flux
81           DOUBLE PRECISION, DIMENSION(:), ALLOCATABLE ::  Flux_nT
82     
83     ! end GHD Theory modification
84     
85     
86           END MODULE mflux
87