Changeset 597

Timestamp:

Aug 12, 2008, 5:33:05 PM (15 years ago)

Author:

gerson bicca

Message:

added expander model (testing)

Location:

branches/gui

Files:

• eml/pressure_changers/turbine.mso (modified) (2 diffs)
• sample/pressure_changers/sample_compressor.mso (modified) (1 diff)

branches/gui/eml/pressure_changers/turbine.mso

@@ -98 +98 @@
 end
 
-
 Model HidraulicTurbineGenerator as HidraulicTurbine
         ATTRIBUTES
@@ -126 +125 @@
         EPower = BPower * Meff;
 end
+
+Model expander
+        
+        ATTRIBUTES
+        Pallete         = true;
+        Icon            = "icon/HidraulicTurbine"; 
+        Brief           = "Testing Model of a expander.";
+        
+PARAMETERS
+
+outer NComp     as Integer                      (Brief = "Number of chemical components", Lower = 1);
+outer PP                as Plugin                       (Brief = "External Physical Properties", Type="PP");
+        Mw(NComp)       as molweight            (Brief = "Molar Weight");
+        
+VARIABLES
+        Eff     as efficiency           (Brief = "Turbine efficiency");
+        Meff    as efficiency           (Brief = "Brake efficiency");
+        Beta    as positive             (Brief = "Volumetric expansivity", Unit = '1/K');
+        Head    as head                         (Brief = "Head Developed");
+        FPower  as power                        (Brief = "Fluid Power");
+        BPower  as power                        (Brief = "Brake Power");
+        Pratio  as positive                     (Brief = "Pressure Ratio");
+        Pdrop   as press_delta          (Brief = "Pressure Drop", DisplayUnit = 'kPa', Symbol ="\Delta P");
+        Mwm     as molweight            (Brief = "Mixture Molar Weight");
+        rho             as dens_mass            (Brief = "Specific Mass");
+        Cp              as cp_mol                       (Brief = "Heat Capacity");
+        
+in      Inlet           as stream       (Brief = "Inlet stream", PosX=0.05, PosY=0.0, Symbol="_{in}");
+out     Outlet          as stream       (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();      
+        
+EQUATIONS
+
+"Calculate Mwm for Inlet Mixture"
+        Mwm = sum(Mw*Inlet.z);
+
+"Calculate rho using a External Physical Properties Routine"
+        rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
+        
+"Calculate Outlet Vapour Fraction"
+        Outlet.v = PP.VapourFraction(Outlet.T, Outlet.P, Outlet.z);
+        
+"Calculate Cp Using a External Physical Properties Routine"
+        Cp = PP.LiquidCp(Inlet.T,Inlet.P,Inlet.z);
+        
+"Pressure Ratio"
+        Outlet.P = Inlet.P * Pratio;
+
+"Pressure Drop"
+        Outlet.P  = Inlet.P - Pdrop;
+
+"Calculate Fluid Power"
+        FPower * rho = -Pdrop * Inlet.F * Mwm;
+
+"Calculate Brake Power"
+        BPower = FPower * Eff;
+
+        BPower = WorkOut;
+
+"Calculate Outlet Temperature"
+        (Outlet.T - Inlet.T) * rho * Cp = (Outlet.h - Inlet.h) * rho
+        + Pdrop * Mwm * (1-Beta*Inlet.T);
+
+"Calculate Outlet Enthalpy"
+        (Outlet.h - Inlet.h) * rho =  -Pdrop * Mwm;
+
+"Molar Balance"
+        Outlet.F = Inlet.F;
+"Outlet Composition"
+        Outlet.z = Inlet.z;
+
+"Calculate Head"
+        Head = Outlet.h - Inlet.h;
+
+end

branches/gui/sample/pressure_changers/sample_compressor.mso

@@ -76 +76 @@
         DEVICES
         S1 as source;
+        W1 as work_source;
         C1 as centrifugal_compressor;
         
         CONNECTIONS
         S1.Outlet to C1.Inlet;
+        W1.Work to C1.WorkIn;