Авторефераты по всем темам  >>  Авторефераты по техническим специальностям

На правах рукописи

Лашина Елена Александровна

Анализ структурно устойчивых периодических решений кинетических моделей каталитических реакций

05.13.18 - математическое моделирование, численные методы и комплексы программ

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Новосибирск - 2012

Работа выполнена в Федеральном государственном бюджетном учреждении науки Институте катализа им. Г.К. Борескова Сибирского отделения Российской академии наук

Научный консультант: кандидат технических наук, Чумакова Наталия Алексеевна

Официальные оппоненты: Лаевский Юрий Миронович, доктор физико-математических наук, профессор, Институт вычислительной математики и математической геофизики Сибирского отделения Российской академии наук, зав. лабораторией Волокитин Евгений Павлович, кандидат физико-математических наук, доцент, Институт математики им. С.Л. Соболева Сибирского отделения Российской академии наук, с.н.с.

Ведущая организация: Московский государственный университет им. М.В. Ломоносова, факультет вычислительной математики и кибернетики

Защита состоится 20 марта 2012 года в 16.30 на заседании диссертационного совета Д 003.061.02 при Учреждении Российской академии наук Институте вычислительной математики и математической геофизики Сибирского отделения РАН по адресу: 630090, г. Новосибирск, проспект Академика Лаврентьева, 6.

С диссертацией можно ознакомиться в библиотеке Федерального государственного бюджетного учреждения науки Института вычислительной математики и математической геофизики Сибирского отделения Российской академии наук.

Автореферат разослан 17 февраля 2012 года.

Ученый секретарь диссертационного совета Д 003.061.02, доктор физико-математических наук Сорокин С.Б.

Общая характеристика работы

Актуальность. При экспериментальном исследовании кинетики и динамики гетерогенных каталитических реакций обнаруживаются автоколебания и нерегулярная динамика скорости реакции. Для теоретического описания таких явлений рассматривают кинетические модели, являющиеся системами нелинейных обыкновенных дифференциальных уравнений (ОДУ). Математическим образом автоколебаний служат грубые (структурно устойчивые) устойчивые периодические решения. Структурно устойчивое периодическое решение системы нелинейных ОДУ, содержащей параметр, можно продолжить по параметру на максимально возможный интервал. Построенное таким образом семейство периодических решений является максимальным. Для значения параметра на границе этого интервала в системе происходит бифуркация периодического решения.

Если стадии реакции различаются по масштабу скоростей и автоколебания имеют релаксационный характер, то в соответствующей модели возможно выделение быстрых и медленных переменных. В работе [1] предложен принцип генерирования нерегулярной динамики в системе трех нелинейных ОДУ с одной медленной переменной, вида = f(x, y, z), = g(x, y, z), = h(x, y, z), где 0 < 1 - малый параметр. Рассматривается случай, когда вырожденная система = f(x, y, z), = g(x, y, z) с параметром z имеет гистерезис на кривой стационарных состояний, и существуют два максимальных семейства грубых устойчивых периодических решений для z (a1, b1) и z (a2, b2), b1 < a2. При z = a1 и z = b2 в вырожденной системе происходит бифуркация Андронова-Хопфа, при z = b1 и z = a2 периодические решения вырождаются в петли сепаратрис седловых особых точек. Нерегулярная динамика системы с тремя переменными связана с возможностью изображающей точки в фазовом пространстве последовательно возвращаться в окрестность гомоклинической траектории вырожденной системы, где решения обладают высокой параметрической чувствительностью.

В большинстве случаев построение решения системы ОДУ проводят с применением численных методов. Если решение обладает высокой параметрической чувствительностью, то важным является анализ глобальной ошибки численного интегрирования.

Таким образом, актуальными являются разработка и исследование кинетических моделей гетерогенных каталитических реакций, описывающих автоколебания и нерегулярную динамику, с учетом выделения в системе быстрых, умеренных и медленных переменных.

Цель и задачи исследования. Целью работы является развитие качественных и численных методов анализа максимальных семейств структурно устойчивых периодических решений, а также сложной и нерегулярной динамики в кинетических моделях гетерогенных каталитических реакций.

Для достижения поставленной цели рассмотрены следующие задачи:

1. Разработка алгоритма и комплекса программ для уточнения петли сепаратрисы седла в системе двух нелинейных ОДУ, использующего понятие проектора на инвариантное многообразие седлового стационарного состояния.

2. Построение и параметрический анализ простейшей кинетической модели, описывающей релаксационные автоколебания и нерегулярную динамику скорости реакции окисления оксида углерода на металлах платиновой группы и учитывающей резкое изменение адсорбционной и реакционной способности катализатора в условиях реакции.

3. Анализ влияния параметра k2, линейно зависящего от парциального давления кислорода в газовой фазе, на структуру максимальных семейств грубых периодических решений вырожденной системы с параметром z в кинетической модели реакции окисления водорода. Построение оценки глобальной погрешности численного интегрирования уток-циклов.

4. Исследование влияния малого параметра на динамику кинетической модели реакции окисления водорода на никеле, являющейся системой трех нелинейных ОДУ с одной медленной переменной.

Научная новизна работы 1. Разработаны алгоритм и комплекс программ для уточнения петли сепаратрисы седла системы двух нелинейных ОДУ, использующие понятие проектора на инвариантное подпространство седлового стационарного состояния.

Доказано, что оператор, определяющий линейную часть возмущения проектора, является вещественным.

2. Предложена кинетическая модель реакции окисления оксида углерода на металлах платиновой группы, которая является системой трех нелинейных ОДУ и учитывает ступенчатые зависимости параметров k2 = k2(z) и E3 = E3(y), характеризующих асорбционные свойства и реакционную способность катализатора, от переменных y и z. Определены достаточные условия, при которых в вырожденной системе двух ОДУ с параметром z существует устойчивый предельный цикл. В численном эксперименте определены значения параметра (характеризующего скорость изменения медленной переменной z), при которых происходят бифуркации удвоения периода, а также значение , при котором динамика системы является нерегулярной.

3. Выполнен параметрический анализ вырожденной системы в кинетической модели гетерогенной каталитической реакции окисления водорода.

Указаны значения параметров z и k2, при которых существуют грубые устойчивое и неустойчивое периодические решения-утки. Построена оценка v(t) - главного члена асимтотического разложения глобальной погрешности численного интегрирования периодических решений-уток, на временном интервале, равном нескольким периодам, и показана зависимость динамики ||v(t)|| от мультипликаторов утки-цикла.

4. С применением численных методов в кинетической модели реакции окисления водорода получена нерегулярная динамика, которая зарождается в результате бифуркации удвоения периода при увеличении малого параметра .

Основные результаты, выносимые на защиту:

1. Алгоритм уточнения петли сепаратрисы седла системы двух нелинейных ОДУ и его реализация в виде комплекса программ. Доказательство вещественности оператора, определяющего линейную часть возмущения проектора на локальное устойчивое инвариантное подпространство седлового стационарного состояния.

2. Кинетическая модель реакции окисления оксида углерода, которая учитывает ступенчатые зависимости параметров E3(y) и k2(z), где y и z - переменные модели. Доказательство достаточности такой зависимости E3(y) для существования в вырожденной системе двух нелинейных ОДУ устойчивого предельного цикла.

3. Результаты анализа семейств структурно устойчивых периодических решений кинетической модели реакции окисления водорода с параметрами k2 и z. Оценка главного члена асимтотического разложения глобальной погрешности численного интегрирования периодических решений-уток на временном интервале, равном нескольким периодам.

4. Результаты исследования нерегулярной динамики кинетических моделей гетерогенных каталитических реакций окисления водорода и оксида углерода.

Практическая ценность. Показано, что предложенные зависимости констант скоростей отдельных стадий от концентраций промежуточных веществ являются достаточными для существования релаксационных автоколебаний и нерегулярной динамики в кинетических моделях гетерогенных каталитических реакций. С помощью разработанного алгоритма уточнения петли сепаратрисы седла определены значения параметра, при которых автоколебания с большими периодом и амплитудой исчезают.

Достоверность результатов. Исследование динамики кинетических моделей проводилось с применением методов качественной теории систем ОДУ, а также методов вычислительной математики. Построение оценки глобальной погрешности позволило обосновать точность результатов численного интегрирования.

Апробация работы. Основные результаты работы были представлены на следующих конференциях: международная школа-конференция молодых ученых по катализу Каталитический дизайн - от исследований на молекулярном уровне к практической реализации, Новосибирск, 2002; международная конференция молодых ученых по математическому моделированию и информатике, Новосибирск, 2002; международная конференция по вычислительной математике, Новосибирск, 2004; международная конференция по химическим реакторам Химреактор-18, Мальта, 2008; международная конференция Дифференциальные уравнения. Функциональные пространства. Теория приближений, Новосибирск, 2008; всероссийская конференция Математика в приложениях, Новосибирск, 2009; международная школасеминар Нелинейный анализ и экстремальные задачи, Иркутск, 2010; на семинаре ИМ СО РАН Избранные вопросы математического анализа под руководством д.ф.-м.н. Г.В. Демиденко; на Объединенном семинаре кафедры вычислительной математики НГУ и ИВМиМГ СО РАН, руководитель д.ф.-м.н. В.П. Ильин; общеинститутском семинаре ИК СО РАН Технология каталитических процессов.

ичный вклад соискателя. Автор принимал участие в обсуждении постановок задач, построении кинетической модели гетерогенной каталитической реакции окисления оксида углерода на металлах платиновой группы.

Автором выполнена реализация алгоритма уточнения петли сепаратрисы седла для системы двух нелинейных ОДУ в виде комплекса программ.

Доказательства выносимых на защиту теорем и результаты исследования структурно устойчивых периодических решений и нерегулярной динамики в рассматриваемых кинетических моделях получены автором.

Публикации. По теме диссертации опубликовано 9 печатных работ, из них 3 в ведущих рецензируемых журналах и изданиях, перечень которых определен Высшей аттестационной комиссией.

Структура и объем диссертации. Диссертационная работа состоит из введения, трех глав, заключения и списка литературы. Работа содержит 34 рисунка, 1 таблицу; список литературы состоит из 77 источников. Общий объем работы составляет 153 страницы.

Краткое содержание работы Во введении обоснована актуальность темы диссертации, сформулированы цель и задачи работы, показана научная новизна.

Первая глава содержит краткий обзор существующих кинетических моделей гетерогенных каталитических реакций окисления оксида углерода и водорода, являющихся системами нелинейных ОДУ. На основании обзора теоретических и экспериментальных данных делается вывод о необходимости создания простейшей математической модели, которая учитывает изменение свойств поверхности катализатора в ходе реакции и описывает регулярные и нерегулярные автоколебания скорости реакции.

В работе предлагается кинетическая модель реакции окисления оксида углерода на металлах платиновой группы. Модель описывает изменения безразмерных концентраций оксида углерода (x) и кислорода (y), адсорбированных на поверхности катализатора, а также концентрации кислорода, внедренного в приповерхностный слой металла (z):

= k1(1 - x - y) - k-1x - k3xy - xz, = 2k2(1 - x - y)2 - k3xy - y(1 - z), (1) = (y(1 - z) - xz), где параметр 0 < 1 является малым, и переменная z - медленная. Параметры ki, i = 1, 2, 3, и положительны.

Предполагается, что существуют критические значения y, при которых происходит изменение реакционной способности поверхности катализатора так, что параметр k3 = k3(y) является функцией от y, имеющей непрерывные первые производные:

k3(y) k30 exp(-E3(y)/(RT )), (2) где T [K] - температура катализатора, R = 1.987 [кал/(мольК)] - универсальная газовая постоянная, E3 = E3(y) [кал/моль] - энергия активации взаимодействия веществ, адсорбированных на поверхности металла. Функция E3(y) > 0 определена для y [0, 1], является ступенчатой и зависит от параметров yc и , так что E3(y) E31 для 0 y yc - , E3(y) E32 для dE3(y) yc+ y 1 и > 0 при yc- < y < yc+, E31 и E32 - положительные dy постоянные, E31 < E32.

Кроме того, предполагается, что при внедрении кислорода в приповерхностный слой металла происходит реконструкция поверхности катализатора и изменяются ее адсорбционные свойства. В модели это отражено в виде зависимости k2 = k2(z), где k2(z) C1 и имеет вид:

k2(z) k20(z)PO2, PO2 [Торр] - парциальное давление кислорода в газовой фазе. Функция k20(z) является ступенчатой и зависит от параметров zc и z, так что k20(z) kdk20(z) для 0 z zc - z, k20(z) k22 для zc + z z 1 и < 0 при dz zc - z < z < zc + z, k20(zc) = (k21 + k22)/2, k21 и k22 - положительные постоянные, k21 > k22.

Расчеты проводились для (E32 - E31)((1 + 2) arctan(y - yc) - (y - yc)) E31 + EE3(y) = +, 2((1 + 2) arctan - ) где |y - yc| < и (k22 - k21)((1 + 2) arctan(z - zc) - (z - zc)) k22 + Ek20(z) = +, 2((1 + 2) arctan - ) где |z - zc| < z.

Параметр k1 линейно зависит от PCO [Торр] - парциального давления СО в газовой фазе, и k1 = k10PCO.

В работе также рассматривается кинетическая модель реакции окисления водорода на никеле, которая описывает изменения безразмерных концентраций адсорбированных водорода (x) и кислорода (y) и кислорода, внедренного в приповерхностный слой металла (z) [1]:

= k1(1 - x - y)2 - k-1x2 - 2k30e- yx2y, 4 = k2(1 - x - y)2 - k40e- y+5zy - k30e- yx2y, (3) = (y(1 - z) - z(1 - x - y)).

Предполагается, что 0< 1 и z является медленной переменной. Значения всех параметров положительны.

Системы (1) и (3) рассматриваются в области = {(x, y, z) : x, y 0, x + y 1, 0 z 1} с границей , и 0 = \ - внутренность области .

Во второй главе приведены методы качественной теории динамических систем и вычислительной математики, используемые в работе для анализа семейств структурно устойчивых периодических решений и нерегулярной динамики кинетических моделей (1) и (3).

Третья глава содержит описание результатов исследования кинетических моделей (1) и (3). В работе доказано, что задача Коши для системы (1) или (3) с начальным условием x(0) = x0, y(0) = y0, z(0) = z0, где (x0, y0, z0) , имеет единственное решение x = x(t), y = y(t), z = z(t) и (x(t), y(t), z(t)) для всех 0 t < +.

В п.3.1 рассматривается вырожденная система относительно переменных x и y, которая получается из (1) при = 0 и содержит z в качестве параметра. Для случая, когда k3(y) const, в работе предложены достаточные условия существования глобально устойчивого стационарного состояния.

Кроме того, доказано, что если = 0 и параметры ki, i = 1, 21, 22, и z таковы, что вырожденная система (1) для каждого k3 = const, k3 [k32, k31], имеет единственное глобально устойчивое стационарное состояние (xs(k3), ys(k3)), и выполнено неравенство 4k-1k2(1 - ys(k3)) - k1(k1 + k-1) + 2k31ys(k3)(2k2(1 - ys(k3)) - k1) > 0, то существуют yc, ys(k32) < yc < ys(k31), и = такое, что для всех c < единственное стационарное состояние (xc, ys) системы (1)-(2), расs сматриваемой при = 0 в области 0, является неустойчивым, и существует замкнутая траектория, окружающая это стационарное состояние. Здесь c k32 = k30 exp(-E32/(RT )), k31 = k30 exp(-E31/(RT )) и ys (yc - , yc + ).

В работе выполнен параметрический анализ вырожденной системы (1) при = 0 с помощью методов теории бифуркаций. А именно, при фиксированных значениях ki, i = -1, 10, 20, 30, yc и , определены явные выражения PCO = PCO(y) и PO2 = PO2(y), где y [0, 1], для значений PCO и PO2, при которых в системе (1) при = 0 происходит бифуркация Андронова-Хопфа или существует негрубое стационарное состояние - седло-узел. Затем, рассматривая кривые, задаваемые этими выражениями в плоскости (PCO, PO2), определены значения параметров PCO и PO2, при которых существует три грубых стационарных состояния или грубый устойчивый предельный цикл.

Для однопараметрического семейства вырожденных систем (1) с параметром z построены максимальные семейства грубых периодических решений и стационарных состояний. Показано, что существуют значения параметров, при которых в пространстве (x, y, z) гладкая кривая ABCD стационарных состояний однопараметрического семейства вырожденных систем (1) при = 0 с параметром z имеет S-образный вид. Ветвь AB является максимальным семейством грубых устойчивых стационарных состояний. Стационарные состояния, лежащие на BC, неустойчивы и имеют тип седло. Точки B и C соответствуют стационарным состояниям, имеющим тип седло-узел. Стационарные состояния, образующие кривую CD, имеют тип узел или фокус.

s s Кроме того, существуют два максимальных семейства S1 и S2 грубых устойmin max чивых периодических решений вырожденной системы для z (z1, z1 ) и min max max min min max z (z2, z2 ), соответственно, где z1 < z2. При z = z1 и z = zпериодические решения зарождаются в результате бифуркации Андроноваmax min Хопфа. При z = z1 и z = z2 имеет место глобальная бифуркация вырождения устойчивого периодического решения в гомоклиническую траекторию - петлю сепаратрисы седла с отрицательной седловой величиной. Стаmin max ционарные состояния, лежащие на кривой CD для z1 < z < z2 (для min max z [0, 1]\(z1, z2 )), неустойчивы (устойчивы).

0.A 0.50.50.50.C 0.Ss 0.50.50.0.5Ss B 0.50.0.0.50.D 0.50.0.0.0.0.8 0 0.05 0.1 0.15 0.2 0.25 0.y Y X s Рис. 1: Слева: Кривая ABCD стационарных состояний и максимальные семейства Ss и S2 устойчивых периодических решений однопараметрического семейства вырожденных систем (1) при = 0. (1 и 2 - петли сепаратрис седловых стационарных состояний). Справа: Проекция притягивающего множества A системы (1) на плоскость (y, z) при = 4.и = 3.37 10-4. Здесь k1 = 0.0396, k-1 = 0.005, k21 = 9 104, k22 = 1.1 104, PO2 = 9 10-7, k30 = 1013, E31 = 28 103, E32 = 33 103, T = 500, zc = 0.5, = z = 0.1, yc = 0.101.

В п.3.2 на основании теоремы о существовании устойчивого предельного цикла в системе трех ОДУ с малым параметром (см. [2]) доказано, что существуют значения и такие, что соответствующая система (1) имеет Z z грубый устойчивый предельный цикл. Численно показано, что существуют значения параметров, при которых при увеличении имеет место каскад бифуркаций удвоения периода. Последовательность бифуркационных значений сходится к некоторому критическому = , и при > в фазовом пространстве системы существует притягивающее множество A, обладающее сложной структурой (см. рис.С1). Рассмотрение трансверсального сечения S множества A и отображения Пуанкаре P, которое переводит точки пересечения траекторий системы с S в их последующие, сводит исследование динамики системы в окрестности множества A к анализу некоторого одномерного отображения. Это отображение можно приблизить непрерывным унимодальным, имеющим отрицательную производную Шварца. Более того, оно топологически эквивалентно кусочно-линейному отображению которое имеет всюду плотную траекторию. Тем самым, отображение Пункаре P и система (1) в окрестности множества A могут обладать схожими свойствами.

В п.3.3 представлены результаты исследования влияния параметра k2, пропорционального парциальному давлению кислорода в газовой фазе, на структуру максимальных семейств структурно устойчивых периодических решений однопараметрического семейства вырожденных систем (3) с параметром z, то есть в случае = 0. С помощью метода продолжения периодического решения системы двух нелинейных ОДУ по параметру и формулы для вычисления первой ляпуновской величины негрубого стационарного состояния - сложного фокуса, описанных во второй главе, определен интервал значений k2, при которых существует максимальное семейство устойчивых периодических решений такое, что на обеих его границах периодическое решение вырождается в сложный устойчивый фокус в результате бифуркации Андронова-Хопфа. При увеличении k2 зарождается семейство неустойчивых периодических решений, существующих для z из интервала (z, zmax). Устойчивые периодические решения существуют для z из интервала (zmin, zmax), zmin < z. При z = z и z = zmin в системе происходит бифуркация Андронова-Хопфа. При z = zmax устойчивый и неустойчивый циклы сливаются в полуустойчивый (см. рис. 2). При дальнейшем увеличении k2 максимальное семейство устойчивых периодических решений распадается на два так, что на одной из границ каждого семейства происходит бифуркация вырождения периодического решения в петлю сепатрисы седла.

Отметим, что при рассматриваемых значениях ki и j в случае, когда k2 >> k1, в фазовом пространстве системы существует область, где переменные x и y вырожденной системы (3) являются медленной и быстрой, -x A 5.0.380.30.0.5 y 0.C x 0.4.Ss 0.4.0.4.4.0.0.B 0.D 4.0.0.0.8 0.0.0.0.4 0.0.0.0.2 0 0.x 0.x y y Рис. 2: Слева: Максимальное семейство Ss устойчивых периодических решений и семейство ABCD стационарных состояний однопараметрического семейства систем (3) при = 0 с параметром z (Вкладка: Семейства устойчивых (сплошные линии) и неустойчивых (пунктирные линии) уток-циклов). Справа: Притягивающее множество A системы (3) при = 1 и = 8.837 10-7. Здесь k1 = 0.2, k-1 = 0.01, k2 = 8, k30 = 100, k40 = 2, 3 = 30, 4 = 12, 5 = 10.

соответственно. В этом случае медленная кривая, которая является множеством точек фазовой плоскости, где скорость изменения быстрой переменной обращается в нуль, имеет S-образный вид. Две ее крайних ветви состоят из устойчивых, а средняя - из неустойчивых стационарных состояний уравнения быстрых движений. Более того, если z таково, что вырожденная система (3) имеет предельный цикл, то он содержит дуги, лежащие в окрестности устойчивой и в окрестности неустойчивой части медленной кривой, и называется уткой-циклом. Такие решения обладают высокой параметрической чувствительностью.

Для оценки глобальной погрешности численного интегрирования уток-циклов в работе рассмотрен метод, основанный на анализе главного члена v(t) асимптотического разложения по малому параметру глобальной погрешности численного интегрирования ( определяется порядком метода и величиной шага сетки). Построение оценки проводилось на временном интервале, равном нескольким периодам, путем вычисления функции v(t) на последовательности вложенных сеток. При интегрировании устойчивого уткицикла максимум и минимум нормы v(t) на периоде линейно возрастают с увеличением числа периодов. На неустойчивом утке-цикле эта зависимость аппроксимируется степенной функцией с основанием равным мультипликатору периодического решения, отличному от единицы.

В п.3.4 приведены результаты исследования влияния структуры решений вырожденной системы на динамику системы (3) при = 0. Рассмотрен z z z - 0.382случай, когда однопараметрическое семейство вырожденных систем (3) с параметром z имеет семейства грубых устойчивых и неустойчивых уток-циклов.

Численно показано, что существуют значения параметра , при которых выполнены условия теоремы о существования устойчивого предельного цикла в системе трех уравнений с малым параметром (см. [2]). Более того, при таких получены значения параметра , при которых динамика системы является нерегулярной, и динамика системы усложняется в результате каскада бифуркаций удвоения периода. При этих значениях в фазовом пространстве системы (3) существует притягивающее множество A (например, см. рис. 2).

Путем построения отображения Пункаре показано, что исследование динамики системы в окрестности этого множества можно свести к анализу одномерного отображения, которое можно приблизить непрерывным унимодальным отображением с отрицательной производной Шварца. Показано, что это отображение имеет цикл длины 3. На основании теоремы Шарковского, отсюда следует, что отображение имеет циклы сколь угодно больших периодов.

Для уточнения гомоклинической траектории - петли сепаратрисы седла в системе двух нелинейных ОДУ в п.3.5 предлагается следующий алгоритм. Пусть дана система двух нелинейных ОДУ с параметром z R:

= f(x, z), x R2 (4) где вектор-функция f является гладкой. Предположим, что система (4) имеет стационарное состояние - седло p(z) для значений z I, так что собственные значения матрицы Якоби 1(z) < 0 и 2(z) > 0. Пусть при z = z0 I состояние равновесия p0 = p(z0) имеет гомоклиническую траекторию L0, стремящуюся к седлу p0 как при t -, так и при t +. Для нахождения бифуркационного значения z0 и уточнения петли L0 будем решать следующую краевую задачу: При заданных 0 < 1 и собственных векторах vj(z), соответствующих j(z), j = 1, 2, необходимо найти такие параметры z и tj и функции j(t), что :

j = (-1)jf(j, z), j(0) = p(z) + vj(z), 1(t1) = 2(t2) l, где ||vj(z)|| = 1 и l - заданное трансверсальное сечение.

Для решения поставленной краевой задачи рассматривается следующий итерационный алгоритм. Пусть z близко к z0 и xj(tj, aj(z), z) решение системы j = (-1)jf(xj, z), удовлетворяющее начальному условию xj(0, aj(z), z) = aj(z), где aj(z) = p(z) + vj(z) и xj(tj, aj(z), z) l, j = 1, 2.

Предположим, что - возмущение параметра z и j - возмущения моментов времени tj, j = 1, 2, такие что x1(t1 + 1, a1(z + ), z + ) = x2(t2 + 2, a2(z + ), z + ) l.

Раскладывая последнее равенство в ряд по и j и оставляя линейную часть, определим следующее приближение для бифуркационного значения параметра. Для вычисления vj(z + ), j = 1, 2, рассмотрим проектор P (z + ) на устойчивое инвариантное многообразие седла p(z + ). Учитывая разложе(1) ние в ряд P (z + ) = P (z) + P (z) + O(2) получим, что vj(z + ) = (1) vj(z) + (-1)j+1P (z)vj(z) + O(2), j = 1, 2. Для вычисления проекторов P (z) и P (z + ) в работе рассмотрено их представление в виде интегралов по некоторому контуру на комплексной плоскости. Доказано, что оператор (1) P (z) является вещественным и для него получено явное выражение.

Данный алгоритм реализован в виде комплекса программ на языке Фортран, его применение при исследовании кинетических моделей позволило определить гомоклинические траектории и бифуркационные значения параметров с достаточной точностью.

Выводы 1. Для системы двух нелинейных ОДУ с параметром разработаны алгоритм и комплекс программ для уточнения гомоклинической траектории - петли сепаратрисы седла, в котором используются проекторы на собственные подпространства линеаризованной задачи в окрестности стационарного состояния. Доказано, что оператор, определяющий линейную часть возмущения проектора, является вещественным.

2. Для описания динамики каталитической реакции окисления оксида углерода на металлах платиновой группы предложена и исследована кинетическая модель с тремя переменными, одна из которых медленная.

Модель учитывает резкое изменение энергии активации взаимодействия адсорбированных веществ, а также константы скорости адсорбции кислорода в некотором узком интервале значений концентраций адсорбированного и внедренного кислорода, соответственно. Предложены достаточные условия, при которых вырожденная система относительно быстрых переменных имеет единственное стационарное состояние и замкнутую траекторию. Численно показано, что в полной системе возможно развитие нерегулярной динамики в результате каскада бифуркаций удвоения периода при увеличении константы скорости внедрения кислорода в приповерхностный слой катализатора.

3. Исследованы структурно устойчивые периодические решения кинетической модели окисления водорода на никеле с двумя переменными: Получены значения параметров, при которых система имеет устойчивый и неустойчивый предельные циклы, обладающие высокой параметрической чувствительностью (утки-циклы). При численном интегрировании уток-циклов на одном периоде получена оценка глобальной погрешности дискретизации. Показана асимптотическая зависимость глобальной погрешности от числа периодов для устойчивого и неустойчивого решения при интегрировании системы на длительном интервале времени.

4. Показано, что нерегулярные колебания в кинетической модели окисления водорода на никеле, являющейся системой трех нелинейных ОДУ с одной медленной переменной, зарождаются в результате каскада бифуркаций удвоения периода при изменении параметра (структурные характеристики катализатора).

Список литературы С [1] Чумаков Г.А., Слинько М.Г. Кинетическая турбулентность (хаос) скорости реакции взаимодействия водорода с кислородом на металлических катализаторах. Докл. АН СССР, 1982, 266, 5, 1194-1198.

[2] Понтрягин Л.С., Родыгин Л.В. Периодическое решение одной системы обыкновенных дифференциальных уравнений с малым параметром при производных. Докл. АН СССР, 1960, 132, 3, 537-540.

Список публикаций в рецензируемых журналах 1. Лашина Е.А., Чумаков Г.А., Чумакова Н.А. Максимальные семейства периодических решений кинетической модели гетерогенной каталитической реакции. Вестник НГУ, серия: математика, механика, информатика,т. V, вып. 4, 2005, с. 3-20.

2. Ivanova (Lashina) E.A., Chumakova N.A., Chumakov G.A., Boronin A.I.

Modeling of relaxation oscillations in CO oxidation on metallic catalysts with consideration of reconstructive heterogeneity of the surface. Chem. Eng. J., v. 107, 2005, 191-197.

3. Lashina E.A., Chumakova N.A., Chumakov G.A., Boronin A.I. Chaotic dynamics in the three-variable kinetic model of CO oxidation on platinum group metals. Chem. Eng. J., v. 154, 2009, 82-87.

Список публикаций в трудах конференций 1. Иванова (Лашина) Е.А., Чумакова Н.А. Оценка глобальной ошибки дискретизации на периодических решениях и решениях-утках одной кинетической модели. Тезисы докладов международной конференции молодых ученых по математическому моделированию и информационным технологиям. Новосибирск: ИВТ СО РАН, 2002, с. 29-30.

2. Иванова (Лашина) Е.А., Чумакова Н.А. Максимальные семейства периодических решений кинетической модели каталитического окисления водорода. Тезисы докладов международной школы-конференции молодых ученых по катализу Каталитический дизайн - от исследований на молекулярном уровне к практической реализации. Новосибирск: ИК СО РАН, 2002, с. 207-208.

3. Ivanova (Lashina) E.A., Chumakova N.A., Chumakov G.A. An algorithm for the saddle-loop homoclinic orbit finding in two-dimensional kinetic model.

Труды международной конференции по вычислительной математике МКВМ-2004. Ч.II, Новосибирск: Изд. ИВМиМГ СО РАН, 2004, 870-875.

4. Лашина Е.А., Чумаков Г.А., Чумакова Н.А. Об одном алгоритме уточнения петли сепаратрисы седла. Тезисы докладов Международной конференции, посвященной 100-летию со дня рождения С.Л. Соболева Дифференциальные уравнения, функциональные пространства, теория приближений, Новосибирск: ИМ СО РАН, 2008, с. 517.

5. Лашина Е.А., Чумакова Н.А., Чумаков Г.А. Оценка глобальной погрешности численного интегрирования на периодических решениях-утках.

Тезисы докладов всероссийской конференции, приуроченной к 80-летию академика С.К. Годунова Математика в приложениях, Новосибирск:

ИМ СО РАН, 2009, с. 170-171.

6. Лашина Е.А., Чумакова Н.А., Чумаков Г.А. Хаотическая динамика одной кинетической модели гетерогенной каталитической реакции. Тезисы докладов II Международной школы-семинара Нелинейный анализ и экстремальные задачи, Иркутск: ИДСТУ СО РАН, 2010, с. 45.

ашина Елена Александровна Анализ структурно устойчивых периодических решений кинетических моделей каталитических реакций.

Автореф. дисс. на соискание учёной степени кандидата физико-математических наук.

Подписано в печать 15.02.2012. Заказ №. Формат 60х84/16. Усл. печ. л. 1. Тираж 100 экз.

Отпечатано на полиграфическом участке издательского отдела Института катализа СО РАН 630090, Новосибирск, пр-т Академика Лаврентьева,    Авторефераты по всем темам  >>  Авторефераты по техническим специальностям