CC BY
53
MSC 8035M
К ВОПРОСУ О ТЕПЛООБМЕНЕ СФЕРИЧЕСКОЙ ЧАСТИЦЫ
В ГАЗООБРАЗНОЙ СРЕДЕ
Н.В. Малай, С.И. Цыбульников
Белгородский государственный университет, ул. Студенческая, 14, Белгород, 308007, Россия, e-mail: e-mail: malay@bsu.edu.ru
Аннотация. Получено приближенное решение задачи о теплообмене движущейся твердой нагретой сферической частицы со средой при малых числах Пекле и Рейнольдса. Рассмотрен случай, когда внутри частицы действуют неравномерно распределенные источники тепла произвольной природы. Задача решается методом сращиваемых асимптотических разложений по числу Рейнольдса. При рассмотрении теплообмена предполагалось, что средняя температура поверхности частицы может существенно отличается от температуры окружающей ее газообразной среды. При решении уравнения конвективного теплопереноса учитывался степенной вид зависимости коэффициентов теплопроводности и плотности газообразной среды от температуры.
Ключевые слова: поле температуры, метод сращиваемых асимптотических разложений.
Введение. Исследование многочисленных энергетических процессов связано с решением задач переноса теплоты и вещества. Перенос этих субстанций в твердых телах, жидкостях и газах подчиняется условно принятым линейным зависимостям, например, перенос теплоты - закону Фурье: плотность теплового потока (удельный тепловой поток) пропорционален температурному градиенту; молекулярный перенос вещества -закону диффузии Фика: плотность потока вещества пропорциональна градиенту концентраций (или разности диффузионных химических потенциалов). На основании этих линеаризованных законов выводятся соответствующие дифференциальные уравнения.
В дальнейшем, мы будем рассматривать случай однородной среды без диффузии и химических реакций, протекающих с конечной скоростью. Если скорость движения сплошной среды мала по сравнению со скоростью звука, то возникающие в результате движения изменения давления настолько малы, что вызываемым ими изменением плотности (и других термодинамических величин) можно пренебречь. Общие уравнения закона сохранения энергии вне и внутри частицы упрощаются и принимают следующий вид [1]:
РесРе (ИеУ) Те = (АеУТе) , (\iVTi) = -qi (1)
с краевыми условиями в сферической системе координат т,в,<р
дТ дТ
г=Н,Т, = ъ, К-д^ = А(Т? - ТЦ , (2)
r ^ то, Те = Te^,
r ^ 0, Т = то,
(3)
(4)
где Я - радиус частицы, аа - постоянная Стефана-Больцмана; а\ - интегральная степень черноты; ре, ие и сре - плотность, массовая скорость и удельная теплоемкость жидкости; qi - плотность тепловых источников, неоднородно распределенных в объеме частицы, за счет которых происходит нагрев ее поверхности. В уравнение (1) входит массовая скорость газообразной среды ие. Для определения этой величины использовались результаты работы [2].
В литературе первое уравнение в выражении (2) называется конвективным уравнением переноса тепла. Оно описывает переноса тепла в газообразной среде за счет движения самой среды левая часть этого уравнения и за счет теплопроводности - правая часть. Наличие левой части делает это уравнение существенно нелинейным, что приводит к большим математическим трудностям при нахождении его решений.
Значимость процесса теплообмена, как в природе, так и в технике определяется, прежде всего, тем, что свойства тел самым существенным образом зависят от их теплового состояния, которое в свою очередь само определяется условиями теплообмена [3]. Эти условия оказывают существенное влияние на процессы изменения состояния вещества, механические, магнитные и другие свойства тел. Кроме того при проектировании экспериментальных установок, в которых необходимо обеспечить направленное движение нагретых частиц; при разработке методов тонкой очистки газов от аэрозольных частиц; при математическом моделировании процесса осаждения частиц в разнотемпературных плоскопараллельных каналах и т.п. необходимо знать поле температуры в их окрестности.
Таким образом, многие задачи по тепло-и массопереносу с которыми сталкиваются сегодня физики, инженеры и специалисты по прикладной математике, не поддаются точному решению. Среди причин, затрудняющих точное решение, можно указать, например, нелинейные уравнения движения. Переменные коэффициенты и нелинейные граничные условия на известных или неизвестных границах сложной формы. Для решения подобных задач мы вынуждены пользоваться различного рода приближениями, комбинируя аналитические и численные методы. Среди аналитических методов весьма мощными являются методы возмущений (асимптотических разложений) по большим или малым значениям параметра или координаты [4].
1. Постановка задачи. Метод решения. Рассматривается установившейся процесс теплообмена в потоке вязкой неизотермической газообразной среде, обтекающей твердую сферическую частицу радиусом Я. На большом расстоянии от сферы скорость потока равна ите. При описании свойств газообразной среды и частицы рассматривается степенной вид зависимости динамической вязкости и теплопроводности от температуры [5], таким образом
ре рете^е , Ае Аете^е , ре рете^е> А Аiтеti ,
рете ре (Тете) , рете ре (Тете) , Аете Ае (Тете) , ^те А (Тете) ,
tk = Тк/Тете, к = е, г, 0, 5 < а, в < 1, —1 < ш < 1-
Здесь и далее индекс е указывает на газообразную среду, индекс г — на принадлежность частице, а индекс то означает параметры газообразной среды на бесконечности, т.е. вдали от частицы.
Определяющими параметрами задачи являются коэффициенты рете, ^ете, Аете и сохраняющиеся в процессе движения частицы величины - Я, ите и Тете. Из этих параметров составить следующие безразмерные комбинации: число Рейнольдса И.ете =
(ретеЦтеЯ)/рете, число Пекле Рвте ^ете^Ргте? где Ргте (сререте/Аете число Прандтля,
ите — величина характерной скорости. Обезразмерим уравнения и граничные условия: Уе — ие/Цте) утхт/Я.
При е =И,ете ^ 1 решение уравнений гидродинамики находятся в виде ряда по е [2]. Перейдем теперь к решению уравнений теплопроводности. Если перейти к безразмерным величинам, то конвективное уравнение теплопроводности принимает вид
Конвективное уравнение переноса теплоты является нелинейным и для его решения применим метод сращиваемых асимптотических разложений [4]. В работе будем рассматриваться задача с возмущениями, действующими в очень узких областях, или зонах, в которых зависимые переменные испытывают достаточно резкие изменения. Ввиду наличия малого параметра при старшей производной эти узкие зоны часто оказываются лежащими вблизи границы области, в которой решается задача. Поэтому в задачах механики жидкостей и газов такие зоны получили название пограничными слоями, в механике твердого тела - областями краевого эффекта, в электрических приложениях - поверхностными, или скин-слоями.
Следовательно, во многих физических задачах, где резкие изменения зависимых переменных часто происходят внутри интересующих нас областей эти узкие зоны называются — ударными слоями (скачками уплотнения), точками перехода или стоксовскими линиями или поверхностями. Указанные быстрые изменения мы не можем исследовать с помощью обычных медленных масштабов; это приводит к необходимости вводить новые - быстрые, увеличенные или растянутые переменные [4]. Таким образом, исследование обыкновенных дифференциальных уравнений показали, что получить равномерно пригодные разложения в случаях, когда в некоторых областях изменения независимых переменных, где зависимые переменные испытывают резкие изменения обычными методами невозможно.
Один из методов, связанных с этой проблемой, заключается в построении прямых разложений (называемых внешними разложениями) с использованием исходных переменных и в построении разложений (называемых внутренними разложениями), описывающих эти резкие изменения и использующих увеличенные масштабы. Внешние разложения становятся непригодными в областях резких изменений, в то время как пригодность внутренних разложений нарушается при выходе из этих областей. Чтобы связать эти разложения, используют так называемую процедуру сращивания [4].
Внутренние и внешние асимптотические разложения обезразмеренной температуры представим как
ГС
«**(е>«) = Е .шс«,«),
п=0
(7)
где £ = еу - «сжатая» радиальная координата [4].
При этом требуется, чтобы /га+1//га ^ 0, /П+1//П ^ 0 при е ^ 0. Недостающие граничные условия для внутреннего и внешнего разложений вытекают из условия тождественности продолжения асимптотических разложений того и другого в некоторую промежуточную область, т.е.
(8)
Асимптотическое разложение решения внутри частицы, как показывают граничные условия на поверхности частицы, следует искать в аналогичном виде
(9)
п=0
Относительно функций /*(е) и /п(е) предполагается, что их порядок малости по е увеличивается с ростом п.
С учетом сжатой радиальной координаты имеем следующее уравнение для температуры **
Рг
1 1 о
+*
/ я** V * я** \ / \
\W-qI + у , г: 1 при е
(10)
(11)
Здесь ** = **(£, $), П - единичный вектор в направлении оси г.
2. Поля температур в окрестности нагретой аэрозольной частицы. При нахождении распределений температур в окрестности аэрозольной частицы мы ограничимся поправками первого порядка малости по е. Они определяются последовательно с учетом условия сращивания. Ввиду громоздкости полученных формул мы приведем нулевые и первые члены разложения в случае малых относительных перепадов температуры
Г
*е(£>0) = *еО + £*еИ Ъ(у,в) = *е0 +£*е1, и(у,6)=и о+еЬп, ^ео(у) = 1 + — ,
У
**0 = 1 *г0(у) = В0 +
1
1
4 п Я2ТетеАетеу2,/у
дгс1У+ — йу------------ фо &у ,
'1 У У2 Л
*е1(£>0) = уехр^РГоо^Ж- 1)|, *е1(у,0) = ^(1 -у) +
*г1 (у, ^) = СОБ % В1У +
СОй
Го /1 А1 А2
— + ш[ ----------------+ —
у2
2 4у3 2у
11
————-—qizdV + -
4 жИ, ТегсАегсу У У 3
У I 2
. Л у2
У
R2y2 2П +1 f+1
фп(у) = -т—Tf,— — Яг Рп(cos в) d(cos в) (п >0), х = cos 9, и = РгооГо .
AircTerc 2 J — 1
Постоянные интегрирования, входящие в выражения для полей температур определяются из граничных условий на поверхности частицы. Что же касается постоянных A1, A2 — они определяются из граничных условий для массовой скорости [2].
Заключение. Получены выражения для распределения температур вне и внутри аэрозольной частицы с учетом влияния движения среды (т.е. учтено влияние конвективного члена в уравнении теплопроводности) при произвольных относительных перепадах температуры в окрестности частицы. Поскольку частица нагрета, то при решении уравнений газовой динамики использовался степенной вид зависимости коэффициентов молекулярного переноса (вязкости и теплопроводности) и плотности газообразной среды от температуры. При нахождении полей температур вне и внутри аэрозольной частицы в случае значительных перепадов температуры предполагалось, что коэффци-ент теплопроводности частицы по величине много больше коэффициента теплопроводности газа, т.е. Aj ^ Ае. При выполнении этого условия в коэффициенте динамической вязкости ^e(r, в) можно пренебречь зависимостью по углу в в системе частица-газ и, считать, что ^e(te) ~ pe(teo) (предполагалась слабая угловая асимметрия распределения температуры). Это позволило рассматривать гидродинамическую часть отдельно от тепловой части, а связь между ними осуществляется через граничные условия.
Литература
1. Ландау Л.Д., Лифшиц Е.М. Механика сплошных сред / М.: ГИТ-ТЛ, 1954. - 795 с.
2. Малай Н.В., Щукин Е.Р., Стукалов А.А., Рязанов К.С. К вопросу о гравитационном движении равномерно нагретой частицы в газообразной среде // ПМТФ. - 2008. - 49;1. -С.74-80.
3. Брюханов О.Н., Шевченко С.Н. Тепломассообмен / М.: Ассоциация строительных вузов. 2005. - 460 с.
4. Ван-Дайк М. Методы возмущений в механике жидкости / М.: Мир, 1967. - 310 с.
5. Бретшнайдер C. Свойства газов и жидкостей. Инженерные методы расчета / М.: Химия. 1966. - 536 с.
ON HEAT EXCHANGE BETWEEN SPHERICAL PARTICLE IN GASEOUS MEDIUM
N.V. Malay, S.I. Tsybulnikov
Belgorod State University,
Studencheskaya St., 14, Belgorod, 308007, Russia, e-mail: e-mail: malay@bsu.edu.ru
Abstract. Approximated solution of the heat exchange between moving hard spherical particle at small Pecle’s and Reynolds’ numbers is obtained. The case when inhomogeneous distributed heat sources acted into particle is under consideration. The problem are solved by the method of joined asymptotic expansions on Reynolds’ number powers when averaged temperature on the particle surface may be essentially differed from the temperature of surround gaseous medium.
Key words: temperature distribution, method of joined asymptotic expansions.
