CC BY
32
УДК 533.72
К ВОПРОСУ О ГИДРОДИНАМИЧЕСКОМ СОПРОТИВЛЕНИИ СФЕРИЧЕСКОЙ ЧАСТИЦЫ С ОДНОРОДНЫМ ВНУТРЕННИМ ТЕПЛОВЫДЕЛЕНИЕМ
А.А. Стукалов
Белгородский государственный университет 308015, Белгород, ул. Победы, 85 e-mail: stykalov@bsu.edu.ru
В приближении Стокса проведено теоретическое описание стационарного движения нагретой аэрозольной частицы сферической формы в неизотермической газообразной среде, внутри которой действуют равномерно распределенные источники (стоки) тепла постоянной мощности. При рассмотрении движения предполагалось, что средняя температура поверхности частицы может существенно отличаться от температуры окружающей ее газообразной среды. В процессе решения гидродинамических уравнений получено аналитическое выражение для гидродинамической силы, действующей на равномерно нагретую аэрозольную частицу, с учетом зависимости вязкости и плотности от температуры.
Ключевые слова: гидродинамическое сопротивление, движение нагретых частиц в газе.
1. Постановка задачи. Движение нагретых частиц в вязких газообразных средах рассматривалось в ряде работ [1-5]. Под нагретой частицей понимают частицу, средняя температура поверхности которой по величине значительно превышает температуру окружающей среды. Нагрев поверхности частицы может быть обусловлен, например, протеканием объемной химической реакции, процессом радиоактивного распада вещества частицы и т.д.
Нагретая поверхность оказывает значительное влияние на теплофизические характеристики окружающей среды и тем самым может существенно повлиять на распределение полей скорости и давления в ее окрестности.
В настоящее время достаточно подробно рассмотрено движение твердых сферических частиц как при малых относительных перепадах в их окрестности, так и при значительных [1-5].
Рассмотрим движение твердой аэрозольной частицы в вязкой неизотермической газообразной среде, внутри которой действуют равномерно распределенные источники (стоки) тепла постоянной мощности при малых числах Рейнольдса под действием некоторой силы, например, электромагнитной, термофоретической, диффузиофоретиче-ской и т.д. Если перейти в систему координат, связанной с частицей, то задача по существу сводится к задаче обтекания нагретой твердой частицы сферической формы плоскопараллельным потоком газа скоростью U » (U » || OZ).
При рассмотрении обтекания нагретой сферы вводятся следующие допущения: 1) все процессы, происходящие в системе частица - газ, рассматриваются в квазистацио-нарном приближении, что возможно в силу малости времени тепловой релаксации частицы; 2) определяющими параметрами задачи являются коэффициенты сре ,
Рею ,№еоЛ» и сохраняющиеся в процессе движения частицы величины - R, U» . Из этих параметров можно составить две безразмерные комбинации: число Рейнольдса
Re^ = (Ре» U»R)/ ^е» <<1 и тепловое число Пекле - Ре» = (среиюR Рео)/ Лео <<1,
где R — радиус частицы [6-7], U» = U» - характерная скорость (величина скорости
268 НАУЧНЫЕ ВЕДОМОСТИ № 9 (49) 2008
набегающего потока); 3) движение частицы рассматривается при значительных (больших) относительных перепадах температуры. Под относительным перепадом температуры понимают разность температур между поверхностью частицы и областью вдали от нее. Относительный перепад температуры считается малым, если выполняется неравенство (T - TJ /Тет<<1, где TS — средняя температура поверхности частицы. При выполнении этого неравенства коэффициенты молекулярного переноса (вязкости и теплопроводности) и плотность можно считать постоянными величинами. Если (Ts — Tj /Тет ~O(1), то относительный перепад температуры считается значительным. В этом случае уже нельзя считать эти величины постоянными. В работе при описании свойств газообразной среды рассматривается степенной вид зависимости динамической вязкости и теплопроводности от температуры: jue = Meo t^ , Ae = Aeo t“ ,
Ai = , где Meo = Me (Tj , ^ = Ae (TJ , Aio = Al (TJ , tk = Tk/Teo , k = e- 1
[8]; 4) коэффициент теплопроводности частицы по величине много больше коэффициента теплопроводности газа. Это ограничение приводит к тому, что в коэффициенте вязкости можно пренебречь зависимостью по углу д в системе «частица - газ» и считать, что вязкость связана только с температурой te0 (r), т.е. Me (te (r,d)) -Me (te0 (r )) . Считая, при этом что te (r, д) = te0 (r Me (r,d) , где Ste (r,d)<< te0 (r). Явный вид функций te0 (r) и Ste (r ,d) определяются из решения конвективного уравнения теплопроводности с соответствующими граничными условиями; 5) частица образована однородным и изотропным по своим свойством веществом; 6) при описании нагретой частицы во внешних заданных полях используется гидродинамический подход [6-7].
В рамках сформулированных выше допущений уравнения для скорости и давления запишутся в виде [6-7]:
д P д ( Г dUj ЗЩ 2 _ дЦ;П Л
J. +-------х m
д xk д x j
Me
v
дxk дxj 3 j дxn
k j n
(m, k, j = 1,2,3) (1.1)
jk
div( Pe U e) = 0, (1.2)
где (Х1, Х2, Х3) декартовы координаты
Эта система гидродинамических уравнений решается со следующими граничными условиями в сферической системе координат (Г ,6,ф)
Г = R, и = 0, Щ = 0. (1.3)
г ие ^ и„ соБ^ег -и„ Бш0е0, Ре ^ Рех. (1.4)
Помимо линеаризованного по скорости уравнения Навье-Стокса (1.1) и непрерывности (1.2) с граничными условиями (1.3) - (1.4), необходимо также решить уравнение теплопроводности, описывающее распределение температуры вне и внутри частицы с соответствующими граничными условиями на ее поверхности. Это связано с тем, что частица нагрета. Разность температур между частицей и газообразной средой может стационарно поддерживаться, например, за счет выделения тепла при химических реакциях на поверхности частицы, за счет радиоактивного распада вещества частицы, внешним излучением и т.п. В частности, если на частицу падает поток монохроматического излучения (длина волны Л0) интенсивностью /о, то
поглощаемая ею энергия равна ПR2/0Кф, где Я - радиус частицы, Кф — фактор
поглощения [9], и распределяется по объему сферы равномерно с плотностью qi.
Это допущение справедливо в том случае, если коэффициент теплопроводности частицы по величине много больше коэффициента теплопроводности внешней среды и Ла >> R .
Следовательно, к уравнениям (1.1) - (1.2) необходимо добавить уравнение тепло-переноса вне и внутри частицы
div(лsv те) = о, (1.5)
(1.6)
div(A1VTi) = -qi
с граничными условиями (1.7)
r = R, T = T, 4 дТ^ = 4д T 1 -- (T4 T4
д r
д r
+ a>a(T4 - С)
(1.7)
отражающие равенство температур и непрерывности радиальных потоков тепла с учетом излучения на поверхности частицы; У0 — постоянная Стефана - Больцмана, У — интегральная степень черноты тела [10]; q = const — постоянная плотность тепловых источников, действующая внутри частицы [9].
Сила, действующая на частицу со стороны потока, определяется по формуле [6-7]
Fz = J(- Pe cos в + yrr cos в - yr0sin é)r2 sin в dв dр. (1.8)
(S)
Здесь угг , Уre, ue и ue - компоненты тензора полных напряжений, радиальная и касательная компоненты массовой скорости U в сферической системе координат
a
Me
2
dU
dy
—divU 3 e
a
гв
Me
due
+ ■
Ue
U
dy y д в y
y
r / R
обезразмеренная радиальная координата.
2. Поле скорости и распределение температуры. Определение силы сопротивления. Чтобы найти силу, действующую со стороны газа на твердую нагретую сферу, нужно знать распределения температуры, массовой скорости и давления в ее окрестности. Интегрируя уравнения (1.5) - (1.6) с соответствующими граничными условиями, получаем
te (y> в)= *e0 (y ) =
1 + Г0
y
1+а
(2.1)
ti (y, в ) = ti 0 (y )= B0 + — + - jVo dy - dy
V
yy
y
1+ю
(2.2)
/
где постоянные интегрирования определяются из соответствующих граничных усло-
вий на поверхности неравномерно нагретой частицы; Г0 = t]+a — 1;
/0
24,
y
+1
T
J qidx , —0
1 + ю
ex —1
4^R 4iJex V
J q^V
T / T
TS ' Tec
T
средняя
температура поверхности неравномерно нагретой частицы, определяемая из решения трансцендентного уравнения (2.3.3)
1
t
S
270
НАУЧНЫЕ ВЕДОМОСТИ
№ 9 (49) 2008
ti:a-1
R2
rtI
- G0G1
(tSí-1)-
1 + а 3 А Г " 1 А ' (2'3)
е<& е^ е^
Перейдем теперь к нахождению полей скорости и давления. В сферической системе координат линеаризованное по скорости уравнение Навье-Стокса, описывающее распределение скорости и давления вне частицы (1.1) имеет вид [7]:
дР„ да 2
ду ду y
+ —&ГТ +
1 да
1 дре даг в
3
у д в ду у
+ ~ат в +
уд 1 да,
гв + сМ.
а
у
аг -
вв
а
рр
у
уд
вв + ctge
у
(аев-а?Л
(2.4)
(2.5)
Здесь агг,агв,авв и <г( координат и они равны
рр
компоненты тензора напряжений в сферической системе
' 2 дШ- - 2 dlvuе ду 3 е
г2 дЩ 2тте 2,.ттл
----- + -ие --divU е
у д в у г 3 е
а
рр
2 2 2
- и + - ctg в Щ - - divU е у у 3
а
г
j
ди 1 ди
в + '
и
и
ду у д у
Решая систему уравнений (2.4) - (2.5) методом разделения переменных и, учитывая допущение 4), в конечном итоге получаем уравнений аналогичное [5]. В [5] проведен подробный анализ решения линеаризованного по скорости уравнения Навье-Стокса и выражения для полей скорости и давления.
Зная выражения для полей скорости и давления, подставляя в (1.8), после интегрирования получаем следующее выражение для силы сопротивления сферы, внутри которой действуют источники тепла постоянной мощности:
^ = 6пR ^ п,, ^ (2.6)
где N | у=1 = G1 0!2 — G2G1/ , Ы2 | у=1 = G1 Gз — G3G1/ , 0[ — первая производная от соответствующих функций. Явный вид которых приведен в [5].
Чтобы оценить какой вклад внутренние источника тепла (нагрев поверхности) оказывает на силу сопротивления, можно рассмотреть наиболее простой случай, когда частица поглощает излучение как черное тело [8]. В этом случае поглощение происходит в тонком слое толщиной ¿Я << Я, прилегающем к нагреваемой части поверхности частицы. При этом плотность тепловых источников внутри слоя толщиной ¿Я определяется с помощью формулы
q,(г, в) =
J п
-----— cos в, — < в <п, R -5R < г < R
SR 2
0, 0 < в <п 2
где /0 - интенсивность падающего излучения.
На рис. 1 приведены значения функции / от интенсивности падающего излучения /о. Численные оценки проводились для частиц меди радиусом 100 мкм, взвешенных в воздухе при нормальных условиях.
Рис.1. График зависимости функции f от интенсивности падающего излучения 10
Литература
1. Kassoy D.R., Adamson T.C., Messiter A.F. Compressible Low Reynolds Number around a Sphere //Phys. Fluids. 1966. Vol. 9. № 4. P. 671-681.
2. Щукин Е.Р., Малай Н.В. Фотофоретическое и термодиффузиофоретическое движение нагретых нелетучих аэрозольных частиц //ИФЖ. 1988. Т. 54, № 4. С. 630-635.
3. Щукин Е.Р., Малай Н.В, Яламов Ю.И. Движение нагреваемых внутренними источниками тепла капель в бинарных газовых смесях //ТВТ. 1988. Т. 25, № 5. С. 1020-1024.
4. Берковский Б.М., Краков М.С., Никифоров И.В., Полевиков В.К. Гидродинамическое сопротивление эллипсоидальной капли при малых числах Рейнольдса //МЖГ. 1987. № 3. С. 4-8.
5. Малай Н.В., Щукин Е.Р., Стукалов А.А., Рязанов К.С. Гравитационное движение равномерно нагретой твердой частицы в газообразной среде //СО РАН ПМТФ. 2008. № 1. С. 74 - 80.
6. Хаппель Дж., Бреннер Г. Гидродинамика при малых числах Рейнольдса. М.: Мир. 1976. 630 с.
7. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т.6. Гидромеханика. М.: Наука. 1986. 736 с.
8. Брейтшгайдер Ст. Свойства газов и жидкостей. Инженерные методы расчета. М.: Химия. 1966. 535 с.
9. Борен К., Хафмен Д. Поглощение и рассеяние света малыми частицами. М.: Мир. 1986. 660 с.
10. Шейндлин А.Е. Излучательные свойства твердых материалов. М.: Энергия. 1974. 471 с.
TO THE QUESTION ON HYDRODYNAMICAL RESISTANCE OF THE SPHERICAL PARTICLE WITH THE HOMOGENEOUS INTERNAL THERMAL EMISSION
А.А. Stukalov
Belgorod State University Pobedy St., 85, Belgorod, 308015, Russia e-mail: stykalov@bsu.edu.ru
The theoretical description of stationary movement of the heated up aerosol particle of the spherical form in not isothermal gaseous environment with intervals distributed heat sources of constant capacity inside is studied in the Stokes approximation.
While consideration of movement it was supposed, that the average temperature of a surface of a particle can differ from temperature of the gaseous environment surrounding it essentially. During the decision of the hydrodynamic equations analytical expression for the hydrodynamic force working on in regular intervals heated up aerosol particle is received, in view of dependence of viscosity and density from temperature.
Key words: hydrodynamic resistance, movement of the heated up particles in gas.
