| 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 | |
|---|
| 22 | using "water_steam/power_plant"; |
|---|
| 23 | |
|---|
| 24 | # Usina termoeletrica |
|---|
| 25 | # Note: model w/o using mixers and splitters, for simplicity. |
|---|
| 26 | FlowSheet USINA |
|---|
| 27 | PARAMETERS |
|---|
| 28 | PP2 as Plugin(Brief="Steam tables", Type="water"); |
|---|
| 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_Simples; |
|---|
| 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 | COND.Q_COND to TORRE.Qin; |
|---|
| 84 | |
|---|
| 85 | VARIABLES |
|---|
| 86 | EF_CIC as Eficiencia; |
|---|
| 87 | POT_TAP as Potencia; |
|---|
| 88 | POT_TMP as Potencia; |
|---|
| 89 | POT_TBP as Potencia; |
|---|
| 90 | POT_TURB as Potencia; |
|---|
| 91 | |
|---|
| 92 | EQUATIONS |
|---|
| 93 | "Trocadores de Calor" |
|---|
| 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 | GE.EF_GE = 0.96; |
|---|
| 122 | TORRE.cpa = 4.186e-3*'MJ/kg/K'; |
|---|
| 123 | |
|---|
| 124 | SPECIFY |
|---|
| 125 | # EF_CIC = 0.27; |
|---|
| 126 | |
|---|
| 127 | GV.Fin_a.F = (357/3.6)*'kg/s'; |
|---|
| 128 | # GV.Fin_a.P = 134 * 0.101325 * 'MPa'; |
|---|
| 129 | # GV.Fin_a.T = (197 + 273.15) * 'K'; |
|---|
| 130 | |
|---|
| 131 | GV.Feco.P = 120 * 'atm'; |
|---|
| 132 | GV.Feco.T = (314 + 273.15)*'K'; |
|---|
| 133 | GV.Fvap.P = 115 * 0.101325*'MPa'; |
|---|
| 134 | GV.Fvap.T = (400 + 273.15)*'K'; |
|---|
| 135 | # GV.Fout_sa.F = 357/3.6 * 'kg/s'; |
|---|
| 136 | GV.Fout_sa.P = 110 * 0.101325*'MPa'; |
|---|
| 137 | GV.Fout_sa.T = (480 + 273.15)*'K'; |
|---|
| 138 | # GV.Fout_ra.F = 344/3.6 * 'kg/s'; |
|---|
| 139 | GV.Fout_ra.P = 28 * 0.101325*'MPa'; |
|---|
| 140 | GV.Fout_ra.T = (480 + 273.15)*'K'; |
|---|
| 141 | # GV.EF_GV = 0.8; |
|---|
| 142 | GV.Q_GV = 360 * 'MJ/s'; |
|---|
| 143 | |
|---|
| 144 | # POT_TURB = 100 * 'MW'; |
|---|
| 145 | TAP1.EF_T = 0.8; |
|---|
| 146 | TAP1.Fout.P = 35 * 0.101325*'MPa'; |
|---|
| 147 | TAP1.y = 0.035; |
|---|
| 148 | TAP2.Fout.P = 29 * 0.101325*'MPa'; |
|---|
| 149 | # TAP2.Fout.T = (345 + 273.15) * 'K'; |
|---|
| 150 | # TAP2.EF_T = 0.8; |
|---|
| 151 | TMP1.Fout.P = 26 * 0.101325*'MPa'; |
|---|
| 152 | TMP1.y = 0.035; |
|---|
| 153 | TMP2.Fout.P = 23 * 0.101325*'MPa'; |
|---|
| 154 | TMP2.y = 0.036; |
|---|
| 155 | TMP3.Fout.P = 20 * 0.101325*'MPa'; |
|---|
| 156 | TMP3.y = 0.037; |
|---|
| 157 | TMP4.Fout.P = 15 * 0.101325*'MPa'; |
|---|
| 158 | TBP1.Fout.P = 8 * 0.101325*'MPa'; |
|---|
| 159 | TBP1.y = 0.03; |
|---|
| 160 | TBP2.Fout.P = 0.2 * 0.101325*'MPa'; |
|---|
| 161 | |
|---|
| 162 | TANQUE.Fout.P = 20 * 0.101325*'MPa'; |
|---|
| 163 | |
|---|
| 164 | # COND.Fout.T = (34 + 273.15)* 'K'; |
|---|
| 165 | # COND.Q_COND = 125 * 'MJ/s'; |
|---|
| 166 | COND.G_S = 0*'K'; |
|---|
| 167 | |
|---|
| 168 | BBP.EF_B = 0.95; |
|---|
| 169 | # BBP.POT_BMB = 2 * 0.105 * 'MW'; # duas bombas |
|---|
| 170 | BBP.Fout.P = 10 * 0.101325*'MPa'; |
|---|
| 171 | BAP.EF_B = 0.95; |
|---|
| 172 | # BAP.POT_BMB = 2 * 6.84 * 'MW'; # duas bombas |
|---|
| 173 | BAP.Fout.P = 130 * 0.101325*'MPa'; |
|---|
| 174 | |
|---|
| 175 | RAP1.G_S = 0*'K'; |
|---|
| 176 | RAP2.G_S = 0*'K'; |
|---|
| 177 | RBP1.G_S = 0*'K'; |
|---|
| 178 | RBP2.G_S = 0*'K'; |
|---|
| 179 | |
|---|
| 180 | # TR_AP1.DP = 5 * 0.101325 * 'MPa'; |
|---|
| 181 | # TR_AP2.DP = 5 * 0.101325 * 'MPa'; |
|---|
| 182 | # TR_BP1.DP = 2 * 0.101325 * 'MPa'; |
|---|
| 183 | # TR_BP2.DP = 2 * 0.101325 * 'MPa'; |
|---|
| 184 | |
|---|
| 185 | TR_AP1.Fout.P = 125 * 0.101325*'MPa'; |
|---|
| 186 | TR_AP2.Fout.P = 120 * 0.101325*'MPa'; |
|---|
| 187 | # TR_AP2.Fout.T = (197 + 273.15) * 'K'; |
|---|
| 188 | TR_BP1.Fout.P = 8 * 0.101325*'MPa'; |
|---|
| 189 | TR_BP2.Fout.P = 6 * 0.101325*'MPa'; |
|---|
| 190 | |
|---|
| 191 | TORRE.Th = (35 + 273.15)*'K'; |
|---|
| 192 | # TORRE.Tc = (26 + 273.15) * 'K'; |
|---|
| 193 | TORRE.Tar_c = (20 + 273.15)*'K'; |
|---|
| 194 | TORRE.DTar = 10*'K'; |
|---|
| 195 | TORRE.F = 7000*'kg/s'; |
|---|
| 196 | |
|---|
| 197 | OPTIONS |
|---|
| 198 | # GuessFile = "USINA"; |
|---|
| 199 | Dynamic = false; |
|---|
| 200 | end |
|---|
| 201 | |
|---|
| 202 | # Somente a caldeira |
|---|
| 203 | FlowSheet Caldeira |
|---|
| 204 | DEVICES |
|---|
| 205 | GV as Gerador_Vapor; |
|---|
| 206 | FEED as Corrente; |
|---|
| 207 | F_RA as Corrente; |
|---|
| 208 | |
|---|
| 209 | PARAMETERS |
|---|
| 210 | PP2 as Plugin(Brief="Steam tables", Type="water"); |
|---|
| 211 | |
|---|
| 212 | CONNECTIONS |
|---|
| 213 | FEED to GV.Fin_a; |
|---|
| 214 | F_RA to GV.Fin_ra; |
|---|
| 215 | |
|---|
| 216 | EQUATIONS |
|---|
| 217 | [FEED.S,FEED.H] = PP2.propPTl(FEED.P,FEED.T); |
|---|
| 218 | GV.Feco.F = GV.Fin_a.F; |
|---|
| 219 | GV.Fin_ra.T = GV.Fout_ra.T-180*'K'; |
|---|
| 220 | GV.Fin_ra.P = GV.Fout_ra.P+1*'MPa'; |
|---|
| 221 | [GV.Fin_ra.S,GV.Fin_ra.H] = PP2.propPTv(GV.Fin_ra.P,GV.Fin_ra.T); |
|---|
| 222 | |
|---|
| 223 | SPECIFY |
|---|
| 224 | FEED.F = 357/3.6 * 'kg/s'; |
|---|
| 225 | FEED.P = 134 * 0.101325 * 'MPa'; |
|---|
| 226 | FEED.T = (197 + 273.15) * 'K'; |
|---|
| 227 | GV.Feco.P = 130 * 0.101325 * 'MPa'; |
|---|
| 228 | GV.Feco.T = (314 + 273.15) * 'K'; |
|---|
| 229 | GV.Fvap.P = 126 * 0.101325 * 'MPa'; |
|---|
| 230 | GV.Fvap.T = (400 + 273.15) * 'K'; |
|---|
| 231 | GV.Fout_sa.P = 121 * 0.101325 * 'MPa'; |
|---|
| 232 | GV.Fout_sa.T = (525 + 273.15) * 'K'; |
|---|
| 233 | GV.Fin_ra.F = 344/3.6 * 'kg/s'; |
|---|
| 234 | GV.Fout_ra.P = 28 * 0.101325 * 'MPa'; |
|---|
| 235 | GV.Fout_ra.T = (526 + 273.15) * 'K'; |
|---|
| 236 | GV.EF_GV = 0.8; |
|---|
| 237 | # GV.Q_GV = 230 * 'MJ/s'; |
|---|
| 238 | |
|---|
| 239 | OPTIONS |
|---|
| 240 | Dynamic = false; |
|---|
| 241 | end |
|---|
