Changeset 742 for branches/gui/sample/pressure_changers/sample_pump.mso
- Timestamp:
- Mar 3, 2009, 11:49:26 PM (14 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/gui/sample/pressure_changers/sample_pump.mso
r740 r742 1 #*------------------------------------------------------------------- 2 * EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC. 3 * 4 * This LIBRARY is free software; you can distribute it and/or modify 5 * it under the therms of the ALSOC FREE LICENSE as available at 6 * http://www.enq.ufrgs.br/alsoc. 7 * 8 * EMSO Copyright (C) 2004 - 2007 ALSOC, original code 9 * from http://www.rps.eng.br Copyright (C) 2002-2004. 10 * All rights reserved. 11 * 12 * EMSO is distributed under the therms of the ALSOC LICENSE as 13 * available at http://www.enq.ufrgs.br/alsoc. 14 * 15 *-------------------------------------------------------------------- 16 * Sample file for model centrifugal_pump and pump (simplified) 17 *-------------------------------------------------------------------- 18 * 19 * This sample file needs VRTherm DEMO (www.vrtech.com.br) to run. 20 * 21 *---------------------------------------------------------------------- 22 * Author: Estefane S. Horn, Paula B. Staudt 23 * $Id$ 24 *--------------------------------------------------------------------*# 1 #*---------------------------------------------- 2 * FlowSheet generated automaticaly by EMSO-GUI 3 *----------------------------------------------*# 4 5 using "streams"; 25 6 using "pressure_changers/pump"; 26 7 27 FlowSheet Pump_P1 28 29 DEVICES 30 S1 as source; 31 P1 as centrifugal_pump; 32 33 CONNECTIONS 34 S1.Outlet to P1.Inlet; 35 36 PARAMETERS 37 PP as Plugin (Brief="External Physical Properties", 8 FlowSheet Sample_Pump_Diagram 9 PARAMETERS 10 PP as Plugin(Brief="Physical Properties", 38 11 Type="PP", 39 12 Components = ["water"], … … 41 14 VapourModel = "PR" 42 15 ); 43 Mw as positive (Brief = "Molar Mass", Unit = 'kg/kmol'); 44 NComp as Integer (Default= 1); 16 NComp as Integer; 17 SET 18 NComp = PP.NumberOfComponents; 19 20 DEVICES 21 suction as simple_source; 45 22 46 23 SET 47 NComp = PP.NumberOfComponents; 48 P1.MechanicalEff = 0.95; 49 P1.PumpEfficiency = 0.72; 50 P1.SuctionArea = 0.001*'m^2'; 51 52 S1.ValidPhases = "Vapour-Liquid"; 53 S1.CompositionBasis = "Molar"; 54 24 55 25 SPECIFY 56 S1.F = 1000 * 'kmol/h'; 57 S1.P = 10 * 'atm'; 58 S1.T = 298 * 'K' ; 59 S1.Composition = [1]; 60 61 P1.StaticHead = 10 * 'm'; 62 P1.Outlet.P = 18 * 'atm'; 63 64 end 26 suction.MolarComposition = 1; 27 suction.F = 1 * 'kmol/h'; 28 suction.T = 298 * 'K'; 29 suction.P = 2 * 'atm'; 65 30 66 FlowSheet pump_Test 67 PARAMETERS 68 PP as Plugin (Brief="External Physical Properties", 69 Type="PP", 70 Components = ["water", "benzene"], 71 LiquidModel = "PR", 72 VapourModel = "PR" 73 ); 74 NComp as Integer; 31 INITIAL 75 32 76 DEVICES 77 p1 as pump; 78 s1 as source; 33 GUESS 34 35 DEVICES 36 discharge as simple_sink; 79 37 80 38 SET 81 NComp = PP.NumberOfComponents;82 83 s1.ValidPhases = "Vapour-Liquid";84 s1.CompositionBasis = "Molar";85 86 87 88 CONNECTIONS89 s1.Outlet to p1.Inlet;90 39 91 40 SPECIFY 92 s1.F = 826 * 'kmol/h'; 93 s1.P = 549 * 'kPa'; 94 s1.T = 313 * 'K'; 95 s1.Composition = 1/NComp; 96 #s1.Outlet.v = 0.68; 97 98 p1.dP = 90 * 'kPa'; 41 42 INITIAL 43 44 GUESS 45 46 DEVICES 47 Pump_Operation as centrifugal_pump; 48 49 SET 50 Pump_Operation.PumpEfficiency = 0.10; 51 Pump_Operation.MechanicalEff = 0.95; 52 Pump_Operation.NPSH_Options = "Include Kinetic Head"; 53 Pump_Operation.SuctionArea = 0.01 * 'm^2'; 54 55 SPECIFY 56 Pump_Operation.Pratio = 7; 57 Pump_Operation.StaticHead = 12 * 'm'; 58 59 INITIAL 60 61 GUESS 62 63 DEVICES 64 work_inlet as work_source; 65 66 SET 67 68 SPECIFY 69 70 INITIAL 71 72 GUESS 73 74 CONNECTIONS 75 suction.Outlet to Pump_Operation.Inlet; 76 Pump_Operation.Outlet to discharge.Inlet; 77 work_inlet.Work to Pump_Operation.WorkIn; 99 78 100 79 OPTIONS 101 80 Dynamic = false; 81 Integration = "original"; 82 NLASolver( 83 File = "sundials", 84 RelativeAccuracy = 1e-3, 85 AbsoluteAccuracy = 1e-6, 86 MaxIterations = 100 87 ); 88 DAESolver( 89 File = "sundials", 90 RelativeAccuracy = 1e-3, 91 AbsoluteAccuracy = 1e-6, 92 EventAccuracy = 1e-2 93 ); 102 94 end
Note: See TracChangeset
for help on using the changeset viewer.