Changeset 595 for branches/gui/eml/pressure_changers/compressor.mso

Timestamp:

Aug 8, 2008, 6:55:17 PM (15 years ago)

Author:

gerson bicca

Message:

updated compressor model (testing)

File:

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

branches/gui/eml/pressure_changers/compressor.mso

@@ -21 +21 @@
 
 Model centrifugal_compressor
-        ATTRIBUTES
+
+ATTRIBUTES
         Pallete         = true;
         Icon            = "icon/CentrifugalCompressor"; 
-        Brief           = "Model of a centrifugal compressor.";
-        Info            =
-"== Assumptions ==
-* Steady State;
-* Only Vapor;
-* Adiabatic.
+        Brief           = "Testing Model of a centrifugal compressor.";
+
         
-== Specify == 
-* the inlet stream;
-* the outlet pressure (Outlet.P);
-* the Isentropic efficiency (Effs).
-";
-        
-        PARAMETERS
-outer PP                as Plugin               (Brief = "External Physical Properties", Type="PP");
-outer NComp     as Integer              (Brief = "Number of chemical components", Lower = 1);
-        R                       as positive     (Default = 8.31451, Brief = "Constant of Gases", Unit= 'kJ/kmol/K');
+PARAMETERS
+
+outer PP                        as Plugin                       (Brief = "External Physical Properties", Type="PP");
+outer NComp     as Integer                      (Brief = "Number of chemical components", Lower = 1);
+        R                                       as positive             (Brief = "Constant of Gases", Unit= 'kJ/kmol/K', Default = 8.31451,Hidden=true);
         Mw(NComp)       as molweight    (Brief = "Molar Weight");
 
-        VARIABLES
-        n                       as positive             (Brief = "Politropic Coefficient", Lower=0);
-        k                       as positive     (Brief = "Isentropic Coefficient", Lower=1e-3); 
-        Cp              as cp_mol               (Brief = "Heat Capacity");
-        Cv                      as cv_mol               (Brief = "Heat Capacity");
+VARIABLES
+
+        n                               as positive             (Brief = "Politropic Coefficient", Lower=0);
+        k                               as positive             (Brief = "Isentropic Coefficient", Lower=1e-3); 
+        Cp              as cp_mol                       (Brief = "Heat Capacity");
+        Cv                      as cv_mol                       (Brief = "Heat Capacity");
         Pratio          as positive             (Brief = "Pressure Ratio", Symbol ="P_{ratio}");        
         Pdrop           as press_delta  (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
         Wp                      as energy_mol   (Brief = "Politropic Head");
         Ws                      as energy_mol   (Brief = "Isentropic Head");
-        Tiso            as temperature  (Brief = "Isentropic Temperature"); 
-        Effp            as positive     (Brief = "Politropic efficiency");
-        Effs            as efficiency   (Brief = "Isentropic efficiency");
-        FPower          as power                (Brief = "Fluid Power");
-        Mwm                     as molweight    (Brief = "Mixture Molar Weight");
-in      Inlet           as stream               (Brief = "Inlet stream", PosX=0, PosY=0.5086, Symbol="_{in}");
-out     Outlet          as streamPH             (Brief = "Outlet stream", PosX=1, PosY=0.5022, Symbol="_{out}");
+        Tiso                    as temperature          (Brief = "Isentropic Temperature"); 
+        Effp                    as positive             (Brief = "Politropic efficiency");
+        Effs                    as efficiency           (Brief = "Isentropic efficiency");
+        FPower          as power                        (Brief = "Fluid Power");
+        Mwm                     as molweight            (Brief = "Mixture Molar Weight");
 
-        SET
+in              Inlet           as stream               (Brief = "Inlet stream", PosX=0.35, PosY=0, Symbol="_{in}");
+out     Outlet          as streamPH     (Brief = "Outlet stream", PosX=1, PosY=0.85, Symbol="_{out}");
+
+in      WorkIn          as power                (Brief = "Work Inlet", PosX=0, PosY=0.55);
+
+SET
+
         Mw = PP.MolecularWeight();
-        
-        EQUATIONS
-        
-        "Calculate Mwm for Inlet Mixture"
+
+        R       = 8.31451*'kJ/kmol/K';
+
+EQUATIONS
+
+"Calculate Mwm for Inlet Mixture"
         Mwm = sum(Mw*Inlet.z);
 
-        "Pressure Ratio"
+"Pressure Ratio"
         Outlet.P = Inlet.P * Pratio;
 
-        "Pressure Drop"
+"Pressure Drop"
         Outlet.P  = Inlet.P - Pdrop;
 
-        "Calculate Cp Using a External Physical Properties Routine"
+"Calculate Cp Using a External Physical Properties Routine"
         Cp = PP.VapourCp(Inlet.T,Inlet.P,Inlet.z);
         
-        "Calculate Cv Using a External Physical Properties Routine"
+"Calculate Cv Using a External Physical Properties Routine"
         Cv = PP.VapourCv(Inlet.T,Inlet.P,Inlet.z);
         
-        "Calculate Isentropic Coeficient"
+"Calculate Isentropic Coeficient"
         k * Cv = Cp;
         
-        "Calculate Isentropic Head"
+"Calculate Isentropic Head"
         Ws = (k/(k-1))*R*Inlet.T*((Outlet.P/Inlet.P)^((k-1)/k) - 1);
         
-        "Calculate Isentropic Outlet Temperature"
+"Calculate Isentropic Outlet Temperature"
 #       Tiso = Inlet.T * (Outlet.P/Inlet.P)^((k-1)/k);
-        PP.VapourEntropy(Tiso, Outlet.P, Outlet.z) = 
-                PP.VapourEntropy(Inlet.T, Inlet.P, Inlet.z);
+        PP.VapourEntropy(Tiso, Outlet.P, Outlet.z) = PP.VapourEntropy(Inlet.T, Inlet.P, Inlet.z);
 
-        "Calculate Real Outlet Temperature"
+"Calculate Real Outlet Temperature"
         Effs * (Outlet.T- Inlet.T) = (Tiso - Inlet.T);
         
-        "Calculate Politropic Coefficient"
+"Calculate Politropic Coefficient"
         n*(ln(Outlet.T/Inlet.T)) = (n-1)*(ln(Outlet.P/Inlet.P));
         
-        "Calculate Politropic Efficiency"
+"Calculate Politropic Efficiency"
         Effp * (n-1) * k = n * (k-1);
         
-        "Calculate Politropic Head"
+"Calculate Politropic Head"
         Ws*Effp = Wp*Effs;
 
-        "Calculate Fluid Power"
-        FPower*Effs = Inlet.F*Ws;
-        
-        "Overall Molar Balance"
+"Overall Molar Balance"
         Outlet.F = Inlet.F;
 
-        "Component Molar Balance"
+"Component Molar Balance"
         Outlet.z = Inlet.z;
+
+# Testing Equations
+
+"Fluid Power"
+        FPower*Effs = Inlet.F*Ws;
+
+"Fluid Power"
+        FPower = WorkIn;
+
 end