source: mso/eml/pressure_changers/pump.mso @ 57

Last change on this file since 57 was 57, checked in by Paula Bettio Staudt, 15 years ago

Added new pressure_changers models

  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 4.7 KB
Line 
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: pump.mso 57 2006-11-16 16:46:48Z paula $                                         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
118#*-------------------------------------------------------------------
119* Model of a pump
120*--------------------------------------------------------------------
121*
122*       Streams:
123*               * an  inlet stream
124*               * an  outlet stream
125*
126*       Specify:
127*               * the inlet stream
128*               * the pump press delta
129*
130*----------------------------------------------------------------------
131* Author: Paula B. Staudt
132* $Id: pump.mso 57 2006-11-16 16:46:48Z paula $
133*--------------------------------------------------------------------*#
134
135
136Model pump
137        PARAMETERS
138ext PP as CalcObject;
139ext NComp as Integer;
140       
141        VARIABLES
142in      Inlet as stream;
143out     Outlet as stream_therm;
144       
145        dP as press_delta (Brief="Pump head");
146       
147        EQUATIONS
148        "Molar Balance"
149        Inlet.F = Outlet.F;
150        Inlet.z = Outlet.z;
151       
152        "Pump head"
153        Outlet.P = Inlet.P + dP;
154       
155        "Temperature"
156        Outlet.T = Inlet.T;
157
158        "Vapourisation Fraction"
159        Outlet.v = Inlet.v;
160end
Note: See TracBrowser for help on using the repository browser.