- Timestamp:
- Jul 28, 2007, 9:41:04 PM (16 years ago)
- Location:
- trunk/eml
- Files:
-
- 1 deleted
- 27 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/eml/electrical/electrical.mso
r303 r325 32 32 end 33 33 34 Model electrical_basic35 ATTRIBUTES36 Pallete = false;37 Brief = "Basic model for electrical devices with one input and one output wire.";38 39 VARIABLES40 in inlet as wire;41 out outlet as wire;42 end43 34 44 Model electrical as electrical_basic 45 ATTRIBUTES 46 Pallete = false; 47 Brief = "Basic electrical device"; 48 Info = 49 "Model for an electrical device in which the inlet current 50 is equal to the outlet one 51 "; 52 53 EQUATIONS 54 outlet.i = inlet.i; 55 end 56 57 Model Resistor as electrical 35 Model Resistor 58 36 ATTRIBUTES 59 37 Pallete = true; … … 63 41 PARAMETERS 64 42 R as resistance; 43 44 VARIABLES 45 in inlet as wire (Brief = "Inlet", PosX=0.3831, PosY=0); 46 out outlet as wire (Brief = "Outlet", PosX=0.3529, PosY=1); 47 65 48 EQUATIONS 66 49 inlet.V - outlet.V = R * outlet.i; 50 outlet.i = inlet.i; 67 51 end 68 52 69 Model Capacitor as electrical53 Model Capacitor 70 54 ATTRIBUTES 71 55 Pallete = true; … … 75 59 PARAMETERS 76 60 C as capacitance; 61 77 62 VARIABLES 78 63 q as charge; 64 in inlet as wire (Brief = "Inlet", PosX=0.3978, PosY=0); 65 out outlet as wire (Brief = "Outlet", PosX=0.3965, PosY=1); 66 79 67 EQUATIONS 80 68 diff(q) = inlet.i; 81 69 inlet.V - outlet.V = (1/C) * q; 70 outlet.i = inlet.i; 82 71 end 83 72 84 Model Indutor as electrical73 Model Indutor 85 74 ATTRIBUTES 86 75 Pallete = true; … … 90 79 PARAMETERS 91 80 L as indutance; 81 82 VARIABLES 83 in inlet as wire (Brief = "Inlet", PosX=0.4638, PosY=0); 84 out outlet as wire (Brief = "Outlet", PosX=0.4638, PosY=1); 85 92 86 EQUATIONS 93 87 inlet.V - outlet.V = L * diff(inlet.i); 88 outlet.i = inlet.i; 94 89 end 95 90 96 Model Supply as electrical_basic91 Model Supply 97 92 ATTRIBUTES 98 93 Pallete = true; … … 103 98 V as voltage; 104 99 V0 as voltage(Default = 0); 100 101 VARIABLES 102 in inlet as wire (Brief = "Inlet", PosX=0.3923, PosY=0); 103 out outlet as wire (Brief = "Outlet", PosX=0.3984, PosY=1); 104 105 105 EQUATIONS 106 106 outlet.V = V0; -
trunk/eml/heat_exchangers/DoublePipe.mso
r301 r325 29 29 PARAMETERS 30 30 31 outer PP 32 outer NComp 31 outer PP as Plugin (Brief="External Physical Properties", Type="PP"); 32 outer NComp as Integer (Brief="Number of Components"); 33 33 34 34 M(NComp) as molweight (Brief="Component Mol Weight"); 35 35 36 HotSide 36 HotSide as Switcher (Brief="Flag for Fluid Alocation ",Valid=["outer","inner"],Default="outer"); 37 37 innerFlowRegime as Switcher (Brief="Inner Flow Regime ",Valid=["laminar","transition","turbulent"],Default="laminar"); 38 38 outerFlowRegime as Switcher (Brief="Outer Flow Regime ",Valid=["laminar","transition","turbulent"],Default="laminar"); … … 40 40 InnerLaminarCorrelation as Switcher (Brief="Heat Transfer Correlation in Laminar Flow for the Inner Side",Valid=["Hausen","Schlunder"],Default="Hausen"); 41 41 InnerTransitionCorrelation as Switcher (Brief="Heat Transfer Correlation in Transition Flow for the Inner Side",Valid=["Gnielinski","ESDU"],Default="Gnielinski"); 42 InnerTurbulentCorrelation as Switcher(Brief="Heat Transfer Correlation in Turbulent Flow for the Inner Side",Valid=["Petukhov","SiederTate"],Default="Petukhov");42 InnerTurbulentCorrelation as Switcher (Brief="Heat Transfer Correlation in Turbulent Flow for the Inner Side",Valid=["Petukhov","SiederTate"],Default="Petukhov"); 43 43 44 44 OuterLaminarCorrelation as Switcher (Brief="Heat Transfer Correlation in Laminar Flow for the Outer Side",Valid=["Hausen","Schlunder"],Default="Hausen"); 45 45 OuterTransitionCorrelation as Switcher (Brief="Heat Transfer Correlation in Transition Flow for the OuterSide",Valid=["Gnielinski","ESDU"],Default="Gnielinski"); 46 OuterTurbulentCorrelation as Switcher (Brief="Heat Transfer Correlation in Turbulent Flow for the Outer Side",Valid=["Petukhov","SiederTate"],Default="Petukhov");47 48 Pi 49 DoInner 46 OuterTurbulentCorrelation as Switcher (Brief="Heat Transfer Correlation in Turbulent Flow for the Outer Side",Valid=["Petukhov","SiederTate"],Default="Petukhov"); 47 48 Pi as constant (Brief="Pi Number",Default=3.14159265); 49 DoInner as length (Brief="Outside Diameter of Inner Pipe",Lower=1e-6); 50 50 DiInner as length (Brief="Inside Diameter of Inner Pipe",Lower=1e-10); 51 51 DiOuter as length (Brief="Inside Diameter of Outer pipe",Lower=1e-10); 52 Lpipe 53 Kwall 54 Rfi 52 Lpipe as length (Brief="Effective Tube Length",Lower=0.1); 53 Kwall as conductivity (Brief="Tube Wall Material Thermal Conductivity",Default=1.0); 54 Rfi as positive (Brief="Inside Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 55 55 Rfo as positive (Brief="Outside Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 56 56 57 57 VARIABLES 58 58 59 in InletInner as stream (Brief="Inlet Inner Stream");60 in InletOuter as stream (Brief="Inlet Outer Stream");61 out OutletInner as streamPH (Brief="Outlet Inner Stream");62 out OutletOuter as streamPH (Brief="Outlet Outer Stream");63 64 Details as Details_Main(Brief="Some Details in the Heat Exchanger");65 Inner as Main_DoublePipe 66 Outer as Main_DoublePipe 59 in InletInner as stream (Brief="Inlet Inner Stream", PosX=0, PosY=0.5225); 60 in InletOuter as stream (Brief="Inlet Outer Stream", PosX=0.2805, PosY=0); 61 out OutletInner as streamPH (Brief="Outlet Inner Stream", PosX=1, PosY=0.5225); 62 out OutletOuter as streamPH (Brief="Outlet Outer Stream", PosX=0.7264, PosY=1); 63 64 Details as Details_Main (Brief="Some Details in the Heat Exchanger"); 65 Inner as Main_DoublePipe (Brief="Inner Side of the Heat Exchanger"); 66 Outer as Main_DoublePipe (Brief="Outer Side of the Heat Exchanger"); 67 67 68 68 SET -
trunk/eml/heat_exchangers/HeatExchangerDetailed.mso
r308 r325 28 28 PARAMETERS 29 29 30 HotSide 30 HotSide as Switcher (Brief="Hot Side in the Exchanger",Valid=["shell","tubes"],Default="shell"); 31 31 ShellType as Switcher (Brief="TEMA Designation",Valid=["Eshell","Fshell"],Default="Eshell"); 32 32 33 33 VARIABLES 34 34 35 in InletTube as stream (Brief="Inlet Tube Stream");36 out OutletTube as streamPH (Brief="Outlet Tube Stream");37 in InletShell as stream (Brief="Inlet Shell Stream");38 out OutletShell as streamPH (Brief="Outlet Shell Stream");39 40 Details as Details_Main(Brief="Details in Heat Exchanger");41 Tubes 42 Shell 43 Baffles as Baffles_Main(Brief="Baffles");35 in InletTube as stream (Brief="Inlet Tube Stream", PosX=0, PosY=0.4915); 36 out OutletTube as streamPH (Brief="Outlet Tube Stream", PosX=1, PosY=0.4915); 37 in InletShell as stream (Brief="Inlet Shell Stream", PosX=0.5237, PosY=1); 38 out OutletShell as streamPH (Brief="Outlet Shell Stream", PosX=0.5237, PosY=0); 39 40 Details as Details_Main (Brief="Details in Heat Exchanger"); 41 Tubes as Tube_Side_Main (Brief="Tube Side"); 42 Shell as Shell_Side_Main (Brief="Shell Side"); 43 Baffles as Baffles_Main (Brief="Baffles"); 44 44 45 45 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# 46 46 # Auxiliar Variables - Must be hidden 47 47 #+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++# 48 Nc 48 Nc as Real (Brief = "Number of Tube rows Crossed in one Crossflow Section",Lower=1); 49 49 Ncw as Real (Brief = "Number of Effective Crossflow rows in Each Window",Lower=1); 50 50 a as Real (Brief = "Variable for calculating Ji heat transfer correction Factor",Lower=1e-3); 51 51 b as Real (Brief = "Variable for calculating shell side pressure drop friction Factor",Lower=1e-3); 52 Rb 52 Rb as Real (Brief = "ByPass Correction Factor for Pressure Drop",Lower=1e-3); 53 53 Rss as Real (Brief = "Correction Factor for Pressure Drop",Lower=1e-3); 54 54 Rspd as Real (Brief = "Pressure Drop Correction Factor for Unequal Baffle Spacing",Lower=1e-3); 55 mw as Real (Brief = "Mass Velocity in Window Zone", Unit='kg/m^2/s');55 mw as Real (Brief = "Mass Velocity in Window Zone", Unit='kg/m^2/s'); 56 56 57 57 PARAMETERS 58 outer PP 58 outer PP as Plugin (Brief="External Physical Properties",Type = "PP"); 59 59 outer NComp as Integer (Brief="Number of Components"); 60 60 61 Pi 62 M(NComp) 63 64 TubeFlowRegime as Switcher (Brief="Tube Side Flow Regime ",Valid=["laminar","transition","turbulent"],Default="laminar");65 ShellFlowRegime 66 ShellRange 67 Side 68 LaminarCorrelation 69 TransitionCorrelation 61 Pi as constant (Brief="Pi Number",Default=3.14159265); 62 M(NComp) as molweight (Brief="Component Mol Weight"); 63 64 TubeFlowRegime as Switcher (Brief="Tube Side Flow Regime ",Valid=["laminar","transition","turbulent"],Default="laminar"); 65 ShellFlowRegime as Switcher (Brief="Shell Side Flow Regime ",Valid=["deep laminar","laminar","turbulent"],Default="deep laminar"); 66 ShellRange as Switcher (Brief="Shell Side Flow Regime Range for Correction Factor",Valid=["range1","range2","range3", "range4","range5"],Default="range1"); 67 Side as Switcher (Brief="Flag for Fluid Alocation ",Valid=["shell","tubes"],Default="shell"); 68 LaminarCorrelation as Switcher (Brief="Tube Heat Transfer Correlation in Laminar Flow",Valid=["Hausen","Schlunder"],Default="Hausen"); 69 TransitionCorrelation as Switcher (Brief="Tube Heat Transfer Correlation in Transition Flow",Valid=["Gnielinski","ESDU"],Default="Gnielinski"); 70 70 TurbulentCorrelation as Switcher (Brief="Tube Heat Transfer Correlation in Turbulent Flow",Valid=["Petukhov","SiederTate"],Default="Petukhov"); 71 71 -
trunk/eml/heat_exchangers/HeatExchangerSimplified.mso
r307 r325 29 29 30 30 PARAMETERS 31 outer PP as Plugin(Brief="External Physical Properties", Type="PP");32 outer NComp as Integer 31 outer PP as Plugin (Brief="External Physical Properties", Type="PP"); 32 outer NComp as Integer (Brief="Number of Components"); 33 33 34 34 M(NComp) as molweight (Brief="Component Mol Weight"); … … 36 36 VARIABLES 37 37 38 in InletHot as stream (Brief="Inlet Hot Stream");39 out OutletHot as streamPH (Brief="Outlet Hot Stream");40 in InletCold as stream (Brief="Inlet Cold Stream");41 out OutletCold as streamPH (Brief="Outlet Cold Stream");42 43 xh(NComp) as fraction (Brief = "Liquid Molar Fraction in Hot Side");44 yh(NComp) as fraction(Brief = "Vapour Molar Fraction in Hot Side");45 vh as fraction (Brief = "Vapour Molar Fraction in Hot Side");46 47 xc(NComp) as fraction (Brief = "Liquid Molar Fraction in Cold Side");48 yc(NComp) as fraction(Brief = "Vapour Molar Fraction in Cold Side");49 vc as fraction (Brief = "Vapour Molar Fraction in Cold Side");50 51 Details as Details_Main(Brief="Heat Exchanger Details");52 HotSide as Main_Simplified 53 ColdSide as Main_Simplified 38 in InletHot as stream (Brief="Inlet Hot Stream", PosX=0, PosY=0.4915); 39 out OutletHot as streamPH (Brief="Outlet Hot Stream", PosX=1, PosY=0.4915); 40 in InletCold as stream (Brief="Inlet Cold Stream", PosX=0.5237, PosY=1); 41 out OutletCold as streamPH (Brief="Outlet Cold Stream", PosX=0.5237, PosY=0); 42 43 xh(NComp) as fraction (Brief = "Liquid Molar Fraction in Hot Side"); 44 yh(NComp) as fraction (Brief = "Vapour Molar Fraction in Hot Side"); 45 vh as fraction (Brief = "Vapour Molar Fraction in Hot Side"); 46 47 xc(NComp) as fraction (Brief = "Liquid Molar Fraction in Cold Side"); 48 yc(NComp) as fraction (Brief = "Vapour Molar Fraction in Cold Side"); 49 vc as fraction (Brief = "Vapour Molar Fraction in Cold Side"); 50 51 Details as Details_Main (Brief="Heat Exchanger Details"); 52 HotSide as Main_Simplified (Brief="Heat Exchanger Hot Side"); 53 ColdSide as Main_Simplified (Brief="Heat Exchanger Cold Side"); 54 54 55 55 SET -
trunk/eml/heat_exchangers/Mheatex.mso
r307 r325 30 30 PARAMETERS 31 31 32 outer PP 33 outer 34 FlowDirection as Switcher(Brief="Flow Direction",Valid=["counter","cocurrent"],Default="cocurrent");35 Ncold 36 Nhot 32 outer PP as Plugin (Brief="Physical Properties", Type="PP"); 33 outer NComp as Integer (Brief="Number of Components"); 34 FlowDirection as Switcher (Brief="Flow Direction",Valid=["counter","cocurrent"],Default="cocurrent"); 35 Ncold as Integer (Brief="Number of Inlet Cold Streams",Lower=1); 36 Nhot as Integer (Brief="Number of Inlet Hot Streams",Lower=1); 37 37 38 38 VARIABLES 39 39 40 in InletHot(Nhot) as stream (Brief="Inlet Hot Streams");41 out OutletHot(Nhot) as streamPH (Brief="Outlet Hot Streams");42 in InletCold(Ncold) as stream (Brief="Inlet Cold Streams");43 out OutletCold(Ncold) as streamPH (Brief="Outlet Cold Streams");40 in InletHot(Nhot) as stream (Brief="Inlet Hot Streams", PosX=0, PosY=0.7156); 41 out OutletHot(Nhot) as streamPH (Brief="Outlet Hot Streams", PosX=1, PosY=0.7156); 42 in InletCold(Ncold) as stream (Brief="Inlet Cold Streams", PosX=1, PosY=0.2793); 43 out OutletCold(Ncold) as streamPH (Brief="Outlet Cold Streams", PosX=0, PosY=0.2793); 44 44 45 45 Method as LMTD_Basic (Brief="Log Mean Temperature Difference Method"); 46 Q as power 47 UA as Real 46 Q as power (Brief="Heat Transfer", Default=7000, Lower=1e-6, Upper=1e10); 47 UA as Real (Brief="UA product",Unit='W/K',Lower=1e-8); 48 48 49 49 EQUATIONS -
trunk/eml/heat_exchangers/PHE.mso
r307 r325 20 20 Model PHE 21 21 22 ATTRIBUTES22 ATTRIBUTES 23 23 Icon = "icon/phe"; 24 24 Pallete = true; … … 45 45 "; 46 46 47 PARAMETERS48 49 outer PP as Plugin(Brief="External Physical Properties", Type="PP");50 outer NComp as Integer 51 Pi as constant(Brief="Pi Number",Default=3.14159265);52 Kp1(15) as constant(Brief="First constant in Kumar calculation for Pressure Drop");53 Kp2(15) as constant(Brief="Second constant in Kumar calculation for Pressure Drop");54 Kc1(14) as constant(Brief="First constant in Kumar calculation for Heat Transfer");55 Kc2(14) as constant(Brief="Second constant Kumar calculation for Heat Transfer");47 PARAMETERS 48 49 outer PP as Plugin (Brief="External Physical Properties", Type="PP"); 50 outer NComp as Integer (Brief="Number of Chemical Components"); 51 Pi as constant (Brief="Pi Number",Default=3.14159265); 52 Kp1(15) as constant (Brief="First constant in Kumar calculation for Pressure Drop"); 53 Kp2(15) as constant (Brief="Second constant in Kumar calculation for Pressure Drop"); 54 Kc1(14) as constant (Brief="First constant in Kumar calculation for Heat Transfer"); 55 Kc2(14) as constant (Brief="Second constant Kumar calculation for Heat Transfer"); 56 56 M(NComp) as molweight (Brief="Component Mol Weight"); 57 57 58 ChevronAngle as Switcher 59 Method as Switcher(Brief="Method of Thermal Calculation",Valid=["NTU","LMTD"],Default="NTU");60 SideOne as Switcher(Brief="Fluid Alocation in the Side I - (The odd channels)",Valid=["hot","cold"],Default="hot");61 Nchannels as Integer(Brief="Total Number of Channels in The Whole Heat Exchanger");62 Nplates as Integer(Brief="Total Number of Plates in The Whole Heat Exchanger",Default=25);63 NpassHot as Integer(Brief="Number of Passes for Hot Side");64 NpassCold as Integer(Brief="Number of Passes for Cold Side");65 Dports as length(Brief="Ports Diameter",Lower=1e-6);66 Atotal as area(Brief="Total Effective Area",Lower=1e-6);67 Aports as area(Brief="Port Opening Area of Plate",Lower=1e-6);68 Achannel as area(Brief="Cross-Sectional Area for Channel Flow",Lower=1e-6);69 Dh as length(Brief="Equivalent Diameter of Channel",Lower=1e-6);70 Depth as length(Brief="Corrugation Depth",Lower=1e-6);71 PhiFactor as Real(Brief="Enlargement Factor",Lower=1e-6);72 Lp as length(Brief="Plate Vertical Distance between Port Centers",Lower=0.1);73 Lpack as length(Brief="Compact Plate Pack Length",Lower=0.1);74 Lv as length(Brief="Vertical Ports Distance",Lower=0.1);75 Lh as length(Brief="Plate Horizontal Distance between Port Centers",Lower=0.1);76 Lw as length(Brief="Plate Width",Lower=0.1);77 pitch as length(Brief="Plate Pitch",Lower=0.1);78 pt as length(Brief="Plate Thickness",Lower=0.1);79 Kwall as conductivity(Brief="Plate Thermal Conductivity",Default=1.0);80 Rfh as positive(Brief="Hot Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);81 Rfc as positive(Brief="Cold Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0);82 83 VARIABLES84 85 in InletHot as stream (Brief="Inlet Hot Stream");86 in InletCold as stream (Brief="Inlet Cold Stream");87 out OutletHot as streamPH (Brief="Outlet Hot Stream");88 out OutletCold as streamPH (Brief="Outlet Cold Stream");89 90 HotSide as Main_PHE(Brief="Plate Heat Exchanger Hot Side");91 ColdSide as Main_PHE(Brief="Plate Heat Exchanger Cold Side");92 Thermal as Thermal_PHE(Brief="Thermal Results");93 94 SET58 ChevronAngle as Switcher (Brief="Chevron Corrugation Inclination Angle in Degrees ",Valid=["30","45","50","60","65"],Default="30"); 59 Method as Switcher (Brief="Method of Thermal Calculation",Valid=["NTU","LMTD"],Default="NTU"); 60 SideOne as Switcher (Brief="Fluid Alocation in the Side I - (The odd channels)",Valid=["hot","cold"],Default="hot"); 61 Nchannels as Integer (Brief="Total Number of Channels in The Whole Heat Exchanger"); 62 Nplates as Integer (Brief="Total Number of Plates in The Whole Heat Exchanger",Default=25); 63 NpassHot as Integer (Brief="Number of Passes for Hot Side"); 64 NpassCold as Integer (Brief="Number of Passes for Cold Side"); 65 Dports as length (Brief="Ports Diameter",Lower=1e-6); 66 Atotal as area (Brief="Total Effective Area",Lower=1e-6); 67 Aports as area (Brief="Port Opening Area of Plate",Lower=1e-6); 68 Achannel as area (Brief="Cross-Sectional Area for Channel Flow",Lower=1e-6); 69 Dh as length (Brief="Equivalent Diameter of Channel",Lower=1e-6); 70 Depth as length (Brief="Corrugation Depth",Lower=1e-6); 71 PhiFactor as Real (Brief="Enlargement Factor",Lower=1e-6); 72 Lp as length (Brief="Plate Vertical Distance between Port Centers",Lower=0.1); 73 Lpack as length (Brief="Compact Plate Pack Length",Lower=0.1); 74 Lv as length (Brief="Vertical Ports Distance",Lower=0.1); 75 Lh as length (Brief="Plate Horizontal Distance between Port Centers",Lower=0.1); 76 Lw as length (Brief="Plate Width",Lower=0.1); 77 pitch as length (Brief="Plate Pitch",Lower=0.1); 78 pt as length (Brief="Plate Thickness",Lower=0.1); 79 Kwall as conductivity (Brief="Plate Thermal Conductivity",Default=1.0); 80 Rfh as positive (Brief="Hot Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 81 Rfc as positive (Brief="Cold Side Fouling Resistance",Unit='m^2*K/kW',Default=1e-6,Lower=0); 82 83 VARIABLES 84 85 in InletHot as stream (Brief="Inlet Hot Stream", PosX=0, PosY=0.7156); 86 in InletCold as stream (Brief="Inlet Cold Stream", PosX=1, PosY=0.7156); 87 out OutletHot as streamPH (Brief="Outlet Hot Stream", PosX=0, PosY=0.2793); 88 out OutletCold as streamPH (Brief="Outlet Cold Stream", PosX=1, PosY=0.2793); 89 90 HotSide as Main_PHE (Brief="Plate Heat Exchanger Hot Side"); 91 ColdSide as Main_PHE (Brief="Plate Heat Exchanger Cold Side"); 92 Thermal as Thermal_PHE (Brief="Thermal Results"); 93 94 SET 95 95 96 96 #"First constant for Kumar's calculating Pressure Drop" … … 139 139 Atotal =(Nplates-2)*Lw*Lp*PhiFactor; 140 140 141 EQUATIONS141 EQUATIONS 142 142 143 143 "Hot Stream Average Temperature" … … 165 165 ColdSide.Properties.Average.Mw = sum(M*InletCold.z); 166 166 167 if InletCold.v equal 0167 if InletCold.v equal 0 168 168 169 169 then … … 249 249 ColdSide.Properties.Wall.Mu = PP.VapourViscosity(ColdSide.Properties.Wall.Twall,ColdSide.Properties.Average.P,InletCold.z); 250 250 251 end252 253 if InletHot.v equal 0251 end 252 253 if InletHot.v equal 0 254 254 255 255 then … … 336 336 HotSide.Properties.Wall.Mu = PP.VapourViscosity(HotSide.Properties.Wall.Twall,HotSide.Properties.Average.P,InletHot.z); 337 337 338 end338 end 339 339 340 340 "Energy Balance Hot Stream" … … 368 368 OutletCold.z=InletCold.z; 369 369 370 switch SideOne370 switch SideOne 371 371 372 372 case "cold": … … 386 386 ColdSide.PressureDrop.Npassage = (2*Nchannels-1+(-1)^(Nchannels))/(4*NpassCold); 387 387 388 end388 end 389 389 390 390 "Hot Stream Mass Flux in the Channel" … … 418 418 ColdSide.PressureDrop.Pdrop =ColdSide.PressureDrop.DPchannel+ColdSide.PressureDrop.DPports; 419 419 420 switch ChevronAngle #Pressure Drop Friction Factor According to kumar's (1984)420 switch ChevronAngle #Pressure Drop Friction Factor According to kumar's (1984) 421 421 422 422 case "30": # ChevronAngle <= 30 … … 510 510 end 511 511 512 end513 514 switch ChevronAngle # Heat Transfer Coefficient According to kumar's (1984)512 end 513 514 switch ChevronAngle # Heat Transfer Coefficient According to kumar's (1984) 515 515 516 516 case "30": # ChevronAngle <= 30 … … 593 593 end 594 594 595 end595 end 596 596 597 597 "Hot Stream Velocity in Channels" … … 637 637 Thermal.Ud*(1/HotSide.HeatTransfer.hcoeff +pt/Kwall+1/ColdSide.HeatTransfer.hcoeff + Rfc + Rfh)=1; 638 638 639 switch Method639 switch Method 640 640 641 641 case "LMTD": … … 649 649 Thermal.Q = Thermal.Eft*Thermal.Cmin*(InletHot.T-InletCold.T); 650 650 651 end651 end 652 652 653 653 "Temperature Difference at Inlet - Counter Flow" … … 681 681 ColdSide.HeatTransfer.NTU*ColdSide.HeatTransfer.WCp = Thermal.Ud*Atotal; 682 682 683 if Thermal.Cr equal 1 # To be Fixed: Effectiveness in true counter flow !683 if Thermal.Cr equal 1 # To be Fixed: Effectiveness in true counter flow ! 684 684 685 685 then … … 697 697 Thermal.Fc =(ln(abs(1-Thermal.Eft*Thermal.Cr))-ln(abs(1-Thermal.Eft)))/(Thermal.NTU*(1-Thermal.Cr)); 698 698 699 end 699 700 end 700 701 702 end -
trunk/eml/heat_exchangers/heater.mso
r319 r325 17 17 *--------------------------------------------------------------------*# 18 18 19 using "streams .mso";19 using "streams"; 20 20 21 Model heater 21 Model heater_basic 22 22 23 ATTRIBUTES 24 Pallete = true; 25 Icon = "icon/heater"; 26 Brief = "Heater or Cooler"; 23 ATTRIBUTES 24 Pallete = false; 25 Brief = "Basic model for Heater or Cooler"; 27 26 Info = 28 27 "Determines thermal and phase conditions of an outlet stream"; 28 29 PARAMETERS 30 outer PP as Plugin (Brief="Physical Properties", Type="PP"); 31 outer NComp as Integer (Brief="Number of Components"); 32 Ninlet as Integer (Brief="Number of Inlet Streams"); 33 Kvalues as Switcher (Brief="Option for Display Phase Equilibrium K-values",Valid=["yes","no"], Default="yes"); 34 35 VARIABLES 36 QDuty as power (Brief = "Actual Duty"); 37 Vfrac as fraction (Brief = "Vapor fraction Outlet Stream"); 38 Lfrac as fraction (Brief = "Liquid fraction Outlet Stream"); 39 Kvalue(NComp) as Real (Brief = "Phase Equilibrium K-values"); 40 in Inlet(Ninlet ) as stream (Brief="Inlet Streams", PosX=0, PosY=0.4833); 41 out Outlet as streamPH (Brief="Outlet Stream", PosX=1, PosY=0.4782); 29 42 30 PARAMETERS 31 outer PP as Plugin (Brief="Physical Properties", Type="PP"); 32 outer NComp as Integer (Brief="Number of Components"); 33 Ninlet as Integer (Brief="Number of Inlet Streams"); 34 Kvalues as Switcher (Brief="Option for Display Phase Equilibrium K-values",Valid=["yes","no"], Default="yes"); 35 36 SUBMODELS 37 in Inlet(Ninlet ) as stream (Brief="Inlet Streams"); 38 out Outlet as streamPH (Brief="Outlet Stream"); 39 in InletHeat as energy_stream (Brief ="Inlet Heat Stream"); 40 out OutletHeat as energy_stream (Brief ="Outlet Heat Stream"); 43 EQUATIONS 41 44 42 VARIABLES 43 QDuty as power (Brief = "Actual Duty"); 44 Vfrac as fraction (Brief = "Vapor fraction Outlet Stream"); 45 Lfrac as fraction (Brief = "Liquid fraction Outlet Stream"); 46 Kvalue(NComp) as Real (Brief = "Phase Equilibrium K-values"); 47 48 EQUATIONS 49 50 "Flow" 45 "Flow" 51 46 Outlet.F = sum(Inlet.F); 52 47 53 for j in [1 : NComp]48 for j in [1 : NComp] 54 49 55 "Composition"50 "Composition" 56 51 Outlet.F*Outlet.z(j) = sum(Inlet.F*Inlet.z(j)); 57 52 58 end53 end 59 54 60 "Vapor fraction Outlet Stream"55 "Vapor fraction Outlet Stream" 61 56 Vfrac = Outlet.v; 62 57 63 "Liquid fraction Outlet Stream"58 "Liquid fraction Outlet Stream" 64 59 Lfrac = 1-Vfrac; 65 60 66 "Heat Duty"61 "Heat Duty" 67 62 QDuty = Outlet.F*Outlet.h - sum(Inlet.F*Inlet.h); 68 63 69 "Net Duty" 70 OutletHeat.Q = InletHeat.Q - QDuty; 71 72 switch Kvalues # Fix for better convergence !!! 64 switch Kvalues # Fix for better convergence !!! 73 65 case "yes": 74 "K-values Phase Equilibrium"66 "K-values Phase Equilibrium" 75 67 Kvalue*(Outlet.x) = Outlet.y; 76 68 77 69 case "no": 78 "K-values Phase Equilibrium"70 "K-values Phase Equilibrium" 79 71 Kvalue = 1; 80 72 81 end 82 73 end 83 74 end 84 75 76 77 Model heater as heater_basic 78 79 ATTRIBUTES 80 Pallete = true; 81 Icon = "icon/heater"; 82 Brief = "Heater"; 83 Info = 84 "Determines thermal and phase conditions of an outlet stream"; 85 86 VARIABLES 87 in InletQ as energy_stream (Brief ="Inlet Heat Stream", PosX=0.1662, PosY=1); 88 out OutletQ as energy_stream (Brief ="Outlet Heat Stream", PosX=0.8245, PosY=0); 89 90 EQUATIONS 91 "Net Duty" 92 OutletQ.Q = InletQ.Q - QDuty; 93 end 94 95 Model cooler as heater_basic 96 97 ATTRIBUTES 98 Pallete = true; 99 Icon = "icon/cooler"; 100 Brief = "Cooler"; 101 Info = 102 "Determines thermal and phase conditions of an outlet stream"; 103 104 VARIABLES 105 in InletQ as energy_stream (Brief ="Inlet Heat Stream", PosX=0.8245, PosY=1); 106 out OutletQ as energy_stream (Brief ="Outlet Heat Stream", PosX=0.1662, PosY=0); 107 108 EQUATIONS 109 "Net Duty" 110 OutletQ.Q = InletQ.Q - QDuty; 111 end -
trunk/eml/heat_exchangers/icon/cooler.svg
r96 r325 12 12 id="svg10172" 13 13 sodipodi:version="0.32" 14 inkscape:version="0.4 4.1"14 inkscape:version="0.45.1" 15 15 width="210mm" 16 16 height="297mm" 17 sodipodi:docbase="C:\Documents and Settings\Administrador\Meus documentos\Bolsa Simulação\Desenhos2" 18 sodipodi:docname="cooler.svg"> 17 sodipodi:docbase="D:\User\Arge\PROJETOS\Alsoc\EMSO\mso\eml\heat_exchangers\icon" 18 sodipodi:docname="cooler.svg" 19 inkscape:output_extension="org.inkscape.output.svg.inkscape" 20 inkscape:export-filename="D:\User\Arge\PROJETOS\Alsoc\EMSO\mso\eml\heat_exchangers\icon\cooler.png" 21 inkscape:export-xdpi="90" 22 inkscape:export-ydpi="90"> 19 23 <metadata 20 24 id="metadata10177"> … … 65 69 </defs> 66 70 <sodipodi:namedview 67 inkscape:window-height=" 900"68 inkscape:window-width="1 280"71 inkscape:window-height="753" 72 inkscape:window-width="1024" 69 73 inkscape:pageshadow="2" 70 74 inkscape:pageopacity="0.0" … … 78 82 showgrid="true" 79 83 inkscape:zoom="1.4329667" 80 inkscape:cx="3 49.04822"81 inkscape:cy="6 50.17997"82 inkscape:window-x=" -4"83 inkscape:window-y=" -4"84 inkscape:cx="311.3198" 85 inkscape:cy="609.68654" 86 inkscape:window-x="0" 87 inkscape:window-y="0" 84 88 inkscape:current-layer="svg10172" /> 85 89 <path 86 style="fill:none;fill-rule:evenodd;stroke: black;stroke-width:3.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1"87 d="M 4 00,524.36218 C 400,574.36218 400,574.36218 400,574.36218"90 style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:3.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1" 91 d="M 490.77562,434.77475 C 540.77196,435.38039 540.77196,435.38039 540.77196,435.38039" 88 92 id="path13752" /> 89 93 <path 90 94 sodipodi:type="arc" 91 style="fill:url(#linearGradient14653);fill-opacity:1;fill-rule:nonzero;stroke: black;stroke-width:4.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1"95 style="fill:url(#linearGradient14653);fill-opacity:1;fill-rule:nonzero;stroke:#000000;stroke-width:4.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1" 92 96 id="path10179" 93 97 sodipodi:cx="400.19324" … … 96 100 sodipodi:ry="100.17168" 97 101 d="M 500.36491 352.14737 A 100.17168 100.17168 0 1 1 300.02156,352.14737 A 100.17168 100.17168 0 1 1 500.36491 352.14737 z" 98 transform=" translate(0,74)" />102 transform="matrix(-1.2112799e-2,0.9999266,0.9999266,1.2112799e-2,45.291601,29.348951)" /> 99 103 <path 100 style="fill:none;fill-rule:evenodd;stroke: black;stroke-width:3.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1"101 d="M 288,542.36218 L 378,392.36218 L 408,462.36218 L 498,312.36218"104 style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:3.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1" 105 d="M 510.13094,323.001 L 359.05179,411.17748 L 428.68327,442.02317 L 277.60412,530.19965" 102 106 id="path11973" /> 103 107 <path 104 style="fill: black;fill-opacity:1;fill-rule:evenodd;stroke:black;stroke-width:1.29729998px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"105 d="M 490.64865,314.82308 L 501.57869,321.71353 L 506.79992,301.90858 L 490.64865,314.82308z "108 style="fill:#000000;fill-opacity:1;fill-rule:evenodd;stroke:#000000;stroke-width:1.29729998px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1" 109 d="M 280.15389,522.87865 L 286.91144,533.89135 L 267.0447,538.8723 L 280.15389,522.87865 z " 106 110 id="path10086" 107 inkscape:transform-center-x="- 8.486542"108 inkscape:transform-center-y=" -50.636822" />111 inkscape:transform-center-x="-7.806366" 112 inkscape:transform-center-y="50.694172" /> 109 113 <path 110 style="fill:none;fill-rule:evenodd;stroke: black;stroke-width:3.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1"111 d="M 402,274.36218 C 402,324.36218 402,324.36218 402,324.36218 L 402,324.36218 L 402,324.36218 L 402,274.36218z "114 style="fill:none;fill-rule:evenodd;stroke:#000000;stroke-width:3.5;stroke-linecap:butt;stroke-linejoin:miter;stroke-miterlimit:4;stroke-dasharray:none;stroke-opacity:1" 115 d="M 232.76974,433.74641 C 282.76607,434.35205 282.76607,434.35205 282.76607,434.35205 L 282.76607,434.35205 L 282.76607,434.35205 L 232.76974,433.74641 z " 112 116 id="path12861" /> 113 117 <path 114 style="fill: black;fill-opacity:1;fill-rule:evenodd;stroke:black;stroke-width:1.00203896px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"115 d="M 406.48459,306.49581 L 396.50102,307.22055 L 402.63462,321.86116 L 406.48459,306.49581z "118 style="fill:#000000;fill-opacity:1;fill-rule:evenodd;stroke:#000000;stroke-width:1.00203896px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1" 119 d="M 274.50102,438.86116 L 274.56591,428.85153 L 289.57909,434.00644 L 274.50102,438.86116 z " 116 120 id="path14641" 117 inkscape:transform-center-x="-1 7.236487"118 inkscape:transform-center-y=" 37.345131" />121 inkscape:transform-center-x="-14.938788" 122 inkscape:transform-center-y="-37.917701" /> 119 123 <path 120 style="fill: black;fill-opacity:1;fill-rule:evenodd;stroke:black;stroke-width:1.00203896px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1"121 d="M 404.51069,572.7855 L 394.50102,572.84548 L 399.64858,587.86116 L 404.51069,572.7855 z "124 style="fill:#000000;fill-opacity:1;fill-rule:evenodd;stroke:#000000;stroke-width:1.00203896px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1" 125 d="M 539.14075,439.87165 L 539.32198,429.86344 L 554.2742,435.19251 L 539.14075,439.87165 z " 122 126 id="path14643" 123 inkscape:transform-center-x="-15. 204017"124 inkscape:transform-center-y=" 38.611301" />127 inkscape:transform-center-x="-15.639411" 128 inkscape:transform-center-y="-38.515247" /> 125 129 </svg> -
trunk/eml/heat_exchangers/icon/heater.svg
r316 r325 15 15 width="210mm" 16 16 height="297mm" 17 sodipodi:docbase=" E:\Documentos Pessoais\Gerson\repositorio\trunk\eml\heat_exchangers\icon"17 sodipodi:docbase="D:\User\Arge\PROJETOS\Alsoc\EMSO\mso\eml\heat_exchangers\icon" 18 18 sodipodi:docname="heater.svg" 19 inkscape:output_extension="org.inkscape.output.svg.inkscape"> 19 inkscape:output_extension="org.inkscape.output.svg.inkscape" 20 inkscape:export-filename="D:\User\Arge\PROJETOS\Alsoc\EMSO\mso\eml\heat_exchangers\icon\heater.png" 21 inkscape:export-xdpi="90" 22 inkscape:export-ydpi="90"> 20 23 <metadata 21 24 id="metadata10177"> … … 78 81 id="base" 79 82 showgrid="true" 80 inkscape:zoom=" 0.71648337"81 inkscape:cx=" 285.99484"82 inkscape:cy=" 609.30613"83 inkscape:zoom="1.0132605" 84 inkscape:cx="439.12625" 85 inkscape:cy="717.22114" 83 86 inkscape:window-x="0" 84 87 inkscape:window-y="22" … … 119 122 inkscape:transform-center-x="-38.963517" 120 123 inkscape:transform-center-y="-14.381614" /> 124 <path 125 style="fill:#000000;fill-opacity:1;fill-rule:evenodd;stroke:#000000;stroke-width:1.00203896px;stroke-linecap:butt;stroke-linejoin:miter;stroke-opacity:1" 126 d="M 504.79979,310.76656 L 511.24729,318.42336 L 519.49898,304.8632 L 504.79979,310.76656 z " 127 id="path2168" 128 inkscape:transform-center-x="-18.314027" 129 inkscape:transform-center-y="-34.531454" /> 121 130 </svg> -
trunk/eml/mixers_splitters/mixer.mso
r313 r325 38 38 39 39 VARIABLES 40 in Inlet_mixer(Ninlet) as stream (Brief = "Inlet streams" );41 out Outlet as stream (Brief = "Outlet stream" );40 in Inlet_mixer(Ninlet) as stream (Brief = "Inlet streams", PosX=0.5, PosY=0); 41 out Outlet as stream (Brief = "Outlet stream", PosX=1, PosY=0.5059); 42 42 out Out_int as stream (Brief = "Intermediate Outlet stream"); 43 43 zeroQ as energy_source (Brief="No Heat rate supplied"); -
trunk/eml/mixers_splitters/splitter.mso
r304 r325 43 43 44 44 VARIABLES 45 in Inlet as stream ;46 out Outlet(NOutlet) as stream ;45 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.5001); 46 out Outlet(NOutlet) as stream (Brief = "Outlet streams", PosX=1, PosY=0.5); 47 47 frac(NOutlet) as fraction (Brief = "Distribution of Outlets", Default=0.5); 48 48 … … 94 94 95 95 VARIABLES 96 in Inlet as stream ;97 out Outlet1 as stream ;98 out Outlet2 as stream ;96 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.5069); 97 out Outlet1 as stream (Brief = "Outlet stream 1", PosX=1, PosY=0.3027); 98 out Outlet2 as stream (Brief = "Outlet stream 2", PosX=1, PosY=0.7141); 99 99 frac as fraction (Brief = "Fraction to Outlet 1"); 100 100 -
trunk/eml/pressure_changers/compressor.mso
r305 r325 18 18 *--------------------------------------------------------------------*# 19 19 20 using " pressure_changers/flux_machine_basic";20 using "streams"; 21 21 22 Model centrifugal_compressor as flux_machine_basic_PH22 Model centrifugal_compressor 23 23 ATTRIBUTES 24 24 Pallete = true; … … 53 53 Ws as energy_mol (Brief = "Isentropic Head"); 54 54 Tiso as temperature (Brief = "Isentropic Temperature"); 55 Effp as efficiency(Brief = "Politropic efficiency");55 Effp as positive (Brief = "Politropic efficiency"); 56 56 Effs as efficiency (Brief = "Isentropic efficiency"); 57 57 FPower as power (Brief = "Fluid Power"); 58 58 Mwm as molweight (Brief = "Mixture Molar Weight"); 59 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.5086); 60 out Outlet as streamPH (Brief = "Outlet stream", PosX=1, PosY=0.5022); 59 61 60 62 SET … … 86 88 "Calculate Isentropic Outlet Temperature" 87 89 Tiso = Inlet.T * (Outlet.P/Inlet.P)^((k-1)/k); 88 90 # PP.VapourEntropy(Tiso, Outlet.P, Outlet.z) = 91 # PP.VapourEntropy(Inlet.T, Inlet.P, Inlet.z); 92 89 93 "Calculate Real Outlet Temperature" 90 94 Effs * (Outlet.T- Inlet.T) = (Tiso - Inlet.T); -
trunk/eml/pressure_changers/pump.mso
r305 r325 18 18 *--------------------------------------------------------------------*# 19 19 20 using " pressure_changers/flux_machine_basic";20 using "streams"; 21 21 22 Model centrifugal_pump as flux_machine_basic22 Model centrifugal_pump 23 23 ATTRIBUTES 24 24 Pallete = true; … … 61 61 pvm as pressure (Brief = "Mixture Vapour Pressure"); 62 62 NPSHa as length (Brief = "Available Net Positive Suction Head"); 63 NS as positive (Brief = "Specific Speed" , Unit = '(rpm*(gal/min)^0.5)/(m^3/4)');63 NS as positive (Brief = "Specific Speed"); 64 64 Q as flow_vol (Brief = "Volumetric Flow Rate"); 65 65 vm as vol_mol (Brief = "Mixture Molar Volume"); 66 66 in Inlet as stream (Brief = "Inlet stream", PosX=1, PosY=0.4025); 67 out Outlet as stream (Brief = "Outlet stream", PosX=0, PosY=0); 68 67 69 SET 68 70 Mw = PP.MolecularWeight(); … … 127 129 "Calculate Volumetric Flow Rate" 128 130 Q = Inlet.F*vm; 129 131 130 132 "Calculate Specific Speed" 131 NS = N*(Q^0.5)/( Head^3/4);132 133 NS = N*(Q^0.5)/((Head/Mwm)^(3/4)); 134 133 135 end 134 135 136 136 137 #*------------------------------------------------------------------- … … 156 157 157 158 VARIABLES 158 in Inlet as stream;159 out Outlet as streamPH;159 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.4727); 160 out Outlet as streamPH (Brief = "Outlet stream", PosX=1, PosY=0.1859); 160 161 161 162 dP as press_delta (Brief="Pump head"); -
trunk/eml/pressure_changers/turbine.mso
r305 r325 18 18 *--------------------------------------------------------------------*# 19 19 20 using " pressure_changers/flux_machine_basic";20 using "streams"; 21 21 22 Model Hidraulic_Turbine as flux_machine_basic22 Model Hidraulic_Turbine 23 23 ATTRIBUTES 24 24 Pallete = true; … … 58 58 rho as dens_mass (Brief = "Specific Mass"); 59 59 Cp as cp_mol (Brief = "Heat Capacity"); 60 60 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.5086); 61 out Outlet as stream (Brief = "Outlet stream", PosX=1, PosY=0.5022); 62 61 63 SET 62 64 Mw = PP.MolecularWeight(); … … 106 108 107 109 "Calculate Head" 108 Head * rho = Pdiff;110 Head = Outlet.h - Inlet.h; 109 111 end -
trunk/eml/pressure_changers/valve.mso
r305 r325 18 18 *-------------------------------------------------------------------*# 19 19 20 using " pressure_changers/flux_machine_basic";20 using "streams"; 21 21 22 22 23 Model valve as flux_machine_basic_PH23 Model valve 24 24 ATTRIBUTES 25 25 Pallete = true; … … 62 62 vm as vol_mol (Brief = "Mixture Molar Volume"); 63 63 x as fraction (Brief = "Opening"); 64 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.7365); 65 out Outlet as streamPH (Brief = "Outlet stream", PosX=1, PosY=0.7365); 64 66 65 67 SET … … 162 164 163 165 VARIABLES 164 in Inlet as stream;165 out Outlet as streamPH;166 in Inlet as stream (Brief = "Inlet stream", PosX=0, PosY=0.7365); 167 out Outlet as streamPH (Brief = "Outlet stream", PosX=1, PosY=0.7365); 166 168 x as fraction (Brief="Plug Position"); 167 169 rho as dens_mass (Brief="Fluid Density", Default=1e3); -
trunk/eml/reactors/batch.mso
r313 r325 38 38 VARIABLES 39 39 40 in Inlet as stream ;41 C(NComp) as conc_mol (Brief="Components concentration");42 r(NReac) as reaction_mol (Brief = "Reaction rates");43 T as temperature (Brief="Reactor temperature");44 P as pressure (Brief="Reactor pressure");40 in Inlet as stream (Brief="Inlet stream", PosX=0, PosY=0); 41 C(NComp) as conc_mol (Brief="Components concentration"); 42 r(NReac) as reaction_mol (Brief = "Reaction rates"); 43 T as temperature (Brief="Reactor temperature"); 44 P as pressure (Brief="Reactor pressure"); 45 45 Vr as volume (Brief = "Reacting Volume"); 46 46 -
trunk/eml/reactors/cstr.mso
r302 r325 37 37 VARIABLES 38 38 39 in Inlet as stream (Brief="Inlet Stream");40 out Outlet as streamPH (Brief="Outlet Stream" );39 in Inlet as stream (Brief="Inlet Stream", PosX=0, PosY=0); 40 out Outlet as streamPH (Brief="Outlet Stream", PosX=1, PosY=1); 41 41 42 42 q as heat_rate (Brief="Heat"); -
trunk/eml/reactors/pfr.mso
r302 r325 50 50 VARIABLES 51 51 52 in Inlet as stream (Brief = "Inlet Stream");53 out Outlet as stream (Brief = "Outlet Stream" );52 in Inlet as stream (Brief = "Inlet Stream", PosX=0, PosY=0.5076); 53 out Outlet as stream (Brief = "Outlet Stream", PosX=1, PosY=0.5236); 54 54 55 55 str(NDisc+1) as streamPH; -
trunk/eml/stage_separators/batch_dist.mso
r318 r325 48 48 49 49 VARIABLES 50 in Inlet as stream (Brief="Feed stream" );51 in InletL as stream (Brief="Liquid inlet stream" );52 out OutletV as vapour_stream (Brief="Vapour outlet stream" );53 in InletQ as energy_stream (Brief="Heat supplied" );50 in Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.9385); 51 in InletL as stream (Brief="Liquid inlet stream", PosX=0.5, PosY=0.1984); # FIXME 52 out OutletV as vapour_stream (Brief="Vapour outlet stream", PosX=1, PosY=0.1984); 53 in InletQ as energy_stream (Brief="Heat supplied", PosX=1, PosY=0.9578); 54 54 55 55 M(NComp) as mol (Brief="Molar Holdup in the distillator"); -
trunk/eml/stage_separators/condenser.mso
r310 r325 48 48 49 49 VARIABLES 50 in InletV as stream(Brief="Vapour inlet stream" );51 out OutletL as liquid_stream(Brief="Liquid outlet stream" );52 out OutletV as vapour_stream(Brief="Vapour outlet stream" );53 in InletQ as energy_stream (Brief=" Heat supplied");50 in InletV as stream(Brief="Vapour inlet stream", PosX=0.1164, PosY=0); 51 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4513, PosY=1); 52 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4723, PosY=0); 53 in InletQ as energy_stream (Brief="Cold supplied", PosX=1, PosY=0.6311); 54 54 55 55 M(NComp) as mol (Brief="Molar Holdup in the tray"); … … 127 127 128 128 VARIABLES 129 in InletV as stream(Brief="Vapour inlet stream" );130 out OutletL as liquid_stream(Brief="Liquid outlet stream" );131 in InletQ as energy_stream (Brief=" Heat supplied");129 in InletV as stream(Brief="Vapour inlet stream", PosX=0.3431, PosY=0); 130 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.34375, PosY=1); 131 in InletQ as energy_stream (Brief="Cold supplied", PosX=1, PosY=0.5974); 132 132 DP as press_delta (Brief="Pressure Drop in the condenser"); 133 133 … … 181 181 182 182 VARIABLES 183 in InletV as stream(Brief="Vapour inlet stream" );184 out OutletL as liquid_stream(Brief="Liquid outlet stream" );185 out OutletV as vapour_stream(Brief="Vapour outlet stream" );186 in InletQ as energy_stream (Brief=" Heat supplied");183 in InletV as stream(Brief="Vapour inlet stream", PosX=0.1164, PosY=0); 184 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4513, PosY=1); 185 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4723, PosY=0); 186 in InletQ as energy_stream (Brief="Cold supplied", PosX=1, PosY=0.6311); 187 187 188 188 M(NComp) as mol (Brief="Molar Holdup in the tray"); -
trunk/eml/stage_separators/flash.mso
r323 r325 50 50 51 51 VARIABLES 52 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5 );53 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0. 5, PosY=0.5);54 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0. 5, PosY=0);55 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.75 );52 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5421); 53 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4790, PosY=1); 54 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4877, PosY=0); 55 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7559); 56 56 57 57 M(NComp) as mol (Brief="Molar Holdup in the tray"); … … 140 140 141 141 VARIABLES 142 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5 );143 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0. 5, PosY=0.5);144 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0. 5, PosY=0);145 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.75 );142 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5421); 143 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4790, PosY=1); 144 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4877, PosY=0); 145 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7559); 146 146 vfrac as fraction; 147 147 … … 191 191 192 192 VARIABLES 193 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5 );194 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0. 5, PosY=0.5);195 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0. 5, PosY=0);196 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.75 );193 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5421); 194 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4790, PosY=1); 195 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4877, PosY=0); 196 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7559); 197 197 vfrac as fraction(Brief="Real vaporization fraction"); 198 198 h as enth_mol(Brief="Mixture enthalpy"); … … 249 249 250 250 VARIABLES 251 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5 );252 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0. 5, PosY=0.5);253 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0. 5, PosY=0);254 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.75 );251 in Inlet as stream(Brief="Feed Stream", PosX=0, PosY=0.5421); 252 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.4790, PosY=1); 253 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.4877, PosY=0); 254 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7559); 255 255 vfrac as fraction(Brief="Real vaporization fraction"); 256 256 vsat as Real(Lower=-0.1, Upper=1.1, Brief="Vaporization fraction if saturated"); -
trunk/eml/stage_separators/reboiler.mso
r310 r325 56 56 57 57 VARIABLES 58 in Inlet as stream(Brief="Feed Stream" );59 in InletL as stream(Brief="Liquid inlet stream" );60 out OutletL as liquid_stream(Brief="Liquid outlet stream" );61 out OutletV as vapour_stream(Brief="Vapour outlet stream" );62 in InletQ as energy_stream (Brief="Heat supplied" );58 in Inlet as stream(Brief="Feed Stream", PosX=0.8127, PosY=0); 59 in InletL as stream(Brief="Liquid inlet stream", PosX=0, PosY=0.5254); 60 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.2413, PosY=1); 61 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.5079, PosY=0); 62 in InletQ as energy_stream (Brief="Heat supplied", PosX=1, PosY=0.6123); 63 63 64 64 M(NComp) as mol (Brief="Molar Holdup in the tray"); … … 141 141 142 142 VARIABLES 143 in InletL as stream(Brief="Liquid inlet stream" );144 out OutletV as vapour_stream(Brief="Vapour outlet stream" );145 in InletQ as energy_stream (Brief="Heat supplied" );143 in InletL as stream(Brief="Liquid inlet stream", PosX=0.3345, PosY=1); 144 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.3369, PosY=0); 145 in InletQ as energy_stream (Brief="Heat supplied", PosX=1, PosY=0.6111); 146 146 vV as volume_mol (Brief="Vapour Molar volume"); 147 147 rhoV as dens_mass (Brief="Vapour Density"); … … 188 188 189 189 VARIABLES 190 in InletL as stream(Brief="Liquid inlet stream" );191 out OutletV as stream(Brief="Vapour outlet stream");192 in InletQ as energy_stream (Brief="Heat supplied" );190 in InletL as stream(Brief="Liquid inlet stream", PosX=0.3345, PosY=1); 191 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.3369, PosY=0); 192 in InletQ as energy_stream (Brief="Heat supplied", PosX=1, PosY=0.6111); 193 193 194 194 EQUATIONS … … 252 252 253 253 VARIABLES 254 in Inlet as stream(Brief="Feed Stream" );255 in InletL as stream(Brief="Liquid inlet stream" );256 out OutletL as liquid_stream(Brief="Liquid outlet stream" );257 out OutletV as vapour_stream(Brief="Vapour outlet stream" );258 259 in InletQ as energy_stream (Brief="Heat supplied"); 254 in Inlet as stream(Brief="Feed Stream", PosX=0.8127, PosY=0); 255 in InletL as stream(Brief="Liquid inlet stream", PosX=0, PosY=0.5254); 256 out OutletL as liquid_stream(Brief="Liquid outlet stream", PosX=0.2413, PosY=1); 257 out OutletV as vapour_stream(Brief="Vapour outlet stream", PosX=0.5079, PosY=0); 258 in InletQ as energy_stream (Brief="Heat supplied", PosX=1, PosY=0.6123); 259 260 260 M(NComp) as mol (Brief="Molar Holdup in the tray"); 261 261 ML as mol (Brief="Molar liquid holdup"); -
trunk/eml/stage_separators/tank.mso
r313 r325 59 59 60 60 VARIABLES 61 in Inlet as stream; 62 out Outlet as liquid_stream; 63 64 in InletQ as energy_stream (Brief="Rate of heat supply"); 61 in Inlet as stream (Brief = "Inlet stream", PosX=0.3037, PosY=0); 62 out Outlet as liquid_stream (Brief = "Outlet liquid stream", PosX=1, PosY=1); 63 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7859); 65 64 Level as length(Brief="Tank level"); 66 65 M(NComp) as mol (Brief="Molar Holdup in the tank"); … … 120 119 121 120 VARIABLES 122 in Inlet as stream; 123 out Outlet as liquid_stream; 124 125 in InletQ as energy_stream (Brief="Rate of heat supply"); 121 in Inlet as stream (Brief = "Inlet stream", PosX=0.1825, PosY=0); 122 out Outlet as liquid_stream (Brief = "Outlet liquid stream", PosX=1, PosY=1); 123 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.6160); 126 124 Level as length(Brief="Tank level"); 127 125 Across as area (Brief="Tank cross section area", Default=2); … … 176 174 VARIABLES 177 175 Level as length(Brief="Tank level"); 178 in Fin as flow_vol(Brief="Input flow" );179 out Fout as flow_vol(Brief="Output flow" );176 in Fin as flow_vol(Brief="Input flow", PosX=0.3037, PosY=0); 177 out Fout as flow_vol(Brief="Output flow", PosX=1, PosY=1); 180 178 181 179 EQUATIONS … … 211 209 212 210 VARIABLES 213 in Feed as stream ;214 in Inlet as stream;215 out Outlet as liquid_stream ;216 217 in InletQ as energy_stream (Brief="Rate of heat supply"); 211 in Feed as stream (Brief = "Feed stream", PosX=0.32, PosY=0); 212 in Inlet as stream (Brief = "Inlet stream", PosX=0.3037, PosY=0); 213 out Outlet as liquid_stream (Brief = "Outlet liquid stream", PosX=1, PosY=1); 214 in InletQ as energy_stream (Brief="Rate of heat supply", PosX=1, PosY=0.7859); 215 218 216 Level as length(Brief="Tank level"); 219 217 M(NComp) as mol (Brief="Molar Holdup in the tank"); -
trunk/eml/stage_separators/tray.mso
r300 r325 45 45 46 46 VARIABLES 47 in Inlet as stream ;48 in InletL as stream ;49 in InletV as stream ;50 out OutletL as liquid_stream ;51 out OutletV as vapour_stream ;47 in Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.4932); 48 in InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0); 49 in InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1); 50 out OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1); 51 out OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0); 52 52 53 53 M(NComp) as mol (Brief="Molar Holdup in the tray"); … … 202 202 203 203 VARIABLES 204 in Inlet as stream ;205 in InletL as stream ;206 in InletV as stream ;207 out OutletL as liquid_stream ;208 out OutletV as vapour_stream ;204 in Inlet as stream (Brief="Feed stream", PosX=0, PosY=0.4932); 205 in InletL as stream (Brief="Inlet liquid stream", PosX=0.5195, PosY=0); 206 in InletV as stream (Brief="Inlet vapour stream", PosX=0.4994, PosY=1); 207 out OutletL as liquid_stream (Brief="Outlet liquid stream", PosX=0.8277, PosY=1); 208 out OutletV as vapour_stream (Brief="Outlet vapour stream", PosX=0.8043, PosY=0); 209 209 210 210 yideal(NComp) as fraction; -
trunk/eml/streams.mso
r323 r325 40 40 z(NComp) as fraction(Brief = "Overall Molar Fraction"); 41 41 h as enth_mol; 42 v as fraction(Brief = "Vapouri sation fraction");42 v as fraction(Brief = "Vapourization fraction"); 43 43 end 44 44 … … 150 150 151 151 VARIABLES 152 out Outlet as stream (Brief = "Outlet stream", PosX=1, PosY=0.5 );152 out Outlet as stream (Brief = "Outlet stream", PosX=1, PosY=0.5256); 153 153 x(NComp) as fraction (Brief = "Liquid Molar Fraction"); 154 154 y(NComp) as fraction (Brief = "Vapour Molar Fraction"); … … 226 226 227 227 VARIABLES 228 out Outlet as stream (Brief = "Outlet stream", PosX=1, PosY=0.5 );228 out Outlet as stream (Brief = "Outlet stream", PosX=1, PosY=0.5256); 229 229 x(NComp) as fraction (Brief = "Liquid Molar Fraction"); 230 230 y(NComp) as fraction (Brief = "Vapour Molar Fraction"); … … 274 274 275 275 VARIABLES 276 in Inlet as stream ;277 v as fraction ;276 in Inlet as stream (Brief = "Inlet Stream", PosX=0, PosY=0.5308); 277 v as fraction (Brief = "Vapourization fraction"); 278 278 x(NComp) as fraction (Brief = "Liquid Molar Fraction"); 279 279 y(NComp) as fraction (Brief = "Vapour Molar Fraction"); … … 331 331 332 332 VARIABLES 333 in Inlet as stream (Brief = "Inlet stream", PosX=0.5, PosY=0);333 in Inlet as stream (Brief = "Inlet Stream", PosX=0, PosY=0.5308); 334 334 end 335 335 … … 354 354 355 355 VARIABLES 356 out OutletQ as energy_stream (Brief = "Outlet energy stream", PosX=1, PosY=0.5 );357 end 356 out OutletQ as energy_stream (Brief = "Outlet energy stream", PosX=1, PosY=0.5349); 357 end
Note: See TracChangeset
for help on using the changeset viewer.