Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

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

Максимов В.И. 1 Нагорнова Т.А. 1
1 Национальный исследовательский Томский политехнический университет
В статье рассмотрены возможные варианты расчета турбулентного режима в геометрически простой воздушной области. Представлены результаты математического моделирования свободно-конвективного теплопереноса в газовой среде, окруженной со всех сторон ограждающими конструкциями. На нижней границе раздела сред учитывался источник тепловыделения (движение нагреваемого воздуха осуществлялось только за счет подъемной силы). Решена система нестационарных дифференциальных уравнений Навье-Стокса в приближении Буссинеска для газа и теплопроводности для твердых стенок с соответствующими начальными и граничными условиями. Процесс конвективного теплопереноса рассчитан с использованием моделей турбулентности Прандтля и Прандтля-Райхарда. Учтены процессы теплообмена рассматриваемой области с окружающей средой. Получены распределения средних температур воздуха от времени. Проведен анализ безразмерного коэффициента теплоотдачи на внутренней границе раздела сред «воздух – твердая стенка».
конвективный теплоперенос
турбулентность
источник локального нагрева
математическое моделирование
средняя температура
уравнения Навье-Стокса
1. Белов И.А., Исаев С.А. Моделирование турбулентных течений: Учебное пособие. – СПб.: Балтийский государственный технический университет, 2001. – 108 с.
2. Лапин Ю.В., Стрелец М.Х. Внутренние течения газовых смесей. – М.: Наука, 1989. – 370 с.
3. Мовчан В.Т., Шквар Е.А. Математическое моделирование пограничных слоев // Прикладная гидромеханика. – 2005. – Т. 79, № 7. – С. 73–85.
4. Соловьев С.В., Хаенко Л.С. Моделирование естественной конвекции в наклонной полости // Вестник ТОГУ. – Физико-математические науки. – 2010. – Т. 16, № 1. – С. 27–34.
5. Ткаченко Г.В., Урюков Б.А. Модель турбулентной естественной конвекции около вертикальной стенки на основе гипотезы С.С. Кутателадзе // РНКТ-4: Труды Четвертой Российской национальой конференции по теплообмену (Москва, 2006). – М., 2006. – Т. 3. – С. 160–163.
6. Шлихнинг Г. Теория пограничного слоя. – М.: Наука, 1974. – 712 с.
7. Heat transfer under heating of a local region of a large production area by gas infrared radiators / Kuznetsov G.V., Kurilenko N.I., Maksimov V.I., Mamontov G.Ya., Nagornova T.A. // Journal of Engineering Physics and Thermophysics. – 2013. – V. 86, № 3. – P. 519–524.
8. Kuznetsov G.V., Sheremet M.A. Conjugate natural convection in an enclosure with local heat sources // Computational Thermal Sciences. – 2009. – Is. 1. – № 3. – Р. 341–360.
9. Maksimov V.I., Nagornova T.A. Influence of heatsink from upper boundary on the industrial premises thermal conditions at gas infrared emitter operation // EPJ Web of Conferences. – 2014. Vol. 76. – Article number 01006.
10. Schmidt W. Turbulente Ausbreitung eines Stormes erhitzter Luft. – ZAMM, 1941. – V. 21, № 5/6. – P. 265–271.

В последнее время все большее внимание уделяется изучению процессов естественной конвекции [4]. Исследование конвективного теплопереноса позволяет проводить выбор оптимальных режимов работы в любой области и соответствующих технических решений. Но при моделировании конвективных течений в областях больших размеров достаточно часто приходится учитывать турбулентность [5]. Одним из важных вопросов при решении задачи движения воздуха или жидкости является модель турбулентности [3]. В монографии Шлихтинга Г. [6] говорится об обоснованности применения полуэмпирических формул для расчета свободной турбулентности (турбулентные течения без ограничивающих стенок) [10]. Поэтому вопрос о возможных вариантах моделей конвективного теплопереноса в замкнутых объемах при локальном нагреве весьма актуален.

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

macs1.tif

Рис. 1. Область решения задачи: 1 – воздух; 2 – ограждающие конструкции; L, Н – длина и ширина области решения; Индексы: L, R, UP, LB – левая, правая, верхняя и нижняя границы раздела сред

Постановка задачи

Для физической постановки задачи выбрана достаточно простая по геометрии замкнутая прямоугольная область, включающая воздух (1) и ограждающие конструкции (2). На трех внешних границах приняты условия теплоизоляции, на левой внешней границе X = 0 реализованы условия теплообмена с окружающей средой. По всему периметру границы раздела сред заданы граничные условия четвертого рода (равенство тепловых потоков и температур) за исключением участка на нижней границе, где находится источник тепловыделения (L1 < X < L2).

Математическая постановка задачи

Численное моделирование естественной конвекции в замкнутой области с учетом турбулентности проведено двумя способами. Постановка задачи включала систему дифференциальных уравнений Навье-Стокса для газа и теплопроводности для твердых стенок. Турбулентная составляющая учитывалась реализацией алгебраических моделей Прандтля-Райхарда и Прандтля [1]. Моделирование было проведено в системе MathLab.

Математическая постановка принималась аналогичной [7, 9] и имела вид:

macs01.wmf, (1)

macs02.wmf, (2)

macs03.wmf, (3)

macs04.wmf. (4)

Начальные условия для системы уравнений (1–4)

macs05.wmf macs06.wmf (5)

Граничные условия на трех внешних границах области решения

macs07.wmf. (6)

На левой внешней границе

macs08.wmf. (7)

На границах раздела ограждающих конструкций и газа, выполняются условия:

macs09.wmf, macs10.wmf, (8)

macs11.wmf.

macs12.wmf, macs13.wmf, (9)

где X, Y – безразмерные декартовы координаты; Θ – безразмерная температура; Ω – безразмерный аналог вихря скорости; macs14.wmf – число Фурье; Gr = gβL3(Tit – T0)v2 – число Грасгофа; g – ускорение, создаваемое массовыми силами, м/с2; а – коэффициент температуропроводности, м2/с; β – температурный коэффициент объемного расширения, К-1; Bi = αL/λ – число Био; α – коэффициент теплоотдачи между внешней средой и рассматриваемой областью решения, Вт/(м2∙К); Re = U∙L/vt – число Рейнольдса; Pr = v/a – число Прандтля; T – температура, К; T0 – температура в начальный момент времени, К; Tit – масштаб температуры, К; V0 – масштаб скорости (скорость конвекции), м2/с; U, V – безразмерные скорости; t – время, с; t0 – масштаб времени, с; λ – коэффициент теплопроводности твердой стенки, Вт/(м⋅К); λ1,2 – относительный коэффициент теплопроводности; νt – коэффициент турбулентной кинематической вязкости, м2/с; τ – безразмерное время; ψ – безразмерный аналог функции тока.

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

macs19.wmf

где Umax и Umin – максимальная и минимальная величины скорости в слое, δ – полуширина слоя смешения, Χ – эмпирический безразмерный параметр, постоянный по толщине слоя, х – расстояние, измеренное в направлении потока.

Модели Прандтля [1] соответствует определение коэффициента турбулентной вязкости по формуле:

macs20.wmf, macs21.wmf,

где lm – путь смешения, x, y – координаты, м; υ – компонента скорости, м/с; k – универсальный коэффициент пропорциональности, не зависящий от числа Рейнольдса.

Результаты исследования и их обсуждение

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

Численные исследования проведены при значениях температур: начальной – Θ0 = 0, источника тепловыделения – Θit = 1, окружающей среды – Θe = 0.

macs2a.tif macs2b.tif

а) б)

macs2c.tif macs2d.tif

в) г)

Рис. 2. Поля температур (a, в) и изолинии функции тока (б, г) при τ = 500 для модели Прандтля-Райхарда (а, б) и модели Прандтля (в, г)

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

При малых временах, например (τ= 50), теплоотвод с левой внешней границы еще не оказывает заметного воздействия и в воздушной области формируется симметричное температурное поле. Распределение температуры (рис. 2) наглядно иллюстрирует влияние подъемной силы, которая обусловлена неоднородностью поля температур. Анализ средних по газовой области температур позволяет говорить об их существенно неравномерном распределении во времени из-за влияния архимедовой силы (таблица).

На рис. 4–5 приведено изменение среднего числа Нуссельта по вертикальным (X = LL, X = LR) границам раздела сред во времени, определенного по выражениям [9]:

macs22.wmf;

macs23.wmf.

При увеличении временного параметра становится заметным вихреобразование. Основной поток нагретого воздуха смещается к правой границе раздела сред, в верхней части которой происходит интенсивный теплообмен между нагретым газом и холодной ограждающей конструкцией. Следствием этого является (рис. 4) увеличение среднего числа Нуссельта на границе X = LR.

Температурные поля, полученные при реализации двух разных подходов к моделированию турбулентности достаточно сопоставимы, направления движения нагретого воздуха совпадают.

Зависимость средней температуры воздуха от времени: 1 – модель Прандтля-Райхарда, 2 – модель Прандтля

t, c

50

100

150

200

250

300

400

450

500

модель Прандтля-Райхарда

278,7

278,8

278,57

278,52

278,66

279,5

279,19

279,17

278,51

модель Прандтля

277,33

278,53

276,91

276,89

276,9

277,23

277,44

277,6

277,64

macs3.tif

Рис. 3. Зависимость среднего числа Нуссельта от времени для границы X = LL: 1 – модель Прандтля-Райхарда, 2 – модель Прандтля

macs4.tif

Рис. 4. Зависимость среднего числа Нуссельта от времени для границы X = LR: 1 – модель Прандтля-Райхарда, 2 – модель Прандтля

Можно отметить, что полученные результаты иллюстрируют стабилизацию режимов теплообмена на границах раздела «воздух – твердая стенка» при относительно больших временах (t = 500), за исключением правой границы раздела «воздух – ограждающая конструкция». При этом средние числа Нуссельта с ростом времени растут по мере увеличения t более чем в 2 раза. Результаты численного моделирования показывают возможность применения разных подходов к описанию тепловых режимов турбулентной конвекции в замкнутых областях с локальными источниками тепловыделения.

Заключение

Проведено теоретическое исследование турбулентной естественной конвекции в системе с источником тепловыделения. Полученные результаты позволяют оценить характер течения и интенсивность анализируемого процесса теплопереноса, а также новую информацию, которая не только характеризует конвективный режим течения, но и является дополнительной для построения и апробации белее сложных моделей турбулентности, учитывающих особенности конвективного теплопереноса в больших замкнутых объемах при локальном нагреве. Результаты математического моделирования позволяют сделать выводы о возможности применения алгебраической модели Прандтля для рассмотрения турбулентных течений при естественной конвекции в геометрически простых областях [2]. Применение данной модели целесообразно в случае, когда движение нагреваемого источником тепловыделения воздуха осуществляется за счет подъемной силы [10]. Следует отметить, что применение алгебраической модели при расчете свободно-конвективных течений в турбулентном режиме значительно проще и требует меньшего времени вычислений, чем более сложные модели турбулентности.

Работа выполнена в рамках НИР Госзадания «Наука» (Шифр ФЦП 2.1321.2014).

Рецензенты:

Борисов Б.В., д.ф.-м.н., профессор кафедры «Теоретическая и промышленная теплотехника» НИ ТПУ, г. Томск;

Голдаев С.В., д.ф.-м.н., старший научный сотрудник, профессор кафедры «Теоретическая и промышленная теплотехника» НИ ТПУ, г. Томск.

Работа поступила в редакцию 29.12.2014.


Библиографическая ссылка

Максимов В.И., Нагорнова Т.А. ПРИМЕНЕНИЕ ПРОСТЫХ МОДЕЛЕЙ ТУРБУЛЕНТНОГО КОНВЕКТИВНОГО ТЕПЛОПЕРЕНОСА В ЗАМКНУТЫХ ОБЪЕМАХ ПРИ ЛОКАЛЬНОМ НАГРЕВЕ // Фундаментальные исследования. – 2014. – № 12-10. – С. 2108-2113;
URL: http://fundamental-research.ru/ru/article/view?id=36534 (дата обращения: 14.12.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074