ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ
«БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ
ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ»
( Н И У « Б е л Г У » )
ИНСТИТУТ ИНЖЕНЕРНЫХ ТЕХНОЛОГИЙ И ЕСТЕСТВЕННЫХ НАУК
КАФЕДРА ТЕОРЕТИЧЕСКОЙ И МАТЕМАТИЧЕСКОЙ ФИЗИКИ
ОСОБЕННОСТИ ТЕПЛОПЕРЕНОСА ПРИ ЗНАЧИТЕЛЬНЫХ ПЕРЕПАДАХ ТЕМПЕРАТУРЫ В ОКРЕСТНОСТИ АЭРОЗОЛЬНОЙ ЧАСТИЦЫ СФЕРИЧЕСКОЙ ФОРМЫ
Магистерская диссертация
обучающегося по направлению подготовки
03.04.02 Физика
очной формы обучения
группы 07001637
Трубаева Алексея Андреевича
Научный руководитель:
д.ф.-м.н., профессор
Малай Н.В.
Рецензент:
к.ф.-м.н., доцент кафедры
естественнонаучных
дисциплин, БУКЭП
Стукалов А.А
БЕЛГОРОД – 2018
ОГЛАВЛЕНИЕ
Введение
3
Глава I. Постановка задачи
9
Глава II. Расчет поля температуры внутри крупной аэрозольной частицы
15
Глава III. Метод сращиваемых асимптотических разложений и поле температуры вне частицы
20
Глава IV. Поток тепла на поверхность частицы. Анализ полученных результатов
31
Заключение
35
Список использованной литературы
36
ВВЕДЕНИЕ
Тема исследования: влияние на теплообмен нагрева поверхности частицы и движения среды при значительных относительных перепадах температуры между поверхностью частицы и окружающей ее газообразной средой
вдали от нее.
Актуальность исследования: явления переноса представляют собой
неравновесные процессы, в результате которых в физической системе происходит пространственный процесс переноса вещества, импульса, энергии, энтропии или какой-либо другой физической величины. Причина процессов
переноса – возмущения, нарушающие состояние термодинамического равновесия: наличие пространственных неоднородностей состава, температуры
или средней скорости движения частиц системы. Перенос физической величины происходит в направлении, обратном ее градиенту, в результате чего
изолированная от внешних воздействий система приближается к состоянию
термодинамического равновесия. Явления переноса протекают стационарно,
если внешние воздействия поддерживаются постоянными,
При взаимодействии конденсированных сред, у которых различные
температуры, происходит обмен энергией. Интенсивность движения частиц
среды существенно зависит от ее температуры. Например, если температура
тела, имевшего меньшую температуру, увеличивается, то интенсивность
движения частиц тела с более высокой температурой уменьшается. Вследствие такого взаимодействия одно из тел остывает, а другое тело нагревается.
Перенос теплоты есть перенос энергии. Поток энергии, который передается
частицами тела с более высокой температурой частицам тела с более низкой
температурой в литературе принято называть тепловым потоком. Отсюда
следует, чтобы возник тепловой поток, т.е. возник процесс теплообмена
между различными областями пространства, заполненного вещественной
средой, необходимо и достаточно, чтобы в этих областях, имели место раз3
ные температуры. Перенос энергии осуществляется одним или несколькими
носителями. Носителями теплоты могут быть любые вещественные системы:
атомы, молекулярные комплексы, вихревые образования, гидродинамические
потоки и т.д.
Значимость процесса теплообмена в производстве, природе и т.д. определяется тем, что свойства тел самым существенным образом зависят от их
теплового состояния, которое в свою очередь само определяется условиями
теплообмена. Эти условия оказывают существенное влияние на процессы
изменения состояния вещества, механические, тепловые, магнитные и другие
свойства тел. Именно этим и объясняется интенсивное развитие теории теплообмена и то исключительно важное значение, которое ей уделяется в энергетике, химической технологии, авиастроении, медицине, сельском хозяйстве
и природе.
В различных областях химических технологий, гидрометеорологии,
сельском хозяйстве, медицине, охраны окружающей среды, т.е. в современной науке и технике нашли очень широкое применение многофазные смеси.
Наибольший интерес из многофазных смесей представляют так называемые
аэродисперсные системы. Аэродисперсные системы состоят из двух фаз, одна из которых есть частицы, а вторая – вязкая газообразная среда. Аэрозолями называют частицы, взвешенные в газообразной среде, а сами частицы –
аэрозольными.
Следует также отметить, что среднее расстояние между аэрозольными
частицами, встречающихся в аэродисперсных системах, намного больше характерного размера аэрозольной частицы. Это означает, что в таких системах
учет влияния аэрозоля на развитие теплообмена в аэродисперсной системе
можно проводить, основываясь на знание законов теплообмена и динамики
движения с бесконечной окружающей средой отдельных аэрозольных частиц. Без знания закономерностей такого теплообмена невозможно математическое моделирование эволюции аэродисперсных систем и решение такого
важного и интересного вопроса как целенаправленное воздействие на аэро4
золи. Поэтому изучение закономерностей теплообмена отдельных частиц с
бесконечной окружающей вязкой неизотермической газообразной средой является важной актуальной задачей, представляющей значительный теоретический и практический интерес.
В научной литературе аэрозольные частицы по размерам условно делятся на крупные, умеренно крупные и мелкие. Для классификации частиц
по размерам применяют критерий Кнудсена [1]:
Kn
cp
R
,
где ср средняя длина свободного пробега молекул газообразной среды;
R характерный размер частицы.
Если Kn 0.01, то такие частицы называются крупными, при выполнении неравенства 0.01 Kn 0.3 - умеренно крупными и при Kn 1
- мелкими.
При описании процесса теплообмена будем использовать термин «относительный перепад температуры». Под относительным перепадом температуры понимают отношение разности между средней температурой поверхности частицы TiS и температурой окружающей частицу вязкой газообразной среды вдали от нее T к последней. Относительный перепад температуры считается значительным, если имеет место следующая оценка
TiS T / T ~ 1 . В этом случае необходимо учитывать зависимость коэффициентов молекулярного переноса (вязкости, теплопроводности) и плотности газообразной среды от температуры, что существенно осложняет анализ
системы газодинамических уравнений и сама вязкая газообразная среда
называется неизотермической. Здесь и далее индексы «e» и «i» будем относить к газу и частице, индексом «s» – обозначены значения физических величин, взятых при средней температуре поверхности частицы равной TiS , а
индексом « » – обозначены средние значения физических величин, характеризующие газовую среду в невозмущенном потоке.
5
Нагрев поверхности частицы до температуры TiS осуществляется за
счет внутренних источников тепла неоднородно распределенных в ее объеме,
плотность этих тепловых источников будем обозначать через qi . В работе
считается, что плотность тепловых источников задана.
Внутренние источники тепла − это модельное представление, которое
удобное для описания реальных процессов, сопровождающихся выделением
тепла в объеме аэрозольной частицы. Так можно моделировать, например,
нагрев поверхности аэрозольной частицы под действием химической реакции, радиоактивного распада вещества частицы, поглощением электромагнитного излучения и т.д.
В работе мы не конкретизируем физическую природу плотности
тепловых источников qi ( r ). Это позволяет рассмотреть процесс теплообмена в общей постановке.
В работе рассматривается стационарный теплообмен между крупной
неравномерно нагретой аэрозольной частицы сферической формы и газообразной средой при значительных относительных перепадах температуры в ее
окрестности.
Объектом исследования является изучение явления теплообмена при
значительных относительных перепадах температуры.
Предметом исследования: влияние на теплоперенос нагрева поверхности аэрозольной частицы и движения газообразной среды с учетом степенного вида зависимости коэффициентов молекулярного переноса (вязкости, теплопроводности) и плотности газообразной среды от температуры.
Цель исследования – получить выражения для полей температур вне и
внутри крупной неравномерно нагретой частицы с учетом движения среды и
степенного вида коэффициентов молекулярного переноса (вязкости, теплопроводности) и плотности газообразной среды от температуры., выражения
для локального теплового потока и числа Нуссельта.
6
Исходя из поставленной цели, были сформулированы следующие задачи
исследования:
- ознакомление с историей развития теплообмена крупных аэрозольных
сферической формы в газообразных средах;
- изучить математические методы решения дифференциальных уравнений
второго порядка в частных производных в сферической системе координат;
- получить выражения для полей температур вне и внутри крупной аэрозольной частицы сферической формы, внутри которой действуют источники тепла с учетом степенного вида зависимости коэффициентов переноса (вязкости, теплопроводности) и плотности газообразной среды от
температуры при значительных относительных перепадах температуры в
ее окрестности;
- провести качественный анализ влияния нагрева поверхности частицы и
движения среды на распределения полей температур вне и внутри неравномерно нагретой аэрозольной частицы сферической формы, а также на
локальный тепловой поток и тепловое число Нуссельта.
Научная новизна исследования и практическая значимость исследования:
проведено теоретическое описание влияния на теплоперенос
нагрева поверхности крупной аэрозольной частицы и движения среды, получены выражения для полей температур вне и внутри частицы с учетом степенного вида зависимости вязкости, теплопроводности и плотности газообразной среды от температуры, а также локальный тепловой поток на частицу
и число Нуссельта.
Результаты выпускной квалификационной работы могут быть использованы при проектировании экспериментальных установок, в которых необходимо обеспечить направленное движение аэрозольных частиц; при разработке методов тонкой очистки газов от аэрозольных частиц; при оценке каналов просветления и т.д.
7
Апробация исследования. Основные результаты выпускной квалификационной
работы
практической
докладывались
конференции
на
III
«Современные
Международной
проблемы
научнофизико-
математических наук» ( 23 – 26 ноября 2017 г., Орел) и опубликованы в
сборнике «Современные проблемы физико-математических наук» (Орел
– 2017).
Структура работы.
Дипломная работа состоит из введения, четырех глав, заключения, списка литературы.
8
ГЛАВА I. ПОСТАНОВКА ЗАДАЧИ
Рассмотрим неравномерно нагретую крупную твердую сферическую
частицу (т.е. частицу, на поверхности которой не происходит фазовый переход) радиуса R , взвешенную в газообразной среде с температурой Te , плотностью e и вязкостью e . Поверхность частицы нагревается за счет источников тепла плотностью qi ( r ), которые неравномерно распределены в ее
объёме. (см. Рис. 1.). Процесс теплообмена будем описывать в сферической
системе координат r, , , начало которой связанно с центром масс аэрозольной частицы. Крупная аэрозольная частица обтекается бесконечным
плоскопараллельным потоком, со скоростью U ( U || OZ ). На большом расстоянии от частицы температура потока равна T , плотность e и вязкость
e . Задача осесимметрична. В окрестности частицы скорость потока и тем-
пература изменяются: скорость из-за возмущающего действия сферы, а температура – из-за неоднородного нагрева поверхности сферы. Поле скорости
считается известным из решения соответствующей гидродинамической задачи.
Рис. 1.
9
При рассмотрении теплообмена сделаны следующие физические допущения:
- характерные значения времен установления распределения полей температур и скорости течения в среде малы по сравнению с характерным
временем нагрева ее до максимальной температуры. Это означает, что
при теоретическом описании теплообмена в силу малости времен тепловой релаксации процессы теплообмена в системе частица-газ протекают квазистационарно;
- предполагается, что примеси в частице отсутствуют, т.е. она образована однородным и изотропным по своим свойством веществом;
- описание теплообмена рассматривается при значительных (больших)
относительных перепадах температуры. При описании свойств газообразной среды и частицы используется степенной вид зависимости вязкости, теплопроводности и плотности от температуры [2]: e e te ,
e e te , e e / te ,
i i T ,
i iti , где e e T , e e T ,
e e T ,
te Te / T ,
ti Ti / T ,
0,5 , , 1 ,
1 1;
- коэффициент теплопроводности частицы по величине много больше
коэффициента теплопроводности газообразной среды, что имеет место
для большинства газообразных сред. Это допущение приводит к тому,
что в коэффициенте вязкости можно пренебречь зависимостью по углу
в системе «частицы–газообразная среда» (предполагается слабая уг-
ловая асимметрия распределения температуры) и, следовательно, вязкость связана только с температурой te0 r , т.е. e te r, e te0 , При
этом te r, te0 r te r, , te r, te0 r , te r,, te0 r определяются из решения тепловой задачи. При таком допущении можно рассматривать гидродинамическую часть отдельно от тепловой части, а
связь между ними осуществляется с помощью граничных условий.
10
Поскольку идет нагрев частицы, т.е. поглощаемая энергия переходит в
тепловую, то через функцию qi ( r ) обозначим плотность неоднородно распределенных тепловых источников внутри объёма частицы. Так как на
решение задачи влияет только температура частицы, а не способ нагрева,
то считаем эту функцию заданной.
Таким образом, с учетом всех выше описанных допущений, распределение полей температуры вне Te и внутри частицы Ti , описывается
следующей системой уравнений [3,4]:
e c p U eTe div(eTe ) ,
div i Ti qi .
(1.1)
(1.2)
Здесь с p удельная теплоемкость при постоянном давлении; i , e коэффициенты теплопроводности газа и частицы; U e − массовая скорость газа
Система дифференциальных уравнений (1.1) ‒(1.2) решается со следующими граничными условиями в сферической системе координат
y r / R, , [1,3]. На поверхности частицы
y 1 , учитываются равенство
температур (1.3) и равенство потоков тепла с учетом излучения (1.4):
Te Ti ,
e
(1.3)
Te
T
i i o 1 RTi 4 T4 .
y
y
(1.4)
Здесь 0 постоянная Стефана-Больцмана, интегральная степень черноты вещества частицы.
На большом расстоянии от частицы ( y ) температура и давление
газообразной среды равна температуре и давлению невозмущенного потока:
11
Te T , Pe P .
(1.5)
В центре частицы ( y 0 ) нужно учесть конечность температуры:
Ti .
(1.6)
Граничные условия для компонент массовой скорости мы не приводим.
Как было отмечено выше, что поле скорости считается известным из решения соответствующей гидродинамической задачи. Это позволяет рассматривать процесс теплообмена для общего случая, т.е. для установившегося движения крупной твердой неравномерно нагретой частицы сферической формы. Установившееся движение означает, что на частицу кроме силы вязкого
сопротивления среды (силы Стокса) действует другая сила, направленная
противоположно силе сопротивления среды. Это может быть силы тяжести
(гравитационное осаждение), термофоретическая сила (термофорез), фотофоретическая сила (фотофорез) и т.д.
С учетом выше сказанного, в дальнейшем нам потребуется выражения
для компонент массовой скорости U e и U re , которые взяты из работы [5,6],
удовлетворяющие условию ограниченности решения при y , имеют вид:
Здесь f
U re U cos G( y), G( y) A1G1 A2G2 G3 ,
(1.7)
Ue U sin g ( y), g ( y) A1G4 A2G5 G6 .
(1.8)
1
Gk 3 yGk 3
, Gk 1
2
y 1
21
k 4,5,6,
f I , f II , G1 ,
G2 , G3 первые, вторые производные по y от соответствующих функций.
G1 y
1
n
C
G
y
C3,n n 3 ln y G1 y ,
,
3
3 1,n
y n 0
n 0
12
1
G2 y C2,n n 2 ln y G1 y .
y n 0
Значения коэффициентов С1,n n 1 , С2,n n 3 и C3,n n 4 определяются с помощью следующих рекуррентных соотношений:
С1,n
1
n 13n 2 13n 8 1 n 2n 3 2 n 2 C1,n 1
nn 3n 5
n 1n 23n 5 2 1 n 2 4 2 n 2 3 n 3 C1,n 2
n 2n 1n 3 1 n 3 3 C 1,n 3
1
n 13n2 n 6 1nn 1 n 2 C2,n 1 3 n 1
n 1n 3n 2
n 1n 23n 1 21nn 2 2 n 2 C2, n 2 n 2n 1n 3 3
C2, n
1 n 3 C2, n 3
4 1 4 ...n 1 4
2 n 2
n
k
1
6
k
02 k 0
n!
1
n 1 3n2 5n 4 1n 2 C3,n 1 n 1n 23n 4
nn 2n 3
21 n 1n 2 2 n 2 n 3 C3, n 2 n 2n 1n 3 1 n 3 3 C3, n 3
C3, n
3 n 3
3 n k 2n k 1 k
20 k 0
k 3k 2 16k 15C1,k k 16k 13 1 2k 5 2 C1,k 1
3k 1k 2 21k 2 3 C1, k 2 , 1
2 2
1
1
, 2 2
, 3
.
1
1
1 2
При вычислении коэффициентов С1,n , C2,n и C3,n по рекуррентным формулам необходимо учитывать, что C1,0 1 , C2,0 1 , C3,1 0 , С2, 2 1 ,
C2,1
1
2 1 2 6 4 , 33 3 10 3 1 2 , 4 /1 ,
20
60
8
13
2
2
0
1 1
2 1 2 6 4 4 3 1 2 3 3 3 4 4 1 , C3, 2 1 3 ,
4
15 4
C3,3 1 , C3,0 1, C1,n , C2,n и C3,n при n 0 равны нулю.
Постоянные интегрирования A1 и A2 определяются из граничных
условий на поверхности аэрозольной частицы.
Перейдем к безразмерным переменным в уравнениях и граничных
условиях, введя безразмерные координату, скорость и температуру следующим образом: yk xk R , t T T , Ve = U e /U , Vi = U i / U .
Материальные постоянные e , e , e и не меняющиеся в процессе
массообмена параметры Te , R и U являются определяющими. Из этих параметров можно составить безразмерную комбинацию, которую в научной
литературе называют числом Рейнольдса: Re e R U / e .
При 1 решение уравнений гидродинамики будем искать в виде
Ve =Ve0 Ve1 +…, Pe Pe(0) ε Pe(1) ... ( Re ).
(1.9)
Для нахождения локального теплового потока и теплового числа Нуссельта необходимо знать поля температур вне и внутри неравномерно нагретой частицы и при их нахождении мы ограничимся до первого порядка малости.
14
ГЛАВА II. РАСЧЕТ ПОЛЯ ТЕМПЕРАТУРЫ ВНУТРИ КРУПНОЙ
АЭРОЗОЛЬНОЙ ЧАСТИЦЫ
Найдем распределение температуры внутри частицы. Для этого необходимо решить следующее уравнение
divi Ti qi , i iT ,
(2.1)
с граничным условием
Ti .
(2.2)
В безразмерном виде оно принимает вид:
qi R 2
.
div ti ti
T i
Здесь ti
(2.3)
Ti
, i iti .
T
Будем искать решение уравнения (2.3) методом разделения переменных по полиномам Лежандра [7]. Для этого сначала разложим правую часть
уравнения (2.3) по полиномам Лежандра:
qi R 2
qin y Pn x , x cos .
λiT n 0
(2.4)
Неизвестная функция qin ( y) определяется следующим образом:
умножим выражение (2.4) на полином Pm x dx и проинтегрируем от -1 до
+1
qi R 2
Pm x
dx qin y P m x Pn x dx ,
λiT
n 0
1
1
1
1
и используя, свойство ортогональности полиномов Лежандра [7]
0; n 0
1 P m x Pn x dx 2 ; n m ,
2n 1
1
15
имеем следующее выражение для функция qin y
2n 1 R 2
qin y
qi Pn x dx .
2 λiT 1
1
(2.5)
В дальнейшем нам потребуются явный вид функций qi0 y и qi1 y
1
1
1 R2
3 R2
,
qi 0 y
q
r
,
dx
q
y
qi r, dx .
i
i1
2 λiT 1
2 λiT 1
(2.6)
Функцию ti у, будем искать в виде:
ti у, ti 0 у ti1 у, .
(2.7)
Линеаризуем выражение стоящее в левой части уравнения (2.3). До
первого порядка малости по имеем:
ti ti 0 ti ` ti 0
ti1
t
1
ti 0 1 i1 ti 0 ti 0 1ti1 .
ti 0
ti 0
Подставляя это выражение в уравнение (2.3), получаем систему из двух
уравнений для нахождения функций ti0 у и ti1 у,
div ti0ti 0 qi 0 ,
div ti 0 1ti1ti 0 ti1ti 0 qi1 cos .
(2.8)
(2.9)
Найдем решение уравнения (2.8). Введем новую функцию
f 0 ( y) ti10 ( y) ,
(2.10)
и тогда наше уравнение (2.8) принимает вид:
f 0 1 qi 0 .
(2.11)
В сферической системе координат, учитывая, что f 0 f 0 ( y ) имеем:
16
d 2 f 0 2 df 0
1 qi0 y .
d y2 y2 d y
(2.12)
Таким образом, мы получили неоднородное обыкновенное дифференциальное уравнение второго порядка. Из теории решения обыкновенных
дифференциальных уравнений известно, что общее решение обыкновенного
неоднородного уравнения второго порядка [8]
y '' p(x) y ' q(x) y f (x)
имеет вид [8]
( x) f ( x)
( x) f ( x)
y ( x ) C11 ( x ) C22 ( x ) 2 ( x ) 1
dx 1 ( x ) 2
dx ,
W
(
x
)
W
(
x
)
x0
x0
x
x
(2.13)
где функции p(x), q(x), f (x) определены и непрерывны на некотором интервале a , b , 1 , 2 функции, являющие решением соответствующего однородного уравнения, а W определитель Вронского. Точка x0 a, b и выбирается таким образом, чтобы соответствующие интегралы легко вычислялись.
В нашем случае функции 1 1; 2
W
1
; определитель Вронского
y
1
. Следовательно, для функции f 0 y имеем:
y2
y
1y 2
f 0 y B0
y (1 )qi 0 dy y 1 qi 0 dy .
y1
y
1
C0
С учетом выражения (2.10) получаем следующее выражение для функции ti0 у :
17
C
1
t0 y B0 0 0 dy 0 dy
y
y1
y
1
y
y
1
1
,
(2.14 )
R 2 y 2 1
где 0
qi dx .
2 iT 1
1
Константа C0 находится из граничного условия конечности температуры (2.2 ) и она имеет вид:
0
C0 y 0 dy ,
1
которую можно представить в следующем виде:
0
0 1
1 R2
1 qi y 2 dydx
C0 y 1 qi 0 y dy
2 T i
1
1 1
2
1 R 2 1
1 R 2 1
r 2 dr
2
qi y dx ( sin )d
q0
( sin )d
2 T i 1
2 T i R R 2 R
0 0
0 0
1 R 2 1
1 R 2 1
2
qi r sin drd
qi r 2 sin ddrd
3
3
2 R T i 0 0
4 R T i 0 0 0
R
2 R
1 R 2 1 1
qi dV ,
3 T i V V
dV r 2 sin drdd , x cos .
Таким образом,
С0
1
qi dV .
4R iT V
(2.15)
Перейдем к нахождению первого приближения для поля температуры
внутри частицы. Исходя из вида уравнения (2.9 ), будем искать решение для
функции ti1 у, в виде:
ti1 ( y, )
1
f 1 y , .
t ( y)
i0
Подставляя (2.16) в уравнение (2.9), получаем
18
(2.16)
f1 qi1 y cos .
(2.17)
Решение уравнения (2.17) ищем в виде
f1 y, i y cos ,
(2.18)
и, подставляя (2.18) в уравнение (2.17) в сферической системе координат получаем следующее уравнение
d 2 i 2 d i
i
2
qi1 .
d y2 y d y
y2
(2.19)
Применим к уравнению (2.19) формулу (2.13). Заметим, что в нашем
случае 1 y , 2 1 / y 2 , W 3 / y 2 .
Таким образом, имеем следующее выражение для функции ti1 у, :
y
y
N 0 cos
C1
1
1
ti1 y B1 y 2 y 2 dy 2 y1dy .
ti 0
ti 0
y
y
y 1
1
(2.20)
Константа C1 также находится, как и C0 (из граничного условия конечности температуры) и равна:
С1
где J 1
R
J1 ,
3 iT
(2.21)
1
qi z dV , z r cos , qi z dV дипольный момент плотности тепV V
V
1
R2 y2
ловых источников, Ф1
qi x d x .
2 iT 1
19
ГЛАВА III. МЕТОД СРАЩИВАЕМЫХ АСИМПТОТИЧЕСКИХ РАЗЛОЖЕНИЙ И ПОЛЕ ТЕМПЕРАТУРЫ ВНЕ ЧАСТИЦЫ
Исследование обыкновенных дифференциальных уравнений показали,
что получить равномерно пригодные разложения в случаях, когда в некоторых областях изменения независимых переменных, где зависимые переменные испытывают резкие изменения обычными методами (например, метод
возмущения, метод растягивания координат и т.п.) невозможно [9].
Один из методов, связанных с этой проблемой, заключается в построении прямых разложений (называемых внешними разложениями) с использованием исходных переменных и в построении разложений (называемых
внутренними разложениями), описывающих эти резкие изменения и использующих увеличенные масштабы. Внешние разложения становятся непригодными в областях резких изменений, в то время как пригодность внутренних
разложений нарушается при выходе из этих областей. Чтобы связать эти разложения, используют так называемую процедуру сращивания [9].
Сращивание представляет собой основную черту метода сращиваемых
асимптотических разложений [9]. Возможность сращивания основана на существовании области перекрытия, в которой пригодны как внутреннее, так
и внешнее разложения. Используя это перекрытие, можно получить точное
соотношение между конечными частными суммами. Реализация этой замечательной возможности осуществима только для возмущения параметра, которое неоднородно в координатах, или для возмущения координаты, которое
неоднородно по другим координатам. Нельзя срастить два различных параметрических разложения, таких, как разложение для больших и малых значений числа Рейнольдса, числа Маха и т.д. Такие ряды могут перекрываться в
том смысле, что они имеют общую область сходимости, но процесс аналитического продолжения дает только приближенное соотношение для некоторого конечного числа членов [9].
20
Воспользуемся методом сращиваемых асимптотических разложений
для нахождения поля температуры в окрестности неравномерно нагретой частицы.
Если перейти к безразмерным величинам, то конвективное уравнение
теплопроводности (1.1) принимает вид:
ε
Pr
Ve te divtete , te Te /T ,
te
(3.1)
с граничным условием
te y , 1 .
(3.2)
Решать это уравнение будем методом сращиваемых асимптотических
разложений. Внутренние и внешние асимптотические разложения обезразмеренной температуры t e будем искать в виде:
– внутренние асимптотические разложения
te ( y, ) f n ( )t en ( y, ) , ( f 0 ( ) 1),
(3.3)
n 0
– внешние асимптотические разложения
*
t e* ( , ) f n* ( )t en
( , ) ,
(3.4)
n 0
где y – "сжатая" радиальная координата [9].
При этом требуется, чтобы
f n 1
f n*1
0 , * 0 , при 0 .
fn
fn
(3.5)
Недостающие граничные условия для внутреннего и внешнего разложений вытекают из условия тождественности продолжения асимптотических
разложений того и другого в некоторую промежуточную область, т.е.
te ( y , ) te* ( 0, ) .
21
(3.6)
С учетом сжатой радиальной координаты имеем следующее уравнение
для температуры te* :
Pr * * *
Ve te divte*te* ,
*
te
te* 1
при .
Pe* , P при .
и Ve* , = nz + Ve1* , + …,
(3.7)
(3.8)
Здесь Vre* Vre* , , Ve* Ve* , ; te* te* , ; Pr e c pe e число
Прандтля, nz – единичный вектор в направлении оси Оz.
Построение решений начинается с определения нулевого члена внешнего разложения. В данном случае, очевидно, задаче удовлетворяет решение:
te* te*0 1.
(3.9)
Найдем нулевой член внутреннего разложения. При 0 имеем следующее уравнение:
div te te 0 .
(3.10)
Из (3.3) имеем:
te ( y, ) te0 ( y, ) .
Недостающие постоянные интегрирования определяются из условия
сращивания (3.6). В случае нулевого приближения для внутреннего разложения оно имеет вид:
te0 ( y , ) 1 .
(3.11)
te0 ( y, ) teo y .
(3.12)
Отсюда следует, что
Подставляя (3.12) в уравнение (3.10) в сферической системе координат
имеем:
1
1 d 2 dteo
y
0.
2
y dy
d y
22
(3.13)
Интегрируя уравнение (3.13) два раза, получаем:
te 0 ( y ) D0 0
y
1
1
.
Постоянная интегрирования D0 определяется из условия сращивания
(3.11). В нашем случае, очевидно, что D0 1 .
Таким образом, мы нашли нулевое приближение для внутреннего разложения и оно имеет вид:
te 0 ( y ) 1 0
y
1
1
.
(3.14)
Постоянная интегрирования 0 определяется из граничных условий на
поверхности неравномерно нагретой частицы (1.3) ‒ (1.4).
Найдем первое приближение для внешнего разложения. До членов первого приближения внешнего разложения из (3.4) имеем:
te* , 1 f1* () te*1 (, ) .
(3.15)
Из (3.14) видно, что для нахождения первого приближения внешнего
разложения необходимо сначала определить явный вид коэффициента f1* () .
Для этого в решении (3.14) перейдем к внешней переменной y . Тогда из
(3.14) следует, что f1* () и, следовательно
te* , 1 te*1 (, ) .
(3.16)
Подставляя выражение (3.16) в уравнение (3.7), учитывая (3.8) и удерживая члены порядка , получаем:
te*1 1 x 2 te*1
Pr ( x
) *te*1
x
,
te*1 0
23
(3.17)
(3.18)
В уравнении (3.17) сделаем замену переменной – z Pr и решение
для функции t e*1 будем искать в следующем виде:
1
te*1 f ( z, x ) exp zx .
2
(3.19)
Подставляя (3.19) в уравнение (3.17) и после преобразований получаем
уравнение Гельмгольца [7] для нахождения функции f z, x :
1
f 0.
4
f
(3.20)
Решение уравнения (3.20) будем искать методом разделения переменным (по полиномам Лежандра):
f ( z, x ) f n z Pn x .
(3.21)
n 0
После подстановки (3.21) в уравнение (3.20) и , учитывая свойства ортогональности полиномов Лежандра, получаем:
2 fn
f n
z2
z
2z
n(n 1) f n 0 .
z 2
z
4
2
(3.22 )
Если в уравнении (3.22) перейти к новой переменной z / 2 , то уравнение (3.22) принимает вид:
2 fn
f n
2
2
n(n 1) f n 0 .
2
4
2
(3.23 )
Таким образом, мы получили дифференциальное, решением которого
являются модифицированные сферические функции Бесселя [10]:
I 1 ,
2 n 2
K 1 .
2 n 2
24
(3.24)
Заметим, что при
функция
I 1
2 n 2
I 1 ,
2 n 2
K 1 0 , т.е.
2 n 2
не удовлетворяет граничному условию (3.18).
Общее решение уравнения (3.17), удовлетворяющее граничному условию (3.18) имеет вид:
1
2
Pr
Pr
*
te1 (, )
exp
x
L
K
Pn ( x ) ,
n
1
n
2
2
n 0
Pr
2
1
2
(3.25)
1
2
(n m)!
Pr
Pr n
K 1
exp
2 m 0 (n m)!m!(Pr ) m
Pr n 2 2 Pr
Произвольные постоянные Ln должны быть определены в результате
сращивания, которое в данном случае заключается в сравнении поведения
функции (3.25) при 0 и функции (3.14) при y , т.е. условие сращивания имеет вид – te ( y , ) te* ( 0, ) . Нетрудно установить, что
1
1
0
,
te ( y , ) 1 0 1
(
1
)
Pr
2n ! P ( x ) ,
1
t ( 0, ) 1
( x 1) ... Ln
m n
2
n 0 n!(Pr )
Pr
*
e
и видим: L0
Pr 0
, Ln 0 ( n 1 ). Следовательно,
(1 )
te*1 ,
0
Pr
exp x 1 .
(1 )
2
(3.26)
Найдем первое приближение для внутреннего разложения. Из (3.3)
имеем:
te ( y, ) te0 ( y) f1 ( )te1 ( y, )
25
(3.27)
Для определения явного вида коэффициента f1 ( ) поступим аналогичным образом, что и для f1* ( ) . Видим – f1 ( ) . Таким образом, первое
приближение для внутреннего разложения равно:
te ( y, ) te0 ( y ) te1 ( y, ) .
(3.28)
Таким образом, мы видим, что до первого порядка малости имеем следующие выражения:
te ( y, ) te0 ( y ) te1 ( y, ) ,
(3.29)
ti ( y, ) ti 0 ( y ) ti1 ( y, ) .
(3.30)
Для нахождения локального теплового потока и теплового числа Нуссельта постоянные интегрирования, входящие в выражения для полей температур вне и внутри частицы должны быть определены. Они определяются из
граничных условий на поверхности частицы. С учетом выражений (3.29) и
(3.30) линеаризуем граничные условия на поверхности частицы:
До членов ~
получаем следующие линеаризованные граничные
условия.
− нулевое приближение:
eS te 0 ti 0
RTe3 4
t e 0 ti 0 ,
t 1
0 1
iS y y
iS i 0
(3.31)
− первое приближение:
eS te1 te0 ti1 ti 0
RTe3 3
te1 ti1 ,
t
t
4 0 1
t t
iS y te0 y e1 y ti 0 i1 y
iS i 0 i1
(3.32)
Постоянные интегрирования 0 , B0 , входящие в нулевые приближения определяются из граничных условий на поверхности частицы (3.31) и
определяются из системы уравнений (3.33).
26
Среднее значение температуры поверхности частицы TiS определяется
из решения следующей трансцендентной системы уравнений:
tiS teS ,
S
2
RTe3 4
eS t R
t 1,
J 0 1
1 iS eS 3iS Te 0
iS iS
в котором eS e teS
,
TeS teSTe , J 0
iS i tiS ,
tiS ti 0
y 1
,
(3.33)
teS te 0 y 1 , TiS tiS Te ,
0
1
4 3 S
1
q
r
,
dV
,
,
, 0 teS
1.
V
R
i
VV
1 0
3
В (3.33) интегрирование ведется по всему объему частицы.
Для te1 y, в двучленном внутреннем разложении получаем из (3.1)
(оставляя члены ~ ε) следующую задачу:
Pr
(Ve)te 0 te0 te1 ,
te 0
(3.34)
eS te1 te0 ti1 ti 0
RTe3 3
te1 ti1 ,
t
t
4 0 1
t t . (3.35)
iS y te0 y e1 y ti 0 i1 y
iS i 0 i1
Чтобы определить поведение функции te1 y , , срастим двучленные внутренние и внешние разложения:
te ( y, ) te0 ( y) te1 ( y, ) , te* , 1
0
Pr
exp x 1 ,
(1 )
2
в результате сращивания получаем:
te1 y ,
3
2
( x 1) , 3
0 Pr
.
1
(3.36)
Таким образом, граничное условие для первого приближения поля
температуры вдали от частицы имеет вид (3.36).
27
Из (3.36) видим, что для нахождения te1 y, необходимо сначала
определить поле скорости, т.е. решить гидродинамическую задачу. Как мы
отмечали в первой главе оно считается заданным и имеет вид (1.7) ‒ (1.8).
Выражения (1.7) – (1.8) можно записать в компактном виде:
Vre y, G y cos , Ve y, g y sin .
(3.37)
Поскольку поле скорости нам известно, то мы можем найти первое
приближение для внутреннего разложения. Решение (3.34) ищем в виде:
te1
1
y e ( y) cos ,
te0
(3.38)
с краевыми условиями
y
3
2
y const ,
,
e y
3
2
e y const
при y ,
(3.39)
при y 1 .
(3.40)
Подставляя (3.38) в (3.34), убеждаемся, что переменные разделяются и
y
3 N1
1 ,
2 y
(3.41)
а e y удовлетворяет уравнению
d 2 e
d y2
(1 l )
2 d e e
2 3 2 G y .
ydy y
y
(3.42)
Таким образом, мы получили неоднородное обыкновенное дифференциальное уравнение второго порядка для нахождения функции e y . Интегрируя, получаем:
e y D1 y
1
y2
2
1
A
A
3
2
1 3 ,
3
где
28
y
y
(3.43)
n1l n 1 l 4 n2 l n
0 Pr
;
1 ( y) 1 l
6 n 0 n 4
1
n0 n 1
C 2 l n 2
3l n
1
n 1
2
2
n2
2 ( y)
l 1 l n
1
1 l ln l lC1 l ln l
1 l
n 1 n 2
n0
n0 n 2
2 1 l
ln
1 l 3 l n
1
1
n ln y n
n2 ln y n2
6 n0 n 4
n 0 n 1
y2
1
1 l 2
3 ( y )
6
1
n
l 2 Cn33l n 3
n 1
1
2
2 3
2
l
l
ln
l
l
C
2
l
ln
l
1
2
l
2
2
2
n
1
n
2
n 0
n 1 2
l 1 l
n 3
3
1 l
n4 l n 3 1 l l n
n1 ln y n1
3
y
n 0 n 1
n 0 n 1
ln
n2 ln y n2 ,
n 0 n 4
n
C 1
k
k 0
, n1
n
n 2
k1
2
, n n , n4 Ck3 ,
k 0
k 0 k 1
k 0 k 4
n
n
n
k 0
k 0
n2 n k 1n k 2n k 3Ck1 , n3 n k 1Ck2 .
Поскольку l 1, то при интегрировании использовались следующие
формулы:
1 l
k
k k 1 2
k!
1
1
n
n
n
1 kl
l ...
l ,
l ,
1 l n ,
2
1 l n 0
1 l n 0
n 0 k n ! n!
1 l
k
k k 1 2
k!
1
n
n
n
1 kl
l ... 1
l ,
1 l 2 n ,
2
k n !n! 1 l n 0
2!
n 0
l n 1
ln 1 l 1
,
n 1
n 0
1 l
k
n
k k 1 2
l n 1
n k n 1! n
1 kl
l ... 1
l , ln 1 l
,
k 1!n!
2!
n 0
n 0 n 1
1 l k 1 kl k k 1 l 2 ... k n 1!l n , ln l
2!
n 0 k 1! n!
n 1
n
n l 1
n 1 l 1
1
1
,
n 0
n 1
n
n 1
29
n
n
n
n
Cn x bn x x Ck bn k .
n 0
n 0
n 0 k 0
Следовательно, первое приближение для внутреннего разложения, удовлетворяющее граничному условию (3.39) имеет вид:
te1 y,
1 3 N1 cos 1 3
1 2 3 A2 2 A1 13 .
te 0 2 y
3
y
y
te 0 y
(3.44)
Найдем постоянные интегрирования, входящие в (3.44) и (2.20). Подставляя (4.44) и (2.20) в граничные условия (3.35), имеем:
N 0 4 0 1
RTe3
es
3
tiS
tiS
3
RTe3 3
,
4
t
2
1 , N1
0
1
iS
2
teS
es
es
RTe3 3
1
1 3 es
1
1
1
4
t
A
3
3
0
1
iS
2
2 1 2 1
teS
teS
3 is
is
is
RT 3
RT 3
2 11 4 0 1 e tiS3 A1 es 1 1 31 1 1 11 4 0 1 e tiS3
is
is
is
2 es
RTe3 3
3C1
, 1
4 0 1
tiS .
tiS
is
is
Здесь индекс "S" означает, что физические величины берутся при средней
относительной температуре поверхности частицы TiS , определяемой из выражения (3.33); функции 11 y , 12 y и т.д. берутся при y 1 , а 11' y и
т.д. первая производная от соответствующих функций.
30
ГЛАВА IV. ПОТОК ТЕПЛА НА ПОВЕРХНОСТЬ ЧАСТИЦЫ. АНАЛИЗ
ПОЛУЧЕННЫХ РЕЗУЛЬТАТОВ
В третьей главе получены общие выражения для распределения полей
температур вне и внутри неравномерно нагретой крупной твердой частицы
сферической формы. При получении этих выражений было учтено: степенной вид зависимости коэффициентов молекулярного переноса (вязкости,
теплопроводности) и плотности газообразной среды от температуры и влияние движение среды (учет конвективных членов в уравнении теплопроводности) и имеют следующий вид:
где te 0 ( y ) 1 0
y
te1 y,
1
1
te ( y, ) te0 ( y ) te1 ( y, ) ,
(4.1)
ti ( y, ) ti 0 ( y ) ti1 ( y, ) ,
(4.2)
1 R J 1 Ф d y Ф0 d
, ti 0 y B0
0
y 0
y y
3
T
y
y
i
1
2
1
y
1
1
2
1
1 3 N1 cos 1 3
1
A
A
.
2
1
t y2 3 3
te0 2 y
y
y 3
e0
y
y
N 0 cos
C1
1
1
ti1 y B1 y 2 y 2 dy 2 y1dy .
ti 0
ti 0
y
y
y 1
1
N 0 4 0 1
С1
RTe3
es
R
J1 ,
3 iT
3
tiS
tiS
3
RTe3 3
,
4
t
2
1 , N1
0 1
es iS
teS 2
J1
1
qi z dV , z r cos ,
V V
плотности тепловых источников,
31
q
i
V
z dV дипольный момент
RTe3 3
1
1 3 es
tiS A2 es 2 1 2 1
3 1 3 11 4 0 1
teS teS 3 is
is
is
es
RTe3 3
RTe3 3
2 11 4 0 1
tiS A1 1 1 31 1 1 11 4 0 1
tiS
is
is
is
2 es
RTe3 3
3C1
,
1
4
tiS .
0 1
tiS
is
is
Поскольку выражения для полей температур вне и внутри частицы
нами получены, то мы можем найти локальный тепловой поток. Если температура вязкой среды не постоянна вдоль ее объема, то перенос тепла будет
происходить за счет теплопроводности – непосредственный молекулярный
перенос энергии из мест с более высокой в места с более низкой температурой. Он не связан с макроскопическим движением и происходит также и в
неподвижной среде.
Если обозначить через jT плотность локального потока тепла, переносимого посредством теплопроводности и ограничиться случаем когда градиент температуры в среде не слишком велик, то мы можем разложить jT в ряд
по степеням градиента температуры и ограничиться первыми членами разложения. Постоянный член в этом разложении, очевидно, исчезает, поскольку jT должно обращаться в нуль вместе с T . Таким образом, получаем [3]:
jT T .
(4.3)
В литературе выражение (4.3) называется законом Фурье. Поскольку
локальный поток энергии должен быть направлен из мест с более высокой в
места с более низкой температурой jT и T должны иметь противоположные направления.
Найдем jT для нашего случая. Поскольку мы рассматриваем теплообмен при больших перепадах температуры, т.е. коэффициент теплопроводности не постоянная величина, а является функцией температуры имеем:
e e te e te0 te1 e te0 te01te1 ,
32
jT e Te e
T te
1 te
er
e
R y
y .
До членов порядка получаем следующее выражение для локального
потока тепла:
T te0
jT e
R
t e1
te1 t e 0
1 t e1
t e0
e
e
e
.
r
r
y
y
t
y
y
e0
(4.4)
Далее нам потребуются:
te 0
te 0
t
, e1 3
1 y y
y
2 yt e 0
N1
1 y 1 1
3
A2 ' 2
cos 1
3 '
2
1
3
2
te0 y 3 1
y 1 y
y 1
3
A1 ' 1
3
,
1
3
y
y 1
te1
sin
12 3 3 A2 2 A1 33 ,
te 0 y
3
y
y
и подставляя полученные выражения в (4.4), получаем следующее выражение для локального потока тепла:
jT e
2
T te01
T
A
er e 32 N1 cos 3 1 3 '3 2
R y 1
R
3
y
y
2 y
3
1
2
sin 1 3
' 2 A1 '
2 3 1 3 er
2 3 A2 A1 3 e . (4.5)
y y
y
y y
3
y
y
Зная локальный поток тепла можно найти: молекулярный поток тепла
QT M и тепловое число Нуссельта NuT :
QT M
j
T
dS P ,
(4.6)
Sp
1
t
NuT e d x .
y y 1
1
33
(4.7)
где dS P дифференциальный векторный элемент поверхности частицы S P ,
направление которого совпадает с направлением внешней нормали.
Из полученных выражений видно, на физические величины, которые
характеризуют теплообмен существенное влияние оказывает нагрев поверхности частицы, зависимость коэффициентов молекулярного переноса (вязкости, теплопроводности) и плотности газообразной среды от температуры и
движение среды (учет конвективных членов в уравнении теплопроводности).
Следует отметить, что в полученные выше выражения входят постоянные интегрирования A1 и A2 . Эти постоянные входят в поле скорости, которое в работе считается известным из решения соответствующей гидродинамической задачи. Это позволяет рассматривать процесс теплообмена для общего случая, т.е. для установившегося движения крупной твердой неравномерно нагретой частицы сферической формы. Установившееся движение
означает, что на частицу кроме силы вязкого сопротивления среды (силы
Стокса) действует другая сила, направленная противоположно силе сопротивления среды. Это может быть силы тяжести (гравитационное осаждение),
термофоретическая сила (термофорез), фотофоретическая сила (фотофорез) и
т.д.
34
ЗАКЛЮЧЕНИЕ
В магистерской диссертации получены аналитические выражения для
полей температур вне и внутри крупной неравномерно нагретой частицы с
учетом степенного вида зависимости коэффициентов молекулярного переноса (вязкости, теплопроводности) и плотности газообразной среды от температуры и влияния движении среды. При нахождении распределения температуры в окрестности частицы использовался метод сращиваемых асимптотических выражений; выражение для локального теплового потока тепла, которое позволяет найти физические величины, характеризующие теплообмен, а
именно: молекулярный поток тепла QT M и тепловое число Нуссельта NuT .
35
СПИСОК ИСПОЛЬЗОВАННОЙ ЛИТЕРАТУРЫ
1. Ю.И. Яламов, В.С. Галоян Динамика капель в неоднородных вязких
средах. Ереван.: Луйс. 1985. 207 с.
2. Ст. Бретшнайдер Свойства газов и жидкостей. Инженерные методы
расчета. М.: Химия. 1966. 535 с.
3. Л.Д, Ландау, Е.М. Лифшиц Гидродинамика. М.: Наука. 1986. 736 с.
4. Дж. Хаппель, Г. Бреннер Гидродинамика при малых числах Рейнольдса. М.: Мир. 1976. 630 с.
5. Малай Н.В., Щукин Е.Р., Стукалов А.А., Рязанов К.С. Гравитационное
движение равномерно нагретой твердой частицы в газообразной среде
// ПМТФ. 2008. № 1. С. 74 – 80.
6. Малай Н.В., Лиманская А.В, Щукин Е.Р. Термофоретическое движение
нагретых крупных аэрозольных частиц сферической формы //ПМТФ.
2016. Т. 57, № 2(336). С. 164-171.
7. А.Н. Тихонов, А.А. Самарский Уравнения математической физики. М.:
Наука. 1972. 735 с.
8. Матвеев А.Н. Методы интегрирования обыкновенных дифференциальных уравнений. М.: Высшая школа. 1967. 564 с.
9. Найфе А. Введение в методы возмущения. М.: Мир. 1984. 535 с.
10. Абрамовица М., Стиган И. Справочник по специальным функциям.
М.: Наука. 1979. 830 с.
36
Отзывы:
Авторизуйтесь, чтобы оставить отзыв