Changeset 57 for mso/eml/pressure_changers/pump.mso

Timestamp:

Nov 16, 2006, 2:46:48 PM (16 years ago)

Author:

Paula Bettio Staudt

Message:

Added new pressure_changers models

File:

• mso/eml/pressure_changers/pump.mso (modified) (2 diffs)

mso/eml/pressure_changers/pump.mso

@@ +1 @@
+#*------------------------------------------------------------------------
+* This file is property of the author and cannot be used, copyed
+* or modified without permission.
+*
+* Copyright (C) 2002-2006  the author
+*-------------------------------------------------------------------------
+* Authors:      Andrey Copat            Estefane da Silveira Horn
+* $Id$                                          Marcos Lovato Alencastro        
+*                                                                                                               Date:   20/02/2006
+*-------------------------------------------------------------------------
+* -> Steady State
+* -> Only Liquid
+* -> Adiabatic
+* -> Isentropic
+*-------------------------------------------------------------------------*#
+using "streams";
+using "pressure_changers/flux_machine_basic";
+
+Model centrifugal_pump as flux_machine_basic
+        
+        PARAMETERS
+ext NComp       as Integer                      (Brief = "Number of chemical components", Lower = 1);
+ext PP                  as CalcObject           (Brief = "External Physical Properties");
+        Mw(NComp)       as molweight            (Brief = "Molar Weight");
+        Eff             as positive             (Default = 0.72, Brief = "Pump Efficiency");
+        Meff            as positive             (Default = 0.95, Brief = "Brake Efficiency");
+        Beta            as positive             (Default = 0, Brief = "Volumetric Expansivity", Unit = "1/K");
+        g                       as acceleration         (Brief = "Gravity Acceleration", Default = 9.81);
+        N                       as vel_angular          (Brief = "Rotation", Default = 100);
+        Lev                     as length                       (Brief = "Loss Friction", Default = 0);
+        
+        VARIABLES
+        rho               as dens_mass          (Brief = "Specific Mass", Unit="kg/m^3");
+        Cp        as cp_mol                     (Brief = "Heat Capacity", Unit="kJ/kmol/K");
+        FPower    as power                      (Brief = "Fluid Power", Unit="kW");
+        BPower    as power                      (Brief = "Brake Power",Unit="kW");
+        EPower    as power                      (Brief = "Eletrical Potency", Unit="kW");
+        Pdiff     as press_delta        (Brief = "Pressure Increase", Unit="kPa");
+        Pratio    as positive           (Brief = "Pressure Ratio");     
+        Head      as head                       (Brief = "Head Developed", Unit="kJ/kmol"); 
+        Head_is   as head                       (Brief = "Isoentripic Head", Unit="kJ/kmol"); 
+        Mwm       as molweight          (Brief = "Mixture Molar Weight");
+        pvm               as pressure           (Brief = "Mixture Vapour Pressure", Unit = "kPa");                      
+        NPSHa     as length                     (Brief = "Available Net Positive Suction Head");
+        NS                as positive           (Brief = "Specific Speed", Unit = "(rpm*(gal/min)^0.5)/(m^3/4)");
+        Q                 as flow_vol           (Brief = "Volumetric Flow Rate");
+        vm                as vol_mol            (Brief = "Mixture Molar Volume", Unit = "m^3/kmol");
+        
+        SET
+        Mw = PP.MolecularWeight();
+        
+        EQUATIONS
+        #Mixtures Properties
+        "Calculate Mwm for Inlet Mixture"
+        Mwm = sum(Mw([1:NComp])*Inlet.z([1:NComp]));
+
+        "Calculate Cp Using a External Physical Properties Routine"
+        Cp = PP.LiquidCp(Inlet.T,Inlet.P,Inlet.z);
+        
+        "Calculate rho using a External Physical Properties Routine"
+        rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
+
+        "Calculate Mixture Vapour Pressure"
+        [pvm] = PP.BubbleP(Inlet.T,Inlet.z); 
+        
+        "Calculate Outlet Vapour Fraction"
+        Outlet.v = PP.VapourFraction(Outlet.T, Outlet.P, Outlet.z);
+        
+        "Calculate Liquid Molar Volume"
+        vm = PP.LiquidVolume(Inlet.T,Inlet.P,Inlet.z);
+        
+        #Mass and Energy Balance and Pump Equations     
+        "Calculate Outlet Stream Pressure"
+        Outlet.P = Inlet.P + Pdiff;
+        
+        "Pratio Definition"
+        Outlet.P = Inlet.P * Pratio;
+        
+        "Calculate Isentropic Head"
+        Head_is = Pdiff * Mwm/rho;
+
+        "Calculate Real Head"
+        Head = Head_is/(Eff*Meff); 
+        
+        "Calculate Outlet Enthalpy"
+        Outlet.h - Inlet.h = Head;
+        
+        "Calculate Fluid Power"
+        FPower = Head_is * Inlet.F;
+
+        "Calculate Brake Power"
+        BPower * Eff = FPower;
+
+        "Calculate Eletric Power"
+        BPower = EPower * Meff;
+        
+        "Calculate Outlet Temperature"
+        (Outlet.T - Inlet.T) * Cp = (Outlet.h - Inlet.h) -  Pdiff * Mwm / rho * (1 - Beta * Inlet.T);
+        
+        "Molar Balance"
+        Outlet.F = Inlet.F;
+
+        "Calculate Outlet Composition"
+        Outlet.z = Inlet.z;
+        
+        "Calculate Net Positive Suction Head"
+        NPSHa = - Lev + (Inlet.P - pvm)/(g*rho); #If Inlet.P is the suction pump pressure, Lev is 0.
+        
+        "Calculate Volumetric Flow Rate"
+        Q = Inlet.F*vm;
+        
+        "Calculate Specific Speed"
+        NS = N*(Q^0.5)/(Head^3/4);
+        
+end
+
+
 #*-------------------------------------------------------------------
 * Model of a pump
@@ -15 +132 @@
 * $Id$
 *--------------------------------------------------------------------*#
-using "streams";
+
 
 Model pump