Исследование процессов испарения и конденсации жидких капель (стр. 7 из 11)

(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]. Предполагалось, что процесс стационарный, испаряющиеся компоненты химически инертны, пары представляют собой идеальный газ. Для переходного режима использовалась формула Фукса - Сутугина. По сути, этот подход представлял собой применение ранее разработанных моделей для бинарной смеси. Сопоставление модельных расчетов с экспериментальными результатами испарения смеси азотной кислоты с водой показало, что при различных внешних условиях (соотношениях компонент и относительной влажности) большинство моделей можно применять, если подогнать соответствующим образом модельные параметры, например, вероятность прилипания.

С точки зрения корректности постановки задачи и нахождения ее решения следует отметить работу Сахни [24]. Вообще говоря, эта работа относится к расчету потока нейтронов - в замедлителе реактора с учетом поглощения их черными сферами. Тем не менее, она полностью применима для испарения и конденсации частиц сферической формы. По сути дела эта работа - решение задачи Милна для сферической геометрии. Эта задача была сведена к интегральному уравнению первого рода, затем преобразовано к сингулярному уравнению типа Коши, которое было затем решено численными методами. Эти результаты, в частности, были использованы Фуксом и Сутугиным [8] для получения аппроксимационной формулы переходного режима. Очевидно, что для получения точных соотношений для потока молекул на поверхность частицы, необходимо решать уравнение для функции распределения Больцмана. Все определяется тем, в каком виде брать правую часть этого уравнения, ниже мы вернемся к этому.

2.2 Постановка задачи.

Наука об аэрозольных частицах началась cрешения проблемы испарения в газообразной среде.

Предположим, нам известен радиус частицы – a, и концентрация молекул пара вокруг этой частицы – n(r), гдеr– радиус-вектор, построенный из точки начала координат (центр частицы). Требуется найти зависимость радиуса частицы от времени – a(t). Итак, массу частицы можно выразить как

, где ρ – плотность частицы, площадь поверхности частицы - , тогда будет справедливо уравнение:

, (2.39)