CC BY
82
УДК 621.643.001:536.2
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ТЕПЛОВЫХ РЕЖИМОВ КАНАЛЬНЫХ ТЕПЛОТРУБОПРОВОДОВ В УСЛОВИЯХ ВЗАИМОДЕЙСТВИЯ С ВЛАЖНЫМ ВОЗДУХОМ
В.С. Логинов, В.Ю. Половников
Томский политехнический университет E-mail: [email protected]
Проведено численное исследование тепломассопереноса в тепловой изоляции теплотрубопровода, проложенного канальным способом, в условиях взаимодействия с влажным воздухом. Установленыi закономерности тепломассообмена и масштабыi тепловых потерь канального теплотрубопровода с влажным воздухом.
Ключевые слова:
Математическое моделирование, теплотрубопровод, тепловые потери, диффузия, массообмен.
Введение
Численное моделирование тепловых режимов теплопроводов является актуальной задачей при разработке энергосберегающих систем передачи тепловой энергии. При этом характеристика те-плотрубопроводов по величине тепловых потерь является особо важным показателем, отражающим действительное состояние тепловых сетей.
В настоящее время основной причиной роста тепловых потерь трубопроводов является [1, 2] изменение термовлажностных условий работы теплосетей и в первую очередь - эксплуатация теплопроводов, имеющих влагонасыщенную теплоизоляцию. При этом наиболее типичным механизмом увлажнения тепловой изоляции теплопроводов является [2] насыщение теплоизоляционного материала влагой, содержащейся в воздухе канального пространства.
Целью данной работы является математическое моделирование тепломассопереноса в тепловой изоляции трубопровода канальной прокладки и численный анализ тепловых потерь при его эксплуатации в условиях увлажнения изоляции влагой, содержащейся в воздухе канального пространства.
Постановка задачи
Рассматривается широко распространенная конфигурация теплотрубопровода [1], уложенного в одноячейковый железобетонный канал и теплоизолированного минеральной ватой. Предполагается, что влага, содержащаяся во влажном воздухе, окружающем трубопровод, проникает в пористую структуру теплоизоляции под действием диффузионного механизма массопереноса и вытесняет сухой воздух, содержащийся в слое минераловат-ной теплоизоляции. Следствием влагонасыщения тепловой изоляции является увеличение эффективной теплопроводности теплоизоляционного слоя, а наличие потока массы на внешней поверхности теплоизоляции приводит к изменению механизмов тепломассообмена теплопровода с окружающей средой.
С учетом рассматриваемых факторов анализ тепловых потерь сводится к решению одномерной нестационарной задачи тепло- и массопереноса для тепловой изоляции трубопровода.
При постановке задачи на внешней границе рассматриваемой системы «слой теплоизоляции -окружающая среда» выставлялись граничные условия третьего рода, а на внутренней - граничные условия первого рода для задачи теплопроводности и второго рода для задачи влагопереноса. В начальный момент времени считалось, что влага в пористой структуре тепловой изоляции трубопровода отсутствует, а значение температуры в рассматриваемой области равно температуре наружного влажного воздуха.
Математическая модель
Задача решалась в цилиндрической системе координат, начало которой связано с осью симметрии трубы. Математическая постановка задачи тепло-массопереноса в тепловой изоляции трубопровода описана при помощи уравнений (1)-(9): т> 0, Я1 < г < Я2;
дТ 1 д ( г , ч дТ \
^ иг=гдгК^-дг(1)
т = 0, R <г < Д2; Ги = Гвв; т > 0, г = R,; Ти = T = const;
дТ
(2) (3)
т> 0, г = R2; = а[Т,(т, R2) - TJ; (4)
дг
дСв „ I д Св 1 дСв
т> 0, R, < г < R2; -т = Dl -f- + - I; (5)
дт I дг г дг I
т = 0, R, < г < R2; Св = 0; дС
т> 0, г = R,; —^ = 0; 1 дг
дСв
(6) (7)
т> 0, г = R2; -= р[Св(т,R2)- С ]. (8) дг
Энергетика
Обозначения: T - температура, К; т - время, с; г - текущий радиус, м; с - теплоемкость, Дж/(кг-К); р - плотность, кг/м3; Я - коэффициент теплопроводности, Вт/(м-К); R - радиус границ области решения, м; а - коэффициент теплоотдачи, Вт/(м2-К); C - относительная объемная концентрация; D - коэффициент диффузии, м2/с; в -коэффициент массоотдачи, м/с.
Индексы: 1, 2 - соответственно внутренняя и внешняя поверхность изоляции; и - изоляция; эф - эффективный, вв - влажный воздух; в - влага.
Теплофизические свойства тепловой изоляции при насыщении ее влагой рассчитывались как эффективные с учетом объемных долей каждой компоненты. Так, например, эффективный коэффициент теплопроводности Яф рассчитывался из выражения [3]:
(9)
Ав =
0,083-105
Яф =ЯС +ЯС +Я с
эф и и в в св св
в
= СввРв,
г _ Л-
V Ав у
где а - коэффициент температуропроводности, м2/с.
Коэффициент диффузии влажного воздуха Dвв вычислялся по формуле [8]:
Р„
Т
вв
273
здесь Р - давление, Па.
Объемная доля влаги во влажном воздухе Св,, определялась из соотношения:
Яв,вв ё. Л „ фР,
С =■
ё = 0,622-
при условии что
С + С + С = 1,
и в св >
где Ясв, Ссв - соответственно коэффициент теплопроводности и относительная объемная концентрация сухого воздуха, содержащегося в структуре теплоизоляционного слоя.
Метод решения и исходные данные
Система уравнений (1)-(9) решена методом конечных разностей [4] с использованием неявной разностной схемы, а разностные аналоги уравнений (1)-(9) разрешены методом «прогонки» [4].
Исследования проводились для трубопровода с диаметром условного прохода 600 мм и тепловой изоляцией из минеральной ваты (толщина 70 мм). Температура влажного воздуха, в соответствии с [5], принималась равной Твв=298 К, а температура внутренней поверхности изоляции теплопровода -Т1=373 К.
При моделировании работы теплотрубопровода предполагалось, что на его внешнем контуре теплообмен осуществляется в условиях естественной конвекции. Коэффициенты теплоотдачи а для режимов естественной конвекции определялись по критериальным уравнениям [6]:
104 < вгРг < 109; Ми = 0,47(вгРг)1/4, вгРг > 109; Ми = 0,1(вгРг)1/3,
здесь Gr, Рг, Nu - числа Грасгофа, Прандтля и Нус-сельта.
Коэффициент массоотдачи в рассчитывался из отношения [7]:
Яв рвв -ФР,
где ¡л - мольная масса, й - массовое влагосодержа-ние воздуха, ф - относительная влажность, Р, - давление насыщения.
Линейные тепловые потери теплотрубопро-вода рассчитывались по следующему выражению:
дь =а[Т,(т,-Твв]^>
здесь Р - площадь поверхности, м2; Ь - единица длины, м.
Результаты исследований
Численный анализ проводился для периода времени, соответствующего выходу процессов тепло- и массопереноса на стационарный режим. Обоснованность и достоверность результатов численных исследований следует из проведенных проверок используемых методов решения задач на сходимость и устойчивость, а также выполнения условий баланса энергии и массы на границах области расчета. Погрешность по балансу энергии и массы во всех вариантах численного анализа не превысила 0,5 %, что можно считать приемлемым при моделировании те-пломассопереноса в изоляции теплопроводов.
В таблице приведены типичные результаты численных исследований.
Таблица. Результаты численных исследований тепловых потерь трубопроводов
ф, % 0 25 50 75 100
41, Вт/м 106,7 116,1 121,9 126,9 131,5
ф, Вт/м (Св=0) 106,7 113,3 116,3 118,2 119,7
В результате численного моделирования температурных полей в теплоизоляции трубопровода, проведенного на базе математической модели (1)-(9) было установлено, что тепловые потери теплопровода (табл.) в условиях увлажнения изоляции влагой, содержащейся во влажном воздухе канального пространства, в целом достаточно близки к нормированному значению теплопотерь (по данным [9] нормативная величина линейных тепловых потерь для трубопровода с диаметром условного прохода 600 мм составляет 114 Вт/м).
Рост тепловых потерь трубопровода (табл.) объясняется с одной стороны увеличением эффективной теплопроводности материала изоляции вследствие его насыщения влагой из влажного воздуха, а с другой - увеличением коэффициентов те-
ос
плоотдачи для воздуха, имеющего большую относительную влажность р.
В результате проведения серии численных экспериментов было выявлено, что для диапазона варьируемых в данной работе параметров величина объемной доли влаги в пористой структуре тепловой изоляции изменяется в пределах от 0,0031 до 0,0126, а расчетные значения коэффициентов теплоотдачи а=2,64...5,19 Вт/(м2-К). Здесь следует отметить, что величина объемного влагосодержания изоляции лимитирована значением относительной влажности воздуха и в рамках предложенной модели не может превышать объемной доли влаги Св,вв.
Сравнение результатов численных исследований тепловых потерь теплотрубопровода (табл.) в условиях взаимодействия с влажным воздухом и отсутствия влаги в пористой структуре теплоизоляции (Св=0) позволяет сделать вывод о том, что реализация режима работы теплопровода, исключающего возможность проникновения влаги в слой тепловой изоляции, позволит минимизировать потери тепловой энергии. На практике данное положение может быть реализовано нанесением на внешнюю поверхность тепловой изоляции слоя надежного гидроизоляционного материала. Анализ расчетных значений тепловых потерь, приведенных в таблице, свидетельствует о том, что снижение потерь тепловой энергии теплопроводов в условиях контакта с влажным воздухом при применении надежно гидроизолированных трубопроводов может составлять величину до 9 %.
Анализ нестационарности влагонасыщения тепловой изоляции позволяет говорить о том, что длительность этого процесса является достаточно незначительной, имеет значение от 7,5-103 до 104 с, и определяется значением коэффициента диффузии материала изоляции Dи и объемным влагосо-держанием воздуха Св,вв. Следовательно, нестационарностью процесса насыщения изоляции влагой, содержащейся в воздухе канального пространства, можно обоснованно пренебречь и в практических расчетах использовать стационарную модель рассматриваемой системы.
Выводы
Разработаны математическая модель тепловла-гопереноса и методика численного анализа величины тепловых потерь при эксплуатации теплотру-бопровода в условиях увлажнения изоляции влагой, содержащейся в воздухе канального пространства.
Показано, что применение гидроизоляции позволяет снизить до 9 % тепловые потери канального теплопровода, работающего в условиях контакта с влажным воздухом.
Установлено, что в рамках предложенной модели период выхода процесса влагонасыщения на стационарный режим может быть исключен из рассмотрения.
Работа выполнена при финансовой поддержке РФФИ (грант № 08-08-00143-а).
СПИСОК ЛИТЕРАТУРЫ
1. Байбаков С.А., Тимошкин А.С. Основные направления повышения эффективности тепловых сетей // Электрические станции. - 2004. - № 7. - С. 19-25.
2. Извеков А.В., Коновальцев С.И., Данилов О.Л., Новиков А.В., Альбертинский Л.И. Потери тепла в вентилируемых каналах тепловых сетей // Теплоэнергетика. - 1994. - № 12. - С. 37-42.
3. Чудновский А.Ф. Теплофизические характеристики дисперсных материалов. - М.: Физматгиз, 1962. - 456 с.
4. Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло- и массообмена. - М.: Наука, 1984. - 288 с.
5. Соколов Е.Я. Теплофикация и тепловые сети. - М.: МЭИ, 1999. - 472 с.
6. Уонг Х. Основные формулы и данные по теплообмену для инженеров. - М.: Атомиздат, 1979. - 216 с.
7. Каст В., Кришер О., Райнике Г., Винтермантель К. Конвективный тепло- и массоперенос. - М.: Энергия, 1980. - 49 с.
8. Петров-Денисов В.Г., Масленников Л.А. Процессы тепло- и влагообмена в промышленной изоляции. - М.: Энергоатомиз-дат, 1983. - 192 с.
9. СНиП 41-03-2003. Тепловая изоляция оборудования и трубопроводов. - М.: Издательство стандартов, 2004. - 28 с.
Поступила 05.09.2008г.
