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

Прикладные задачи

нелинейной теории колебании и волн

УДК 517:957 https://doi.org/10.18500/0869-6632-2019-27-4-85-98

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

Ю. А. Хазова, В. А. Лукьяненко

ФГАОУ ВО «Крымский федеральный университет им. В.И. Вернадского» Россия, Республика Крым, 295007 Симферополь, просп. Академика Вернадского, д. 4 E-mail: yuliya@hotmail.com, art-inf@yandex.ru Автор для переписки Юлия Александровна Хазова, yuliya@hotmail.com Поступила в редакцию 20.03.2019, принята к публикации 2.07.2019

Тема. Работа посвящена изучению математической модели, описывающей оптическую систему с двумерной обратной связью. Примером такой оптической системы может быть нелинейный интерферометр с зеркальным отражением поля. Математической моделью выступает нелинейное функционально-дифференциальное параболическое уравнение с преобразованием отражения пространственной переменной и условиями на круге. Цель работы состоит в исследовании условий возникновения пространственно неоднородных стационарных решений. Предполагается дать ответ на вопрос об асимптотической форме рождающихся решений и определении их устойчивости. Методы. Исследование проводится методами теоретического анализа, а именно используются метод центральных многообразий, метод Фурье и метод сведения к интегральному уравнению. Используя метод разделения переменных, доказана лемма о собственных значениях и собственных функциях соответствующей линеаризованной задачи. Для определения асимптотической формы решения для линеаризованной и соответствующей нелинейной параболической задачи применялся метод сведения к интегральному уравнению. Приведены необходимые выкладки по доказательству единственности и непрерывной зависимости решения от начальных условий. Результаты и обсуждение. На основе метода центральных многообразий доказана теорема о существовании и устойчивости пространственно неоднородного стационарного решения. Получено представление для неоднородной структуры, рождающейся в результате бифуркации типа «вилка» при переходе бифуркационного параметра через критическое значение. Согласно с доказанной теоремой это решение рождается асимптотически устойчивым. Эта теорема носит локальный характер и работает в окрестности бифуркационного значения коэффициента диффузии. Результаты, приведенные в данной работе, являются продолжением исследований в области нелинейной оптики. Возможность анализа рождающихся структур не только в окрестности бифуркационного значения параметра, но и на всем промежутке изменения выбранного параметра, остается актуальной. Представленные в работе результаты могут быть применены как в теоретическом анализе задач нелинейной оптики, так и в практической обработке и интерпретации информации, полученной при постановке вычислительных или физических экспериментов.

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

Образец цитирования: Хазова Ю.А., Лукьяненко В.А. Применение интегральных методов для исследования одной параболической задачи//Известия вузов. ПНД. 2019. T. 27, № 4. С. 85-98. https://doi.org/10.18500/0869-6632-2019-27-4-85-98

Финансовая поддержка. Исследование выполнено при поддержке Программы развития федерального государственного автономного образовательного учреждения высшего образования «Крымский федеральный университет имени В.И. Вернадского» на 2015-2024 годы по проекту «Сеть академической мобильности «Академическая мобильность молодых ученых России» в 2017 году на базе ФГБУН Институт математики имени С.Л. Соболева СО РАН, лаборатория дифференциальных и разностных уравнений.

https://doi.org/10.18500/0869-6632-2019-27-4-85-98 Application of integral methods for the study of the parabolic problem

Yu.A. Khazova, V.A. Lukianenko

V.I. Vernadsky Crimean Federal University 4 Vernadskogo Prospekt, Simferopol 295007, Republic of Crimea, Russia E-mail: yuliya@hotmail.com, art-inf@yandex.ru Correspondence should be addressed to Yu. A. Khazova, yuliya@hotmail.com Received 20.03.2019, accepted for publication 2.07.2019

Topic. The work is devoted to the study of a mathematical model describing an optical system with two-dimensional feedback. An example of such an optical system could be a nonlinear interferometer with a specular reflection of the field. The mathematical model is a nonlinear functional differential parabolic equation with the transformation of the spatial variable reflection and conditions on the disk. Aim of the work is to study the conditions for the occurrence of spatially inhomogeneous stationary solutions. It is supposed to answer the question of the asymptotic form of the solutions being born and the determination of their stability. Methods. The study is carried out by methods of theoretical analysis, namely, the method of central manifolds, the Fourier method and the method of reducing to an integral equation are used. Using the method of separation of variables, a lemma on eigenvalues and eigenfunctions of the corresponding linearized problem is proved. To determine the asymptotic form of the solution for the linearized and corresponding nonlinear parabolic problems, the method of reduction to an integral equation was used. The necessary calculations on the proof of uniqueness and the continuous dependence of the solution on the initial conditions are given. Results and discussion. Based on the method of central manifolds, a theorem on the existence and stability of a spatially inhomogeneous stationary solution is proved. A representation is obtained for an inhomogeneous structure that is born as a result of a bifurcation of the «fork» type when the bifurcation parameter passes through a critical value. According to the proved theorem, this solution is born asymptotically stable. This theorem is local in nature and works in the vicinity of the bifurcation value of the diffusion coefficient. The results presented in this paper are a continuation of research in the field of nonlinear optics. The possibility of analyzing the structures being born not only in the vicinity of the bifurcation parameter value, but also throughout the change interval of the selected parameter, remains relevant. The results presented in this paper can be applied both in the theoretical analysis of the problems of nonlinear optics and in the practical processing and interpretation of information obtained in the formulation of computational or physical experiments.

Key words: parabolic equation, transformation of the reflection, integral operator.

Reference: Khazova Yu.A., Lukianenko V.A. Application of integral methods for the study of the parabolic problem. Izvestiya VUZ. Applied Nonlinear Dynamics, 2019, vol. 27, no. 4, pp. 85-98. https://doi.org/10.18500/0869-6632-2019-27-4-85-98

Acknowledgements. The research was carried out with the support of the Academic Mobility Network «Academic Mobility of Young Scientists of Russia» by V.I. Vernadsky Crimean Federal University on the basis by Sobolev Institute of Mathematics SB RAS, Laboratory of Differential and Difference Equations in 2017 year.

Введение

Наблюдаемый в последнее время интерес к нелинейным оптическим системам с пространственными преобразованиями в двумерной обратной связи [1,2] вызван использованием этих систем в задачах оптической обработки информации. Богатый спектр локальных и нелокальных преобразований позволяет генерировать и визуально наблюдать процессы самоорганизации: пульсирующие волны, диссипативные структуры, ведущие центры, оптическую турбулентность. Такие структуры можно использовать для кодирования и хранения информации. Возникновение световых структур в оптических системах с нелокальной обратной связью является следствием потери устойчивости некоторого пространственно-однородного светового режима в результате внесения возмущений в условия протекания процесса. Математической моделью этих систем является квазилинейное параболическое уравнение с преобразованием пространственных переменных [3]. Анализ таких систем нелинейной оптики возможен как в рамках метода квазинормальных форм [4], так и методом малого параметра, методом центрального многообразия [5] или методом Галеркина [6].

В работах [7-10] исследовалось нелинейное параболическое уравнение с преобразованием отражения на окружности и отрезке. В данной статье в качестве области, на которой рассматривается уравнение, выбран круг.

и периодичности

ограниченности по r в нуле

и начальным условием

1. Нелинейное параболическое уравнение

1.1. Постановка задачи. Рассматривается смешанная краевая задача для нелинейного параболического уравнения в круге

ut + u = DAu + K

с преобразованием отражения пространственной переменной

Qu

с условиями Неймана

5u%M = о

dr

u

u

Здесь Au = 1 dr (rfr) + дфи - оператор Лапласа в полярной системе координат. Задача (1)-(6) моделирует динамику фазовой модуляции u(r, ф, t), световой волны, прошедшей тонкий слой нелинейной среды керровского типа с преобразованием отражения координат в двумерной обратной связи. Здесь D - коэффициент диффузии нелинейной среды, положительный коэффициент K пропорционален интенсивности светового потока, Y - видность (контрастность) интерференционной картины, 0 < Y < 1.

В данной работе рассматриваются вопросы существования, формы и устойчивости в метрике H1 пространственно неоднородных стационарных решений, бифурцирующих из пространственно однородных стационарных решений, то есть решений u(r, ф, t) = w определяемых из уравнения

w = K

С ростом K количество сосуществующих корней этого уравнения неограниченно растет, причем при K ^ то их состав постоянно меняется: появляются новые и исчезают старые. Поэтому фиксируем гладкую ветвь решений уравнения

w = w

1.2. Линеаризация уравнения в окрестности стационарного решения. Линеаризуем

vt + v = DAv — Ky sin w Qv + S(Qv), где S(v) = Ky((cos w(cos v — 1) — sin w(sin v — v)).

Таким образом, уравнение для V с выделенным линейным оператором Ь имеет вид

vt + Lv = K y((cos w(cos Qv — 1) — sin w(sin Qv — Qv)), Lv = v — DAv + Ky sin w Qv.

Оператор L представим в виде суммы двух операторов

Lv = (Lo + Li)v, Lo = I — DA, Li = —ЛQ, Л = Л(К, у) = —Ky sin w,

где оператор отражения Q определен согласно равенству Qv(r, ф, t) = v(r, п — ф, t) и обладает свойством Q2 = I.

Используя разложение для функций cos v и sin v, получаем

vt + v = DAv + К yQ

v2 v4

COs w[ — 21 + 4! — "J — sin wT — 3! + 5! —

v3 v5

что позволяет записать варианты моделей с выделенной линейной частью: квадратичное уравнение vt = DAv — v — Куsinw Qv — Ky^OfpQv2; кубическое уравнение vt = DAv — v— —К y sin wQv — Ky ^Ofp Qv2 + ky fifp Qv3; а также уравнение с выделенной линеаризованной частью и слагаемым, содержащим разложение v по нечетным степеням или по четным.

Заметим, что условие cos w = 0 приводит к разложению нелинейного оператора только по нечетным степеням. В этом случае первое условие в ( ) имеет вид w = K и не зависит от y, K = w = п + пп = ^п, п = 0,1, 2 ..., sin (^п) = (—1)n, п = 0,1, 2 ... А второе условие в ( ) 1 + Ky sin w = 1 + Ky(—1)n = 0 для нечетных n (ky = 1) и для четных (Ky = —1) выполняется автоматически, так как K> 0, 0 < y < 1.

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

vt + Lv = vt + v — DAv + Ky sin w Qv = 0,

где v = v(r, ф, t) - искомая функция, 0 < r < ri, 0 < ф < 2п, t > 0, Av = 1 (rf^) + Г2-оператор Лапласа в полярной системе координат.

В соответствии с ( )-( ) для v получаем условия Неймана на границе при r = ri: = 0; условие периодичности v(r, ф, t) = v(r, 2п + ф, t); условие ограниченности в начале координат |v(0, ф, t)| < c < то и начальное условие v(r, ф, 0) = u0(r, ф) — w = v0.

Решение (10) будем строить, используя метод разделения переменных

v = X(r, ф)Т(t),

XT' + XT = DAXT — Ky sin w  QXT,

dX (r1, ф) dr

= 0, X(r, ф) = X(r, 2п + ф).

Для нетривиальных решений имеем

T' lo. nAX ~ • QX . T = —1 + dx — K y sin w "X = —

Тогда для функции X получаем задачу Штурма-Лиувилля

' DAX - Kу sin w  QX + (-1 + X)X = 0,

= 0,

dX(n, ф) dr

X(r, ф) = X(r, 2п + ф), [ |X(0, ф)| < то,

а для функций T:

T'(t) + XT (t) = 0.

Лемма 1. Задача Штурма - Лиувилля ( ) имеет собственные функции вида

Xfcm(r, ф) = (Jfc(A,£mr) cos кф, Jk(^kmr) sin кф}, которым соответствуют собственные значения

Km = D

Xfcm = D

Mkm

ri

+ (-1)kKy sin w + 1,

2 + (-l)k+1 Ky sin w + 1,

к = 0,1,2,..., т = 1,2,..., где Зк (ж) - функция Бесселя, цкт - корни уравнения

•к

я(г)Ф(ф), в результате

получаем

D

1 д 1

___

- Ky sin wR(r)^fa) + (X - 1)д(г)Ф(ф) = 0,

Ф(ф) = Ф(2п + ф), R'(ri) = 0, |R(0)| < то.

В представлении

r2R''(r) + rR'(r) Ky sin w 2 £Ф(ф) (-1+ X)r2 _ Ф''(ф) _

R

D

Ф

D

Ф

выражение является постоянным на собственных функциях Фк(ф) задачи Штурма-Лиувилля

( Ф''(ф) + уф(ф) = 0, \ Ф(ф) = Ф(2п + ф),

Фк (ф) = {cos кф, sin кф}, vk = к2, к = 0,1____

Оператор Q на системе функций Фк(ф) принимает следующие значения

Q cos кф = cos к(п — ф) = (—1)k cos кф, к = 0,1..., Q sin кф = sin к(п — ф) = (—1)fc+1 sin кф, к = 0,1...

2

Уравнение Бесселя для функций R(r) также принимает различные значения на системе функций {cos кф, sin кф}

r2R''(r) + rR'(r) + (Ir2 - Л2) R(r) = 0, к = 0,1...,

, Г (-1+i)-(-1)kKY sin w  , Г (-1+^)-(-1)fc+1 K Y sin w

где для cos кф Л = ^— D—-и для sin кф Л = ^— D--•

В уравнении Бесселя сделаем замену Vi r = ж, R(r) = Rix /Л) = y (ж). Получаем классическую форму уравнения Бесселя

ж2у''(ж) + жу'(ж) + (ж2 - к2)у(ж) =0, к = 0,1...

Из условия ограниченности решения R(r) в нуле следует выбрать ограниченные решения уравнения Бесселя у(ж) = Jk(ж), к = 0,1... Условие Неймана на границе приводит к соотношению R'(r1) = у' (л/ir^ = 0, J' (л/ir^ =0, к = 0,1... Для каждого к уравнение ( ) имеет счетное число корней : J') =0, к = 0,1..., m = 1,2 ...

Для собственных значений Л получаем л/Лг1 = или ikm = (^kp) , к = 0,1..., m = 1, 2 ... Тогда для функций cos кф находим

~ / 4 2 Dikm = -1 + ikm - (-1)'KY sin W = D f ^

ikm = d(2 + (-1)'KY sin W + 1,

для sin кф

Dikm = -1 + i'm - (-1)k+1KY sin W = D (^Tf)' ) 2 + (-1)k+1KY sin w + 1.

is _ n ¡ rkm Afcm = D

Для ri = 1 собственные значения X задачи Штурма-Лиувилля ( ) на функциях cos кф равны Xkm = + (-1)kKy sin w + 1, а для sin кф: = + (-1)fc+1KY sin w + 1.

Им отвечают собственные функции

Xfcm(r, ф) = (Jfc(Xkmr)cos кф, Jfc(Xf,mr)sin кф}, к = 0,1,2,..., m = 1, 2,.... (15) Лемма доказана.

2. Нелинейное интегральное уравнение

2.1. Представление решения линеаризованной задачи. Решение задачи (10) предста-вимо в виде ряда по функциям (15), то есть

те те

v

fc=Q m=1

где X£m, A|m определяются формулами (13). Коэффициенты cfcm, dfc m определяются из начальных условий.

Перепишем (16) в виде

a (r 0) ^

v

k=i

где

a0(r 0) = ¿ Come-X0mtJo(XOmr), Ofc(r, t) = J] Ckme-^ Jk(X^),

m= 1 m= 1

те

bk(r,t) = ^ dfcme-XkmiJfc(XjLr).

m=1

Если ряд (17) сходится, то из начального условия следует

a (r 0) ^

q0(r, ф) — w = —2--+ ^(ak(r, 0) cos кф + bk(r, 0) sin кф).

k=1

Отсюда

2п

1

ak(r, 0) = — i[q0(r, ф) — w] cos кф ^ф, к = 0,1,2,... п У o

2п

bk(r, 0) = — [q0(r, ф) — w] sin кф ^ф, к = 1,2,____

п

0

Очевидно, коэффициенты Ckm, dkm определяются разложениями

a0

те

ak(r, 0) = ^ CkmJk(Xkmr), bk(r, 0) = ^ dkmJk(^ImO, -^- = C0m J0(X0mr).

m=1 m=1 m=1

Откуда следует, что

rr 1 !• ri

/ ak(r, 0) Jk(Xkmr)r dr / bk(r, 0)Jk(A|mOrdr

c = ^0_ d = _

ckm í-ri , dkm í-ri .

/ Jfc(Xfcmr)rdr / Jfc(Xkmr)rdr

00

2.2. Нелинейное интегральное уравнение. Полученное разложение используем для сведения решения задачи для исходного нелинейного уравнения (9) к интегральному уравнению типа Вольтерры II рода. Для этого правую часть уравнения — + = / разложим в ряд по собственным функциям

те те

/(Г Ф^) = ^ X) [/кт(^)Зк(Хктг) СОв кф + /кт(*)Зк(А|тг) 8Ш кф] .

к=0 т=1

те

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

т\ т\ 2п

Ят^У З2(^ктг)г^г = П///(Г Ф, Зк (^ктг) г СОЙ кФ ^Ф ^

00

r1 r1 2п

/fcm(t^ Jfc(Xkmr)rdr = П/ j / (r, (Xkmr) r sin кф ¿ф dr

п

0 0 0

r1 r1 2п

/0)m(tW T2n c m\mjm _ 1

J02(X0mr)rdr = П/j / (r, ф^Ш^О r ¿ф dr.

2 1 0 v 0m / _

п

0 0 0

Для определения функций Т(¿), в отличии от (12), получаем неоднородные дифференциальные уравнения

Ткт№ + Хкт Ткт(^) = /кт(^), Ткт(^) + ХктТкт = /кт(^).

Отсюда следует, что

t

,-Xf t ^ -Xc

TCm(t) = Ckme-^ + | e-Xkm(t-T)/kCm(T) dT, 0

t

TSm(t) = dkme-Xkmt + | e-Xkm(t-T)/kSm(x) dT.

0

Получим представление, в котором v0

Ckrne-Akmt + J e-Xkm(t-T)/km(T) ¿T Jk(X^r) cos кф+ 0

+ ^dkme-^ + ^ e-Xkm(t-T)/km(T) d^ Jk(X^r) sin кф v(r, ф,t) = v0(r, ф,t)+ / c(/)+ / s(/),

v(r, ф, t) =

k=0 m=1

где

t r1 2п

Ic,s(/) = П/// (/(P,0, T)) EE e-XC;m(t-T)Jk(XCk'mP)cos кeJk(XCk'mr)p cos кф ¿ф dp dT.

п

000 k=0 m=1

Обозначим ядра интегральных операторов /с, /5 через

тете

л c,s

Е Е e-Xkm(t-T)Jk(Xck',mr) Jk(Х££P)P cos кф cos ке

кс"^ — T, ф, e,r, p) = k=0 m=1-^r^-.

/ Jfc(xk',m)rdr 0

0

Подставляя в г

г

Интегральное уравнение для Ог будет иметь вид

Ог = Ого + (Ог) + О/(Ог),

Ь п 2п

О/с5(Ог) = П//(Ог)Ле(£ - т,Оф, е,г, р) ^ффйт,

ооо

Ь г1 2п

О/55(Ог) = !// ! 5 (Ог)к^ - т,Оф, е,г, р) ^ф ф йт. ооо

Ядро кс(£ — т, Оф, е, г, р) содержит под знаком суммы (—1)к, а — т, Оф, е, г, р) содержит (—1)к+1, то есть являются знакопеременными.

2.3. Корректность начально-краевой задачи. При исследовании существования и единственности начально краевой задачи

г4 + ¿ог = / (Ог), г(0) = го, (/ = ¿1 + 5).

Доказательство существования и единственности начально-краевой задачи функционально-дифференциального уравнения диффузии проводится по следующей схеме. Исходная задача записывается в виде задачи Коши для нелинейного операторного функционально-дифференциального уравнения в соответствующем функциональном пространстве с учетом граничных условий. Далее используются результаты по разрешимости операторных дифференциальных уравнений [11,12]. Сводим задачу к операторному уравнению, для которого выполняются условия теоремы Лерэ-Шаудера о разрешимости нелинейных операторных уравнений ([13], с. 84-85). Эти результаты при фиксированных параметрах позволяют установить единственность и непрерывную зависимость решения от начальных условий. Приведем соответствующие выкладки.

Используя свойство полугруппы, порожденной уравнением ( ) на  = {0 < г < Г1, 0 < ф < 2п} € Е2, уравнение рассматриваем при условии ( ). обозначим через Н = £2(П) гильбертово пространство со скалярным произведением (и, г) = иг^ж, Нг(П), I € пространство Соболева измеримых на  функций со скалярным произведением

(и,г)г = ^ / даи(ж)даг(ж)^ж, о<|а|<г ^

где а = (а1, а2) € , да = д"1 д£2, |а| = а1 + а2. Обозначим через Н1 = Н1 () П =

= ди\г=п =0, I € с нормой ||  ||г. Пусть Н-1 - пространство, сопряженное Н1; норма определяется равенством ||и||н-1 = зир{(и, г)/||г||Н1, г € Н1, г = 0}. Имеет место вложение Н1 С Н С Н-1 (по теореме Соболева), вложение Н1 С Н вполне непрерывно.

Пусть В - банахово пространство. Обозначим через С (В) банахово пространство непрерывных и ограниченных по вещественной оси функций со значениями в пространстве В и нормой ||/||с(в) = зир;ед||/(¿)||б• Введем пространство и = С(Н) Р|М2(Н-1) с нормой II/||и = ||/||с(я) + ||/||м2(Н1), где М2(В) банахово пространство измеримых функций / : Я ^ В с нормой

1

2

II/IIm2(B) = «телу ||/(t + s)||Bds 0

Отображение F : v ^ cosQv = Ф(ф^) из U в M2(H-1) дифференцируемо по Фреше. Действительно, Ф € C2(R) - вещественнозначная функция и |Ф(га)(^)| < M, n = 1, 2 для всех ^ € R. Из [14], неравенства Гельдера и теоремы Соболева о вложениях следует, что

Так как

1 1

0

Ф следует:

+ v)) - ф^и) - ф'^И =

1

1(Ф (д(и + ей)) - Ф(ди) - Ф (ои)^)^ 0

то есть получаем ||Ф(ф(и + V)) — Ф(^и) — Ф'(фи)^||м2(я-!) < С11V|и• Следовательно, дифференциал Т' задается формулой Т'(и)-и = Ф'(^и)^ и имеют место оценки

где постоянная С не зависит от и, V, ад.

Таким образом, Т : V ^ Ф(0>) из Н1 в Н-1 является дифференцируемым по Фреше в каждой точке пространства Н1. Значит, для Т > 0 из [14] получаем следующий результат.

Утверждение 1. Задача ( )-( ) с начальным условием ( ) имеет единственное решение и, принадлежащее классу £те([0, Т], Н) П £2([0, Т], Н1).

Докажем справедливость оценок для решений, которые получены на основании следующих ниже выкладок. Умножим (1) в Н на и. Используя интегрирование по частям и условие (6), получаем

Q

Далее, из неравенства

1 i 1 1 -K (1 + Ycos Qu)^ — -(K(1 + y))2 + -||u||2 n J 2 2

Q

следует dt|u|2 + ||u||2 — (K(1 + y))2. Откуда получаем

0

Следовательно, решение задачи (1)-(5) с начальным условием (6) продолжимо на положительную ось.

Утверждение 2. Полугруппа {St} порождается задачей (1)-(6) и обладает следующими свойствами:

1) {St}(H, H) - равномерно непрерывна;

2) {St}(H, H) - равномерно ограничена при t > 0;

3) полугруппа {St} имеет поглощающее множество, ограниченное в H.

Действительно, если равенство

dt

умножить на u — v в H, то получим

Откуда

Здесь а-1 - значение нормы оператора О. Непрерывность {5Ь} из Н в Н следует из неравенства ||и(£) — г(£)||2 < ||и(0) — г(0)||2е2(аКу-1)Ь. Справедливость 2) следует из (20). В качестве поглощающего множества {5Ь} можно взять шар {и € Н : ||и|| < К(1 + у) + 1}.

2.4. Устойчивость структуры решения. Для анализа устойчивости структуры решения

с и те кт и лкт

в зависимости от параметра D необходимо оценить собственные значения A|,m и X|m

Xkm = D (^ 2 + ( —1)kЛ + 1, 4m = D (2 + ( —1)k+1A + 1,

где w есть корень уравнения w = K(1 + у cos w) при условии, что 1 + Ky sin w(K, y) = 0, а M-fcm - корни уравнения J(^fcm) = 0, k = 0,1..., m = 1, 2 ....

Как правило, фиксируется K и выбирается D ^ 1 так, чтобы при выполнении условия 1: Л = Л(К, у) < —1, Л = —Ky sin w решение было асимптотически устойчивым. Например, k = 1. Тогда

X1m = D (2 — Ky sin w + 1, X?m = D (^ 2 — (—KY sin w) + 1,

А,Цт может менять знак при уменьшении D. Критические собственные значения

m = D (М^) 2 + ( —1)21+1Л + 1, ^m = D (^Tf) 2 + ( —1)21 Л + 1, которым отвечают критические собственные функции

J21+1 (Ья+1,тг) cos(21 + 1)ф, J21 (X^r) sin(2^, l = 0,1, 2, . . . , при изменении параметра D меняют знак на противоположный.

Теорема 1. Существует So > 0 такое, что, если 0 < D — Di < So, то уравнение имеет два асимптотически устойчивых решения:

/D — Di\ 1/2 (r, Ф, D) ^ ) J Ji(X1ir)cos ф+

+ 2 ("dcdd^) Лcigw((xio — 2Xii)-1 + (X12 — 2A.ii)-1 ^ф^^ ± К^)3/2 (Ai3 — 3Aii)-^Л — 4Л2ctg2 w(Ai2 — 2Aii)-^ x

x J2(Aiir)J3(Aiir) cos 3ф, где ci(D) = [Л — 4(Acig w)2 ((Aio — 2Aii)-i + 2(A^ — 2Aii)-i)] J2(A^r) < 0.

Необходимые выкладки для доказательства теоремы можно найти в работе [5].

Заключение

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

Библиографический список

1. Ахманов С.А., Воронцов М.А., Иванов В.Ю. Генерация структур в оптических системах с двумерной обратной связью: на пути к созданию нелинейно-оптических аналогов нейронных сетей // в кн. Новые физические принципы оптической обработки информации. М.: Наука, 1990. С. 263-325.

2. Разгулин А.В., Романенко Т.Е. Вращающиеся волны в параболическом функционально-дифференциальном уравнении с поворотом пространственного аргумента и запаздыванием // Ж. вычисл. матем. и матем. физ. 2013. № 53. C. 1804-1821.

3. Budzinskiy S.S., Larichev V.A., Razgulin V.A. Reducing dimensionality to model 2D rotating and standing waves in a delayed nonlinear optical system with thin annulus aperture // Nonlinear Anal.-Real World Appl. 2018. No. 44. P. 559-572.

4. Колесов А.Ю., Розов Н.Х. Оптическая буферность и механизмы ее возникновения // ТМФ. 2004. Т. 140, № 1. C. 14-28.

5. Хазова Ю.А., Шиян О.В. Теорема о существовании и устойчивости решения одного параболического уравнения // Динамические системы. 2018. Т. 8(36), № 3. C. 275-280.

6. Хазова Ю.А. Решение типа «бегущие волны» в параболической задаче с преобразованием поворота // Известия вузов. Прикладная нелинейная динамика. 2017. Т. 25, № 6. C. 57-69.

7. Белан Е.П., Хазова Ю.А. Динамика стационарных структур в параболической задаче на окружности с отражением пространственной переменной //Динамические системы. 2014. Т. 4(32), № 1-2. C. 43-57.

8. Хазова Ю.А. Динамика стационарных структур в параболической задаче на отрезке с отражением пространственной переменной // Динамические системы. 2014. Т. 4(32), № 3-4. C. 245-257.

9. Хазова Ю.А. Стационарные структуры в параболической задаче с отражением пространственной переменной // Таврический вестник информатики и математики. 2015. № 3(28). C. 82-95.

10. Хазова Ю.А. Интегральное представление приближенных решений параболического уравнения // Актуальные направления научных исследований XXI века: Теория и практика. 2018. Т. 6, № 6 (42). C. 384-386.

11. Лионс Ж.-Л., Мадженес Э. Неоднородные граничные задачи и их приложения. М.: Мир, 1971.

12. Temam R. Infinite-dimensional dynamical systems in mechanics and physics. Springer-Verlag, New York, 1997.

13. Лерэ Ж., Шаудер Ю. Топология и функциональные уравнения //Успехи математических наук. 1946. Т. 1, № 3-4. C. 71-95.

14. Бабин А.В., Вишик М.И. Аттракторы эволюционных уравнений. М.: Наука, 1989.

References

1. Akhmanov S.A., Vorontsov M.A., Ivanov V.Yu. Controlling transverse-wave interactions in nonlinear optics generations of spatiotemporal structures. J. Optical Soc. Amer. Ser. B, 1992, vol. 9, no. 1, pp. 78-90.

2. Razgulin A.V., Romanenko T.E. Rotating waves in parabolic functional differential equations with rotation of spatial argument and time delay. Computational Mathematics and Mathematical Physics, 2013, vol. 53, pp. 1804-1821 (in Russian).

3. Budzinskiy S.S., Larichev V.A., Razgulin V.A. Reducing dimensionality to model 2D rotating and standing waves in a delayed nonlinear optical system with thin annulus aperture. Nonlinear Anal.-Real World Appl., 2018, no. 44, pp. 559-572.

4. Kolesov A.Yu., Rozov N.Kh. Optical buffering and mechanisms for its occurrence. Theoret. and Math. Phys., 2004, vol. 140, no. 1, pp . 14-28 (in Russian).

5. Khazova Yu.A., Shiyan O.V. Theorem on the existence and stability of the solution of one parabolic equation. Dinamicheskie Sistemy, 2018, vol. 8(36), no. 3, pp. 275-280 (in Russian).

6. Khazova Yu.A. Traveling waves solution in parabolic problem with a rotation. Izvestiya VUZ. Applied Nonlinear Dynamics, 2017, vol. 25, no. 6, pp. 57-69 (in Russian).

7. Belan E.P., Khazova Yu.A. Dynamics of stationary structures in a parabolic problem with reflection spatial variable in the case of a circle. Dinamicheskie Sistemy (Dynamical Systems), 2014, vol. 4(32), no. 1-2, pp. 43-57 (in Russian).

8. Khazova Yu.A. Dynamics of stationary structures in a parabolic problem with reflection spatial variable in the case of a segment. Dinamicheskie Sistemy (Dynamical Systems), 2014, vol. 4(32), no. 3-4, pp. 245-257 (in Russian).

9. Khazova Yu.A. Stationary structures in a parabolic problem with reflection spatial variable. Tavricheskiy vestnik matematiki i informatiki, 2015, no. 3(28), p. 82-95 (in Russian).

10. Khazova Yu.A. Integral representation of approximate solutions for a parabolic equation. Aktualniye napravleniya nauchnih issledovaniy XXI veka: Teoriya i praktika, 2018, vol. 6, pp. 384-386 (in Russian).

11. Lions J.L., Magenes E. Problemes aux Limites Non Homogenes et Applications. Vol. 1, 2. Paris: Dunod, 1968.

12. Temam R. Infinite-Dimensional Dynamical Systems in Mechanics and Physics. Springer-Verlag, New York, 1997.

13. Leray J., Schauder J. Topology and functional equations. Uspehi Matem. Nauk, 1946, vol. 1, no. 3-4, pp. 71-95.

14. Babin A.V., Vishik M.I. Attractors of Evolution Equations. North-Holland, Amsterdam, 1992.

Хазова Юлия Александровна - родилась в Севастополе (1989), окончила Таврический национальный университет имени В.И. Вернадского (2011) по специальности математика. Работает старшим преподавателем кафедры дифференциальных уравнений и геометрии факультета математики и информатики Таврической академии (структурное подразделение) Крымского федерального университета имени В.И. Вернадского. Опубликовала 14 научных статей по направлениям нелинейная оптика, параболические уравнения с преобразованиями пространственной переменной.

Россия, Республика Крым, 295007 Симферополь, просп. Академика Вернадского, д. 4, ФГАОУ ВО «Крымский федеральный университет им. В.И. Вернадского» E-mail: hazova.yuliya@hotmail.com

Лукьяненко Владимир Андреевич - родился в с. Юрковка Винницкой обл. (1949), окончил Одесский государственный университет имени И.И. Мечникова (1972) по специальности математика. Работает доцентом кафедры дифференциальных уравнений и геометрии факультета математики и информатики Таврической академии (структурное подразделение), Крымский федеральный университет имени В.И. Вернадского. Опубликовал свыше 130 научных статей по направлениям уравнения типа свертки, краевые задачи теории аналитических функций и уравнения в частных производных, линейные и нелинейные уравнения 1-го рода, матмоделирование, интеллектуализация обработки данных.

Россия, Республика Крым, 295007 Симферополь, просп. Академика Вернадского, д. 4 ФГАОУ ВО «Крымский федеральный университет им. В.И. Вернадского» E-mail: art-inf@yandex.ru