Данное выражение может быть заложено в отдельную подпрограмму, которая после небольшого числа циклов производит расчёт теплообменника по входным температурам и расходам сред, давлениям пара на входе и выходе, составу уходящих газов и комплексу kiFi. Пример такой подпрограммы приведён в приложении 5. Теплофизические свойства веществ в программе определяются по полиномам из [2] (для газов) и формуляции IF‑97 (для воды и пара). Аналогичные выражения для испарительного элемента с известной входной температурой газов и элемента с одной фазой, у которого заданы температуры сред на горячем конце, принимают вид соответственно
, .Расход греющего пара на деаэратор по отношению к расходу конденсата можно определить как
Кроме того, для улучшения сходимости можно использовать переходы, отличающиеся от простой подстановки нового значения (при этом последовательность таких переходов должна сходиться к требуемой величине). Например, итерации по недогреву в экономайзере
выполняют переход в зависимости не только от значения, определённого по расчёту элементов, но и от значений, принятого изначально и полученного на предыдущей итерации (см. алгоритм). Циклы с мало зависящими друг от друга переменными для упрощения алгоритма можно объединять; в самом внешнем (нулевого уровня) цикле все критерии сходимости опять проверяются, и если хотя бы один из них не выполнен, расчёт продолжается (циклы остальных уровней реализованы с постусловием).Список переменных, определяемых для данного алгоритма методом последовательных приближений, представлен в таблице 4; сам алгоритм приведён в приложении 2. После 51 цикла нижнего уровня (по t0) алгоритм сошёлся с требуемой точностью, результаты приведены в таблице 5.
Таблица 3. Организация расчёта на сходимость: циклы и переменные
Уровень цикла | Наименование величины | Обозначение | ЕИ | Критерий сходимости | Новое значение | Число шагов |
1 | Удельный расход греющего пара | gгп | – | 17 | ||
Недогрев в экономайзере до tб | δtнед | град. | ||||
Температура в конденсаторе | tк | град. | ||||
2 | Расход свежего пара | G0 | кг/с | 48 | ||
3 | Температура свежего пара | t0 | град. | 51 |
2.3. Проверка расчёта КПД ПГУ по обратному методу
КПД ПГУ при расчёте определялся по прямому методу
Потери энергии в ПГУ включают следующие компоненты:
· теплота уходящих газов, за вычетом теплоты воздуха и топлива,
Qух = 2(GухIух–GвIв–10–Gтопhтоп)=2·(226,47·472,98–222,02·225,63–4,46·1579,7) =
= 79 083 кДж/с;
· теплота, отводимая в конденсаторе, Qк=2G0(hк–hк´)=2·28,57·(2323,49–142,16)=
= 128 093 кДж/с;
· потери в ГТУ в генераторе, на трение в подшипниках и редукторе и наружное охлаждение
· потери в генераторе и на трение в утилизационной установке
ΔNпт = Nг(1/(ηмηг)–1) = 64540·(1/(0,995·0,985)–1) = 1303 кДж/с;
· потери с рабочим телом (выпар, продувка) с отсчётом от энтальпии подпитки
Таблица 4. Результаты итеративного расчёта одноконтурной ПГУ
Наименование величины | Обозначение | ЕИ | Значение | Расхождение с ном. режимом |
Температура в конденсаторе | °C | 33,9 | -3,1 | |
Давление там же | кПа | 5,3 | -16% | |
Энтальпия насыщенной воды там же | кДж/кг | 142,16 | -12,84 | |
Энтальпия насыщенного пара там же | кДж/кг | 2562,66 | -5,51 | |
Температура пара перед турбиной | °C | 525 | -10 | |
Давление там же | МПа | 5,4 | -0,30% | |
Энтальпия там же | кДж/кг | 3488,41 | -0,66% | |
Давление в барабане | МПа | 6,12 | -0,30% | |
Температура там же | °C | 277 | -0,07% | |
Энтальпия насыщенной воды там же | кДж/кг | 1220,71 | -0,08% | |
Энтальпия насыщенного пара там же | кДж/кг | 2783,18 | 0,01% | |
Энтальпия на выходе из экономайзера | кДж/кг | 1176,74 | 0,55% | |
Объём воды в циркуляционном насосе | м³/кг | 0,00132 | -0,05% | |
Приращение энтальпии там же | кДж/кг | 1,4 | -0,35% | |
Температура газов на выходе из испарителя (pinch-point) | °C | 287 | -0,05% | |
Энтальпия газов там же | кДж/кг | 590,03 | 0,00% | |
Энтальпия воды на выходе из ГПК | кДж/кг | 331,44 | -0,80% | |
Расход греющего пара на обогрев 1 кг конденсата | – | 0,05 | 2,44% | |
Удельный расход конденсата | – | 1,023 | 0,00% | |
Удельный расход греющего пара | – | 0,049 | 2,44% | |
" – " – " выпара | – | 0,00153 | 0,00% | |
" – " – " продувки | – | 0,021 | 0,11% | |
" – " – " добавочной воды | – | 0,0225 | 0,10% | |
" – " – " питательной воды | – | 1,069 | 0,11% | |
Тепловая мощность испарителя | кДж/с | 48581 | 0,04% | |
Удельный расход пара | кг/кг | 0,128 | -3,51% | |
Расход свежего пара | кг/с | 28,92 | 0,32% | |
т/ч | 8,03 | |||
Температура на выходе из ГПК | °C | 79 | -0,79% | |
Тепловая мощность ГПК | кДж/с | 5567 | 5,45% | |
Приращение энтальпии в питательном насосе | кДж/кг | 9,32 | -0,31% | |
Тепловая мощность экономайзера | кДж/с | 22523 | 1,34% | |
Энтальпия уходящих газов | кДж/кг | 466 | 0,48% | |
Температура уходящих газов | °C | 173 | 1,98 | |
Энтальпия газов за пароперегревателем | кДж/кг | 804,54 | -1,04% | |
Температура газов там же | °C | 478 | -1,55% | |
Энтальпия газов за экономайзером | кДж/кг | 490,58 | 0,52% | |
Температура газов там же | °C | 196 | 1,19% | |
Температура на выходе из экономайзера | °C | 268 | 0,48% | |
Давление пара за клапанами | МПа | 5,13 | -0,30% | |
Энтропия свежего пара | 7,0121 | -0,39% | ||
Энтропия пара за клапанами | 7,0349 | -0,39% | ||
Энтальпия конца адиабатного процесса расширения | кДж/кг | 2275,01 | -0,49% | |
Энтальпия в конденсаторе | кДж/кг | 2357,57 | 0,28% | |
Энтропия там же | 7,7 | 1,14% | ||
Влажность там же | – | 0,084934 | -0,005 | |
Внутренняя мощность паровой турбины | МВт | 65415 | -2,26% | |
Электрическая мощность ПСУ | МВт | 64112 | -2,26% | |
Электрическая мощность ПГУ | МВт | 222969 | 2,74% | |
КПД ПГУ | – | 0,500 | -0,0090 | |
Минимальный температурный напор | град. | 10,07 | 0,07 | |
Недогрев до кипения в экономайзере | град. | 8,52 | -1,48 |
3. Расчёт двухконтурной утилизационной ПСУ