Changeset 346 for trunk/eml/streams.mso

Timestamp:

Aug 19, 2007, 10:32:18 PM (15 years ago)

Author:

Argimiro Resende Secchi

Message:

Adding entropy in streamPH, sources and sinks.

File:

• trunk/eml/streams.mso (modified) (8 diffs)

trunk/eml/streams.mso

@@ -35 +35 @@
 
         VARIABLES
-        F as flow_mol;
-        T as temperature;
-        P as pressure;
-        z(NComp) as fraction(Brief = "Overall Molar Fraction");
-        h as enth_mol;
-        v as fraction(Brief = "Vapourization fraction");
+        F as flow_mol                   (Brief = "Stream Molar Flow Rate");
+        T as temperature                (Brief = "Stream Temperature");
+        P as pressure                   (Brief = "Stream Pressure");
+        z(NComp) as fraction    (Brief = "Stream Molar Fraction");
+        h as enth_mol                   (Brief = "Stream Enthalpy");
+        v as fraction                   (Brief = "Vapourization fraction");
 end
 
@@ -101 +101 @@
         
         VARIABLES
-        x(NComp) as fraction(Brief = "Liquid Molar Fraction");
-        y(NComp) as fraction(Brief = "Vapour Molar Fraction");
-        
+        x(NComp) as fraction    (Brief = "Liquid Molar Fraction");
+        y(NComp) as fraction    (Brief = "Vapour Molar Fraction");
+        s as entr_mol                   (Brief = "Stream Entropy");
+
         EQUATIONS
         "Flash Calculation"
         [v, x, y] = PP.FlashPH(P, h, z);
+        
         "Enthalpy"
         h = (1-v)*PP.LiquidEnthalpy(T, P, x) +
                 v*PP.VapourEnthalpy(T, P, y);
+        
+        "Entropy"
+        s = (1-v)*PP.LiquidEntropy(T, P, x) +
+                v*PP.VapourEntropy(T, P, y);
 end
 
@@ -153 +159 @@
         x(NComp)                        as fraction                     (Brief = "Liquid Molar Fraction");
         y(NComp)                        as fraction                     (Brief = "Vapour Molar Fraction");
-        hl                                      as enth_mol;
-        hv                                      as enth_mol;
+        hl                                      as enth_mol                     (Brief = "Liquid Enthalpy");
+        hv                                      as enth_mol                     (Brief = "Vapour Enthalpy");
+        s                                       as entr_mol                     (Brief = "Stream Entropy");
+        sl                                      as entr_mol                     (Brief = "Liquid Entropy");
+        sv                                      as entr_mol                     (Brief = "Vapour Entropy");     
         zmass(NComp)            as fraction                     (Brief = "Mass Fraction");
         Mw                                      as molweight            (Brief = "Average Mol Weight");
@@ -168 +177 @@
         
         "Overall Enthalpy"
-        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, x) +
-                Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
-        
+        Outlet.h = (1-Outlet.v)*hl + Outlet.v*hv;
+
+        "Liquid Enthalpy"
         hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, x);
+
+        "Vapour Enthalpy"
         hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
-        
+
+        "Overall Entropy"
+        s = (1-Outlet.v)*sl + Outlet.v*sv;
+
+        "Liquid Entropy"
+        sl = PP.LiquidEntropy(Outlet.T, Outlet.P, x);
+        
+        "Vapour Entropy"
+        sv = PP.VapourEntropy(Outlet.T, Outlet.P, y);
+
         "Average Molecular Weight"
         Mw = sum(M*Outlet.z);
@@ -229 +249 @@
         x(NComp)                        as fraction             (Brief = "Liquid Molar Fraction");
         y(NComp)                        as fraction             (Brief = "Vapour Molar Fraction");
-        hl                                      as enth_mol;
-        hv                                      as enth_mol;
+        hl                                      as enth_mol             (Brief = "Liquid Enthalpy");
+        hv                                      as enth_mol             (Brief = "Vapour Enthalpy");
+        s                                       as entr_mol             (Brief = "Stream Entropy");
+        sl                                      as entr_mol             (Brief = "Liquid Entropy");
+        sv                                      as entr_mol             (Brief = "Vapour Entropy");     
 
         EQUATIONS
@@ -237 +260 @@
         
         "Overall Enthalpy"
-        Outlet.h = (1-Outlet.v)*PP.LiquidEnthalpy(Outlet.T, Outlet.P, x) +
-                Outlet.v*PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
-
+        Outlet.h = (1-Outlet.v)*hl + Outlet.v*hv;
+
+        "Liquid Enthalpy"
         hl = PP.LiquidEnthalpy(Outlet.T, Outlet.P, x);
+
+        "Vapour Enthalpy"
         hv = PP.VapourEnthalpy(Outlet.T, Outlet.P, y);
+        
+        "Overall Entropy"
+        s = (1-Outlet.v)*sl + Outlet.v*sv;
+
+        "Liquid Entropy"
+        sl = PP.LiquidEntropy(Outlet.T, Outlet.P, x);
+        
+        "Vapour Entropy"
+        sv = PP.VapourEntropy(Outlet.T, Outlet.P, y);
 end
 
@@ -285 +319 @@
         Fw                              as flow_mass    (Brief = "Stream Mass Flow");
         Fvol            as flow_vol     (Brief = "Volumetric Flow");
-        
+        s                               as entr_mol             (Brief = "Stream Entropy");
+
         EQUATIONS
         "Flash Calculation"
@@ -318 +353 @@
         zmass = M*Inlet.z / Mw;
         
+        "Overall Entropy"
+        s = (1-v)*PP.LiquidEntropy(Inlet.T, Inlet.P, x) +
+                v*PP.VapourEntropy(Inlet.T, Inlet.P, y);
 end