Changeset 453 for branches/packed
- Timestamp:
- Feb 7, 2008, 5:26:42 PM (15 years ago)
- Location:
- branches/packed
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/packed/eml/stage_separators/column.mso
r450 r453 1306 1306 stage([top+topdown:topdown:bot]).OutletV to stage([top:topdown:bot-topdown]).InletV; 1307 1307 stage([top:topdown:bot-topdown]).OutletL to stage([top+topdown:topdown:bot]).InletL; 1308 1309 1310 end 1311 1308 end 1309 1310 Model PackedDistillation_kettle_cond 1311 PARAMETERS 1312 outer PP as Plugin(Brief = "External Physical Properties", Type="PP"); 1313 outer NComp as Integer; 1314 NStages as Integer(Brief="Number of stages", Default=2); 1315 topdown as Integer(Brief="Trays counting (1=top-down, -1=bottom-up)", Default=1); 1316 top as Integer(Brief="Number of top tray"); 1317 bot as Integer(Brief="Number of bottom tray"); 1318 1319 SET 1320 top = (NStages-1)*(1-topdown)/2+1; 1321 bot = NStages/top; 1322 1323 VARIABLES 1324 stage(NStages) as packedStage; 1325 cond as condenser; 1326 reb as reboiler; 1327 sptop as splitter; 1328 pump1 as pump; 1329 1330 CONNECTIONS 1331 #vapor 1332 reb.OutletV to stage(bot).InletV; 1333 stage([top+topdown:topdown:bot]).OutletV to stage([top:topdown:bot-topdown]).InletV; 1334 stage(top).OutletV to cond.InletV; 1335 1336 #liquid 1337 cond.OutletL to sptop.Inlet; 1338 sptop.Outlet2 to pump1.Inlet; 1339 pump1.Outlet to stage(top).InletL; 1340 stage([top:topdown:bot-topdown]).OutletL to stage([top+topdown:topdown:bot]).InletL; 1341 stage(bot).OutletL to reb.InletL; 1342 1343 EQUATIONS 1344 reb.OutletV.F = stage(bot).OutletV.F; 1345 end -
branches/packed/eml/stage_separators/tray.mso
r451 r453 286 286 d as length (Brief="size of an element of packing"); 287 287 #h as length (Brief="Height of packing"); 288 C as Real (Brief="Constant for resitance factor equation"); 289 Cp as Real (Brief="Constant for resitance at loading point factor equation"); 288 C as Real (Brief="Constant for resitance factor equation"); #Billet pp 80 (Cf+Cfl)/2 289 Cp as Real (Brief="Constant for resitance at loading point factor equation"); # Billet pp94 290 290 Mw(NComp) as molweight (Brief = "Component Mol Weight"); 291 291 … … 329 329 uV * Across = OutletV.F * vV; 330 330 331 "Coefficient of Resistance" 332 ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m'; 331 # "Coefficient of Resistance" 332 # ksi * C^2 * (uL/uV * sqrt(rhoV/rhoL) * (niL/niV)^5.8)^3 = g/1*'s^2/m'; 333 #ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) * 333 334 ksil = Cp * (exp(uL*rhoL/a/niL/200)*(hL/hLs)^0.3) * (64/Rev+(1.8/Rev)) * 334 ((e-hL)/e);#^1.5;335 336 "Liquid holdup and Liquid flow"337 a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2));338 339 "Pressure drop and Vapor flow"340 (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3);335 ((e-hL)/e); 336 337 # "Liquid holdup and Liquid flow" 338 # a^2 * niL * uL = hL^1 *(g*rhoL/3 - ksi*a*rhoV*uV^2/(4*hL*(e-hL)^2)); 339 340 # "Pressure drop and Vapor flow" 341 # (InletV.P - OutletV.P)/hs = ksil *(a/2 + 2/ds)*(uV^2*rhoV/(e-hL)^3); 341 342 342 343 "Vapor Reynolds number" … … 344 345 345 346 "Theoretical Liquid Holdup" 346 hLs = (12*a^2*niL*uL/g/rhoL)^0.333; 347 #hLs = (12*a^2*niL*uL/g/rhoL)^0.333; 348 hLs = (12*a^2*niL*uL/g/rhoL); 347 349 end 348 349 350 350 351 #*------------------------------------------------------------------- -
branches/packed/sample/stage_separators/sample_column.mso
r448 r453 582 582 end 583 583 584 FlowSheet Packed_kettle_cond_Test 585 PARAMETERS 586 PP as Plugin(Brief="Physical Properties", 587 Type="PP", 588 Components = [ "methanol", "water"], 589 LiquidModel = "IdealLiquid", 590 VapourModel = "Ideal" 591 ); 592 NComp as Integer; 593 594 VARIABLES 595 Qc as energy_source (Brief="Heat rate removed from condenser"); 596 Qr as energy_source (Brief="Heat rate supplied to reboiler"); 597 598 SET 599 NComp = PP.NumberOfComponents; 600 601 DEVICES 602 col as PackedDistillation_kettle_cond; 603 feed as source; 604 zero as stream; 605 606 CONNECTIONS 607 feed.Outlet to col.reb.Inlet; 608 zero to col.stage([1:col.NStages]).Inlet; 609 Qc.OutletQ to col.cond.InletQ; 610 Qr.OutletQ to col.reb.InletQ; 611 612 SPECIFY 613 feed.Outlet.F = 0.98 * 'mol/min'; 614 feed.Outlet.T = (50+273.15) * 'K'; 615 feed.Outlet.P = 0.91 * 'kPa'; 616 feed.Outlet.z = [0.16, 0.84]; 617 618 zero.F = 0 * 'kmol/h'; 619 zero.T = 300 * 'K'; 620 zero.P = 1 * 'atm'; 621 zero.z = 1/NComp; 622 zero.v = 0; 623 zero.h = 0 * 'J/mol'; 624 625 col.sptop.frac = 0.4; 626 col.sptop.Outlet1.F = 0.10 * 'mol/min'; 627 col.reb.OutletL.F = 0.88 * 'mol/min'; 628 col.cond.OutletV.F = 0 * 'kmol/h'; 629 Qr.OutletQ.Q = 700 * 'cal/min'; 630 Qc.OutletQ.Q = -700 * 'cal/min'; 631 col.pump1.dP = 0.1 * 'atm'; 632 col.stage.Emv = 1; 633 col.stage.hs = 100 * 'mm'; 634 635 SET 636 col.NStages = 10; 637 col.cond.V = 1 * 'l'; 638 col.cond.Across = 100 * 'cm^2'; 639 col.reb.V = 2 * 'l'; 640 col.reb.Across = 200 * 'cm^2'; 641 col.stage.V = 576.4 * 'cm^3'; 642 col.stage.Q = 0 * 'kW'; 643 col.stage.Ap = 56.74 * 'cm^2'; 644 col.stage.ds = 85 * 'mm'; 645 col.stage.d = 20 * 'mm'; 646 col.stage.C = 2.18; 647 col.stage.Cp = 1.329; 648 col.stage.e = 0.662; 649 col.stage.a = 185.4 * 'm^2/m^3'; 650 651 INITIAL 652 # condenser 653 col.cond.OutletL.T = (63.5+273.15) *'K'; 654 col.cond.Level = 2 * 'cm'; 655 col.cond.OutletL.z(1) = 0.16; 656 657 # reboiler 658 col.reb.OutletL.T = (82+273.15) *'K'; 659 col.reb.Level = 2 * 'cm'; 660 col.reb.OutletL.z(1) = 0.16; 661 662 # column stages 663 col.stage.OutletL.T = [(63.5+273.15):((82+273.15)-(63.5+273.15))/(col.NStages-1):(82+273.15)] * 'K'; 664 col.stage.Level = 1 * 'cm'; 665 col.stage.OutletL.z(1) = 0.16; 666 667 OPTIONS 668 TimeStep = 0.1; 669 TimeEnd = 50; 670 #GuessFile = "/home/paula/SectionColumn_Test_with8tray.rlt"; 671 #GuessFile="Distillation_kettle_cond_Test.rlt"; 672 #Dynamic = false; 673 end -
branches/packed/sample/stage_separators/sample_tray.mso
r451 r453 143 143 144 144 t1.Emv = 1; 145 t1.hs = 0. 075 * 'm';146 #t1.OutletV.F = 147.1 * 'kmol/h';147 #t1.OutletL.F = 229.5 * 'kmol/h';148 #t1.ksi = 0.805;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 149 150 150 SET … … 152 152 t1.Q = 0 * 'kW'; 153 153 t1.Across = 0.8 * 'm^2'; 154 t1.V = 0.8 * 'm^2' * 0. 075* 'm';154 t1.V = 0.8 * 'm^2' * 0.4 * 'm'; 155 155 t1.ds = 1.009 * 'm'; 156 156 t1.d = 50 * 'mm'; … … 167 167 OPTIONS 168 168 #InitialFile = "/home/paula/tray_Test.rlt"; 169 InitialFile = "/home/paula/packedStage_Test.rlt"; 169 170 #DAESolver(File="dassl"); 170 TimeStep = 0.1;171 TimeEnd = 10;171 TimeStep = 2; 172 TimeEnd = 50; 172 173 end
Note: See TracChangeset
for help on using the changeset viewer.