# Ticket #78: Reator.mso

File Reator.mso, 3.3 KB (added by , 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 | *# |

16 | using "correntes"; |

17 | |

18 | |

19 | Model Reator |

20 | |

21 | |

22 | PARAMETERS |

23 | outer NComp as Integer (Brief = "Nº de componentes", Lower = 1); |

24 | outer 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 | |

136 | end |