Changeset 477 for branches/packed/sample/stage_separators/sample_tray.mso
 Timestamp:
 Mar 7, 2008, 3:46:44 PM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/packed/sample/stage_separators/sample_tray.mso
r470 r477 321 321 TimeEnd = 60; 322 322 end 323 324 #to compare with tray_Test_1 325 FlowSheet packedStage_Test_1 326 PARAMETERS 327 PP as Plugin(Brief="Physical Properties", 328 Type="PP", 329 Components = [ "npentane", "benzene"], 330 LiquidModel = "PR", 331 VapourModel = "PR" 332 ); 333 NComp as Integer; 334 335 SET 336 NComp = PP.NumberOfComponents; 337 338 VARIABLES 339 deltaP as Real (Unit='atm/m'); 340 341 DEVICES 342 t1 as packedStage; 343 feed as source; 344 inL as liquid_stream; 345 inV as vapour_stream; 346 347 CONNECTIONS 348 feed.Outlet to t1.Inlet; 349 inL to t1.InletL; 350 inV to t1.InletV; 351 352 EQUATIONS 353 deltaP = (t1.InletV.P  t1.OutletV.P)/t1.hs; 354 355 SPECIFY 356 feed.Outlet.F = 113.4 * 'kmol/h'; 357 feed.Outlet.T = 291 * 'K'; 358 feed.Outlet.P = 1.66 * 'atm'; 359 feed.Outlet.z = [0.5, 0.5]; 360 361 inL.P = 165 * 'kPa'; 362 inL.T = 315 * 'K'; #310 * 'K'; 363 inL.F = 61.99 * 'kmol/h'; 364 inL.z = [0.1641, 0.8359];#[0.5, 0.5];# 365 366 inV.F = 201.25 * 'kmol/h'; 367 inV.P = 150 * 'kPa'; 368 inV.T = 315 * 'K'; #321 * 'K'; 369 inV.z = [0.0584, 0.9416];#[0.5, 0.5];# 370 371 t1.OutletV.P = 145 * 'kPa'; 372 373 SET 374 #Metal Pall Ring  nominal packing size 50 mm  Billet and Schultes, 1999. 375 t1.Q = 0 * 'kW'; 376 t1.Ap = 0.8 * 'm^2'; 377 t1.V = 0.8 * 'm^2' * 0.4 * 'm'; 378 t1.ds = 1.009 * 'm'; 379 t1.Cpo = 0.763; 380 t1.e = 0.951; 381 t1.a = 112.6 * 'm^2/m^3'; 382 t1.hs = 0.4 * 'm'; 383 t1.Qsio = 1; 384 385 INITIAL 386 t1.OutletL.T = 315 *'K'; 387 t1.ML = 0.25 * 'kmol'; 388 t1.OutletL.z(1) = 0.1641; 389 390 OPTIONS 391 DAESolver(File="sundials"); 392 TimeStep = 0.001; 393 TimeEnd = 0.1; 394 end 395 396 #to compare with tray_Test_2 397 FlowSheet packedStage_Test_2 398 PARAMETERS 399 PP as Plugin(Brief="Physical Properties", 400 Type="PP", 401 Components = [ "isobutane", "npentane", "propylene", 402 "benzene", "isobutene" ], 403 LiquidModel = "PR", 404 VapourModel = "PR" 405 ); 406 NComp as Integer; 407 408 SET 409 NComp = PP.NumberOfComponents; 410 411 VARIABLES 412 deltaP as Real (Unit='inH2O/m*mm/in'); 413 phiL as Real; 414 415 DEVICES 416 t1 as packedStage; 417 feed as source; 418 inL as liquid_stream; 419 inV as vapour_stream; 420 421 CONNECTIONS 422 feed.Outlet to t1.Inlet; 423 inL to t1.InletL; 424 inV to t1.InletV; 425 426 EQUATIONS 427 deltaP = (t1.InletV.P  t1.OutletV.P)/t1.hs; 428 phiL = t1.vL * t1.ML/t1.V; 429 430 SPECIFY 431 feed.Outlet.F = 0 * 'kmol/h'; 432 feed.Outlet.T = 300 * 'K'; 433 feed.Outlet.P = 1.66 * 'atm'; 434 feed.Outlet.z = [0.226, 0.425, 0.035, 0.025, 0.289]; 435 436 inL.F = 71.21 * 'kmol/h'; 437 inL.P = 2.22 * 'atm'; 438 inL.T = 297.6 * 'K'; 439 inL.z = [0.226, 0.425, 0.035, 0.025, 0.289]; 440 441 inV.F = 175.3 * 'kmol/h'; 442 inV.P = 2.3062 * 'atm'; 443 inV.T = 308.3 * 'K'; 444 inV.z = [0.265, 0.233, 0.150, 0.014, 0.338]; 445 446 #t1.OutletV.P = 2.305 * 'atm'; 447 #t1.OutletV.F = 165 * 'kmol/h'; 448 449 SET 450 #Metal Pall Ring  nominal packing size 50 mm  Billet and Schultes, 1999. 451 t1.Q = 0 * 'kW'; 452 #t1.Ap = 0.8 * 'm^2'; 453 #t1.V = 0.8 * 'm^2' * 0.4 * 'm'; 454 # t1.ds = 1.009 * 'm'; 455 t1.Cpo = 0.763; 456 t1.e = 0.951; 457 t1.a = 112.6 * 'm^2/m^3'; 458 # t1.hs = 0.4 * 'm'; 459 460 t1.V = 4 * 'ft^2' * 1 * 'ft'; 461 t1.Ap = 4 * 'ft^2'; 462 t1.hs = 1 * 'ft'; 463 t1.ds = 2.26 * 'ft'; 464 t1.Qsio = 100; 465 466 INITIAL 467 t1.OutletL.T = 290 *'K'; 468 t1.ML = 0.02 * 'kmol'; 469 t1.OutletL.z([1:4]) = [0.226, 0.425, 0.035, 0.025]; #[0.022, 0.0425, 0.0035, 0.95]; 470 471 OPTIONS 472 DAESolver(File="sundials"); 473 #InitialFile = "/home/paula/tray_Test_2.rlt"; 474 TimeStep = 0.01; 475 TimeEnd = 10; 476 end 477 478 FlowSheet packedStageAirWater_Test 479 PARAMETERS 480 PP as Plugin(Brief="Physical Properties", 481 Type="PP", 482 Components = [ "nitrogen", "water"], 483 LiquidModel = "IdealLiquid", 484 VapourModel = "Ideal" 485 ); 486 NComp as Integer; 487 488 SET 489 NComp = PP.NumberOfComponents; 490 491 VARIABLES 492 deltaP as Real (Unit='inH2O/m*mm/in'); 493 phiL as Real; 494 495 DEVICES 496 t1 as packedStage_AirWater; 497 feed as source; 498 inL as liquid_stream; 499 inV as vapour_stream; 500 501 CONNECTIONS 502 feed.Outlet to t1.Inlet; 503 inL to t1.InletL; 504 inV to t1.InletV; 505 506 EQUATIONS 507 deltaP = (t1.InletV.P  t1.OutletV.P)/t1.hs; 508 phiL = t1.vL * t1.ML/t1.V; 509 510 SPECIFY 511 feed.Outlet.F = 0 * 'kmol/h'; 512 feed.Outlet.T = 300 * 'K'; 513 feed.Outlet.P = 1 * 'atm'; 514 feed.Outlet.z = [0 1]; 515 516 inL.F = 71.21 * 'kmol/h'; 517 inL.P = 1 * 'atm'; 518 inL.T = 293 * 'K'; 519 inL.z = [0 1]; 520 521 inV.F = 175.3 * 'kmol/h'; 522 inV.P = 1.06 * 'atm'; 523 inV.T = 295 * 'K'; 524 inV.z = [1 0]; 525 526 t1.OutletV.P = 1.05 * 'atm'; 527 528 SET 529 #Metal Bialecki Ring  nominal packing size 25 mm  Billet and Schultes, 1999. 530 t1.Cpo = 0.891; 531 t1.Ch = 0.692; 532 t1.e = 0.956; 533 t1.a = 210 * 'm^2/m^3'; 534 535 t1.Q = 0 * 'kW'; 536 t1.V = 0.018 * 'm^2' * 1.4 * 'm'; 537 t1.Ap = 0.018 * 'm^2'; 538 t1.hs = 1.4 * 'm'; 539 t1.ds = 0.15 * 'm'; 540 541 INITIAL 542 t1.OutletL.T = 294 *'K'; 543 t1.ML = 0.0002 * 'kmol'; 544 t1.OutletL.z(1) = 0.01; 545 546 OPTIONS 547 #InitialFile = "/home/paula/packedStageAirWater_Test.rlt"; 548 DAESolver(File="dassl"); 549 TimeStep = 10; 550 TimeEnd = 100; 551 end
Note: See TracChangeset
for help on using the changeset viewer.