diff options
Diffstat (limited to 'Simulator/Simulator/UnitOperations/CentrifugalPump.mo')
-rw-r--r-- | Simulator/Simulator/UnitOperations/CentrifugalPump.mo | 76 |
1 files changed, 0 insertions, 76 deletions
diff --git a/Simulator/Simulator/UnitOperations/CentrifugalPump.mo b/Simulator/Simulator/UnitOperations/CentrifugalPump.mo deleted file mode 100644 index 6b5934b..0000000 --- a/Simulator/Simulator/UnitOperations/CentrifugalPump.mo +++ /dev/null @@ -1,76 +0,0 @@ -within Simulator.UnitOperations; - -model CentrifugalPump "Model of a centrifugal pump to provide energy to liquid stream in form of pressure" - //=========================================================================== - //Header files and Parameters - extends Simulator.Files.Icons.CentrifugalPump; - parameter Simulator.Files.ChemsepDatabase.GeneralProperties C[Nc]; - parameter Integer Nc = 2 "Number of components"; - parameter Real Eff(unit = "-") "Efficiency"; - //=========================================================================== - //Model Variables - Real Pin(unit = "Pa", min = 0, start = Pg) "Inlet stream pressure"; - Real Tin(unit = "K", min = 0, start = Tg) "Inlet stream temperature"; - Real Hin(unit = "kJ/kmol",start=Htotg) "Inlet stream molar enthalpy"; - Real Fin(unit = "mol/s", min = 0, start = Fg) "Inlet stream molar flow"; - Real xin_c[Nc](each unit = "-", each min = 0, each max = 1, start=xg) "Inlet stream components molar fraction"; - Real Tdel(unit = "K") "Temperature increase"; - Real Pdel(unit = "K") "Pressure increase"; - Real Q(unit = "W") "Power required"; - Real rho_c[Nc](each unit = "kmol/m3", each min = 0) "Component molar density"; - Real rho(unit = "kmol/m3", min = 0) "Density"; - Real Pvap(unit = "Pa", min = 0, start = Pg) "Vapor pressure of mixture at Outlet temperature"; - Real NPSH(unit = "m") "Net Positive Suction Head"; - Real Pout(unit = "Pa", min = 0, start = Pg) "Outlet stream pressure"; - Real Tout(unit = "K", min = 0, start = Tg) "Outlet stream temperature"; - Real Hout(unit = "kJ/kmol",start=Htotg) "Outlet stream molar enthalpy"; - Real Fout(unit = "mol/s", min = 0, start = Fg) "Outlet stream molar flow"; - Real xout_c[Nc](each unit = "-", each min = 0, each max = 1, start=xg) "Outlet stream molar fraction"; - //============================================================================ - //Instantiation of Connectors - Simulator.Files.Interfaces.matConn In(Nc = Nc) annotation( - Placement(visible = true, transformation(origin = {-100, -2}, extent = {{-10, -10}, {10, 10}}, rotation = 0), iconTransformation(origin = {-100, 16}, extent = {{-10, -10}, {10, 10}}, rotation = 0))); - Simulator.Files.Interfaces.matConn Out(Nc = Nc) annotation( - Placement(visible = true, transformation(origin = {102, 0}, extent = {{-10, -10}, {10, 10}}, rotation = 0), iconTransformation(origin = {100, 100}, extent = {{-10, -10}, {10, 10}}, rotation = 0))); - Simulator.Files.Interfaces.enConn En annotation( - Placement(visible = true, transformation(origin = {2, -100}, extent = {{-10, -10}, {10, 10}}, rotation = 0), iconTransformation(origin = {0, -70}, extent = {{-10, -10}, {10, 10}}, rotation = 0))); - - extends GuessModels.InitialGuess; -equation -//============================================================================ -//Connector equation - In.P = Pin; - In.T = Tin; - In.F = Fin; - In.H = Hin; - In.x_pc[1, :] = xin_c[:]; - Out.P = Pout; - Out.T = Tout; - Out.F = Fout; - Out.H = Hout; - Out.x_pc[1, :] = xout_c[:]; - En.Q = Q; -//============================================================================= -//Pump equations - Fin = Fout; - xin_c = xout_c; - Pin + Pdel = Pout; - Tin + Tdel = Tout; -//============================================================================= -//Calculation of Density - for i in 1:Nc loop - rho_c[i] = Simulator.Files.ThermodynamicFunctions.Dens(C[i].LiqDen, C[i].Tc, Tin, Pin); - end for; - rho = 1 / sum(xin_c ./ rho_c); -//============================================================================== -//Energy Balance and NPSH Calculation - Hout = Hin + Pdel / rho; - Q = Fin * (Hout - Hin) / Eff; - NPSH = (Pin - Pvap) / rho; -//=============================================================================== -//Vapor Pressure of mixture at Outlet Temperature - Pvap = sum(xin_c .* exp(C[:].VP[2] + C[:].VP[3] / Tout + C[:].VP[4] * log(Tout) + C[:].VP[5] .* Tout .^ C[:].VP[6])); - annotation( - Documentation(info = "<html><head></head><body><div>Centrifugal pump is generally used to provide energy to a liquid material stream. The energy supplied is in form of pressure.</div><div><br></div>To simulate a centrifugal pump, Efficiency of the pump should be provided as calculation parameter. Additional one of the following input variables must be defined:<div><ol><li>Outlet Pressure</li><li>Pressure Increase</li><li>Power Required</li></ol><div><br></div></div><div>For example on simulating a centrifual pum, go to <i><b>Examples</b></i> >> <i><b>Pump</b></i></div></body></html>")); - -end CentrifugalPump; |