Changeset 744

Timestamp:

Mar 6, 2009, 7:56:26 PM (15 years ago)

Author:

gerson bicca

Message:

improved SepComp? model (testing)

Location:

branches/gui

Files:

• eml/mixers_splitters/icon/SepComp.png (added)
• eml/mixers_splitters/icon/SepComp.svg (added)
• eml/mixers_splitters/sepComp.mso (modified) (1 diff)
• sample/mixers_splitters/Sample_SepComp_Diagram.pfd (added)
• sample/mixers_splitters/sample_sepComp.mso (deleted)

branches/gui/eml/mixers_splitters/sepComp.mso

@@ -114 +114 @@
         ATTRIBUTES
         Pallete         = true;
-        Icon            = "icon/splitter"; 
+        Icon            = "icon/SepComp"; 
         Brief           = "Model of a separator of components";
         Info            =
 "== Assumptions ==
 * thermodynamics equilibrium
-* adiabatic
         
 == Specify ==
 * the inlet stream
-* (NComp - 1) molar fractions to 1 of the outlet streams
-* the fraction of split of the outlet streams
+* (NComp) Overhead_Splits or (NComp) Bottom_Splits 
+* the Pressure and Temperature of the outlet streams
 ";
 
 PARAMETERS
 
-outer PP                        as Plugin       (Brief = "External Physical Properties", Type="PP");
-outer   NComp           as Integer      (Brief = "Number of chemical components", Lower = 1); 
-                mainComp        as Integer      (Brief = "Component specified", Default = 1, Lower = 1);
+outer   PP                      as Plugin       (Brief = "External Physical Properties", Type="PP");
+outer   NComp   as Integer (Brief = "Number of chemical components"); 
         
 VARIABLES
 
-in              Inlet           as stream (Brief = "Inlet stream", PosX=0, PosY=0.5001, Symbol="_{in}");
-out     Outlet1         as stream (Brief = "Outlet stream 1", PosX=1, PosY=0.3027, Symbol="_{out1}");
-out     Outlet2         as stream (Brief = "Outlet stream 2", PosX=1, PosY=0.7141, Symbol="_{out2}");
+in              Inlet                                           as stream               (Brief = "Inlet stream", PosX=0, PosY=0.5001, Symbol="_{Inlet}");
+out     Overhead_Outlet         as streamPH     (Brief = "Overhead_Outlet stream", PosX=0.5, PosY=0, Symbol="_{Overhead}");
+out     Bottom_Outlet           as streamPH     (Brief = "Bottom_Outlet stream", PosX=0.5, PosY=1, Symbol="_{Bottom}");
+in              InletQ                                  as power                (Brief="Rate of heat supply", PosX=1, PosY=0.7559, Symbol="_{in}");
+
+        Overhead_Splits(NComp)  as fraction     (Brief = "Fraction to Overhead_Outlet", Symbol="\phi_{overhead}");
+        Bottom_Splits(NComp)            as fraction     (Brief = "Fraction to Bottom_Outlet", Symbol="\phi_{bottom}");
         
-        frac                    as fraction (Brief = "Fraction to Outlet 1", Symbol="\phi");
-        recovery                as fraction (Brief = "Recovery of the component specified", Symbol="\eta");
+EQUATIONS
 
-        EQUATIONS
-        
-"Flow"
-        Outlet1.F = Inlet.F * frac;
-        Outlet1.F + Outlet2.F = Inlet.F;
-        
-        recovery*Inlet.z(mainComp) = frac*Outlet1.z(mainComp);
-        
-        sum(Outlet1.z) = 1;
-        
-for i in [1:NComp] do
-        
- "Composition"
-        Outlet1.F*Outlet1.z(i) + Outlet2.F*Outlet2.z(i) = Inlet.F*Inlet.z(i);
+"Composition Overhead"
+        Overhead_Outlet.F*Overhead_Outlet.z = Inlet.F * Inlet.z*Overhead_Splits;
 
-end
-        
-"Pressure"
-        Outlet1.P = Inlet.P;
-        Outlet2.P = Inlet.P;
-        
-"Enthalpy"
-        Outlet1.h = (1-Outlet1.v)*PP.LiquidEnthalpy(Outlet1.T, Outlet1.P, Outlet1.z) + 
-                                Outlet1.v*PP.VapourEnthalpy(Outlet1.T, Outlet1.P, Outlet1.z);
-        Outlet2.h = (1-Outlet2.v)*PP.LiquidEnthalpy(Outlet2.T, Outlet2.P, Outlet2.z) + 
-                                Outlet2.v*PP.VapourEnthalpy(Outlet2.T, Outlet2.P, Outlet2.z);
-        
-"Temperature"
-        Outlet1.T = Inlet.T;
-        Outlet2.T = Inlet.T;
-        
-"Vapourization Fraction"
-        Outlet1.v = PP.VapourFraction(Outlet1.T, Outlet1.P, Outlet1.z);
-        Outlet2.v = PP.VapourFraction(Outlet2.T, Outlet2.P, Outlet2.z);
+"Bottom Composition Constraints"
+        sum(Bottom_Outlet.z) = 1;
+
+"Sum Of Splits"
+        Overhead_Splits+Bottom_Splits = 1;
+
+"Overhead Composition Constraints"
+        sum(Overhead_Outlet.z) = 1;
+
+ "Global Composition"
+        Overhead_Outlet.F*Overhead_Outlet.z + Bottom_Outlet.F*Bottom_Outlet.z = Inlet.F*Inlet.z;
+
+"Energy Balance"
+        Inlet.F*Inlet.h = Overhead_Outlet.F*Overhead_Outlet.h + Bottom_Outlet.F*Bottom_Outlet.h - InletQ;
+
+#"Overhead Pressure"
+#       Overhead_Outlet.P = Inlet.P;
+
+#"Bottom Pressure"
+#       Bottom_Outlet.P = Inlet.P;
+
+#"Overhead Temperature"
+#       Overhead_Outlet.T = Inlet.T;
+
+#"Bottom Temperature"
+#       Bottom_Outlet.T = Inlet.T;
 
 end