Ticket #78: Reator.mso

File Reator.mso, 3.3 KB (added by Antonio José, 15 years ago)
Line 
1#*
2 * EMSO template file for new models.
3 *
4 * The user should edit the contents of this file in order
5 * to develop a new model.
6 *
7 * $Id: model.mso 114 2007-01-15 14:20:51Z rafael $
8 *#
9
10
11#*
12 * The using command makes available entities declared in another files.
13 *
14 * With this command the user can make use of the EML types
15 *#
16using "correntes";
17
18
19Model Reator
20
21       
22        PARAMETERS
23outer NComp as Integer (Brief = "Nº de componentes", Lower = 1);
24outer PP   as Plugin(Brief="Physical Properties", Type="PP");
25    An as coefficient(Unit='m^1.5*s');
26        A as area(Brief="Cross section area");
27    D as length(Brief="diametro do reator");
28        #B as Real(Unit='(atm*m^8*kg)/(K^1.5)/mol^3/s');
29       
30        VARIABLES
31        Ç(NComp) as corrente1;
32        Emf as Real(Brief="Bed voidage at minimum fluidization",Default=0.4);
33        Ar as Real;
34        rop as dens_mass(Brief="densidade do catalisador");
35        ro as dens_mass(Brief="densidade da mistura");
36        mi as viscosity(Brief="viscosidade da mistura gasosa");
37        g as acceleration(Brief="aceleração da gravidade");
38        dp as length(Brief="diametro da partícula");
39        umf as velocity(Brief="velocidade supeficial na minima fluidização",Default=0.003);
40        db as length(Brief="diametro da bolha");
41        dbm as length;
42        dbo as length;
43        uo as velocity(Brief="velocidade superficial");
44        ub as velocity(Brief="velocidade da bolha");
45        S as fraction(Brief="fração volumar da fase bolha");
46        Qe as flow_vol(Brief="vazao volumetrica da fase densa");
47        Qb as flow_vol(Brief="vazao volumetrica da fase bolha");
48    H as length(Brief="altura do leito");
49        Hmf as length(Brief="Altura do leito na minima fluidização");
50        #Kbe as coefficient(Brief="Coeficiente total de Transferencia de Massa");
51        #Kb as coefficient;
52        #Ke as coefficient;
53        Dm(NComp) as diffusivity(Brief="Difusividade do componente i na mistura gasosa");
54        Dj(NComp,NComp) as diffusivity(Brief="Difusividade do componente i no componente j");
55    x as length;
56        V(NComp) as vol_mol(Brief="Volume molar dos componentes");
57        Mw(NComp) as molweight;
58        T as temperature;
59        Te(NComp) as temperature(Brief="Temperatura normal de ebulição dos componentes",Lower=0);
60        P as pressure;
61        z(NComp) as fraction;
62        #R as Real;
63       
64        EQUATIONS
65        "Bed voidage at minimum fluidization"
66        Emf=0.586*((1/Ar)^0.029)*((ro/rop)^0.021);
67        Ar=ro*(rop-ro)*g*(dp^3)/(mi^2);
68       
69        "velocidade supeficial na minima fluidização"
70        umf=(mi/(ro*dp))*(((25.25)^2+0.0651*Ar)^0.5-25.25);
71       
72        "diametro da bolha"
73        db=dbm-(dbm-dbo)*exp(-0.3*x/D);
74        dbm=0.652*(An*abs(uo-umf))^0.4;
75        dbo=0.347*(An*abs(uo-umf))^0.4;
76               
77        "velocidade da bolha"
78        ub=uo-umf+0.711*abs(g*db)^0.5;
79       
80       
81        "fração volumar da fase bolha"
82        S=(uo-umf)/ub;
83       
84        "vazao volumetrica da fase densa"
85        Qe=umf*A;
86       
87        "vazao volumetrica da fase bolha"
88        Qb=(uo-umf)*A;
89       
90        "altura do leito"
91        H=Hmf/(1-S);
92        H/D=40;
93       
94        for i in [1:NComp]
95               
96                for j in [1:NComp]
97                       
98                        if i equal j then
99                       
100                        Dj(i,j)=0 * 'm^2/s';
101                       
102                        else
103                       
104                        Dj(i,j) = 4.357e-9*(T^1.5/(P*((V(i)*1e3)^1/3+(V(j)*1e3)^1/3)^2))*(1/Mw(i)+1/Mw(j))^1/2;
105                       
106                        end
107                       
108                end
109               
110        end
111       
112        for i in [1:NComp]
113                V(i)= PP.VapourVolume(Te(i), P, Ç(i).z);
114        end     
115       
116       
117        for i in [1:NComp]
118               
119                for j in [1:NComp]
120                       
121                        if i equal j then
122                               
123                        Dm(i)=0 * 'm^2/s';
124                       
125                        else
126                       
127                        Dm(i)=(1-z(i))/(sum(z(j)/Dj(i,j)));
128                       
129                        end
130               
131                end
132       
133        end
134       
135               
136end