Прямая и обратная задачи для сингулярной системы с медленными и быстрыми переменными в химической кинетике Текст научной статьи по специальности «Математика»

Владикавказский математический журнал 2015, Том 17, Выпуск 1, С. 39-46

УДК 541.124:541.126:517.9

ПРЯМАЯ И ОБРАТНАЯ ЗАДАЧИ ДЛЯ СИНГУЛЯРНОЙ СИСТЕМЫ С МЕДЛЕННЫМИ И БЫСТРЫМИ ПЕРЕМЕННЫМИ В ХИМИЧЕСКОЙ КИНЕТИКЕ1

Л. И. Кононенко

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

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

Введение

Моделирование многих процессов в различных областях научной и практической деятельности (геофизика, химическая кинетика, медицина и т. д.) приводит к решению прямых и обратных задач и существует давно. Но как раздел современной прикладной математики — обратные задачи естествознания — появились сравнительно недавно.

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

Задача отыскания решений системы по некоторым исходным данным при известных функциях в правых частях представляет собой так называемую прямую задачу для дифференциальных уравнений. Основная цель исследований в данной работе — постановка и анализ задачи, обратной к этой (в химической кинетике обратные задачи называются задачами идентификации). Обратная, задача сводится к нахождению неизвестных правых частей по некоторым данным о решении прямой задачи. Сформулирована и обоснована гипотеза о том, что по заданию решения на медленной поверхности можно восстановить неизвестные правые части, в частности, получены условия существования и единственности коэффициентов в правой части медленной подсистемы вырожденной системы, заданной полиномом. Начнем с краткого обзора по теории систем дифференциальных уравнений с медленными и быстрыми переменными. Мы будем иметь дело с прямыми и обратными задачами для обыкновенных дифференциальных уравнений [1—9] в отличие от изучения теории обратных задач для уравнений с частными производными, которые представлены, например, в [10-16], причем в системах, рассматриваемых нами, присутствует малый параметр.

© 2015 Кононенко Л. И.

1 Работа поддержана Российским фондом фундаментальных исследований, проект № 15-01-00745, и Сибирским отделением РАН, междисциплинарный интеграционный проект № 80.

В [14], где собраны почти все основные направления теории обратных и некорректных задач, рассматривается также система линейных обыкновенных дифференциальных уравнений, решение задачи Коши которой описывает процесс химической кинетики

^ = qí1Ul{t) + ^2и2(£) Н-----Ь <?гпи„(Г),

«¿(0) =Тг, г = 1,2

где щ (Ь) — концентрация г-го вещества в момент времени Ь.

Постоянные параметры характеризуют зависимость скорости изменения Пг(Ь) от концентраций веществ, участвующих в процессе.

Прямая задача формулируется следующим образом: определить щ(Ь), зная параметры и^в начальный момент времени.

Обратная задача сводится к нахождению параметров по решению системы Пг(Ь), которое соответствует начальному условию Заметим, что иногда тоже неизвестны и их нужно определить вместе с .

Поскольку величины, участвующие в процессах химической кинетики, разномасштабны, мы будем рассматривать системы с малым параметром, описывающие эти процессы и представленные в таких работах, как [17-23].

Рассматривается сингулярно возмущенная система обыкновенных дифференциальных уравнений

х = / (х,у,г,е), еу = д(х,у,Ь,е),

где х £ Жт — медленные, у £ Ж™ — быстрые переменные; /, д — достаточно гладкие функции, Ь £ Ж — переменная, имеющая смысл времени, X, у — производные по времени, е — положительный малый параметр.

Система рассматривается в ограниченной выпуклой инвариантной притягивающей области Ш С Жт х Ж™.

1. О прямой задаче для системы (1)

Прямая задача для системы дифференциальных уравнений с малым параметром (1)

сводится к отысканию пары функций х(Ь), у(Ь), удовлетворяющих системе (1), по неко-

/д

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

Гладкая поверхность Б в Жт х Ж™ х Ж называется интегральным многообразием системы (1), если любая траектория этой системы, имеющая хотя бы одну общую точку с Б, целиком принадлежит поверхности Б. Формально, если при Ь = Ьо точка (х(Ьо),у(Ь0),Ь0) £ Б, то траектория (х(Ь), у(Ь), Ь) целиком принадлежит Б.

Уравнение х = /(х,у,Ь,е) в системе (1) называется медленной подсистемой , а уравнение еу = д(х, у, Ь, е) — быстрой подсистемой системы (1). Если в (1) положить е = 0, получим порождающую или вырожденную систему

х = /(х,у,Ь, 0), (2)

0 = д(х,у,Ь, 0). (3)

Уравнение д(х, у, 0) = 0 задает медленную поверхность. Это уравнение медленной поверхности может иметь одно или несколько решений, каждое из которых задает лист медленной поверхности,.

Опишем подробнее. Поверхность, задаваемая уравнением (3), называется медленной поверхностью. Возьмем какое-нибудь решение х(£о), у (¿о) ¿о уравнения д(х,у,£, 0) = 0, т. е. выберем точку (х(£о),у(^о), ¿о) на медленной поверхности Г. Если ¿о, 0)) ф 0, то в некоторой окрестности Уо точки (ж(£о),^о) £ х Ж,

С Мт, существует по теореме о неявной функции [24] вектор-функция у = Но(х,£), уо = Но(хо,Ьо^ у ^ С Ж", которая является решением уравнения (3), т. е.

д(х, Но(х, ■£),■£, 0) = 0 при всех (х,£) £ Уо- Пересечение поверхности Г с поверхностью 0)) = 0 является поверхностью (кривой) Г1 на единицу меньшей размерности, чем Г Она делит медленную поверхность на части, в каждой из которых 0)) не меняет знака. Каждую из этих частей называют листом медленной поверхности,.

Листы интегрального многообразия медленных движений (или медленного интегрального многообразия) являются уточнением при учете малого параметра е листов медленной поверхности и получаются из них с помощью асимптотического разложения

е

Н(х, ¿,е) = Но (х, ¿) + еН^х, ¿) +-----+ екНк(х, ¿) + ..., (4)

где коэффициенты разложения Нк (х,£) подсчитываются по рекуррентной формуле, приведенной, например, в [9],

hk = -B

k-1

-1

(5)

Jk) _ ahk-1 _ dhp f(k-l-p) / dt dx J

p=0

В = det(^(x,ho(x,t),t,0)) /0.

Среди интегральных многообразий системы (1) нас интересуют m-мерные интегральные многообразия (размерность медленных переменных), которые представимы в виде графика вектор-функции y = h(x,t, е). Выполняется соотношение

lim h(x, t, е) = h0(x, t),

где ho(x,t) — функция, график которой является листом медленной поверхности.

е

к траекториям вырожденной системы.

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

Согласно методу интегральных многообразий мы должны сделать следующее:

1) исследовать строение медленной поверхности, т. е. найти количество и форму листов;

2) найти границы листов;

3) проверить характер устойчивости листов;

4) провести качественный анализ медленных подсистем на устойчивых листах (в частности, найти стационары, их классификацию; особый интерес вызывают колебания различных типов и решения-утки);

5) провести анализ системы в целом.

Нахождение решения системы (1) сводится к отысканию решения вырожденной сие

нулю. Этот факт следует из работ А. Н. Тихонова [4, 5], в которых доказаны теоремы о предельном переходе к решению вырожденной задачи при стремлении малого параметра к нулю. Правые части системы (1) являются достаточно гладкими функциями, поэтому удовлетворяют требуемым условиям, в частности, обеспечивают единственность решения.

При использовании метода интегральных многообразий для решения конкретных задач центральным становится вопрос о вычислении функции /г(х, Ь, е), описывающей многообразие. Как правило, точное вычисление не является возможным и используются различные виды приближенных вычислений. Мы будем использовать для приближенного вычисления асимптотическое разложение функции /г(х, Ь, е) по степеням малого параметра, приведенное в формулах (4), (5).

Существование интегрального многообразия было доказано в [3,9]. Приведем соответствующую теорему и сформулируем постановку прямой задачи. (Заметим, что термин «прямая» задача обычно употребляется в контексте с «обратной» задачей, а в остальных случаях говорят просто «задача».)

Постановка прямой задачи. Пусть для системы (1) выполнены условия:

I. Уравнение д(х,у,Ь, 0) = 0 имеет изолированное решение у = До(х,Ь) при Ь £ Ж, х £ Жт

II. В области ^о = |(х,у,Ь,е) : х £ Жт, ||у — Л,о(х,Ь)|| < р, Ь £ Ж, 0 ^ е ^ ео | функции /, д и До равномерно непрерывны и ограничены вместе с частными производными по переменным до (к + 2)-го порядка включительно (к ^ 0).

III. Собственные значения (г = 1,...,«) матрицы /го(ж, 0) подчиняются неравенству Яе А^(х,Ь) ^ —27 < 0.

Требуется по заданным функциям /(х,у,Ь, е), д(х,у,Ь,е) в правой части системы (1) найти х(Ь), у(Ь) в области ^о-

Теорема. Пусть выполняются условия 1-111. Тогда существует такое е1 (0 < е1 ^ ео), что для каждого е £ (0,е1 ] система (1) имеет интегральное многообразие медленных движений у = Д(х, Ь, е), движение по которому описывается уравнением

х = /(х, Д(х, Ь, е), Ь, е). (6)

Если х(Ь) — решение этого уравнения, то пара х(Ь), у(Ь), где у(Ь) = Д(х, Ь, е), является решением исходной системы (1), т. е. пара х(Ь), у(Ь) есть решение прямой задачи.

В качестве примеров прямой задачи были рассмотрены две модели из химической кинетики.

I. Математическая модель реактора идеального смешения.

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

х 1 = а — х1 — а[ш1 + (шз — ш1 — ш2)х1], х 2 = Ь — х2 — + Ш4 + (шз — Ш1 — Ш2)х2 ], (7)

у1 = в(2Ш1 — Шз — Ш4), у/2 = в (Ш2 — шз),

ГД6Ш1 = К^(1 — у1 —у2)2 — К— 1 у2, Ш2 = К2х2(1 —у1 — у2) — К—2у2, Шз = у^, Ш4 = К4х2у1 — обезразмеренные скорости четырех стадий реакции.

Областью изменения переменных является множество

W = |(хьх2,уъу2) : 0 ^ х1 ^ а, 0 ^ х2 ^ Ь, 0 ^ у1, 0 ^ у2, у1 + у2 ^ 1}.

Система (7) изучалась в [6, 7].

Скорость реакции на поверхности катализатора существенно выше, чем скорости адсорбции. Предполагается, что основным механизмом реакции является адсорбционный, а ударный механизм учитывается как дополнительный. Поэтому мы используем при анализе модели следующую иерархию параметров: к-2, к-1, к4 <С К1, к2 <С 1. Константы десорбции предполагаются малыми сравнительно с константами адсорбции. Кроме того, а ^ ^ие = 1/в ^ к-2, к-1; к4.

II. Математическая, модель каталитической реакции окисления.

Рассматривается детальный механизм реакции СО + О2 на иридии. Кинетической схеме этой реакции соответствует система дифференциальных уравнений с безразмерными переменными

х 1 = 2Ь1х22 — Ь2 х6х1 — 6вх1х2, х 2 = Ь4 х7 — Ь5х2 — Ьв х1х2 — 6дх2хз — Ь12 х2х4, х 3 = Ь2х6х1 — 2Ьз х3 — Ьбхз + х5 — Ьдх2хз + 2Ью х4х5 + &12х2х4, (8)

х4 = 263х2 — Ь1о х4х5 — Ь12х2х4,

х 5 = Ьб хз — &7х5 — 6юх4х5 — ЬИх5,

где хб = 1 — хз — х4 — х5, х7 = 1 — х1 — х2 — хз — х4 — х5.

Система (8) исследовалась в [8, 17, 18]. Там же приведены выражения для коэффициентов Ь (г = 1, 2,... , 12). Областью изменения переменных является множество

W = ,...,х5) : 0 < хг < 1, ^ < 1, г = 1,..., Л.

I ¿=1 J

При анализе модели используем следующую иерархию параметров: Ью > Ьв » Ь7 > Ь1,Ь2,Ьз ,Ь4,Ьб, Ьц ,Ь12 » Ь5, Ьд.

В [19] также проведен подробный анализ систем (7) и (8) с применением метода интегральных многообразий, который позволяет провести исследование согласно шагам, описанным ранее. В работах [20-23] так же рассматривались системы с малым параметром, в частности, в [20-22] — системы, описывающие задачи химической кинетики.

2. Обратная задача для системы

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

е=0

системы (2);

2) функция / в правой части медленной подсистемы системы (1) задается в виде многочлена р-й степени / = ^Ьцхгу^, так как в задачах химической кинетики правые части системы часто являются полиномами, более того будем рассматривать многочлен первой степени;

3) рассматриваются системы с одной медленной и одной быстрой переменными, т. е. т = п = 1;

4) функцию д(х, у, Ь, е) считаем заданной и удовлетворяющей всем условиям теоремы о неявной функции в каждой точке области, в частности, ф о, следовательно, при е = 0 медленная поверхность, уравнение которой д(х, у, Ь, 0) = 0, задана;

5) медленная поверхность состоит из одного листа.

В рамках сделанных ограничений вырожденная система

х = ^ Ь^ хУ

¿+КР (9)

0 = д(х,у,Ь, 0)

имеет вид

х = Ро + Р1х + Р2у, (10)

0 = д(х,у,Ь, 0).

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

Учитывая перечисленные ограничения 1)-5), получим следующую постановку обратной задачи для сингулярно возмущенной системы.

Постановка обратной задачи. Рассматривается система дифференциальных уравнений со следующими данными

х = Ро + Р1х + Р2у,

0 = д(х, у, Ь, 0), х £ Ж, у £ Ж, Ь £ Ж, (х,у) £ Ш, (11)

х(Ь%) = а%, х(^) = в%, г = 1,2,3,

где функция д удовлетворяет условиям 4), 5). Для данных а%, в% требуется найти коэффициенты Ро, р>1, Р2, удовлетворяющие системе (11).

Из второго уравнения системы (10) при условии ф о выразим быструю

у х д

влетворяет всем условиям теоремы о неявной функции. Имеем у = ^>(х,Ь). Подставив это выражение в первое уравнение системы (11), имеем

х(Ь) = ро + Р1х(Ь)+ р2^(х(Ь),Ь). (12)

Используя данные х(Ь%) = а%, х(Ь%) = в% и вводя обозначения 7% = у% = ^>(а%,Ь%), г = 1, 2, 3, имеем систему трех уравнений с тремя неизвестными Ро, р1, Р2:

Ро + Р1«1 + Р2 71 = в1,

Ро + Р1«2 + Р272 = в2, (13)

Ро + Р1«з + Р2 7з = вз.

Значит, на медленной поверхности (е = 0) искомые коэффициенты р%, г = 0,1, 2, вы-

д0

дГ-

числяются по следующим формулам: р^ = Здесь А0 — определитель системы (13):

1 «1 71

Д° = 1 а2 72

1 аз 7з

а АР0, Др!; Др2 — определители, получаемые из А0 заменой у'-го столбца, у = 1,2,3, соответствующим столбцом коэффициентов при ро, Ръ Р2:

А0

ро

в1 а1 71 1 в1 71 1 а1 в1

в2 а2 72 , Д0 = ' 0! 1 в2 72 , д0 = ' 02 1 а2 в2

вз аз 7з 1 вз 7з 1 аз вз

(15)

Нуль в верхнем индексе определителей означает, что коэффициенты рг, Я = 0,1, 2, подсчитываются на медленной поверхности (при е = 0).

Из теории линейных алгебраических систем [26, 27] вытекает следующее условие существования и единственности коэффициентов рг;. До = 0.

Имеет место следующее достаточное условие существования и единственности коэффициентов системы (13), которая была получена из вырожденной системы (10).

Предложение. Пусть данные и, а^ вг, Я = 1, 2, 3, таковы, что выполнены условия

1)-5) и определитель системы (13) при этих значениях Я = 1, 2, 3, не равен нулю.

Тоща обратная задача имеет единственное решение, которое определяется равенствами до

(М = -дтг, где А и Д04, я = 1,2,3, вычисляются по формулам (14) и (15).

В дальнейших исследованиях мы предполагаем включить в рассмотрение случай

е=0

чество листов медленной поверхности.

Автор выражает благодарность В. Г. Романову за давнюю совместную статью [28], которая инициировала данную, коллегам за помощь в работе и рецензенту за справедливую критику и ценные замечания.

Литература

1. Митрополъский К). А., Лыкова О. В. Интегральные многообразия в нелинейной механике.—М.: Наука, 1963.

2. Васильева А. В., Бутузов В. Ф. Сингулярно возмущенные уравнения в критических случаях.—М.: Изд-во МГУ, 1978.

3. Голъдштейн В. М., Соболев В. А. Качественный анализ сингулярно возмущенных систем.— Новосибирск: Изд-во Ин-та математики СО АН СССР, 1988.

4. Тихонов А. Н. Системы дифференциальных уравнений, содержащие малые параметры при производных // Мат. сб.—1952.—Т. 31(73), № З.-С. 575-586.

5. Тихонов А. Н. О зависимости решений диффернциальных уравнений от малого параметра // Мат. сб.—1948.—Т. 22(64), № 2.-С. 193-204.

6. Голъдштейн В. М., Кононенко Л. И., Лазман М. 3., Соболев В. А., Яблонский Г. С. Качественный анализ динамических свойств каталитического изотермического реактора идеального смешения // Математические проблемы химической кинетики.—Новосибирск: Наука, 1989.—С. 176-204.

7. Кононенко Л. И. О гладкости медленной поверхности сингулярно возмущенных систем // Сиб. журн. индустр. математики.—2002.—Т. 5, № 2(10).—С. 109-125.

8. Кононенко Л. И., Волокитин Е. 17. Параметризация и качественый анализ сингулярной системы в математической модели реакции каталитического окисления // Сиб. журн. индустр. математики.—2012,—Т. 15, № 1(49).—С. 43-52.

9. Воропаева Н. В., Соболев В. А. Геометрическая декомпозиция сингулярно возмущенных систем.— М.: Физматлит, 2009.

10. Романов В. Г. Обратные задачи для гиперболических систем // Вычислительные методы в математической физике, геофизике и оптимальном управлении.—Новосибирск: Наука, 1978.—С. 128-142.

11. Лаврентьев М. М., Романов В. Г., Шишатский С. П. Некорректные задачи математической физики и анализа. —М: Наука, 1980.—286 с.

12. Аниконов К). Е. Несколько вопросов теории обратных задач для кинетических уравнений // Обратные задачи математической физики.—Новосибирск, 1985.—С. 28-41.

13. Алексеев А. С. Некоторые обратные задачи теории распространения волн 1, 2 // Изв. АН СССР. Сер. геофизика.—1962,—№ 11,—С. 1514-1531.

14. Кабанихин С. И. Обратные и некорректные задачи.—Новосибирск: Сибирское научное изд-во, 2009.-456 с.

15. Голубятников В. 17. Обратная задача для уравнения Гамильтона — Якоби на замкнутом многообразии // Сиб. мат. журн.—1997.—Т. 38, № 2.—С. 276-279.

16. Кожанов А. И. Нелинейные нагруженные уравнения и обратные задачи // Журн. вычисл. математики и мат. физики.—2004.—Т. 44, № 4.—С. 694-716.

17. Гайнова И. А., Фадеев С. И., Елохнн В. И., Воронин А. И. Реакция окисления СО на поликристаллической фольге иридия. Моделирование кинетики поверхностных процессов //Междунар. конф. по вычислительной математике. Труды конференции. Часть I.—Новосибирск, 2004.—С. 449-454.

18. Кононенко Л. И. Интегральные многообразия в математической модели реакции каталитического окисления // Изв. РАЕН. Сер. MMMIIY. 1997. Т. 1, № 4.-С. 53-59.

19. Кононенко Л. И. Медленные поверхности в задачах химической кинетики // Мат. заметки ЯГУ.— 2012.—Т. 19, вып. 2.-С. 49-67.

20. Chumakov G. A., Chumakova N. A. Relaxation oscillations in a kinetic model of catalytic hydrogen oxidation involving a chase on canards // Chem. Eng. J.—2003.—Vol. 91, № 2-3.—P. 151-158.

21. Кононенко Л. И. Релаксации в сингулярно возмущенных системах на плоскости // Вести. НГУ. Сер. Математика, механика, информатика.—2009.—Т. 9, вып. 4.—С. 45-50.

22. Соболев В. А., Щенакина Е. А. Интегральные поверхности со сменой устойчивости и траектории-утки // Изв. РАЕН. Сер. МММИУ.-1997.-Т. 1, № З.-С. 176-187.

23. Вобкова А. С., Колесов А. Ю., Розов Н. X. Проблема «выживания уток» в трехмерных сингулярно возмущенных системах с двумя медленными переменными // Мат. заметки.—2002.—Т. 71, вып. 6.— С. 818-831.

24. Решетник Ю. Г. Курс математического анализа. Ч. 1. Кн. 2: Интегральное исчисление функций одной переменной. Дифференциальное исчисление функций многих переменных.—Новосибирск: Изд-во Института математики им. С. Л. Соболева СО РАН, 1999.—512 с.

25. Ермакова А. Методы макрокинетики, применяемые при математическом моделировании химических процессов и реакторов. Институт катализа СО РАН им. Г. К. Борескова.—Новосибирск,

2001.-188 с.

26. Годунов С. К. Современные аспекты линейной алгебры.—Новосибирск: Научная книга, 1997.— 388 с.

27. Годунов С. К. Лекции по современным аспектам линейной алгебры. Новосибирск: Научная книга,

2002.-216 с.

28. Романов В. Г., Слинючева Л. И. Обратная задача для линейных гиперболических систем первого порядка // Матем. проблемы геофизики.—Новосибирск: Издательство ВЦ СО АН СССР, 1972.— Вып. З.-С. 187-215.

Статья поступила 28 сентября 2014 г.

Кононенко Лариса Ивановна Институт математики им. С. Л. Соболева СО РАН старший научный сотрудник лаб. прикладного анализа РОССИЯ, 630090, Новосибирск, пр. Ак. Коптюга, 4;

Новосибирский государственный университет, доцент РОССИЯ, 630090, Новосибирск, Академгородок, ул. Пирогова, 2 E-mail: larakon2@gmail.ru, larak@math.nsc.ru

DIRECT AND INVERSE PROBLEMS FOR A SINGULAR SYSTEM WITH SLOW AND FAST VARIABLES IN CHEMICAL KINETICS

Kononenko L. I.

Direct and inverse problems for singular systems with small parameter are stated, which describe catalytic reactions in chemical kinetics. The solution of the direct problem is based on the method of integral manifolds. The inverse problem reduces to finding the coefficients of the polynomial in the right-hand part of the slow equation according to the solution given on the slow surface of the system. The above arguments make it possible to obtain existence and uniqueness conditions for the coefficients in the right-hand part of the slow subsystem of the degenerate system.

Key words: mathematical modeling, singularly perturbed system, integral manifold, slow surface, inverse problem.