CC BY
25
УДК 551.578.71:551.509.616
О ФРАКЦИОННОЙ МОДЕЛИ НАКОПЛЕНИЯ ЗАГРЯЗНЯЮЩИХ ВЕЩЕСТВ В ТУМАНЕ
Харченко А. В., Бояршинов М. Г., Трусов П. В. (Пермь)
Abstract
The fractional model descnbing processes of accumulation of impurities in drops of a fog for arbitrary speeds of condensation, absorption of impurities, nucléation etc. is constructed. The model admits possibility of irregidar splitting on a fraction, that allows essentiaUy to reduce their amoimt.
Введение
Накопление газообразных загрязняющих веществ в жидких аэрозольных частицах может оказывать сильное влияние на характер распространения примеси в атмосфере и осаждения на подстилающую поверхность. Построение модели, учитывающей процессы взаимодействия аэрозолей и газообразных примесей, представляет собой актуальную задачу при исследовании переноса и рассеивания атмосферной примеси.
Моделирование процессов, происходящих в полидисперсной среде, строится на основе решения кинетического уравнения для функции распределения частиц по размерам [2, 3]. Для изучения накопления примеси в каплях аэрозоля в [3] получено аналитическое решение кинетического уравнения в предположении, что скорость изменения массы воды в каплях находится из формулы Максвелла для квазистатического испарения сферических капель, справедливой только для крупных капель, а скорость накопления примеси — по равновесному закону Генри. При этом в [3] пренебрегается нуклеацией, коагуляцией и дроблением капель. Обобщение полученного аналитического решения на случай произвольных скоростей конденсации и растворения примеси, а также учет взаимодействия капель вызывает затруднения. Использование метода моментов для решения кинетического уравнения даже в случае отсутствия примеси, но с учетом зависимости скорости конденсации от размера капель, приводит к появлению заметной погрешности [2]. Поэтому представляется актуальным построение модели, способной описать накопление газообразной примеси в каплях аэрозоля в случае сильной зависимости скорости массообмена капель от их размеров, которая может быть легко обобщена на случай нескольких примесей.
Хорошо зарекомендовали себя при изучении процессов в полидисперсных средах фракционные модели (ФМ), когда весь спектр капель разбивается на конечное число фракций. При этом выделяют два подхода к построению ФМ — эйлеров и лагранжев [2]. С одноименными подходами в механике сплошных сред эти названия связаны только по принципу описания эволюции размеров капель. Эйлеров подход к построению ФМ предполагает постоянные размеры капель одной фракции (функция распределения капель по размерам аппроксимируется кусочно-постоянной функцией с неизменными интервалами однородности), при изменении размеров частица переходит в другую
фракцию независимо от причины, вызвавшей это изменение. При лагранжевом описании размеры капель внутри одной фракции могут меняться за счет фазовых переходов на поверхности частиц, коагуляции или. дробления (функция распределения капель аппроксимируется суммой конечного числа дельта-функций, которые могут смещаться по шкале размеров). Лагранжев подход позволяет более естественно описывать непрерывное изменение размеров капель, однако затрудняет учет нуклеации, в то время как эйлеров подход дает возможность простого учета зародышеобразования, но размер капель при этом меняется дискретно. В ФМ используется гипотеза о взаимопроникающих континуумах, что позволяет ввести для дисперсной фазы распределенную плотность и считать, что в каждой точке иолидисперсная среда имеет несколько скоростей, температур и т.д. ФМ используются для изучения движения газа с частицами конденсата в соплах [2], динамики облаков и туманов.
В представляемой работе предлагается эйлерова ФМ накопления гаюобразных примесей в каплях водного аэрозоля. Модель построена для произвольного разбиения капель на фракции, что позволяет учитывать коагуляцию капель без потери массы крупных капель даже при сравнительно небольшом количестве фракций. Одновременный учет накопления нескольких примесей позволяет исследовать химические превращения загрязняющих веществ в каплях туманов. В предложенной ФМ отсутствует ограничение малости массы примеси в капле, что расширяет рамки применения модели на высот«• концентрированные примеси.
Цостроение фракционной модели
Введем функцию распределен!¡и капель f\t,x. u‘''j по массе содержащихся в них
воды д'0-’ и загрязняющих веществ |Ди\ i-\, S Б каждый момент времени t В произвольной точке рассматриваемой области х Кинетическое уравнение баланса в фазовом пространстве, аналогичное уравнению Лиувнлдя, можно записать в виде
-4......»‘"'Ьо
сх •-о Щ 1=1
Здесь второе слагаемое в левой части описывает перенос капель всех размеров со скоростью которая предполагается одинаковой для капель одной массы, третье
слагаемое описывает увеличение или уменьшение капель за счет фазовых переходов на поверхности ИЛИ сорбцию ИЛИ десорбцию Примеси ЖИДКОСТЬЮ, При ЭТОМ j —
скорость изменения массы воды в капле; ¡к'"{х, — изменения массы /-той примеси в
капле; первое слагаемое в правой частя описывает юмогашую конденсацию, J(x) —
(0)
скорость нуклеации капель, )я, ' — масса критического зародыша, произведение дельтафункций выражает предположение, что концентрация примесей в газовой смеси много меньше насыщающей, т е. вновь образовавшиеся капли полностью состоят из воды, функционал СГ описывает взаимодействие капель —• исчезновение мелких капель и увеличение количества крупных капель при коагуляции, обратные процессы при
дроблении и т.п, функция определяет химические превращения веществ внутри
капель. Аргумент х указывает на зависимость скоростей процессов от полей температур (несушей среды и капель), концентраций (пара и примеси) и т. п. Принимается гипотеза об однородном распределении температуры и концентрации примеси по всему объему капель, при этом считается, что капли одной массы имеют равную температуру.
сС ос
Величина и: = [соответствует количественной концентрации
о о
капель (количество капель в единице объема) всех размеров, а величина
СС СО
р,, -]•••[ — распределенную плотность всей дисперсной фазы, при этом
о о
5
т ~ Ё ' — полная масса капли
Аппроксимируем функцию /(?, кусочно-постоянной функцией по
и) /-
аргу ментам и следующим ооразом:
.V«1» Л'>*> 5 Г, V./- 4,-1
/0,-о /л,=о *=о
где Цд ’ = 0 (ц<0> вводится для упрощения записи, сохраняет смысл массы
критического зародыша, / 0) ¡5, = 0), А(-) — ступенчатая функция Хевисайда. Будем
‘ 0_/' '
считать, что капли, удовлетворяющие неравенствам |л(*(,/) < принадлежат к
одной фракции. Температура и скорость всех капель, принадлежащих к одной фракции, предполагаются одинаковыми, равными соответственно И и Т0) ^.Функционал
взаимодействия капель и функция химического взаимодействия веществ предполагаются также независимыми от массы в пределах каждой фракции, а скорости изменения массы
воды и примесей /уо; в пределах одной фракции завися! от |ди ' линейно.
Далее для краткости аргументы функций опускаются. Количественная концентрация капель каждой фракции определяется выражением
10)
Л’
I |.АЛ1'-...ф-' /
а распределенная плотность соотношением Р
", • (Л)
..О! ,.<*'> ' " 1=0
(3)
"г" (о) (*Г ' ' "У-г- .2/,=0
■л0) .Ц;;'.Г:
— 11^ — И*'1-1. П1 , ™ —(и-). , 4-15-1 , ) Ж ... — ^ м^.
Используя символ Кронекера
5„ =
V
(?)
и интегрируя но массе воды и примесей в капле в пределах одной фракции уравнение (1), умноженное на массу капли, после использования аппроксимации (2) можно получить уравнение для распределенной плотности каждой фракции
и $ ^
Т.'Р. о ~ ^ ¡?)Р I о; , - Л,,., Пб д,) т<5 Ю! ..¡х) + % .ы- ,'Я ’|"
р/ - ■■■/ (;х ./ Г '-I ,= 1 'Ь ! ■■■.: : •
М--:)б
Р.% ,о- 5
(6)
+
+
-и
Здесь при д - -1 сумма в правой части описывает увеличение числа капель в рассматриваемой фракції» за счет конденсации водяного пара и поглощения примеси каплями более малого размера, <] = 1 - испарения жидкосчи и выделения примеси каплями большего размера, при у = 0 сумма описывает уменьшение количества капель за счет перехода в другие фракпиг., последнее слагаемое обусловлено изменением
распределенной плотности за счет масеообмена капель без перехода из одной фракции в др> тую
Определении) вида функций J, !" , % и функционала а посвящено множество работ, среди которых можно упомянуть классические теории: гомогенной нуклеации Френкеля- Зельдовича [1]; квазистатического испарения Максвелла [5]; коагуляции Смолуховекого [4] и т. д. Не следует забывать, что величины J, а (о: ,5,, /\’ и
у .. ,. являются функциями температуры капель и газовой смеси, концентрации пара и
;ш рязняющего вещества, поэтому дія полноты модель должна включать уравнения дкижекия и энергии всей полидисперсной среды и каждой отдельной фракции, аналогичные приведенным в [2]. а также уравнения диффузии пара и примеси с соответствующими источниковыми членами, которые догъкны описывать массообмен с
каплями.
Решение полученной системы дифференциальных уравнений в частных
производных вида (6), дополненной упомянутыми уравнениями, в обшем виде
представляет трудности с вычислительной точки зрения лаже при небольших значениях -V ". Поэтому необходимы дополнительные упрощающие предположения. В качестве
последних может выступать гипотеза об односкоростной среде [2] или более сильная гипотеза о неподвижной и однородной среде [3]. Считая концентрацию капель невысокой, можно пренебречь функционалом взаимодействия [2].
Коагуляцию, происходящую под действием броуновских сил, можно описать с помошью модели Смолуховского [4], в которой интенсивность взаимодействия капель зависит только от их размеров и диффузии, состав капель при этом может быть произвольным. Для капель, состоящих из нескольких веществ, можно записать
ц!«> ц!5)
(7)
00 зо
о о
^/ = 0
5
V
І-
1 = 9
где К — коэффициент коагуляции, который определяется как вероятность столкновения капель различных размеров с последующим слипанием. Первый интеграл в правой части описывает увеличение числа крупных капель за счет слияния более мелких, а второй — исчезновение капель за счет коагуляции с каплями всех размеров.
Для простоты предположим, что
д/;о>
VI я
I - І І - £ [V-
/(5>=о*ІОІ=о я- 51=о ' ..... (8)
< п<(>> - фОО#« - - № Ш. -
то есть для любых двух фракций коэффициент коагуляции остаегся постоянным в пределах каждой из фракций.
Теперь, умножая (7) на массу капли, используя аппроксимации (2) и (8), интегрируя по массе в пределах каждой фракция, можно получить
,{¡1 /(01
СТ,о.
Р>(С|.../'5'1Рд:,0| <<*>*{<>)
=0Я"' ,=0
П" *!) ¡/"-^‘4)4^4) ДЛ^ ДИ^
4 II
тшГ I
ІГ
1 + 5„
И' и, +К0
' {■;) (і) ]/ (і)і+8й ОН+л., і
И/// +^,)Д^ #| +
(9)
к': /
,(>)1+8й
.И+з.
'Ил
2 + 5,
{¡}2-Ь.г. (/)2+5„
■ и
())2+5 !/)2+8
/1
'г1///
+ и Г
’Р
у*. 2, ••• Л ---------------
/о:-0 /5)=о "Ь'О^М и*І,
п\-" V.) (,'
и)
где ц/ = тіпі тах^ц/(і),ц.(„ пу,: ,КУ-г:)+і) ; ц# - тах: ц;(\'„т
= тіп[ тахіні1;'/, , тіщ Ц(,) +1 + ^'
ц(,у = шіп( тах( и1,',' \ ч' ; , та:4 Цу, .и”!
Громоздкий вид первой суммы в правой части обусловлен произвольностью выбора Ац|:,’ , что позволяет сократить количество фракций, чго, в свою очередь, приводит к
повышению эффективности расчетов.
Расчет простейшего случая
Для иллюстрации применения разработанной модели рассмотрим простейший модельный случай накопления одного загрязняющего вещества в водном тумане без химических реакций. Среда предполагается неподвижной и однородной, температурными эффектами пренебрегаем. Скорость конденсации пара определяется формулой Максвелла
!{(" -- 4ттуг>(и,г(с^)‘ - с °‘). 0°)
где и" - - коэффициент диффузии пара, г - ЦЪт / 4трж — радиус капли, форма которой предполагается сферической, р.*. -• плотность жидкости капли, с'с< — концентрация пара. с'г ' — равновесная концентрация пара при данной температуре. Тогда для /'¡01 ;и можно записать выражение
\
ЛЩ ' г' г ; Л,(0) {0;\
/до , 4яр^ I [сх с. ). (11)
Скорость накопления примеси как и в [1], определяется равновесным законом Генри
1! ! 1 1' '
>•"’ д . (12)
V 4 7Т7* Н
Iі' - 4л:Л)аУ
где 1У" — коэффициент диффузии примеси, с " — концентрация примеси; Н безразмерная константа Генри. Отсюда можно получить выражения для Л1,,:
Ът
/‘ji-.ro --4тс/;(1)
V
4~Р :
■ І і \
К Р. !
' т:^НУ (к,)
і 3 /
Коэффициент коагуляции в неподвижной среде для сферических частиц
определяется соотношением [4]
. . Ч 2кТ(\ 1'. .
Кип т) -■-----: — -г 7 л/-?/■;. гм,
V V?' /V - ' (Т4)
где У .г - радиусы капель с массами т, т, к постоянная Больцмана, Т.\ —
температура і; вязкость среды.
Далее будем считать, что происходит только конденсация и поглощение примеси /'|С) ,5, > 0, процессом нуклеации пренебрегаем (,/ = 0)
Как указывалось выше, полученную систему дифференциальных уравнений необходимо дополнить уравнениями диффузии пара и примесей, которые с учетом высказанных предположений об однородности и неподвижности среды могут быть записаны в виде
СІС
■м
„.(0)
(ІІ
1т,
(15)
При расчетах принималось, что вязкость среды и коэффициенты диффузии не зависят от состава и температуры. В начальный момент все капли имеют размер критического зародыша, при этом пх о = 2 • 105 см ', с'1'1 =5-1СГ6 г/смг, с1" = 1СГ9 г/см', Т - 237.15 К. Подобные условия могут реализовываться при образовании радиационного тумана, который возникает при более быстром ночном охлаждении поверхности земли за счет лучеиспускания по сравнению с малоподвижным воздухом над ней [1]. Рассматривалось 400 фракций (.М1-01 = 40, Аг{: =10). Полученная система обыкновенных дифференциальных уравнений решалась методом Рунге-Кутта четвертого порядка аппроксимации.
у
а) б) в)
Рис. 1. Зависимость концентрации пара и примеси от времени (а) и распределенная плотность дисперсной фазы, полученная без учета коагуляции (б)
и с учетом коагуляции (б)
На рис. 1(6,в) представлена распределенная плотность фракций капель в зависимости от содержания в каплях воды и примеси в момент установления равновесной концентрации примеси и пара На графике видно, что первоначально практически монодисперсная примесь с течением времени стала полидисперсной, и при этом содержание примеси в каплях одинаковой массы может быть различно. На рис Ка) приведена зависимость концентрации пара с'"" и примеси с'0 от времени. Для установления равновесной концентрации при постоянной температуре требуется время порядка половины секунды.
Выводы
Построена эйлерова фракционная модель, описывающая процессы накопления примеси в каплях туманов для произвольных скоростей конденсации, поглощения примеси, нуклеации и т.и Модель допускает возможность неравномерного разбиения на фракции, что позволяет существенно сокращать их количество.
В отличие от решения, полученного в [3], в представленной модели не требуется априорного знания функции распределения капель по полной массе: напротив, это распределение получается в результате решения. В отличие от [3] отсутствует также ограничение на массу примеси в капле.
Библиографический список
1 Амелин Л.Г Теоретические основы образования гумана при конденсации пара. - VI.
Химия, 1972. - 304 с.
2. Болдарев A.C., Гасило» В. А , Зайчик Л И., Ольховская О.Г. Численное моделирование
квазиодномерных и двумерных течений спонтанно конденсирующегося пара в трансзвуковых соплах Н Теплофизика высоких температур, - 1998. - Т. 36. - №1. - С. 135 - ¡40.
3. Степанов A.C.. Захарова И.М., Новикова Л.Д. Моделирование процессов накопления загрязняющих веществ в каплях туманов // Метеорология и гидрология. - 1997 - №4. -С. 25-36.
4 Туницкий H.H., Каминский В.Л., Тимашев С.Ф, Методы физико-химической кинетики - М.: Химия, 1972. - >98 с
5 Фукс Н А Испарение и рост капель в газообразной среде. - М.: Изд-во АН СССР, ! °5S - 91 с.
Получено 15.06,99.
