Changeset 571 for branches/gui/eml/streams.mso

Timestamp:

Jul 24, 2008, 8:02:39 PM (15 years ago)

Author:

gerson bicca

Message:

testing a new way to specify a stream composition

File:

• branches/gui/eml/streams.mso (modified) (1 diff)

branches/gui/eml/streams.mso

@@ -531 +531 @@
 
 end
+
+Model source_testing
+
+#Testing a new way to specify stream composition !!!!
+# Must be added a set of equations to solve the  "Switch CompostionBasis" !!!!!!
+
+        ATTRIBUTES
+        Pallete = true;
+        Icon = "icon/Source";
+        Brief = "Material stream source";
+        Info = "
+        This model should be used for boundary streams.
+        Usually these streams are known and come from another process
+        units.
+
+        The user should specify:
+         * Total molar (mass or volumetric) flow
+         * Temperature
+         * Pressure
+         * Molar (mass or volumetric) composition
+        
+        No matter the specification set, the model will calculate some
+        additional properties:
+         * Mass density
+         * Mass flow
+         * Mass compostions
+         * Specific volume
+         * Vapour fraction
+         * Volumetric flow
+         * Liquid and Vapour compositions
+        ";
+
+        PARAMETERS
+        outer PP                        as Plugin               (Brief = "External Physical Properties", Type="PP");
+        outer NComp             as Integer              (Brief = "Number of chemical components", Lower = 1); 
+                  M(NComp)      as molweight    (Brief = "Component Mol Weight");
+                  rhoModel              as Switcher             (Brief = "Density model", Valid = ["volume", "correlation"], Default="volume");
+                  CompostionBasis               as Switcher             (Brief = "Molar or Mass Compostion", Valid = ["Molar", "Mass"], Default="Molar");
+        
+        Composition(NComp) as fraction  (Brief = "Stream Composition");
+        
+        SET
+
+        M   = PP.MolecularWeight();
+
+        VARIABLES
+        out Outlet                      as stream                       (Brief = "Outlet stream", PosX=1, PosY=0.5256, Symbol="_{out}");
+        x(NComp)                        as fraction                     (Brief = "Liquid Molar Fraction",Hidden=true);
+        y(NComp)                        as fraction                     (Brief = "Vapour Molar Fraction",Hidden=true);
+        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",Protected=true);
+        Mw                                      as molweight            (Brief = "Average Mol Weight");
+        vm                                      as volume_mol           (Brief = "Molar Volume");       
+        rho                                     as dens_mass            (Brief = "Stream Mass Density");
+        rhom                            as dens_mol                     (Brief = "Stream Molar Density");
+        Fw                                      as flow_mass            (Brief = "Stream Mass Flow");
+        Fvol                    as flow_vol         (Brief = "Volumetric Flow");
+        T_Cdeg                          as temperature          (Brief = "Temperature in °C", Lower=-200);
+        
+        EQUATIONS
+
+switch CompostionBasis
+
+        case "Molar":
+"Stream Molar Composition"
+        Outlet.z = Composition/sum(Composition); 
+
+"Stream Mass Composition"
+        zmass = M*Outlet.z / Mw;
+
+        case "Mass":
+"Stream Mass Composition"
+        zmass = Composition/sum(Composition);
+
+"Stream Molar Composition"
+        Outlet.z = (Composition/sum(Composition)/M)*sum(M*Outlet.z);# do not converge fixme !!!!!
+
+end
+
+        "Flash Calculation"
+        [Outlet.v, x, y] = PP.Flash(Outlet.T, Outlet.P, Outlet.z);
+        
+        "Overall Enthalpy"
+        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);
+
+        switch rhoModel
+                case "volume":
+        "Molar Density"
+                rhom * vm = 1;
+                
+                case "correlation":
+        "Mass Density"
+                rho*((1-Outlet.v)/PP.LiquidDensity(Outlet.T,Outlet.P,x) + Outlet.v/PP.VapourDensity(Outlet.T,Outlet.P,y)) = 1;
+        end
+        
+        "Mass or Molar Density"
+        rhom * Mw = rho;
+
+        "Flow Mass"
+        Fw      =  Mw*Outlet.F;
+
+        "Molar Volume"
+        vm = (1-Outlet.v)*PP.LiquidVolume(Outlet.T, Outlet.P, x) + Outlet.v*PP.VapourVolume(Outlet.T,Outlet.P,y);
+        
+        "Volumetric Flow"
+        Fvol = Outlet.F*vm ;
+        
+        "Temperature in °C"
+        T_Cdeg = Outlet.T - 273.15 * 'K';
+
+end