ТЕПЛОПЕРЕНОС В ЗАМКНУТЫХ ВОЗДУШНЫХ ПОЛОСТЯХ СТРОИТЕЛЬНЫХ МАТЕРИАЛОВ И ИЗДЕЛИЙ

Юркевич А. А.

Кандидат технических наук, доцент, Ижевский государственный технический университет имени М.Т.Калашникова

ТЕПЛОПЕРЕНОС В ЗАМКНУТЫХ ВОЗДУШНЫХ ПОЛОСТЯХ СТРОИТЕЛЬНЫХ МАТЕРИАЛОВ И ИЗДЕЛИЙ

Аннотация

Авторами приводятся результаты численного эксперимента по определению составляющих сложного теплообмена (теплопроводность, конвекция и излучение), протекающем в полостях и порах неоднородных строительных  материалах и изделиях. Сложный теплообмен в полости и поре определяется на основе совместного решения задачи конвективного теплообмена в объемной постановке и теплообмена излучением в замкнутой системе.  

Ключевые слова: сложный теплообмен, теплопроводность, конвекция, излучение, пустотные строительные материалы и изделия.

Yurkevich A.A.

PhD of Technical Sciences, assosiate professor, Izhevsk State Technical University named MT Kalashnikov
HEAT TRANSFER IN CLOSED AIR CAVITY CONSTRUCTION MATERIALS AND PRODUCTS

Abstract

The authors present the results of the numerical experiment to determine the components of the complex heat transfer (conduction, convection and radiation) flowing into the cavities and pores of non-homogeneous materials and products. Complex heat transfer in the cavity and the pore depends on the basis of the joint solution of the problem of convective heat transfer in a three-dimensional setting and radiation heat transfer in a closed system.

Keywords: complex heat transfer, conduction, convection, radiation, hollow building materials and products.

Введение. Одним из способов повышения термического сопротивления штучных строительных материалов (кирпичи, камни керамические, шлакоблоки и т. п.), является увеличение их пустотности - соотношения между объемом занимаемым воздухом и твердой массой. Пустотность увеличивается за счет изготовления пустот цилиндрической или прямоугольной формы (в основном вертикальными). Количество которых варьируется от 2 до 60. Полости замкнутые, так как с четырех сторон ограничены стенками изделия, а с двух (верхней и нижней) - цементно-песчаным раствором.

Большое количество пустот практически исключает возможность решения сопряженной задачи конвективного теплообмена с целью оптимизации размеров пустот, их формы, количества и относительного расположения, так как количество одновременно решаемых задач равно числу пустот. Поэтому задача оптимизации формируется как краевая задача теплопроводности многосвязной области с заменой процесса конвективного теплообмена в полостях, сопровождающегося переносом теплоты излучением, равновеликой (эффективной) по тепловому балансу процессом теплопроводности [4, 5], который определяется по уравнению:

Физико-математическая модель и постановка задачи. Плотность теплового потока  определяется в результате решения задачи конвективного теплообмена в полости или поре:

где   

Рис. 1. Схема для расчета конвективного теплообмена  в полости или поре

Термогравитационная конвекция несжимаемой вязкой жидкости в приближении Буссинеска описывается следующей системой уравнений:

• уравнения вихря

• уравнение энергии, описывающее температурное поле в движущейся  жидкости

 

•  уравнения составляющих вектора скорости   

• уравнения векторного потенциала (уравнения Пуассона)

• граничные условия для уравнений вихря

при X = 0 и X = H в плоскостях перпендикулярных оси X

при Y = 0 и Y = L  в плоскостях перпендикулярных оси Y 

при Z = 0 и Z = M в плоскостях перпендикулярных оси Z

• граничные условия для уравнений Пуассона

при X = 0 и X = H в плоскостях перпендикулярных оси X

при Y = 0 и Y = L  в плоскостях перпендикулярных оси Y 

при Z = 0 и Z = M в плоскостях перпендикулярных оси Z

• на всех поверхностях поры или полости  (условие прилипания);
• граничные условия для уравнения энергии

в плоскостях перпендикулярных оси 

в плоскостях перпендикулярных оси 

в плоскостях перпендикулярных оси Z , при Z = 0 и Z = M  q = 0 .

Плотность потока результирующего излучения   на тепловоспринимающую поверхность, определяется по формуле

Здесь: R - плотность потока эффективного излучения поверхности; F -  угловой коэффициент излучения; i,j - индексы поверхности. Для расчетной схемы на рис. 1 N = 6, то есть 1 £ i, j £ 6.

Плотность потока эффективного излучения поверхности определяется по уравнению

здесь  - постоянная Стефана-Больцмана, (5,6699´10^-8 Вт/(м² К⁴));  - степень черноты i-ой поверхности; - абсолютная температура i-ой поверхности. Средние температуры i-х поверхностей вычислялись по формулам

где  - температура определяемая в результате решения задачи конвективного теплообмена.

Метод решения и результаты численного моделирования. Решение задачи сложного теплообмена находилось численно в такой последовательности.  В начальный момент времени газ находится в состоянии покоя, распределение температуры в полости соответствовало среднему значению между «холодной» и «горячей» поверхностей.

Стационарное решение достигается методом установления. Трехмерные уравнения вихря, три уравнения Пуассона  и уравнение энергии с соответствующими граничными условиями, численно реализованы методом конечных разностей по неявной схеме  переменных направлений (схема стабилизационной поправки) с разностями против потока. Исходная система уравнений считается решенной, когда значения рассматриваемых функций в расчетных точках не меняется в течении 10-100 шагов по времени и максимальная относительная погрешность  становилась меньше 10^-3. Были проведены численные исследования по сходимости результатов к стационарному  решению на равномерной сетке с разным количеством узлов и  Анализ результатов расчетов показал, что равномерная сетка с количеством узлов обеспечивает требуемую точность.

Адекватность учета свободной конвекции была проверена сравнением результатов численного расчета с опытными данными Михеева М. А. [1] и Эккерта Э. Р. [2] для эффективного коэффициента теплопроводности вертикальных воздушных прослоек. Среднее расхождение соответственно составляет 5,4 % и 3,6%. Максимальное расхождение не превышает 10 % и 6,6%.

В результате решения задачи конвективного теплообмена были получены значения температур на поверхностях полости или поры. После этого, для задачи теплообмена излучением находились средние значения температуры для каждой поверхности. Задача теплообмена излучением в замкнутой системе численно решена зональным методом [3].

Тестирование численной реализации задачи теплообмена излучением в замкнутой системе проводилась на корректных контрольных примерах.

В воздушных полостях и порах при сложном теплообмене передача теплоты осуществляется тремя способами: теплопроводностью, конвекцией и излучением. Для выяснения доли каждой из этих составляющих в сложном теплообмене в газовых полостях и порах были проведены численные исследования. Рассматривались кубические полости и поры с размерами от 0,001 до 0,1 м, результаты представлены в виде графиков на рис. 2 для разных значений перепадов температур . Степень черноты всех поверхностей   (большинство строительных материалов имеют данную степень черноты).

Рис. 2. Зависимости составляющих сложного теплообмена от размеров воздушных полостей и пор при =20°С.

Рис.3. Зависимости составляющих сложного теплообмена от размеров воздушных полостей и пор при =10°С.

Рис. 4. Зависимости составляющих сложного теплообмена от размеров воздушных полостей и пор при  =5°С.

Рис. 5. Зависимости составляющих сложного теплообмена от размеров воздушных полостей и пор при  =3°С.

Рис. 6. Зависимости составляющих сложного теплообмена от размеров воздушных полостей и пор при  =1°С.

По результатам численных исследований (рис. 3-7) сделаны следующие выводы:

1. для полостей и пор малых размеров (H<0,002 м основная доля теплоты передается теплопроводностью (до 90 %);
2. для полостей и пор с размерами H»0,01 метра доли теплоты передаваемой излучением и теплопроводностью становятся одинаковыми и составляют   50 %  каждая;
3. в полостях и порах с размерами более 0,01 метра доля теплоты передаваемая излучением, постепенно увеличивается и при H»0,05 метра, составляет 60 - 80 %;
4. доля теплоты, передаваемая конвекцией, в полостях и порах с размером H»0,05 метра составляет в зависимости от разности температур 10 - 25 % от общего количества теплоты.

Обозначения. - коэффициента теплопроводности газа; - составляющая  , учитывающая влияние конвекции; -составляющая ,учитывающая влияние излучения; - составляющая , учитывающая теплопроводность и конвекцию в полости или поре; - плотность теплового потока за счет теплопроводности и конвекции;   - плотность потока результирующего излучения на тепловоспринимающую (охлажденной) поверхность; - температура на нагреваемой и охлажденной поверхности газовой полости или поры; Q_К - количество теплоты предоваемого конвекцией; Q_ИЗЛ - количество теплоты предоваемого излучением; Q - полное (суммарное) количество теплоты предоваемого теплопроводностью, конвекцией и излучением; L, M, H - эквивалентная высота, ширина и глубина газовой полости или поры; Здесь - составляющие скорости по координатам  - составляющие  вихря по координатам  - составляющие соленоидального векторного потенциала по координатам  - коэффициент кинематической вязкости, м²/с;  - ускорение свободного падения, м/с²; - коэффициент температурного расширения газа, 1/К;  - коэффициент температуропроводности, м²/с; t - температура, определяемая на основе численного решения, °С;  - плотность потока результирующего излучения для тепловоспринимающей поверхности; R - плотность потока эффективного излучения поверхности; F - диффузионные угловые коэффициенты излучения; N- количество поверхностей;  - постоянная Стефана-Больцмана, (5,6699´10^-8 Вт/(м² К⁴));  - степень черноты i-ой поверхности;  - абсолютная температура поверхности,  - средняя температура поверхности; - температура определяемая в результате решения задачи конвективного теплообмена; - количество узловых точек, соответственно по оси 

Литература

1. Михеев М. А.  Основы теплопередачи  - М.-Л.: Государственное энергетическое издательство, 1949. - С. 83-87
2. Эккерт Э. Р., Дрейк Р. М. Теория тепло- и массообмена. - М.: Госэнергоиздат, 1961. -
3. Оцисик М. Н. Сложный теплообмен. М.: Мир, 1976. -  С. 171-182
4. Юркевич А. А., Диденко В. Н., Корепанов Е. В. Метод прогнозирования коэффициента теплопроводности неоднородных (пористых, пустотных) строительных материалов с учетом конвективного теплообмена и излучения //Вестник ИжГТУ: Периодический научно-теоретический журнал Ижевского государственного технического университета. -     Вып. 2. - Ижевск: Изд-во ИжГТУ, 1998. - С. 15 - 18.
5. Юркевич А.А. Метод прогнозирования теплоизоляционных свойств строительных материалов и изделий: дис. канд.техн.наук. – Ижевск, 1999. – 125 с.