ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ СОПРЯЖЕННОГО СВОБОДНО-КОНВЕКТИВНОГО ТЕПЛООБМЕНА В ЗАМКНУТОЙ ПОЛОСТИ С ТРЕУГОЛЬНЫМ ЭЛЕМЕНТОМ ОБЪЕМНОГО ТЕПЛОВЫДЕЛЕНИЯ Текст научной статьи по специальности «Физика»

ВЕСТНИК ПЕРМСКОГО УНИВЕРСИТЕТА

2018 • ФИЗИКА • Вып. 4 (42)

УДК 536.2

PACS 02.60.Cb; 44.25.+f

Численное исследование сопряженного свободно - конвективного теплообмена в замкнутой полости с треугольным элементом объемного тепловыделения

Н. С. Гибановa, М. А. Шереметb

Томский государственный университет 634050, Томск, пр. Ленина, 36 a email: [email protected] b email: [email protected]

Выполнено численное исследование нестационарных ламинарных режимов сопряженного свободно-конвективного теплопереноса в замкнутой квадратной полости с тепловыделяющим элементом треугольной формы, расположенным в центре нижней границы рассматриваемой полости. Вертикальные стенки полости считались изотермическими и выполняли функцию теплоотвода. Процессы переноса импульса, массы и энергии в рассматриваемой замкнутой области описывались системой дифференциальных уравнений в частных производных с использованием преобразованных переменных «функция тока - завихренность». Краевая задача была решена методом конечных разностей на структурированной сетке. Особое внимание уделялось анализу эволюции изолиний функции тока и температуры, характеризующих особенности гидродинамики и теплопереноса. Выполнено также исследование влияния степени черноты поверхности локального источника тепловыделения и внутренних границ замкнутой полости на гидродинамику и теплоперенос. Показаны особенности эволюции температурного факела, проанализированы распределения локальных и интегральных характеристик для различных чисел Рэлея.

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

Поступила в редакцию 29.10.2018; принята к опубликованию 19.11.2018

Numerical investigation of conjugate free convective heat transfer inside an enclosure with triangular heat generating source

N. S. Gibanova, M. A. Sheremetb

Tomsk State University, Lenin Avenue 36, 634050, Tomsk a email: [email protected] b email: [email protected]

Numerical investigation of unsteady laminar conjugate natural convection in a square enclosure with a triangular heat generating source located in a center of bottom wall has been carried out. Vertical walls of the cavity were isothermal and they can be considered as heat removal system. The process of heat, mass and momentum transfer inside the enclosure has been described by the system of partial differential equations using the modified variables "stream function - vorticity". The boundary-value problem has been solved by the finite difference method using the structured grid.

© Гибанов Н.С., Шеремет М.А., 2018

L@ ®

распространяется на условиях лицензии J Creative Commons Attribution 4.0 International (CC BY 4.0).

The special focus has been given to analysis of evolution of streamlines and isotherms, which describe the behavior of hydrodynamics and heat transfer. Also the study of impact of the surface emissivity of heater walls and cavity walls has been performed. The features of unsteady thermal plume, distributions of both local and integral characteristics at different Ra numbers have been shown.

Keywords: natural convection; surface radiation; laminar flow; triangular shape of energy source; finite difference method.

Received 29.10.2018; accepted 19.11.2018 doi: 10.17072/1994-3598-2018-4-52-59

1.Введение

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

Процессы переноса массы и энергии вследствие естественной конвекции в замкнутой полости исследованы достаточно широко. Имеется множество различных работ по данной тематике. В работе [1] проводились исследования турбулентной естественной конвекции в двумерном, трехмерном и квазидвумерном приближениях. Авторы оценивали возможность применения данных двумерной или квазидвумерной задач для описания пространственных процессов. В результате было отмечено, что двумерные расчеты слабо сопоставимы с трехмерными, а данные квазидвумерного приближения качественно воспроизводят зависимость числа Нуссельта от толщины полости, но значения чисел Нуссельта меньше на 30%. В [2] проведено численное исследование влияния пористого слоя на естественную конвекцию наножид-кости в открытой полости с горячей вертикальной стенкой. Авторы отмечают возможность интенсификации теплообмена с ростом концентрации на-ночастиц при условии близости пористого слоя и греющей стенки. Экспериментальное и численное исследование режимов естественной конвекции воздуха внутри наклонной полости проведено в [3]. В качестве рабочей области была рассмотрена замкнутая квадратная полость с боковыми изотермическими стенками-нагревателями. В середине нижней стенки располагался охладитель, установ-

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

Относительно сопряженных задач конвективно-радиационного теплопереноса можно отметить, что такие работы являются наиболее интересными вследствие взаимодействия различных видов передачи энергии. В настоящее время имеются исследования, посвященные сопряженной турбулентной естественной конвекции и поверхностному излучению [5, 6], а также сложному теплообмену в случае ламинарных режимов в двумерном и трехмерном приближении с локальными источниками энергии [7-9].

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

2. Математическая модель и метод решения

Область решения рассматриваемой задачи представляет собой замкнутую квадратную полость с вертикальными изотермическими стенками постоянной минимальной температуры Тс (рис.1). Горизонтальные стенки считались адиабатическими. В центральной части нижней стенки рассматриваемой полости находился локальный источник объемного тепловыделения Q треугольной формы.

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

дО дТ ЗО ЗТ ЗО

дт ЗУ дХ дХ ЗУ '

уд2 О д2 5© дХ2 + дУ2 I+ дХ'

З2Т З2Т -Т +-Т = "О,

дХ ЗУ2

д© д¥д&_д¥д& дт дУ дХ дХ дУ

1

з2© З2© +

^зх 2 зу

з©

~дг

а а

/

у/Рг-Яа

З2© З2© ^ —7 + —т + О

ЗХ2 ЗУ2

Для определения плотности радиационного потока использовался подход, основанный на вычислении плотности потока эффективного излучения, заключающийся в решении следующих уравнений [12]:

вгаЛ Як ^ ^к-!Я! ■

1=1

N

Як = (1 -бк )Е рк-Л

г=1

(2.5)

(2.6)

Рис. 1. Область исследования

Для описания процессов переноса массы, импульса и тепловой энергии используется система нестационарных двумерных уравнений Навье-Стокса, записанная в приближении Буссинеска [10]. После введения преобразованных переменных, таких как функция тока у (и = ду/ду, V = -ду/дх) и завихренность ю

(ю = ¿V/ Зх -ди/ Зу), а также обезразмеривания система уравнений примет следующий вид [11]:

(2.1)

(2.2)

(2.3)

Для описания процесса теплопереноса в источнике тепловыделения используется следующее уравнение теплопроводности:

(2.4)

+ 8к (1 -д)4 |©к + 0.51+|

Для учета радиационной составляющей в процессе теплообмена между поверхностями важно учитывать всё падающее на поверхность и отраженное от нее излучение. Поверхности считаются диффузно-серыми и разбиваются на N малых участков, для которых формулируется условие теплового баланса. Отраженное и собственное излучение объединяются в одно эффективное излучение Як. Важную роль в распределении тепла играет геометрия системы. Для оценки этого параметра вычисляются угловые коэффициенты излучения. Угловой коэффициент ¥ы определяет долю энергии излучения, которая передается с к-й поверхности на /-ю. Плотность потока результирующего излучения определяется как разность эффективного излучения участка "к" и суммы долей потоков эффективного излучения всех участков системы, достигающих данной к-й поверхности.

В уравнениях (2.1)-(2.6) использовались следующие параметры: X, У - безразмерные декартовы координаты; т - безразмерное время; © - безразмерная температура; Т - безразмерная функция тока; и, V - безразмерные составляющие вектора скорости в проекциях на оси X, У, соответственно; О - безразмерная завихренность скорости; а/ - коэффициент температуропроводности среды, заполняющей полость [м2/с]; а„ - коэффициент температуропроводности материала источника энергии [м2/с]; Qrad - безразмерная плотность радиационного потока; ¥ы - угловые коэффициенты; гк - степень черноты поверхности стенок и источника тепла; д = Т0 / Т - температурный параметр; Яа = gPATI3/- число Рэлея; V - коэффициент кинематической вязкости [м2/с]; в -температурный коэффициент объемного расширения [К-1]; g - ускорение свободного падения [м/с2]; Ь - характерный размер полости [м]; АТ - перепад температур [К]; Рг = v|а^ - число Прандтля;

Os = QI}/()^АТ) - число Остроградского; -коэффициент теплопроводности материала источника энергии [Вт/(м-К)]; Q - плотность объемного тепловыделения источника энергии [Вт/м3].

Граничные условия для представленной задачи

+

(2.1)-(2.6) имели следующий вид:

- на вертикальных стенках x = 0, X = 1:

дТ

Т = 0, — = 0, 0 = 0;

дХ

- на горизонтальных стенках У = 0, У = 1:

dW d®

W = 0, — = 0, — = ±NmdQmd ; dY dY

на поверхности источника энергии:

dW

W = 0, — = 0, dn

А д® д® f

Xf dn

dn

Nrad Qra

где Nrad - радиационно-кондуктивный параметр; X/ - коэффициент теплопроводности внутренней среды [Вт/(м-К)].

В начальный момент времени предполагалось, что среда внутри рассматриваемой области находится в состоянии покоя (Т = О = 0), начальная температура внутри области считалась равной 0.5. В следующий момент времени мгновенно задавались представленные выше граничные условия.

т—'—I—1—I—1—г

Рис. 2. Влияние сеточных параметров на изолинии функции тока W, изолинии температуры ®, среднего конвективного (Nucon) и среднего полного (Nutotai) чисел Нуссельта

Представленная краевая задача (2.1)-(2.6) решалась с помощью метода конечных разностей на равномерной сетке. Для аппроксимации конвективных слагаемых использовалась монотонная аппроксимация Самарского [13], для диффузионных слагаемых - центральные разности. Уравнения параболического типа (2.1), (2.3) и (2.4) решались на основе локально-одномерной схемы А.А. Самар-

ского, позволяющей перейти от двумерной задачи к двум одномерным задачам. Полученная в результате аппроксимации система линейных алгебраических уравнений решалась методом прогонки. Разностные уравнения эллиптического типа (2.2), (2.5) и (2.6) были решены методом последовательной верхней релаксации. Для определения угловых коэффициентов ¥ы был использован метод натянутых нитей Хоттеля [14, 15].

В целях оптимизации вычислений были проведены тестовые расчеты при Яа = 105, О5 = 1, Рг = 0.7, е = 0.6 на нескольких сетках с различным количеством узлов. В результате самой оптимальной была выбрана сетка размерностью 150x100. На рис. 2 приведены сравнения локальных и интегральных параметров на разных стеках. Можно заметить, что отличия выбранной сетки от наиболее подробной сетки являются незначительными.

3. Результаты численного моделирования

Представленная краевая задача была численно решена при следующих значениях безразмерных параметров: 0 < е < 0.9; 0 < т < 6000; Яа = 105; Рг = 0.7; X, = 150 Вт/(м-К); X/ = 0.026 Вт/(м-К); Nгad = 77.87. Особое внимание было уделено анализу влияния фактора нестационарности, а также приведенной степени черноты на структуру течения и теплообмен в полости.

На рис. 3 изображены изолинии функции тока Т и изотермы 0 для различных значений степени черноты е.

Можно отметить, что в рассматриваемой полости формируется конвективное течение симметричное относительно центра области вследствие центрального расположения источника тепловыделения. В верхней части области зарождаются ядра конвективных ячеек, ориентированные под углом, соответствующим наклону стенок нагревательного элемента. Как видно из рисунка, изотермы достаточно наглядно отражают распределение тепловой энергии внутри области. При отсутствии излучения (рис 3, а) область прогревается наиболее интенсивно. Можно проследить характер распределения энергии в полости, когда тепловой факел, отходя от нагревателя, достигает верхней адиабатической стенки и начинается более интенсивное распределение тепла в горизонтальном направлении. Далее происходит взаимодействие с холодным фронтом, исходящим от боковых охлаждающих стенок, и массы прогретого воздуха начинают опускаться вниз вдоль холодных стенок. При введении дополнительного механизма переноса тепла - излучения температурный рост внутри полости снижается при постепенном увеличении степени черноты (рис. 3, б-г). Такие изменения можно объяснить способностью стенок поглощать и отдавать тепло более интенсивно.

Рис. 3. Изолинии функции тока Т и изотермы 0 при Яа = 105, т= 6000: а - е = 0, б - е = 0.3, в - е = 0.6, г -, е= 0.9

Рис. 4. Временная зависимость средней температуры в области 0ср от приведенной степени черноты стенок

Рис. 5. Временная зависимость средней температуры внутри источника энергии 0ср от приведенной степени черноты стенок

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

На рис. 6 представлена эволюция изолиний функции тока Т и температуры 0. Как можно заметить из рисунка, в начальный момент времени (рис. 6, а) происходит кондуктивный нагрев зоны вблизи локального источника энергии. Около стенок нагревателя начинают формироваться конвективные ячейки. Структура вихрей симметрична относительно центральной вертикальной оси, и циркуляция среды внутри этих зон отражает зарождение восходящих течений около нагретых поверхностей источника энергии. При дальнейшем прогреве области (рис. 6, б) форма вихрей, а также структура ядер конвективных течений изменяются. Над нагревателем начинает формироваться тепловой факел. В момент времени т = 10 (рис. 6, в) вихревое течение заполняет большую часть полости, а ядра конвективных ячеек смещаются в верхнюю ее часть. Толщина теплового факела уменьшается с ростом вертикальной зоны проникновения. При дальнейшем увеличении времени (рис. 6, в-е) происходит значительный прогрев области решения. Ядра конвективных ячеек вытягиваются, приобретая сначала овальную форму, а затем смещаются к верхней части полости и располагаются под некоторым углом к вертикальной оси. Структура теплового факела также модифицируется, достигая верхней границы области, тепловые потоки начинают распределяться вдоль верхней стенки и затем опускаются в нижнюю часть полости.

ч» 0

Рис. 6. Изолинии функции тока Т и температуры © при Яа = 105, е = 0.6 в различные моменты времени: а - т = 1, б - т = 3, в -х = 10, г - т = 100, д - т = 1000, е - т = 6000

4.Заключение

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

Было отмечено влияние излучения на характеристики теплообмена в полости и в источнике энергии. Изменение степени черноты полости не только позволяет управлять температурными режимами, но и, главное, уменьшать температуру внутри тепловыделяющего элемента.

Показана динамика развития полей функции тока и температуры. Установлены зависимости формы теплового факела и структуры течения от времени.

Исследование выполнено при финансовой поддержке гранта Президента для молодых ученых МД-2819.2017.8.

Список литературы

1. Теймуразов А. С. Турбулентный конвективный теплоперенос в тонких вертикальных слоях жидкости // Вестник Пермского университета. Физика. 2016. Вып. 1 (32). С. 18-24.

DOI: 10.17072/1994-3598-2016-1-18-24.

2. Miroshnichenko I. V., Sheremet M. A., Oztop H. F. Abu-Hamdeh N. Natural convection of alumina-water nanofluid in an open cavity having multiple porous layers // International Journal of Heat and Mass Transfer. 2018. Vol. 125. P. 648-657.

DOI: 10.1016/j.ijheatmasstransfer.2018.04.108

3. Leporini M., Corvaro F., Marchetti B., Polo-nara F., Benucci M. Experimental and numerical investigation of natural convection in tilted square cavity filled with air // Experimental Thermal and Fluid Science. 2018. Vol. 99. P. 572-583.

DOI: 10.1016/j.expthermflusci.2018.08.023

4. Михайленко С. А., Шеремет М. А. Моделирование конвективного теплопереноса во вращающейся замкнутой полости с локальным источником энергии // Вестник Пермского университета. Физика. 2017. Вып. 1 (35). С. 1925. DOI: 10.17072/1994-3598-2017-1-19-25

5. Мирошниченко И. В., Пахомов М. А., Шеремет М. А. Численный анализ турбулентных режимов сопряженного конвективно-радиационного теплопереноса в замкнутой области со стеклянной стенкой // Вестник Пермского университета. Физика. 2018. Вып. 1 (39). С. 17-25. DOI: 10.17072/1994-3598-2018-1-1725

6. Miroshnichenko I. V., Sheremet M. A. Radiation effect on conjugate turbulent natural convection in a cavity with a discrete heater // Applied Mathematics and Computation. 2018. Vol. 321. P. 358371. DOI: 10.1016/j.amc.2017.11.010

7. Kuznetsov G. V., Sheremet M. A. Conjugate natural convection with radiation in an enclosure // International Journal of Heat and Mass Transfer. 2009. Vol. 52. P. 2215-2223.

DOI: 10.1016/j.ijheatmasstransfer.2008.12.006

8. Singh D. K, Singh S. N. Conjugate free convection with surface radiation in open top cavity // International Journal of Heat and Mass Transfer. 2015. Vol. 89. P. 444-453.

DOI: 10.1016/j.ijheatmasstransfer.2015.05.038

9. Patil S., Sharma A. K., Velusamy K. Conjugate laminar natural convection and surface radiation in enclosures: Effects of protrusion shape and position // International Journal of Heat and Mass Transfer. 2016. Vol. 76. P. 139-146.

DOI: 10.1016/j. icheatmasstransfer.2016.05.006

10. Gibanov N. S., Sheremet M. A. Natural convection in a cubical cavity with different heat source configurations // Thermal Science and Engineering Progress. 2018. Vol. 7. P. 138-145.

DOI: 10.1016/j.tsep.2018.06.004

11. Кузнецов Г. В., Шеремет М. А. Численное моделирование температурных полей узлов и блоков радиоэлектронной аппаратуры и электронной техники // Микроэлектроника. 2009. Т. 38. № 5. С. 344-352.

12. Martyushev S. G., Sheremet M. A. Numerical analysis of conjugate natural convection and surface radiation in an enclosure with local heat source // Computational Thermal Sciences. 2013. V. 5. P. 11-25.

DOI: 10.1615/ComputThermalScien.2012006040

13. Алешкова И. А., Шеремет М. А. Математическое моделирование сопряженной термогравитационной конвекции в пористой среде // Вестник Удмуртского университета. Математика. Механика. Компьютерные науки. 2010. Вып. 2. С. 49-56.

14. Siegel R., Howell J. R. Thermal radiation heat transfer. London: Taylor and Francis, 2002. 868 p.

15. Спэрроу Э. М., Сесс Р. Д. Теплообмен излучением. Л.: Энергия, 1971. 296 с.

References

1. Teimurazov A. S. Turbulent convective heat transfer in thin vertical fluid layers. Bulletin of Perm University. Physics, 2016, no. 1 (32), pp. 18-24 (In Russian). DOI: 10.17072/1994-3598-2016-118-24

2. Miroshnichenko I. V., Sheremet M. A., Oztop H. F. Abu-Hamdeh N. Natural convection of alumina-water nanofluid in an open cavity having multiple porous layers. International Journal

of Heat and Mass Transfer, 2018, vol. 125, pp. 648-657.

DOI: 10.1016/j.ijheatmasstransfer.2018.04.108

3. Leporini M., Corvaro F., Marchetti B., Polo-nara F., Benucci M. Experimental and numerical investigation of natural convection in tilted square cavity filled with air. Experimental Thermal and Fluid Science, 2018, vol. 99, pp. 572-583.

DOI: 10.1016/j.expthermflusci.2018.08.023

4. 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;. DOI: 10.17072/1994-3598-2017-1-19-25

5. 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 (In Russian). DOI: 10.17072/1994-3598-2018-1-17-25

6. Miroshnichenko I. V., Sheremet M. A. Radiation effect on conjugate turbulent natural convection in a cavity with a discrete heater. Applied Mathematics and Computation, 2018. vol. 321, pp. 358-371. DOI: 10.1016/j.amc.2017.11.010

7. Kuznetsov G. V., Sheremet M. A. Conjugate natural convection with radiation in an enclosure. International. Journal of Heat Mass Transfer, 2009, vol. 52, pp. 2215-2223.

DOI: 10.1016/j.ij heatmasstransfer.2008.12.006

8. Singh D. K, Singh S. N. Conjugate free convection with surface radiation in open top cavity. International Journal of Heat and Mass Transfer, 2015, vol. 89, pp. 444-453.

DOI: 10.1016/j.ijheatmasstransfer.2015.05.038

9. Patil S., Sharma A. K., Velusamy K. Conjugate laminar natural convection and surface radiation in enclosures: Effects of protrusion shape and position. International Journal of Heat and Mass Transfer, 2016, vol. 76, pp. 139-146.

DOI: 10.1016/j.icheatmasstransfer.2016.05.006

10. Gibanov N. S., Sheremet M. A. Natural convection in a cubical cavity with different heat source configurations. Thermal Science and Engineering Progress, 2018, vol. 7, pp. 138-145.

DOI: 10.1016/j.tsep.2018.06.004

11. Kuznetsov G. V., Sheremet M. A. Chislennoe modelirovanie temperaturnih polej uzlov i blokov radioelektronnoi apparaturi i elektronnoi tehniki (Numerical simulation of temperature fields of units and blocks of radio-electronic equipment and electronic equipment). Microelectronics, 2009, vol. 38, pp. 344-352 (In Russian).

12. Martyushev S.G., Sheremet M.A. Numerical analysis of conjugate natural convection and surface radiation in an enclosure with local heat source.

Computational Thermal Sciences, 2013, vol. 5, pp. 11-25.

DOI: 10.1615/ComputThermalScien.2012006040 13. Aleshkova I. A., Sheremet M. A. Mathematical simulation of conjugate natural convection in a porous medium. Vestnik Udmurtskogo Universiteta

Matematika. Mekhanika. Komputernie Nauki, vol. 2010, no. 2, pp. 49-56 (In Russian).

14. Siegel R., Howell J. R. Thermal radiation heat transfer. London: Taylor and Francis, 2002. 868 p.

15. Sparrow E.M, Cess R. D. Radiation heat transfer. Brooks Pub. Co., 1966, 322 p.

Просьба ссылаться на эту статью в русскоязычных источниках следующим образом: Гибанов Н. С., Шеремет М. А. Численное исследование сопряженного свободно-конвективного теплообмена в замкнутой полости с треугольным элементом объемного тепловыделения // Вестник Пермского университета. Физика. 2018. № 4 (42). С. 52-59. doi: 10.17072/1994-3598-2018-4-52-59

Please cite this article in English as:

Gibanov N. S., Sheremet M. A. Numerical investigation of conjugate free convective heat transfer inside an enclosure with triangular heat generating source. Bulletin of Perm University. Physics, 2018, no. 4 (42), pp. 52-59. doi: 10.17072/1994-3598-2018-4-52-59