Сохрани и опубликуйсвоё исследование
О проекте | Cоглашение | Партнёры
Дипломная работа 03.03.02 Физика
Источник: Белгородский государственный университет - национальный исследовательский университет (НИУ «БелГУ»)
Комментировать 0
Рецензировать 0
Скачать - 763,4 КБ
Enter the password to open this PDF file:
-
ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ АВТОНОМНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ОБРАЗОВАНИЯ «БЕЛГОРОДСКИЙ ГОСУДАРСТВЕННЫЙ НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ УНИВЕРСИТЕТ» ( Н И У « Б е л Г У » ) ИНСТИТУТ ИНЖЕНЕРНЫХ ТЕХНОЛОГИЙ И ЕСТЕСТВЕННЫХ НАУК Кафедра теоретической и математической физики ТЕРМОФОРЕЗ АЭРОЗОЛЬНОЙ ЧАСТИЦЫ ЦИЛИНДРИЧЕСКОЙ ФОРМЫ Дипломная работа студента очной формы обучения направления подготовки 03.03.02 Физика 4 курса группы 07001210 Лифинцев Иван Сергеевич Научный руководитель: д.ф.-м.н., профессор Малай Н.В. БЕЛГОРОД - 2016
Реферат Рассмотрен термофорез крупной твердой аэрозольной частицы цилиндрической формы, внутри которой действуют неоднородно распределенные по ее объему тепловые источники. Получены аналитические выражения для полей температуры вне и внутри аэрозольной частицы, а также выражения для силы и скорости термофореза цилиндрической частицы.
Содержание. Введение.………………………………………………………………….…….4 I. Постановка задачи. …………………………………………………….…….9 II. Решение уравнения теплопроводности. ……………………………….…16 III. Решение уравнений гидродинамики. Нахождение полей скорости и давления…………………………………………………………….…….. …...23 IV. Вывод выражения для силы термофореза. ……………………………...34 V. Анализ полученныхрезультатов. …………………………………………36 Заключение. …………………………………………………………………...39 Список использованной литературы………………………………………....40 3
ВВЕДЕНИЕ Актуальность исследования.В дипломной работе исследуется поведение аэрозольных частиц несферической формы в вязких неоднородных газообразных средах, внутри которых действуют неравномерно распределенные источники тепла. В настоящее время имеет место огромный интерес установление комплекса физических закономерностей, играющих решающую роль в предсказании характера поведения аэрозольных частиц в газах и газовых смесях. Это обстоятельство неслучайно, так как с каждым годом увеличивается использование аэрозолей в практических целях– промышленности, технике, сельском хозяйстве, медицине и т.д. Образующиеся в результате производственной деятельности человека аэрозоли могут, с одной стороны, содержать ценные вещества, с другой,–оказывать вредное влияние на людей и окружающую среду. В связи с обострением экологической ситуации необходимо решать вопросы очистки промышленных отходов от аэрозольных частиц, природа образования которых может быть произвольной. Однокомпонентную или многокомпонентную газообразную среду, вместе со взвешенными в ней частицами называют аэродисперсной системой. Аэродисперсные системы находят широкое применение в практической деятельности человека, например, в технике, медицине, сельском хозяйстве и т.д. В связи с этим всестороннее изучение их свойств представляет значительный научный и практический интерес. На входящие в состав аэродисперсных систем аэрозольные частицы могут действовать силы различной природы, вызывающие упорядоченное движение частиц относительно центра инерции газообразной среды. Так, например, седиментация и электрофорез происходят в поле гравитации и электрической силы. В газообразных средах с неоднородным распределением температуры возникновение упорядоченного движения частиц обусловлено действием сил молекулярного происхождения. Их появление вызвано передачей некомпенсированного импульса частицам молекулами газообразной среды. При этом 4
движение частиц, обусловленное внешним градиентом температуры, называется термофорезом. Неоднородное распределение температуры в окрестности частиц может возникнуть при их нагреве и охлаждении источниками или стоками тепла произвольной физической или химической природы. Возникающее при этом движение частиц называется теплофоретическим. Частным случаем теплофоретического движения является фотофорез. Фотофоретическое движение возникает при нагреве частиц световым излучением. При неоднородном распределении температуры движение аэрозольных частиц происходит, например, в каналах тепло – и массообменников, в зонах просветления облаков и туманов, устройствах предназначенных для нанесения тонких покрытий из аэрозольных частиц, химических реакторах и т.д. Аэрозольные частицы могут оказывать значительное влияние на протекание физических и физико-химических процессов различного вида в аэродисперсных системах (например, процессов массо – и теплообмена). Это обстоятельство нужно учитывать при нахождении физических величин внутри рассматриваемых систем. Проблема вынужденного переноса аэрозольных частиц под действием градиентов температуры (термофорез) в газах имеет большое научное и прикладное значение. Явление термофореза всегда сопутствует термодинамически неравновесным системам. Термофоретический перенос аэрозольных частиц играет важную роль в установках по улавливанию аэрозольных примесей; для нанесения тонких покрытий; для описания атмосферных процессов, при получении веществ с заданными свойствами и ряде других областей деятельности человека [1,2]. Твердые и жидкие частицы, находящиеся во взвешенном состоянии в газообразной среде, называются аэрозольными частицами. Для классификации аэрозольных частиц по размерам применяют критерий Кнудсена: Kn R , где средняядлинасвободного пробега молекул газообразной среды; R радиус частицы. 5
Частицы называют крупными, если Kn 0.01, умеренно крупными при 0.01 Kn 0.3 и мелкими при Kn 1. Если на поверхности вещества происходит испарение или конденсация частиц вещества, то такие частицы называются летучими. Нелетучими называются частицы, на поверхности которых фазовый переход не происходит В настоящей дипломной работе рассматриваются твердые крупные частицы. Среднее расстояние между аэрозольными частицами у значительной части встречающейся на практике аэродисперсных систем намного больше характерного размера частиц. В таких условиях учет влияния аэрозоля на развитие физического процесса можно проводить, основываясь на знание законов динамики движения и тепло- и массообмена с бесконечной окружающей средой отдельных аэрозольных частиц. Поэтому изучение закономерностей движения отдельных частиц в газообразных средах является важной актуальной задачей, представляющей значительный теоретический и практический интерес. Рассмотрение термофоретического движения аэрозольных частиц является достаточно сложной задачей. Это связано с тем, что в газе движение конкретной частицы определяется как поверхностными явлениями (обусловленные непосредственным взаимодействием молекул газообразной среды с поверхностью частицы), так и с объемными эффектами, возникающими из-за неоднородных распределений гидродинамического и температурного полей. Одна из задач, возникающая при практических приложениях вызвана тем, что многие частицы, встречающиеся в промышленных установках и в природе, имеют форму поверхности отличную от сферической, например, цилиндрическую. В связи c этим возникает как теоретический, так и практический интерес раccмотрение термофоретического движения цилиндрической аэрозольной частицы. 6
Тема исследования. Термофорез твердой крупной цилиндрической аэрозольной чаcтицы, внутри которой действуют неравномерно распределенные в ее объеме источники тепла. Объектом исследования является изучение явления термофореза в газообразной среде. Предметом исследования:влияние формы поверхности и внутренних источников тепла, неоднородно распределенных в объеме частицы на термофорез крупной аэрозольной частицы цилиндрической формы. Цель исследования –получить аналитические выражения для термофоретическойсилы, действующей на твердую крупную аэрозольную частицу цилиндрической формы и скорости термофореза с учетом внутренних источников тепла. Исходя из поставленной цели, были сформулированы следующие задачи исследования: - ознакомление с историей развития термофореза цилиндрических частиц; - изучить математические методы решения дифференциальных уравнений второго порядка в частных производных в цилиндрической системе координат; - получить выражения для поля температуры внутри частицы, с учетом тепловых источников неоднородно распределенных в ее объеме; - получить выражения для силы и скорости термофореза с учетом внутренних источников тепла; - провести качественный анализ влияния на силу и скорость термофореза внутренних источников тепла. Научная новизна исследования. В настоящей дипломной работе делается попытка проанализировать явление термофореза в вязкой газообразной среде для твердой крупной аэрозольной частицы цилиндрической формы с учетом внутренних источников тепла. 7
Практическая значимость исследованиязаключается в том, что его материалы и выводы дополняют и углубляют исследования по данной проблеме и могут быть использованы при разработке общих курсов и факультативов по экологическим вопросам для школ, гимназий, лицеев. Апробация исследования. Основные результаты выпускной квалификационной работы докладывались на Международной научно-практической и научно-методической конференции «Современные проблемы математики и механики: теория и практика» (г. Белгород, БУКЭиП, 2016) и по теме выпускной квалификационной работы опубликованы тезисы. Структура работы. Дипломная работа состоит из введения, пяти глав, заключения, списка литературы. 8
Глава I. Постановка задачи Пусть в неограниченной неподвижной изотропной бинарной газовой смеси создан и поддерживается с помощью внешних источников постоянный градиент температуры T . В такую неоднородно нагретую среду помещаетсякрупнаяцилиндрическая частица радиуса R . Причем R значительно меньшего ее длины L . Предположим, что ось частицы перпендикулярна направлению градиента температуры T . Свяжем систему отсчета с частицей, причем так, чтобы ось z совпадала с осью частицы, а полярная ось x совпадала с направлением T (см. рис.1). При таком выборе системы координат наша задача сводится к анализу обтекания частицы плоскопараллельным потоком газа со скоростью U . Определенная в такой системе координат скорость газа на бесконечности, равна с обратным знаком величине скорости термофореза, U th U 9
Здесь и далее индексы «е» и «i» будем относить соответственно к газообразной среде и частице; индексом « » – обозначены параметры газа на бесконечности, т.е. вдали от частицы. Общая система уравнений газовой динамики нелинейная [3-4], и при решении поставленной выше задачи были сделаны следующие предположения: 1. все процессы, происходящие в системе, будем рассматривать в квазистационарном приближении, что возможно в силу малости времен тепловой и диффузионной релаксации умеренно крупных и крупных аэрозольных частиц; 2. движение аэрозольных частиц во внешнем заданном поле градиента температуры будем считать происходящим при числах Пекле и Рейнольдса много меньших единицы; 3. частица образована однородным и изотропным по своим свойством веществом. 4. задача рассматривается при малых относительных перепадах температуры в окрестности частицы, т.е. выполняется следующее условие Te TS / T 1. Здесь Te температура газа, TS средняя температура поверхности частицы и T температура газообразной среды вдали от частицы. При выполнении этого условия коэффициенты молекулярного переноса (вязкости и теплопроводности) можно считать постоянными величинами. При введенных выше предположениях решалась следующая система газодинамических уравнений, описывающая распределения полей скорости Ve , давления Pe и температур вне Te и внутри Ti частицы [4]: e Ve Ve Pe e Ve (1.1) div Ve 0 (1.2) Te 0 (1.3) 10
Ti qi (1.4) i Запишем граничные условия, удовлетворяющие нашей задачи в полярной системе координат. – на поверхности частицы (при r R ) справедливы граничные условия (1.6) – (1.9): Vr 0 , V e (1.6) K TS e Te , (1.7) Te e R Te Ti i , r r Te Ti , (1.8) (1.9) – в дали от цилиндрической частицы (при r ) они имеют вид (1.10)- (1.13) : Vr U cos , (1.10) V U sin , (1.11) Te T | T | r cos , Pe P . (1.12) (1.13) – и внутри частицы при r 0 , имеем: i , (1.14) где e , e , e плотность, теплопроводность и динамическая вязкость газа; i – теплопроводность частицы; r , полярные координаты; Vr , V компоненты вектора скорости в полярной системе координат; | T | величина заданного градиента температуры, KTS коэффициент теплового скольжения, q i плотность тепловых источников, неоднородно распределенных в объеме частицы. 11
Наличие источников тепла внутри частицы приводит к тому, что средняя температура ее поверхности TS может отличаться от температуры окружающей газообразной среды вдали от нее ( T ). Нагрев поверхности частицы оказывает влияние на теплофизические характеристики окружающей газообразной среды (вязкости, теплопроводности, плотности) и в конечном итоге на распределение полей скорости, давления и температуры в окрестности. В нашем случае это влияет на термофорез. Внутренние источники тепла − это модельное представление, удобное для описания реальных процессов, сопровождающихся выделением тепла в объеме аэрозольной частицы. Так можно моделировать нагрев поверхности под действием химической реакции, радиоактивного распада вещества частицы, поглощением электромагнитного излучения и т.д. Поскольку поглощенная частицей энергия превращается в тепловую энергию (идет на нагрев частицы), то обозначим через qi ( r ) плотность тепловыделения в объеме частицы и будем называть qi ( r ) плотностью тепловых источников, неоднородно распределенных в ее объеме. Считается, что эта функция задана (см. ниже). Если, например, неоднородный нагрев обусловлен поглощением электромагнитного излучения, то qi x , y , z определяется из решения электродинамической задачи. Электродинамическая задача, позволяет рассчитать функцию источников тепла в объеме частицы при заданных значениях длины волны, интенсивности излучения, размеров частицы и комплексного показателя преломления[5,6]. При макроскопическом рассмотрении объемная плотность источников тепла равняется с обратным знаком дивергенции вектора плотности потока излучения I 0 (вектора Умова – Пойтинга) [7]: qi ( r )= I 0 , 12
где I 0 c Re E , H , E , H – векторы напряженности электрического и 4 магнитного поля внутри частицы, которые описываются уравнениями Максвелла: , E i k H , , H i m2k E . Здесь k 2 / 0 волновое число, 0 длина волны падающего излучения, m n1 i k1 комплексный показатель преломления вещества частицы, n1 , k1 действительная и мнимая части комплексного показателя преломления вещества частицы. В результате функция qi r принимает следующий вид [5,7]: qi r 2 kk1I 0 Br , где Br E r / E02 безразмерная функция, описывающая распределение 2 электрического поля в объеме частицы, E0 модуль напряженности электрического поля падающей волны. Функция Br может быть рассчитана на основе строгой электродинамической теории Ми [5]. Следует отметить, что возможностьсуществования неоднородности распределения поглощенной энергии излучения в частице зависит от оптических констант материала частицы и параметра дифракции 2 R0 / 0 , R0 эффективный радиус частицы. Таким образом, выражение для плотности энергии излучения в частице, трансформируемой в тепло, можно записать в виде [5,7]: qi 4 n1k1 I0B , ns где ns показатель преломления среды. На поверхности аэрозольной частицы (1.6) – (1.9) учтено: непроницаемость поверхности частицы для нормальной компоненты и тепловое скольжения для касательной компоненты массовой скорости Ve ; непрерывность радиального потока тепла. Вдали от частицы ( r ) справедливы граничные 13
условия (1.10) – (1.13) и конечность физических величин, характеризующих аэрозольную частицу при r 0 учтено в (1.14). Коэффициент K TS определяются методами кинетической теории газов(из решения уравнения Больцмана в слое Кнудсена[8]). При коэффициентах аккомодации тангенциального импульса и энергии, равных единицы газокинетические коэффициенты (в случае плоской поверхности) равен соответственно: KTS 1.152 [8]. Определяющими параметрами задачи являются коэффициенты e , e , e и сохраняющиеся в процессе движения цилиндрической частицы величины – R, T , T . Из этих параметров можно составить три безразмерные комби- нации: R T / T 1 – параметр, характеризующий перепад температуры на размере частицы, числа Рейнольдса Re e U R / e 1 и Пекле Pe c peU R e / e 1 , c pe удельная теплоемкость при постоянном давле- нии. Обезразмерим уравнения и граничные условия следующим образом: U e Ve / U , tk Tk / T ( k e, i ), pe Pe / P . При 1 набегающий поток оказывает лишь возмущающее влияние и поэтому решение уравнений газовой динамики можно использовать метод теории возмущения, т.е. U e U e0 U e1 ..., pe pe0 pe1 ..., t t 0 t 1 ... (1.15) При нахождении силы и скорости термофореза мы ограничимся поправками первого порядка малости по . Общая сила, действующая на единицу длины аэрозольной цилиндрической частицы, определяется интегрированием тензора вязких напряжений по поверхности частицы и имеет следующий вид [3,4] 2 F R prr cos pr sin d , (1.16) 0 14
где prr , pr компоненты тензора вязких напряжений и в цилиндрической системе координат равны: p rr Pe 2e V 1 Vr V Vr , pr e . r r r r 15
Глава II.Решение уравнения теплопроводности Чтобы найти силу и скорость термофореза необходимо знать поле температуры в окрестности аэрозольной частицы. Для этого необходимо решить уравнения (1.3) – (1.4). В общем случае уравнение теплопереноса имеет следующий вид [4]: Te (Ve )Te , (2.1) где коэффициент температуропроводности[4]. Скорость аэрозольной частицы в поле градиента температуры по порядку величины равна e T e | T | . eTe L0 e Te Здесь (2.2) T средний градиент температуры вдали от частицы. L0 Этот градиент равен отношению перепада температуры T на конеч- ном отрезке L0 к величине этого отрезка. Скорость термофореза равна взятой с обратным знаком скорости центра инерции среды на большом расстоянии от частицы, поэтому | U th |~ e | T | e Te (2.3) Отношение правой части уравнения (2.1) к его левой части равно | UT Te | R R ~ | T | ~ | T | . Te Te Te В (2.4) учтено, что для большинства газов (2.4) e ~ 1. e Покажем, что для реальных условий величина R | T | 1. To Найдем, при каких значениях | T | нарушается последнее соотношение. Для этого запишем: 16
| T | 1 . Te R (2.5) Верхняя граница твердых гидрозольных частиц порядка 300 микрон. Отсюда для | T | следует: Te 3 102 K K . | T | 104 6 R 3 10 м см (2.6) Такой градиент температуры никогда не реализуется в промышленных и природных условиях. Не наблюдается даже градиенты температуры порядка 3000 К . см Следовательно, в уравнении (2.1) можно пренебречь правой нелинейной частью, и мы получаем следующие уравнения, описывающие распределения полей температур внутри и вне крупной аэрозольной частицы цилиндрической формы: Ti qi i , Te 0 . (2.7) Найдем сначала общее решение уравнения t 0 (2.8) в полярной системе координат ( y, ) . Уравнение (2.8) решаем методом разделение переменной, т.е. t ( y, ) R ( y) F( ) , или в развернутом виде 1 t ( y, ) 1 2 t ( y, ) y 0 y y y y2 2 Подставляя сюда t ( y, ) R ( y) F( ) , разделяя переменные (разделим за- y2 ) и, обозначив константу разделения через n 2 , получаем сиRF стему из двух уравнений: тем все на 17
1 d 2F n2 ; 2 F d y d dR (y ) n2 ; R dy dy Решая уравнение : d 2F Fn2 0. 2 d методом Эйлера, получаем следующее решение: F ( ) An cos(n ) Bn sin( n ), гдеn=0,1,2,3… Найдем решение уравнения y d dR (y ) n2 ; R dy dy Здесь необходимо рассмотреть два случая. а).при n=0 оно принимает следующий вид: y d dR (y ) 0; dy dy иего решением является функция: R( y) N 0 ln( y) D. б).при n>1 имеем: y2 d2 R dR y n2 R 0 2 dy dy Решая это уравнение методом Эйлера, получаем: R( y ) (N n 0 n yn Dn ). yn Таким образом, общее решение уравнения Лапласа в полярных системах координат ∆𝑡 = 0 имеет вид: 18
t ( y , ) N 0 ln( y ) D0 (Nn yn n 0 Dn )( An cos(n ) Bn sin( n ) .) yn Константы интегрирования N0, D0, Nn, Dn, An, Bn – определяются из граничных условий. В частности, мы рассматриваем термофорез цилиндрической частицы с радиусом R и для граничного условия вдали от частицы Te T cos ( ) получаем: t e 1 cos( )( где N1 Te R D 1) T y (2.9) Te R ; 𝑁𝑛 = 0, (𝑛 ≥ 2), N0=0; D0=1; An=0, n≥2; Bn=0, T n≥1. Найдем температуру внутри частицы, для этого необходимо решить следующее уравнение: Ti qi i (2.10) или в развернутом виде: ti qi R 2 1 1 2 ti (y ) 2 . (2.11) y y y y 2 i T Это есть неоднородное дифференциальное уравнение в частных производных второго порядка и его решение будем искать метод разделения переменных. Для этого сначала правую часть уравнения (2.11) представим в виде: qi R 2 q0 ( y ) q1 ( y ) cos( ). i T (2.12) а).найдем сначала 𝑞0 . Для этого умножим (2.12) на 𝑑ф и проинтегрируем от нуля 0 до 2𝜋 , получаем: 19
2 qi R 2 0 ( i T )d q0 ( y )2 , и, следовательно, 2 R2 q0 ( y ) ( ) qi ( y )d 2 iTi 0 б).найдем 𝑞1 (𝑦). Поступая аналогичным образом, получаем 2 qi R 2 0 ( i T ) cos( )d 2 q 1 2 q 0 ( y ) cos( ) d 0 ( y ) cos2 ( ) d , 0 qi R 2 i T 2 2 1 cos 2 d 2 0 qi ( y ) cos( )d q1 ( y ) 0 q1 ( 2 ) q1 , 2 2 R2 q1 ( y ) q ( y ) cos( ) d . (2.13) iT 0 i Поскольку функции 𝑞0 (𝑦)и 𝑞1 (𝑦) нами получены, то мы можем уже и решить само уравнение (2.11): ti 1 1 2 ti (y ) 2 q0 ( y ) q1 ( y ) cos( ). y y y y 2 Подставляем сюда выражение ti ( y, ) ti0 ( y ) ti1 ( y ) cos( ). (2.14) и после сокращения на cos(ф) получаем систему из двух уравнений: dti0 1 d y ( ) q0 ; y dy dy dti1 ti1 1 d y ( ) q1 ; 2 y dy dy y 20
Таким образом, мы получили два неоднородных уравнения, которые мы будем решать в виде: t n ( y ) t nодн t nчаст . (2.15) Заметим, что общим решением обыкновенного неоднородного дифференциального уравнения второго порядка [10,11] y '' p(x) y ' q(x) y f (x) , является функция x y C11 C 2 2 2 x0 1f W x dx 1 x0 2f W dx , где 1 , 2 функции, являющие решением однородного уравнения, а W определитель Вронского. С учетом выше сказанного, получаем: 0 tio ( y ) M 0 ln( y ) M 1 ln( y ) 0 ydy 1 0 (2.16) 0 ( y ) ln( y ) dy , 1 Постоянные интегрирования определяются из граничных условий. Таким образом, мы нашли общие решенияуравнений теплопроводности вне и внутри частицы: q Ti i , Te 0 , i которые удовлетворяют граничным условиям конечности температуры внутри частицы и вдали от частицы (1.12), (1.14): T D y 1 ), (2.18) T y ti ( y, ) ti0 ( y ) ti1 ( y ) cos( ), (2.19) t e ( y, ) 1 cos( )( где 21
0 tio ( y ) M 0 N 0 ln( y ) ln( y ) 0 ( y ) dy 1 0 0 ln( y ) dy , 1 K 1 1 ti1 ( y ) Гy y 1 dy y 2 y y 1 y y y dy , 1 1 0 0 1 1 0 yqi , 1 yq1 , N 0 0 dy, K 1 ydy , R2 q0 ( y ) 2 iT q1 ( y ) 2 R iT 2 q i cos( ) d , 0 2 q i cos( )d . 0 Оставшиеся постоянные интегрирования определяем из граничных условий на поверхности частицы и они имеют следующий вид: e 1 T R i 2 K , M 0 1 , N0 0 , D1 e e T 1 1 i i e T R e i 1 Г 2 K T i 1 e 1 e i i 1 22
Глава III. Решение уравнений гидродинамики. Нахождение полей скорости и давления Наибольшую трудность представляет решение уравнения гидродинамики – уравнение Навье-Стокса. e Ve Ve Pe e Ve (3.1) В это уравнение входит нелинейный (конвективный) член Ve Ve . Если число Рейнольдса много меньше единицы, то этот член квадратичен по скорости и Стокс пренебрег им при решении задачи обтекания твердой частицы сферической формы [12,13,14]. В литературе такой способ решения уравнения Навье-Стокса получил название «метод Стокса». Однако особенностью обтекания цилиндрической частицы в гидродинамике заключается в том, что обычный метод решения уравнений гидродинамики (разложение полей скорости и давления в ряды по малому параметру как это сделал Стокс) не проходит. В литературе это получило название «парадокс Стокса»[12,13,14]. Специально для случая обтекания цилиндрической частицы в 1910 году Озееном был разработан новый метод решения уравнений гидродинамики, который в последствии в литературе был назван «метод Озеена» [12,13,14]. Исследования показали, что если в задаче об обтекании кругового цилиндра плоскопараллельным потоком вязкой жидкостью в основных уравнениях отбросить полностью инерционные члены, то задача не имеет решения. Форма цилиндра при этом не имеет никакого значения. Высказанное утверждение, справедливо для цилиндра произвольной формы. Озеен показал на примере движения сферы в вязкой жидкости, что мы получим гораздо лучшие результаты, если в уравнениях движения оставим только важнейшие из инерционных членов, отбросив остальные инерционные члены. 23
Рассмотрим задачу об обтекании цилиндром потоком, имеющим на бесконечности скорость, параллельную оси Ox и равную U .Озеен предложил заменить скорость V на U , т.е. в далеких от цилиндра точках имеем Vx U Vx' , Vy U Vy' , где Vx' , Vy' малые величины. Рассматривая теперь инерционные члены, стоящие в левой части уравнения (3.1) мы увидим, что они отличаются от U Vx , x U Vy y , (3.2) малыми членами второго порядка, если считать Vx' , Vy' и инерционные члены за малые члены первого порядка. Поэтому мы получим гораздо лучшее приближение в далеких от цилиндра областях, если заменим уравнения (3.1) следующими обобщенными уравнениями Стокса: (уравнения Озеена) Vx 1 Pe e Vx , x e x Vy 1 Pe U e V y , (3.3) x e y Vx V y 0, x y U где e e / e кинематическая вязкость газовой среды. Уравнения (3.3) имеют очень простой вид, если их записать в векторной форме: U V 1 grad Pe e Ve , x e div Ve 0 (3.3) Нужно отметить, что для области газа, непосредственно примыкающей к цилиндру, замена инерционных членов величинами (3.2) ничуть не лучше за- 24
мены этих членов нулями, так как в этой области Vx , Vy малы (на самой поверхности Vx , Vy обращаются в нуль), и мы не можем использовать факт малости Vx' по сравнению с U . Однако, поскольку мы рассматриваем движения с малыми числами Рейнольдса, то, как полные инерционные члены, так и заменяющие их в наших уравнениях выражения (3.2) будут малы по сравнению с членами, происходящими от сил вязкости. Следовательно, в области, примыкающей к цилиндру, уравнения (3.3) и уравнения Стокса являются в одинаковой мере хорошими приближениями к полной системе дифференциальных уравнений (3.1). Граничные условия, которыми должно удовлетворять решение системы (3.2): Vx Vy 0 при r R Vx U , Vy 0 при r (3.4) Для решения этой задачи мы применим метод Ламба[12,13,14]. Образуем, расхождение и вихрь от обеих частей первого из уравнений (3.3), принимая во внимание второе из этих уравнений, придем к формулам: Pe 0 , (3.5) U (3.6) rot V . где , x Используем теперь симметрию движения относительно оси Ox . x 0 . (3.7) тогда условие div x y z 0 x y z приводит к равенству y y z 0, z 25
откуда следует, что y , z z . y (3.8) И так мы имеем формулы: Vy Vz 0, y z Vx Vz , z x z Vy x Vx . y y Сразу видно, что частным решением этой системы является V x , Vy 0 , и общим же решением будет Vx , x Vy 0. (3.9) Подстановка значений (3.8) в уравнение (3.6) приводит к равенствам, которые будут удовлетворены, если e . x U С другой стороны, подставляя значения (3.9) в уравнение неразрывности, находим: x и, на основании предыдущего равенства, e U , откуда видно, что следует принять e U , где удовлетворяет уравнению Лапласа 0 . (3.10) Введем для краткости обозначение k U / 2 e , тогда выражения для составляющих скорости примут следующий окончательный вид: 26
Vx 1 1 . ; Vy x 2k x y 2k y (3.11) причем функция χ соответствует уравнению 2k 0. x (3.12) Подставляя значения (3.11) в уравнения (3.2),можем решить эти уравнения относительно Pe , в результате получим: Pe P eU . x Уравнению (3.10) удовлетворяет функция lnr, а также (3.13) ln r 2 ln r , и т.д. 2 x x Поэтому мы принимаем A0 ln r A1 ln r 2 ln r A2 ... x x 2 (3.14) В уравнении (3.12) делаем стандартную подстановку e kx , (3.15) тогда получаем для определения ψ уравнение типа Гельмгольца 2 2 k 2 0 2 2 x y имеющее в полярных координатах (r, θ) вид: 2 1 1 2 2 k 2 0 2 2 r r r r (3.16) Отыскиваем те решение этого уравнения, которое зависит только от r и, следовательно, удовлетворяет обыкновенному дифференциальному уравнению 2 1 k 2 0 2 r r r (3.17) 27
Решениями этого последнего уравнения являются функции Бесселя с чисто мнимым аргументом[10,11]: I 0 (kr) J 0 (ikr) и K 0 (kr) i 2 H 01 (ikr). Однако I 0 (kr) безгранично возрастает вместе с r, в то время как K 0 (kr) стремится к 0 при r ; поэтому единственным приемлемым для нас реше2 нием уравнения (3.17) является функция K 0 (kr) .Функции K 0 k r , K 0 2k r x x ,…также будут решениями уравнения (3.15). Мы можем по этому принять: kx K 0 (kr) 2 K 0 (kr) U e C0 K 0 (kr) C1 C2 ... . 2 x x (3.18) Так как ln r x cos 2 , x r r 2 ln r cos 2 ,…, 2 x r2 то для φ находим следующее разложение: A0 ln r A1 cos cos 2 A2 ... r r2 (3.19) С другой стороны, мы имеем при малых kr разложения: I 0 (kr) 1 1 (kr) 4 2 ( kr ) ... 22 (2 4) 2 2 2 1 k 4r 4 1 k r K 0 (kr) I 0 (kr) ln kr 1 ..., 2 (2 4) 2 22 2 где - постояннаяМаскароны ( =1,7811; ln =0,57722). Ограничиваясь самым первым приближением, мы можем принять, что 1 K 0 ( kr) ln kr , 2 K 0 (kr) 1x cos ; x r r r поэтому вблизи цилиндра считая k r малым, будем иметь приближенно 28
U C0 ln kr kr cos ln kr . C1 1 2 1 2 cos . r (3.20) Теперь мы можем вычислить по формулам 1 V cos r r 2 r (3.21) 1 1 sin V r 2 kr Найдем частные производные от функций и по r и по : A A A 0 1 cos 2 32 cos 2 ... r r r r (3.22) A A 1 sin 2 22 sin 2 ... r r (3.23) C 1 1 C1 cos 0 U C0 ln kr kr cos ln kr r r r r 2 2 1 C C0 k cos ln kr 21 cos 2 r 3.24) C1 1 1 C1 U C ln kr kr cos ln kr cos sin 0 2 2 r r 1 C0 kr sin ln kr 2 (3.25) Подставляя (3.22), (3.23), (3.24), (3.25) в систему уравнений (3.21), получим: Vr A0 A1 A C C C 1 C 2 cos 2 32 cos 0 0 cos ln kr 0 cos 1 2 r r r 2kr 2 2kr 2 2 C 1 1 U cos C0 ln kr kr cos ln kr cos 1 cos2 r 2 2 29
V A1 A C 1 C sin 2 32 sin 2 0 kr sin ln kr 21 sin U sin r r 2kr 2 r 1 1 C C0 ln kr kr cos ln sin 1 cos sin r 2kr 2 при вычислениях мы отбрасываем некоторые члены, малые в сравнении с оставляемыми, а именно, вычисляя cos и sin , мы пользуемся более упрощенным выражением для : 1 2 U C0 ln kr ; кроме того, в выражении для берутся только два первых члена, тогда получаем окончательные выражения для Vr и V : A0 C0 C C 1 A C U 0 ln kr 21 0 1 2 cos Vr r 2kr 2 2 2 2kr r (3.26) C1 A1 C0 1 V 2 ln kr U sin 2 r 2 2 2 kr Полагая в полученных формулах r = R, составляя равенства Vr ( R, ) V ( R, ) 0 и приравнивая 0 коэффициенты при 1, cos и sin приходим к трем уравнениям для определения четырех коэффициентов: A0 , A1 , С0 , С1 : A0 C0 0 R 2 kR C0 C1 A0 1 1 ln kR 0 2 U 2 R 2 k 2 kR C0 A1 C1 1 U ln kR 0 2 2 R 2 2 2 kR Решим эту систему уравнений: A0 C U 0 2 R 2 C0 1 C1 A1 C1 1 1 ln k kR 2kR2 R 2 U 2 ln 2 kR 2kR2 0 30
2U C0 1 C 1 C ln kR 0 ln kR 0 0 2 2 2 2 2 1 C0 1 2 ln kR 4U , 2 C0 A0 4U , 1 1 2 ln kR 2 2U 4 e , 1 1 2 ln 1 kR k 1 2 ln kR 2 2 A1 C1 2k U R2 1 1 2 ln kR 2 При принятой степени приближения коэффициенты A1 и C1 по отдельности определены быть не могут. Подставляя найденные значения коэффициентов в формулы (3.26), находим выражения для проекции скоростей, пригодные в области вблизи цилиндра: Vr U cos R3 r 1 2 ln 2 r R 1 1 2 ln kR 2 V U sin R r 1 2 2 ln r R 1 1 2 ln kR 2 (3.27) 2 На больших расстояниях от цилиндра, беря только члены, содержащие A0 и C 0 , и рассматривая абсолютное движение цилиндра в газе, покоящейся на бесконечности, будем иметь: A0 1 C0 e kr cos K 0 (kr) cosK 0 (kr), r 2 1 V С0 e kr cos K 0 (kr) sin . 2 Vr (3.28) Приведем теперь формулу для давления; вследствие (3.5) и (3.6) имеем: 31
Pe P eU A0 cos cos 2 eU A1 ... r r2 и, ограничиваясь необходимым приближением, получаем Pe P eU A0 cos r (3.29) Вычислим силу, действующую на цилиндр; общая формула [4] F prr cos pr sin ds S в данном случае упрощается; 2 F R prr cos pr sin d , (3.30) 0 где prr и pr компоненты тензора вязких напряжений [4]. Подставляя (3.30) полученные выше выражения и после интегрирования, получаем F 2 e A0 (3.31) Подставляя найденное выше значение A0 , получаем формулу Ламба для величины сопротивления, испытываемого цилиндром при его движении в газе, причем эта сила, отнесенная к единице длины цилиндра равна F 8 eU 1 1 2 ln kR 2 (3.32) Таким образом, используя метод Ламба, нами получены выражения для полей скорости и давления в окрестности цилиндрической частицы, и они имеют следующий вид: 32
A0 C0 C C 1 A C U 0 ln kr 21 0 1 2 cos Vr r 2kr 2 2 2 2kr r C A C 1 V 21 0 ln kr 1 2 U sin (3.33) 2 2 2kr r cos Pe P eU A0 r где A0 , C0 , A1 , C1 постоянные интегрирования, k U e / 2 e , ln 0.5772 . 33
Глава IV. Вывод выражения для силы и скорости термофореза. В главе 3 используя метод Ламба, нами были получены выражения для полей скорости и давления в окрестности цилиндрической частицы, и они имеют следующий вид: A0 C0 C C 1 A C U 0 ln kr 21 0 1 2 cos Vr r 2kr 2 2 2 2kr r C A C 1 V 21 0 ln kr 1 2 U sin (4.1) 2 2 2kr r cos Pe P eU A0 r где A0 , C0 , A1 , C1 постоянные интегрирования, k U e / 2 e , ln 0.5772 . Учитывая (4.1) можно найти полную силу, действующую на аэрозольную умеренно крупную частицу цилиндрической формы. Полная сила, как мы отмечали в главе 1, определяется интегрированием тензора вязких напряжений по поверхности частицы и имеет следующий вид [4] 2 F R prr cos pr sin d , (4.2) 0 где prr , pr компоненты тензора вязких напряжений и в цилиндрической системе координат равны: prr Pe 2e V 1 Vr V Vr , pr e . r r r r Вычисляя компоненты тензора вязких напряжений с использованием (4.1) и после интегрирования, получаем F 2 e A0 (4.3) Постоянная интегрирования A0 определяется из граничных условий на поверхности твердой крупной аэрозольной частицы цилиндрической формы и равна 34
e K TS e i U T e T e 1 i A0 2 (4.4) 1 1 1 ln kR 2 2 Подставляя A0 в (4.3) получаем выражение для полной силы, действующей на единицу длины крупной аэрозольной цилиндрической частицы e 4 e K TS e i F U T T e 1 e 1 1 1 ln kR i 2 2 Из этой формулы видим, что полученная сила, действующая на единицу цилиндра, аддитивно складывается из термофоретической силы Fth и силы сопротивления среды FC F Fth Fc , (4.5) где e K TS e2 i Fth 4 T (4.6) T e 1 1 e 1 ln kR 1 2 2 i Fc 4 e U 1 1 1 ln kR 2 2 , (4.7) Как мы выше отмечали, поскольку частица движется равномерно, то это означает, что полная сила, действующая на нее равна нулю и скорость термофореза равна с обратным знаком величине скорости U . Приравнивая выражение (4.5) к нулю получаем выражение для скорости термофореза умеренно крупной аэрозольной частицы цилиндрической формы 35
U th e i K TS e T e 1 e i T (4.8) 36
Глава 5. Анализ полученных результатов Используя метод Ламба для гидродинамической части и уравнения Лапласа для теплопереноса (уравнение теплопроводности) в дипломной работе получены выражения для термофоретической силы и скорости твердой умеренно крупной аэрозольной частицы цилиндрической формы, которые имеют следующий вид: e K TS e2 i Fth 4 T , T e 1 1 e 1 ln kR 1 2 2 i U th e i K TS e T e 1 e i T (5.1) (5.2) При получении аналитических выражений для силы и скорости термофореза мы учли, что внутри крупной аэрозольной частицы цилиндрической формы действуют неоднородно распределенные источники тепла. Сравним скорости термофореза твердой крупной цилиндрической и сферической частиц. Скорость термофореза сферической частицы равна[12]: U сф 2 e i K TS e T e 1 2 e i T (5.3) При сопоставлении формул (5.2) и (5.3) видно, что при предельном переходе скорость термофореза крупной цилиндрической частицы будет переходить в скорость термофореза крупной сферической частицы. Коэффициент перехода будет равен: 37
U th K TS e T T e e 1 2 i 1 1 e 1 2 e i K TS e i T T e 1 e i 2 U сф e i В частном случае, когда i e и 0 , скорость термофореза сфериR ческой частицы почти в два раза больше скорости частицы цилиндрической формы. Скорости термофореза для сферической и цилиндрической частиц почти совпадают, если e i . 38
Заключение Рассмотрен термофорез крупной твердой аэрозольной частицы цилиндрической формы, внутри которой действуют неоднородно распределенные по ее объему тепловые источники. Получены аналитические выражения для полей температуры вне и внутри аэрозольной частицы, а также выражения для силы и скорости термофореза цилиндрической частицы. 39
Список использованной литературы 1. Страус В. Промышленная очистка газов. М.: Химия. 1981. 616 с. 2. Вальдберг А.Ю., Исянов Л.М., Яламов Ю.И. Теоретические основы охраны атмосферного воздуха от загрязнения промышленными аэрозолями. Санкт-Петербург: НИИОГАЗ-ФИЛЬТР. 1993. 235 с. 3. Хаппель Дж,, Бреннер Г. Гидродинамика при малых числах Рейнольдса М.: Мир. 1976. 630 с. 4. Ландау Л.Д., Лифшиц Е.М. Теоретическая физика. Т. 6. Гидродинамика. М.: Наука. 1986. 736 с. 5. Борен К., Хафмен Д. Поглощение и рассеяние света малыми частицами М.: Мир. 1986. 660 с. 6. Волковицкий О.А., Седунов Ю.С., Семенов Л.П. Распространение интенсивного лазерного излучения в облаках. Ленинград: Гидрометеоиздат. 1982. 312 с. 7. Береснев С.А., Кочнева Л.Б. Фактор асимметрии поглощения излучения и фотофорез аэрозолей// Физика атмосферы и океана. 2003. Т. 16. № 2. С.134-141 с. 8. Поддоскин А.Б., Юшканов А.А., Яламов Ю.И. Теория термофореза умеренно крупных аэрозольных частиц //ЖТФ. 1982. Т. 52. Вып. 11. С. 22532261 9. С. Чепмен, Т. Каулинг Математическая теория неоднородных газов М.: Иностр. литература. 1960. 510 с. 10. А.Н. Тихонов, А.А. Смарский Уравнения математической физики М.: Наука. 1972. 735 с. 11.Камке Э. Справочник по обыкновенным дифференциальным уравнениям. М.: Физ.-мат. лит-ра. 1961.703 с. 12.Ю.И. Яламов, В.С. Галоян Динамика капель в неоднородных вязких средах Ереван: Луйс. 1985. 89 с.
13. Н.Е. Кочин, И.А. Кибель, Н.В. Розе Теоретическая гидромеханика. М.: Наука. 1963. Ч. II. С. 727 с. 14. Г. Ламб Гидродинамика. М.: ОГИЗ. 1947. 928 с. 41
Отзывы:
Авторизуйтесь, чтобы оставить отзыв