УДК 519.6
Н. Н. БЕЛЯЕВ, О. А. ВЕРГУН, Ю. В. КОВТУН (ДИИТ)
МОДЕЛИРОВАНИЕ ПРОЦЕССА НЕЙТРАЛИЗАЦИИ ОКСИДА ФОСФОРА В АТМОСФЕРЕ
Запропоновано метод розрахунку процесу нейтралiзацiï токсичного газу в атмосфера В якосл математи-чно1 моделi використовуеться двовимiрне рiвняння переносу домiшок. Розрахунок здiйснюeться з викорис-танням неявно1 рiзницевоï схеми.
Предложен метод расчета процесса нейтрализации токсичного газа в атмосфере.. Метод основывается на численном интегрировании двухмерного уравнения конвективно-диффузионного переноса примеси. Для расчета используется неявная разностная схема. Приводятся результаты вычислительного эксперимента.
A numerical model was used to calculate the process of toxic gas neutralization. The 2D transport model is used to simulate the process. The implicit difference scheme was used. The rsults of numerical experiment are presented.
Одним из эффективным методов защиты атмосферы от загрязнения при миграции токсичных веществ является метод нейтрализации. При разработке плана ликвидации аварийной ситуации весьма важно иметь простые и эффективные методы расчета процесса нейтрализации, чтобы оперативно принять научно обоснованное решение по размещению технических средств, подающих нейтрализатор в шлейф токсичного газа. Следует иметь в виду, что данная задача является достаточно сложной с точки зрения теоретического решения, поскольку в атмосфере происходит процесс переноса токсичного газа, нейтрализатора и одновременно процесс их химического взаимодействия. В этой связи основным методом решения становится численное моделирование (СББ модели). В рамках данного направления, для практики, имеет смысл разработка регулятор-ных численных моделей, т.е. моделей, учитывающих, с одной стороны, основные физические процессы переноса, а с другой стороны, требующих небольших затрат машинного времени при практической реализации и незначительного количества исходных данных (вполне доступных) для моделирования. Поэтому целью данной работы явилась разработка эффективного метода расчета процесса нейтрализации токсичного вещества в атмосфере при его аварийной утечке с помощью подачи нейтрализатора от вертолета (движущегося или неподвижного).
Математическая модель
Будем рассматривать аварийную ситуацию, в результате которой в атмосферу происходит эмиссия токсичного газа - Р2О5 (продукт горе-
ния фосфора) Для нейтрализации токсичного газа используется подача воды в шлейф Р2О5. Процесс взаимодействия токсичного газа и нейтрализатора рассчитывается на основе уравнения реакции
Р2О5 + ЗН2О ^ 2Н3РО4.
Чтобы рассчитать это взаимодействие, а значит определить эффективность процесса нейтрализации, необходимо знать концентрацию воды и Р2О5 в каждой точке пространства, где возможно это взаимодействие. Вода и токсичный газ рассеиваются в атмосфере под действием ветра, атмосферной диффузии, поэтому их концентрация изменяется под действием этих факторов в каждой точке. Для расчета переноса этих веществ в атмосфере воспользуемся моделью градиентного типа, осредненной по высоте переноса примеси (токсичный газ, нейтрализатор) [3]:
dC duC dvC
dt dx
dy
- cC = div (| gradC ) -
N
+1 & (Щх - хг (0)5(у - у (Г)), (1)
1=1
где С - концентрация примеси в атмосфере; и,V - компоненты вектора скорости воздушной
среды; ц = х, цу) - коэффициенты турбулентной диффузии; & - интенсивность выброса примеси (токсичного газа, нейтрализатора); 5(х - хг (^))5(у - уг (^)) - дельта-функция Дирака; хг, уг - координаты источника эмиссии примеси, зависящие от времени; с - коэффициент, учитывающий оседание примеси, химический
© Беляев Н. Н., Вергун О. А., Ковтун Ю. В., 2009
распад загрязнителя, характер взаимодействия с подстилающей поверхностью [3]; t - время.
Постановка краевых условий для задачи переноса рассмотрены в [2, 3].
Метод решения
Для численного интегрирования уравнения (1) используется прямоугольная разностная сетка. Уравнение (1) интегрируется с использованием попеременно-треугольной неявной разностной схемы расщепления. Выполним следующее физическое расщепление уравнения переноса примеси при интегрировании на временном шаге &:
д С д иС д уС
д t д х д х
+ с С _ 0;
дС дС^
д t д х
д х
д С ц-
д у У д у
д С д t
_! а М Щг - г),
где гг = (хг , уг) - координаты источника эмиссии примеси.
Первое уравнение из данной системы описывает процесс переноса примеси под действием ветра, второе уравнение - под действием диффузии, а третье - изменение концентрации примеси под действием источников эмиссии.
Выполним следующую замену производных:
4 К СП+1 - СП; дt ~ At '
диС ди+С ди - С
дх дх дх дуС ду+С сУ~С
ду ду ду и + \и\ и - и
V + V
где и _ -
2
и _-
2
V _■
V- V
V _ ■
2
д^+С ^ V+1-V +г]Сг,]-1 _ 1+Сп+1
ду Ау у '
ЬЦС «Д 1 +1Сг,1 +1 -1 _ Ь-СП+1 и т.д. ду Ау у
Вторые производные аппроксимируем так:
У^П+1 У^П+1
4+1,1 - Су
дСЛ 1 ~
дх у |*Ц х1
^п+1 Ч1 ~ у^п+1 ■ Сг -1,1
Ах2
Ах2
_ М:хСп+1 + М+хСп+1;
дС_
ду У у ду
/~т+\ _+1
Сг, 1+1 ~ Су
Ау2
СС п+1 сс п+1
-Ц г■1 -2 г-1 -1 _МууСп+1 + М+Сп+\ 2 Ау 2
~ 2ЦгЦг+1 где ц _-и т.д.
1 Цг + Ц+1
Здесь Ь+, Ь-, М+, Мх и т.д. - обозначения разностных операторов.
Решение исходного уравнения разбивается на ряд шагов:
, 1
- на первом шаге к _ — :
4
п+к п л
Ч - Ч 1
А t
+ 2 (Ь+Ск + 1ЬуСк ) сп _ 0;
- на втором шаге к _ п +—; с _ п +— :
2 4
Са - си 1
А г
+ - (4 + Ь-уСк ) + - Ск _ 0;
- на третьем шаге к _ п +— ; с _ п +— :
4 2
С к 4с 1
_1 (сс + М+Ск + МуСс + М+уСк );
- на четвертом шаге к _ п + 1; с _ п
Ск - Сс 1 At 2
(Ск + М+ Сс + М~ Ск + М+ сс)
\ хх хх уу уу /
На пятом шаге расщепления учитывается влияние источников на изменение концентрации примеси и расчетная зависимость на этом шаге имеет вид:
5 п+1 5 п
С г, 1 - С г, 1 ( 1 ) 2
At
ы АхАу
2
В дискретном виде дельта-функция Дирака «распределяется» по объему разностной ячейки с учётом сохранения суммарного количества выбрасываемой примеси. Функции 5г тождественно равняются нулю, кроме разностных ячеек, где расположен источник выброса примеси.
Результаты моделирования
Построенная численная модель реализована в виде пакета программ «AIR-2D» (программы записаны на алгоритмическом языке FORTRAN). Применение этого пакета программ для решения конкретной задачи рассматривается ниже.
Рассматривается следующий сценарий: происходит аварийный выброс P2O5 в количестве Q = 1 (этот и другие параметры задачи приводятся далее в безразмерном виде). На месте выброса образовался шлейф (рис. 1). В момент времени t = 15 начинается подача воды (нейтрализатора) в шлейф в количестве Qn = 22.2. Подача воды заканчивается в момент времени t = 25. Скорость ветра - 7 м/с; коэффициент диффузии - 4.7 м2/с; коэффициент с = 0.01. Размеры расчетной области 1500 х 800 м.
На рис. 1 показана зона загрязнения атмосферы токсичным газом непосредственно перед началом нейтрализации. Как видно из данного рисунка, возле места аварии образовался шлейф, вытянутый в направлении движения воздушного потока. На последующих трех рисунках показано, как изменяется зона загрязнения атмосферы при подаче нейтрализатора от неподвижного вертолета. Видно, что первоначально происходит «разрыв» шлейфа в том месте, где осуществляется подача нейтрализатора, т.е. образуется зона свободная от токсичного газа (рис. 2). Затем зона загрязнения локализуется только возле источника эмиссии токсичного газа (t = 25, рис. 3). Но т.к. подача нейтрализатора прекращается в момент времени t = 25, а источник эмиссии токсичного газа продолжает действовать, то, как видно из рис. 4, шлейф токсичного газа практически в течение короткого промежутка времени после прекращения подачи воды «восстанавливается» в размерах.
На последующих рисунках (рис. 5 - 7) представлена динамика изменения зоны загрязнения при подаче нейтрализатора из движущегося вертолета. Представленные результаты моделируют «неудачный» маршрут вертолета, когда он движется вдоль границы шлейфа. Тем не менее, как видно из результатов моделирования, и в этом случае происходит изменение, хо-
тя и незначительное, формы зоны загрязнения атмосферы вследствие взаимодействия токсичного газа и нейтрализатора.
Рис. 1. Зона загрязнения атмосферы после аварии (вид сверху), t = 15 (нет подачи нейтрализатора)
Рис. 2. Зона загрязнения после аварии (вид сверху), t = 20 (идет подача нейтрализатора, вертолет не движется)
Рис. 3. Зона загрязнения после аварии (вид сверху), t = 25 (идет подача нейтрализатора, вертолет не движется)
Выводы
Рассмотренный в данной работе метод расчета нейтрализации токсичного газа в атмосфере является эффективным инструментом решения комплекса задач при при разработке плана ликвидации аварийной ситуации. На решение задачи с помощью модели требуется около 3.. .5 секунд машинного времени, что позволяет в течение короткого промежутка времени выполнить расчет десятков вариантов организа-
ции подачи нейтрализатора. При практическом применении разработанного пакета программ не требуется каких-либо дополнительных специальных сервисных программ. Разработанный пакет программ имеет удобную для пользователя систему ввода исходной информации и вывода получаемых результатов. Дальнейшее совершенствование предложенного метода необходимо вести в направлении создания пространственной модели переноса примеси.
Рис. 4. Зона загрязнения после аварии (вид сверху), t = 30 (подача нейтрализатора прекращена)
Рис. 5. Зона загрязнения после аварии (вид сверху),
t = 16 (идет подача нейтрализатора, вертолет движется)
Рис. 7. Зона загрязнения после аварии (вид сверху), t = 24 (идет подача нейтрализатора, вертолет движется)
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Аварии и катастрофы. Предупреждение и ликвидация последствий [Текст] : учеб. пособие в 5 кн. / под ред. В. А. Котляревского и А. В. Забе-гаева. - М.: Изд-во АСВ, 2001. - 200 с.
2. Численное моделирование распространения загрязнения в окружающей среде [Текст] / М. З. Згуровский и др. - К.: Наук. думка, 1997. - 368 с.
3. Марчук, Г. И. Математическое моделирование в проблеме окружающей среды [Текст] / Г. И. Марчук. - М.: Наука. - 1982. - 320 с.
4. Самарский, А. А. Теория разностных схем [Текст]. - 2-е изд., испр. / А. А. Самарский. -М.: Наука, 1983. - 616 с.
Поступила в редколлегию 18.03.2009.
Принята к печати 02.04.2009.
Рис. 6. Зона загрязнения после аварии (вид сверху), t = 18 (идет подача нейтрализатора, вертолет движется)