Цикл Ренкина – термодинамический цикл, по которому работают паросиловые или паротурбинные установки (ПТУ), являющиеся в настоящее время основными тепловыми агрегатами, производящими электрическую энергию в мире [1]. Соответственно, разработка любой паротурбинной установки начинается с нахождения расположения основных (характерных) точек цикла.
Чаще всего построение цикла Ренкина производится в так называемой i-s диаграмме, по оси абсцисс которой откладывается энтропия системы, а по оси ординат – энтальпия (рис. 1).
Цикл состоит из нескольких термодинамических процессов, протекающих последовательно и отделяемых точками их начала и окончания [2]:
1–2 – адиабатное расширение пара в турбине (изменение при этом энтальпии в турбине H0 называют располагаемым теплоперепадом);
1–2д – действительное расширение пара в турбине (изменение энтальпии при этом в турбине Hд называют использованным или действительным теплоперепадом);
2–3 – изобарный (и изотермический, если он происходит в области влажного пара) отвод теплоты в конденсаторе;
Рис. 1. Цикл Ренкина с характерными точками
3–4 – адиабатное повышение давления в питательном насосе;
4–1 – изобарный подвод теплоты в паровом котле, который состоит из трёх участков:
4–5 – нагрев воды до кипения;
5–6 – парообразование;
6–1 – перегрев пара.
Буквой К обозначена критическая точка, выше которой водяная фаза не может существовать и водяной пар ведёт себя как газ – заполняет весь выделенный объём [3]. Такой пар называют перегретым.
Левая пограничная линия (левее точки К) разделяет области воды и пара со степенью сухости x = 0. Правая пограничная линия (правее точки К) разделяет области перегретого пара и пара со степенью сухости x = 1.
С помощью q1 обозначена теплота, подведённая к рабочему телу (воде) в котле, q2 – теплота, отведённая в конденсаторе.
Изображение цикла Ренкина в i-s диаграмме позволяет не только наглядно показать протекание основных процессов в паросиловой установке, но и найти параметры, характеризующие эффективность цикла. Вычислив влияние координат характерных точек цикла на эффективность, можно целенаправленно оптимизировать параметры цикла.
Таким образом, для повышения качества проектирования паросиловых установок, работающих по циклу Ренкина, необходим быстрый и малозатратный способ нахождения характерных точек данного цикла [4]. В связи с этим целью данной работы является автоматизация построения необходимых зависимостей и характерных точек с помощью специально разработанной компьютерной программы.
Материалы и методы исследования
Для создания расчётной программы использовался алгоритм определения положения точек начала и окончания термодинамических процессов ПТУ, а также различных типов коэффициентов полезного действия (КПД) установки, приведённый в работе [5].
Средством разработки программного продукта являлся получивший широкое распространение в настоящее время язык программирования Python. В данном случае он был выбран благодаря своим неоспоримым достоинствам:
- лёгкость написания программ научной направленности;
– большой набор средств для создания графического интерфейса пользователя (GUI);
- возможность вывода графиков, обладающих высоким качеством, а также автоматическое их сохранение в виде файлов распространённых графических форматов во время выполнения программы;
- возможность создания программы в виде одного исполняемого файла.
Но главным критерием выбора языка Python являлось то, что для него существует мощный программный модуль CoolProp, позволяющий получать физические и, в частности, термодинамические свойства многих веществ. Текущая версия содержит свойства 122 веществ, среди которых в данном исследовании было использовано вещество вода-пар. Для аппроксимации свойств воды и водяного пара в модуле CoolProp используется новейшая на данный момент методика IF97 [6].
Рис. 2. Интерфейс программы для расчёта цикла Ренкина
Основное окно программы для расчёта цикла Ренкина было создано с помощью подключённого модуля – библиотеки tkinter языка Python. Командами данного модуля созданы пояснительные надписи, поля для ввода исходных данных, кнопка запуска расчёта, [7]. Вывод результатов расчёта цикла реализован в виде графиков с использованием модуля matplotlib [8].
Созданный интерфейс программы показан на рис. 2.
Как видно из рисунка, для расчёта цикла необходимо ввести лишь начальные и конечные параметры водяного пара, а также КПД элементов паротурбинной установки.
Результаты исследования и их обсуждение
Для отработки функций программы использовался приведённый в [5] пример расчёта основных точек и характеристик цикла Ренкина. Исходные данные согласно данному примеру введены в главное окно программы и запущен расчёт.
На рис. 3 приведён построенный в i-s диаграмме с помощью программы цикл.
В верхней части изображения также выводятся два КПД: термический КПД цикла, определяющий его термодинамическое совершенство, и эффективный КПД действительного цикла, полученный с учётом дополнительных потерь в котле, паропроводах, в электрогенераторе, в насосе, механических потерь в турбине.
Первый, термический, КПД цикла зависит в большой степени от выбора положения основных точек цикла, которое в некоторых пределах можно изменять. Для демонстрации этого влияния написан ещё один вариант рассматриваемой программы. В нём осуществляется вывод графика зависимости термического КПД от начального давления (в точке 1 – перед турбиной).
На рис. 4 показана данная зависимость при начальной температуре t1 = 550 °C (такой же, как и при построении цикла в вышеприведённом примере). Давление p1 изменялось от 5 до 50 МПа. Остальные исходные данные те же, что и на рис. 1.
Рисунок демонстрирует, что максимального значения термический КПД цикла Ренкина в данном случае достигает при давлении p1, примерно равном 38 МПа.
Подобный график можно вывести и для любых других исходных данных. Если снизить начальную температуру t1 до 450 °C, то программа выдаст картину, изображённую на рис. 5.
Из рисунка видно, что при снижении начальной температуры максимальный термический КПД цикла получается при меньшем значении начального давления, равном примерно 24 МПа. Максимально достижимый термический КПД цикла при этом снижается на 2,5 %.
Данное явление влияния начальных параметров цикла на КПД описано в литературе по термодинамике, например в [9].
Подобным образом с помощью программы на языке Python с подключением описанных выше модулей-библиотек можно рассчитать и построить любые другие характеристики цикла Ренкина.
Рис. 3. Точки цикла Ренкина и КПД, полученные в программе
Рис. 4. Зависимость термического КПД цикла от давления пара перед турбиной при начальной температуре 550 °C
Рис. 5. Зависимость термического КПД цикла от давления пара перед турбиной при начальной температуре 450 °C
Процесс написания и использования компьютерного приложения доказал удобство применения языка Python для реализации алгоритмов научных расчётов. Расчёт точек и характеристик цикла Ренкина, выполненный с помощью программы в сравнении с примером, описанным в статье [5], показал полное соответствие полученных результатов.
Следует отметить уникальную возможность использования модуля CoolProp при нахождении физических и термодинамических свойств веществ, и, пожалуй, нет больше ни одного свободно распространяемого программного обеспечения, реализующего в полной мере эту возможность.
Основное предназначение разработанной программы видится в оптимизации цикла паротурбинной установки с целью получения наибольшей её эффективности при дальнейшей детальной разработке.
Заключение
Таким образом, в рамках описанной работы с использованием приведённого в литературе алгоритма расчёта цикла Ренкина разработана программа на языке Python с использованием дополнительных программных модулей, позволяющая решить следующие задачи:
- построить цикл существующей паротурбинной установки с указанием точек начала и окончания протекания всех термодинамических процессов;
- подобрать характеристики разрабатываемого цикла Ренкина с целью получения наивысшей экономичности паросиловой установки;
- использовать программу для наглядной демонстрации взаимного влияния исходных факторов и результирующих характеристик цикла.
С использованием разработанной программы можно производить как проверку характеристик уже используемой установки при каком-либо изменении исходных параметров, так и проектировать новые паротурбинные установки, работающие по циклу Ренкина. Такая автоматизация процесса проектирования приведёт к повышению надёжности и эффективности работы ПТУ [10]. Также данную программу можно рекомендовать использовать для обучения студентов, персонала тепловых станций основам протекания термодинамических процессов и как иллюстрацию применения языка Python в создании научно ориентированных приложений.
Библиографическая ссылка
Ильичев В.Ю., Юрик Е.А. РАЗРАБОТКА ПРОГРАММЫ ДЛЯ ИССЛЕДОВАНИЯ ТЕРМОДИНАМИЧЕСКОГО ЦИКЛА РЕНКИНА // Научное обозрение. Технические науки. – 2020. – № 2. – С. 32-36;URL: https://science-engineering.ru/ru/article/view?id=1284 (дата обращения: 06.11.2024).