МАТЕМАТИКА
MATHEMATICS
УДК 517.927.2
doi: 10.21685/2072-3040-2023-3-1
О решении нелинейного интегрального уравнения Липпмана - Швингера методом сжимающих отображений
Ю. Г. Смирнов1, Д. А. Лабуткина2
1,2Пензенский государственный университет, Пенза, Россия 1mmm@pnzgu.ru, 2labutkinadara@gmail.com
Аннотация. Актуальность и цели. Рассматривается скалярная трехмерная краевая задача дифракции волны на неоднородном рассеивателе для уравнения Гельмгольца с нелинейной зависимостью волнового числа рассеивателя от поля. Материалы и методы. Исходная задача сводится к объемному нелинейному интегральному уравнению Липпмана - Швингера по рассеивателю. Для исследования интегрального уравнения применен метод сжимающих отображений. Результаты. Доказаны существование и единственность решения в пространстве непрерывных функций при некоторых ограничениях на параметры задачи. Выводы. Доказана сходимость итерационного процесса в методе сжимающих отображений и представлена оценка скорости сходимости.
Ключевые слова: уравнение Гельмгольца, интегральные уравнения, метод сжимающих отображений, разрешимость краевой задачи, численный метод
Финансирование: исследование выполнено за счет гранта Российского научного фонда № 20-11-20087.
Для цитирования: Смирнов Ю. Г., Лабуткина Д. А. О решении нелинейного интегрального уравнения Липпмана - Швингера методом сжимающих отображений // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. 2023. № 3. С. 3-10. doi: 10.21685/2072-3040-2023-3-1
On the solution of the nonlinear Lippmann - Schwinger integral equation by the method of contracting maps
Yu.G. Smirnov1, D.A. Labutkina2
1,2Penza State University, Penza, Russia 1mmm@pnzgu.ru, 2labutkinadara@gmail.com
Abstract. Background. A scalar three-dimensional boundary value problem of wave diffraction on an inhomogeneous scatterer for the Helmholtz equation with a nonlinear dependence of the scattering wavenumber on the field is considered. Materials and methods. The boundary value problem is reduced to the volume nonlinear Lippmann - Schwinger integral equation on the scatterer. The method of contracting maps is used to study the integral equation. Results. The existence and uniqueness of the solution in the space of continuous functions under certain conditions on the parameters of the problem are proved. Con-
© Смирнов Ю. Г., Лабуткина Д. А., 2023. Контент доступен по лицензии Creative Commons Attribution 4.0 License / This work is licensed under a Creative Commons Attribution 4.0 License.
elusions. The convergence of the iterative process in the method of contracting maps is proved and an estimate of the convergence rate is presented.
Keywords: Helmholtz equation, integral equations, method of contracting maps, solvability of boundary value problem, numerical method
Financing: the research was financed by the grant No. 20-11-20087 of the Russian Science Foundation.
For citation: Smirnov Yu.G., Labutkina D.A. On the solution of the nonlinear Lippmann -Schwinger integral equation by the method of contracting maps. Izvestiya vysshikh uehebnykh zavedeniy. Povolzhskiy region. Fiziko-matematieheskie nauki = University proceedings. Volga region. Physical and mathematical sciences. 2023;(3):3-10. (In Russ.). doi: 10.21685/2072-3040-2023-3-1
Введение
Краевые задачи сопряжения для уравнения Гельмгольца встречаются во многих разделах математической физики: акустике, механике, электродинамике. В частности, это задачи дифракции акустических или электромагнитных волн на рассеивателях. Основные типы этих задач исследованы достаточно подробно [1-4]. Однако в последнее время появился интерес к краевым задачам, в которых волновая функция нелинейно зависит от искомого решения. В качестве примера можно рассмотреть случай нелинейного мета-материала [5].
Одним из наиболее эффективных методов решения краевых задач сопряжения является метод интегральных уравнений. Краевая задача сопряжения сводится к объемному интегральному уравнению Липпмана - Швингера по области рассеивателя. Такой подход позволяет не только исследовать свойства и разрешимость задачи, но и получать ее численное решение. Вычислительные алгоритмы, построенные для решения интегральных уравнений, легко распараллеливаются, что позволяет использовать суперкомпьютеры для их решения.
В настоящей статье мы рассмотрим вопросы существования и единственности решения для краевой задачи дифракции и интегрального уравнения, а также докажем сходимость итерационной процедуры для метода сжимающих отображений в пространстве непрерывных функций и получим оценку скорости сходимости.
1. Постановка задачи
Рассмотрим задачу дифракции на рассеивателе (в области) Q с гладкой
2 3
границей dQ класса C , расположенном в однородном пространстве R , характеризующемся постоянным волновым числом k0 .
Пусть U0 (x) - падающая волна - известная функция, удовлетворяющая
в R3 однородному уравнению Гельмгольца
Au0 + k) U0 = 0,
описывающая источник поля. Например, это может быть плоская волна U0(x) = e'k°Х3,x = (,X2,X3)e R3.
В области Q среда неоднородна и характеризуется волновой функцией к(х;и) . На дQ будем определять только предельные значения волновой
функции и волнового числа с разных сторон поверхности. Требуется определить решение задачи сопряжения (полное поле), где искомая (комплексно-значная) функция и( х) должна удовлетворять:
- условиям гладкости:
ие С2^пС1 ^)пС1 (3\Q); (1)
- уравнению Гельмгольца:
Лм + к 2( x; u)u = 0, (2)
где
к 2( x; u) = <
к02,xе R3\Q,
02 ,2 (3)
к (x) + alu\ ,x е Q;
- условию излучения Зоммерфельда для рассеянного поля и3 = и — и^ :
— ¡к0и8 = о|1|, г = |х| (4)
дг I г )
- условиям сопряжения на границе dQ :
du
[u ]9Q = 0
dn
= 0, (5)
dQ
где [ • - разность следов функции с разных сторон дQ .
2
Здесь п обозначает внешнюю нормаль к области Q ; функция к (х; и) описывает вид нелинейной зависимости волновой функции в области Q ; не-
2 2
прерывная функция к (х) > к0 определяет характер неоднородности рассеи-вателя; а> 0 - коэффициент нелинейности - положительная постоянная, характеризующая уровень нелинейности поля. Зависимость вида (3) для волновой функции является аналогом закона Керра [6].
2. Сведение краевой задачи к интегральному уравнению
Перепишем уравнения (2), (3) в следующем виде:
Ди(х) + к02и(х) = 0, х е Я3 \ Q,
2 /99 \
Ди(х) + к0и(х) = (( — к (х;и) )и(х), хе Q . (6)
Обозначим
G = G (x, y) =
eik0\x-у\
4n|x - y| '
Применяя вторую формулу Грина к функциям u(y) и G(x,y) в R , с учетом условий излучения (4), аналогично [7, с. 55], получаем интегральное уравнение в Q :
u(x) = JG(x,y)(k2(y;u)-k0)u(y)dy + u0(x), xe Q , (7)
Q
где x = (xj,x2,x3) и y = (yj,y2,Уз)e R . Полученное уравнение естественно называть нелинейным уравнением Липпмана - Швингера. Заметим, что если коэффициент нелинейности а = 0, то (7) совпадает с известным (линейным) уравнением Липпмана - Швингера [7, 8].
Учитывая свойства объемного потенциала [9, с. 395], будем искать решения уравнения (7) в пространстве непрерывных функций C(Q). Запишем уравнение (7) в операторном виде:
u = A(u), (8)
где оператор A(u): C(Q) ^ C(Q) определяется формулой
A(u)(x) := J G(x,y) (k2(y;u) - k0 )u(y)dy + uo(x).
Q
Введем также линейный оператор Aou : C(Q) ^ C(Q) по формуле
Aou (x) = J G(x, y)u(y)dy .
Q
Тот факт, что эти операторы непрерывны в пространстве C (Q), является следствием вида волновой функции (3) и свойств объемного потенциала [9, с. 395].
3. Существование и единственность решения интегрального уравнения
Пусть R > 0 - фиксированное число такое, что ||uo|C < R . Будем рассматривать уравнение (8) в шаре радиуса R: u e Cr (Q):=
:={u e C(Q): ||u||c ^ R} . Введем обозначение K := max k2(x) - k0 , K > 0.
xeQ
Пусть также М :=||Сд)^с(д) • Тогда справедлива Лемма 1. Если выполняется условие
М(КЯ + аЯ3) и01С < Я, (9)
то (нелинейный) оператор Л действует в шаре Ся (б): Л( ): Ся (б) ^ Ся (б) • Доказательство. Достаточно проверить следующую цепочку неравенств при условии, что ||и||с < Я :
||A(u)||C < M (к2(x;u) - £o)u
+
C
C
<
< M ( K\u\c + a
\u\|C 1 +
c
< M ((
KR + aR3 )+ uj < R.
Теперь можно доказать следующее утверждение. Теорема 1. Если выполнено условие
М (К + 3aR2) < 1,
(10)
то оператор А() : Cr (Q) ^ Cr (Q) сжимающий.
Доказательство. Пусть щ е Cr (Q), ^ е Cr (Q). Рассмотрим разность
А(щ) - A(«2) = JG(x,y)(к2 (y;щ) - ко )щ (y)dy -
Q
-J G(x, y)(k2(y; u2) - к,2 u (y)dy =
Q
= J G(x,y) ((к2 (y) - к,2Xu (y) - u2 (y)) + a(u (y) Ц (y)|2 - u2 (y) |u2 (y)|2 ))dy
Q
и оценим ее по норме:
||A(u ) - A(u2 )||C < M (к2(y) - ko ) ) - u2) + a(ux \ux |2 - u2 |u21
C
<M^K||u1 -«2||C +a (ui \щ\2 -U2 |mi|2 ) + (U2 |mi|2 -U2 |м21 <M(к\щ -U21|C + aR2 Щ -U21|C + a|Щ2(|щ\ -|u2щ\ + |u2|C ) < < M ( К U - U21|C + aR2 ||uj - U21|C + 2aR21\u1 - u2 ||C ) <
< M (к + 3aR2 ) - u2||C ,
откуда следует, что при выполнении условия (10) оператор будет сжимающим. При оценке использовалось простое, справедливое для любых комплексных чисел zj и ^2 неравенство - <- z^ , которое легко получить, записав эти числа в тригонометрической форме.
Из принципа сжимающих отображений [10, с. 351], леммы 1 и теоремы 1 следует, что решение уравнения (8) существует и единственно при выполнении условий (9), (10). Таким образом, верна
Теорема 2. Пусть выполнены условия (9), (10). Тогда решение уравнения (8) существует и единственно в шаре Cr (Q).
Оценки (9), (10) выполняются, в частности, при малых значениях коэффициента нелинейности a и, одновременно, при малых величинах K .
4. Сходимость итерационного метода
В качестве следствия из теоремы 1 и принципа сжимающих отображений [10, с. 351] получаем
Утверждение 1. При выполнении условий (9), (10) итерационный метод
и(и+1) = А(и(п)), п = 0,1,2,...,
(11)
сходится
ком, что
и — и
(п)
С
^ 0, п ^ ^ при любом начальном приближении и(0) та-
С
< Я, со скоростью геометрической прогрессии с показателем
q < 1, и верна оценка скорости сходимости
и — и
( п)
С 1 — <
А(и(0)) — и(0)
С
где q := М(К + 3аЯ ), а и е Ся (б) - точное решение уравнения (8).
В качестве хорошего начального приближения можно взять решение линейного уравнения при а = 0, которое подробно исследовано, численные методы для его решения разработаны [7].
5. Разрешимость краевой задачи
Пусть и е Ся (б) - (точное) решение интегрального уравнения (7) (или (8)). Рассмотрим интегральное представление
и(х) = |С(х,у)(к2(у;и) — к%)иШу + и0(х), хе Я3 \б . (12)
б
Докажем, что функция и е Ся (б), определенная в б и продолженная 3 —
в Я \ б по формуле (12), является решением краевой задачи (1)-(5). В области Я3 \ б функция и е С(Я3 \ б) удовлетворяет уравнению (2) за счет выбора ядра G(х,у), поскольку в (12) хй б и можно дифференцировать (12) под знаком интеграла. Условие излучения (4) также выполняется автоматически за счет выбора ядра G(х, у).
В силу свойств объемного потенциала [9, с. 395] правая часть уравнения (7) есть непрерывно дифференцируемая функция в Я3 . Тогда левая часть
уравнения (7) также непрерывно дифференцируема в б и, по формуле (12),
3
в Я \ б и условия сопряжения (4) выполняются. Поскольку теперь
и е С1 (б), то, снова применяя утверждение о гладкости объемного потенциа-
2
ла [9, с. 395], получаем, что и е С (б), и условия гладкости (1) также выполняются.
Остается проверить справедливость уравнения (2) в области б. Имеем
п
Au (x) + £q u( x) = J (A x + k0 )G( x, y) [k2 (y; u) - k0 )u( y)dy + (A x + k0 )u0 (x) =
Q
- J 5(x - y) [ k2(y; u) - £2) u (y )dy = - [ k2(x; u) - £q ) u (x),
Q
2
откуда следует, что Au (x) + k (x; u(x))u (x) = 0. Таким образом, доказана
Теорема 3. Решение u( x) интегрального уравнения (7) (или (8)) вместе с интегральным представлением (12) дает единственное решение краевой задачи (l)-(5).
Заключение
В работе рассмотрена нелинейная краевая задача дифракции волны на неоднородном рассеивателе для уравнения Гельмгольца с нелинейной зависимостью волнового числа рассеивателя от интенсивности поля. С помощью объемного потенциала исходная задача сведена к исследованию нелинейного интегрального уравнения Липпмана - Швингера по области рассеивателя. Доказаны существование и единственность решения как интегрального уравнения, так и краевой задачи дифракции при некоторых ограничениях на параметры задачи. Для доказательства применен принцип сжимающих отображений.
Для численного решения интегрального уравнения предложен итерационный метод сжимающих отображений. Доказана сходимость метода при произвольном начальном приближении при некоторых ограничениях на параметры задачи. Получена оценка скорости сходимости.
Список литературы
1. Ладыженская О. А. Краевые задачи математической физики. М. : Наука, 1973. 540 с.
2. Санчес-Паленсия Э. Неоднородные среды и теория колебаний. М. : Мир, 1984. 472 с.
3. Nedelec J-C. Acoustic and Electromagnetic Equations. Integral Representations for Harmonic Problems. Springer Science+Business Media, 2001. Vol. 144. 316 p.
4. Колтон Д., Кресс Р. Методы интегральных уравнений в теории рассеяния. М. : Мир, 1987. 311 с.
5. Нелинейности в периодических структурах и метаматериалах / под ред. Ю. С. Кив-шаря, Н. Н. Розанова. М. : Физматлит, 2015. 384 с.
6. Ахмедиев Н. Н., Анкевич А. Солитоны, нелинейные импульсы и пучки. М. : Физматлит, 2003. 304 с.
7. Смирнов Ю. Г., Цупак А. А. Математическая теория дифракции акустических и электромагнитных волн на системе экранов и неоднородных тел. М. : Русайнс, 2016. 226 с.
8. Colton D., Kress R., Inverse Acoustic and Electromagnetic Scattering Theory. Springer Science+Business Media, 2013. 405 p.
9. Владимиров В. С. Уравнения математической физики. М. : Наука, 1981. 512 с.
10. Треногин В. А. Функциональный анализ. М. : Наука, 1993. 440 с.
References
1. Ladyzhenskaya O.A. Kraevye zadachi matematicheskoy fiziki = Boundary value problems of mathematical physics. Moscow: Nauka, 1973:540. (In Russ.)
2. Sanches-Palensiya E. Neodnorodnye sredy i teoriya kolebaniy = Inhomogeneous media and the theory of vibrations. Moscow: Mir, 1984:472. (In Russ.)
3. Nedelec J-C. Acoustic and Electromagnetic Equations. Integral Representations for Harmonic Problems. Springer Science+Business Media, 2001;144:316.
4. Kolton D., Kress R. Metody integral'nykh uravneniy v teorii rasseyaniya = Methods of integral equations in scattering theory. Moscow: Mir, 1987:311. (In Russ.)
5. Kivshar Yu.S., Rozanova N.N. (eds.). Nelineynosti vperiodicheskikh strukturakh i met-amaterialakh = Nonlinearities in periodic structures and metamaterials. Moscow: Fiz-matlit, 2015:384. (In Russ.)
6. Akhmediev N.N., Ankevich A. Solitony, nelineynye impul'sy i puchki = Solitons, nonlinear pulses and beams. Moscow: Fizmatlit, 2003:304. (In Russ.)
7. Smirnov Yu.G., Tsupak A.A. Matematicheskaya teoriya difraktsii akusticheskikh i el-ektromagnitnykh voln na sisteme ekranov i neodnorodnykh tel = Mathematical theory of diffraction of acoustic and electromagnetic waves on a system of screens and inhomogeneous bodies. Moscow: Rusayns, 2016:226. (In Russ.)
8. Colton D., Kress R., Inverse Acoustic and Electromagnetic Scattering Theory. Springer Science+Business Media, 2013:405.
9. Vladimirov V.S. Uravneniya matematicheskoy fiziki = Equations of mathematical physics. Moscow: Nauka, 1981:512. (In Russ.)
10. Trenogin V.A. Funktsional'nyy analiz = Functional analysis. Moscow: Nauka, 1993:440. (In Russ.)
E-mail: labutkinadara@gmail.com
Авторы заявляют об отсутствии конфликта интересов / The authors declare no conflicts of interests.
Поступила в редакцию / Received 24.04.2023
Поступила после рецензирования и доработки / Revised 10.06.2023 Принята к публикации / Accepted 20.08.2023
Информация об авторах / Information about the authors
Юрий Геннадьевич Смирнов
доктор физико-математических наук, профессор, заведующий кафедрой математики и суперкомпьютерного моделирования, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
Yuriy G. Smirnov
Doctor of physical and mathematical sciences, professor, head of the sub-department of mathematics and supercomputer modeling, Penza State University (40 Krasnaya street, Penza, Russia)
E-mail: mmm@pnzgu.ru
Дарья Андреевна Лабуткина
аспирант, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
Dar'ya A. Labutkina
Postgraduate student, Penza State University (40 Krasnaya street, Penza, Russia)