УДК: 533.72
МОЛЕКУЛЯРНЫЙ ТЕПЛООБМЕН С ГАЗООБРАЗНОЙ СРЕДОЙ СИЛЬНО НАГРЕТОЙ НЕПОДВИЖНОЙ ТВЕРДОЙ УМЕРЕННО КРУПНОЙ СФЕРИЧЕСКОЙ ЧАСТИЦЫ
Е.Р. Щукин, Н.В. Малай, З.Л. Шулиманова
* Белгородский государственный университет, ул. Студенческая, 14, Белгород, 308007, Россия, e-mail: e-mail: malayObsu.edu.ru **Институт высоких температур РАН, Москва, 127412
Аннотация. Проведено математическое моделирование процесса молекулярного теплообмена с окружающей средой неподвижной умеренно крупной твердой сферической аэрозольной частицы при значительных перепадах температуры в её окрестности. Полученные формулы позволяют, с учетом скачка температуры и зависимости коэффициента теплопроводности от температуры, непосредственно находить распределение температуры в окрестности частицы и величину молекулярного потока тепла, отводимого от поверхности частицы. Анализ теоретических результатов показал, что увеличение температуры поверхности частицы приводит к монотонному возрастанию скачка температуры газа у её поверхности. Это в случае умеренно крупной частицы может привести к сильному уменьшению величины молекулярного потока тепла, отводимого от её поверхности.
Ключевые слова: газообразная среда, распределение температуры, твердые аэрозольные частицы, молекулярный теплообмен.
Введение. Температура поверхности твердых аэрозольных частиц, нагреваемых тепловыми источниками электромагнитной или химической природы [1-5], может значительно превышать температуру несущей газообразной среды [1-5]. Такой нагрев аэрозольные частицы могут испытывать, например, в зонах прохождения лазерного излучения через аэрозоли при
их просветлении, диагностике и сгорании [2-5].
Нагретые частицы могут оказывать значительное влияние на характер распределения температуры в аэрозоле. В этом случае при проведении математического моделирования, протекающих в аэрозоле, зависящих от температуры физических процессов, необходимо учитывать и молекулярный теплообмен аэрозольных частиц с газовой средой.
В большей части встречающихся на практике аэродисперсных систем среднее расстояние между частицами значительно превосходит их характерные размеры, а числа Рейнольдса и Пекле частиц много меньше единицы [4-7]. В таких системах математическое моделирование, зависящих от температуры физических процессов, можно проводить, основываясь на знании закономерностей теплообмена с бесконечной средой одиночных неподвижных аэрозольных частиц [5,8]. Поэтому изучение закономерностей и молекулярного теплообмена с бесконечной газообразной средой одиночных сильно нагретых неподвижных аэрозольных частиц представляет значительный научный и практический интерес.
Большое влияние на теплоперенос и распределение температуры в аэрозоле могут оказывать нагреваемые тепловыми источниками твердые аэрозольные частицы с числом Кнудсена Кп < 0.3 [1-5]. В случае сферических частиц Кп = А^/К [14-17], где К — радиус частицы, А^ — средняя длина свободного пробега молекул газа у поверхности частицы.
Математическое моделирование процесса молекулярного теплообмена с газовой средой твердых частиц с Кп < 0, 3 можно проводить, используя нелинейное гидродинамическое уравнение баланса тепла [10,13]. При этом на поверхности частицы учитывают газокинетическое граничное условие для скачка температуры [9,11-14]. Условие для скачка температуры позволяет учитывать влияние на процесс молекулярного теплообмена, окружающего каждую частицу, тонкого слоя Кнудсена [9]. Если при оценке теплообмена температуру газа у поверхности частицы можно считать равной температуре поверхности частицы, то аэрозольную частицу называют крупной (в случае твердых сферических частиц с коэффициентами аккомодации равными единице к крупным относят частицы с числом Кп < 0, 01). Если скачок температуры оказывает значительное влияние на процесс теплообмена, то частицу называют умеренно
крупной (при коэффициентах аккомодации равных единиц. Для таких частиц выполняются условия 0, 01 < Кп < 0, 3).
Формулы, приведенные в опубликованных до настоящего времени теоретических работах [15- 17], позволяют непосредственно оценивать молекулярный теплообмен с газообразной средой только сильно нагретых неподвижных крупных твердых, в частности, сферических аэрозольных частиц. Ниже в квазистационарном приближении проведено математическое моделирование процесса молекулярного теплообмена с газообразной средой неподвижной умеренно крупной твердой сферической частицы. Найденные при этом формулы позволяют при заданной температуре поверхности крупной или умеренно крупной частицы непосредственно оценивать распределение температуры в её окрестности и величину, отводимого от поверхности частицы молекулярного потока тепла. Численный анализ, полученных формул показал, что в случае сильно нагретых умеренно крупных твердых сферических частиц скачок температуры может оказать заметное влияние на процесс молекулярного теплообмена.
Моделирование теплообмена с учетом скачка температуры. Пусть в однокомпо-нентном газе с температурой Тете и давлением рнаходится неподвижная умеренно крупная твердая сферическая частица с радиусом Л. Внутри частицы действуют тепловые источники, которые могут вызвать её сильный нагрев. Коэффициент теплопроводности вещества частицы значительно больше коэффициента теплопроводности газообразной среды. При этом распределение температуры Тр вдоль поверхности частицы близко к однородному и в связи с этим температуру Тр можно считать постоянной величиной. В связи с малыми временами тепловой релаксации, процесс теплопереноса в окрестности частицы протекает квазистационарно. Радиус частицы достаточно мал для того, чтобы можно было пренебречь влиянием гравитационной конвекции на процесс переноса тепла в её окрестности. Температура поверхности частицы считается известной.
Коэффициент теплопроводности газа ке зависит от температуры газа Те. У недиссоции-рованных газов зависимость коэффициента ке от Те близка к степенной [18]. Поэтому, оценивая молекулярный теплообмен частицы с недиссоциированной газообразной средой, можно использовать следующую степенную зависимость ке от Те:
ке — кете (Те/Тете) , (1)
где кете - значение коэффициента теплопроводности при температуре газа Те — Тете. Например, при Тете — 293К, р — 1 бар формула (1) в интервале от 150 К до 2000К при ш — 0.85, кете — 0, 0255 Вт/м-К описывает с точностью до 5% зависимость от коэффициента ке воздуха и с точностью до 3% в интервале от 100К до 6000К зависимость от Те коэффициента ке гелия при ш — 0, 697, кете — 0,149 Вт/м-К [18]. Это достаточно хорошо показывают данные табл. 1.
Таблица 1
Экспериментальные [18] и найденные с помощью формулы (1) зависимости коэффициента ке воздуха и гелия от температуры Те при давлении р — 1 бар и Тето — 293 К.
Воздух Те, К 150 220 290 500 800 1500 1700 1900 2000
ке ■ Ю3 Вт/м-К [181 13.8 19.8 25.5 40.7 57.3 100 113 128 137
ке • 10й Вт/м-К; и = 0.85, ке ос = 0.0255 Вт/м-К 14.4 20.0 25.3 40.2 59.9 102.2 113.6 124.9 130.5
Гелий Те, К 100 300 600 1000 2000 3500 4500 5500 6000
ке ■ 10й Вт/м-К [181 72.0 151 250 354 579 826 970 1180 1200
ке • 10й Вт/м-К; и = 0.697, ке со = 0.149 Вт/м-К 70.4 151.5 245.6 350.6 568.3 839.5 1000 1150 1222
Тр — Tes = —Кт , dr
При рассмотренных выше условиях процесс молекулярного теплообмена твердой частицы с газообразной средой протекает сферически симметрично. Математическое моделирование сферически симметричного квазистационарного процесса теплообмена целесообразно проводить в сферической системе координат с началом в центре частицы. В этой системе координат распределение температуры Те в окрестности одиночных крупных и умеренно крупных частиц описывается уравнением переноса тепла (2) [6,10,16,17]
с1 ЛТе ^ Г 21
с1г е Лт '
в котором г радиальная координата; ке = ке^ (Те/Те^)ш — коэффициент теплопроводности газа; ¿е = (Те/Тете) — обезразмеренная температура газа. При заданной величине температуры поверхности частицы Тр решение (2) нужно проводить совместно с граничными условиями (3) [11-13]:
¿Те
) Тез — Те\г=л, Те\г—^оо — Теоо , (3)
т=Я
где Те8 — значение у поверхности частицы интерполированной из объема температуры газа [14-18]; разность ДТе5 = Тр — Те8 называют скачком температуры газа у поверхности частицы. Найденные в первом приближении по длине Л^ аналитические выражения для коэффициента Кт в случае однокомпонентного газа имеют следующий вид [11-14]:
Кт = СтЛs . (4)
В выражении (4) Ст — коэффициент скачка температуры, зависящий от коэффициентов аккомодации тангенциального импульса дт и энергии де [11-14]; Лs —- средняя длина свободного пробега молекул газа при температуре Тез у поверхности частицы. У большей части встречающихся на практике аэрозольных частиц коэффициент цт близок к единице [19], а значения коэффициента де лежат в пределах от 0,8 до 1 [20]. При коэффициентах аккомодации равных единице значения скачка температуры близки к Ст = 2, 2 [11-14]. Значения Лs в (4) можно оценивать по формуле
^ = Л<x/tes , (5)
где — средняя длина свободного пробега молекул газа при температуре Те(Х, t es — (Tes/Te'x).
Перейдя в граничной задаче (2) - (3) к безразмерной температуре с учетом (4), (5) получаем:
^^^ ^^в 0 /
(¿Г 6 (¿Г '
Aíeí¡ — —Ст/^oote ~Т~
dr
j tes — te\r=R , te\r—— 1 j (7)
r=R
где А^ = tp — tes - обезразмеренный скачок температуры ¿е, tp = (Тр/Тете). Аналитическое решение граничной задачи (6) - (7), позволяющее при (Лте/Я) tes < 0, 3 с относительной точностью до 0,02% находить значения имеет следующий вид:
te —
В выражении (8)
R
1 + - " 1)
r
1/(1+")
(8)
tes — tp ^tes j (9)
где скачок обезразмеренной температуры Ates равен:
Ate, = (^i - \J А\ — 4A0A-^j j 2A2
(10)
Ao = e(tp - , Ai =[1 + e (2tp - (1 - w) í-w)] , A2 = e
1 + ÍÜÜ^) r (i+O + 2 p
e =
Ct (A,
1 + w \ R
Выражение для отводимого от поверхности частицы молекулярного потока тепла (т равно:
Я{тм) = -4тгК2ке ^ т ат
(M)
т (11)
r=R
Подставляя в (11) функцию (8), получаем где коэффициент
Q(TM) = -inRk^Tf ,
f(M) = (4+ш - 1) / (1 + w) .
(12) (13)
R = 1Д мкм
R = 5 мкм
R = 15 мкм
7 tr
Рис. 1. Кривые зависимости скачка Ates обезразмеренной температуры воздуха te от обезразмеренной температуры поверхности частиц tp при R =1,1 мкм; 5 мкм; 15 мкм.
Полученные в процессе решения задачи о молекулярном теплообмене формулы (8) и (12) позволяют непосредственно, при известных температуре и радиусе поверхности умеренно крупной твердой сферической аэрозольной частицы, находить распределение безразмерной температуры te в её окрестности и величину молекулярного потока тепла, отводимого от поверхности частицы. При выводе этих формул было учтено влияние, оказываемое на процесс молекулярного теплообмена частицы с газообразной средой, скачка температуры и зависимости
зо
коэффициента теплопроводности газа от температуры. Найденные формулы позволяют проводить оценки как при малых, так и больших перепадах температуры в окрестности частицы.
Рис. 2. Кривые зависимостей интерполированной температуры Ь, от числа Кнудевна Кп — Ад/Д частиц при Ьр — 3; 5; 7.
Проведенный с помощью формул (8) и (12) численный анализ показал, что увеличение температуры поверхности крупных и умеренно крупных частиц приводит к монотонному возрастанию скачка температуры газа. Крые, показывающие зависимость обезразмеренного скачка температуры АЬез от обезразмеренной температуры поверхности частицы Ьр при Д — 1,1 мкм, Д — 5 мкм, Д — 15 мкм, приведены на рис. 1.
Максимальные скачки температуры возникают у поверхности умеренно крупных частиц, что может привести к сильному уменьшению величины молекулярного потока тепла, отводимого от их поверхности. Это достаточно хорошо показывает ход приведенных на рис. 2 и рис. 3 кривых зависимостей от числа Кнудсена интерполированной температуры Ьез и отношения
¡ТМ) — 0м )/4пДЙе^ Те^ (14)
при Тр — 3, Тр — 5, Тр — 7.
Температуру газа Ьез у поверхности крупных частиц можно считать равной температуре поверхности частиц с точностью до 1,5%. Допускаемая при этом при вычислении молекулярных потоков тепла 0ТМ) ошибка не превышает 2,5%. При построении кривых на рис. 1-3 расчеты были проведены для находящихся в воздухе с Тете — 293 К и давлением р — 1 бар частиц с коэффициентами аккомодации равными единице.
О 0.05 0,1 0,15 0,2 0,25 0,3 Xs jR
Рис. 3. Кривые зависимостей отношения /yM) = Qj^^/4nRkecxTt от числа Кнудевна при tp = 3; 5; 7.
Зависимости от радиуса частиц чисел Кнудсена Кп = Л^/Я в воздухе при Тр = 3, Тр = 5, Тр = 7 приведены в табл. 2.
Таблица 2
Зависимость от радиуса частицы Я числа Кнудсена в воздухе
с Теж = 293 К, = 1 бар при Тр = 3, Тр = 5, Тр = 7.
Tp = 3 R, мкм 0.485 1 3 5 7 9 15 30
A s/R 0.300 0.161 0.059 0.036 0.026 0.20 0.012 0.006
Tp = 5 R, мкм 0.780 1 3 5 7 9 15 30
A s/R 0.300 0.245 0.094 0.058 0.042 0.033 0.020 0.010
Tp = 7 R, мкм 1.08 3 5 7 9 15 30 43
A s/R 0.300 0.126 0.080 0.058 0.046 0.028 0.014 0.010
Литература
1. Кузиковский А.В. Динамика сферической частицы в мощном оптическом поле // Изв. ВУЗов. Физика. - 1970. - № 5. - С.89-94.
2. Зуев В.Е., Землянин А.А., Копытин Ю.Д., Кузиковский А.В. Мощное лазерное излучение в атмосферном аэрозоле / Новосибирск: Наука, 1984. - 224 с.
3. Bennet U.S., Rosasco G.I. Heating microscopic particles with laser beams //J. Appl. Phys. -1978. - 49;2. - P.640-647.
4. Букатый В.И., Копытин Ю.Д., Погодаев В.А. Горение углеродных частиц инициированное лазерным излучением // Изв. ВУЗов. Физика. - 1983. - № 2. - С.14-22.
5. Букатый В.И., Суторихин И.А., Краснопевцев В.Н., Шайдук А.М. Воздействие лазерного излучения на твердый аэрозоль / Барнаул: АГУ, 1994. - 196 с.
6. Фукс Н.А. Механика аэрозолей / М.: Изд-во АН СССР, 1955. - 352 с.
7. Спурный К., Йех К., Седлачек Б., Шторх О. Аэрозоли / М.: Атомиздат, 1964. - 360 с.
8. Фукс Н.А. Испарение и рост капель в газообразной среде / М.:Изд-во АН СССР, 1958. -92 с.
9. Коган М.Н. Динамика разреженного газа / М.: Наука, 1967. - 250 с.
10. Лойцянский Л.Г. Механика жидкости и газа / М.: Наука, 1970. - 904 с.
11. Loyalka S.K., Sielvert C.E., Thomas I.R. Temperature - jump problem with arbitrary accommodation // Phys. Fluids. - 1978. - 21;5. - P.854-855.
12. Яламов Ю.И., Поддоскин А.Б., Юшканов А.А. О граничных условиях при обтекании неоднородно нагретым газом сферической поверхности малой кривизны // ДАН СССР. -1980. - 254. - С.343-346.
13. Маясов Е.Г., Юшканов А.А., Яламов Ю.И. О термофорезе нелетучей сферической частицы в разреженном газе при малых числах Кнудсена // Письма в ЖТФ. - 1988. -14;6. - С.498-502.
14. Латышев А.В., Юшканов А.А. Аналитическое решение модельного БГК- уравнения Больцмана в задаче о температурном скачке с учетом аккомодации энергии // Математическое моделирование. - 1992. - 4;10. - С.61-66.
15. Щукин Е.Р., Бахтилов В.И. Нагрев и испарение крупных сфероидальных частиц под действием внутренних источников тепла // Физика аэродисперсных систем и физическая кинетика: Сборник/ МОПИ им.Н.К.Крупской. - 1979. - Вып. 4. Ч.1. - С.144-165 / Деп. в ВИНИТИ, №3828-79.
16. Щукин Е.Р. Квазистационарное испарение и рост капель чистых веществ при значительных перепадах температуры в их окрестности //M.: Объединенный институт высоких температур РАН, 1995. - 88 с. Деп. в ВИНИТИ, №412-В95.
17. Shchukin E.R. Solution of some non-linear problems in the theory of heating, vaporization and burning of solid particles and drops / Mathematical Modeling Problems, Methods, Applications. (Edited by Uvarova L.A. et al.) / New York: Kluwer Academic, Plenum Publishers, 2001. - P.255-267.
18. Варгафтик Н.Б. Справочник по теплофизическим свойствам газов и жидкостей / М.: Наука, 1972. - 720 с.
19. Баранцев Р.Г. Взаимодействие разреженных газов с обтекаемыми поверхностями / М.: Наука, 1975. - 344 с.
20. Коленчиц О.А. Тепловая аккомодация систем газ-твердое тело / Минск: Наука и техника, 1977. - 126 с.
MOLECULAR THERMAL EXCHANGE WITH GAS MEDIUM OF STRONG HEATED IMMOVABLE HARD INTERMEDIATELY GREAT
SPHERICAL PARTICLE
E.R. Shchukin, N.V. Malay, Z.L. Shulimanova
* Belgorod State University, Studencheskaya St., 14, Belgorod, 308007, Russia, e-mail: e-mail: malay@bsu.edu.ru **High Temperature Institute of RAS, Moscow, 127412, Russia
Abstract. It is proposed the mathematical modeling of molecular heat exchange of intermediately great hard spherical aerosol particle with surround medium at essentially large temperature difference in its neighborhood. Obtained formulas permits to find directly the temperature distribution in the particle neighborhood and the value of molecular heat flux from the particle surface. It may be done with the account of the temperature step and the thermal conductivity coefficient dependence on temperature. Analysis of theoretical results has shown that the increase of the particle surface temperature leads to monotone increase of gas temperature step near the surface. In the case of intermediately large particle, this may be to lead to strong decrease of molecular heat flux from its surface.
Key words: gaseous medium, temperature distribution, hard aerosol particles, molecular thermal exchange.