CC BY
36
УДК 551.578.71:551.509.616
А. В. Харченко Пермский государственный технический университет
ОПРЕДЕЛЕНИЕ ОБЛАСТИ ОСАЖДЕНИЯ АЭРОЗОЛЯ С ПОМОЩЬЮ ФРАКЦИОННОЙ МОДЕЛИ
Abstract
In the present work the fractional model for research of the formed toxic precipitation from cm gaseous or dispersible impurity and water atmospheric aerosol is offered. The model is illustrated by the solution of a model problem on one-dimensional propagation of a impurity aerosols, which enters interaction with an atmospheric aerosol with formation, at the end\ of toxic precipitation.
Введение
Проблема определения зоны выпадения токсичных осадков вокруг промышленных зон в настоящее время весьма актуальна. Многие токсичные вещества выбрасываются в атмосферу промышленными объектами в дисперсном или в газообразном состоянии, после чего они могут взаимодействовать с каплями туманов и облаков, образуя токсичный аэрозоль, который переносится воздушными массами и может выпадать в виде осадков.
Накопление примеси в каплях аэрозоля, перенос капель в атмосфере и их выпадение сопровождаются различными процессами, которые можно условно разделить на четыре основные группы (рис. 1). Однако, если ограничиться невысокими скоростями ветра, положительной температурой воздуха, а также отсутствием химических реакций в каплях, то количество процессов, рассмотрение которых необходимо для моделирования образования и выпадения токсичных осадков, можно существенно уменьшить.
Рис. 1. Процессы происходящие с атмосферным аэрозолем
В настоящее время исследованию процесса накопления и переноса примеси атмосферным аэрозолем уделяется много внимания. Однако известные в настоящее время
работы, посвященные данному вопросу, концентрируются на моделировании одного или нескольких процессов и поэтому не обладают необходимой полнотой. Так, в [1] наряду с учетом поверхностных фазовых переходов игнорируется процесс коагуляции, являющийся, по существу, определяющим фактором роста капель. В работах [2, 3, 4], напротив, рассматривается только процесс коагуляции. Исследование [5] наиболее полно охватывает процессы массообмена, однако, как и во всех упомянутых ранее источниках, исследуются только моменты распределения капель по массе составляющих их веществ, чем достигается высокая эффективность методов, при этом теряются детали распределения, которые важны при исследовании выпадения токсичных осадков. Кроме того, во всех указанных работах рассматривается пространственно однородный случай В [6] исследуются перенос и диффузия токсичного аэрозоля с учетом коагуляции, однако не учитывается наличие в атмосфере чистого водного аэрозоля, не рассматривается также и массообмен на поверхности аэрозольных частиц.
Таким образом, существует необходимость построения модели, описывающей образование, перенос, диффузию и выпадение токсичных осадков с учетом коагуляции, нуклеации, поверхностных фазовых переходов и гравитационной седиментации.
Постановка задачи
Для простоты рассматриваются однородные по объему капли, которые находятся
в газовой смеси, состоящей из воздуха, паров воды (с концентрацией с^) и Я загряз-
(') —
няющих веществ (с концентрацией с ,/ = 1,5). Вводится функция распределения капель /(?,х.по массе содержащихся в них воды /г''- и загрязняющих веществ /и ' , / = 1,5 в каждый момент времени / в произвольной точке рассматриваемой области х Кинетическое уравнение баланса в фазовом пространстве записывается в виде
■у 5 /у ^ ^ ^ у'
^. + у.и/ + £-----^- = V■DVf+J + a-G. (1.1)
^ 1=о дц
Здесь второе слагаемое в левой части описывает перенос капель со скоростью и(/,х,/и|^,...,//|-А)); третье слагаемое описывает увеличение или уменьшение размера
капель за счет фазовых переходов на поверхности капель и сорбции или десорбции
Л‘)(, т (0) (5) (0) (.?)\
примеси ЖИДКОСТЬЮ, при ЭТОМ 1 | /, X. / ,С ,...,С ,/Л — скорость измене-
ния массы /-го вещества в капле; первое слагаемое в правой части описывает диффузию капель как пассивной примеси, /}(/, х, Т, т) — коэффициент диффузии; второе слагае-
( Т (0) (5) (0) (^)'ч
мое описывает гомогенную конденсацию, ./ /,х,Т,с ,// —скорость
нуклеации капель; функционал а(/,х,/,/и^),...>(ц^)) описывает взаимодействие капель
— исчезновение мелких капель и увеличение количества крупных капель при коагуляции, функция С(!.х,т) описывает гравитационную седиментацию.
Решение уравнения (1.1) совместно с уравнениями движения, а также переноса водяного пара и газообразных примесей, позволит описать образование, перенос и выпадение токсичных осадков.
Фракционная модель
Далее, во избежание громоздких формул, рассматривается случай одной примеси (5 = 1), рассуждения для произвольного количества примесей 5 в пространственно од-
нородном случае без учета седиментации приведены в [7]. Следуя этой работе, функция / аппроксимируется кусочно-постоянной функцией по аргументам /и"' следующим образом:
/ (о) (1)Л (о) (о)у (0) (ОЛЦУ (1) (1)У (1) (1)
(2.1)
где /л':'} = 0 (р^'— масса вещества в критическом зародыше), Н() — ступенчатая фикция Хевисайда. Будем считать, что капли, удовлетворяющие неравенствам
(0) ^ (0) (0) (1) ^ (1) (1) „ , р < и <М у1 и р. - И < М^+}> принадлежат к одной фракции, при этом все капли, ~ ^ (1) принадлежащие к одной фракции, имеют одинаковую массу т - р + р . .
Используя гипотезу о взаимопроникающих континуумах и рассматривая каждую фракцию как «газ» одинаковых частиц, можно ввести распределенную плотность фракции, которая определяется соотношением
(о) 0)
Г Г , , (о), (1) , л ^
^=11УМ“ М = /*"м)И)‘
(') (0 (О
(°)А (О
(2.2)
(о) (1)
где Д и'-и' -и', т, .V V Лт
-/+1 О ’ ИХМ) У ('+1Х./+0 'У
Jj2 Далее
для простоты принимается
гипотеза об односкоростной однотемпературной среде, когда несущая газовая смесь и все «газы» частиц в каждой точке имеют равную скорость и температуру.
Умножая обе части уравнения (1.1) на массу капли т и и/ггегрируя по массе в пределах одной фракции, используя аппроксимацию (2.1), можно получить уравнение баланса распределенной плотности для каждой фракции:
р Р
— р + V • и р - О Ар +./ + сг - С н------------
3 и и ч ч ч ч ч т
(,+їХ^+ї) (,+іХ/+0/
(2.3)
1 ^ Ч’1 (:~чЫ))Р{'+^ Го (0) Ло)
ж
?=-1 Ар т, ч , і-('+9+ї)0+І)
о ^1(лА0) (1) ^ г(0) (1) И0) ^ V.
'Ы'+? (<+?)0'+і) 6\ (<+я1ЬН)^иЦ) ('+?Х>+і)^ч>)1
I
их.
І+Ч.
ж
\(лг( 1) (о) (1) (0) (1) (0)У
_ т!4/' - ^ -+/' ^ +/' ■ >" I
9=-1 Ар ' т, у .
.1+4 ('+]Х/+'3+0
Функция 2{х,у) = (- 1)Л Ь(1 - 2*^п(у))з§п(>') определяет направление процессов массо-обмена. Эта функция, в зависимости от аргументов, принимает значения
у < 0 ^ = 0 у > 0
х - -1 1 0 0
х = 0 -1 0 1
X = 1 0 0 -1,
что позволяет ей «включать» и «выключать» соответствующие слагаемые, описывающие обмен массой с другими фракциями.
Выражение в правой части (2.3) описывает увеличение числа капель в рассматриваемой фракции, в первой и второй сумме соответственно, за счет конденсации водяно-
го пара и поглощения примеси каплями меньшего размера (# = -1); испарение жидкости и выделение примеси каплями большего размера (ц - 1); при с/ = 0 обе суммы описывают уменьшение количества капель за счет перехода в другие фракции.
Как и в [7], коагуляция описывается с помощью модели Смолуховского. Для капель, состоящих из двух веществ, можно записать:
( (о) (1)\ ( (о) (1Пгг (0) (1) ~(о) -ОЛ ,/~(°) ,~0)
сги,М ,М 1 = -./1 ,М ,М ,М ,М ЩИ Ф +
/ \ ^0 0
: о) (II
1 МС ^Г~(°) ~(!) (°) ~(°) М ~(1'Л ,( (о) ~(о) (0 ,~(о) ,~(1)
] ] Киї ,ц ,/и -ц -ц ууі ,ц у\ц -/и ,ц -ц уц <1ц ,
где К — ядро коагуляции, которое определяется величиной вероятности столкновения капель различных размеров с последующим слипанием. Для простоты предположим, что для любых двух фракций ядро коагуляции остается постоянным в пределах каждой из фракций. Тогда для случая одной примеси [7]
V ( (о) (0)
Л' м р К 1 p+\q+\p+\q+\ р р К у/ <р +<р Ш \
_ ^ ^ ц pqiJ 1 V ' ^ ^ ^ V М ‘/^ \ ,}1 .Л J (2.4)
°~рц Ррч т 2 А (0)А ^ ^ ’
“" ИМ '■0^""т(«хЛАчН)4/‘' ^
где
ы 1 ц/ = — р* 2
^ (Ч) («)У (?)2 І?)2^ ( (?) , («)У (чр (д)2
■ (а ^ а ^2 ^2
+ пші(Ди. ,АМк ІМ3рЛ-М2рЛ
1 ( (<у)3 (Ч)3 (?)3 (?}3
+
(?) ( (?) (?)У (?) (?) ^ ( (<?) («)у (?) (?)
Рр1к=^ш+мк+1]{м4р1Гм3р1к)-1м, + мк ){м2р1к-м1р,к)
■ Г Л А ("ЛГ ^ ^ 1 ( № (>?)2 (?)2
+ шш А{1 ,Д// ы -и ---------\и -и -и +и
^ / к %р!к 2р\к у \pilc 2р:к Ъргк Ар1
(') ■ Г ^ (О (О О')^ (') ^ (0 ( (') • (>) (') (<■) у,
^ЗШП[1Ш^ '**1 >^^^4№"1ПаХ^, ’ИШ1[^+1+^+1’Л+1^
0) • ( ( (О . ( 0) (0 (') (')У (<п
и = тт тая и ,тт и +и , и + и \\,и 1.
21]к ^ ’ 'уО+1 к к+\ )) ^1+1 у
(') . ( ( (') ( (0 (>) (') (')у О)4^
^"тш[таЧ'а' ,тахК+1+^
Громоздкий вид второй суммы в правой части (2.4) обусловлен произвольностью выбора А// , однако поскольку ц не зависят от времени, выражение в скобках является постоянным множителем. Кроме того, большинство слагаемых во второй сумме остается нулевым на протяжении всего процесса.
Как указывалось выше, полученную систему дифференциальных уравнений необходимо дополнить уравнениями диффузии пара и примеси, которые с учетом высказанных предположений об односкоростной многоконтинуальной среде могут быть записаны в виде
^£_ +у. и с{,) = У-£>!і)Ус(,) +В{,) -11 У')ф(0,ф(1),
(') ■ п(')
где с - концентрация 1-го вещества, В - источники и стоки 1-го вещества
Используя аппроксимацию (2.1) для случая одной примеси и считая коэффициент диффузии постоянным по пространству, можно получить
Теперь, задаваясь конкретными зависимостями для скоростей массообмена, седиментации и нуклеации, а также граничными и начальными условиями, можно, пользуясь построенной фракционной моделью, получить решение рассматриваемой задачи.
Простейший случай
Для примера применения фракционной модели рассматривается решение одномерной задачи распространения токсичного аэрозоля, первоначально состоящего только из загрязняющего вещества, во влажной атмосфере в постоянном поле скорости В данном случае рассматривается плоский вертикальный слой высотой Н и длиной / , предполагается перенос и диффузия аэрозоля и пара только в горизонтальном направлении, совпадающем с направлением ветра, горизонтальные скорости аэрозольных частиц всех размеров, пара и воздуха равны и, однако вертикальная скорость оседания частиц V зависит от их массы. Концентрация пара и распределенные плотности фракций предполагаются однородными по высоте слоя. Нуклеация для простоты пренебре-гается.
Для определения скорости конденсации пара используется закон Максвелла,
ной температуре.
Поскольку используется односкоростное приближение, т е. относительные горизонтальные скорости отсутствуют, можно рассматривать только гравитационную коагуляцию. В этом случае [8]
сфере. Скорость оседания частиц каждой фракции под действием силы тяжести определяется законом Стокса,
где т] - вязкость воздуха, g - ускорение свободного падения.
Скорость уменьшения распределенной плотности фракции за счет осаждения частиц, с учетом однородности по высоте слоя, определяется по формуле [9]:
Поскольку при моделировании токсичных осадков наибольший интерес представляют область их выпадения и количество выпавшего вещества в различных точках
V * У
где р — плотность жидкости капли, с'г ' — равновесная концентрация пара при дан-
gm
V
У
V
подстилающей поверхности, имеет смысл ввести переменные р (х), описывающие
плотность выпавшей ')-той примеси на подстилающую поверхность в точке х . Соответствующее уравнение для единственной примеси в одномерном случае будет иметь вид
,(1) »г „ ..(О
а
К М
-=1Е<л
/=0 ;=0
т„
Построенная система уравнений дополняется начальными условиями
А)У(4 с(0)(х,О)=4О)(4 У1)(ж,0)=^1)(дс).
Для определенности предполагается,что и > 0, тогда граничные условия на входной границе (х = 0 ) задаются в виде
Р,) (о,о = Ру {Л с10)(0, ?) = с(О)(0,
а на выходной границе {х-Ь)
д п д (о)
-р= О,
дх
0.
Результаты
Для численного решения задача разбивается на две подзадачи. Первая, о переносе газов и частиц аэрозоля в пространстве, решается с помощью конечно-разностного метода с абсолютно устойчивой схемой. Вторая, о перераспределении массы между фракциями, газами и об осаждении примеси в каждой точке, решается с помощью метода Рунге-Кутта.
В качестве иллюстрации полученных результатов приводится решение задачи с
/ \ (0 { \ (о)
граничными условиями /?0|(0,/)-М е , рхо[0,1) = М ,а для остальных фракций
р., (0. /) — 0. Такое задание условий соответствует кратковременному выбросу монодисперсно й аэрозольной примеси во влажную атмосферу от наземного источника ( Н = 1).
На рис. 2 представлен график эволюции зависимости плотности выпавшей примеси от расстояния. Отчетливо просматривается конкуренция двух пиков, вызванная ростом частиц. Левый пик, ярко выраженный на ранних стадиях процесса, обусловлен высокой плотностью мелких примесных частиц, поступающих через входную границу. С течением времени поступление частиц снижается, однако за счет конденсации и коагуляции появляются более крупные частицы, оседание которых протекает более интенсивно, что приводит в конечном счете к росту правого пика
Рис. 2. График плотности выпавшей примеси в различные моменты времени: а) при t -5с ; б) при ( = 10 с, в) при < = 15 с; г) при / = 20 с
На рис. 3 приведены распределенные плотности фракций в различных точках области. В начальной области (рис. 3,а) присутствуют чисто водяные капли, но уже к середине области (рис. 3,6) они исчезают образуя, смешанные капли, при этом спектр частиц остается практически неизменным до конца области (рис. 3,в).
; 1*7 и'!’
Р
81ИШМ У *'
V/,
«г
□«
щи1-***-11
Ї4ІЗ«.11
Рис. 3. Фракционный состав при I = 15 с в точках, а) х = Ь/4\б) х = Ь/2; в) х = ЗЫ4
Выводы
Предложена модификация фракционной модели, которая позволяет описывать процессы осаждения токсичного аэрозоля на подстилающую поверхность. Рассмотрен достаточно простой случай распространения аэрозольной примеси в атмосфере. Несмотря на простоту задачи, решение позволило выявить влияние деформирующегося спектра частиц на вид зависимости плотности выпавшей примеси от расстояния. В начальной стадии процесса более интенсивно осаждаются мелкие частицы, фракция которых имеет большую распределенную плотность. На более поздних этапах доминирующим становится осаждение крупных смешанных частиц, появление которых обусловлено коагуляцией и конденсацией водяного пара на поверхности частиц.
Библиографический список
1. Степанов А С,, Захарова И.М., Новикова Л.Д. Моделирование процессов накопления загрязняющих веществ в каплях туманов // Метеорология и гидрология. - 1997 -№4 - С. 25-36.
2. Московченко А.В., Гладков С О. Динамика химического состава аэрозолей в условиях коагуляции // Доклады АН. - 1993. - Т. 328. - №3. - С. 325-329,
3. Пискунов В.Н., О кинетике коагуляции композитных частиц // Доклады АН. - 1996 -Т. 350. -№1. - С. 70-73.
4. Trautmann Т., Wanner С. A fast and efficient modified sectional method for simulating multicomponent collision kinetics // Atmospheric Environment. - 1999. - Vol. 33. - P. 1631-1640.
5. Russel L.M., Seinfeld J.H. Size- and Composition-Resolved Externally Mixed Aerosol Model // Aerosol Science and Technology. - 1998, - Vol. 28. - №5. - P. 403—416.
6. Алоян A.E. Зайганов В.А. Лушников A.A. Макаренко С.В. Перенос трансформирующегося аэрозоля в атмосфере // Известия АН СССР. Физика атмосферы и океана. - 1991. -Т. 27. -№11. -С. 1232-1240.
7 Харченко А. В., Бояршинов М. Г., Трусов П В., О фракционной модели накопления загрязняющих веществ в тумане // Математическое моделирование систем и процессов: Сб. науч. тр./ - Перм. гос. техн.ун-т; Пермь, 1999. - №7 - С. 92-99.
8. Ломая В.А., Мазин И.П., Неизвестный А.И. Влияние турбулентности на эффективность коагуляции облачных капель // Известия АН СССР. Физика атмосферы и океана. - 1990. - Т. 26. -№8. - С. 813-819.
9, Райст П. Аэрозоли: Введение в теорию - М.: Мир, 1987 - 280 с.
Получено 20.03 ,2000
