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

Timestamp:

Oct 19, 2009, 5:41:09 PM (14 years ago)

Author:

gerson bicca

Message:

revised compressor model

File:

• branches/gui/eml/pressure_changers/compressor.mso (modified) (6 diffs)

branches/gui/eml/pressure_changers/compressor.mso

@@ -42 +42 @@
 outer NComp                     as Integer              (Brief = "Number of chemical components", Lower = 1);
         Rgas                                            as positive             (Brief = "Constant of Gases", Unit= 'kJ/kmol/K', Default = 8.31451,Hidden=true);
-        Mw(NComp)                       as molweight    (Brief = "Molar Weight");
-        CompressorType  as Switcher             (Brief = "Compressor Model Type",Valid=["Polytropic With GPSA Method","Isentropic With GPSA Method","Isentropic With ASME Method","Polytropic With ASME Method"], Default="Isentropic With GPSA Method");
+        Mw(NComp)                       as molweight    (Brief = "Molar Weight",Hidden=true);
+        CompressorOperation     as Switcher             (Brief = "Compressor Operation Model",Valid=["Polytropic","Isentropic"], Default="Isentropic");
+        Method  as Switcher             (Brief = "Method of Calculation",Valid=["GPSA Method","ASME Method"], Default="GPSA Method");
 
 VARIABLES
@@ -148 +149 @@
         FluidPower = Head*Mwm*Inlet.F;
 
-switch CompressorType
-
-        case "Isentropic With GPSA Method":
+switch CompressorOperation
+
+        case "Isentropic":
+        
+switch Method
+        
+        case "ASME Method":
+
+"Efficiency"
+        EfficiencyOperation = IsentropicEff;
+
+"Isentropic Coefficient"
+        IseCoeff*ln(rho_ise/rho_in) = ln(Outlet.P/Inlet.P);
+
+"Polytropic Coefficient"
+        PolyCoeff*ln(rho_out/rho_in) = ln(Outlet.P/Inlet.P);
+
+"Isentropic Head"
+        HeadIsentropic*rho_in = (IseCoeff/(IseCoeff-1.001))*Inlet.P*HeadCorrection*((Pratio)^((IseCoeff-1.001)/IseCoeff) - 1);
+
+"Polytropic Head"
+        HeadPolytropic*rho_in = (PolyCoeff/(PolyCoeff-1.001))*Inlet.P*HeadCorrection*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
+
+"Schultz Polytropic Head Correction"
+        HeadCorrection*Mwm*(IseCoeff/(IseCoeff-1.001))*(Outlet.P/rho_ise -Inlet.P/rho_in) = (hise-Inlet.h);
+        
+        case "GPSA Method":
 
 "Efficiency"
@@ -169 +194 @@
 "Isentropic Head"
         HeadIsentropic = Head*IsentropicEff;
-        
-        case "Polytropic With GPSA Method":
+
+end
+
+        case "Polytropic":
+
+switch Method
+
+        case "GPSA Method":
 
 "Efficiency"
@@ -190 +221 @@
         HeadIsentropic = Head*IsentropicEff;
         
-        case "Isentropic With ASME Method":
-
-"Efficiency"
-        EfficiencyOperation = IsentropicEff;
+        case "ASME Method":
+
+"Efficiency"
+        EfficiencyOperation = PolytropicEff;
 
 "Isentropic Coefficient"
@@ -206 +237 @@
 "Polytropic Head"
         HeadPolytropic*rho_in = (PolyCoeff/(PolyCoeff-1.001))*Inlet.P*HeadCorrection*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
-
-"Schultz Polytropic Head Correction"
-        HeadCorrection*Mwm*(IseCoeff/(IseCoeff-1.001))*(Outlet.P/rho_ise -Inlet.P/rho_in) = (hise-Inlet.h);
-
-        case "Polytropic With ASME Method":
-
-"Efficiency"
-        EfficiencyOperation = PolytropicEff;
-
-"Isentropic Coefficient"
-        IseCoeff*ln(rho_ise/rho_in) = ln(Outlet.P/Inlet.P);
-
-"Polytropic Coefficient"
-        PolyCoeff*ln(rho_out/rho_in) = ln(Outlet.P/Inlet.P);
-
-"Isentropic Head"
-        HeadIsentropic*rho_in = (IseCoeff/(IseCoeff-1.001))*Inlet.P*HeadCorrection*((Pratio)^((IseCoeff-1.001)/IseCoeff) - 1);
-
-"Polytropic Head"
-        HeadPolytropic*rho_in = (PolyCoeff/(PolyCoeff-1.001))*Inlet.P*HeadCorrection*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
         
 "Schultz Polytropic Head Correction"
@@ -232 +243 @@
 
 end
+
+end