source: mso/eml/pressure_changers/turbine.mso @ 75

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

Updated pressure_changers files header

File size: 3.1 KB
RevLine 
[75]1#*-------------------------------------------------------------------
2* EMSO Model Library (EML) Copyright (C) 2004 - 2007 ALSOC.
[57]3*
[75]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* Model of a Hidraulic Turbine
17*--------------------------------------------------------------------
18*
19*       - Assumptions
20*               * Steady State
21*               * Only Liquid
22*               * Adiabatic
23*               * Isentropic
24*
25*----------------------------------------------------------------------
26* Authors: Andrey Copat, Estefane S. Horn, Marcos L. Alencastro
27* $Id$
28*--------------------------------------------------------------------*#
29
[57]30using "streams";
31using "pressure_changers/flux_machine_basic";
32
33Model Hidraulic_Turbine as flux_machine_basic
34
35        PARAMETERS
36ext NComp       as Integer                      (Brief = "Number of chemical components", Lower = 1);
37ext PP                  as CalcObject           (Brief = "External Physical Properties");
38        Mw(NComp)       as molweight            (Brief = "Molar Weight");
39        Eff     as positive                     (Default = 0.72, Brief = "Pump efficiency");
40        Meff    as positive                     (Default = 1.0, Brief = "Brake efficiency");
41        Beta    as positive                     (Default = 0, Brief = "Volumetric expansivity", Unit = "1/K");
42       
43       
44        VARIABLES
45        Head    as head                         (Brief = "Head Developed");
46        FPower  as power                        (Brief = "Fluid Power");
47        BPower  as power                        (Brief = "Brake Power");
48        EPower  as power                        (Brief = "Eletrical Potency");
49        Pdiff   as press_delta          (Brief = "Pressure Increase");
50        Pratio  as positive                     (Brief = "Pressure Ratio");     
51        Mwm     as molweight            (Brief = "Mixture Molar Weight");
52        rho             as dens_mass            (Brief = "Specific Mass", Unit="kg/m^3");
53        Cp              as cp_mol                       (Brief = "Heat Capacity");
54       
55        SET
56        Mw = PP.MolecularWeight();     
57       
58        EQUATIONS
59        #Mixtures Properties
60        "Calculate Mwm for Inlet Mixture"
61        Mwm = sum(Mw([1:NComp])*Inlet.z([1:NComp]));
62       
63        "Calculate rho using a External Physical Properties Routine"
64        rho = PP.LiquidDensity(Inlet.T,Inlet.P,Inlet.z);
65       
66        "Calculate Outlet Vapour Fraction"
67        Outlet.v = PP.VapourFraction(Outlet.T, Outlet.P, Outlet.z);
68       
69        "Calculate Cp Using a External Physical Properties Routine"
70        Cp = PP.LiquidCp(Inlet.T,Inlet.P,Inlet.z);
71       
72        #Mass and Energy Balance and Turbine 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 Fluid Power"
80        FPower = Pdiff * Inlet.F / (rho/Mwm);
81       
82        "Calculate Brake Power"
83        BPower = FPower * Eff;
84       
85        "Calculate Eletric Power"
86        EPower = BPower * Meff;
87       
88        "Calculate Outlet Temperature"
89        (Outlet.T - Inlet.T)*Cp = (Outlet.h - Inlet.h) -  Pdiff /(rho/Mwm) * (1-Beta*Inlet.T);
90       
91        "Calculate Outlet Enthalpy"
92        Outlet.h - Inlet.h =  Pdiff / (rho/Mwm);
93       
94        "Molar Balance"
95        Outlet.F = Inlet.F;
96       
97        "Calculate Outlet Composition"
98        Outlet.z = Inlet.z;
99
100        "Calculate Head"
101        Head = Pdiff/rho;
102
103end
Note: See TracBrowser for help on using the repository browser.