Математическая модель эволюции аэрозолей при ультразвуковом воздействии Текст научной статьи по специальности «Физика»

УДК 544.733.422:519.87

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ЭВОЛЮЦИИ АЭРОЗОЛЕЙ ПРИ УЛЬТРАЗВУКОВОМ ВОЗДЕЙСТВИИ

О.Б. Кудряшова, А.А. Антонникова

Институт проблем химико-энергетических технологий СО РАН, г. Бийск E-mail: olgakudr@inbox.ru

Математическая модель основана на уравнении Смолуховского, описывающем динамику изменения функции распределения частиц аэрозолей по размерам сучетом ультразвукового воздействия, испарения (для жидкокапельных аэрозолей) и осаждения. Проведено исследование асимптотического поведения функции вероятности столкновений частиц аэрозоля от частоты ультразвука.

Ключевые слова:

Аэрозоль, распределение по размерам, ультразвуковое воздействие, испарение, коагуляция, осаждение.

Key words:

Aerosol, size distribution, uitrasonic influence, evaporation, coagulation, sedimentation.

Процессу ультразвуковой (УЗ) коагуляции аэрозолей посвящено много работ. Многие авторы рассматривают процессы взаимодействия частиц теоретически, но результаты эксперимента, тем не менее, описываются эмпирическими формулами.

В работах [1, 2] приводятся математические модели, использующие балансовый подход Смолу-ховского в сочетании с выкладками ортокинети-ческой и гидродинамической гипотез взаимодействия частиц. При этом проведено сравнение с экспериментом по оптимальным частотам и амплитудам звукового воздействия. Однако полученные выражения не наглядны, не позволяют просмотреть асимптотические закономерности, и для прогнозирования процессов коагуляции аэрозолей приходится либо проводить численный эксперимент, либо пользоваться более наглядными, но не обоснованными физически интерполяционными формулами, описывающими физический эксперимент.

В работе предложены выражения для кинетики коагуляции в зависимости от основных параметров воздействия, свойств аэрозоля и среды: частоты и амплитуды звуковых колебаний, концентрации и дисперсного состава исходного аэрозоля, вязкости и температуры среды, физико-химических параметров материала частиц.

Рассмотрим трансформацию распределения частиц в произвольном облаке по размерам с течением времени ¡. Следуя [3-5], запишем балансовое уравнение (интегральный вариант уравнения Смо-луховского), описывающее изменение со временем функции распределения частиц по размерам: д/ (Б, г)

dt

/- +12 +1 з,

(1)

где К(Б,Б1) - вероятность столкновений частиц;

Н - верхняя граница облака; g -

Аф _

Рчї

ускорение свободного падения; п - динамическая вязкость среды; рч - плотность частицы. Все частицы, масса которых превышает критическое значение Дф(0, выпадают из облака и не принимают дальнейшего участия в коагуляции; спектр частиц на каждый момент времени / будет обрезан справа за счет седиментации крупных частиц, причем, постепенно эта граница будет смещаться в сторону все более малых частиц.

Член 12 описывает возникновение частиц диаметра Б за счет столкновения капель с диаметрами Д и Б-Д:

/2 _

1 а

- JK(D -Dj,Dj)f (Dj,t)f (D -Dj,t)dDj.

Жидкие капли испаряются тем быстрее, чем меньше их размер, за счет кривизны поверхности. Член 13 описывает уменьшение массы частиц за счет их испарения с учетом уравнения Кельвина для капиллярного эффекта:

'■-£ (£/ (°)

_д_

dm

2nDf M(pdrop - ppl )f (Dy

RT

где 11 описывает убыль частиц с диаметром Б за единицу времени в единице объема за счет столкновения капли диаметра Б с любой каплей диаметра Д:

Аф(' )

I, = -/(А г) | к (А Д)/(Д, г )оА,

где т - масса капли; Б/ - коэффициент диффузии; М - молекулярная масса жидкой капли; Я - универсальная газовая постоянная; Т - абсолютная температура; р1щ и рр! - парциальное давление над каплей и плоской поверхностью.

Член уравнения 13 имеет смысл для жидкокапельных аэрозолей, состоящих из микронных и суб-микронных капель, испарение для которых существенно из-за кривизны их поверхности, причем при невысокой влажности окружающей среды (р1щ>>^рр)

Начальные условия для уравнения (1): при = /(Б,/0)=/0(Б) - начальное распределение частиц

U

по размерам. Для описания функции распределения частиц по размерам обычно применяют гамма-распределение: f(D)=aD*exp (—bD), где b, а- параметры распределения; a - нормировочный коэффициент. В качестве характерного размера данного распределения можно выбрать модальный диаметр D0=a/b.

Относительная концентрация аэрозоля уменьшается по сравнению с начальной п0 за счет испарения и седиментации, и на момент времени t это уменьшение составит:

Д t (Аф(<)/Do ® ^

— =J J I3dx + J f (D, t )dx

Пг\ n n

D„(t)/ D„

dt. (2)

Важным вопросом является определение вероятности столкновений частиц; чем она выше, тем быстрее произойдет коагуляция и осаждение аэрозоля. В отсутствие какого-либо воздействия эта величина обуславливается броуновским движением [3]. В работе [4] вероятность столкновений частиц аэрозоля без дополнительного воздействия считалась пропорциональной сумме масс частиц:

K (D, Д) = ^ D3 + Д3),

(3)

где кь - коэффициент пропорциональности; V -кинематический коэффициент вязкости среды.

В условиях ультразвукового воздействия повышается вероятность столкновений частиц. В выражение (3) должны войти параметры, характеризующие ультразвуковое воздействие, прежде всего амплитуда и частота излучения.

Согласно [6], коагуляцию можно рассматривать как образование агрегатов в звуковом поле в результате взаимодействия частиц, вызванного акустическими течениями вокруг них в звуковом поле. Следуя этой работе, рассчитаем число встреч частиц в звуковом поле, полагая, что частицы встретятся в том случае, если линия тока потока, возникающего около одной из частиц и увлекающего вторую, пройдет в «трубке» диаметром, равным двум диаметрам частицы 2Б. Число встреч на единицу длины выделенной частицы за единицу времени равно и0, где - площадка, огра-

ниченная линиями тока; и0 - скорость движения частиц. При невысоких уровнях звукового давления в [6] получено, что число встреч пропорционально диаметру частиц, квадрату скорости их движения (которая, в свою очередь, определяется амплитудой звуковых колебаний), концентрации частиц, коэффициенту обтекания кобт, и обратно пропорционально вязкости среды:

U0 По к026т (D3 + D3)

(4)

Взвешенная в газе частица под действием сил звукового поля вовлекается в колебательное движение («ортокинетическая» гипотеза). В зависимости от свойств среды, размеров и плотности частицы она может увлекаться средой лучше или хуже,

что определяется коэффициентом увлечения крл -отношение амплитуды скорости взвешенной частицы к амплитуде скорости частицы газа. Считая, что между частицей и средой действует сила Стокса, получим формулу для коэффициента увлечения в виде [7]:

1

у™1 _ Г, 2 2 ’ (5)

VI + о т

где о - частота акустического воздействия; т=рчБ2/18п - время релаксации частицы; п - динамический коэффициент вязкости среды. Анализ этого выражения показывает, что амплитуда колебания частицы тем больше отличается от амплитуды колебания окружающей среды, чем больше размер и плотность частицы, чем выше частота звука и меньше вязкость.

Учитывая выражение для коэффициента увлечения частицы в звуковом поле (4), можно считать, что частица увлекается звуковым полем, повышая вероятность столкновений с фиксированной частицей в к2обтка(1—кувл)2 раз, где ка - коэффициент пропорциональности. Тогда вероятность столкновений частиц с диаметрами Б и Б1 пропорциональна Ж(1-ку1л)2, или, с учетом (2)-(5):

к (д А) =

knL( D3 + Dj3)

1 -

1

ч2\

Vi + а2т2

(6)

Анализируя выражение (6), заметим, что:

• при отсутствии акустического поля (каио2=0) вероятность столкновений сводится к броуновской (3);

• с повышением амплитуды воздействия (а значит, скорости и0) вероятность повышается;

• акустическое воздействие с относительно низкими частотами (ю2!2«!) неэффективно

(крл—1 (!—крл)—0);

• существует оптимальная частота воздействия (ю2!2»!), наиболее эффективная с точки зрения коагуляции аэрозолей (крл—0, (!—кувл)—>!). Такое поведение функции (6) соответствует экспериментальным и теоретическим данным [1, 2, 4]. Пользуясь выражением (6), найдем предельно низкие (неэффективные) и оптимальные частоты ультразвукового воздействия:

1

® ■ = —

mm

т

1/

1-

Ку1 kbn0D03 ) UoJk

-1

.0(1)

— 1

(7)

(8)

Обозначим повышение вероятности столкновений частиц при ультразвуковом воздействии через

( / Л2 Л

K{. K =

1+KUо2

1 -

1

-\А + а2т2

. Поведение по-

VJ =

max

V

вышающего коагуляцию коэффициента в зависимости от частоты иллюстрирует рис. 1 (водный аэрозоль, каЩ= 1).

К1

Рис. 1. Зависимость коэффициента /С, от частоты УЗ-воздей-ствия для водного аэрозоля: 1) D=5 мкм; 2) D=1 мкм

Расчет для водного аэрозоля минимальной и оптимальной частот воздействия по формулам (7) и (8) приведен на рис. 2.

(а), кГц

Рис. 2. Оптимальная (1) и минимальная частота (2) для водного аэрозоля

Рис. 3 иллюстрирует уменьшение массы субми-кронного аэрозоля во времени в процессе коагуляции и испарения под действием ультразвука и без УЗ-воздействия. Параметры начального распределения частиц аэрозоля по размерам: Ь=1, а=0,1. Неиспаряющийся аэрозоль (кривая 1) коагулирует и осаждается медленно; ультразвуковое воздействие существенно ускоряет этот процесс. ((ют)2=2, кривая 2). Быстро испаряющийся субми-кронный водный аэрозоль теряет массу в результате испарения (кривая 3); ультразвуковое воздействие еще увеличивает скорость убыли массы аэрозоля (кривая 4). Эти расчеты соответствуют многочисленным экспериментальным данным, в том числе полученным авторами ранее [4]. В опытах, результаты которых приведены в работе [4], время полного осаждения аэрозоля (мука, вода) уменьшалось при УЗ-воздействии в 1,5-2,0 раза.

Начальный спектр распределения частиц по размерам претерпевает деформацию в течение времени. На рис. 4 приведено массовое распределение частиц по размерам для водного аэрозоля в начальный момент времени (кривая 1) и через 1 минуту при УЗ-воздействии для неиспаряемого (кривая 3) и испаряемого водного аэрозоля (кривая 2). Массовая

функция распределения по размерам g(Б) связана со счетной соотношением g(Б)=m/m10/(Б), где т10 -среднеарифметическая масса частиц:

да

Що = |щ/(Б)^’ т - масса частицы диаметра Б. о

т/т0

О 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 100 , с

Рис. 3. Зависимость относительной массы аэрозоля от времени: 1) без испарения и без УЗ-воздействия; 2) без испарения, с УЗ-воздействием; 3) с испарением (вода), без УЗ-воздействия; 4) с испарением и с УЗ-воз-действием

г(Я)

0,5 1 1,5 2 2,5 3 3,5 4 4,5 5 5,5 6 6,5 7 7,5 8 8,5 9 9,5 В, мкм

Рис. 4. Массовое распределение частиц по размерам при УЗ-воздействии: 1) начальный момент времени; 2) и 3) неиспаряемый аэрозоль через 1мин

Спектр распределения неиспаряемого аэрозоля за этот промежуток времени практически не изменился; водяные капли из-за быстрого испарения сильно уменьшились в размерах, сместив максимум функции распределения в область ~ 10-7 микрона. Качественный вид кривой 2, рис. 4, подтвержден в экспериментах по динамическим измерениям функции распределения субмикронных частиц водного аэрозоля [8]; количественное расхождение теоретических и экспериментальных результатов по средним диаметрам и концентрации частиц не превысило 20 %, что говорит об адекватности предлагаемой модели.

Выводы

На основе модели коагуляции аэрозоля под действием ультразвука в виде варианта интегрального уравнения Смолуховского, со стоком (испарение) и обрезанием спектра (осаждение), получе-

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

СПИСОК ЛИТЕРАТУРЫ

1. Хмелев В.Н., Шалунова К.В., Цыганок С.Н., Барсуков Р.В., Сливин А.Н., Шалунов А.В. Ультразвуковая коагуляция аэрозолей. - Бийск: АлтГТУ, 2010. - 228 с.

2. Хмелев В.Н., Шалунов А.В., Голых Р.Н., Шалунова К.В. Теоретическое исследование процесса акустической коагуляции газодисперсных систем // Измерения, автоматизация и моделирование в промышленности и научных исследованиях (ИАМП-2010): Матер. VII Всеросс. научно-техн. конф. -Бийск, 6-7 октября 2010 г. - Бийск: Изд-во Алт. гос. техн. ун-та, 2010. - С. 222-227.

3. Волощук В.М. Кинетическая теория коагуляции. - Л.: Гидро-метеоиздат, 1984. - 284 с.

4. Антонникова А.А., Кудряшова О.Б., Хмелев М.В., Шалунов А.В. Исследование эффективности применения ультразвуковых колебаний для осаждения мелкодисперсных аэрозолей // Научно-технический вестник Поволжья. - 2011. - №6. -С. 96-99.

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

5. Кудряшова О.Б. Математическая модель эволюции жидкокапельных аэрозолей // Известия Томского политехнического университета. - 2012. - Т. 320. - № 2. - С. 129-133.

6. Розенберг Л.Д. Физические основы ультразвуковой технологии. - М.: Наука, 1970. - 689 с.

7. Денисов А.С., Подольский А.А., Турубаров В.И. Об увлечении аэрозольных частиц в звуковом поле при числах Рейнольдса <1 // Акустический журнал. - 1965. - Т. 11. - Вып. 1. -С. 146-155.

8. Титов С.С. Турбидиметрический высокоселективный метод и быстродействующий измерительный комплекс определения параметров нестационарных многофазных сред: дис. ... канд. техн. наук. - Бийск, 2011. - 153 с.

Поступила 10.04.2012 г.