Changeset 57 for mso/eml/pressure_changers/pump.mso
- Timestamp:
- Nov 16, 2006, 2:46:48 PM (16 years ago)
- 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 *-------------------------------------------------------------------------*# 16 using "streams"; 17 using "pressure_changers/flux_machine_basic"; 18 19 Model centrifugal_pump as flux_machine_basic 20 21 PARAMETERS 22 ext NComp as Integer (Brief = "Number of chemical components", Lower = 1); 23 ext 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 115 end 116 117 1 118 #*------------------------------------------------------------------- 2 119 * Model of a pump … … 15 132 * $Id$ 16 133 *--------------------------------------------------------------------*# 17 using "streams"; 134 18 135 19 136 Model pump
Note: See TracChangeset
for help on using the changeset viewer.