CC BY
15
УДК 517.9
DOI: 10.14529/mmp 170108
СТАЦИОНАРНЫЕ ТОЧКИ УРАВНЕНИЯ «РЕАКЦИЯ-ДИФФУЗИЯ» И ПЕРЕХОДЫ В СТАБИЛЬНЫЕ СОСТОЯНИЯ
A.C. Коротких
Рассмотрена бесконечномерная динамическая система, заданная уравнением «реакция-диффузия» с кубической нелинейностью при краевом условии Неймана и фиксированном значении средней величины. Изложена методика приближенного вычисления бифурцирующих решений при малых и конечных значениях закритическо-го приращения параметра. Предложена также методика «трассировки» траекторий спуска из произвольного состояния (с произвольной концентрацией) в стабильное состояние (с концентрацией, реализующей минимум функционала энергии). Методика основана на вычислении сужения функционала энергии на линейную оболочку основных собственных функций (мод) оператора Лапласа и приближенном построении трассы спуска в виде последовательности точек, сопровождающих траекторию динамической системы. В случае малого закритического приращения бифуркационного параметра вычислены асимптотические представления бифурцирующих решений. В случае конечного закритического приращения бифуркационного параметра приведены примеры вычисления трассы спуска в точки минимума функционала энергии.
Ключевые слова: уравнение «реакция-диффузия»; гладкий функционал; экстремаль; бифуркация стационарных состояний; моды бифуркации; вариационный метод Ляпунова - Шмидта; траектории спуска в точки минимума.
Введение и постановка задачи
Как известно, структурные изменения вещества изучаются посредством специальным образом подобранных м одел ь н ых дифференциальных уравнений при кра-бвых условиях Пуассона. Например, при изучении динамики спонтанного разделения фаз (бинарного) вещества, (сплава), при исследовании стабильных распределении концентраций вещества и т.п. часто используется градиентная динамическая система типа «реакция-диффузия> [1-3]
V = -^аЛ (V (т) - С И) = Ат + Хт - т3 - С, (1)
в которой т = т(х,Ь) - относительная концентрация компоненты вещества,
хе и с к2,
V(w) := JJ i i-^- - A— + -jj dx1dx2
и
- интеграл энергии, и - область, занятая веществом, Ь - время. Как правило, разыскиваются такие решения т = т(х,Ь), для которых выполнено граничное условие Неймана
дт
7Т =0' дп ди
где ди - граница области и, п векторное поле нормалей к границе.
Первоочередной интерес представляют стационарные состояния - решения w = w(x) стационарного уравнения
Aw + Xw - w3 - C, C = const, (2)
или, что одно и то же, экстремали (критические точки) функционала
V(w) := V(w) - C [w] ,
где [w] := ff w dxidx2 - среднее значение w. и _
Поиск экстремалей V(w) можно осуществлять посредством редукции Ляпунова-Шмидта к конечномерной задаче grad W(£) = 0, £ Е Rn, где
W(£) := inf V(w)
(w,ej )=íj
- ключевая функция, построенная по начальным собственным функциям (модам бифуркациями) ej , j = 1, ... , m оператора Лапласа на области U (при заданных краевых условиях Пуассона) в соответствующим образом подобранном функциональном пространстве состояний [2, 3].
Обратимся к функциональным пространствам Hi0 , Hi1 , Hi2 , где H0 = L2(Q), Hík -замыкание по норме \\w\\k (к = 1, 2) множества тригонометрических полиномов вида
m
£p,qep,q , P,q — 0 , ep,q(x1, x2) - собственная функция оператора A (при краевом
p+q=1
условии Неймана ^ =0):
J дп ди '
ep,o(xi,x2) = V2cos(pnxi), eo,q(xi,x2) = V2cos(qnx2), ep,q (x1,x2) = 2 cos(pnx1) cos(qnx2), pq = 0 .
Собственная функция e<p q отвечает собственному значению Xp,q = —n2(p2 + q2), совокупность функций {ep,q} является ортонормированной системой векторов в L2(Q).
w \ w\ k
шениями:
\\w\1 := — JJ w A(w) dx1 dx2 = JJ |V(w) |2 dx1dx2 , n Q
\\w\2 := j j w A2(w) dx1dx2 = jj(A(w))2 dx1dx2 .
nn
Линейный оператор A := —A : Hi2 ——> Hi0 является изоморфизмом. При этом оператор A положительный и диагонализируемый:
A : qepq 1—y ixpq1 £pqep>q> P,q - 0' P + q - 1.
p,q p,q
Пусть B - положительный и симметричный оператор, для которого B2 = A (то есть B - положительный квадратный корень из A). Оператор B изоморфно действует из Hi2 на 'H1 и из H^ на Hi0 :
B : £p,qep,q1—у \jIXp,q1 £pqep,q = n £p,qeP,q, p + q - 1.
p q p q p q
Из этого соотношения следует, что ||ш||2 = \\Б (ш)^. Так как уравнение «реакция-диффузия> допускает операторное представление
ш = — grad V(ш) , ш еН2 , то заменой ш = Вш оно сводится к эквивалентному, но неградиентному, уравнению
—д __^
— = —В grad V(В-1 ш), гШ еН-2 . (3)
—Ь
Если рассмотреть градиент функционала д(гш) := V(В-1 ги) в тройке Н1 С Н1 С Н0 (см. определение градиента в [2, 3]), то можно рассмотреть градиентное уравнение
—шд _
— = В-1 grad V(В-1 ш), ш еН-1. (4)
Стационарные решения уравнений (1) и (4) совпадают. Формально, динамические уравнения (1) и (4) не эквивалентны, но эти уравнения являются в некотором смысле родственными. Во всяком случае, траектории решений уравнения (4) дают представление о динамике концентраций при переходе вещества в стабильное состояние (по модельному уравнению (1)). Математическое же преимущество исследования решений уравнения в форме (4) состоит в том, что правая часть (4) является нелинейным оператором класса Лере - Шаудера, то есть допускающим представление в виде «единица плюс вполне непрерывный >:
—В-1 grad V(В-1 ш) = ш + В-1 д(В-1 ш) = ш + д(ш), д(А) := А А-2 А — В-1 (В-1 А)3 — 3К В-1 (В-1 А)2 + а , а = [[ (В-1 (В-1 А)3 + 3К В-1 (В-1 А)2) —Х1—Х2.
Для такого класса уравнений нетрудно провести обоснование применимости метода конечномерной редукции, основанного на комбинации двух редукций - ритцевской аппроксимации функционала энергии и редукции Пуанкаре - Ляпунова - Шмидта [2, 3].
Прямое обоснование применимости редукции Пуанкаре - Ляпунова - Шмидта к уравнению (1) или (3), по-видимому, можно провести, опираясь на формулу Дюамеля для нелинейных полугрупп (см., например, [4, с. 158]).
Заметим, что в недавней опубликованной работе A.B. Казарникова и C.B. Ре-виной [5] получены формулы асимптотических приближений к бифурцирующему из нуля периодическому решению в системе Релея с диффузией, имеющей более общий вид, чем система (1) данной работы. Получние закритической ветви автоколебаний проведено в ней на основе (невариационной) схемы Ляпунова - Шмидта, ранее предложенной В.И. Юдовичем.
Представленный в данной статье подход допускает применение в задаче, рассмот-репной в работе [5]. Более точное и эффективное приближение вычисления трассы спуска можно получить, используя подходы, разработанные P.A. Свиридюком и его учениками на основе теории полугрупп [6-8].
1. Локальный анализ бифурцирующих стационарных решений уравнения «реакция-диффузия»
Рассмотрим вновь уравнение диффузии с кубической нелинейностью, записанное в форме
гг = —grad (V (т) — С [ад]) = Дт + Хт — г3 — С,
где т = т(х,Ь) - концентрация изучаемого компонента, х = х(х1;х2), х Е и = [0,1] х [0,1] С К2,
V(w) :--
(
\Vw\2 Aw2 w . , —2--A— + — ) dxidx2
-)
интеграл энергии по области U = [0,1] х [0,1] С R2. Будем предполагать выполнение
граничного условия Неймана fr
dU
0
концентрацию вещества в целом (см., например, [1]):
[w] = ii w(xi-x2)dxidx2 = к>0
п
(5)
При исследовании локальных бифуркаций экстремалей часто используется локальная конечномерная ритцевская аппроксимация функционала энергии
Ж(£) := V(К + 6б1 + 6е2 + ... + Сиви), построенная по «главным> собственным функциям (модам) в^ оператора Лапласа на
и
ном функциональном пространстве состояний):
во = 1, в1 = л/2еов(пх1), в2 = л/2ео8(пх2), в3 = 2еов(пх1) еов(пх2)), в4 = л/2еов(2пх1), в5 = л/2еов(2пх2), в6 = 2 еов(2пх1) еов(пх2)),
в7 = 2 еов(пх1) еов(2пх2)), в8 = 2 еов(2пх1) еов(2пх2)), ....
В нелокальных задачах также можно использовать ритцевскую аппроксимацию, но при этом для достижения требуемой точности решения необходимо использовать большое количество мод, что приводит к большой размерности аппроксимирующей системы. Снизить ее размерность можно за счет использования нелинейной ритцев-ской аппроксимации, например, в виде нелокально продолженной ключевой функции. Дальнейший анализ можно осуществить переходом к конечномерной задаче
grad W(С) = 0 , С е Rn
W(С) :=
inf V (w)
(w, ej )=Zj
нелокально продолженная ключевая функция Ляпунова - Шмидта. Запишем искомую функцию концентрации т в виде
w = K + w,
w dxldx2 = 0.
(6)
Сужение функционала V на подпространство функций с фиксированным средним (равным К) осуществляется подстановкой (2), в результате которой получаем V(т) := V(К + на П2-
= Л (^1-+ ^ + Кй3^ ¿х^ , А = Л - 2 К2, (7)
и, соответственно,
-grad(V)(w) = Д^ + А га - V3 - 3К й2 - а, (8)
а = // (й3 + 3К й2) ¿х14х2 (в тройке пространств И2 С Ц0 С Ц0 [2-3]). Для отыс-п
Л
(в точке V = 0 Дй + Л гй = 0. Минимальным критическим значением параметра Л является число п (соответственно Л — Л* = п2 + 8). Модами бифуркации при таком Л ЯВЛЯЮТСЯ функции
е\ = л/2 С08(пХ1) , &2 = л/2 С0в(пХ2).
Теорема 1. При малых значениях концентрации и при малых 8 := Л - Л* главной частью ключевой функции, соответствующей функционалу (3), является многочлен (4-ой степени)
U (6,6) = - 2 (£2 + Û ) + 8(^4 + £) + 2 &
Доказательство. Доказательство вытекает из определяющего ключевую функцию соотношения ^
W(£i,6):= х inf е V(V) ,
И ИЗ того, что
Ws(6,6) = Vfée + ^2) + О (1612 + Ы2)2 =
= - 2 (£2 + £2) + 4(^4 + £4) + l&g + О ш2 + Ы2)2 (9)
(см. [9 - 10]). Здесь учтено, что ff w3dx2dx2 = 0, если V = £2e2 + £2e2. Несложные
п
вычисления показывают, что
a = JJ e2 dx2dx2 = JJ e4 dx2dx2 = 3 , п п
2
Ь = JJ е1е2йх1йх2 = I J е1йх1 I ^ е^х2 | =1
П \0 / \0
Такой вид главной части ключевой функции связан, во-первых, с ее симметрией относительно преобразования (6,6) ^ (6,6)
и, во-вторых, с тем, что кубическая часть в тейлоровском разложении функционала V обнуляется на 6е1 + £2е2. Из теоремы 1 вытекает следующее утверждение.
Теорема 2. Для функционала (3) вблизи нуля при малых концентрациях имеется ветвь устойчивых функций концентраций вида
и = К + е(в! + 62) + о(е),
где е = с 51 - малый параметр при с = .
Таким образом, удалось получить первое асимптотическое слагаемое «(критической ветви бифурцирующих экстремалей - по информации о главной части ключевой функции.
2. Построение нелокальной ключевой функции методом градиентного спуска
При построении нелокальной ключевой функции можно также воспользоваться прямой процедурой кратчайшего спуска (см. [10-13]) в точку минимума V при условии, что направление спуска зависит от параметра. Первым шагом этой процедуры является выбор величины сдвига вдоль градиента из начальной (порождающей) точки, с целью уменьшения значения функционала энергии.
Пусть б\,... ,6П фиксированный базис ритцевской аппроксимации (базис Ритца), составленный из собственных функций оператора Лапласа (в порядке возрастания номеров собственных функций без пропусков отдельных функций), и пусть б1,б2 -основные моды, по которым допускается вырождение. Пусть при этом VRitz(С) :=
V ( Ск бк)- ритцевская аппроксимация функционала энергии (по базису Ритца). \к=1 )
В качестве нулевого приближения к функции
Ж(С) := 1п£ VШz(и) , С =(С1,С2) ,
{т,вз )=(,} ,,з=1,2
рассматривается функция
^с(С) := V(С1 61 + С2 62)
61, 62
Первый шаг заключен в выборе «поправки:» к дающей первое приближение к Ж (С) в виде
Wl(С) := ^(а(С)) ,
"'Л
Ь112
a(£)) = ao - So go , ao = (£1,6, 0,..., 0), go := grad?3 in VrUz (ao), so = ,
\Go go, go)
где Go = hess^...^ VRuz(ao) - матрица Гессе (в нулевой порождающей точке ao) функции Vmtz то перемениым ... ,£n.
a1 = ao -
so go , и т.д. На шаге с номером к + 1 получаем
afc+i = ak - sk gk , gk := g(ak).
Выбор функциональной величины сдвига sk = sk (£) вдоль антиградиента делается посредством формулы
sk = ■ (10)
\Gk gk ,gk)
I on Bulletin of the South Ural State University. Ser. Mathematical Modelling, Programming
& Computer Software (Bulletin SUSU MMCS), 2017, vol. 10, no. 1, pp. 125-137
где Ск = Ьевв^,...,^ Ушх (а0) _ матрица Гессе (в точке ак) функции УКцх по переменным
К сожалению, выбор функциональной величины вк = 8к(С)
сдвига (вдоль антиградиента) в виде (6) приводит к повышенному росту информационного потока, сопровождающего вычисления, и, как следствие, к существенному замедлению работы алгоритма и к быстрому достижению предела возможностей вычислительных устройств. Это препятствие можно преодолеть, заменив функциональный множитель (6) числовым множителем а, служащим универсальной оценкой снизу всех функциональных множителей вида (6). Подбор такого ограничителя снизу можно осуществить, используя следующие легко проверяемые неравенства (для произвольной симметричной матрицы С = (дз,к)):
\\С\\г <№ <\\С\\з,
:= шах^ : » е 8рее(С)} , \\С\\2 := ^г (СТС) , \\С\\3 := ^ \д3,к\.
3, к
Ниже изображены границы кластеров концентрации при соответствующих коэф-
ак
выделением нулевых линий уровня соответствующих промежуточных и финальной функций). Вычисления, проводимые на основании изложенных вьтттте теоретических положений, дают возможность визуализировать критические точки (вычислить сколь угодно точно координаты, определить их характер и расположение), и с любой наперед заданной точностью локализовать притягивающие точки: сначала вычислить точки минимума ключевой функции, а затем по связывающей формуле определить притягивающие точки уравнения (1). Прямая графическая локализация позволяет сколь угодно точно определить координаты точек минимума для ключевой функции (рис. 1 3).
В пределах заданного уровня точности получаем координатное изображение (0, 093; 0, 093) для точки минимума, лежащей на биссектрисе координатного угла. Этой
точке соответствует стабильная концентрация (рис. 4).
Рис. 1. Линии уровня ключевой функции и локализация ее точки минимума при К = 0, 5 А = п2 + 0, 4
Рис. 2. Линии уровня и 3Д-график функции концентрации при K = 0, 5 A = п2 + 0,4
3. Трассировка спуска из произвольного состояния в стабильное
Если т - начальная концентрация, находящаяся в общем положении, то выходящая из нее траектория динамической системы (1) стремится к концентрации, реализующей минимум функционала энергии. С помощью пошагового метода Коши -Пикара из теории обыкновенных дифференциальных уравнений и его варианта для градиентных уравнении (в виде градиентного спуска) можно наблюдать изменения функции концентрации посредством визуального изучения строения линий уровня функции концентрации в «контрольных> точках траектории спуска. Это позволяет получить представление о динамике изменения концентрации.
Итак, пусть задана основная начально-краевая задача
т = Ат + Ат — т3 — С = 0,
т = т(х, ¿) - функция концентрации, х = х(х^х2), х е и = [0,1] х [0,1] С К2.
Рассмотрим собственные функции оператора А = — А :
е0 = 1, е1 = л/2со^(пх1), е2 = л/~2со$,(пх2), е3 = 2сов^х^ сов(пх2)),
е4 = л/2сов(2пх1), е5 = л/2 сов(2пх2), е6 = 2сов(2пх1) сов(пх2)), е7 = 2 сов(пх1) сов(2пх2)), ез = 2 сов(2пх1) сов(пх2)), ед = л/2 сов(3пх1), е10 = л/2 сов(3пх2) .... и соответствующие собственные А
А1 = А1 = п2 , Аз = 2 п2, А4 = А5 = 4 п2 ,
Аб = А7 = 5 п2, Аз = 8 п2, Ад = А10 = 9 п2 ....
Функцию концентрации запишем в виде
т = К + т, II т ¿х1йх2 = 0. (11)
п
Основные итерационные а0
ленности положим К = 0, 5,
а0 = 0, 5 + 0, 3 сов(3пх1) + 0,4 сов(3пх2) + сов(пх1) + 2 сов(пх2). Далее имеем
ап+1 = ап — а /п, и = /(ап)). (12)
Для ритцевской аппроксимации функционала V имеем представление
Уяих (С) = V К +
[к + Ск е^ .
Суть первого тттага процедуры градиентного спуска заключена в выборе сдвига вдоль антиградиента из нач^ьной (порождающей) точки ао в точку, дающую уменьшение значения функционала энергии
т = ао — в/о, /о := /(ао) .
Второй тттаг процедуры градиентного спуска является повторением первого тттага в а1
в точке а2 := а1 — а /1; а2 = а1 — а /1, /1 := /(а1), и так далее.
После к-го шага получаем ак+1 = ак — а /к, /к := /(аи). Последовательность точек ак, к = 1, 2,... называется трассой прямого спуска. Процесс построения трассы называть трассировкой траектории спуска.
ак
трацией (рис. 5).
х[1] <1]
х[1] <1]
Рис. 5. Изменение линий уровня функции концентрации вдоль трассьт спуска (из начальной точки ао = 0, 5 + в\ + 2 e2 + 0, 3 e3 + 0, 4 e4 до третьего шara a3, K = 0, 5, A = 15)
Литература
1. Михлитт, С.Г. Вариационные методы в математической физике ./' С.Г. Михлитт. - М.: Наука, 1970. - 512 с.
2. Даринский, Б.M. Бифуркации экстремалей фредгольмовых функционалов / Б.М. Да-рннскнй, Ю.И. Сапронов, С.Л. Царев // Современная математика. Фундаментальные направления. - 2004. - Т. 12. - С. 3-140.
3. Марри, Дж. Нелинейные дифференциальные уравнения в биологии. Лекции о моделях / Дж. Марри. - М.: Мир, 1983. - 397 с.
4. Хэссард, Б. Теория и приложения бифуркации рождения цикла / Б. Хэссард, Н. Каза-ринов, И. Вэн. - М.: Мир, 1985. - 280 с.
5. Казарников, A.B. Возникновение автоколебаний в системе Рэлея с диффузией / A.B. Ка-зарнпков, C.B. Ревпна // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2016. - Т. 9, № 2. - С. 16-28.
6. Свиридюк, Г.А. Квазистационарные траектории полулинейных динамических уравнений типа Соболева / Г.А. Свиридюк // Известия АН СССР. Серия: Математическая. -1993. - Т. 57, № 3. - С. 192-207.
7. Свиридюк, Г.А. О задаче Веригина для обобщенного фильтрационного уравнения Бус-синеска / Г.А. Свиридюк, С.А. Загребина // Известия вузов. Математика. - 2003. -№ 7. - С. 54-58.
8. Загребина, С.А. О задаче Шоуолтера - Сидорова / С.А. Загребина // Известия вузов. Математика. - 2007. - № 3. - С. 22-28.
9. Костина, Т.И. Нелокальное вычисление ключевых функций в задаче о периодических решениях вариационных уравнений / Т.И. Костина // Вестник Воронежского государственного университета. Серия: Физика. Математика. - 2011. - № 1. - С. 181-186.
10. Сапронов, Ю.И. Моделирование диффузорных течений жидкости посредством редуцированных уравнений / Ю.И. Сапронов // Вестник ЮУрГУ. Серия: Математическое моделирование и программирование. - 2014. - Т. 7, № 2. - С. 74-86.
11. Красноселький, М.А. Итерационный процесс с минимальными невязками / М.А. Крас-носелький, С.Г. Крейн // Математический сборник. - 1952. - Т. 31 (73), № 2. - С. 315-334.
12. Лемешко, A.A. О равномерной сходимости с производными галеркинских приближений к решениям уравнений с параметрами / A.A. Лемешко // Математические модели и операторные уравнения. - 2003. - Т. 2. - С. 94-103.
13. Лемешко, A.A. О равномерной сходимости ньютоновских приближений к решениям уравнений с параметрами / A.A. Лемешко // Сборник трудов молодых ученых математического факультета Воронежского государственного университета. - 2003. - С. 74-83.
14. Ковалева, М.И. Огибающие кривые, точки возврата и бифуркационный анализ нелинейных задач // М.И. Ковалева, Т.И. Костина, Ю.И. Сапронов. - Воронеж: Военный учебно-научный центр Военно-воздушных сил «Военно-воздушная академия имени профессора Н.Е. Жуковского и Ю.А. Гагарина:», 2015. - 242 с.
Андрей Сергеевич Коротких, аспирант, кафедра математического моделирования, Воронежский государственный университет (г. Воронеж, Госсийская Федерация), [email protected].
Поступила в редакцию 20 сентября 2016 г.
MSC 90C30, 90C90 DOI: 10.14529/mmp 170108
STATIONARY POINTS OF THE "REACTION-DIFFUSION" EQUATION AND TRANSITIONS TO STABLE STATES
A.S. Korotkih, Voronezh State University, Voronezh, Russian Federation,
Of concern is a stationary " reaction-diffusion" equation with cubic non-linearity is Neumann boundary conditions and fixed average value of the desired bifurcating solutions. A method of approximate calculation of bifurca-ting solutions for small and finite values of supercritical parameter increment are presented. Computing is based on the Lyapunov-Schmidt reducing procedure and is leaning on key functions Ritz' approximation of the set of eigenfunctions (modes) of main linear part of gradient energy functional. A technique of evaluating of a functional space size, where Lyapunov - Schmidt reduction can be applied is performed. In case of local reduction the main part of the key function has been found and asymptotic presentation of bifurcating solutions for small supercritical increment of bifurcation parameter is calculated. The relation between solutions search procedures for "reaction-diffusion" equations and Cahn - Hilliard equation (with extended Neumann boundary conditions) is also performed. Graphs are presented.
Keywords: continuously differentia,ble functional; extremal; bifurcation; Lyapunov -Shmidt method.
References
1. Mikhlin S.G. Variational Methods in Mathematical Physics. Oxford, London, Edinburgh, N.Y., Paris, Frankfurt, Pergamon Press, 1964. 584 p. DOI: 10.1002/zamm. 19650450439
2. Darinskii B.M., Sapronov Y.I., Tsarev S.L. Bifurcations of Extremals Fredholm Functional. Journal of Mathematical Science, 2007, vol. 145, issue 6, pp. 5311-5453. DOI: 10.1007/sl0958-007-0356-2
3. Murray J.D. Lectures on Nonlinear Differential-Equation. Models in Biology About the Models. Oxford, Clarendon Press, 1977.
4. Hassard B.D., Kazarinoff N.D., Wan Y.-H. Theory and Applications of Hopf Bifurcation. Cambridge, Cambridge University Press, 1981. 320 p.
5. Kazarnikov A.V., Revina S.V. The Onset of Auto-Oscillations in Raileigh Systems with Diffusion. Bulletin of South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2016, vol. 9, no. 2, pp. 16-28. (in Russian) DOI: 10.14529/mmpl60202
6. Sviridyuk G.A. Quasistationary Trajectories of Semilinear Dynamical Equations of Sobolev Type. Russian Academy of Sciences. Lzvestiya Mathematics, 1994, vol. 42, no. 3, pp. 601-614. (in Russian) DOI: 10.1070/IM1994v042n03ABEH001547
7. Sviridyuk G.A., Zagrebina S.A. On the Verigin Problem for the Generalized Boussinesq Filtration Equation. Russian Mathematics (Lzvestiya VUZ. Matematika), 2003, vol. 47, no. 7, pp. 55-59. (in Russian)
8. Zagrebina S.A. On the Showalter - Sidorov Problem. Russian Mathematics (Lzvestiya VUZ. Matematika), 2007, vol. 51, no. 3, pp. 19-24. (in Russian) DOI: 10.3103/S1066369X07030036
9. Kostina T.I. [Nonlocal Calculation of Key Functions in the Problem of Periodic Solutions of Variational Equations]. Proceeding of Voronezh State University. Series: Physics. Mathematics, 2011, no. 1, pp. 181-186. (in Russian)
10. Sapronov Yu.I. Modelling Fiquid Flows in Diffusers by Reduced. Bulletin of South Ural State University. Series: Mathematical Modelling, Programming and Computer Software, 2014, vol. 7, no. 2, pp. 74-86. (in Russian) DOI: 10.14529/mmpl40207
11. Krasnosel'skij M.A., Krejn S.G. [The Iterative Process Minimal Residuals]. Sbornik: Mathematics, 1952, vol. 31 (73), no. 2, pp. 315-334. (in Russian)
12. Femeshko A.A. [Uniform Convergence with Derivatives Galerkin Approximations to the Solutions of Equations with Parameters]. Matematicheskij rnodeli i operatornye uravnenija [Mathematical Models and Operator Equations], 2003, vol. 2, pp. 94-103. (in Russian)
13. Femeshko A.A. [Uniform Convergence of Newton's Approximations to Solutions of Equations with Parameters]. Sbornik trudov molodyh uchenyh matematicheskogo fakul'teta Voronezhskogo gosudarstvennogo Universiteta [Proceedings of Young Scientists Mathematical Faculty of Voronezh State University], 2003, pp. 74-83. (in Russian)
14. Kovaleva M.I., Kostina T.I., Sapronov Y.I. Ogibayushchie krivye, tochki vozvrata i bifurkatsionnyy analiz nelineynykh zadach [The Envelope Curve, the Point of Return and Bifurcation Analysis of Nonlinear Problems]. Voronezh, Military Training and Research Center of the Air Force "Air Force Academy named after Professor N.E. Zhukovsky and Y.A. Gagarin", 2015. 242 p.
Received September 20, 2016
