ЧИСЛЕННЫЙ АНАЛИЗ СОПРЯЖЕННОГОКОНВЕКТИВНО-РАДИАЦИОННОГО ТЕПЛОПЕРЕНОСАВ ЗАМКНУТОЙ ОБЛАСТИ
Проведено математическое моделирование сложного теплопереноса в замкнутой прямоугольной области с теплопроводными стенками конечной толщины при наличии локального источника тепловыделения. В качестве механизмов переноса тепла в газовой полости рассматривались естественнаяконвекция и излучение, а в элементах ограждающей твердой оболочки - теплопроводность. Получены поля скорости и температуры, характеризующиеосновные закономерности исследуемого процесса в широком диапазоне изменения определяющих параметров.
Numerical analysis of conjugate convective-radiativeheat transfer in an enclosure.pdf Математическое моделирование сопряженного конвективно-радиационноготеплопереноса в газовой полости и теплопроводности в элементах твердого мате-риала имеет достаточно широкие области приложения [1, 2]. С точки зрения тео-ретического значения большой интерес вызывает наиболее полный анализ воз-можных режимов взаимодействия термогравитационной конвекции и излучения вгазовой полости. Исследование процессов теплопереноса в областях при наличиитеплопроводных стенок приобретает существенные связи с практикой (анализ те-плового состояния радиоэлектронной аппаратуры; исследование технологическихрежимов в машиностроении; поиск оптимальных параметров по созданию ком-фортных условий в жилых помещениях) [1 - 3].Целью настоящей работы является математическое моделирование нестацио-нарной сопряженной естественной конвекции в замкнутой области при наличиилокального источника тепловыделения в условиях радиационного теплопереносав газовой полости и неоднородного теплообмена с внешней средой.Постановка задачиРассматривается краевая задача нестационарного конвективно-радиационноготеплопереноса в прямоугольной области, представленной на рис 1. Область реше-ния представляет собой газовую полость, ограниченную теплопроводными стен-ками конечной толщины. На внутренней стороне левой стенки находится источ-ник тепловыделения, имеющий постоянную температуру в течение всего процес-са. Считается, что в начальный момент времени стенки полости и заполняющая еенесжимаемая жидкость имеют постоянную и одинаковую во всех точках темпера-туру, причем жидкость неподвижна. Горизонтальные стенки y = 0, y = 2h1+H ивертикальная стенка x = 2l1+L предполагаются теплоизолированными с наружной1 Работа выполнена в рамках реализации ФЦП «Научные и научно-педагогические кадры инновацион-ной России» на 2009 - 2013 годы (ГК № П357).Численный анализ сопряженного конвективно-радиационного теплопереноса 97стороны. На границе x = 0 осуществляется конвективно-радиационный теплооб-мен с окружающей средой. Газ, занимающий внутреннюю полость, считается вяз-кой, ньютоновской, несжимаемой, излучающей жидкостью, удовлетворяющейприближению Буссинеска.В такой постановке процесс переноса тепла в анализируемой области решения(рис. 1) описывается системой нестационарных двумерных уравнений конвекции вприближении Буссинеска в газо-вой полости [4], где слагаемое вуравнении энергии, характери-зующее излучение, описываетсяна основе приближения Россе-ланда [5] . Для анализа поля тем-пературы в твердых стенках ис-пользуется нестационарное дву-мерное уравнение теплопровод-ности [2] с нелинейными гранич-ными условиями. Основным вприближении Буссинеска являет-ся предположение о том, что рас-сматривается «слабая» конвек-ция: вызванные неоднородно-стью температуры отклоненияплотности от среднего значенияпредполагаются настолько ма-лыми, что ими можно пренеб-речь во всех уравнениях, кромеуравнения движения, где это от-клонение учитывается в члене с подъемной силой. В приближении Росселанда [5]излучающую среду можно рассматривать, как некоторый континуум фотонов, т.е.можно считать, что на каждый элемент среды, как и в случае молекулярной прово-димости, непосредственно влияют только его соседние элементы. В таких условияхперенос энергии излучения в среде можно уподобить диффузионному переносу.Если пренебрегать вязкой диссипацией энергии, то уравнения неразрывности,движения и энергии в газовой полости для рассматриваемой задачи будут иметьвидu v 0x y + = ; (1)2 22 2u u u v u p u ut x y x x y⎛ ⎞ ⎛ ⎞ ³⎜ + + ⎟ = − + ®⎜ + ⎟ ⎝ ⎠ ⎝ ⎠; (2)( )2 22 2 y г 0v u v v v p v v g T Tt x y y x y⎛ ⎞ ⎛ ⎞ ³⎜ + + ⎟ = − + ®⎜ + ⎟ + ³¤ − ⎝ ⎠ ⎝ ⎠; (3)2 2г г г г гг г г 2 2 p div r ,C T u T v T T T qt x y x y⎛ ⎞ ⎛ ⎞ ³ ⎜ + + ⎟ = ⎜ + ⎟ − ⎝ ⎠ ⎝ ⎠(4)где div rx ryrq qqx y = + .1213x1L 1 0yHhитTитlитT еgqh2h1l1h1l1Рис. 1. Область решения: 1 - твердые стенки,2 - газовая полость, 3 - источник тепла g
Ключевые слова
сопряженный теплоперенос,
естественная конвекция,
излучение,
математическое моделирование,
conjugate heat transfer,
natural convection,
radiation,
mathematical simulationАвторы
Мартюшев Семен Григорьевич | Томский государственный университет | аспирант кафедры теоретической механики | |
Шеремет Михаил Александрович | Томский государственный университет | кандидат физико-математических наук, доценткафедры теоретической механики | sheremet@math.tsu.ru |
Всего: 2
Ссылки
Соковишин Ю.А., Мартыненко О.Г. Введение в теорию свободно-конвективного теп- лообмена. Л.: Изд-во Ленингр. ун-та, 1982. 224 с.
Андрижиевский А.А., Лукашевич А.Г., Трифонов А.Г. Пространственное моделирование тепловых сборосов в водные системы // V Минский Международный форум по тепло- и массообмену (электронная версия докладов). Минск: Институт тепло- и массообмена им. А.В. Лыкова, 2004.
Лыков А.В. Теория теплопроводности. М.: Высшая школа, 1967. 600 с.
Жукаускас А.А. Конвективный перенос в теплообменниках. М.: Наука, 1982. 472 с.
Спэрроу Э.М., Сесс Р.Д. Теплообмен излучением. Л.: Энергия, 1971. 296 с
Вержбицкий Г.В. Основы численных методов. М.: Высшая школа, 2002. 840 с.
Пасконов В.М., Полежаев В.И., Чудов Л.А. Численное моделирование процессов тепло и массообмена. М.: Наука, 1984. 288 c.
Роуч П. Вычислительная гидродинамика. М.: Мир, 1980. 616 c.
Тарунин Е.Л. Вычислительный эксперимент в задачах свободной конвекции. Иркутск: Изд-во Иркут. ун-та, 1990. 225 с.
Самарский А.А. Теория разностных схем. М.: Наука, 1977. 656 с.
De Vahl Davis G. Natural convection of air in a square cavity: a bench mark numerical solution // Int. J. Numer. Methods Fluids. 1983. V. 3. Р. 249 - 264.
Hortmann M., Peric M., Sheuerer G. Finite volume multigrid prediction of laminar natural convection: benchmark solutions // Int. J. Numer. Methods Fluids. 1990. V. 11. Р. 189 - 207.
Kalita J.C., Dalal D.C., Dass A.K. Fully compact higherorder computation of steady-state natural convection in a square cavity // Phys. Rev. E64. 2001. Р. 1 - 13.
Bejan A. Convection Heat Transfer. John Wiley & Sons Inc., 1995.
Merrikh A.A., Lage J.L. Natural convection in an enclosure with disconnected and conducting solid blocks // Int. J. Heat and Mass Transfer. 2004. V. 46. P. 995 - 1008.
House J.M., Beckermann C., Smith T.F. Effect of a centered conducting body on natural convection heat transfer in an enclosure // Numer. Heat Transfer. Part A. 1990. V. 18. Р. 213 - 225.