[315] | 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 | 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; |
---|
| 202 | end |
---|
| 203 | |
---|
| 204 | # Somente a caldeira |
---|
| 205 | FlowSheet 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; |
---|
| 243 | end |
---|