DOI: 10.18454/IRJ.2016.51.161 Икрамов Р. Д.1, Мустафина С.А.2
1ORCID: 0000-0003-0325-7859, Аспирант, 2Доктор физико-математических наук, Стерлитамакский филиал Башкирского государственного университета ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ДИНАМИКИ РАСШИРЕННОЙ МОДЕЛИ РЕАКЦИИ
БЕЛОУСОВА-ЖАБОТИНСКОГО
Аннотация
В статье рассмотрена одна из моделей колебательной реакции Белоусова-Жаботинского - 6-ти стадийная модель Орегонатора. При проведении численного исследования использован L-устойчивый метод на основе АВС-схемы с действительными коэффициентами, подходящий для решения систем дифференциальных уравнений с высоким коэффициентом жесткости. Для модели Орегонатора проведен анализ колебаний, найдены периоды колебаний. Реакция Белоусова-Жаботинского является ярким примером колебательных химических реакций.
Ключевые слова: реакция Белоусова-Жаботинского, Орегонатор, колебания, осцилляции.
Ikramov R.D.1, Mustafina S.A.2
1ORCID: 0000-0003-0325-7859, Postgraduate student, 2PhD in Physics and Mathematics,
Sterlitamak Branch of Bashkir State University NUMERICAL INVESTIGATION OF DYNAMICS OF THE EXTENDED MODEL OF THE BELOUSOV-ZHABOTINSKI REACTION
Abstract
6-stage Oregonator model, one of the models of the oscillating Belousov-Zhabotinski reaction, was considered in the article. L-stable method based on ABC-scheme with real coefficients is used for solving differential equations system. Method is appropriate for systems with high stiffness coefficient. Article considered analysis of oscillations, values of the periods of oscillations. Belousov-Zhabotinski's reacton is the one of the greatest chemical oscillating reaction.
Keywords: Belousov-Zhabotinski's reacton, Oregonator, oscillations.
Большинство химических превращений протекает через огромное количество элементарных стадий с большим числом различных интермедиантов, которые могут включать атомы и радикалы. При этом их концентрации не всегда описываются возрастающими или убывающими кривыми. При этом стоит учитывать трудность экспериментальной регистрации малых концентраций промежуточных веществ, а большое количество перекрестных взаимодействий между веществами лишь усложняет исследование реакции.
Колебательными называются реакции, в которых происходит периодическое изменений концентраций веществ, участвующих в самой реакции.
Колебательные процессы играют важную роль не только в химии, но также и биологии, биохимии и в других отраслях наук. Моделирование колебательных процессов позволяет понять, как протекают те или иные явления, происходящие в организме человека (колебания сердечной мышцы, биоритмы), в каталитических реакциях и т.д. Таким образом, исследование колебательных процессов и реакций является на сегодняшний день актуальной задачей.
Основная трудность при решении задач, описывающих колебательные процессы, заключается в необходимости обеспечения правильных значений амплитуд и фаз на протяжении многих периодов. Основная причина этому заключается в большом значении коэффициента жесткости и высокой размерности системы дифференциальных уравнений (СДУ) математических моделей. Поскольку характер поведения и тип СДУ может меняться в ходе интегрирования, то подбор численного метода интегрирования является важным этапом в процессе исследования колебательных процессов.
В начале XXI-го века в связи с бурным развитием компьютерных и информационных технологий вновь возник интерес к нестационарным и нелинейным явлениям в химической кинетике, называемых критическими [1]. Они характеризуются, в основном, множественностью стационарных состояний и автоколебаниями. Автоколебательными системами, в свою очередь, называются такие системы, проявляющиеся, во-первых, в незатухающих колебаниях вне зависимости от начальных условий, во-вторых, в их устойчивости по отношению к отклонениям от стационарного режима. Таким образом, в автоколебательной системе устанавливаются и поддерживаются незатухающие колебания за счет сил, зависящих от состояния системы [2].
Наиболее известной автоколебательной гомогенной реакцией считается реакция Белоусова-Жаботинского, в гетерогенных реакциях автоколебания были найдены Г.К. Боресковым и М.Г. Слинько.
В 1951 году Б.П. Белоусов открыл и экспериментально исследовал химическую реакцию окисления лимонной кислоты броматом, катализируемой ионами металлов в кислой среде. Возникающий в реакции колебательный режим с периодической сменой окраски стал классическим примером колебательных реакций [4]. Таким образом, Б.П. Белоусов внес основной вклад в открытие колебательных химических реакций. Позднее в 1955 г. И.Р. Пригожин доказал возможность существования в открытой системе химических колебаний около стационарного состояния, достаточно удаленного от термодинамического равновесия.
Механизм реакции Белоусова-Жаботинского насчитывает более 80 стадий, очевидно, что исследование закономерностей, а также решение прямой, обратной, оптимизационной задач не является возможным. Ученые из штата Орегон Р. Филд, Е. Кёрёш, Р. Нойес предложили механизм реакции Белоусова-Жаботинского, состоящий из 11 реакций между 15 различными соединениями, причем концентрации соединений достаточно сильно изменяются в ходе колебательного цикла [5]. Поэтому кинетика реакции может быть описана системой нелинейных обыкновенных дифференциальных уравнений первого порядка, полученных применением закона действующих масс к каждой стадии, чаще всего решаемых численно. Для исследования колебаний, бистабильности и бегущих волн аналитически необходимо сведение полной модели к более простой. Р. Филд и Р. Нойес первыми показали, что механизм Филда -
Кереша-Нойеса может быть упрощен до поддающейся анализу математической модели. Поэтому авторами работы [6] была предложена абстрактная и простая модель реакции Белоусова-Жаботинского под названием «Орегонатор», которую впоследствии расширили [7] до шести стадий:
А + У ^Х + Р, Х + У <г-> 2 Р, А + Х <г-> 2И/, С + Ш ^ X + г, 2Х <-» А + Р, г С,
где А = Вг О- (бромат-ион), С = М (п) - ион металла катализатора, Р = Н ОВ г (бромноватистая кислота), Ш = ВгО2 (диоксид брома), X = НВгО2 (бромистая кислота), У = Вг~ (бромид-ион), 1 = М(п + 1) - окисленная форма иона металла катализатора, д - стехиометрический множитель.
Математическая модель Орегонатора, дающего сложный предельный цикл, имеет вид [8]:
С1 = -171 - + с2 = — уг — р2 + 0.462г?6 с3 = -г?4 + г?6
с'4 = Р1-Р2-Р3+Р4- 2 (1)
с5 = у±+2у2 + г?5 с6 = 2г?3 - г?4 Су = р4- г?6
где задаются следующими соотношениями:
р1 = /Сх^Сг — /с_1с4с5,
= ^2С2С4 — к-2С5С5'
= к3с±с4 — к_3с6с6, 1^4 — к4с3с6 к_4с4с7,
= ^5С4С4 — к-5с4с5'
У6 = к6с7.
Концентрации реагентов обозначены следующим образом: с1 = [Вг О-] , с2 = [В г~] , с3 = [М (п) ] , с4 = [НВг О2 ] , с^ = [Н ОВ г] , с^ = [Вг О2] , с7 = [М (п + 1) ] .
Константы скорости принимают следующие значения (М 1с -1):
/сх = 0.084,/е.! = 104, к2 = 4 ■ 108, к_2 = 5 ■ 10"5, к3 = 2 ■ 103, к_3 = 2 ■ 107, к4 = 1.3 ■ 105, к_4 = 2.4-107, к5 = 4 ■ 104, /с_5 = 4 ■ Ю-11, к6 = 0.65 с"1.
Основная трудность в моделировании колебательных химических реакциях заключается в разработке численных алгоритмов, достаточно точных и устойчивых для получения адекватных и качественно точных результатов [9]. Так как константы скорости различаются между собой на многие порядки, то коэффициент жесткости системы (1) принимает большое значение [10]. В силу этого, использование явных численных схем не представляется возможным. Для численного исследования модели (1) использовалась явно-неявная АВС-схема с действительными коэффициентами [11]:
[Е-к-Гу+^-Гу2 ][с (Ь + К)-с (0] = [Е-Н--Гу]к-Г, (2)
где с ( Ь) - численное решение системы (1) на шаге интегрирования в момент времени £ с шагом й, [у - матрица Якоби системы дифференциальных уравнений (1), Е - единичная матрица. Данная схема является Ь-устойчивой и обладает вторым порядком точности.
Результаты интегрирования системы (1) при начальных условиях (М) : с1 = 0.1 387, с2 = 1.5 34 ■ 10 - 7, с3 = 1. 1 7 6 ■ 10-4, с4 = 3. 165 ■ 10-8, с5 = 1.9 5 6 ■ 10-4, с6 = 5 . 814 ■ 10-7, с7 = 6. 3 1 ■ 10-6 представлены на рис.1 (а-ж). Шаг интегрирования .
д)
О 100 200 300 400 500 t, С
е)
О 100 200 300 400 500 t, С
ж)
Рис. 1 - Колебания значений концентрации реагентов в зависимости от времени в модели (1): а) [ВгОз], б) [Вг~], в) [М(п)], г) [НВг02], д) [НОВг], е) [Вг02], ж) [М(п + 1)].
По рис. 1 видно, что модель Орегонатора (1) характеризуется длительными постепенно затухающими колебаниями. Реагент [Вг0^] характеризуется постепенным периодическим уменьшением своей концентрации, поскольку является исходным реагентом. Реагент характеризуется постепенным увеличением своей
концентрации, накапливаясь как продукт реакции. Остальные компоненты системы характеризуются релаксационными режимами колебаний с быстрыми скачками. Период скачков составляет Т = 6 3 . 5 сек и постепенно увеличивается [12]. Суммарная концентрация катализатора и его иона остается постоянной, что подтверждает адекватность данной модели.
а) б)
Рис. 2 - Фазовые портреты модели (1): а) зависимость концентрации [М (п) ] от [Вг ] , б) [М (n + 1 ) ] от [М (п) ] .
На рис. 2 представлены фазовые портреты в логарифмической плоскости, представляющие собой зависимости значений концентраций одного реагента от другого. На рис. 2а представлена зависимости концентрации металла катализатора от концентрации бромид-иона. Фазовый портрет представляет собой предельный цикл, замкнутые круговые области, что говорит, что в данной системе наблюдаются автоколебания. На рис. 2б представлен предельный цикл, вырожденный в прямую линию, что говорит о логарифмической зависимости между концентрациями катализатора и его иона.
Литература
1. Быков В. И. Моделирование критических явлений в химической кинетике. М.: Наука, 1988. 345 с.
2. Рубин А. Б., Пытьева Н. Ф., Ризниченко Г.Ю. Кинетика биологических процессов. Учебное пособие. М.: Изд-во Моск. ун-та, 1977. 330 с.
3. Lotka A. J. Elements of Physical Biology. Baltimore: Williams & Wilkins Company, 1925. 460 p.
4. Жаботинский А.М. Концентрационные колебания. М.: Наука, 1974. 179 с.
5. Жаботинский А.М., Отмер Х., Филд Р. Колебания и бегущие волны в химических системах. М.: Мир, 1988. 720 с.
6. Field R.J., Koros E., Noyes R.M. Oscillations in chemical systems. IV. Limit cycle Behavior in a model of real chemical reaction // J. Chem. Phys. - 1974. - Vol. 60. - P. 1877-1884.
7. Showalter K., Noyes R.M., Ber-eli K. A modified Oregonator exhibiting complicated limit cycle behavior in a flow system // J. Chem. Phys. - 1978. - Vol. 69. - P. 2514-2524.
8. Икрамов Р.Д., Мустафина С.А. Численное исследование моделей реакции Белоусова-Жаботинского на основе двухстадийного метода Розенброка с комплексными коэффициентами // Системы управления и информационные технологии. - 2014. - №2. - С. 11-14.
9. Икрамов Р.Д., Мустафина С.А. Численное исследование колебательных реакций с помощью метода Розенброка с действительными коэффициентами // Журнал Средневолжского математического общества. - 2014. - Т.16. - №1. -С.71-75.
10. Икрамов Р.Д., Мустафина С.А. Численное исследование моделей Орегонатора с использованием двухстадийного метода Розенброка с комплексными коэффициентами // Информационные технологии моделирования и управления. - 2014. - №6. - С. 211-217.
11. Булатов М.В., Тыглиян А.В., Филиппов С.С. Об одном классе одношаговых одностадийных методов для жестких систем обыкновенных дифференциальных уравнений // Журнал вычислительной математики и математической физики. - 2011. - Т.51. - №7. - С.1251-1265.
12. Икрамов Р. Д. Моделирование и численное исследование динамики колебательных химических реакций полунеявными методами: Автореф. дис. канд. физ. -мат. каук. - Уфа, 2016. - 24 с.
References
1. Bykov V. I. Modelirovanie kriticheskih javlenij v himicheskoj kinetike. M.: Nauka, 1988. 345 s.
2. Rubin A. B., Pyt'eva N. F., Riznichenko G.Ju. Kinetika biologicheskih processov. Uchebnoe posobie. M.: Izd-vo Mosk. un-ta, 1977. 330 s.
3. Lotka A. J. Elements of Physical Biology. Baltimore: Williams & Wilkins Company, 1925. 460 p.
4. Zhabotinskij A.M. Koncentracionnye kolebanija. M.: Nauka, 1974. 179 s.
5. Zhabotinskij A.M., Otmer H., Fild R. Kolebanija i begushhie volny v himicheskih sistemah. M.: Mir, 1988. 720 s.
6. Field R.J., Koros E., Noyes R.M. Oscillations in chemical systems. IV. Limit cycle Behavior in a model of real chemical reaction // J. Chem. Phys. - 1974. - Vol. 60. - P. 1877-1884.
7. Showalter K., Noyes R.M., Ber-eli K. A modified Oregonator exhibiting complicated limit cycle behavior in a flow system // J. Chem. Phys. - 1978. - Vol. 69. - P. 2514-2524.
8. Ikramov R.D., Mustafina S.A. Chislennoe issledovanie modelej reakcii Belousova-Zhabotinskogo na osnove dvuhstadijnogo metoda Rozenbroka s kompleksnymi kojefficientami // Sistemy upravlenija i informacionnye tehnologii. -2014. - №2. - S. 11-14.
9. Ikramov R.D., Mustafina S.A. Chislennoe issledovanie kolebatel'nyh reakcij s pomoshh'ju metoda Rozenbroka s dejstvitel'nymi kojefficientami // Zhurnal Srednevolzhskogo matematicheskogo obshhestva. - 2014. - T.16. - №1. - S.71-75.
10. Ikramov R.D., Mustafina S.A. Chislennoe issledovanie modelej Oregonatora s ispol'zovaniem dvuhstadijnogo metoda Rozenbroka s kompleksnymi kojefficientami // Informacionnye tehnologii modelirovanija i upravlenija. - 2014. - №6. - S. 211-217.
11. Bulatov M.V., Tyglijan A.V., Filippov S.S. Ob odnom klasse odnoshagovyh odnostadijnyh metodov dlja zhestkih sistem obyknovennyh differencial'nyh uravnenij // Zhurnal vychislitel'noj matematiki i matematicheskoj fiziki. - 2011. -T.51. - №7. - S.1251-1265.
12. Ikramov R. D. Modelirovanie i chislennoe issledovanie dinamiki kolebatel'nyh himicheskih reakcij polunejavnymi metodami: Avtoref. dis. kand. fiz.-mat. kauk. - Ufa, 2016. - 24 s.
DOI: 10.18454/IRJ.2016.51.139 Магомедов А.М.
Доктор физико-математических наук, профессор, Дагестанский государственный университет Работа выполнена при поддержке проекта № 2588 в рамках базовой части государственного задания Минобрнауки РФ и Отдела математики и информатики ДНЦ РАН О РАСКРАШИВАЕМОСТИ ДВУДОЛЬНЫХ ГРАФОВ СПЕЦИАЛЬНОГО ВИДА
Аннотация
A. S. Asratyan и C. J. Casselgren доказали, что задача об интервальной реберной раскраске бирегулярного двудольного графа со степенями 6 и 3 соответственно («(6,3)-граф») NP-полна. В статье показано, что минимальная мощность n вершин, при которой (6,3)-граф не допускает раскраски требуемого вида, равна 18: любой (6,3)-граф с n< 18 допускает интервальную реберную раскраску; для каждого n, не меньшего, чем 18 (и кратного 3), существует (6,3)-граф с n вершинами, для которого такая раскраска невозможна. Результаты находят применение в вопросах построения оптимальных расписаний учебных занятий.
Ключевые слова: граф, двудольный, раскраска.
Magomedov A.M.
PhD in Physics and Mathematics, Professor, Dagestan State University THE COLORING OF SPECIAL KIND BIPARTITE GRAPHS
Abstract
A. S. Asratyan and C. J. Casselgren proved that the problem of the interval edge-coloring of biregular bipartite graph with degrees 6 and 3 respectively («(6,3)-graph») is NP-complete. The article shows that the minimum capacity of vertices n such that (6,3)-graph prevents coloring of the required form equal to 18: any (6,3)-graph with n< 18 allows interval edge-coloring; for each n not less than 18 (and divided by 3) there is (6,3)-graph with n vertices, for which a coloring impossible. The results are used for constructing an optimal schedules training sessions.
Keywords: graph, bipartite, coloring.
1. Введение
В статье использованы обозначения и определения из книги [1]. Интервальной реберной раскраской графа t цветами будем называть отображение множества ребер графа в множество { 1 , 2 ,. . ,,t} , такое, что: 1) для каждого , , найдется ребро, раскрашенное в цвет , 2) в каждой вершине графа все представленные в ней
цвета попарно различны и образуют целочисленный интервал.
Двудольный граф в котором степени всех вершин равны 6, а степени всех вершин равны 3,
будем называть (6,3)-бирегулярным или, кратко, (6,3)-графом. При | X| = п (6,3)-граф G = (X, У, Е) будем называть (6, 3 ) п-графом. Будем называть (6,3)-граф раскрашиваемым или нераскрашиваемым в зависимости от того, допускает ли он интервальную реберную раскраску шестью цветами или нет.
В [2] показана NP-полнота задачи о раскрашиваемости (6,3)-графа. Пример нераскрашиваемого (6 , 3 ) 2 „-графа был построен еще в 1991 г. [3], в 2015 г. пример нераскрашиваемого (6, 3 ) 7-графа опубликован в [4].
Сформулируем основной результат статьи.
Теорема 1. При п < 5 каждый (6 , 3 ) п-граф раскрашиваем, для любого п > 6 найдется нераскрашиваемый (6, 3 ) п-граф.
2. Нераскрашиваемые (6,3 ) „ -графы
Отображение множества ребер Е (6,3)-графа G = (X, У, Е) в множество из двух цветов, такое, что в каждой вершине цвета всех трех ребер, инцидентных , одинаковы, а каждой вершине инцидентны три ребра каждого из