Исследование турбулентного теплопереноса в замкнутой прямоугольной области с теплопроводными ограждающими конструкциями в условиях лучистого нагрева внутренних границ Текст научной статьи по специальности «Физика»

УДК 536.33:536.244

ИССЛЕДОВАНИЕ ТУРБУЛЕНТНОГО ТЕПЛОПЕРЕНОСА В ЗАМКНУТОЙ ПРЯМОУГОЛЬНОЙ ОБЛАСТИ С ТЕПЛОПРОВОДНЫМИ ОГРАЖДАЮЩИМИ КОНСТРУКЦИЯМИ В УСЛОВИЯХ ЛУЧИСТОГО НАГРЕВА ВНУТРЕННИХ ГРАНИЦ

Г.В. КУЗНЕЦОВ, А.Э. НИ

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

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

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

Введение

В условиях непрерывного истощения ископаемых топлив все более актуальными становятся задачи энергосбережения. Одно из перспективных направлений генерации тепловой и электрической энергии - использование нетрадиционных энергетических ресурсов (солнечные батареи и коллекторы, тепловые насосы, ветроустановки и др.). Но в силу значительных материальных затрат использование альтернативных источников энергии оказывается во многих случаях неэффективным.

Высоким потенциалом энергосбережения обладают промышленные предприятия с крупногабаритными производственными цехами. Установлено [1, 2], что использование газовых инфракрасных излучателей (ГИИ) для локального обогрева помещений позволяет экономить до 40% энергетических ресурсов по сравнению с традиционным конвективным отоплением. По этой причине целесообразно исследование тепловых режимов областей с лучистыми системами теплоснабжения.

© Г.В. Кузнецов, А.Э Ни Проблемы энергетики, 2015, № 7-8

Известные математические модели теплопереноса в условиях радиационного нагрева объектов теплоснабжения, основанные на теории [3, 4], не описывают в полной мере совместно протекающие при работе ГИИ процессы кондукции, конвекции и излучения. Для решения таких сложных задач целесообразно использование аппарата [5, 6], разработанного для моделирования сопряженного теплопереноса в замкнутых областях с локальными источниками энергии в условиях естественной конвекции.

Исследование свободно-конвективного теплопереноса в отапливаемом с использованием ГИИ помещении выполнено в работах [7, 8]. Но численный анализ турбулентной термогравитационной конвекции [7] проведен без учета распределения лучистой энергии по внутренним поверхностям ограждающих конструкций (предполагалось, что вся энергия, поступающая от ГИИ, равномерно распределена только по нижней горизонтальной границе раздела «газ - стенка»). Четыре возможных варианта описания распределения радиационных потоков по внутренним поверхностям области анализа исследованы в работе [8], но в рамках ламинарной модели течения воздуха. В крупногабаритных же производственных помещениях преобладает турбулентная конвекция.

Целью работы является численное моделирование сопряженного турбулентного теплопереноса в замкнутой прямоугольной области, нагреваемой газовым инфракрасным излучателем, с учетом распределения лучистой энергии по закону Ламберта.

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

Рис. 1. Область решения задачи: 1 - воздух; 2 - ограждающие конструкции; 3 - газовый инфракрасный излучатель (условное обозначение)

Исследованы два возможных достаточно реальных варианта описания распределения лучистой энергии по внутренним поверхностям ограждающих конструкций:

1. Тепловой поток равномерно распределен только по нижней горизонтальной границе раздела «газ - стенка».

2. Распределение радиационной энергии задано по закону Ламберта [9]. Исследуемый процесс сопряженного турбулентного теплопереноса описывается

нестационарными уравнениями Рейнольдса и энергии для газа и уравнением теплопроводности для ограждающих конструкций. В безразмерных переменных «вихрь скорости О - функция тока ^ - температура © » эти уравнения имеют вид [10, 11]: © Проблемы энергетики, 2015, № 7-8

ао ттао тгао а2 —+и—+V— = —-аг ах ду дх2

о

У2¥ = -о,

а©1 + и а©1+Vа©1

а2

дг

ах

ау ах2

• Рг Ргг

а у2

©1

О

а©! "ах

а у2

л/Яа • Рг Рг

©

а©^

а^о-

'2

■ = У2©2.

(1) (2)

(3)

(4)

Турбулентная вязкость определялась с использованием модели Болдуина -Ломакса [12]:

2 I I

Vt = ^ , 1т = ку,

где / — путь смешения; у - координата; к - постоянная Кармана (к « 0,4). Начальные условия для уравнений (1) - (4):

¥(Х, У, 0) = 0; О(Х, У, 0) = 0; и(X, У, 0) = 0; V(X, У, 0) = 0;

©! (X, У, 0) = ©2 (X, У, 0) = 0. Граничные условия для уравнений (1) - (4): на внешних границах области решения:

X = 0, X = 1,2,0 <У < 1: а©2(X,У,Г) = 0;

ах

У = 0, У = 1,0<X < 1,2: а©2(X,У,г) = 0.

аУ

на внутренних границах «твердая стенка - воздух», параллельных оси ОХ:

а¥

¥ = 0,-= 0,

аУ

©,. =© ],

а©, ь] а©г-

= —---L + К1

где

аУ аУ

г = 1,2 ]=12

на внутренних границах «твердая стенка - воздух», параллельных оси ОУ:

а¥

¥ = 0, — = 0,

ах

а©/ ^ ] а©/ где

1 - - 1 + К1 ,

ах хг ах

г = 1,2 ] = 12

а • к) ^ „ Я-Р-(Ь — 2 • /1)3 • (Т1 — Т0)

где го =-0—- - число Фурье; Яа =-1-— - число Рэлея;

(Ь — 2 • /1)2 К1=д • (Ь — 2 • /1)

V • а

^ Т — Т0)

число Кирпичева; Рг = —

число Прандтля;

—

Рг? = — - турбулентное число Прандтля; а — коэффициент температуропроводности,

а

м2 /с ; а — турбулентный коэффициент температуропроводности; я — ускорение

свободного падения м/с2; р — коэффициент объёмного расширения, 1/К;

V — кинематическая вязкость, м2 /с; д - плотность теплового потока, Вт/м2; © Проблемы энергетики, 2015, № 7-8

а

X - коэффициент теплопроводности, Вт/(м • К); Х,У - безразмерные координаты, соответствующие ху; (Ь - 2/1) - характерный размер газовой полости по оси х, м; х-безразмерное время; и, V- безразмерные скорости

соответствующие пу; То - температура воздуха, и стенок в начальный момент времени, К; Тц - температура на поверхности излучателя, К;

Уравнения (1) - (4) с соответствующими начальными и граничными условиями решены методом конечных разностей, аналогично работам [5, 13]. Для оценки работоспособности используемого алгоритма решена краевая задача

термогравитационной конвекции в замкнутой квадратной полости [11, 14]. Результаты сравнения средних чисел Нуссельта приведены в табл. 1.

Таблица 1

Сравнение средних чисел Нуссельта_

Число Рэлея Ки ^ Ки ^ [11] Ки^ [14]

Яа = 108 33,92 30,32 31,2

Яа = 109 56,21 49,52 58,1

На основании сравнения полученных Кис результатами [11, 14] можно сделать обоснованный вывод о работоспособности используемого алгоритма решения. Результаты и обсуждения

Численное исследование проведено при значениях безразмерных критериев, соответствующих достаточно типичным условиям теплопереноса при работе ГИИ:

число Рэлея Яа = 10 , число Прандтля Рг=0,71, турбулентное число Прандтля Ргг =1,

число Кирпичева 0<Кл < 30. Результаты моделирования показаны на рис. 2.

г

и-

0.5

О (31 (36 11 12х

а)

г

б)

I ___ о оое--- Г..... д К 71 1 —_____---® Т"!

г;

Рис. 2. Изолинии температур (а, в) и функций тока (б, г) при т = 1000 : а, б) Равномерное распределение теплового потока по нижней горизонтальной границе; в, г) Распределение лучистой энергии и по вертикальным внутренним границам области в соответствии с законом Ламберта

Установлено, что в момент времени т = 1000 формируются вихревые (рис. 2, б, г) структуры в газовой полости. Причиной таких циркуляционных течений является непрерывный подвод лучистой энергии к внутренним поверхностям ограждающих конструкций. В первом случае (рис. 2, а, б) происходит прогрев газа у границы у=Ы, 11<х<14 и продвижение воздушных масс в направлении оси У. Во втором случае (при более равномерном нагреве внутренних границ и соответствующем отводем теплоты в ограждающие конструкции) циркуляционные течения принимают более выраженный характер (рис. 2, г).

С увеличением времени до т = 1500, в случае равномерного распределения теплового потока только по нижней горизонтальной границе (рис.3, а, б), формируется выраженное симметричное относительно оси Х=0,6 циркуляционное течение воздуха. При распределении лучистой энергии по закону Ламберта часть теплового потока поступает к вертикальным ограждающим конструкциям. Повышается температура в области Х=0,1; Х=1,1; 0,08<У<0,8 (рис.4) и происходит образование вторичных вихрей в газовой полости (рис. 3, г). Необходимо отметить, что такое перераспределение энергии за счет теплоотвода в вертикальные ограждающие конструкции ведет к понижению температуры по высоте области анализа в сечении Х=0,6 (рис. 4, б).

О 0.1

1.1 1.2

Рис. 3. Изолинии температур (а, в) и функций тока (б, г) при х = 1500 : (а, б) Равномерное распределение теплового потока по нижней горизонтальной границе; (в, г) Распределение лучистой энергии и по вертикальным внутренним границам области в соответствии с законом Ламберта

На рис. 4 показаны распределения температур по оси У в сечениях Х=0,12 и

Х=0,6.

У

О

У б

Рис. 4. Распределение температур по оси У при х = 1500 в сечении: а) Х=0,12; б) Х=0,6.

1 - равномерное распределение теплового потока по нижней горизонтальной границе;

2 - распределение лучистой энергии по закону Ламберта

Из рис. 4 видно, что при лучистом нагреве внутренних границ замкнутой прямоугольной области, заполненной газом и ограниченной твердыми теплопроводными стенками конечной толщины, формируются существенно неоднородные температурные поля. Распределение лучистой энергии по вертикальным ограждающим конструкциям приводит к увеличению градиента температур в сечении Х=0,6 (рис. 4, б). При этом растет температура воздуха у вертикальных стенок (рис. 4, а). Можно сделать вывод о значимости фактора распределения лучистых потоков по внутренним границам при анализе основных закономерностей теплопереноса в областях с радиационной системой отопления.

Проведена оценка влияния турбулизации течения на процесс теплопереноса в рассматриваемой области решения путем сравнения рассчитанных средних чисел Нуссельта с результатами работы [8]. На границе раздела «воздух - стенка» (у=М, /1<х</4) безразмерный коэффициент теплоотдачи определяется [5] как

Результаты для случая распределения лучистой энергии по закону Ламберта приведены в табл. 2.

Таблица 2

Сравнение средних чисел Нуссельта

Время, т Nu av ( Ra = 108 ) Nu^ [8] ( Ra = 106 )

500 2,99 1,11

1000 3,23 1,14

1500 3,32 1,15

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

Заключение

По результатам численного моделирования установлено, что процесс радиационного нагрева замкнутых прямоугольных помещений с теплопроводными и теплоаккумулирующими твердыми стенками носит существенно нестационарный характер. При этом формируются неоднородные температурные поля. Показано, что распределение тепловых потоков по внутренним поверхностям ограждающих конструкций приводит к увеличению градиента температур по высоте области анализа. Можно сделать вывод, что при постановке задач сопряженного теплопереноса с радиационными источниками энергии целесообразно учитывать распределение лучистой энергии и по вертикальным границам. При анализе тепловых режимов локальных рабочих зон в крупногабаритных помещениях с теплопроводными ограждающими конструкциями в условиях работы газовых инфракрасных излучателей целесообразно использовать модели естественной конвекции.

Работа выполнена в рамках НИР Госзадания «Наука» №13.1339.2014/К (шифр Федеральной целевой научно-технической программы 2.1410.2014).

Summary

Mathematical modeling of conductive - convective heat transfer in a closed rectangular area with thermally conductive enclosure structures in conditions of radiant heating of internal boundaries is passed. The fields of temperature and stream functions for the two possible descriptions of the radiant energy distribution (even distribution of heat flow only along the lower horizontal boundary and the Lambert's law radiant energy distribution along vertical inner surfaces) are obtained. It is established that the absorption of radiant fluxes by vertical enclosure structures radiant flux increases the temperature

gradients in the characteristic sections. Scale of the influence of turbulence on the flow rate of heat transfer is showed.

Keywords: conjugate heat transfer, numerical simulation, turbulence, thermogravitational convection, Baldwin — Lomax model, gas infrared emitter.

Литература

1. Романова Е. А. Энергосберегающие системы газового отопление и вентиляции // Новые технологии и изобретения. 2005. №2. С. 19 - 23.

2. Roth K. Infrared radiant heaters // ASHRAE. Journal. June. 2007.

3. Бухмиров В.В., Крупенников С.А., Солнышкова Ю.С. Алгоритм расчета систем лучистого отопления помещений // Вестник Ивановского государственного энергетического университета. 2010. № 4. С. 23 - 25.

4. Бухмиров В.В., Крупенников С.А., Солнышкова Ю.С. Модификация зонального метода для решения задач радиационного теплообмена: Основные положения // Вестник Ивановского государственного энергетического университета. 2009. № 2. С. 61 - 63.

5. Kuznetsov G.V., Sheremet M.A. Two-dimensional problem of natural convection in a rectangular domain with local heating and heat-conducting boundaries of finite thickness // Fluid Dynamics. 2006. V.41. P. 881 - 890.

6. Kuznetsov G.V., Sheremet M.A. Turbulant regime of thermogravitational convection in a closed cavity // Journal of Engineering Physics and Thermophysics. 2010. V. 83. №2. P. 346 - 357.

7. Кузнецов Г.В., Куриленко Н.И., Максимов В.И., и др. Свободно- конвективный теплоперенос в отапливаемых с использованием газовых инфракрасных излучателей производственных помещений // Известия высших учебных заведений. Проблемы энергетики. 2013. №1 - 2. С. 18 - 25.

8. Kuznetsov G.V., Nagornova T.A., Ni A.E. Computational modeling of conjugate heat transfer in a closed rectangular domain under the conditions of radiant heat supply to the horizontal and vertical surfaces of enclosure structures // Journal of Engineering Physics and Thermophysics. 2015. V.88. № 1. P. 168 - 177.

9. Теплопередача: Учебник для вузов / В.П. Исаченко, В.А. Осипова, А.С. Сукомел. 5-е изд., стереотипное. ООО «ТИД «Арис», 2014. 416 с., ил.

10. Лойцянский Л.Г. Механика жидкости и газа. М.: Наука, 1978. 736 с.

11. Sajjadi H., Kefayatib R. Lattice Bolzman simulation of turbulent natural convection in tall enclosure // Thermal Science. 2015. V.19. №1. P. 155 - 166.

12. Белов И.А., Исаев С.А. Моделирование турбулентных течений: Учебное пособие, СПб: Балтийский государственный технический университет, 2001. 108 с.

13. Kuznetsov G.V., Strizhak P.A. Heat and mass transfer at ignition of liquid fuel droplets spreading over the surface of massive hot bodies // Journal of Engineering Thermophysics. 2010. V.19. № 2. P. 75 - 84.

14. Dixit H.N., Babu V. Simulation of High Rayleigh Number Natural Convection in a Square Cavity using the Lattice Boltzmann Method // Int. J. Heat Mass Transfer. 2006. V. 49. P. 727-739.

Поступила в редакцию 29 апреля 2015 г.

Кузнецов Гений Владимирович - д-р физ.-мат. наук, профессор, заведующий кафедрой ТПТ ЭНИН Томского политехнического университета (ТПУ). Тел: 8(382)2564010. E-mail: [email protected].

Ни Александр Эдуардович - инженер кафедры ТПТ ЭНИН Томского политехнического университета (ТПУ). Тел:8(952)8802713. E-mail: [email protected].