ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА
2018
• ФИЗИКА •
Вып. 1 (39)
УДК 536.2
PACS 02.60.Cb; 44.25 +
Численный анализ турбулентных режимов сопряженного конвективно - радиационного
a Томский государственный университет 634050, Томск, пр. Ленина, 36 email: [email protected] b Институт теплофизики СО РАН
630090,Новосибирск, проспект Академика Лаврентьева,! email: [email protected] c Томский государственный университет 634050, Томск, пр. Ленина, 36 email: [email protected]
Проведено математическое моделирование нестационарной турбулентной термогравитационной конвекции и теплового поверхностного излучения в замкнутой полости с твердыми теплопроводными стенками конечной толщины при наличии источника энергии постоянной температуры. Процесс переноса импульса, массы и энергии описывался системой дифференциальных уравнений в переменных «функция тока - температура - завихренность». Краевая задача решена методом конечных разностей на неравномерной сетке. В качестве замыкающей модели турбулентности использовалась стандартная к -s модель. Проанализировано влияние степени черноты ограждающих поверхностей и фактора нестационарности на режимы течения и теплоперенос. В результате проведенных исследований получены распределения как локальных (изолинии функции тока и температуры), так и интегральных (средние радиационное и конвективное числа Нуссельта) характеристик. Показано, что увеличение степени черноты ограждающих поверхностей проявляется в изменении структуры течения и интенсификации теплопереноса.
Ключевые слова: естественная конвекция; поверхностное излучение; турбулентность; источник энергии; метод конечных разностей
Поступила в редакцию 03.12.2017; принята к опубликованию 11.01.2018
Numerical analysis of conjugate turbulent convective-radiative heat transfer in an enclosure with a glass wall
I. V. Miroshnichenkoa, M. A. Pakhomovb, M. A. Sheremetc
a Tomsk State University, Lenin Avenue 36, 634050, Tomsk email: [email protected]
b Institute of Thermophysics SB RAS, 630090, Novosibirsk, Academician Lavrentyev Avenue, 1 email: [email protected]
c Tomsk State University, Lenin Avenue 36, 634050, Tomsk email: [email protected]
теплопереноса в замкнутой
со стеклянной стенкой
И. В. Мирошниченкоa, М. А. Пахомовb, М. А. Шереметc
c
© Мирошниченко И.В., Пахомов М.А., Шеремет М.А., 2018
17
Mathematical simulation of unsteady turbulent natural convection and thermal surface radiation in a closed cavity with solid heat-conducting walls of finite thickness in the presence of a heat source of constant temperature has been carried out. The governing equations for the conservation of mass, momentum, and energy have been written using a system of partial differential equations in terms of stream function, temperature and vorticity. The boundary-value problem has been solved by the finite difference method. The effect of turbulence has been modeled using the standard k -s model. The effects of the surface emissivity and the nonstationary factor on the fluid flow and heat transfer have been studied. As a result of the studies, distributions of both local (isolines of the stream function and temperature) and integral (average radiative and convective Nusselt numbers) characteristics have been obtained. It has been shown that an increase in the surface emissivity leads to a modification of flow patterns and the heat transfer enhancement.
Keywords: natural convection; surface radiation; turbulence; energy source; finite difference method.
Received 03.12.2017; accepted 11.01.2018 doi: 10.17072/1994-3598-2018-1-17-25
1. Введение
Исследования совместных эффектов естественной конвекции и поверхностного излучения, ввиду первостепенной важности в энергетике и строительстве, бурно развиваются в наше время [1-3]. Интерес исследователей к численному моделированию данных процессов обусловлен, в первую очередь, возможностью определить реальное тепловое состояние исследуемых объектов, в том числе в системах-потребителях тепловой энергии. Более того, стоящая довольно остро у нас в стране проблема теплоэнергосбережения мотивирует продолжать исследования в области конвективно-радиационного теплопереноса в целях повышения энергоэффективности объектов промышленного или жилого назначения.
К настоящему времени проведено большое количество численных и экспериментальных исследований ламинарных режимов естественной конвекции в замкнутых областях с источниками энергии [4-6]. Так, например, влияние размеров и формы локального тепловыделяющего элемента на теплопередачу в квадратной полости детально проанализировано в [7]. Сопряженный конвективно-радиационный теплообмен в замкнутом кубе численно исследован в [8] методом контрольного объема. Установлено, что с ростом коэффициента теплопроводности твердых стенок полости существенно увеличивается роль как радиационного, так и конвективного механизмов теплообмена на поверхности источника энергии. Однако во множестве прикладных задач рассматриваются режимы течения, которые характеризуются достаточно большими скоростями движения среды в полости, либо сама полость имеет большие размеры. В таких задачах естественным образом возникает необходимость моделирования турбулентных режимов переноса массы, импульса и энергии. К наиболее высокоточным методам расчета турбулентных течений относятся DNS (прямое числен-
ное моделирование) и LES (метод моделирования крупных вихрей). Оба метода требуют существенных вычислительных ресурсов, что несколько ограничивает применение данных подходов.
Самым распространенным методом к исследованию подобного типа задач на сегодняшний день является использование RANS подхода, заключающегося в применении осредненных по Рейнольд-су уравнений Навье-Стокса, которые замыкаются с помощью различных полуэмпирических моделей турбулентности. В работе [9] проведен детальный обзор существующих экспериментальных и численных научных работ в области турбулентной естественной конвекции в замкнутых областях. Показано, что исследований, посвященных численному анализу сложного теплообмена (конвекция, кондукция и излучение) в областях различной геометрии не так много. Связано это, прежде всего, с математическими и вычислительными сложностями в описании совместных эффектов рассматриваемых процессов передачи тепла.
В настоящей работе излагаются результаты математического моделирования нестационарных турбулентных режимов сопряженной термогравитационной конвекции и теплового поверхностного излучения в замкнутой области с твердыми стенками, состоящими из материалов с разными теп-лофизическими свойствами, при наличии источника энергии. Основной акцент сделан на изучении процесса эволюции во времени теплового факела и гидродинамических структур, а также на исследовании влияния приведенной степени черноты ограждающих конструкций на режимы течения и теплообмена в рассматриваемой области.
2. Математическая модель и метод решения
Область решения представляет собой замкнутую полость, изображенную на рис. 1. Внутри полости находится ньютоновская диатермичная сре-
да, которая удовлетворяет приближению Бус-синеска. Локальный источник энергии расположен вдоль нижней стенки. Теплофизические характеристики одной из твердых стенок (4 на рис. 1) соответствуют параметрам стекла (коэффициент теплопроводности '/. \ = 1.15 Вт/(м-К), приведенная степень черноты поверхности стенки е4 = 0.9) и не варьируются в задаче. На внешней границе левой стенки задаются условия конвективного теплообмена с окружающей средой, а на оставшихся внешних границах - условие теплоизоляции. При описании радиационной составляющей полного теплового потока использовались следующие допущения: рассматривалось только поверхностное излучение от источника тепловыделения и твердых стенок и не учитывалось излучение, поступающее извне.
Рис. 1. Область решения: 1 - твердые стенки, 2 - воздушная полость, 3 - источник энергии постоянной температуры, 4 - стеклянная стенка
Для сгущения расчетной сетки к твердым стенкам использовалось следующее алгебраическое преобразование координат [10]:
l = a +
1 t пк
А 1 + -
b - a b - a
X -
a + b
2
+ п
tg ук
) = a + -
i * пк
t 1 + tgv-
b - a b - a
Y -
a + b
(2.1)
2
п
tg у к
где к - параметр сжатия, а, Ь -начальная и, соответственно, конечная точки отрезка прямой с предполагаемым сгущением.
Процесс переноса импульса, массы и энергии для данной задачи описывается системой нестационарных уравнений Навье-Стокса. Краевая задача сформулирована в безразмерных переменных «функция тока - завихренность - температура» с использованием преобразования вида (2.1). Система безразмерных уравнений имеет следующий вид:
dl ( d^X
dX2 { dX ) д% + d2) д¥ (d)^2 д2 dY2 д) [dY ) дrf ~ '
dd^Ju-dlX ) КдО+
дт [ dX д% ) dX д$
(у_ d^X ) 7сП = [ dY д) ) dY д)
(2.2)
d) дО. dY д)
dl ду, I dlY д у dX2 ~д1 + {dX ) 1
d2) ду ( d))2 д2у
dY2 д) \ dY) д)
х|О + 2 ) + 41 ) +
dY д) ) dX \dY) дф) д) dl д© dX д1'
д© тт dl д© тйпд© dl
— + U—--+ V —— = —— х
дт dX д1 dY д) dX
'51
У
1
v
л/Ra • Pr Prt
d) д dY д)
dl д© dX д1
.. \
(2.3)
(2.4)
л/Ra • Pr Prt
d) д© dY д)
дК TTdl дК тйпдК
— + U—--+ V—--=
дт dX д1 dY д)
dl дК dX д^
d) дК dY д)
(2.5)
+ P + G - E,
дЕ т^ЕдЕ TTdn дЕ
— + U—--+ V—-— =
дт dX д^ dY д)
(2.6)
Е
Е2
+ C1s (Pk + C3sGk ) Е - C2s —
К
К
+
2
+
д© ~дт
й4 д© йц
л/Яа • Рг
й2ц 5©
йУ2 дц ^ йУ ^ дц
й4
йХ
д2 © д#2
7 =
>и
2 д2 ©
(2.7)
ЦШ ]2 + 2Гц^2
дХ д4
дУ дц
дцди+ д¥
дУ дц дХ д4
дц д© РГ дУ дц
где р - слагаемое, характеризующее порождение турбулентности за счет сдвиговых напряжений; Gi описывает генерацию или диссипацию турбулентной кинетической энергии за счёт выталкивающей силы. Турбулентная вязкость вычислялась по формуле Колмогорова - Прандтля
К2
"'= ^ Т
Уравнения (2.5) и (2.6) - транспортные уравнения низкорейнольдсовой к -е модели турбулентности [10]. Безразмерная координата первого пристеночного узла
Яа гди <
Рг V дУ
Здесь X, У - безразмерные физические координаты; 4, ц -безразмерные координаты в вычислительной области; т - безразмерное время; © - безразмерная температура; Т - безразмерная функция тока (и = дТ/дУ, V = -дТ/дХ) ; и, V - безразмерные составляющие скорости в проекции на оси X, У, соответственно; О - безразмерная завихренность скорости (О = б V/дХ - ди/дУ) ; К - безразмерная кинетическая энергия турбулентности; Е -безразмерная скорость диссипации кинетической энергии турбулентности; у1 - коэффициент турбулентной вязкости; Рг = у1/а1 - турбулентное число Прандтля; аг - коэффициент турбулентной температуропроводности; Рг = Ц(^ра2) - число
Прандтля; ц - динамический коэффициент вязкости воздуха, кг/(м-с); р - плотность воздуха, кг/м3; а2 - коэффициент температуропроводности воздуха, м2/с; а,2 = а/а2 - относительный коэффициент температуропроводности; g - ускорение свободного падения, м/с2; Яа = -Те)Ь3/(ца2) -
число Рэлея; в - температурный коэффициент объемного расширения, 1/К; Ь - характерный размер полости, м (см. рис. 1).
Основываясь на гипотезах Буссинеска и градиентной диффузии имеем, соответственно, следующие соотношения:
Здесь у - расстояние от стенки до ближайшего уз*
ла, и - динамическая скорость, V - кинематический коэффициент вязкости.
Рекомендуемые эмпирические константы, используемые в данной модели, представлены в таблице 1.
Таблица 1. Эмпирические константы для к-е модели турбулентности
Сц С1е С2е С3е СТк сте Рг
0.09 1.44 1.92 0.8 1.0 1.3 1.0
Начальные условия для сформулированной системы уравнений (2.2)-(2.7) следующие:
Т (4, ц,0) = О (4, ц,0) = К (4,ц,0)= Е (4,ц,0) = 0, ©(4,ц,0) = 0.5.
Граничные условия имеют следующий вид:
д4 д©
• на границе 4 = 0:--= Ы •© (условие
дХ д4
конвективного теплообмена с окружающей средой);
• на границе 4 = 1+21/Ь: д©/д4 = 0 (условие теплоизоляции);
• на границе ц = 0: д©/дц= 0 (условие теплоизоляции);
• на границе ц = 1+2//Ь: с©/дц = 0 (условие теплоизоляции);
• на поверхности источника энергии: © = 1;
• на внутренних границах раздела воздух-стенка, параллельных оси 04:
Т = 0, К = 0,
дТ дЕ дТ = 0, — = 0, ©=©, дц дц
. дцд©^=цд^-м о ■
4,2 ^ ^ гайх-'Тай >
дУ дц дУ дц
• на внутренних границах раздела воздух-стенка, параллельных оси 0ц:
Т = 0, К = 0,
дТ дЕ дТ = 0, — = 0, ©.=©,, д4д4 1 2
. 54д©1 = 54д©1 -м 0
'■2 дХ д4 дХ д4 гай0а.
Здесь ^ = стТ£Т - Т )] - радиационно"
кондуктивный параметр; ст - постоянная Стефана-Больцмана, Вт/(м2-К4); Л1 - коэффициент теплопроводности твердых стенок, Вт/(м-К), Х2 - коэффициент теплопроводности воздуха, Вт/(м-К), Х4 -
а
V
2
коэффициент теплопроводности стеклянной стенки, Вт/(м-К), Д2 = Д/Д - относительный коэффициент теплопроводности (в зависимости от материала ограждающей стенки), Б1 = аЬ/Д - число Био материала твердой стенки; а - коэффициент теплообмена между внешней средой и рассматриваемой областью решения, - степень черноты стеклянной стенки, ё - степень черноты ограждающих твердых стенок.
Для нахождения безразмерной плотности радиационного потока 0 применялся метод решения с использованием плотности потока эффективного излучения. Данный метод заключается в разрешении двух разностных дифференциальных уравнений (2.8)-(2.9) методом последовательной верхней релаксации:
0гай ,к = Як ^ ^к-1Я1,
I =1
N
(2.8)
(2.9)
Здесь Як - безразмерная плотность эффективного излучения А-й поверхности; /•'/■, - угловой коэффициент между поверхностями "А" и ";'"; ёк -приведенная степень черноты к-й поверхности; С = Те¡Ты - температурный параметр; ©к - температура к-й поверхности; N - количество поверхностей.
Полученная система дифференциальных уравнений в частных производных (2.2)-(2.7) решена методом конечных разностей. Подробное описание используемых методов решения можно найти
в [10].
Для валидации разработанного метода решения проведено сравнение с экспериментальной работой [11]. На рис. 2 и 3 представлено сравнение профилей температуры и вертикальной компоненты скорости с экспериментом. Разработанный численный метод показал достаточно хорошее согласие с экспериментальными данными.
В таблице 2 представлены значения средних радиационного
N
Ии, = гай
0.9
21/Ь+0.9
I
21/Ь
0гай\ц=цк8 й4
и конвективного
Ии = ■
21/Ь+0.9
1 21! Ьь + 0.
0.9 2{ь
дц д© дУ дц
й4
чисел Нуссельта [12] в зависимости от размера расчетной сетки. В связи с тем, что уменьшение шага разностной сетки отражается на существен-
ном повышении времени счета, по результатам численного эксперимента выбрана сетка размерностью 120x120 (так как различия в значениях среднего конвективного числа Нуссельта для сеток 120x120 и 180x180 составило менее 5 %).
Рис. 2. Профиль температуры в сечении У = 0.5 при Яа = 1.58 • 109 в сравнении с экспериментальными данными [11]
Рис. 3. Профиль вертикальной компоненты скорости в сечении У = 0.5 при Яа = 1.58 • 109 в сравнении с экспериментальными данными [11]
Таблица 2. Средние конвективное и радиационное числа Нусселъта при Яа = 109, г = 6000 , ё = 0.3 в зависимости от размерности вычислительной сетки
Размер сетки разница, % разница, %
60x60 39.79 - 26.25 -
120x120 50.34 20.95 26.16 0.34
180x180 52.85 4.74 26.10 0.23
I=1
3. Результаты численного моделирования
Сформулированная краевая задача решалась численно при следующих значениях безразмерных комплексов: 0 < е < 0.9 ; 0 < т < 10000; Яа = 109; Рг = 0.7; ^ 2 = 6-8, Н-ай = 390.15. Основное внимание было сосредоточено на исследовании влияния приведенной степени черноты ё и фактора нестационарности на структуру течения и теплоперенос в рассматриваемой области.
Изолинии функции тока и температуры, отражающие влияние нестационарности процессов переноса тепла и массы, представлены на рис. 4. В начальный момент времени (рис. 4, а) в рассматриваемой полости формируются две конвективные ячейки, располагающиеся в области основания. Данные вихревые структуры отличаются, в первую очередь, направлением циркуляции газа, и их появление вызвано несколькими причинами:
радиационным теплообменом, постоянным подводом тепла от источника тепловыделения и воздействием окружающей среды посредством кондук-тивного механизма переноса энергии.
С течением времени (рис. 4, в) увеличивается масштаб правой конвективной ячейки (направление циркуляции которой соответствует движению по часовой стрелке) вследствие формирования в области над источником энергии теплового факела, смещенного ближе к правой стенке. Дальнейший рост времени (рис. 4, г - 4, ж) приводит к существенному прогреву внутренней части газовой полости и модификации изолиний функции тока - значительно увеличивается вихрь, расположенный в левой части полости (направление циркуляции которого соответствует движению против часовой стрелки).
Представленные распределения изолиний функции тока и температуры позволяют проследить за эволюцией термогидродинамических структур во времени, что, в конечном итоге, при-
Рис. 4. Изолинии функции тока Ч'и температуры ©при ё = 0.3 в различные моменты времени т. а-5, 6-15, в-30, г - 150, д - 500, е - 3000, ж-10000
водит к лучшему пониманию происходящих процессов переноса тепла и массы.
Данное увеличение связано, в первую очередь, с охлаждением на границе 4 = 0 вследствие воздействия низкой температуры окружающей среды. Также следует отметить более интенсивное проникновение фронта пониженной температуры извне через стеклянную стенку. Это обусловлено её более высокими теплопроводящими свойствами по сравнению с остальными ограждающими стенками.
Рис. 5. Изолинии функции тока ¥ и температуры ©при г = 10000: а - ё = 0, б — е = 0.9
На рис. 5 представлены распределения изолиний функции тока и температуры в зависимости от степени черноты ограждающих поверхностей. Следует отметить, что случай ё = 0 подразумевает рассмотрение только конвективного теплоперено-са без поверхностного излучения.
Как уже было отмечено выше, в рассматриваемой области формируются две конвективные ячейки, одна из которых (меньшей интенсивности) располагается в правом верхнем углу полости. Увеличение приведенной степени черноты поверхностей приводит к существенной трансформации вихревых структур внутри области решения. Эта трансформация обусловлена, в первую очередь, увеличением размеров основного вихря за счет пропорционального уменьшения вихря, расположенного в правом верхнем углу. Все это приводит к росту интенсивности течения. Подтверждением данного факта являются максимальные значения функции тока в ядрах конвективных
1-1-Г
0.7 0,8 0,9 ©
Рис. 6. Профиль температуры © при т = 10000 в сечении X = 0.6 в зависимости от приведенной степени черноты
Более того, учет излучательной способности поверхностей отражается в существенном повышении не только средней температуры в области решения (см. рис. 5), но и локальной температуры в верхней части полости (см. рис. 6). В частности, увеличение степени черноты ограждающих конструкций с 0 до 0.9 приводит к увеличению средней температуры во всей области более чем на 62 %.
Влияние ё на среднее радиационное число Нуссельта при Ка = 109 показано на рис. 7. Видно, что увеличение приведенной степени черноты поверхностей проявляется в интенсификации теплообмена на границе раздела двух сред. Рост ё от 0.3 до 0.9 приводит к увеличению среднего радиационного числа Нуссельта более чем в 2.4 раза.
ячеек
I = 0.029 <
I тах|^=0 I 1Г
= 0.035.
0 2000 4000 6000 8000
Рис. 7. Зависимость среднего радиационного числа Нуссельта на поверхности источника энергии от времени и приведенной степени черноты
4. Заключение
Проведено математическое моделирование турбулентной естественной конвекции и теплового поверхностного излучения в замкнутой полости с теплопроводными стенками конечной толщины при наличии источника энергии. В результате проведенных исследований получены различные распределения интегральных и локальных параметров, характеризующих течение и теплообмен в рассматриваемой области. Моделирование проводилось при следующих значениях безразмерных комплексов: 0 < е < 0.9, 0 < т < 10000, Ra = 109, Рг = 0.7, =6.8, N„,¿= 390.15.
Показано, что с течением времени в результате нагрева от источника энергии и охлаждения от левой теплопроводной стенки в области формируются две вихревые структуры разной интенсивности. Помимо этого, в работе продемонстрирована эволюция теплового факела.
Установлено, что повышение скорости движения и температуры воздушных масс в области решения происходит вследствие увеличения степени черноты поверхностей ограждающих конструкций.
Исследование выполнено при финансовой поддержке РФФИ в рамках научного проекта № 1738-50057 «мол_нр».
Список литературы
1. Wu T., Lei C. On numerical modelling of conjugate turbulent natural convection and radiation in a differentially heated cavity // Int. J. Heat Mass Transfer. 2015. Vol. 91. P. 454-466.
DOI: 10.1016/j.ij heatmas stransfer.2015.07.113.
2. Варапаев В. Н., Голубев С. С. Сравнение точного и приближенного учета теплового излучения границ при численном решении сопряженной задачи естественной конвекции в вертикальном слое оконных ограждений // Вестник МГСУ. 2011. № 8. С. 129-136.
3. Wang Y., Menga X., Yang X., Liua J. Influence of convection and radiation on the thermal environment in an industrial building with buoyancy-driven natural ventilation // Energy and Buildings. 2014. Vol. 75. P. 394-401.
DOI: 10.1016/j.enbuild.2014.02.031.
4. АстанинаМ. С., Шеремет М. А. Моделирование термогравитационной конвекции с переменной вязкостью в замкнутой полости с локальным источником энергии // Вестник Пермского Университета. Серия: Физика. 2015. Вып. 3 (31). С. 52-58.
5. Durgam S., Venkateshan S. P., Sundararajan T. Experimental and numerical investigations on optimal distribution of heat source array under natural and forced convection in a horizontal channel // International Journal of Thermal Sciences. 2017. Vol. 115. P. 125-138.
DOI: 10.1016/j.ijthermalsci.2017.01.017.
6. Михайленко С. А., Шеремет М. А. Моделирование конвективного теплопереноса во вращающейся замкнутой полости с локальным источником энергии // Вестник Пермского Университета. Физика. 2017. Вып. 1 (35). С. 19-25.
7. Гибанов Н. С., Шеремет М. А. Влияние формы и размеров локального источника энергии на режимы конвективного теплопереноса в квадратной полости // Компьютерные исследования и моделирование. 2015. Т. 7, № 2. С. 271-280.
8. Yang G., Wu J. Y. Effects of natural convection, wall thermal conduction, and thermal radiation on heat transfer uniformity at a heated plate located at the bottom of a three-dimensional rectangular enclosure // Numerical Heat Transfer, Part A: Applications. 2016. Vol. 69. P. 589-606.
DOI: 10.1080/10407782.2015.1090238.
9. Miroshnichenko I. V., Sheremet M. A. Turbulent natural convection heat transfer in rectangular enclosures using experimental and numerical approaches: A review // Renewable and Sustainable Energy Reviews. 2018. Vol. 82. P. 40-59.
DOI: 10.1016/j.rser.2017.09.005.
10. Miroshnichenko I. V., Sheremet M. A. Turbulent natural convection combined with thermal surface radiation inside an inclined cavity having local heater // International Journal of Thermal Sciences. 2018. Vol. 124. P. 122-130.
DOI: 10.1016/j.ijthermalsci.2017.09.023.
11. Ampofo F., Karayiannis T. G. Experimental benchmark data for turbulent natural convection in an air filled square cavity // Int. J. Heat Mass Transfer. 2003. Vol. 46. P. 3551-3572.
DOI: 10.1016/S0017-9310(03)00147-9.
12. Siegel R., Howell J. R. Thermal radiation heat transfer. London: Taylor & Francis, 2002. 868 p.
References
1. Wu T., Lei C. On numerical modelling of conjugate turbulent natural convection and radiation in a differentially heated cavity. International Journal of Heat and Mass Transfer, 2015, vol. 91, pp. 454-466.
DOI: 10.1016/j .ijheatmasstransfer.2015.07.113.
2. Varapaev V. N., Golubev S. S. Comparison of the exact and approximate calculation boundary heat radiation in numerical solution of conjugate problem natural convection in window vertical slot // Vestnik MGSU, 2011, no. 8, pp. 129-136 (In Russian).
3. Wang Y., Menga X., Yang X., Liua J. Influence of convection and radiation on the thermal environment in an industrial building with buoyancy-driven natural ventilation. Energy and Buildings, 2014, vol. 75, pp. 394-401.
DOI: 10.1016/j.enbuild.2014.02.031.
4. Astanina M. S., Sheremet M. A. Simulation of natural convection with variable viscosity in an enclosure with a local heat source // Bulletin of Perm University. Series: Physics, 2015, no. 3 (31), pp. 52-58 (In Russian).
5. Durgam S., Venkateshan S. P., Sundararajan T. Experimental and numerical investigations on optimal distribution of heat source array under natural and forced convection in a horizontal channel. International Journal of Thermal Sciences, 2017, Vol. 115, pp. 125-138.
DOI: 10.1016/j.ijthermalsci.2017.01.017.
6. Mihajlenko S. A., Sheremet M. A. Simulation of convective heat transfer inside a rotating enclosure with a local heat source. Bulletin of Perm University. Physics, 2017, vol. 1 (35), pp. 19-25 (In Russian).
7. Gibanov N. S., Sheremet M. A. Effect of shape and sizes of a local heat source on convective heat transfer in a square cavity. Computer Research and Modeling, 2015, vol. 7, no 2, pp. 271-280 (In Russian).
8. Yang G., Wu J. Y. Effects of natural convection, wall thermal conduction, and thermal radiation on
heat transfer uniformity at a heated plate located at the bottom of a three-dimensional rectangular enclosure. Numerical Heat Transfer, Part A: Applications, 2016, vol. 69, pp. 589-606. DOI: 10.1080/10407782.2015.1090238.
9. Miroshnichenko I.V., Sheremet M.A. Turbulent natural convection heat transfer in rectangular enclosures using experimental and numerical approaches: A review. Renewable and Sustainable Energy Reviews, 2018, vol. 82, pp. 40-59.
DOI: 10.1016/j.rser.2017.09.005.
10. Miroshnichenko I.V., Sheremet M.A. Turbulent natural convection combined with thermal surface radiation inside an inclined cavity having local heater. International Journal of Thermal Sciences, 2018, vol. 124, pp. 122-130.
DOI: 10.1016/j.ijthermalsci.2017.09.023.
11. Ampofo F., Karayiannis T.G. Experimental benchmark data for turbulent natural convection in an air filled square cavity. International Journal of Heat and Mass Transfer, 2003, vol. 46, pp. 35513572. DOI: 10.1016/S0017-9310(03)00147-9.
12. Siegel R., Howell J.R. Thermal radiation heat transfer. London: Taylor & Francis, 2002. 868 p.
Просьба ссылаться на эту статью в русскоязычных источниках следующим образом: Мирошниченко И. В., Пахомов М. А., Шеремет М. А. Численный анализ турбулентных режимов сопряженного конвективно-радиационного теплопереноса в замкнутой области со стеклянной стенкой // Вестник Пермского университета. Физика. 2018. № 1 (39). С. 17-25. doi: 10.17072/1994-3598-2018-1-17-25
Please cite this article in English as:
Miroshnichenko I. V., Pakhomov M. A., Sheremet M. A. Numerical analysis of conjugate turbulent convective-radiative heat transfer in an enclosure with a glass wall // Bulletin of Perm University. Physics, 2018, no. 1 (39), pp. 17-25. doi: 10.17072/1994-3598-2018-1-17-25