УДК 519.6
Н. Н. БЕЛЯЕВ, О. П. ШАКИНА, А. Ю. ШМИГЕЛЬ (ДИИТ)
ЗАГРЯЗНЕНИЕ ВОЗДУШНОЙ СРЕДЫ ПРИ ГОРЕНИИ ФОСФОРА
Розроблено тривимiрну чисельну модель розрахунку процесу поширення токсично! рiчовини при ropiHHi фосфору. Модель базуеться на чисельному iнтегруваннi рiвняння конвективно-дифузiйного переносу домишки. Наводяться результати обчислювального експерименту.
Разработана трехмерная численная модель распространения оксида фосфора в атмосфере при горении фосфора. Модель основывается на численном интегрировании уравнения конвективно-диффузионного переноса примеси. Приводятся результаты вычислительного эксперимента.
The 3D numerical model to simulate the toxic gas dispersion during phosphor firing.The model is based on the K-gradient transport model. The results of numerical experiment are presented.
Введение. Авария при транспортировке фосфора железнодорожным транспортом, которая произошла на Львовской железной дороге в июле 2007 года поставила задачу разработки более надежных, чем существующая нормативная методика [5], методов прогноза уровня загрязнения окружающей среды при горении фосфора. Нормативная методика дает возможность только оценить размеры (площадь) зоны загрязнения и только до высоты 10 метров от уровня земли и не позволяет рассчитать концентрацию оксида фосфора (продукта горения фосфора) на различном расстоянии от места аварии. Особую актуальность приобретает вопрос разработки эффективных методов прогноза уровня загрязнения атмосферы при авариях связанных с эмиссией токсичных веществ с целью оперативного получения научно обоснованной информации об уровне загрязнения воздушной среды, подстилающей поверхности для организации эффективной защиты населения от химического поражения. Целью настоящей работы является создание регуляторной математической модели, ориентированной на использование экологами в повседневной практике для экспресс-прогноза загрязнения атмосферы при горении фосфора.
Достоинством предложенной модели является возможность учета основных физических факторов, влияющих на процесс переноса оксида фосфора в атмосфере и при этом требуются небольшие затраты машинного времени при практической реализации модели.
Постановка задачи. Будем считать, что произошла разгерметизация цистерны с фосфором, в результате которой имело место испарение воды, находящейся в цистерне. Это привело к тому, что фосфор вступил в контакт с воздушной средой и начался процесс его горения:
4P + 5O2
•2PA.
В результате горения от цистерны начинает распространяться шлейф оксида фосфора, представляющего угрозу для окружающей среды.
Для моделирования процесса переноса оксида фосфора в атмосфере будем использовать трехмерное уравнение миграции примеси [2, 3]
дС диС д(^ - ^)С — +-+-+ —-—— =,
дt дх ду д2
ду дС_
дх I х дх
дх д
дС
ду \ у ду
f оц f)+1 &(> )«( - r)
(i)
где С - концентрация загрязняющего вещества (Р205), попадающего в атмосферу при горении фосфора;
и, V, ^ - компоненты вектора скорости воздушной среды;
- скорость оседания примеси;
ц = (цх, , ) - коэффициент турбулентной диффузии;
Q - интенсивность выброса токсичного вещества;
5 (г - г ) - дельта-функция Дирака;
Г = (С, у^, ) - координаты источника выброса (цистерны).
В данной модели имеется возможность задавать несколько источников выделения оксида фосфора, т. е. моделировать поступление оксида фосфора в атмосферу при горении от нескольких цистерн.
В разработанной математической модели
используется неравномерный профиль скорости ветра и вертикального коэффициента диффузии
и = и,
I 2 >
V 21 У
Ц 2 = к1
I 2 >
V 21 У
Конвективные производные представим в виде:
диС ди+С ди - С
дх дх дх
Цх Ц
Ц 2 х =-=
где и1 - скорость ветра на высоте 21; к^0,2; п = 0,16; т «1.
Для моделирования процесса загрязнения атмосферы при горении фосфора необходимо знать массовую скорость выгорания фосфора Gг, кг/м2-мин. Как правило, массовая скорость выгорания того или иного вещества определяется экспериментальным путем. Для экспресс-прогноза эту величину можно ориентировочно назначить исходя из имеющихся литературных данных [2]. Так как молекула фосфора имеет слабые связи, то можно считать, что массовая скорость выгорания его близка к массовой скорости выгорания бензина. Тогда, зная геометрические размеры цистерны, площадь ее поперечного сечения S можно рассчитать интенсивность Я поступления оксида фосфора в атмосферу
Я() = 5(г)в,
где 1 - время.
Для защиты атмосферы используют подачу воды в шлейф оксида фосфора. В этом случае процесс взаимодействия описывается уравнением
Р205 + 3Н20 ^ 2Н3Р04.
В разработанном пакете прикладных программ процесс взаимодействия оксида фосфора с водой рассчитывается в отдельном модуле. Процесс рассеивания воды от источника (в данной работе рассматривается от вертолета) моделируется с помощью уравнения виды (1). В модели полагается, что капли воды осаждаются со скоростью .
Метод решения. Рассмотрим построение попеременно-треугольной разностной схемы для численного интегрирования уравнения переноса. Заменим производную по времени разделённой разностью «назад»:
дС сЩ - СП
л и Ах . и Ах
1--1 +-
Ц Ц
дуС ду+С ду~С
где и =
ду ду ду
дл>С дм>+С дн>~ С -=-+-
д2 д2 д2
и + и - и + и + V + IV -—; и =—!—1 , V =■
V - V
V = -
2
w+ =
2 2
w +- w - |w| ■; w =-
2
2
Для аппроксимации конвективных производных используем выражения:
ди С _ и-1, Щк С"+1- Щ,к - СЩ = Т- С П+1 .
д х Ах
ду+С ^ ^ ,Щ+1,кСЩк - V+кСг,Щ-1,к =
ду Ау
дv-C ^ Vi-Щ+1,кС1,Щ+1,к - УЩкСЩк
ду Ау
ды+С ^ К],к+1СЩк - Щ,к-1
д2 А2
1 Wi,],к+1Сг,],к+1 - WijkCi,],к _
д2
А2
г+сп+1
= Т-Сп+1.
Компоненты скорости потока будем задавать на гранях разностной ячейки (контрольного объема). Вторые производные аппроксимируем следующими выражениями:
<ч / Л Г>п+1 пп+1
д I дС ] _ Ц+1,щ,к - сук
дх Г х &х
Ах
у^п+1 у^п+1
и>к - С-1-м = м-хсп+1 + мххсп+1
~Цх Л 2 Ах
хх '
<ч / Г"п+1 Г'п+1 пп+1 пп+1
д I дС \ _ Ц,Щ,к+1 - Сщк _ ЧЩ,к - Су,к-1
&1ц 2 ^ 2 А22--ц 2 А22
= М-Сп+1 + М+Сп+\
дг
Аг
2
В данных выражениях Lx , Ьх , Ly, Ц , Ц, Ц,
M+xx, М—х и т. д. - обозначения разностных операторов. С учетом этих обозначений, разностный аналог трехмерного уравнения переноса примеси будет иметь вид:
С п+1 сп
Ик - г]к + цСпх1 + ц —пх1 + ЬУССпх1 +
А!
У
+Цу Спх1 х Ц Спх1 х Ц Сп+1 хоС^1 = = (м+хСпх1 + Мух—пх1 ххЫ+уСпх1 х +МуУСпх1 хМ+Спх1 хМуСпх1).
УУ 22 22 I
Расщепим решение данного разностного уравнения при интегрировании на временном интервале & так:
к 1
- на первом шаге к = — :
о ^ к ']к
г^пхк _п 1
—————— х1 ЦСк хЦСк хЦСк)х-—
^ 2v х у 2 ' 4 !
/ 1 1 - на втором шаге к = п х—; с = п х— :
2 4
С— — —— х1 (Ск+ЦСк хЦ Ск) + А 2v х у '
о —— = 1 (м—х—к х МХх—с х М—у Ск х
хМУуСс х Муг Ск х МХ2 Сс);
1 3 1
- на третьем шаге к = п х— ; с = п х— ис-
4 2
пользуется формула (4);
1 Л 3
- на четвертом шаге к = п х 1; с = п х используется формула (3).
В данные разностные соотношения входит «возмущенный» коэффициент диффузии Д :
=
Ц х
1х
и Ах Ц
х =•
Ц х
1-
и Ах Ц
и т. д.
На пятом шаге расщепления учитывается влияние источников на изменение концентрации примеси и расчетная зависимость на этом шаге имеет вид:
5 пх1 5
С г, — ,к — С
N „ />пх1/2\ г,—,к =^ ql (t ' ) 5
Аt 1=1 АхАуА2
В дискретном виде дельта-функция Дирака «размазывается» по объему разностной ячейки с учётом сбережения суммарного количества выбрасываемого загрязнения. Функции 51 тождественно равняются нулю, кроме ячеек, где расположен 1-й источник загрязнения. На каждом шаге расщепления расчет неизвестной концентрации оксида фосфора осуществляется по явной формуле бегущего счета.
Практическая реализация. На основе разработанной численной модели создан пакет прикладных программ СНБМ1САЬ-3Б. Разработанная численная модель была применена для решения задачи о загрязнении атмосферы при поступлении оксида фосфора от двух железнодорожных цистерн. Численное моделирование выполнено при следующих параметрах: размеры расчетной области 250^221^55 м, коэффициент диффузии - 2,35 м2/с, скорость ветра - 3,3 м/с (на высоте 10 м), аварийный выброс происходит на высоте 1,3 м, профиль ветра -неравномерный; скорость оседания капель воды - 0,01 м/с. Для защиты атмосферы от загрязнения осуществляется подача воды от вертолета на высоте 9 м, интенсивность подачи воды - 150 кг/сек, вода подается в течение промежутка времени (100...116 с).
На рис. 2.4 показана зона загрязнения атмосферы возле цистерн для различных направлений ветра.
Рис. 1. Выход оксида фосфора из цистерны
0.216Е+03
—0.124Ет с 0.116Е<00 г _и. 1 иг:г I ип л
П.1(№ИП 1 —0.920Е-01 л 0.840Е-01 а —и,7Ь0Е-01 1 0.680Е-01 е —Я,ЬВЙЕ-Й1 0.520Е-01 ——е.440Е-01 у
0.360Е-01 —И.З¥1Е-И1 0.200Е-01 —Й.120Е-01 0.400Е-02
0.526Е+01
0.500Е»01 соогА1па!б х 0.245Е+03
Рис. 2. Зона загрязнения атмосферы при выбросе оксида фосфора из двух цистерн г = 70 с (сечение 2 = 9 м), а = -45°
0.216Е+03
-_0.6 20Е-01 о 0.580Е-01 г —0.540Е-01 ¿1 0.500Е-01 ! —0.460Е-01 и 0.420Е-01 а —0.380Е-01 1 0.340Е-01 е ^^.300Е-01 0.260Е-01 —0.220Е-01 у 0.180Е-01
0.100Е-01 ^^.600Е-02 0.200Е-02
0.526Е+01
0.500Е+01 1:11111'(111м 1.н к 0.245Е+03
Рис. 3. Зона загрязнения атмосферы при выбросе
оксида фосфора из двух цистерн г = 70 с (сечение 2 = 9 м), а = 45° , * - место сброса воды от вертолета
0.216Е+03
_0.878Е*00 □ 0.822Е»00 г _0.7Б5Е*00 А 0.708Е+00 1 _0.&52Е*00 п 0.595Е+00 а —0.538Е*00 г 0.482Е*00 е -0.42СЕ '30
0 ::С.ОЕ |0О
—0.312Е*00 |
0.255Е+00 —в.ЮОЕ'ОО 0.142Е+00 ^.850Е-01 0.283Е-01
0.526Е+01
0.500Е+01 ¡и.11.С х 0.245Е+03
Рис. 4. Зона загрязнения атмосферы при выбросе оксида фосфора из двух цистерн г = 70 с (сечение г = 1,3 м), а = 45°
Из рис. 2, 3 отчетливо видно, как происходит переориентация зоны загрязнения при изменении направления ветра. Так при угле а = -45° отчетливо выделяется две подзоны с высоким градиентом концентрации оксида фосфора, которые соответствуют положению двух железнодорожных цистерн. Однако, при угле а = 45° эти две подзоны сливаются в одну зону загрязнения на высоте 9 м. На высоте 1,3
м (рис. 4) при данной ориентации ветра можно отчетливо видеть положение дымящихся цистерн. На последующих рисунках показана динамика изменения формы, размеров и интенсивности зоны загрязнения при подачи воды от вертолета. Как видно из рис. 5 в первые секунды после окончания подачи воды зона загрязнения резко уменьшается в размерах (отметим, что подача воды прекращается в момент времени г = 116 с). Однако, достаточно быстро зона загрязнения «восстанавливается» в прежних размерах. Это отчетливо видно из рис. 7, который необходимо сравнить с рис.3, где представлена зона загрязнения перед подачи воды с вертолета, т. е. одноразовая подача воды является не слишком эффективным способом защиты атмосферы от загрязнения.
0.21ЬЕ»03
_0.310Е-01 о 0.290Е-01 г _0.270Е-01 А 0.250Е-01 1 _0.230Е-01 п 0.210Е-01 а —а.19еЕ-в11 0.170Е-01 е —Г1.15ЙЕ-(1| 0.130Е-01 —0.110Е-01 у 0.900Е-02
—елтк-ег
0.500Е-02 —Е1.300Е-82 0.100Е-02
0.Е2ЬЕ»01
0.500Е+01 сооНгпа^с х 0.245Е+03
Рис. 5. Зона загрязнения атмосферы при выбросе
оксида фосфора из двух цистерн г = 120 с (сечение 2 = 9 м), а = 45° (после подачи воды)
0.21ЬЕ»03
—0.&20Е-01 □ 0.580Е-01 г —0.540Е-01 а 0.500Е-01 ; —8.41Ж 01 г, 0.420Е-01 а —0.380Е-01 1 0.340Е-01 е ^.300Е-01 0.260Е-01
-0.22аЕ-И1 у
0.180Е-01
-0.140С 01
0.100Е-01 ^.600Е-02 0.200Е-02
0.52ЬЕ*01
0.Б00Е+01 соогД1т1а1е х 0.245Е+03
Рис. 6. Зона загрязнения атмосферы при выбросе
оксида фосфора из двух цистерн г = 160 с (сечение 2 = 9 м), а = 45° (после подачи воды)
ш
0.200Е-02
0.526Е+01 ____
0.500Е+01 г.ппгИ у. 0.245Е+03
Рис. 7. Зона загрязнения атмосферы при выбросе
оксида фосфора из двух цистерн t = 200 с (сечение 2 = 9 м), а = 45° (после подачи воды)
Ниже представлены данные об интенсивности нейтрализации оксида фосфора при подачи воды:
t = 120 с фосфора
t = 140 с фосфора
t = 160 с фосфора
t = 160 с фосфора
t = 200 с фосфора
t = 260 с фосфора
t = 280 с фосфора
Из этой таблицы видно, что процесс нейтрализации закончился к моменту времени t = 260 с.
Выводы. В работе построена трехмерная численная модель расчета загрязнения воздушной среды в случае горения фосфора. Модель основывается на численном интегрировании уравнения переноса примеси в атмосфере. Разработанная модель позволяет выполнить расчет процесса нейтрализации оксида фосфора (продукт горения фосфора) путем подачи воды от вертолета. Модель дает возможность оперативно рассчитать эффективность процесса нейтрализации с учетом расположения дымящихся цистерн, места, интенсивности и режима подачи воды, метеоусловий. Расчет по предложенной модели требует около 5 с на получение прогнозной информации об уровне загрязнения атмосферы и эффективности организованного процесса нейтрализации. Для практической реализации модели необходима стандартная метеорологическая информация. Дальнейшее совершенствование модели необходимо проводить в направлении ее адаптации к моделированию процесса нейтрализации оксида фосфора с учетом рельефа местности.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Котляревского В. А. Аварии и катастрофы. Предупреждение и ликвидация последствий: Учебн. пособ. В 5-ти кн. / Под ред. В. А. Котляревского и А. В. Забегаева. - М.: АСК, 2001. -200 с.
2. Згуровский М. З. Численное моделирование распространения загрязнения в окружающей среде / М. З. Згуровский, В. В. Скопецкий, В. К. Хрущ, Н. Н. Беляев - К.: Наук. думка, 1997. - 368 с.
3. Марчук Г. И. Математическое моделирование в проблеме окружающей среды. - М.: Наука, 1982. - 320 с.
4. Методика прогнозування насладив впливу (ви-киду) небезпечних х1м1чних речовин при авар1ях на промислових об'ектах 1 транспорт!. -К.: Наук. думка, 2001. - 33 с.
Поступила вредакцию 19.09.2007.
нейтрализовано нейтрализовано нейтрализовано нейтрализовано нейтрализовано нейтрализовано нейтрализовано
8,26 кг оксида 12,41 кг оксида 14,38 кг оксида 15,49 кг оксида 16,14 кг оксида 16,45 кг оксида 16,45 кг оксида