Changeset 689 for branches/gui/eml

Timestamp:

Nov 20, 2008, 10:56:29 PM (14 years ago)

Author:

gerson bicca

Message:

updates

Location:

branches/gui/eml/pressure_changers

Files:

• compressor.mso (modified) (2 diffs)
• expander.mso (added)
• icon/expander.png (added)
• icon/expander.svg (added)
• turbine.mso (modified) (1 diff)

branches/gui/eml/pressure_changers/compressor.mso

@@ -159 +159 @@
 
 "Isentropic Coefficient"
-        HeadIsentropic = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(IseCoeff/(IseCoeff-1.001))*Rgas*Inlet.T*((Pratio)^((IseCoeff-1.001)/IseCoeff) - 1);
-
-"Polytropic Coefficient"
-        HeadPolytropic = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(PolyCoeff/(PolyCoeff-1.001))*Rgas*Inlet.T*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
+        HeadIsentropic*Mwm*((IseCoeff-1.001)/IseCoeff) = (0.5*Zfac_in+0.5*Zfac_out)*Rgas*Inlet.T*((Pratio)^((IseCoeff-1.001)/IseCoeff) - 1);
+
+"Polytropic Coefficient"
+        HeadPolytropic*Mwm*((PolyCoeff-1.001)/PolyCoeff) = (0.5*Zfac_in+0.5*Zfac_out)*Rgas*Inlet.T*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
 
 "Head Correction"
@@ -179 +179 @@
 
 "Isentropic Coefficient"
-        HeadIsentropic = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(IseCoeff/(IseCoeff-1.001))*Rgas*Inlet.T*((Pratio)^((IseCoeff-1.001)/IseCoeff) - 1);
-
-"Polytropic Coefficient"
-        HeadPolytropic = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(PolyCoeff/(PolyCoeff-1.001))*Rgas*Inlet.T*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
+        HeadIsentropic*Mwm*((IseCoeff-1.001)/IseCoeff) = (0.5*Zfac_in+0.5*Zfac_out)*Rgas*Inlet.T*((Pratio)^((IseCoeff-1.001)/IseCoeff) - 1);
+
+"Polytropic Coefficient"
+        HeadPolytropic*Mwm*((PolyCoeff-1.001)/PolyCoeff) = (0.5*Zfac_in+0.5*Zfac_out)*Rgas*Inlet.T*((Pratio)^((PolyCoeff-1.001)/PolyCoeff) - 1);
 
 "Head Correction"

branches/gui/eml/pressure_changers/turbine.mso

@@ -126 +126 @@
 end
 
-Model expander
-        
-ATTRIBUTES
-        Pallete         = true;
-        Icon            = "icon/HidraulicTurbine"; 
-        Brief           = "Model of an expansor.";
-        Info            =
-"To be documented";
-        
-PARAMETERS
-
-outer PP                as Plugin               (Brief = "External Physical Properties", Type="PP");
-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");
-
-VARIABLES
-
-        IseCoeff        as positive     (Brief = "Isentropic Coefficient", Lower=0.2); 
-        Pratio          as positive             (Brief = "Pressure Ratio", Symbol ="P_{ratio}");        
-        Pdrop           as press_delta  (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
-        Pdecrease       as press_delta  (Brief = "Pressure Decrease", DisplayUnit = 'kPa', Symbol ="P_{decr}");
-
-        Head                            as energy_mass  (Brief = "Head",Hidden=true);
-        HeadIsentropic          as energy_mass  (Brief = "Isentropic Head");
-        Tisentropic                     as temperature  (Brief = "Isentropic Temperature"); 
-        
-        IsentropicEff   as efficiency   (Brief = "Isentropic efficiency");
-        MechanicalEff   as efficiency   (Brief = "Mechanical efficiency");
-        
-        FluidPower      as power                (Brief = "Fluid Power");
-        BrakePower      as power                (Brief = "Brake Power");
-        PowerLoss       as power                (Brief = "Power Losses",Lower=0);
-        Mwm                     as molweight    (Brief = "Mixture Molar Weight");
-        rho                     as dens_mass    (Brief = "Mass Density");
-        Zfac_in         as fraction     (Brief = "Compressibility factor at inlet");
-        Zfac_out        as fraction     (Brief = "Compressibility factor at outlet");
-
-in      Inlet           as stream       (Brief = "Inlet stream", PosX=0.05, PosY=0.0, Symbol="_{in}");
-out     Outlet          as streamPH     (Brief = "Outlet stream", PosX=0.65, PosY=1, Symbol="_{out}");
-
-out     WorkOut         as power        (Brief = "Work Outlet", PosX=1, PosY=0.46);
-
-SET
-
-        Mw = PP.MolecularWeight();
-
-        Rgas    = 8.31451*'kJ/kmol/K';
-
-EQUATIONS
-
-"Overall Molar Balance"
-        Outlet.F = Inlet.F;
-
-"Component Molar Balance"
-        Outlet.z = Inlet.z;
-
-"Average Molecular Weight"
-        Mwm = sum(Mw*Inlet.z);
-
-"Pressure Ratio"
-        Outlet.P = Inlet.P * Pratio;
-
-"Pressure Drop"
-        Outlet.P  = Inlet.P - Pdrop;
-
-"Pressure Decrease"
-        Outlet.P  = Inlet.P - Pdecrease;
-
-"Mass Density"
-        rho = PP.VapourDensity(Inlet.T, Inlet.P, Inlet.z);
-        
-"Compressibility factor at Inlet Conditions"
-        Zfac_in = PP.VapourCompressibilityFactor(Inlet.T,Inlet.P,Inlet.z);
-        
-"Compressibility factor at Outlet Conditions"
-        Zfac_out = PP.VapourCompressibilityFactor(Outlet.T,Outlet.P,Outlet.z);
-
-"Isentropic Head"
-        HeadIsentropic*Mwm = (PP.VapourEnthalpy(Tisentropic,Outlet.P,Outlet.z)-Inlet.h);
-
-"Actual Head"
-        Head*Mwm = (Outlet.h-Inlet.h);
-
-"Isentropic Coefficient"
-        HeadIsentropic = (0.5*Zfac_in+0.5*Zfac_out)*(1/Mwm)*(IseCoeff/(IseCoeff-1.001))*Rgas*Inlet.T*((Outlet.P/Inlet.P)^((IseCoeff-1.001)/IseCoeff) - 1);
-
-"Isentropic Outlet Temperature"
-        PP.VapourEntropy(Tisentropic, Outlet.P, Outlet.z) = PP.VapourEntropy(Inlet.T, Inlet.P, Inlet.z);
-
-
-if IsentropicEff equal 1
-
-        then
-"Discharge Temperature"
-        Outlet.T = Tisentropic;
-
-        else
-
-"Discharge Temperature"
-        (PP.VapourEnthalpy(Outlet.T,Outlet.P,Outlet.z)-Inlet.h)= (PP.VapourEnthalpy(Tisentropic,Outlet.P,Outlet.z)-Inlet.h)*IsentropicEff;
-
-end
-
-"Fluid Power"
-        FluidPower = IsentropicEff*HeadIsentropic*sum(Mw*Inlet.z)*Inlet.F+PowerLoss;
-
-"Brake Power"
-        BrakePower = WorkOut;
-
-"Brake Power"
-        BrakePower = FluidPower*MechanicalEff;
-
-"Power Loss"
-        PowerLoss = BrakePower - FluidPower;
-
-end