source: branches/new_gui/my_folders/water_steam/sample_rankine.mso @ 897

Last change on this file since 897 was 897, checked in by gerson bicca, 13 years ago

added more folders

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