Changeset 458 for branches/packed
 Timestamp:
 Feb 14, 2008, 5:09:53 PM (15 years ago)
 Location:
 branches/packed
 Files:

 3 edited
Legend:
 Unmodified
 Added
 Removed

branches/packed/eml/stage_separators/tray.mso
r455 r458 366 366 uL as velocity (Brief="volume flow rate of liquid, m^3/m^2/s", Default = 0.007); 367 367 uV as velocity (Brief="volume flow rate of vapor, m^3/m^2/s", Default = 1.14); 368 dp as length (Brief="Particle diameter" );369 K as Real (Brief="Wall factor");370 Rev as Real (Brief="Reynolds number of the vapor stream" );371 Qsio as Real (Brief="Resistance coefficient" );368 dp as length (Brief="Particle diameter", Default=1e3); 369 invK as Real (Brief="Wall factor"); 370 Rev as Real (Brief="Reynolds number of the vapor stream", Lower = 0, Default=100); 371 Qsio as Real (Brief="Resistance coefficient", Lower = 0); 372 372 373 373 SET … … 393 393 394 394 "Particle diameter" 395 dp = 6 * ( (1e)/a);395 dp = 6 * (1e)/a; 396 396 397 397 "Wall Factor" 398 1 = K *(1 + (2*dp/(3*ds*(1e))));398 invK = (1 + (2*dp/(3*ds*(1e)))); 399 399 400 400 "Reynolds number of the vapor stream" 401 Rev = K*dp*uV*rhoV / (miV*(1e)); 402 403 #"Resistance Coefficient" 404 #Qsio = Cpo * (64/Rev + 1.8/Rev^0.08); 405 #Qsio = Cpo * (64/Rev);# + 1.8/Rev^0.08); 401 Rev*invK = dp*uV*rhoV / (miV*(1e)); 402 403 "Resistance Coefficient" 404 Qsio = Cpo * (64/Rev + 1.8/Rev^0.08); 406 405 407 406 "Pressure drop and Vapor flow" 408 (InletV.P  OutletV.P)/hs * 'm'/'Pa' = Qsio*a*uV^2*rhoV / (2*e^3*K) * 's^2'*'m^2'/'kg';407 (InletV.P  OutletV.P)/hs = Qsio*a*uV^2*rhoV*invK / (2*e^3); 409 408 410 409 "Efficiency" … … 414 413 #* 415 414 * Model of a tray with reaction 416 * *#415 **# 417 416 Model trayReact 418 417 ATTRIBUTES 
branches/packed/sample/stage_separators/sample_column.mso
r455 r458 99 99 100 100 OPTIONS 101 TimeStep = 10 ;101 TimeStep = 100; 102 102 TimeEnd = 1000; 103 103 … … 277 277 278 278 OPTIONS 279 TimeStep = 0.1;279 TimeStep = 5; 280 280 TimeEnd = 50; 281 281 #GuessFile="Distillation_kettle_cond_Test.rlt"; … … 508 508 PP as Plugin(Brief="Physical Properties", 509 509 Type="PP", 510 Components = [ "isobutene", "npentane", "propylene", "benzene", "isobutane"], 510 Components = [ "isobutane", "npentane", "propylene", 511 "benzene", "isobutene" ], 511 512 LiquidModel = "PR", 512 513 VapourModel = "PR" … … 519 520 DEVICES 520 521 sec as Section_Column_Packed; 521 feed as source;522 feed as liquid_stream; 522 523 reb as vapour_stream; 523 524 cond as liquid_stream; … … 525 526 526 527 CONNECTIONS 527 feed.Outlet to sec.stage(10).Inlet; 528 529 zero to sec.stage(1:9).Inlet; 530 zero to sec.stage([11:sec.NStages]).Inlet; 531 532 reb to sec.stage(sec.NStages).InletV; 528 feed to sec.stage(5).Inlet; 529 zero to sec.stage([1:4]).Inlet; 530 zero to sec.stage([6:sec.NStages]).Inlet; 531 reb to sec.stage(8).InletV; 533 532 cond to sec.stage(1).InletL; 534 533 535 VARIABLES536 deltaP_TOP as Real (Unit='atm/m');537 deltaP_BOT as Real (Unit='atm/m');538 539 EQUATIONS540 deltaP_TOP = (sec.stage(1).InletV.P  sec.stage(1).OutletV.P)/sec.stage(1).hs;541 deltaP_BOT = (sec.stage(sec.NStages).InletV.P  sec.stage(sec.NStages).OutletV.P)/sec.stage(sec.NStages).hs;542 534 543 535 SPECIFY 544 feed. Outlet.F = 113.4 * 'kmol/h';545 feed. Outlet.T = 291 * 'K';546 feed. Outlet.P = 168.3 * 'kPa';547 feed. Outlet.z = 1/NComp;536 feed.F = 113.4 * 'kmol/h'; 537 feed.T = 291 * 'K'; 538 feed.P = 1.66 * 'atm'; 539 feed.z = [0.2, 0.2, 0.2, 0.2, 0.2]; 548 540 549 541 zero.F = 0 * 'kmol/h'; 550 542 zero.T = 300 * 'K'; 551 543 zero.P = 1 * 'atm'; 552 zero.z = 1/NComp;544 zero.z = [0.2, 0.2, 0.2, 0.2, 0.2]; 553 545 zero.v = 0; 554 546 zero.h = 0 * 'J/mol'; 555 547 556 548 cond.F = 85 * 'kmol/h'; 557 cond.P = 2. 2* 'atm';558 cond.T = 2 73* 'K';559 cond.z = [0. 26682, 0.19950, 0.26682, 0.000051, 0.26680];560 561 reb.F = 1 43.4* 'kmol/h';562 reb.P = 5* 'atm';563 reb.T = 3 61* 'K';564 reb.z = [0. 0001611, 0.51834, 0.00000076, 0.48115, 0.00035];565 566 #sec.stage.Qsio = 1;567 568 SET 569 sec.H = 4 * 'm';570 sec.NStages = 16;549 cond.P = 2.33 * 'atm'; 550 cond.T = 283.5 * 'K'; 551 cond.z = [0.599, 0.044, 0.035, 0.007, 0.315]; 552 553 reb.F = 137.57 * 'kmol/h'; 554 reb.P = 2.46 * 'atm'; 555 reb.T = 325 * 'K'; 556 reb.z = [0.16, 0.542, 0.013, 0.008, 0.277]; 557 558 #sec.stage.Qsio = 0.8; 559 560 SET 561 sec.H = 6 * 'ft'; 562 sec.NStages = 8; 571 563 sec.stage.Q = 0 * 'kW'; 572 sec.stage.Ap = 1.5 * 'm^2';564 sec.stage.Ap = 3.94 * 'ft^2'; 573 565 sec.stage.ds = 1.009 * 'm'; 574 566 sec.stage.Cpo = 0.763; … … 577 569 578 570 INITIAL 579 sec.stage.OutletL.T = [273:(360273)/(sec.NStages1):360] *'K';571 sec.stage.OutletL.T =[283:(325283)/(sec.NStages1):325] *'K'; 580 572 sec.stage.ML = 0.5 * 'kmol'; 581 sec.stage.OutletL.z([1:4]) = 1/NComp;582 573 sec.stage.OutletL.z([1:4]) = [0.2, 0.2, 0.2, 0.2]; 574 583 575 OPTIONS 576 #InitialFile = "/home/paula/SectionColumn_Test_with8tray.rlt"; 584 577 DAESolver(File="dassl"); 585 578 TimeStep = 0.01; … … 637 630 col.stage.Qsio = 0.05; 638 631 639 EQUATIONS 640 col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q); 641 642 SET 632 # EQUATIONS 633 # col.reb.OutletV.F = col.K * sqrt(Qr.OutletQ.Q); 634 635 SET 636 col.H = 1 * 'm'; 643 637 col.NStages = 5; 638 644 639 col.cond.V = 1 * 'l'; 645 640 col.cond.Across = 100 * 'cm^2'; 646 col.K = 0.2 * 'mol*min^0.5/kg^0.5/m'; 641 # col.K = 0.2 * 'mol*min^0.5/kg^0.5/m'; 642 647 643 col.reb.V = 2 * 'l'; 648 644 col.reb.Across = 200 * 'cm^2'; 645 649 646 col.stage.Q = 0 * 'kW'; 650 647 col.stage.Ap = 56.74 * 'cm^2'; … … 652 649 col.stage.e = 0.662; 653 650 col.stage.a = 185.4 * 'm^2/m^3'; 654 col.H = 1 * 'm';655 651 col.stage.Cpo = 0.763; 656 652 
branches/packed/sample/stage_separators/sample_tray.mso
r455 r458 96 96 end 97 97 98 FlowSheet packedStage_Test99 PARAMETERS100 PP as Plugin(Brief="Physical Properties",101 Type="PP",102 Components = [ "npentane", "benzene"],103 LiquidModel = "PR",104 VapourModel = "PR"105 );106 NComp as Integer;107 108 SET109 NComp = PP.NumberOfComponents;110 111 VARIABLES112 deltaP as Real (Unit='atm/m');113 114 DEVICES115 t1 as packedStage_Billet;116 feed as source;117 inL as liquid_stream;118 inV as vapour_stream;119 120 CONNECTIONS121 feed.Outlet to t1.Inlet;122 inL to t1.InletL;123 inV to t1.InletV;124 125 EQUATIONS126 deltaP = (t1.InletV.P  t1.OutletV.P)/t1.hs;127 128 SPECIFY129 feed.Outlet.F = 113.4 * 'kmol/h';130 feed.Outlet.T = 291 * 'K';131 feed.Outlet.P = 1.66 * 'atm';132 feed.Outlet.z = [0.5, 0.5];133 134 inL.F = 61.99 * 'kmol/h';135 inL.P = 1.63 * 'atm';136 inL.T = 310 * 'K';137 inL.z = [0.1641, 0.8359];138 139 inV.F = 201.25 * 'kmol/h';140 inV.P = 1.48 * 'atm';141 inV.T = 321 * 'K';142 inV.z = [0.0584, 0.9416];143 144 t1.Emv = 1;145 t1.hs = 0.4 * 'm';#0.075 * 'm';146 t1.OutletV.F = 147.1 * 'kmol/h';147 # t1.OutletL.F = 229.5 * 'kmol/h';148 #t1.ksi = 0.805;149 150 SET151 #Ceramic Pall Ring  nominal packing size 50 mm152 t1.Q = 0 * 'kW';153 t1.Across = 0.8 * 'm^2';154 t1.V = 0.8 * 'm^2' * 0.4 * 'm';155 t1.ds = 1.009 * 'm';156 t1.d = 50 * 'mm';157 t1.C = 2.37;158 t1.Cp = 0.662;159 t1.e = 0.78;160 t1.a = 120 * 'm^2/m^3';161 162 INITIAL163 t1.OutletL.T = 350 *'K';164 t1.ML = 0.5 * 'mol';165 t1.OutletL.z(1) = 0.2;166 167 OPTIONS168 #InitialFile = "/home/paula/tray_Test.rlt";169 #InitialFile = "/home/paula/packedStage_Test.rlt";170 TimeStep = 2;171 TimeEnd = 50;172 end173 174 98 FlowSheet packedStage_BilletSchultes_Test 175 99 PARAMETERS 176 100 PP as Plugin(Brief="Physical Properties", 177 101 Type="PP", 178 Components = [ "npentane", "benzene"], 102 Components = [ "isobutane", "npentane", "propylene", 103 "benzene", "isobutene" ], 179 104 LiquidModel = "PR", 180 105 VapourModel = "PR" … … 203 128 204 129 SPECIFY 205 feed.Outlet.F = 113.4* 'kmol/h';206 feed.Outlet.T = 291* 'K';207 feed.Outlet.P = 1 .66* 'atm';208 feed.Outlet.z = [0.5, 0.5];130 feed.Outlet.F = 0 * 'kmol/h'; 131 feed.Outlet.T = 300 * 'K'; 132 feed.Outlet.P = 1 * 'atm'; 133 feed.Outlet.z = 1/NComp; 209 134 210 inL.F = 61.99* 'kmol/h';211 inL.P = 1.63* 'atm';212 inL.T = 310* 'K';213 inL.z = [0. 1641, 0.8359];135 inL.F = 7121 * 'kmol/h'; 136 inL.P = 2.22 * 'atm'; 137 inL.T = 297.6 * 'K'; 138 inL.z = [0.226, 0.425, 0.035, 0.025, 0.289]; 214 139 215 inV.F = 201.25 * 'kmol/h'; 216 inV.P = 1.48 * 'atm'; 217 inV.T = 321 * 'K'; 218 inV.z = [0.0584, 0.9416]; 140 inV.F = 175.3 * 'kmol/h'; 141 inV.P = 2.3062 * 'atm'; 142 inV.T = 308.3 * 'K'; 143 inV.z = [0.265, 0.233, 0.150, 0.014, 0.338]; 144 145 #t1.OutletV.F = 177.9 * 'kmol/h'; 146 #t1.Qsio = 0.61; 219 147 220 t1.hs = 0.4 * 'm';#0.075 * 'm';221 # t1.OutletV.F = 147.1 * 'kmol/h';222 # t1.OutletL.F = 229.5 * 'kmol/h';223 # t1.Qsio = 2.2;224 225 148 SET 226 149 #Metal Pall Ring  nominal packing size 50 mm  Billet and Schultes, 1999. … … 232 155 t1.e = 0.951; 233 156 t1.a = 112.6 * 'm^2/m^3'; 157 t1.hs = 0.4 * 'm';#0.075 * 'm'; 234 158 235 159 INITIAL 236 t1.OutletL.T = 310 *'K';237 t1.ML = 0. 5 * 'mol';238 t1.OutletL.z( 1) = 0.05;160 t1.OutletL.T = 290 *'K'; 161 t1.ML = 0.2 * 'kmol'; 162 t1.OutletL.z([1:4]) = [0.286, 0.203, 0.148, 0.006]; 239 163 240 164 OPTIONS 241 InitialFile = "packedStage_BilletSchultes_Test.rlt"; 165 #NLASolver(File="nlasolver"); 166 #InitialFile = "packedStage_BilletSchultes_Test.rlt"; 242 167 TimeStep = 2; 243 168 TimeEnd = 50;
Note: See TracChangeset
for help on using the changeset viewer.