Changeset 605 for branches/gui/eml

Timestamp:

Aug 22, 2008, 1:16:35 AM (15 years ago)

Author:

gerson bicca

Message:

expander model finished (missing documentation)

File:

• branches/gui/eml/pressure_changers/turbine.mso (modified) (4 diffs)

branches/gui/eml/pressure_changers/turbine.mso

@@ -158 +158 @@
         FluidPower      as power                (Brief = "Fluid Power");
         BrakePower      as power                (Brief = "Brake Power");
-        PowerLoss       as power                (Brief = "Power Losses");
+        PowerLoss       as power                (Brief = "Power Losses",Lower=0);
         Mwm                     as molweight    (Brief = "Mixture Molar Weight");
         rho                     as dens_mass    (Brief = "Mass Density");
@@ -207 +207 @@
         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);
@@ -217 +223 @@
         Outlet.T = Tisentropic;
 
-"Head"
-        Head = HeadIsentropic;
-
-"Isentropic Coefficient"
-        IseCoeff= 1;
-
         else
-"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);
-        Head = (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);
 
 "Discharge Temperature"
-        Outlet.T = Inlet.T*((Outlet.P/Inlet.P)^((IseCoeff-1.001)/IseCoeff));
-
-"Head"
-        Head*Mwm = (PP.VapourEnthalpy(Outlet.T,Outlet.P,Outlet.z)-Inlet.h);
+        (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 = Head*sum(Mw*Inlet.z)*Inlet.F;
-        #FluidPower = Head*sum(Mw*Inlet.z)*Inlet.F;
+        FluidPower = Head*sum(Mw*Inlet.z)*Inlet.F;
 
 "Brake Power"
@@ -244 +237 @@
 
 "Brake Power"
-        BrakePower = (FluidPower/MechanicalEff)-PowerLoss;
+        BrakePower = FluidPower*MechanicalEff;
 
 "Power Loss"