Changeset 270 for trunk/sample/stage_separators/sample_column.mso
 Timestamp:
 Jun 16, 2007, 4:18:47 PM (16 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

trunk/sample/stage_separators/sample_column.mso
r247 r270 281 281 end 282 282 283 FlowSheet Distillation_thermosyphon_subcooling_Test 284 PARAMETERS 285 PP as Plugin(Brief="Physical Properties", 286 Type="PP", 287 Components = [ "isobutane", "npentane", "propylene", 288 "benzene", "isobutene" ], 289 LiquidModel = "PR", 290 VapourModel = "PR" 291 ); 292 NComp as Integer; 293 294 VARIABLES 295 Qc as heat_rate (Brief="Heat rate removed from condenser"); 296 Qr as heat_rate (Brief="Heat rate supplied to reboiler"); 297 Qttop as heat_rate (Brief="Heat rate removed from condenser"); 298 Qtbottom as heat_rate (Brief="Heat rate supplied to reboiler"); 299 300 SET 301 NComp = PP.NumberOfComponents; 302 303 DEVICES 304 col as Distillation_thermosyphon_subcooling; 305 feed as source; 306 zero as stream; 307 308 CONNECTIONS 309 feed.Outlet to col.trays(5).Inlet; 310 zero to col.trays([1:4]).Inlet; 311 zero to col.trays([6:col.NTrays]).Inlet; 312 Qttop to col.ttop.Q; 313 Qtbottom to col.tbottom.Q; 314 Qc to col.cond.Q; 315 Qr to col.reb.Q; 316 317 SPECIFY 318 feed.Outlet.F = 113.4 * 'kmol/h'; 319 feed.Outlet.T = 291 * 'K'; 320 feed.Outlet.P = 168.3 * 'kPa'; 321 feed.Outlet.z = 1/NComp; 322 323 zero.F = 0 * 'kmol/h'; 324 zero.T = 300 * 'K'; 325 zero.P = 1 * 'atm'; 326 zero.z = 1/NComp; 327 zero.v = 0; 328 zero.h = 0 * 'J/mol'; 329 330 col.sptop.Outlet2.F = 85 * 'kmol/h'; 331 col.sptop.frac = 0.444445; 332 col.spbottom.Outlet1.F = 100 * 'kmol/h'; 333 334 Qr = 3.7743e6 * 'kJ/h'; 335 Qc = 3.71e6 * 'kJ/h'; 336 Qttop = 0 * 'kJ/h'; 337 Qtbottom = 0 * 'kJ/h'; 338 339 col.pump1.dP = 16 * 'kPa'; 340 col.trays.Emv = 1; 341 col.cond.DP = 0.5 * 'atm'; 342 343 SET 344 col.NTrays = 8; 345 col.trays.V = 4 * 'ft^3'; 346 col.trays.Ah = 0.394 * 'ft^2'; 347 col.trays.lw = 20.94 * 'in'; 348 col.trays.hw = 0.125 * 'ft'; 349 col.trays.Q = 0 * 'kW'; 350 col.trays.beta = 0.6; 351 col.trays.alfa = 4; 352 col.trays.Ap = 3.94 * 'ft^2'; 353 col.reb.DP = 0.5 * 'atm'; 354 col.tbottom.Across = 2.20 * 'm^2'; 355 356 INITIAL 357 # ttop 358 col.ttop.Outlet.T = 260 *'K'; 359 col.ttop.Level = 1 * 'm'; 360 col.ttop.Outlet.z([1:4]) = [0.65, 0.05, 0.01, 0.01]; 361 362 # tbottom 363 col.tbottom.Outlet.T = 330 *'K'; 364 col.tbottom.Level = 1 * 'm'; 365 col.tbottom.Outlet.z([1:4]) = [0.1, 0.7, 0.01, 0.01]; 366 367 # column trays 368 col.trays.OutletL.T = [290:(330290)/(col.NTrays1):330] * 'K'; 369 col.trays.Level = 1.2 * col.trays.hw; 370 col.trays.OutletL.z([1:4]) = [0.5, 0.05, 0.01, 0.01]; 371 372 OPTIONS 373 TimeStep = 0.1; 374 TimeEnd = 50; 375 #time = [0:0.01:1, 2:50]; 376 #GuessFile="Distillation_kettle_cond_Test.rlt"; 377 #Dynamic = false; 378 end 379 283 380 FlowSheet Column_ctrl 284 381 PARAMETERS … … 493 590 end 494 591 592
Note: See TracChangeset
for help on using the changeset viewer.