Changeset 517 for branches/rate/eml
 Timestamp:
 May 19, 2008, 11:29:24 AM (15 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/rate/eml/stage_separators/tray.mso
r516 r517 119 119 == Initial == 120 120 * the plate temperature (OutletL.T) 121 * the liquid height (Level) OR the liquid flow OutletL.F121 * the liquid height (Level) OR the liquid holdup (ML) 122 122 * (NoComps  1) OutletL compositions 123 123 "; … … 494 494 495 495 VARIABLES 496 NL(NComp) as flow_mol 497 NV(NComp) as flow_mol 496 NL(NComp) as flow_mol_delta (Brief = "Stream Molar Rate on Liquid Phase"); 497 NV(NComp) as flow_mol_delta (Brief = "Stream Molar Rate on Vapour Phase"); 498 498 T as temperature (Brief = "Stream Temperature"); 499 499 P as pressure (Brief = "Stream Pressure"); … … 519 519 end 520 520 521 Model trayRate 521 Model trayRateBasic 522 522 ATTRIBUTES 523 523 Pallete = false; 524 524 Icon = "icon/Tray"; 525 Brief = "Basic equations of a tray column model.";525 Brief = "Basic equations of a tray rate column model."; 526 526 Info = 527 527 "This model contains only the main equations of a column tray nonequilibrium model without … … 561 561 Level as length (Brief="Height of clear liquid on plate"); 562 562 interf as interface; 563 563 564 SET 565 NC1=NComp1; 566 564 567 EQUATIONS 565 568 "Component Molar Balance" … … 622 625 interf.NL(1:NC1)=interf.a*sumt(interf.kL*(interf.x(1:NC1)OutletL.z(1:NC1)))/vL+ 623 626 OutletL.z(1:NC1)*sum(interf.NL); 627 628 # interf.NL(1:NC1)=0.01*'kmol/s'; 624 629 625 630 interf.NV(1:NC1)=interf.a*sumt(interf.kV*(OutletV.z(1:NC1)interf.y(1:NC1)))/vV+ … … 629 634 OutletV.P = OutletL.P; 630 635 interf.P=OutletL.P; 631 632 SET633 NC1=NComp1;634 636 end 635 637 638 Model trayRate as trayRateBasic 639 ATTRIBUTES 640 Pallete = false; 641 Icon = "icon/Tray"; 642 Brief = "Complete rate model of a column tray."; 643 Info = 644 "== Specify == 645 * the Feed stream 646 * the Liquid inlet stream 647 * the Vapour inlet stream 648 * the Vapour outlet flow (OutletV.F) 649 650 == Initial == 651 * the plate temperature of both phases (OutletL.T and OutletV.T) 652 * the liquid height (Level) OR the liquid flow holdup (ML) 653 * the vapor holdup (MV) 654 * (NoComps  1) OutletL compositions 655 "; 656 657 PARAMETERS 658 Ah as area (Brief="Total holes area"); 659 lw as length (Brief="Weir length"); 660 g as acceleration (Default=9.81); 661 hw as length (Brief="Weir height"); 662 beta as fraction (Brief="Aeration fraction"); 663 alfa as fraction (Brief="Dry pressure drop coefficient"); 664 665 VapourFlow as Switcher(Valid = ["on", "off"], Default = "on"); 666 LiquidFlow as Switcher(Valid = ["on", "off"], Default = "on"); 667 668 VARIABLES 669 rhoL as dens_mass; 670 rhoV as dens_mass; 671 672 EQUATIONS 673 "Liquid Density" 674 rhoL = PP.LiquidDensity(OutletL.T, OutletL.P, OutletL.z); 675 "Vapour Density" 676 rhoV = PP.VapourDensity(InletV.T, InletV.P, InletV.z); 677 678 switch LiquidFlow 679 case "on": 680 "Francis Equation" 681 # OutletL.F*vL = 1.84*'m^0.5/s'*lw*((Level(beta*hw))/(beta))^1.5; 682 OutletL.F*vL = 1.84*'1/s'*lw*((Level(beta*hw))/(beta))^2; 683 when Level < (beta * hw) switchto "off"; 684 685 case "off": 686 "Low level" 687 OutletL.F = 0 * 'mol/h'; 688 when Level > (beta * hw) + 1e6*'m' switchto "on"; 689 end 690 691 switch VapourFlow 692 case "on": 693 InletV.F*vV = sqrt((InletV.P  OutletV.P)/(rhoV*alfa))*Ah; 694 when InletV.F < 1e6 * 'kmol/h' switchto "off"; 695 696 case "off": 697 InletV.F = 0 * 'mol/s'; 698 when InletV.P > OutletV.P + Level*g*rhoL + 1e1 * 'atm' switchto "on"; 699 end 700 end
Note: See TracChangeset
for help on using the changeset viewer.