Что считает калькулятор
Свойства воды и водяного пара считаются по промышленной формуляции IAPWS-IF97 (International Association for the Properties of Water and Steam, Revised Release R7-97(2012)). Это стандарт, по которому составлены современные таблицы и h-s диаграммы: в том числе таблицы А. А. Александрова и Б. А. Григорьева. Коэффициенты уравнений перенесены из официального текста IAPWS, результат проверен по контрольным таблицам самой формуляции: 72 значения свойств (по 18 в таблицах 5, 15, 33 и 42 для областей 1, 2, 3 и 5) и 6 точек линии насыщения (таблицы 35 и 36) совпадают во всех 9 значащих цифрах.
Области IAPWS-IF97
- Область 1 — вода от 0 до 350 °C при давлении выше давления насыщения, до 100 МПа. Уравнение для энергии Гиббса g(p, T).
- Область 2 — пар от 0 до 800 °C при давлении ниже давления насыщения (и ниже границы B23 при 350…590 °C). Энергия Гиббса: идеально-газовая и остаточная части.
- Область 3 — окрестность критической точки (350…590 °C, высокие давления). Уравнение для энергии Гельмгольца f(ρ, T); плотность по p и T находится итерациями.
- Область 4 — линия насыщения от 0 °C до критической точки (373,946 °C, 22,064 МПа): уравнения p_s(T) и T_s(p). Влажный пар — смесь насыщенных воды и пара со степенью сухости x.
- Область 5 — пар от 800 до 2000 °C при давлении до 50 МПа.
Как решаются обратные задачи
Основные уравнения IF97 заданы в переменных (p, T) или (ρ, T). Для пар (p, h) и (p, s) температура находится итерациями (метод хорд с бисекцией), для пары (h, s) — ещё и давление: при s = const энтальпия монотонно растёт с давлением. В IF97 есть и быстрые обратные уравнения T(p, h), T(p, s), но они отличаются от основных на сотые доли кельвина; итерации дают согласованный с основными уравнениями ответ.
Процессы на h-s диаграмме
- Расширение в турбине. Идеальный процесс — изоэнтропа 1–2s: s₂s = s₁, располагаемый теплоперепад H₀ = h₁ − h₂s. Действительный процесс с внутренним КПД ηoi: h₂ = h₁ − ηoi·H₀, работа l = h₁ − h₂, мощность N = G·l.
- Дросселирование (вентиль, клапан): h₂ = h₁ — горизонталь на h-s диаграмме; давление падает, энтропия растёт, перегретый пар обычно охлаждается.
Пример
Пар 10 МПа, 500 °C: h₁ = 3375,1 кДж/кг, s₁ = 6,5993 кДж/(кг·К). Расширение до 10 кПа: изоэнтропный конец влажный, x₂s = 0,793, h₂s = 2089,6 кДж/кг, H₀ = 1285,4 кДж/кг. С ηoi = 0,85: h₂ = 2282,5 кДж/кг, x₂ = 0,874, l = 1092,6 кДж/кг; при G = 100 кг/с мощность около 109 МВт.
Ограничения
Точность — точность самой IF97 (для большинства свойств доли процента и лучше). Метастабильные состояния (переохлаждённый пар, перегретая жидкость) и лёд не считаются. Вблизи критической точки (в пределах долей кельвина от 373,946 °C) двухфазные значения v′, v″ находятся по уравнению области 3 и менее точны. Вязкость, теплопроводность и поверхностное натяжение — отдельные формуляции IAPWS, здесь их нет.
Источник: IAPWS R7-97(2012), Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam, iapws.org; Wagner W., Kruse A. Properties of Water and Steam. Springer, 1998; Александров А. А., Григорьев Б. А. Таблицы теплофизических свойств воды и водяного пара. М.: МЭИ.