CC BY
26
УДК 621.375.826
К.В. Соломатин
Влияние параметров расширения насыщенного пара на спектр размеров частиц конденсата
Расширение насыщенного пара имеет место в ряде практических технологических задач и может сопровождаться конденсацией в потоке (разлет продуктов испарения при лазерной абляции, движение газа в сопле и т.д.) [1]. Далее рассмотрим именно процессы конденсации, не останавливаясь подробно на газодинамических аспектах в описании явления.
Спектр размеров дисперсной конденсированной фазы представляет несомненный научный интерес, поскольку, с одной стороны, влияет на параметры течения потока, а с другой - определяет вид и состав получаемого в процессе аэрозоля. Последнее означает теоретическую возможность контроля величины среднего размера вторичных частиц и их дисперсии Путем варьирования динамики параметров состояния системы.
При адиабатическом расширении двухфазной системы закон сохранения энергии можно представить в виде (в приближениях равновесного расширения и однородного потока) [2]:
(с;и-х)+<?гх) сіт+яга - х)^г = их,
У
(1)
Г" Г‘
гдеЧ и (г -молярнаятеплоемкость паровой и конденсированной фаз соответственно; % -степень конденсации (доля молекул, относящихся к конденсированной фазе, в общем числе молекул); Т- абсолютная температура потока; Я- универсальная газовая постоянная; V-объем системы; I - скрытая теплота фазового перехода.
Величина степени конденсации может быть представлена в виде интеграла
* = Н>-/(г)Жг)£(/,г >/г,
(2)
где /(/) - скорость зародышеобразования (число частиц конденсата, образующихся за единицу времени, в расчете на одну молекулу паровой фазы); #(/,т) - число молекул во вторичной частице, образовавшейся в момент времени т к моменту времени Л
Равновесная скорость роста отдельной вторичной частицы может быть описана уравнением [3]
Л
(
■ 4л-д,^2
‘-у-ехр
і
ЯГ
_Ь_
ЯГ
(3)
//
Р„ = Р* -ехр -
ЯГ
(4)
где рг - постоянная, то получаем Г, = X/ Я1п
На начальном этапе расширения возможны два механизма возникновения зародышей: гомогенный (в случае чистого пара имеющего низкую температуру) и гетерогенный (в остальных случаях), который обычно и реализуется на практике.
Для случая гетерогенной конденсации в химически однородном паре вокруг каждого иона возникает вторичная частица и далее, вплоть до наступления сильного разрежения, переохлаждение (9=1- Т/Т ) намного меньше единицы. При этом уравнение (3) принимает вид
' А
л 2 2/1 — /-'У 4 1 2/1 —
— ~4тга,.е /IV-------------------= 4 яа.,г гIV
сії ЯГ
ЯГ
- 1п
Р,
<5)
здесь аа~ размер молекулы пара; и и у- - концентрация и средняя скорость молекул газа; Т — температура газа, насыщенного при данной плотности.
Поскольку давление насыщенных паров зависит от температуры по закону
а для степени конденсации Х(0 = ХЛ‘)’
где / - начальная степень конденсации, равная степени ионизации вещества.
Для равновесного состояния эту величину несложно рассчитать, используя формулу Саха [4].
Для решения задачи сделаем несколько упрощений и приближений.
Во-первых, перейдем от временной переменной к объему, поскольку скорость его изменения (и), согласно постановке задачи, считаем заданной. Тогда с1\' = Бис/!, где 5 - сечение цилиндрического сосуда (для одномерного случая расширения).
Во-вторых, будем считать пар идеальным газом, подчиняющимся уравнению /?Г=у(1~х)ЯГ (где V -количество вещества), для того чтобы исключить давление.
В-третьих, пренебрежем различием между Су и Су и примем С, = /?/('/— 1), где у - показатель адиабаты.
В-четвертых, обезразмерим переменные с помощью формул: Т -Т /Тп, V = V / Уа, где индекс «О»
(здесь и всюду далее) означает начальное значение величин.
В результате получаем систему из двух дифференциальных уравнений первого порядка
Влияние параметров расширения насыщенного пара.
, . сії Т 1 <іх
Х^У~ЕТ0с1У'
*£ = А-ХУ\ 1-^) —• с{у У
тгг.
О -Х)Т
(7)
(8)
V / \1 У
где А = —“-Д, "—безразмерный параметр (/0-начальная 372 (я)(|
длина цилиндра; (Л)и - начальная длина свободного пробега в газе).
Начальные условия:
Т = 1, при V = 1.
(9)
Численное решение системы выполним с помощью программы ООРШ8 (модифицированный метод Рунге-Кутты [5]).
Результаты решения представлены на рисунках 1-3.
V/ К
Рис. 1. Зависимость степени конденсации от относительного объема двухфазной смеси
Для численного решения системы (7-8) выберем пар углерода (поскольку давление пара становится заметным при температурах более 3500 К, что соответствует значительной температурной ионизации). Диапазон начальных температур (3,5 4-4,5)-103 К, для него степень ионизации составляет х0= (1 3) • 10-6
и практически линейно зависит от температуры. Для упрощения расчетов и их интерпретации скорость расширения выберем постоянной в диапазоне (1 + 100) м/с.
Согласно результатам численного решения (рис. 1), величина степени конденсации в случае расширения с постоянной скоростью практически не зависит от скорости расширения. Температура (рис. 2) медленно
Рис. 2. Зависимость относительной температуры в системе от относительного объема двухфазной смеси
Рис. 3. Зависимость среднего размера вторичной частицы от начальной температуры насыщенного пара
понижается по линии, примерно соответствующей кривой упругости пара. Оба результата, по-видимому, объясняются допущением равновесности процесса: переохлаждение в системе очень незначительное. При увеличении объема примерно в два раза степень конденсации находится в пределах (0,04-0,05). Поэтому и зависимость от начального состояния системы (температуры или давления) вплоть до момента нарушения теплового равновесия является слабой.
С учетом начального числа центров конденсации зависимость среднего размера (диаметра а) частиц конденсата от температуры при двукратном увеличении объема представлена на рисунке 3. Для температуры 4000 К средний размер составляет около 8 нм. Как видно, для получения более крупных частиц требуется выбирать меньшее значение начальной температуры насыщенного пара (однако при этом массовая концентрация частиц будет существенно снижаться).
Библиографический список
1. Горбунов. В.Н. Неравновесная конденсация в высокоскоростных потоках газа / В.Н. Горбунов, У.Г. Пирумов. -М., 1984.
2. Райзер, Ю.П. О конденсации в облаке испаренного вещества, расширяющегося в пустоту / Ю.П. Райзер // ЖЭТФ. - 1959. - Т.37, вып. 6(12).
3. Хирс, Д. Испарение и конденсация / Д. Хирс, Г. Паунд. - М., 1966.
4. Физический энциклопедический словарь / Гл. ред. А.М. Прохоров. - М., 1995.
5. Хайрер, Э. Решение обыкновенных дифференциальных уравнений. Нежесткие задачи : пер. с англ. / Э. Хайрер, С. Нерсетт, Г. Ваннер. - М., 1990.
