CC BY
9
НАЗАРОВ В. И., ПЕСКОВА Е. Е., ЯЗОВЦЕВА О. С.
ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ ЖЕСТКИХ СИСТЕМ С ИСПОЛЬЗОВАНИЕМ (4,2)-МЕТОДА
Аннотация. Описано применение (4,2)-метода к численному решению задачи Коши для нелинейной системы обыкновенных дифференциальных уравнений, соответствующей математической модели брутто-реакции пиролиза этана. Проведен сравнительный анализ (4,2)-метода и явного метода Эйлера, показана эффективность (4,2)-метода.
Ключевые слова: химическая кинетика, жесткие системы дифференциальных уравнений, явный метод Эйлера, пиролиз этана.
NAZAROV V. I., PESKOVA E. E., YAZOVTSEVA O. S.
NUMERICAL MODELING OF HARD SYSTEMS BY (4,2)-METHOD
Abstract. The article deals with the use of the (4,2)-method for a numerical solution of the Cauchy problem for the nonlinear system of ordinary differential equations, which corresponds to the mathematical model of the ethane pyrolysis brutto-reaction. A comparative analysis of the (4,2)-method and the explicit Euler method is carried out. The efficiency of the (4,2)-method is shown.
Keyword: chemical kinetics, hard systems of differential equations, explicit Euler method, ethane pyrolysis.
Моделирование процесса термического пиролиза этана для получения этилена является актуальной задачей. При решении таких задач требуется найти решение задачи Коши для систем обыкновенных дифференциальных уравнений, составленных на основе схемы реакций. Большая часть этих систем является жесткой по части фазовых переменных. Выполнение свойств жесткой системы влечет за собой большое различие в величине коэффициентов, что, в свою очередь, не позволяет применить к решению задачи Коши стандартные явные методы. В работах [1; 2] подробно описано построение (ш, &)-методов для решения жестких систем дифференциальных уравнений, основанных на методах типа Розенброка. В данной работе описан (4,2)-метод, примененный к решению задачи Коши для нелинейной системы обыкновенных дифференциальных уравнений, соответствующей математической модели брутто-реакции пиролиза этана.
Рассмотрим брутто-реакцию пиролиза этана [3; 4]:
2С2^6 ^ ^2^4 + 2СН4.
Соответствующие кинетические уравнения имеют вид
[С2Н6] ^ [С2Н4] + [Н2], w± = klCl,
2[С2Н6] ^ [С2Н4] + 2[СН4], W2 = к2с2,
где Ci — концентрация -того вещества, kj — константа скорости у'-той прямой стадии, которые вычисляются согласно выражениям Аррениуса при температуре 800 К, таким образом к1 = 0.051, к2 = 0.7.
Составим систему обыкновенных дифференциальных уравнений, в которой фазовые переменные соответствуют концентрациям веществ, участвующих в реакции:
с1 = — к1с1 — 2к2с1
с2 = к1с1 + k2cl ^
¿3 = к1с1
-С4 = 2к2с1
Поскольку начальный момент реакции присутствует только этан, то начальные условия будут иметь вид:
Ci(0) = 1, С2(0) = 0, Сз(0) = 0, С4(0) = 0. (2)
Для нахождения изменений количественных характеристик реакции численно решим задачу Коши (1)-(2) (т, к)-методом и проведем сравнительную характеристику со способом решения методом Эйлера.
Используя (4,2)-метод четвертого порядка точности, разностную схему для каждого компонента вектора решения построенной модели можно представить в виде [1; 2]:
4
,п+1
= УП + X Prti ,°п = Е — ™fn,
У
1=1
ОпЧ1 = ^(Уп), ОпЧ2 = 41' ОпЧз = ^(уп + Р31Ч1 + (З32Ч2) + а3242' ОпЧ4 = 43+ а4242'
где т - шаг интегрирования решаемой задачи, Е — единичная матрица размерности (4x4),
у ей4, Ч1ек4'Уо = (1'0'0'0).
Вещественные коэффициенты, при которых схема обладает устойчивостью и 4-м порядком точности, а, р^ а^2, р3к, 1 < / < 4,3 < } < 4,1 < q <2 , имеют следующие значения [1; 2]:
а = 0.57281606248213, р{ = 1.27836939012447, р2 = —1.00738680980438, р3 = 0.9265539109350, р4 = —0.33396131834691, @31 = 1.00900469029922, Р32 = —0.25900469029921, а32 = —0.49552206416578, а42 = —1.28777648233922.
Для системы (1) /(у) и fn имеют вид
(- - 2^2УЛ
/(У) =
V
&1У1 + ^2У1 ^1У1 2^2У1
/
/п =
5/(у) Зу
/— - 4^2У1 + 2^У1 ^1У1 4^2У1
0 0 0 0
0 0 0 0
0 0 0 0
V
На основе построенной схемы был разработан программный комплекс для решения системы уравнений (1) с начальными данными (2). На рисунке представлены изменения концентраций веществ во времени при температуре 800 К. Для оценки эффективности (4,2)-метода эта же система была решена явным методом Эйлера, разностная схема для которого имеет вид [5]:
Уп+1=Уп + т/(у").
Сравнительные результаты испытаний представлены в таблице 1, за опорное было принято решение, полученное с помощью системы Wolfram Mathematica с шагом т = 0.1.
Таблица 1
Результаты работы методов
Число итераций Время выполнения Точность, £
Явный метод Эйлера (4,2) — метод Явный метод Эйлера (4,2) — метод
10 6 мс 7 мс 0.034 0.001
50 59 мс 49 мс 0.00762 0.00055
93 105 мс 98 мс 0.0043 0.0001
Из таблицы 1 видно, что многостадийный метод имеет преимущество по времени выполнения, точности и корректности по отношению к явному методу.
Рис. 1. Изменение концентраций веществ при 800 К.
Таким образом, сравнительный анализ показал эффективность (т, ^)-метода для решения жестких систем обыкновенных дифференциальных уравнений, в частности, для нелинейной системы обыкновенных дифференциальных уравнений, соответствующей математической модели брутто-реакции пиролиза этана, поскольку заданная точность достигается за меньшее количество итераций по сравнению с методом Эйлера.
ЛИТЕРАТУРА
1. Новиков Е. А. Исследование (т,2)-методов решения жестких систем // Вычислительные технологии. - 2007. - Т. 12.- № 5. - С. 103-115.
2. Галанин М. П., Ходжаева С. Р. Методы решения жестких обыкновенных дифференциальных уравнений. Результаты тестовых расчетов // Препринты ИПМ им. М. В. Келдыша. - 2013. - № 98. - 29 с.
3. Губайдуллин И. М., Пескова Е. Е., Язовцева О. С. Математическая модель динамики многокомпонентного газа на примере брутто-реакции пиролиза этана [Электронный ресурс] // Огарев^Ы^. - 2016. - №20. - Режим доступа: http://joumal.mrsu.ru/arts/matematicheskaya-model-dinamiki-mnogokomponentnogo-gaza-na-primere-brutto-reakcii-piroliza-etana.
4. Мухина Т. Н., Барабанов Н. Л., Бабаш С. Е., Меньщиков В. А., Аврех Г. Л. Пиролиз углеводородного сырья. - М.: Химия, 1987. - 240 с.
5. Самарский А. А., Гулин А. В. Численные методы: учеб. пособие для вузов. - М.: Наука, 1989. - 432 с.
