ВведениеВ последнее время в мировой энергетической политике наблюдается тенденция к уменьшению карбонового следа (или выбросов CO2) от промышленности. Добиться уменьшения такого рода выбросов можно несколькими способами. Первый способ, самый популярный в настоящее время, – увеличение доли возобновляемых источников энергии (ВИЭ) в энергетическом балансе. Способ имеет как положительные, так и отрицательные стороны. К минусам относятся высокие капиталозатраты в проекты возобновляемой энергетики. Также к минусам можно отнести непостоянность выработки энергетического ресурса возобновляемыми источниками, обусловливаемого высокой зависимостью от погодных условий. Второй способ, который на данный момент является не таким популярным как ВИЭ, – повышение энергоэффективности промышленности. Для России второй способ является наиболее приемлемым, т. к. в нашей стране в энергетическом балансе преобладает промышленность, для которой важно постоянство получения энергоресурсов, а соответственно, независимость от погодных условий. В современных условиях существует большое количество способов повышения энергоэффективности различных технологических процессов. В настоящей работе остановимся на таком способе повышения энергоэффективности, как использование установок органического цикла Ренкина (ОЦР) для утилизации теплоты низкопотенциальных источников. Материалы исследованияПростейшая схема ОЦР-установки представляет собой классический цикл Ренкина. Главное отличие заключается в использовании в качестве рабочего тела органических веществ, которые кипят при более низких температурах, чем классическое рабочее тело в цикле Ренкина – вода. Принципиальная схема простейшей ОЦР-установки представлена на рис. 1.   Рис. 1. Принципиальная схема ОЦР-установки Согласно рис. 1 простейшая ОЦР-установка состоит из питательного насоса, в котором рабочее тело сжимается до необходимого давления; испарителя, в котором происходит передача рабочему телу энергии от внешнего низкопотенциального источника; паровой турбины, в которой рабочее тело расширяется, и конденсатора, в котором происходит отвод тепловой энергии от парожидкостной смеси и конденсация рабочего тела. В настоящей работе рассмотрим подробнее испаритель ОЦР-установки.Испаритель чаще всего представляет собой прямоточный котел, в котором происходит парообразование. Парообразование в испарителях ОЦР-установок не изучено в полной мере. Одним из вопросов исследования теплообмена в испарителе, который недостаточно освещен, является влияние физических свойств теплоносителя на коэффициент теплоотдачи.Сегодня большинство исследований посвящено изучению характера изменения энергетической эффективности ОЦР-установок в зависимости от вещества теплоносителя, а соответственно, и от его физических свойств. В работе [1] было установлено, что для диапазона температур низкопотенциального источника теплоты 150–250 °С для установок ОЦР в качестве теплоносителя наиболее выгодно использовать пентан (C5H12), этанол (C2H5OH) и метанол (CH3OH). Эти вещества и будем рассматривать далее в работе. Описание моделиПроведено исследование [2] зависимости коэффициента теплоотдачи от длины участка испарителя ОЦР-установки. Для настоящего исследования модель будет аналогичной исследованию [2]. Данная модель представляет собой прямой участок трубы испарителя ОЦР-установки (длина – 90 см, толщина стенки – 4,5 мм, внутренний диаметр трубы – 40 мм). В качестве материала трубы была задана сталь, в качестве условий однозначности – граничные условия 3-го рода. Граничные условия приняты в соответствии с рекомендациями работ [3–9]. Внешний вид модели представлен на рис. 2.   Рис. 2. Внешний вид исследуемой модели  Для исследования были выбраны рабочие тела, параметры которых представлены в таблице.Параметры рассматриваемых рабочих телВеществоФормулаПоверхностное натяжение, Н/мТемпература кипения при нормальных условиях, КТеплота парообразования, кДж/кгПентанC5H120,0155075309357,5ЭтанолC2H5OH0,022348351855,2МетанолCH3OH0,02226573381 100,3 Каждое из исследуемых рабочих тел в настоящее время используется в качестве рабочего тела в различных ОЦР-установках с температурой низкопотенциального источника теплоты от 100 до 200 °С.Тепловая энергия к рабочему телу передается от стенок трубы, имеющих температуру 440 К. Скорость теплоносителя на входе в трубу испарителя была задана 0,125 м/с. Температура теплоносителя на входе в испаритель задавалась на несколько градусов ниже температуры насыщения. Временной шаг принят 0,01 с. Количество итераций – 500. Таким образом было смоделировано 5 с процесса кипения теплоносителя – от состояния жидкости до установившегося режима. Данные условия сохранялись постоянными при изменении вида теплоносителя. Исследовалась только зависимость коэффициента теплоотдачи различных теплоносителей от их физических свойств. Расчет производился в версии 19R3 программного комплекса ANSYS с помощью модуля Fluent.Результаты исследованияВ качестве первичной обработки результатов были рассмотрены картины кипения теплоносителей. Пример полученной картины кипения для метанола представлен на рис. 3.   a                                                                      б Рис. 3. Полученная картина кипения: а – итерация – 200; б – итерация – 500 В результате моделирования были получены значения коэффициента теплоотдачи в каждой точке участка испарителя. После этого была определена область, в которой происходит установившийся режим течения с помощью анализа изменения коэффициента теплоотдачи во времени. Далее приведены зависимости именно для установившегося режима течения.Также в результате анализа были получены зависимости изменения коэффициента теплоотдачи по длине трубы испарителя для каждого вещества теплоносителя (рис. 4).  l, мα, Вт/м2·К1 600 1 400 1 200 1 000 800 600 400 200 0   Рис. 4. Изменение коэффициента теплоотдачи по длине трубы испарителяОбсуждение результатовПосле получения зависимости изменения коэффициента теплоотдачи по длине трубы испарителя были построены зависимости коэффициента теплоотдачи от удельной теплоты парообразования теплоносителя для участка трубы испарителя 0,7–0,89 м (рис. 5).   Рис. 5. Зависимость коэффициента теплоотдачи от удельной теплоты парообразования Также была получена зависимость изменения коэффициента теплоотдачи от силы поверхностного натяжения теплоносителя (рис. 6).   Рис. 6. Зависимость коэффициента теплоотдачи от силы поверхностного натяжения После получения зависимостей, представленных выше, была получена зависимость коэффициента теплоотдачи от отношения силы поверхностного натяжения к теплоте парообразования материала теплоносителя (рис. 7).  Рис. 7. Зависимость коэффициента теплоотдачи от физических свойств материала теплоносителя На рис. 4 можно заметить, что при увеличении длины трубы коэффициент теплоотдачи падает. Это объясняется тем, что пар заполняет практически полностью конечный участок трубы испарителя. Пар заполняет практически всю пристенную область трубы. Коэффициент теплоотдачи при этом падает из-за низкой теплопроводности пара по сравнению с жидкостью. Также можно заметить, что величина падения коэффициента теплоотдачи для разных видов теплоносителя отличается, что свидетельствует о его зависимости от физических свойств рабочего тела. Отчетливо видно (см. рис. 4), что значение падения коэффициента теплоотдачи по длине теплоносителя для пентана значительно больше аналогичной величины для метанола.Согласно графику зависимости коэффициента теплоотдачи от удельной теплоты парообразования (см. рис. 5) при увеличении значения удельной теплоты парообразования коэффициент теплоотдачи падает. При увеличении силы поверхностного натяжения рабочего тела (см. рис. 6) коэффициент теплоотдачи также снижается.Описанные выше явления можно объяснить следующим образом. При увеличении теплоты парообразования для образования пузырька с паром необходимо осуществить подвод большего количества теплоты. Соответственно, число центров парообразования и размеры пузырьков пара становятся меньше, чем в случае с более низкой теплотой парообразования теплоносителя. Подобное явление наблюдается и в случае с увеличением силы поверхностного натяжения вещества теплоносителя. При ее увеличении для образования пузырька с паром необходимо передать больше энергии, которая затрачивается на преодоление силы поверхностного натяжения. Таким образом, в обоих случаях уменьшается интенсивность образования пузырьков пара. Степень сухости пара становится больше, а значит, коэффициент теплоотдачи начинает падать.В ходе анализа получившихся зависимостей коэффициента теплоотдачи от отношения силы поверхностного натяжения к удельной теплоте парообразования было получено эмпирическое уравнение для расчета коэффициента теплоотдачи:                          Уравнение позволяет вычислять коэффициент теплоотдачи во время процесса кипения органического теплоносителя при диаметре трубы 40 мм и в диапазоне начальной скорости потока 0,1–0,125 м/с. Погрешность значений, полученных с помощью данного уравнения, не превышает 12,1 % от экспериментально полученной величины. ЗаключениеОдним из путей энергосбережения в судовых энергетических установках является утилизация теплоты выхлопных газов. Генерация электроэнергии при относительно низких температурах выхлопных газов возможна только при использовании рабочих тел с высокими давлениями в этом диапазоне температур. Оптимальным техническим решением для этих целей можно считать применение ОЦР. Значительная неравномерность коэффициента теплоотдачи по длине трубы прямоточного утилизационного котла является основной задачей при проектировании вспомогательной судовой энергетической установки на основе ОЦР. В связи с этим в работе рассмотрена проблема влияния теплофизических свойств теплоносителя ОЦР-установки на коэффициент теплоотдачи во время кипения теплоносителя в испарителе. Проиллюстрирована картина кипения органического теплоносителя в трубах испарителя установки ОЦР. В результате исследования было определено, что коэффициент теплоотдачи уменьшается при увеличении длины испарителя. Также была определена зависимость коэффициента теплоотдачи во время процесса кипения от таких физических свойств вещества теплоносителя, как теплота парообразования и коэффициент поверхностного натяжения. Определено, что при увеличении отношения силы поверхностного натяжения к удельной теплоте парообразования коэффициент теплоотдачи увеличивается. После изучения данных зависимостей было получено эмпирическое уравнение для расчета коэффициента теплоотдачи во время кипения органического теплоносителя, учитывающее его изменение по длине трубы испарителя и физические свойства вещества теплоносителя.