CC BY
26
УДК 536.24
Моделирование циклически сопряженного теплообмена в регенеративном воздухоподогревателе
Елин Н.Н., д-р техн. наук, Ометова М.Ю., канд. техн. наук, Рыбкина Г.В., инж.
Предложена математическая модель теплообмена в регенеративном воздухоподогревателе с неподвижной кирпичной насадкой, базирующаяся на системе дифференциальных уравнений в частных производных с граничными условиями для конвективно-радиационного теплообмена, позволяющая рассчитывать нестационарное температурное поле насадки и определять оптимальные режимные параметры процесса.
Ключевые слова: регенератор, теплообменник, нестационарный теплообмен.
Modelling Cyclically Conjugated Heat Exchange in Regenerative Air Heater
N.N. Yelin, Doctor of Engineering, M.Yu. Ometova, Candidate of Engineering, G.B. Rybkina, Engineer
The author considers a heat exchange mathematical model in regenerative air heater with fixed brick nozzle based on differential equations system in partial derivatives with boundary conditions for convective and radiative heat exchange. The model allows to calculate unsteady temperature field of the nozzle and to define optimal operating conditions of the process.
Key words: regenerator, heat exchanger, unsteady heat transfer.
Одна из важнейших проблем при использовании регенеративных теплообменников заключается в определении оптимальной продолжительности цикла, при которой достигается максимальная тепловая мощность аппарата, то есть количество теплоты, передаваемое от уходящих продуктов сгорания к подогреваемому воздуху, отнесенное к общему времени цикла.
Данная задача не может быть решена с помощью известных методов расчета, в которых вместо актуальных значений температур теплоносителей используются значения, ос-редненные по длине канала и времени цикла, а также ряд других допущений. Экспериментальные исследования регенеративных теплообменников представляют большие сложности вследствие большого количества конструктивных и режимных параметров, влияющих на процесс теплообмена. Поэтому основным методом исследования процессов, происходящих в регенеративных теплообменниках, стало их математическое моделирование.
В качестве объекта моделирования рассмотрим элемент насадки регенератора, представленный на рис. 1.
Математическая модель регенеративного воздухоподогревателя включает следующие уравнения:
1) уравнение энергии для насадки
5Г1
~дТ
(
д2Т
д2ТЛ
д2Т
\
(1)
чдх2 ду2 дг* ^
2) уравнение энергии для теплоносителя
дТ2 _ __д_ д ? д у'
где Т - температура; р - плот ность; с - т е пло-емкость; а - температуропроводность; т - ско-
С 2 р 2 ■
-(с2 р2 Ю 2 Т2) >
(2)
рость; индекс «1» относится к насадке, «2» - к теплоносителю;
3) уравнение неразрывности для теплоносителя
р 2 Ю 2 Э _ С 2, (3)
где в2 - массовый расход теплоносителя, постоянный по длине и во времени в течение полуцикла нагрева или охлаждения.
Рис. 1. Схема расчетной области
Для расчета изменения плотности теп лоносителя вдоль потока используется урав нение состояния идеального газа
Р
Р2 _ И Т'
(4)
Расчеты [1] показали, что изменением давления по длине теплообменного аппарата можно пренебречь.
Граничные условия на поверхности насадки при течении продуктов сгорания следующие:
д Т
^ = ак (Т2 — T1w ) + Яизл ,
д n
(5)
где Хн - коэффициент теплопроводности насадки; ак - коэффициент конвективного теплообмена [2]; дизл - плотность теплового потока излучения газа [4]; n = x,z; индекс «w» относится к поверхности насадки.
При движении воздуха вторым слагаемым в правой части (5) можно пренебречь [1].
Начальные условия выбраны следующими: Tiw (,y,z,0) = Tiwо (,y,z); (6)
T2 (0, t) = T2ex ; (7)
T2 (y,0) = T20 (y ), (8)
где Т2вх - температура теплоносителя на входе в канал; индекс «0» относится к распределению температур в начале полуцикла нагрева или охлаждения.
Численная реализация данной модели осуществлялась методом сеток с применением явных схем аппроксимации [3]. Контроль точности вычислений проводился по балансу тепловой энергии.
Численные эксперименты с использованием данной математической модели выполнялись для регенеративного теплообменника, имеющего следующие геометрические размеры насадки, изготовленной из шамотного кирпича: L = 2 м, h1 = 50 мм, h2 = 80 мм. Теплоносителями являются воздух и дымовые газы, образующиеся при сжигании топлива заданного состава. Атмосферный воздух необходимо нагреть от 20 до 4500С дымовыми газами с температурой 9000С. Скорости потоков теплоносителей на входе равны ю0 = 1 м/с. Теплофизические свойства теплоносителей, необходимые для замыкания системы уравнений (1)-(8), рассчитывались по [2].
На рис. 2 представлены результаты расчета процесса прогрева проточной части регенератора, насадка которого в начальный момент времени имеет температуру 200С.
Т, 0С
I, см
Рис. 2. Распределение температуры горячего теплоносителя по длине проточной части регенератора для различных моментов времени: 1 - t = 2 с; 2 - t = 1000 с; 3 - t = 2000 с; 4 - t = 3000 с
В начальный момент распределение температуры горячего теплоносителя по длине проточной части регенератора имеет прямоугольную форму высотой Т2вх = 9000С. По мере поступления теплоносителя в канал форма
профиля становится все более пологой, средняя по длине канала температура теплоносителя быстро уменьшается (тонкие линии на рис. 2 представляют результаты расчетов профиля температуры через каждые 0,2 с). Через 2 с, когда поток теплоносителя достигает выходного сечения канала, температура теплоносителя на выходе становится равной 2300С, а распределение температуры по длине имеет вид, представленный кривой 1 (рис. 2). Такое относительно быстрое охлаждение теплоносителя обусловлено низкой начальной температурой насадки, большими температурными напорами, а следовательно, большой плотностью теплового потока от теплоносителя к насадке.
После того как поток теплоносителя заполнил проточную часть регенератора, начинается относительно медленное увеличение температуры газа до 5000С, соответствующей установившемуся режиму работы, в течение 3000 с (рис. 2, кривые 2, 3, 4).
На рис. 3 представлены результаты расчета температурных полей в поперечных сечениях насадки в различные моменты времени. Анализ полученных результатов показывает, что углы насадки практически не участвуют в теплообмене, а необходимы, скорее всего, для механической прочности конструкции. Этот результат дает основание предположить, что наиболее рациональной формой проходного сечения насадки является щелевая, а наименее рациональной - квадратная.
б)
750 650 550 450 350 250
Рис. 3. Температурное поле в трех поперечных сечениях насадки: а - через 50 с после начала нагрева; б - в конце периода нагрева (через 400 с)
В процессе численного интегрирования системы уравнений (1)-(8) величины коэффициентов теплопереноса рассчитывались на каждом шаге по длине регенератора и по времени. Результаты расчетов распределения коэффициентов теплоотдачи от горячего теплоносителя к насадке по длине проточной части для установившегося режима работы представлены на рис. 4.
На входе в канал доли лучистого и конвективного теплообмена для рассматриваемого примера практически одинаковы. В процессе движения теплоноситель охлаждается, и это влияет на величину коэффициентов теплоотдачи.
а, Вт/(м2С)
, см
Рис. 4. Изменение средних коэффициентов теплоотдачи по длине теплообменного аппарата
Поток тепла радиацией, пропорциональный разнице четвертых степеней температур теплоносителя и насадки, при этом уменьшается наиболее быстро. В выходном сечении проточной части он становится почти в два раза меньше, чем конвективный.
Поток тепла конвекцией заметно уменьшается на участке термической стабилизации потока (равной примерно 1/2 длины проточной части теплообменного аппарата), а затем остается практически постоянным, несмотря на уменьшение скорости потока вследствие охлаждения газа (это объясняется тем, что при этом величина критерия Рейнольдса практически не изменяется вследствие уменьшения вязкости газа).
Важнейшим показателем энергетической эффективности регенеративного теплообменника является его тепловая мощность.
На рис. 5 представлены результаты исследования влияния времени цикла на тепловую мощность регенеративного воздухоподогревателя.
Известно, что существенное влияние на этот показатель эффективности оказывает время перекидки клапанов [3]. Поэтому при исследовании зависимости тепловой мощности регенеративного воздухоподогревателя N от времени цикла в качестве параметра задачи использовалась продолжительность паузы для перекидки клапанов.
N, МВт
т, с
Рис. 5. Зависимости тепловой мощности от длительности периода цикла при различной длительности паузы перекидки клапанов: 1 - тпі = 10 с; 2 - Хщ = 20 с; 3 - Тпз = 30 с
Резкий подъем кривых в области малой длительности цикла объясняется большой долей продолжительности паузы перекидки клапанов в общем времени цикла. При увеличении длительности цикла это влияние уменьшается. Максимальное значение мощности достигается тогда, когда сумма потерь теплоты, обусловленная паузой для перекидки клапанов и потерь теплоты с уходящими дымовыми газами, минимальна.
Предложенная математическая модель позволяет исследовать влияние теплофизических, геометрических и режимных параметров на процессы теплообмена в регенеративном воздухоподогревателе, а также определить оптимальную продолжительность цикла, при которой достигается максимальная тепловая мощность аппарата.
Список литературы:
1. Хавер С.В., Елин Н.Н., Мизонов В.Е., Иванов А.Б.
Моделирование нагрева и охлаждения насадки регенеративного теплообменника // Изв. вузов. Химия и химическая технология. - 2007. - Т. 50. - Вып. 12. - С. 105-107.
2. Тепловой расчет котельных агрегатов (Нормативный метод) / Под ред. Н.В. Кузнецова, В.В. Митора, И.Е. Дубовско-го, Э.С. Карасиной. - М.: Энергия, 1973.
3. Лемлех И.М. Высокотемпературный нагрев воздуха в черной металлургии. - М.: Металлургиздат, 1983.
4. Самарский А.А. Теория разностных схем. - М.: Наука, 1977.
Елин Николай Николаевич,
Ивановский государственный архитектурно-строительный университет,
доктор технических наук, профессор, зав. кафедрой гидравлики и водоснабжения,
адрес: г. Иваново, ул. Громобоя, д. 32, кв. 16
телефон (4932) 41-61-22,
e-mail: [email protected]
Рыбкина Галина Владимировна,
Ивановский государственный архитектурно-строительный университет, инженер ФКПиПКС,
адрес: г. Иваново, ул. Свердлова, д. 10, кв. 2,
телефон 46-15-26,
e-mail: [email protected]
Ометова Мария Юрьевна,
Ивановский государственный архитектурно-строительный университет, кандидат технических наук, доцент каф. гидравлики и водоснабжения, адрес: г. Иваново, ул. 4 Деревенская, д. 62, кв. 47, телефон (4932) 47-20-13, e-mail: ometova_m@ rambler.ru
