source: trunk/sample/water_steam/sample_power_plant.mso @ 488

Last change on this file since 488 was 315, checked in by Argimiro Resende Secchi, 16 years ago

Water and steam models.
File size: 6.1 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 of a Power Plant.
17*--------------------------------------------------------------------
18* Author: Argimiro R. Secchi
19* $Id: sample_power_plant.mso 195 2007-03-07 20:30:12Z arge $
20*-------------------------------------------------------------------*#
21
22using "water_steam/power_plant";
23
24# Usina termoeletrica
25# Note: model w/o using mixers and splitters, for simplicity.
26FlowSheet USINA
27        PARAMETERS
28    propterm as Plugin(Brief="Steam tables", Type="water", File="propterm");
29
30        DEVICES
31        TAP1    as Turbina_sangra;
32        TAP2    as Turbina;
33        TMP1    as Turbina_sangra;
34        TMP2    as Turbina_sangra;
35        TMP3    as Turbina_sangra;
36        TMP4    as Turbina;
37        TBP1    as Turbina_sangra;
38        TBP2    as Turbina;
39        COND    as Condensador_2alim;
40        TANQUE  as Tanque;
41        RAP1    as Condensador_2alim;
42        RAP2    as Condensador;
43        RBP1    as Condensador_2alim;
44        RBP2    as Condensador;
45        BBP             as Bomba;
46        BAP             as Bomba;
47        GV              as Gerador_Vapor;
48        GE              as Gerador_Eletrico;
49        TR_AP2  as Trocador;
50        TR_AP1  as Trocador;
51        TR_BP1  as Trocador;
52        TR_BP2  as Trocador;
53        TORRE   as Torre;
54
55        CONNECTIONS
56        GV.Fout_sa  to  TAP1.Fin;
57        GV.Fout_ra  to  TMP1.Fin;
58        TAP1.Fout   to  TAP2.Fin;
59        TAP1.Fouts  to  RAP2.Fin;
60        TAP2.Fout   to  GV.Fin_ra;
61        TMP1.Fout   to  TMP2.Fin;
62        TMP1.Fouts  to  RAP1.Fin1;
63        TMP2.Fout   to  TMP3.Fin;
64        TMP2.Fouts  to  TANQUE.Fin1;
65        TMP3.Fout   to  TMP4.Fin;
66        TMP3.Fouts  to  RBP2.Fin;
67        TMP4.Fout   to  TBP1.Fin;
68        TBP1.Fout   to  TBP2.Fin;
69        TBP1.Fouts  to  RBP1.Fin1;
70        TBP2.Fout   to  COND.Fin1;
71        COND.Fout   to  BBP.Fin;
72        BBP.Fout    to  TR_BP1.Fin;
73        TR_BP1.Fout to  TR_BP2.Fin;
74        TR_BP2.Fout to  TANQUE.Fin3;
75        TANQUE.Fout to  BAP.Fin;
76        BAP.Fout    to  TR_AP1.Fin;
77        TR_AP1.Fout to  TR_AP2.Fin;
78        TR_AP2.Fout to  GV.Fin_a;
79        RBP2.Fout       to  RBP1.Fin2;
80        RBP1.Fout   to  COND.Fin2;
81        RAP2.Fout       to  RAP1.Fin2;
82        RAP1.Fout   to  TANQUE.Fin2;
83
84        VARIABLES
85        EF_CIC as Eficiencia;
86        POT_TAP as Potencia;
87        POT_TMP as Potencia;
88        POT_TBP as Potencia;
89        POT_TURB as Potencia;
90
91        EQUATIONS
92        "Trocadores de Calor"
93        COND.Q_COND = TORRE.Q;
94        RBP1.Q_COND = TR_BP1.Q;
95        RBP2.Q_COND = TR_BP2.Q;
96        RAP1.Q_COND = TR_AP1.Q;
97        RAP2.Q_COND = TR_AP2.Q;
98
99        "Eficiencias das Turbinas"
100        TAP2.EF_T = TAP1.EF_T;
101        TMP1.EF_T = TAP2.EF_T;
102        TMP2.EF_T = TMP1.EF_T;
103        TMP3.EF_T = TMP2.EF_T;
104        TMP4.EF_T = TMP3.EF_T;
105        TBP1.EF_T = TMP4.EF_T;
106        TBP2.EF_T = TBP1.EF_T;
107
108        "Potencias das Turbinas"
109        POT_TAP = TAP1.POT_TURB + TAP2.POT_TURB;
110        POT_TMP = TMP1.POT_TURB + TMP2.POT_TURB + TMP3.POT_TURB + TMP4.POT_TURB;
111        POT_TBP = TBP1.POT_TURB + TBP2.POT_TURB;
112        POT_TURB = POT_TAP + POT_TMP + POT_TBP;
113
114        "Potencia do Gerador Eletrico"
115        GE.POT_GE = GE.EF_GE * POT_TURB;
116
117        "Eficiencia do Ciclo"
118        EF_CIC * GV.Q_GV = GE.POT_GE - BBP.POT_BMB - BAP.POT_BMB;
119
120        SET
121        BBP.v_esp = 1.01e-3*'m^3/kg';
122        BAP.v_esp = 1.01e-3*'m^3/kg';
123        GE.EF_GE = 0.96;
124        TORRE.cpa = 4.186e-3*'MJ/kg/K';
125
126        SPECIFY
127#       EF_CIC = 0.27;
128
129        GV.Fin_a.F = (357/3.6)*'kg/s';
130#       GV.Fin_a.P = 134 * 0.101325 * 'MPa';
131#       GV.Fin_a.T = (197 + 273.15) * 'K';
132
133        GV.Feco.P = 120 * 'atm';
134        GV.Feco.T = (314 + 273.15)*'K';
135        GV.Fvap.P = 115 * 0.101325*'MPa';
136        GV.Fvap.T = (400 + 273.15)*'K';
137#       GV.Fout_sa.F = 357/3.6 * 'kg/s';
138        GV.Fout_sa.P = 110 * 0.101325*'MPa';
139        GV.Fout_sa.T = (480 + 273.15)*'K';
140#       GV.Fout_ra.F = 344/3.6 * 'kg/s';
141        GV.Fout_ra.P = 28 * 0.101325*'MPa';
142        GV.Fout_ra.T = (480 + 273.15)*'K';
143#       GV.EF_GV = 0.8;
144        GV.Q_GV = 360 * 'MJ/s';
145
146#       POT_TURB = 100 * 'MW';
147        TAP1.EF_T = 0.8;
148        TAP1.Fout.P = 35 * 0.101325*'MPa';
149        TAP1.y = 0.035;
150        TAP2.Fout.P = 29 * 0.101325*'MPa';
151#       TAP2.Fout.T = (345 + 273.15) * 'K';
152#       TAP2.EF_T = 0.8;
153        TMP1.Fout.P = 26 * 0.101325*'MPa';
154        TMP1.y = 0.035;
155        TMP2.Fout.P = 23 * 0.101325*'MPa';
156        TMP2.y = 0.036;
157        TMP3.Fout.P = 20 * 0.101325*'MPa';
158        TMP3.y = 0.037;
159        TMP4.Fout.P = 15 * 0.101325*'MPa';
160        TBP1.Fout.P = 8 * 0.101325*'MPa';
161        TBP1.y = 0.03;
162        TBP2.Fout.P = 0.2 * 0.101325*'MPa';
163
164        TANQUE.Fout.P = 20 * 0.101325*'MPa';
165
166#       COND.Fout.T = (34 + 273.15)* 'K';
167#       COND.Q_COND = 125 * 'MJ/s';
168        COND.G_S = 0*'K';
169
170        BBP.EF_B = 0.95;
171#       BBP.POT_BMB = 2 * 0.105 * 'MW'; # duas bombas
172        BBP.Fout.P = 10 * 0.101325*'MPa';
173        BAP.EF_B = 0.95;
174#       BAP.POT_BMB = 2 * 6.84 * 'MW'; # duas bombas
175        BAP.Fout.P = 130 * 0.101325*'MPa';
176
177        RAP1.G_S = 0*'K';
178        RAP2.G_S = 0*'K';
179        RBP1.G_S = 0*'K';
180        RBP2.G_S = 0*'K';
181
182#       TR_AP1.DP = 5 * 0.101325 * 'MPa';
183#       TR_AP2.DP = 5 * 0.101325 * 'MPa';
184#       TR_BP1.DP = 2 * 0.101325 * 'MPa';
185#       TR_BP2.DP = 2 * 0.101325 * 'MPa';
186
187        TR_AP1.Fout.P = 125 * 0.101325*'MPa';
188        TR_AP2.Fout.P = 120 * 0.101325*'MPa';
189#       TR_AP2.Fout.T = (197 + 273.15) * 'K';
190        TR_BP1.Fout.P = 8 * 0.101325*'MPa';
191        TR_BP2.Fout.P = 6 * 0.101325*'MPa';
192
193        TORRE.Th = (35 + 273.15)*'K';
194#       TORRE.Tc = (26 + 273.15) * 'K';
195        TORRE.Tar_c = (20 + 273.15)*'K';
196        TORRE.DTar = 10*'K';
197        TORRE.F = 7000*'kg/s';
198
199        OPTIONS
200#       GuessFile = "USINA";
201        Dynamic = false;
202end
203
204# Somente a caldeira
205FlowSheet Caldeira
206        DEVICES
207        GV              as Gerador_Vapor;
208        FEED    as Corrente;
209        F_RA    as Corrente;
210
211        PARAMETERS
212    propterm as Plugin(Brief="Steam tables", Type="water", File="propterm");
213
214        CONNECTIONS
215        FEED to GV.Fin_a;
216        F_RA to GV.Fin_ra;
217
218        EQUATIONS
219        [FEED.S,FEED.H] = propterm.propPTl(FEED.P,FEED.T);
220        GV.Feco.F = GV.Fin_a.F;
221        GV.Fin_ra.T = GV.Fout_ra.T-180*'K';
222        GV.Fin_ra.P = GV.Fout_ra.P+1*'MPa';
223        [GV.Fin_ra.S,GV.Fin_ra.H] = propterm.propPTv(GV.Fin_ra.P,GV.Fin_ra.T);
224
225        SPECIFY
226        FEED.F = 357/3.6 * 'kg/s';
227        FEED.P = 134 * 0.101325 * 'MPa';
228        FEED.T = (197 + 273.15) * 'K';
229        GV.Feco.P = 130 * 0.101325 * 'MPa';
230        GV.Feco.T = (314 + 273.15) * 'K';
231        GV.Fvap.P = 126 * 0.101325 * 'MPa';
232        GV.Fvap.T = (400 + 273.15) * 'K';
233        GV.Fout_sa.P = 121 * 0.101325 * 'MPa';
234        GV.Fout_sa.T = (525 + 273.15) * 'K';
235        GV.Fin_ra.F = 344/3.6 * 'kg/s';
236        GV.Fout_ra.P = 28 * 0.101325 * 'MPa';
237        GV.Fout_ra.T = (526 + 273.15) * 'K';
238        GV.EF_GV = 0.8;
239#       GV.Q_GV = 230 * 'MJ/s';
240
241        OPTIONS
242        Dynamic = false;
243end
Note: See TracBrowser for help on using the repository browser.