CC BY
37
удк 544.43
Е. В. Степашина (асп.), А. И. Байтимерова (к.ф.-м.н., доц.), С. А. Мустафина (д.ф.-м.н., проф.)
Программный комплекс автоматизации процедуры уточнения механизма химической реакции на основе DRGEP-метода
Стерлитамакская государственная педагогическая академия им. Зайнаб Биишевой, кафедра математического моделирования 453103, г Стерлитамак, пр. Ленина, 37; тел./факс. (347) 3431039, e-mail: [email protected]
E. V. Stepashina, A. I. Baitimerova, S. A. Mustafina
The software package for the computer-aided procedure of qualified mechanism of chemical reactions based on DRGEP-method
Sterlitamak State Pedagogical Academy named Zaynab Biisheva 37, Lenina pr, 453103, Sterlitamak, Russia; ph./ fax. (347) 3431039, e-mail: [email protected]
Для автоматизации процедуры уточнения механизма химической реакции разработан программный комплекс. Данный программный продукт позволяет осуществить построение кинетической модели, решить прямую задачу химической кинетики и на основе DRGEP-мето-да исключать из механизма вещества, мало влияющие на образование целевых веществ.
Ключевые слова: граф реакции; DRGEP-ме-тод; механизм реакции.
Кинетические модели, основанные на детальных механизмах сложных химических реакций, как правило, представляют собой системы дифференциальных уравнений большой размерности. Размерность определяется списком веществ, который часто превышает 102— 103 реагентов (например, нефтехимические процессы, некоторые реакции горения). Для анализа механизма реакции большой размерности требуется точное описание поведения лишь нескольких веществ, которые называются целевыми. Возникает задача замены исходной системы системой меньшей размерности, в каком-то смысле эквивалентной исходной, сохраняющей при этом изменение концентраций целевых веществ.
За последнее десятилетие предложено большое количество методов, нацеленных на сокращение механизмов реакций. Один из подходов к решению этой задачи — анализ матрицы чувствительности веществ к константам скоростей элементарных реакций. Методами, реализующими данный подход, являются метод анализа главных компонент матрицы чув-
Дата поступления 17.05.11
The software package was designed for the computer-aided procedure of qualified mechanism of chemical reactions. The package allows to carry out construction of kinetic models, to solve a direct problem of chemical kinetics, and to exclude substances a little influencing formation of target substances on the basis of DRGEP-method.
Key words: graph of the reaction; DRGEP-method; reaction mechanism.
ствительности (Principial Components Analysis (PCA)) 1 и прямой анализ чувствительности (Direct Sensitivity Analysis (DSA)) 2. Для того, чтобы применять эти алгоритмы, необходимо рассчитать матрицу чувствительности во всех анализируемых моментах времени реакции. Расчет такой матрицы представляет собой отдельную, достаточно ресурсоемкую задачу.
В настоящее время широко применяется геометрическая трактовка механизма реакции. На основе теории двудольных графов возможен анализ существования и свойств решения системы обыкновенных дифференциальных уравнений кинетической модели схемы реакции .
Одним из методов, реализующих графовый подход к механизму реакции, является метод анализа графа прямых связей (Direct Relation Graph (DRG)) 4. Для того, чтобы применять этот алгоритм, необходимо знать скорости всех реакций в анализируемые моменты времени, то есть достаточно иметь только численное решение кинетической задачи. Как правило, от упрощенного механизма требуется точное описание поведения лишь нескольких
веществ. Поэтому выбирается набор целевых веществ, который является стартовым набором для итерационной процедуры построения упрощенного механизма. Хотя метод DRG рассматривает только прямые связи между веществами, тем не менее, вещества могут быть связаны через промежуточные вещества, то есть косвенно. В связи с этим на основе метода DRG был разработан метод анализа графа прямых связей с распространением ошибки (Direct Relation Graph with Error Propogation (DRGEP)) 5, учитывающий косвенное влияние веществ. По сравнению с методом DRG, метод DRGEP позволяет более точно сократить механизм, так как в схеме реакции связь между косвенно связанными веществами может превосходить прямую зависимость между ними.
Для автоматизации процедуры уточнения механизма химической реакции в среде визуального программирования Delphi на языке Object Pascal был разработан программный комплекс «ChemReductor». Данный программный продукт позволяет осуществить построение кинетической модели, решить прямую задачу химической кинетики и на основе DRGEP-метода исключить из механизма вещества, мало влияющие на образование целевых веществ.
Постановка задачи
DRGEP—метод сокращения схемы химической реакции, использует коэффициент зависимости между веществами rAB .
Е \УАРг$в
Г =
' AB
i=j
Е \УАг®г
ГАС ГАВ Х ГВС ,
(2)
где гАВ — коэффициент зависимости между веществами А и В;
гВс — коэффициент зависимости между веществами В и С.
Если в схеме реакций вещество А связано с веществом С как прямо, так и косвенно, рассчитывается обобщенный коэффициент зависимости между веществами А и С:
RAC = ™ХЛГАСгг } all path г
(3)
Каждое вещество связано с другими веществами с различными обобщенными коэффициентами зависимости. Любое вещество Х будет выбрано в качестве вещества, связанного с целевым веществом А, если
(4)
где е — является определенным пороговым значением (0<Е<1).
Результирующими веществами сокращенного механизма является объединение веществ из всех подмножеств каждого целевого вещества. Остальные вещества на данный момент являются избыточными по отношению к целевым веществам и могут быть безопасно удалены из списка продуктов химической реакции. Следовательно, все стадии, которые потребляют избыточные вещества, могут быть удалены.
(1) Вычислительный эксперимент
где — скорость г-й реакции;
УА{ — стехиометрический коэффициент вещества А в г-й реакции (стехиометрический коэффициент положительный, если А — продукт, и отрицательный, если А — реагент);
Зт =1, если вещество В участвует в г-й реакции, Зт =0, в противном случае;
N — количество реакций в системе.
Коэффициент зависимости между каждой парой веществ является элементом матрицы. Каждый ее элемент удовлетворяет соотношению 0 < гАВ < 1 .
Если вещество А связано с веществом С косвенно (реакция вида А^В, В^С), то коэффициент зависимости между веществами А и С рассчитывается по формуле:
Рассмотрим работу программного комплекса на примере сокращения механизма диме-ризации а-метилстирола в присутствии цеолита МаНТ 6. Схема реакции имеет вид:
2 Х1 о Х2, 2 Х1 о Х3, 2X ^ Х4, X2 о Х3, Х2 ^ Х4,
X3 —' X4'
Xj + X 2 —^ X 5 ,
X + X3 — X5,
Xj + X4 — X5,
i=1
где X — а-метилстирол;
X 2 — 4-метил — 2,4 — дифенилпептен-1 (а-димер);
Х3 — 4-метил — 2,4 —дифенилпептен-2 (¿8-димер);
X 4 — 1,1,3 — триметил-3-фенилиндан (циклический димер);
X 5 — тримеры.
1 этап. Построение кинетической модели. Схеме реакции (5) соответствуют уравнения скоростей стадий:
М>1 = к1П1 5
щ = к2П2,
Wз — кзП1,
Щ4 и к 3
Щ и
Щ = к6П2 5
Щ II к 3
Щ = к8п2 5
II к 9 3
Щ10 = к10п1п2,
Щ11 — кцППз,
Щ12 — к 12 П1 п 4 •
(6)
где щ — концентрация вещества Хг; константы скоростей ks(T), 5=1,...12, зависят от температуры Т, исходя из уравнения Аррениуса:
_ Ел
кг (Т) = к> *Т
(7)
С учетом матрицы стехиометрических коэффициентов дифференциальные уравнения, описывающие схему реакций, имеют вид:
йг йп2 йг йп3 йг йп4 йг йп5 йг
1 = _2ю1 + 2ю2 _ 2ю3 + 2Ю4 _ 2Ю5 _ а>10 _ ю11 _ сохг,
2 этап. Решение прямой задачи химической кинетики. Пользователем вводятся исходные данные: начальные концентрации веществ (щ, мольная доля), опорная температура (T0p, К), значение констант скоростей реакций при этой температуре ^((Т0р), м3/кг-ч, значения энергий активации (Ег, ккал/моль), температура, при которой протекает реакция (Т, К), время протекания реакции ч) и шаг расчета (^ ч). По введенным данным рассчитываются константы скоростей реакций
Е>Т т
кг — ког (Тор) • еКТ Тор ,
решается система дифференциальных уравнений (кинетическая модель схемы реакции) явным методом Рунге — Кутта четвертого порядка. Полученное решение представляется в численном и графическом виде.
3 этап. Сокращение механизма химической реакции на основе ОЯОЕР-ме-тода. Для наглядного представления связей между веществами строится ориентированный граф, отражающий зависимость веществ друг от друга (рис. 1). Вершины графа соответствуют веществам, участвующим в реакции; ребро, направленное от г-й вершины к j-й вершине, соответствует образованию вещества Х}- из X
В качестве целевых веществ были выбраны вещества Х2 и Х3. Расчет проводился с точностью е=0.001 в момент времени t=0.04 ч. Рассчитываются коэффициенты зависимости Гдв для всех веществ. С помощью алгоритма поиска графа в глубину определяются все пути, соединяющие каждое целевое вещество с остальными. По найденным путям рассчитываются обобщенные коэффициенты связи веществ (рис.2).
Рис. 1. Граф связей между веществами
Обобщенный коэффициент связи R:
х1 к2 хЗ к4 к5
к1 1 0,001380 0,000426 0 0
к2 0,949961 1 0,000405 0 0
хЗ 0,651590 0,204660 1 0 0
и4 0,959459 0,035307 0,005035 1 0
х5 1 0,203815 0,985872 0,000850 1
Рис.2. Таблица обобщенных коэффициентов связей веществ
На основе рассчитанных обобщенных коэффициентов связей и точности определяются исключаемые вещества и реакции (рис. 3).
ких и термодинамических моделей реакции (5) можно использовать схему реакции (9), которая является в каком-то смысле эквивалентной
0.8-
-X, -X.
- Х( х. 0.4-
0.6 -
—X»
0.00 0.01
0.02 0,03 0.04 0.05
t
Рис.3. Сокращенный механизм реакции
Результаты и обсуждение
Исходя из полученных значений обобщенных коэффициентов связей и точности расчета, из схемы реакции были исключены вещества Х4 и Х5. В соответствии с этим сокращенная схема реакции имеет вид:
2 Xj о X2, 2 Xj о X3, X2 оX3.
(9)
На рис. 4 представлена динамика концентраций веществ сокращенного механизма и концентраций этих же веществ в исходном механизме.
Как видно из рисунка, уточнение механизма реакции (5) не изменило общую динамику изменения концентраций целевых веществ во времени. Поэтому для анализа кинетичес-
Рис. 4. Изменение концентраций веществ (Xi — вещества исходного механизма, Х'— вещества сокращенного механизма, i=1,2,3)
Таким образом, программный комплекс позволяет уточнить механизм химической реакции при заданной точности и в заданный момент времени протекания реакции; решить прямую задачу для нежестких систем дифференциальных уравнений при заданных единицах измерения температуры, констант скоростей реакции, начальных концентраций, энергий активации. Полученное решение представляется в численном и графическом виде. Для работы с программой необходимо иметь персональный компьютер типа IBM PC Pentium IV с операционной системой Windows 2000 и выше и оперативной памятью от 512 Мб.
Литература
1. Nancy J. Brown, Guoping Li. // International Journal of Chemical Kinetics.— 1997.— V.29.— P.393.
2. Robert P. Dickinson and Robert J. Gelinas. — A Direct Method. // Journal of Computational Physics.- 1976.- V.21.— P.123.
3. Вольперт А. И., Худяев С. И. Анализ в классах разрывных функций и уравнения математической физики.— М: Наука, 1975.— 395 с.
4. T. Lu and C.K. Law. // Proceedings of the Combustion Institute.— V. 30.— 2005.— P. 1333.
5. A. G. Xia, D. V. Michelangeli, P. A. Makar. // Atmospheric Chemistry and Physics.— 2009.— V.9.— P. 4341.
6. Байтимерова А. И. Математическое моделирование и численное исследование каталитических процессов в каскаде реакторов: Дис. ... канд. физ.-мат. наук.— Уфа.— 2009.— 127 с.
