CC BY
114
УДК 544.733.422:519.87
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ВЗРЫВНОЙ ГЕНЕРАЦИИ ЖИДКОКАПЕЛЬНЫХ АЭРОЗОЛЕЙ
О.Б. Кудряшова, Б.И. Ворожцов
Институт проблем химико-энергетических технологий СО РАН, г. Бийск E-mail: olgakudr@inbox.ru
Математическая модель содержит уравнения, описывающие динамику изменения термодинамических параметров в корпусе устройства взрывного генератора жидкокапельных аэрозолей, процесс истечения и размер распыленных капель. Описаны процессы генезиса аэрозольного облака и его дальнейшей эволюции. Найден критерий, характеризующий эффективность процессов кавитации. Получены оценки времени распространения волны сжатия и времени истечения потока; определены зависимости размера кавитационных пузырьков и частиц аэрозоля от параметров конструкции генератора и характеристик веществ; построено распределение частиц аэрозоля по размерам; установлены закономерности изменения функции распределения частиц по размерам с учетом процессов коагуляции, испарения, осаждения.
Ключевые слова:
Жидкокапельный аэрозоль, кавитация, распределение по размерам, испарение, коагуляция. Key words:
Liquid-drop aerosol, cavitation, size distribution, evaporation, coagulation.
При проведении миротворческих операций в XXI в. оружие нелетального действия или просто «гуманное» оружие, действующее в атмосфере, будет основным [1]. Один из видов такого оружия -останавливающие аэрозоли, например, с использованием вытяжки красного жгучего перца. Физико-химическое оружие нелетального действия позволяет с помощью аэрозольных боеприпасов распылять в районе расположения боевой техники противника химические вещества, приводящие к порче и остановке двигателей и электрогенераторов. К настоящему времени уже исследован широкий перечень химических ингибиторов горения. Аналогичные исследования проведены для составов, разлагающих резину.
В таких случаях требуется быстрое распыление высокодисперсных аэрозолей в некоторой, иногда труднодоступной, зоне. При этом, чем выше дисперсность получаемых аэрозолей, тем лучше достигаемый эффект, поскольку высокая удельно-массовая поверхность капельного объема дает большую суммарную площадь испарения мелких капель, что повышает скорость воздействия химических агентов. Таким образом, большой практический интерес представляют аэрозоли с характерным размером частиц порядка одного микрона и менее, причем в ряде задач требуется быстрое (мгновенное) получение таких аэрозолей без изменения физико-химических свойств диспергируемых веществ. Эта проблема до сих пор не была решена.
Именно взрывной способ позволяет достичь высокой скорости получения аэрозолей. Взрывной способ распыления давно используется на практике и хорошо изучен [2, 3]. Однако до сих пор не ставилась задача и, соответственно, не были получены таким способом аэрозоли субмикронных размеров, несмотря на потребности описанных выше практических приложений. Для получения микронных и субмикронных аэрозолей в работе предложено специальное устройство взрывного распылителя.
Настоящая работа посвящена теоретическому и экспериментальному исследованию процесса генерации жидкокапельных аэрозолей с помощью новой конструкции взрывного распылителя (рис. 1), а также дальнейшей эволюции таких аэрозолей. Генератор представляет собой цилиндрический корпус - 1, содержащий контейнер с водой - 3, ограниченный мембранами - 4. На одном из днищ корпуса имеется заряд взрывчатого вещества - 2, другой конец имеет зазор (сопло) - 6, ограниченный отражателем - 5, через который происходит выброс воды после инициирования заряда. При этом создаются условия для кавитации в воде, как будет показано ниже.
Рис. 1. Схема взрывного генератора аэрозоля
Максимальное давление в волне рт определяется давлением газов, образующихся во взрывной камере объемом VI. В предположении мгновенной детонации давление газов составит:
V-11 <
Рт =Г
У
V
где 0 - энергия взрывчатого превращения; у- показатель адиабаты продуктов детонации. Давление пе-
редается через мембрану в столбик жидкости, порождая ударную волну с амплитудой рт, которая возбуждает собственные колебания столбика жидкости. Пусть 2т - амплитуда смещения частиц в возбужденной волне; р„ - плотность жидкости; р0 -плотность воздуха; с - скорость распространения волны в жидкости, равная скорости звука. Тогда
Z =
т
2Q
tpxSxG) С
отсюда D1 =
Р = Pm
1 +1-
Pm
4lL\k pw ((Si / S2)2-1)
пли диаметром, равным толщине водяного слоя, можно оценить количество таких капель. Диаметр кавитационного элемента до истечения из отверстия распылителя обозначим Бь после истечения, в момент разрушения, диаметр элемента - Д, пузырька под атмосферным давлением р0 - Б3. Расширение будем считать происходящим мгновенно
где частота колебаний а>=(пс)/Ь. Время достижения волны противоположной поверхности жидкости равно 4=Ь/с, где Ь - высота столбика жидкости. Тогда, с учетом выражения для массы жидкости М„=Ь8рю где - площадь сечения столбика жидкости, получим безразмерный параметр:
(процесс адиабатический): — =
( \ Р
Ро
Ъу
Учиты-
вая условие равенства объема воды в элементе до и после истечения (Д3-Д3=Д3-Д3), получим выражение для толщины водяного слоя в момент разрушения к:
h = D
L ПС
характеризующий эффективность процессов кавитации - величина плоского разрыва по отношению к высоте столбика жидкости.
Полученное значение максимального смещения частиц жидкости Zm в волне разгрузки приводит к появлению разрыва и испарению жидкости в эту пустоту. Следующая затем фаза сжатия вызывает образование шарообразного пузырька в массе воды диаметром D1, заполненной водяными парами. Плоская волна разрежения порождает пл-оский разрыв с эффективной толщиной: Zef/=Zm/V2. Выделим в этом разрыве элемент с характерным размером Zf. Таких элементов образуется: N=Sx]Z]j!. Элемент будет частично заполнен жидкой водой, частично - парами воды, его плотность составляет около 0,8pw [3]. Это соответствует индексу кавитации k=VyV&0,8, где Ve - объем элемента, Vw -объем жидкости в кавитационном элементе. Объем элемента V^Z^. Тогда Д3=(1-к) Zf.
Процесс расширения кавитационного пузырька адиабатический, поэтому pm((1-k)Zef^r=pDir,
L(1 - k)
\Mw
смесь, представляющая собой водяной каркас с распределенными в нем кавитационными пузырьками, под действием давления за счет дальнейшего расширения продуктов детонации вырывается через отверстие. Каждая единица объема жидкости содержит включения, разрушающие неразрывность потока, размер которых D1 зависит от давления.
Зависимость давления в жидкости определим, решая совместно уравнение неразрывности и уравнение Бернулли. Получим:
- К
1 - k
f p V/y
КГ о У
, \ 1/Ъу
Р
КГ о У
Пароводяная
Если предположить, что в момент истечения кавитированной жидкости каждый пузырек, находящийся под давлением р, раздуется до атмосферного, и в этот момент лопнет, разорвавшись на ка-
Исходный слой I парогазовой смеси площадью и толщиной к разобьется, в конечном итоге, на п
, 65,
аэрозольных частиц диаметром к, ц =—Ц при' пЬ
чем с падением давления и увеличением диаметра частиц с каждым моментом времени их число будет падать. Введем относительное количество ча-
N
стиц: пгеН = п / ^ ц, где п - количество частиц
/=1
слоя ц N - количество слоев жидкости (вплоть до полного истечения). Считая послойно количество выбрасываемых частиц, с учетом постепенного уменьшения давления, можно построить зависимость относительного количества частиц от их размеров, а затем аппроксимировать полученное выражение функцией гамма-распределения.
Проведем численные оценки и расчеты в соответствии с предложенной моделью взрывного генезиса жидкокапельного аэрозоля.
Выберем следующие параметры распылителя: ^=1 см3, толщина слоя жидкости Ь=1 см, площадь гидродинамической трубки ^=1 см2. Пусть вещество обладает следующими характеристиками: 0=50 кДж - энергия взрыва, у=1,3 - показатель адиабаты выделившихся газов. Тогда максимальное давление жидкости составитрт=11,5 МПа, амплитуда смещения 2п^671 мкм. Безразмерный параметр, характеризующий эффективность процессов кавитации в нашем случае имеет величину около 0,067 - максимальный разрыв жидкости, обусловленный ударной волной, составит около 7 % от высоты столбика жидкости. Это достаточно большое значение, которое свидетельствует о том, что процессы кавитации в рассматриваемой задаче существенны.
Следуя изложенному выше алгоритму, получим для наших условий зависимость относительной массы капель от их диаметра, которую можно аппроксимировать с помощью гамма-функции распределения частиц по размерам. Параметры такого распределения составят: а=0,1, Ь=0,9,
m
Рис. 2. Массовая функция распределения частиц по размерам
(Дх)=ох"ехр(-Ьх^); массовая функция распределения по размерам связана со счетной соотношением: g(x)=m/m10Дx), где т10 - среднеарифметиче-
да
ская масса частиц: т10 = |т[(хх); т - масса ча-
0
стицы диаметра х). Параметры функции распределения частиц по размерам, полученные экспериментально, через 1 с после распыления хорошо согласуются с расчетными: «=0,11, ¿=0,6. Модальный диаметр счетного распределения частиц по размерам составляет около 1 мкм; 80 % частиц имеют диаметр менее 20 мкм. Массовая функция распределения частиц по размерам приведена на рис. 2.
Итак, наибольшая масса полученных взрывным кавитирующим способом капель имеет диаметр менее 10 мкм, и для них существенны процессы испарения, обусловленные кривизной их поверхности.
Динамика испарения малых капель определяется уравнением Максвелла:
1 = т = АягР/М(рЛтр - рр1)
Л ЯТ '
где Б} - коэффициент диффузии; М - молекулярный вес жидкой капли; Я - универсальная газовая постоянная; Т- абсолютная температура; г - радиус частицы (в нашем случае г=Н/2); - парциальное давление над каплей; рг1 - парциальное давление над плоской поверхностью. Обратимся к формуле Томсона (Кельвина), определяющей давление насыщенного пара над каплями жидкости:
,, , , 2аМ
1п(РЛгор / РР,) =——>
р р —ЛТг
где а - коэффициент поверхностного натяжения вещества капли.
Капля полностью испарится в момент времени ¡, при котором г=0. Таким образом, время жизни капли определится выражением:
,.2
ь =
г
2МБгрр
ЯТ—
( ( ехр
2аМ
—ЯТГо
-1
0 У
Для водяного аэрозоля: а=0,0727 Н/м, М=18 г/моль, 7=293 К, рр=2486 Па. Учитывая выражение для времени жизни капли и полученные ранее параметры функции гамма-распределения, выясним, какая массовая доля диспергированного аэрозоля испарится в первые секунды. Так, за 1 с испарятся все капли с диаметром меньшим, чем 1,54 мкм (массовая доля т/т0 таких частиц составит 42 %), за 2 с - 1,94 мкм (55 %) и так далее. Результаты полученного расчета (кривая) и эксперимента (точки) динамики отношения конечной массы капель т к начальной массе т0 приведены на рис. 3. Уже через 15 с останется только 10 % от исходной массы аэрозоля, остальные 90 % массы жидкости испарятся. Это хорошо согласуется с данными эксперимента.
Рассмотрим эволюцию распределения частиц по размерам с течением времени. Следуя [4], запишем балансовое уравнение (интегральный вариант уравнения Смолуховского), описывающее изменение со временем функции распределения частиц по размерам:
дg (т, t)
Ы
■ = I + Ь +1,
Здесь 11 описывает убыль капель с массой т за единицу времени в единице объема за счет столкновения капли массой т с любой каплей массы т':
да
11 = -g(т, t) |К(т, т')g(т', {)йт',
где K(m,m) - вероятность столкновения капель с массами m и m' в единицу времени. Определим вероятность столкновения капель (ядро уравнения Смолуховского) как линейную функцию их масс: K(m,m )=b (m+m).
Член I2 описывает возникновение частиц массой m за счет столкновения капель с массой m и m-m ':
1m '
I2 - — J K (m - mm')g (m , t)g (m - mt)dm
член 13 - уменьшение массы капель за счет их испарения, определяется уравнением Максвелла, отнесенным к массе частицы:
7 ■ =дт (^ (т \
д 4лгИ/ М(р^ - рр1)g(т)
дт ЯТ
На рис. 4 приведены результаты численного расчета массовой функции распределения через
10 и 45 с после распыления в соответствии с приведенной выше моделью.
Массовое распределение постепенно «размазывается» за счет коагуляции, и его максимум смещается в сторону меньших размеров за счет испарения мелких капель.
Выводы
Представлены результаты теоретического исследования взрывного способа генерации водяного аэрозоля, который можно рассматривать как предельный случай гидравлического метода распыле-
ния: резкий скачок давления в системе происходит в результате действия взрыва. Показана роль кавитации в процессе формирования высокодисперсного аэрозоля. Предложена модель дальнейшей эволюции полученного аэрозоля с учетом процессов испарения и коагуляции; с помощью численных расчетов получено распределение частиц аэрозоля по размерам в зависимости от времени. Представленные результаты сравнения экспериментальных и теоретических исследований свидетельствуют о физической адекватности предлагаемой математической модели.
СПИСОК ЛИТЕРАТУРЫ
1 Буренок В.М., Корчак В.Ю., Смирнов С.С. Оружие нелетального действия - составная часть системы вооружения будущего // Вестник Академии военных наук. - 2007. - № 4. -С. 117-127.
2 Стебновский С.В. Импульсное диспергирование как предельный режим разрушения жидкого объема // Физика горения и взрыва. - 2008. - Т. 44. - № 2. - С. 117-128.
3 Кедринский В.К. Газодинамика взрыва: эксперимент и модели. - Новосибирск: Изд-во СО РАН, 2000. - 435 с.
4 Волощук В.М. Кинетическая теория коагуляции. - Л.: Гидро-метеоиздат, 1984. - 284 с.
Поступила 17.06.2010 г.
УДК 533.9.01
ОСОБЕННОСТИ ЗАТУХАНИЯ ЭЛЕКТРОМАГНИТНОГО ПОЛЯ В ПЛАЗМЕ ВЫСОКОЧАСТОТНОГО ЕМКОСТНОГО РАЗРЯДА
В.А. Власов, Ю.Ю. Луценко, Н.В. Корепанова, Е.П. Зеленецкая
Томский политехнический университет E-mail: luts@mail.ru
Приведены результаты измерений осевого распределения гармонических составляющих электрического поля емкостного разряда, горящего в среде аргона и воздуха, при изменении частоты основной гармоники. Установлен полирезонансный характер затухания гармоник в плазме разряда. Рассмотрено влияние температуры электронов, и соответственно концентрации электронов плазмы разряда на затухание частотных составляющих. Установлен аномальный рост третьей гармонической составляющей электрического поля при уменьшении электронной температуры плазмы разряда.
Ключевые слова:
Высокочастотный емкостной разряд, гармоники, электромагнитное поле, электронная температура, плазма, параметрический
резонанс.
Key words:
High-frequency capacitive discharge, harmonic components, electromagnetic field, electron temperature, plasma, parametric resonance.
Введение
Высокочастотный разряд емкостного типа, горящий при атмосферном давлении, представляет собой тонкий плазменный шнур, окружённый слабосветящейся диффузионной оболочкой. Горение высокочастотных разрядов емкостного типа осуществляется за счёт [1, 2] диссипации поверхностной электромагнитной волны, распространяющейся вдоль канала разряда. По условиям возбуждения поверхностной электромагнитной волны в плазме разряда данные разряды можно подразделить на высокочастотный факельный разряд и высокочастотный емкостной разряд. Высокочастотный факельный разряд зажигается с конца металлического электро-
да, обычно имеющего цилиндрическую форму. Возбуждение же емкостного разряда осуществляется электродом, охватывающим разрядную трубку. Поэтому высокочастотный емкостной разряд в отличие от высокочастотного факельного разряда является безэлектродным разрядом.
В работе [3] нами было установлено отсутствие затухания первых двух гармонических составляющих электромагнитного поля вдоль оси факельного разряда, горящего в воздухе. Для случая факельного разряда, горящего в среде аргона, было установлено отсутствие затухания первых трёх гармоник. В то же время с увеличением частоты электромагнитной волны её затухание должно возрастать. Для
