Дипломная работа: Исследование процессов испарения и конденсации жидких капель

и  - отношение молекулярных масс A и B. Первый член в знаменателе (8.9) отвечает за столкновения между молекулами А, в то время как второй - за столкновения между молекулами A и молекулами B. Если концентрация молекул A очень низка (в действительности, в атмосферных приложениях, так и есть),  и (8.9) может быть упрощена, фактическим пренебрежением столкновениями между молекулами А:

                                       (2.10)

Отметим, что молекулярная концентрация  может легко быть вычислена из уравнения Клайперона , где p - давление системы. Средний{скупой} путь свободного пробега газа проекции прямой на заднем плане газ не зависит от концентрации непосредственно. Это не удивительно, поскольку мы предположили, что концентрация A настолько низка, что молекулы никогда не сталкиваются и не взаимодействуют друг с другом. Далее мы сосредотачиваемся на взаимодействиях частиц с единственным газом - воздухом, где средняя длина свободного пробега задана формулой (2.6).

2.1.2 Непрерывный режим.

Неустановившаяся диффузия молекул вида A к поверхности частицы радиуса :

                             (2.11)

где c(r,t) концентрация молекул А, а - молярный поток (количество молей падающих на единицу площади в единицу времени) в любом радиальном положении r. Это уравнение - просто выражение массового баланса в бесконечно малой сферической ячейке, вокруг частицы. Молярный поток молекул A дается согласно закону Фика (Бирд и др., 1960),

                     (2.12)

где  - молекулярная фракция частиц A,  - радиальный поток воздуха в положении r, и  - коэффициент диффузии молекул А в воздухе. Так как выделенных направлений нет,  во всех r. Принимая во внимание предположение, применимое в большинстве атмосферных состояний - , (2.12) может быть переписана как

                                        (2.13)

Теперь, комбинируя (2.11) и (2.13), получим:

                                      (2.14)

Если  - концентрация молекул А вдали от частицы, а - концентрация паровой фазы на поверхности частицы, и частица первоначально находится в атмосфере частиц А с концентрацией, равной , начальные, и граничные условия для (11.4) записываются так:

                                   (2.15)

                                             (2.16)

                                             (2.17)

решение (2.14) в граничных условиях (2.5) - (2.17), будет выглядеть так:

    (2.18)

Временная зависимость концентрации в любом радиальном положении r дается третьим членом на правой стороне (2.18). Отметим, что для больших значений t, значение верхнего предела интегрирования приближается к нулю и профиль концентрации приближается к установившемуся состоянию, задаваемому

                                (2.19)

Полный поток молекул А (молей в секунду) к частице обозначен Jc, индекс c показывает, что режим непрерывный (continuum), и задаётся, как:

                                                (2.20)

или, используя (2.19) и (2.13), как

                                   (2.21)

Если , поток молекул A - к частице, а если  - наоборот. Вышеупомянутый результат был впервые получен Максвеллом (1877), и (11.11) часто называется потоком Максвела.

Массовый баланс на растущей или испаряющейся частице:

                                        (2.22)

где  - плотность частицы и  молекулярный вес A. Объединяя (2.21) с (2.22) получим,

                                   (2.23)

Когда  и  постоянны, (2.23) можно проинтегрировать, что даст:

                                 (2.24)

Использование независимого от времени установившегося профиля, заданного (2.19), для вычисления размера частицы во времени (11.24) может казаться противоречивым. Использование установившегося диффузионного потока, для вычисления темпа роста частицы подразумевает, что профиль концентрации пара около частицы достигает установившейся величины прежде, чем произойдёт заметное изменение величины молекулы. Так как рост действительно происходит в сотни раз медленнее чем диффузия, профиль около частицы фактически всегда остается в ее стационарном значении.

2.1.3 Свободно - молекулярный (кинетический) режим.

В трёхмерном случае число столкновений молекул с единицей поверхности в единицу времени равно  (Moore, 1962)

,                                           (2.25)

где - скорость молекул:

   .                                     (2.26)

Учитывая это, молярный поток  (молей в единицу времени) на частицу радиусом :

                                 (2.27)

где  - вероятность прилипания. Отношение молекулярного потока в кинетическом режиме к потоку в непрерывном режиме, равно:

                                        (2.28)

2.1.4 Переходный режим.

Установившийся поток молекул пара к сфере, когда частица является достаточно большой по сравнению со средней длинной свободного пробега молекул пара, задаётся уравнением Максвелла (2.20). Так как это уравнение основано на решении уравнения переноса в непрерывном режиме, оно перестаёт действовать, когда средняя длина свободного пробега молекул пара становится сопоставимым диаметру частицы. В другом случае, выражение, основанное на кинетической теории газов (2.27) также не справедливо в этом случае, где . Когда , явления, как говорят, лежат в переходном режиме.

Распределение концентрации диффузионных молекул и фонового газа в переходном режиме строго описывается уравнением Больцмана. К сожалению, не существует общего решения уравнения Больцмана, справедливого для всего диапазона чисел Кнудсена. Как следствие, большинство исследований явлений переноса избегает решать непосредственно уравнение Больцмана и ограничивают себя подходом, основанным на так называемом методе подгонки потоков. Подгонка потоков предполагает, что кинетические эффекты ограничены областью , а вне этой области имеет место непрерывный режим.

 Расстояние  имеет порядок средней длины свободного пробега . Предполагают, что в пределах этой внутренней области применима простая кинетическая теория газов.

• Теория Фукса Соответствие непрерывных и свободномолекулярных потоков молекулы относится ко времени Николая Альбертовича Фукса, который предложил, что подгонкой двух потоков в , можно получить граничное условие к уравнению диффузии. Предположим, коэффициент прилипания равен единице,

   (2.29)

Тогда решая стационарное уравнение переноса для разбавленной системы,

                                   (2.30)

используя как граничные условия (11.27) и , получаем решение:

                            (2.31)

где поправочный коэффициент :

                     (2.32)

Связав бинарную диффузию и среднюю длину свободного пробега, используя , и , получим:

              (2.33)

Заметим, что определение средней длины свободного пробега  подразумевает, что для a= 1,

                                               (2.34)

и отношение Фукса (2.33) преобразуется, используя (2.34),

                      (2.35)

Значение , используемого в выражениях выше не было определено в теории и должно быть выбрано опытным путем или оценено в соответствии с независимой теорией. Несколько выборов для  были предложены: самое простое, самим Фуксом, =0. Другие предложения по этой теме высказаны Дэвисом , в 1983 году: , .

• Подход Фукса и Сутугина. Фукс и Сутугин в 1971 году последовали решению уравнения Больцмана, данного Сахни в 1966 году, для , где - отношение молекулярного веса диффундирующего вещества и воздуха, для создания следующей интерполяционной формулы переходного режима.

                                     (2.36)

Уравнение (2.36) основано на результатах для  и поэтому непосредственно применимо, чтобы описать молекулы в более тяжелом фоновом газе. Средняя длина свободного пробега, включенная в определение числа Кнудсена в (2.36): .

• Подход Дахнеке. Дахнеке(Dahneke) в 1983 использовал поток, соответствующий подходу Фукса, но, предполагая что -  и , получил,

                                     (2.37)

где . Средняя длина свободного пробега свободного пробега, включенного в определение числа Кнудсена в (2.37):

• Подход Лоялки. Лоялка в 1983 построил улучшенные интерполяционные формулы для переходного режима, решая линеаризованное уравнения Больцмана с помощью БГК модели (Bhatnagar, Gross, Krook - 1954):

                (2.38)

Средняя длина свободного пробега, используемая Лоялкой: , коэффициент скачка концентрации имел значение . Виллиамс и Лоялка в 1991 году указали, что (2.38) не работает вблизи свободно - молекулярного режима.

2.2    Подведение итогов

Для получения возможно более точных результатов по испарению и конденсации частиц применяются самые разные подходы: от полуэмпирических, некоторые из которых перечислены выше, до достаточно обоснованных с математической точки зрения. К этим работам, в частности, относится серия работ [16] - [21]. В основе этих работ лежит расчет потока пара на частицу интегрированием функции распределения, полученной в результате решения линеаризованного уравнения Больцмана:

       (2.39)

здесь:  - одночастичная функция распределения по скоростям и координатам i-ого газового компонента,  - телесный угол,  - вектор относительной скорости,  - сечение столкновений,  - скорость молекулы, которая сталкивается с рассматриваемой молекулой - для которой записывается уравнение Больцмана, суммирование производится по всем газовым составляющим. Вообще говоря, в левой части уравнения (2.39) следует добавить слагаемое - , где  сила, действующая на молекулу, m - ее масса,  - ее ускорение. Предполагается, что силовое поле отсутствует. В таком виде уравнение Больцмана слишком сложное, чтобы для него можно было найти решение, кроме самых простых случаев, например, для равновесного распределения по скоростям. В таком виде оно используется для исследования ее решений. Правая часть этого уравнения называется интегралом столкновения, вся сложность поиска решений связана именно с этим интегралом столкновений. В частности существует принцип Гильберта [26], [27], в соответствии с которым решение уравнения (2.39) можно найти в виде разложения по моментам распределения в начальный момент времени. На этом основан метод моментов Греда. Однако этот метод более применим к задачам гидродинамики, нежели, к проблемам кинетики. Основные приближения, которые используются для получения решения уравнения (2.39) сводятся к тому, чтобы упростить интеграл столкновений. При этом предполагается, что распределение по скоростям мало отличается от равновесного распределения. Таким образом конструируется уравнение для функции, описывающей отклонение распределения от равновесного. Этот подход аналогичен методам, описанным в работе Черчиньяни [21], [22]. В конечном счёте этот метод приводит к интегральному уравнению Фредгольма первого или второго рода - в зависимости от выбранной формы аппроксимации. Дополнительные осложнения возникают при постановке граничных условий. Наибольшие продвижения возможны в этом направлении при сферической форме испаряющихся капель. Попытки получить точное решение приводят к довольно сложным зависимостям, с которыми сложно работать и сопоставлять с экспериментальными данными. Кроме этого, приходится делать предположение скачка концентраций на поверхности частицы. Для диффузионного и около диффузионного режима столкновений молекул пара с частицей, когда задачу можно свести к решению уравнения диффузии, авторам [17] удалось создать метод расчета конденсации и испарения для несферических частиц, используя формализм функций Грина - задача сводится к решению соответствующего интегрального уравнения, при этом могут быть использованы численные методы - аналитические зависимости в этом случае получить не удается. Еще сложнее описать процессы испарения и конденсации частиц, в среде, состоящей из нескольких летучих компонентов [23]. Предполагалось, что процесс стационарный, испаряющиеся компоненты химически инертны, пары представляют собой идеальный газ. Для переходного режима использовалась формула Фукса - Сутугина. По сути, этот подход представлял собой применение ранее разработанных моделей для бинарной смеси. Сопоставление модельных расчетов с экспериментальными результатами испарения смеси азотной кислоты с водой показало, что при различных внешних условиях (соотношениях компонент и относительной влажности) большинство моделей можно применять, если подогнать соответствующим образом модельные параметры, например, вероятность прилипания.

Страницы: 1, 2, 3, 4, 5, 6