Ignore:
Timestamp:
Mar 3, 2009, 11:49:26 PM (13 years ago)
Author:
gerson bicca
Message:

added a work stream into pump model

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
     5using "streams";
    256using "pressure_changers/pump";
    267
    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",
     8FlowSheet Sample_Pump_Diagram
     9        PARAMETERS
     10        PP as Plugin(Brief="Physical Properties",
    3811                Type="PP",
    3912                Components = ["water"],
     
    4114                VapourModel = "PR"
    4215        );
    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;
    4522
    4623        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
    5525        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';
    6530
    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
    7532
    76 DEVICES
    77         p1 as pump;
    78         s1 as source;
     33        GUESS
     34
     35        DEVICES
     36        discharge as simple_sink;
    7937
    8038        SET
    81         NComp = PP.NumberOfComponents;
    82        
    83         s1.ValidPhases = "Vapour-Liquid";
    84         s1.CompositionBasis = "Molar";
    85        
    86 
    87        
    88         CONNECTIONS
    89         s1.Outlet to p1.Inlet;
    9039
    9140        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;
    9978
    10079        OPTIONS
    10180        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        );
    10294end
Note: See TracChangeset for help on using the changeset viewer.