Ignore:
Timestamp:
Nov 16, 2006, 2:46:48 PM (16 years ago)
Author:
Paula Bettio Staudt
Message:

Added new pressure_changers models

File:
1 edited

Legend:

Unmodified
Added
Removed
  • mso/eml/pressure_changers/pump.mso

    r1 r57  
     1#*------------------------------------------------------------------------
     2* This file is property of the author and cannot be used, copyed
     3* or modified without permission.
     4*
     5* Copyright (C) 2002-2006  the author
     6*-------------------------------------------------------------------------
     7* Authors:      Andrey Copat            Estefane da Silveira Horn
     8* $Id$                                          Marcos Lovato Alencastro       
     9*                                                                                                               Date:   20/02/2006
     10*-------------------------------------------------------------------------
     11* -> Steady State
     12* -> Only Liquid
     13* -> Adiabatic
     14* -> Isentropic
     15*-------------------------------------------------------------------------*#
     16using "streams";
     17using "pressure_changers/flux_machine_basic";
     18
     19Model centrifugal_pump as flux_machine_basic
     20       
     21        PARAMETERS
     22ext NComp       as Integer                      (Brief = "Number of chemical components", Lower = 1);
     23ext PP                  as CalcObject           (Brief = "External Physical Properties");
     24        Mw(NComp)       as molweight            (Brief = "Molar Weight");
     25        Eff             as positive             (Default = 0.72, Brief = "Pump Efficiency");
     26        Meff            as positive             (Default = 0.95, Brief = "Brake Efficiency");
     27        Beta            as positive             (Default = 0, Brief = "Volumetric Expansivity", Unit = "1/K");
     28        g                       as acceleration         (Brief = "Gravity Acceleration", Default = 9.81);
     29        N                       as vel_angular          (Brief = "Rotation", Default = 100);
     30        Lev                     as length                       (Brief = "Loss Friction", Default = 0);
     31       
     32        VARIABLES
     33        rho               as dens_mass          (Brief = "Specific Mass", Unit="kg/m^3");
     34        Cp        as cp_mol                     (Brief = "Heat Capacity", Unit="kJ/kmol/K");
     35        FPower    as power                      (Brief = "Fluid Power", Unit="kW");
     36        BPower    as power                      (Brief = "Brake Power",Unit="kW");
     37        EPower    as power                      (Brief = "Eletrical Potency", Unit="kW");
     38        Pdiff     as press_delta        (Brief = "Pressure Increase", Unit="kPa");
     39        Pratio    as positive           (Brief = "Pressure Ratio");     
     40        Head      as head                       (Brief = "Head Developed", Unit="kJ/kmol");
     41        Head_is   as head                       (Brief = "Isoentripic Head", Unit="kJ/kmol");
     42        Mwm       as molweight          (Brief = "Mixture Molar Weight");
     43        pvm               as pressure           (Brief = "Mixture Vapour Pressure", Unit = "kPa");                     
     44        NPSHa     as length                     (Brief = "Available Net Positive Suction Head");
     45        NS                as positive           (Brief = "Specific Speed", Unit = "(rpm*(gal/min)^0.5)/(m^3/4)");
     46        Q                 as flow_vol           (Brief = "Volumetric Flow Rate");
     47        vm                as vol_mol            (Brief = "Mixture Molar Volume", Unit = "m^3/kmol");
     48       
     49        SET
     50        Mw = PP.MolecularWeight();
     51       
     52        EQUATIONS
     53        #Mixtures Properties
     54        "Calculate Mwm for Inlet Mixture"
     55        Mwm = sum(Mw([1:NComp])*Inlet.z([1:NComp]));
     56
     57        "Calculate Cp Using a External Physical Properties Routine"
     58        Cp = PP.LiquidCp(Inlet.T,Inlet.P,Inlet.z);
     59       
     60        "Calculate rho using a External Physical Properties Routine"
     61        rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
     62
     63        "Calculate Mixture Vapour Pressure"
     64        [pvm] = PP.BubbleP(Inlet.T,Inlet.z);
     65       
     66        "Calculate Outlet Vapour Fraction"
     67        Outlet.v = PP.VapourFraction(Outlet.T, Outlet.P, Outlet.z);
     68       
     69        "Calculate Liquid Molar Volume"
     70        vm = PP.LiquidVolume(Inlet.T,Inlet.P,Inlet.z);
     71       
     72        #Mass and Energy Balance and Pump Equations     
     73        "Calculate Outlet Stream Pressure"
     74        Outlet.P = Inlet.P + Pdiff;
     75       
     76        "Pratio Definition"
     77        Outlet.P = Inlet.P * Pratio;
     78       
     79        "Calculate Isentropic Head"
     80        Head_is = Pdiff * Mwm/rho;
     81
     82        "Calculate Real Head"
     83        Head = Head_is/(Eff*Meff);
     84       
     85        "Calculate Outlet Enthalpy"
     86        Outlet.h - Inlet.h = Head;
     87       
     88        "Calculate Fluid Power"
     89        FPower = Head_is * Inlet.F;
     90
     91        "Calculate Brake Power"
     92        BPower * Eff = FPower;
     93
     94        "Calculate Eletric Power"
     95        BPower = EPower * Meff;
     96       
     97        "Calculate Outlet Temperature"
     98        (Outlet.T - Inlet.T) * Cp = (Outlet.h - Inlet.h) -  Pdiff * Mwm / rho * (1 - Beta * Inlet.T);
     99       
     100        "Molar Balance"
     101        Outlet.F = Inlet.F;
     102
     103        "Calculate Outlet Composition"
     104        Outlet.z = Inlet.z;
     105       
     106        "Calculate Net Positive Suction Head"
     107        NPSHa = - Lev + (Inlet.P - pvm)/(g*rho); #If Inlet.P is the suction pump pressure, Lev is 0.
     108       
     109        "Calculate Volumetric Flow Rate"
     110        Q = Inlet.F*vm;
     111       
     112        "Calculate Specific Speed"
     113        NS = N*(Q^0.5)/(Head^3/4);
     114       
     115end
     116
     117
    1118#*-------------------------------------------------------------------
    2119* Model of a pump
     
    15132* $Id$
    16133*--------------------------------------------------------------------*#
    17 using "streams";
     134
    18135
    19136Model pump
Note: See TracChangeset for help on using the changeset viewer.