source: branches/gui/sample/water_steam/sample_rankine.mso @ 963

Last change on this file since 963 was 963, checked in by Argimiro Resende Secchi, 6 years ago

Update water_steam

File size: 4.2 KB
Line 
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 Rankine cicle.
17*--------------------------------------------------------------------
18* Author: Argimiro R. Secchi
19* $Id: sample_rankine.mso 195 2007-03-07 20:30:12Z arge $
20*-------------------------------------------------------------------*#
21
22using "water_steam/power_plant";
23
24# Simple Rankine cycle
25FlowSheet Rankine
26
27    PARAMETERS
28    PP2 as Plugin(Brief="Steam tables", Type="water");
29
30    DEVICES
31        Turb    as Turbina;
32        Cond    as Condensador;
33        Bomba1  as Bomba;
34        GV              as Gerador_Vapor_Simples;
35        GE              as Gerador_Eletrico_Simples;
36
37        CONNECTIONS
38        GV.Fout      to  Turb.Fin;
39        Turb.Fout    to  Cond.Fin;
40        Cond.Fout    to  Bomba1.Fin;
41        Bomba1.Fout  to  GV.Fin;
42
43        VARIABLES
44        EF_CIC as Eficiencia;
45
46        EQUATIONS
47        "Potencia do Gerador Eletrico"
48        GE.POT_GE = GE.EF_GE * Turb.POT_TURB;
49
50        "Eficiencia do Ciclo"
51        EF_CIC * GV.Q_GV = GE.POT_GE - Bomba1.POT_BMB;
52
53        SET
54        GE.EF_GE = 0.97;
55
56        SPECIFY
57#       EF_CIC = 0.275;
58
59        GV.Fout.P = 11.3 *'MPa';
60        GV.Fout.T = (530 + 273.15)*'K';
61#       GV.Fout.F = 100 * 'kg/s';
62        GV.EF_GV = 0.8;
63#       GV.Q_GV = 230 * 'MJ/s';
64
65        Turb.POT_TURB = 63*'MW';
66        Turb.EF_T = 0.8;
67        Turb.Fout.P = 0.007*'MPa';
68
69#       Cond.Fout.T = (34 + 273.15) * 'K';
70#       Cond.Q_COND = 125 * 'MJ/s';
71        Cond.G_S = 5.0*'K';
72
73#       Bomba1.Fout.F = 60 * 'kg/s';
74        Bomba1.EF_B = 0.7;
75#       Bomba1.POT_BMB = 0.9 * 'MW';
76
77        OPTIONS
78        Dynamic = false;
79end
80
81# Simple Rankine cycle with extraction
82FlowSheet Rankine_s
83        DEVICES
84        TAP1    as Turbina_sangra;
85        TAP2    as Turbina;
86        COND    as Condensador_2alim;
87        BBP             as Bomba;
88        RAP2    as Condensador;
89        TR_AP2  as Trocador;
90        GV              as Gerador_Vapor;
91        GE              as Gerador_Eletrico_Simples;
92        F_RA    as Corrente;
93
94    PARAMETERS
95    PP2 as Plugin(Brief="Steam tables", Type="water");
96
97        CONNECTIONS
98        F_RA            to      GV.Fin_ra;
99        GV.Fout_sa  to  TAP1.Fin;
100        TAP1.Fout   to  TAP2.Fin;
101        TAP2.Fout   to  COND.Fin1;
102        TAP1.Fouts  to  RAP2.Fin;
103        RAP2.Fout   to  COND.Fin2;
104        COND.Fout   to  BBP.Fin;
105        BBP.Fout    to  TR_AP2.Fin;
106        TR_AP2.Fout to  GV.Fin_a;
107
108        VARIABLES
109        EF_CIC as Eficiencia;
110        POT_TURB as Potencia;
111
112        EQUATIONS
113        F_RA.T = GV.Fout_ra.T-180*'K';
114        F_RA.P = GV.Fout_ra.P+1*'MPa';
115        [F_RA.S,F_RA.H] = PP2.propPTv(F_RA.P,F_RA.T);
116
117#       GV.Fin_a.P = GV.Fout_sa.P;
118
119        "Trocadores de Calor"
120        RAP2.Q_COND = TR_AP2.Q;
121
122        "Eficiencias das Turbinas"
123        TAP1.EF_T = TAP2.EF_T;
124
125        "Potencias das Turbinas"
126        POT_TURB = TAP1.POT_TURB + TAP2.POT_TURB;
127
128        "Potencia do Gerador Eletrico"
129        GE.POT_GE = GE.EF_GE * POT_TURB;
130
131        "Eficiencia do Ciclo"
132        EF_CIC * GV.Q_GV = GE.POT_GE - BBP.POT_BMB;
133
134        SET
135        GE.EF_GE = 0.97;
136
137        SPECIFY
138#       EF_CIC = 0.275;
139
140        GV.Fin_a.F = (357/3.6)*'kg/s';
141#       GV.Fin_a.P = 134 * 0.101325*'MPa';
142#       GV.Fin_a.T = (197 + 273.15)*'K';
143        GV.Feco.P = 130 * 0.101325*'MPa';
144        GV.Feco.T = (314 + 273.15)*'K';
145        GV.Fvap.P = 126 * 0.101325*'MPa';
146        GV.Fvap.T = (400 + 273.15)*'K';
147        GV.Fout_sa.P = 121 * 0.101325*'MPa';
148        GV.Fout_sa.T = (525 + 273.15)*'K';
149        GV.Fin_ra.F = (344/3.6)*'kg/s';
150        GV.Fout_ra.P = 28 * 0.101325*'MPa';
151        GV.Fout_ra.T = (526 + 273.15)*'K';
152        GV.EF_GV = 0.8;
153#       GV.Q_GV = 400*'MJ/s';
154
155#       TAP1.POT_TURB = 43*'MW';
156#       TAP1.EF_T = 0.8;
157        TAP1.Fout.P = 2*'MPa';
158        TAP1.y = 0.035;
159#       TAP2.POT_TURB = 23*'MW';
160        TAP2.EF_T = 0.8;
161        TAP2.Fout.P = 0.2*'MPa';
162#       TAP2.Fout.T = 315*'K';
163
164#       COND.Fout.T = (34 + 273.15)*'K';
165#       COND.Q_COND = 125*'MJ/s';
166        COND.G_S = 5.0*'K';
167
168        BBP.EF_B = 0.95;
169#       BBP.POT_BMB = 2 * 0.8*'MW'; #duas bombas
170        BBP.Fout.P = 146 * 0.101325*'MPa';
171
172        RAP2.G_S = 1*'K';
173
174#       TR_AP2.Fout.P = 134 * 0.101325*'MPa';
175        TR_AP2.DP = 20 * 0.101325*'MPa';
176
177        OPTIONS
178#       GuessFile = "Rankine_s";
179        Dynamic = false;
180end
Note: See TracBrowser for help on using the repository browser.