Changeset 38 for mso/eml/stage_separators/reboiler.mso
 Timestamp:
 Oct 23, 2006, 5:26:39 PM (17 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

mso/eml/stage_separators/reboiler.mso
r1 r38 141 141 end 142 142 143 #* 144 * Model of a dynamic reboiler with reaction 145 **# 146 Model reboilerReact 147 PARAMETERS 148 ext PP as CalcObject; 149 ext NComp as Integer; 150 Across as area (Brief="Cross Section Area of reboiler"); 151 V as volume (Brief="Total volume of reboiler"); 152 153 stoic(NComp) as Real(Brief="Stoichiometric matrix"); 154 Hr as energy_mol; 155 Pstartup as pressure; 156 157 VARIABLES 158 in Inlet as stream; #(Brief="Feed Stream"); 159 in InletL as stream; #(Brief="Liquid inlet stream"); 160 out OutletL as stream_therm; #(Brief="Liquid outlet stream"); 161 out OutletV as stream_therm; #(Brief="Vapour outlet stream"); 162 163 Q as heat_rate (Brief="Heat supplied"); 164 M(NComp) as mol (Brief="Molar Holdup in the tray"); 165 ML as mol (Brief="Molar liquid holdup"); 166 MV as mol (Brief="Molar vapour holdup"); 167 E as energy (Brief="Total Energy Holdup on tray"); 168 vL as volume_mol (Brief="Liquid Molar Volume"); 169 vV as volume_mol (Brief="Vapour Molar volume"); 170 Level as length (Brief="Level of liquid phase"); 171 Vol as volume; 172 startup as Real; 173 rhoV as dens_mass; 174 r as reaction_mol (Brief = "Reaction resulting ethyl acetate", Unit = "mol/l/s"); 175 C(NComp) as conc_mol (Brief = "Molar concentration", Lower = 1); 176 177 EQUATIONS 178 "Molar Concentration" 179 OutletL.z = vL * C; 180 181 "Component Molar Balance" 182 diff(M)= Inlet.F*Inlet.z + InletL.F*InletL.z 183  OutletL.F*OutletL.z  OutletV.F*OutletV.z + stoic*r*ML*vL; 184 185 "Energy Balance" 186 diff(E) = Inlet.F*Inlet.h + InletL.F*InletL.h 187  OutletL.F*OutletL.h  OutletV.F*OutletV.h + Q + Hr * r * vL*ML; 188 189 "Molar Holdup" 190 M = ML*OutletL.z + MV*OutletV.z; 191 192 "Energy Holdup" 193 E = ML*OutletL.h + MV*OutletV.h  OutletL.P*V; 194 195 "Mol fraction normalisation" 196 sum(OutletL.z)=1.0; 197 198 "Liquid Volume" 199 vL = PP.LiquidVolume(OutletL.T, OutletL.P, OutletL.z); 200 "Vapour Volume" 201 vV = PP.VapourVolume(OutletV.T, OutletV.P, OutletV.z); 202 "Vapour Density" 203 rhoV = PP.VapourDensity(OutletV.T, OutletV.P, OutletV.z); 204 205 "Level of liquid phase" 206 Level = ML*vL/Across; 207 208 Vol = ML*vL; 209 210 "vaporization fraction " 211 OutletV.v = 1.0; 212 OutletL.v = 0.0; 213 214 "Mechanical Equilibrium" 215 OutletL.P = OutletV.P; 216 217 "Thermal Equilibrium" 218 OutletL.T = OutletV.T; 219 220 "Geometry Constraint" 221 V = ML*vL + MV*vV; 222 223 "Chemical Equilibrium" 224 PP.LiquidFugacityCoefficient(OutletL.T, OutletL.P, OutletL.z)*OutletL.z = 225 PP.VapourFugacityCoefficient(OutletV.T, OutletV.P, OutletV.z)*OutletV.z; 226 227 sum(OutletL.z)=sum(OutletV.z); 228 229 end
Note: See TracChangeset
for help on using the changeset viewer.