Нелинейные одно- и двухпараметрические задачи сопряжения на собственные значения для системы уравнений Максвелла в слое тема автореферата и диссертации по математике, 01.01.02 ВАК РФ

Валовик, Дмитрий Викторович АВТОР
доктора физико-математических наук УЧЕНАЯ СТЕПЕНЬ
Пенза МЕСТО ЗАЩИТЫ
2014 ГОД ЗАЩИТЫ
   
01.01.02 КОД ВАК РФ
Автореферат по математике на тему «Нелинейные одно- и двухпараметрические задачи сопряжения на собственные значения для системы уравнений Максвелла в слое»
 
Автореферат диссертации на тему "Нелинейные одно- и двухпараметрические задачи сопряжения на собственные значения для системы уравнений Максвелла в слое"

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Пензенский государственный университет»

Нелинейные одно- и двухпараметрические

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

Специальность 01.01.02 - Дифференциальные уравнения, динамические системы и оптимальное управление

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

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

Валовик Дмитрий Викторович

СЕН 2014

Пенза - 2014 005552704

005552704

Работа выполнена на кафедре математики и супсркомпьютерного моделирования Федерального государственного бюджетного образовательного учреждения высшего профессионального образования «Пензенский государственный университет».

Научный консультант:

Официальные оппоненты:

Ведущая организация:

доктор физико-математических наук, профессор, зав. кафедрой «Математика и суперкомньютсрное моделирование» ФГБОУ ВПО «Пензенский государственный университет» Смирнов Юрий Геннадьевич

Карчевский Евгений Михайлович, доктор физико-математических наук, ФГАОУ ВПО «Казанский (Приволжский) <]>сдсральный университет» (г. Казань), профессор кафедры «Прикладная математика»;

Тихонравов Александр Владимирович, доктор физико-математических наук, профессор, ФГБОУ ВПО «Московский государственный университет им. М. В. Ломоносова» (г. Москва), директор Научно-исследовательского вычислительного центра (НПВЦ);

Шестопалов Юрий Викторович, доктор физико-математических наук, профессор, ФГБОУ ВПО «Московский государственный технический университет радиотехники, электроники и автоматики» (МГТУ МПРЭА) (г. Москва), ведущий научный сотрудник

Институт вычислительной математики (ИВМ) РАН (г. Москва)

Защита состоится 15 октября 2014 г. в 15 часов 30 минут на заседании диссертационного совета Д.501.001.43 ФГБОУ ВПО «Московский государственный университета им. М. В. Ломоносова» по адресу: г. Москва, ГСП-1, Ленинские горы, МГУ, д. 1, стр. 52, 2-й учебный корпус, факультет вычислительной математики и кибернетики.

С диссертацией и авторефератом можно ознакомиться в Научной библиотеке ФГБОУ ВПО «Московский государственный университета им. N1. В. Ломоносова», а также на официальном сайте факультета вычислительной математики и кибернетики МГУ им. М. В. Ломоносова: http://es.msu.ru в разделе «Диссертации».

Автореферат разослан сентября 2014 г.

Ученый секретарь

диссертационного совета ^ у

доктор физ.-мат. наук, профессор Захаров Евгений Владимиров!!

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

Актуальность темы

Задачи сопряжения на собственные значения для системы уравнений Максвелла возникают при изучении распространения электромагнитных волн в неоднородных волноведущих структурах. К таким задачам относится распространение поляризованных ТЕ- и/или TM-волн в плоских диэлектрических слоях и диэлектрических цилиндрических волноводах1 (интерес представляют в том числе и многослойные структуры2).

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

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

Актуальность исследования задач о распространении ТЕ- н/илн TM-волн в нелинейных волноведущих структурах обусловлена двумя обстоятельствами. Во-нервых, разработка новых методов исследования таких нелинейных задач на собственные значения актуальна с математической точки зрения, поскольку отсутствуют общие методы исследования указанного класса задач. В данной диссертации предложен общий метод исследования указанного класса задач. Эти нелинейные задачи сводятся к отысканию собственных значений (значений постоянной распространения), при которых волна может распространяться. Отметим, что рассматриваемые задачи являются нелинейными как но искомым функциям, так и по спектральному параметру (или парс спектральных параметров). Во-вторых, задачи с «простой» геометрией (плоские слои, круглые цилиндрические волноводы) имеют широкие практические приложения1. Центральной проблемой здесь является определение условий существования постоянных распространения. Знание постоянных распространения необходимо при конструировании волноведущих структур.

^хмедиев Н. II., Анкевич А. Солитопы, нелинейные импульсы и пучки. - М.: Физматлит, 2003 ; Шеи И. Р. Принципы нелинейной оптики. - М.: Паука, 19S9 ; Ponatll II.-Б., Stegeman G. I. (editors) Modem problems in condensed matter sciences. Vol. 29. Nonlinear surface electromagnetic phenomena. - North-Holland: Elsevier Science Publishers, 1991.

2Joannopoulos J. D. ct al. Photonic Crystals. Molding the flow of light. - Princctou and Oxford: Princeton University Press, 2008 ; Lourtioz J.-M. ct al. Photonic Crystals. Towards nanoHcalc photonic devices. - Springer-Verlag Berlin Heidelberg, 2005.

3Адамс M. Введение в теорию оптических волноводов. - М.: Мир, 1984 ; Вайнштейн JI. А. Электромагнитные волны. - М.: Радио и связь, 1988 ; Стрэттон Дж. А. Теория электроыагпегиз-ма. - Москва-Ленинград: ГИТТЛ, 1948.

В последние десятилетия предпринимались попытки обобщения известных результатов линейной теории на случай нелинейных задач. Первые фундаментальные результаты были получены В. М. Елеонским, JI. Г. Оганссьянцем и В. П. Силиным4. Первые математически строгие решения задач сопряжения на собственные значения для ТЕ-волн были нолучены в работах В. С. Серова, Ю. Г. Смирнова, Ю. В. Шестоиало-ва, H.-W. Schürmann'1. Также важные результаты нолучены в работах A. D. Boardman, D. N. Christodoulides, К. М. Leung6. Заметим однако, что общего метода исследования предложено не было, так же как не было получено общих результатов о существовании и локализации собственных значений.

Как известно, в линейной среде ТЕ- и TM-волны распространяются не взаимодействуя. Явление нелинейности приводит к новому, принципиально важному результату: существует режим распространения ТЕ- и TM-волн, в котором ТЕ- и TM-волны, распространяясь каждая со своей постоянной распространения и на своей частоте, взаимодействуют, но сохраняют структуру поверхностных волн, образуя связанную волну. Изучение связанных волн, с одной стороны, интересно с физической точки зрения, поскольку они описывают новые режимы распространения волн в волноведущих структурах, которые, в частности, могут оказаться полезными на практике. С другой стороны, возникает новый класс .математических задач на собственные значения. Постоянные распространения в такой задаче существуют дискретными парами, что соответствует парным собственным значениям, пли двухпараметричсской задаче на собственные значения7. Математические методы исследования таких задач пока также не разработаны.

Взаимодействие между ТЕ- и ТМ-волнамп в нелинейной волноведу-щей структуре с керровской нелинейностью рассматривалось как у нас8,

4EIeonskii Р. N., Oganes'yants L. G., Silin V. P. // Soviet Physics JETP. - 1972. - V. 35, ■V 1. - P. 44-47.

'Смирнов Ю. Г., Куприянова С. H. // Жур. выч. мат. и матсм. физ. - 2004. - Т. 44, 10. -С. 1850-1860 ; Schürmann H.-W., Serov V. S., Shestopalov Yu. V. // Phys. Rev. E. - 1998. -V. 58, № 1. - P. 1040-1050 ; Schürmann H.-W., Serov V. S., Shestopalov Yu. V. // .1. Phys. A: Math. Gen. - 2002. - V. 35, K' 50. - P. 10789-10801 ; Schürmann H.-W., Smirnov Yu. G., Shestopalov Yu. V. // Phys. Rev. E. - 2005. - V. 71, № 1. - P. 016614-1-10 ; Serov V. S., Shestopalov Yu. V., Schürmann H.-W. // Dokl. Math. - 1999. - V. 60. - P. 742-744 ; Smirnov Yu. G., Valovik D. V. // ISRN Math. Phys. - 2013. - V. 2013. - P. 1-7.

6См., например, Leung К. M. // Phys. Rev. Ii. - 1985. - V. 32, .Y" 8. - P. 5093-5101 ; Joseph R. I., Christodoulides D. N. // Optics Letters. - 1987. - V. 12, № 10. - P. 826-828.

7 Atkinson F. V., Mingarelli A. B. Multiparameter eigenvalue problems. Sturm - LiouvilJc theory. -NW: CRC Press, 2011.

8Елеонский В. M., Оганесьянц JI. Г., Силин В. П. // Успехи физ. наук. - 1972. - Т. 107, ■V 3. - С. 516-518 ; Eleonskii V. М., Oganes'yants L. G., Silin V. P. // Soviet Phys. JETP. -1973. - V. 36, № 2. P. 282 285.

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

Цель работы

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

Методы исследования

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

Научная новизна работы

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

Основные результаты диссертации

Результаты диссертации состоят в следующем:

• предложен и развит новый математический подход - метод интегральных дисперсионных уравнений, позволяющий исследовать нелинейные одно- и двухнарамстричсские задачи сопряжения на собственные значения для системы уравнений Максвелла;

• введены понятия собственного значения и парных (пли связанных) собственных значений для некоторых классов нелинейных задач сопряжения на собственные значения;

9Воаг(1тап А. Б., Т\уаг<1(т'вИ Т. // Л. Ор1. Бос. Лш. 13. - 1988. - V. 5, .V 2. - Р. 523-528 ; Воаг(1тап А. Б., Ти^агастзЫ Т. // РЬув. Г1ст. А. - 1989. - V. 39, № 5. - Р. 2481-2-191.

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

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

• в результате исследования найдены новые типы нелинейных волн (ТЕ-, ТМ-, ТЕ-ТМ-волн) в изученных волноведущих структурах.

Теоретическая и практическая ценность

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

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

Обоснованность и достоверность результатов

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

Апробация работы

Основные результаты диссертации докладывались на следующих конференциях и семинарах:

• семинаре кафедры физики Osnabrück University, руководитель -проф. H.-W. Schürmann (Германия, г. Оснабрюк, 2010, 2011);

• семинаре но электродинамике факультета ВМК, МГУ им. М. В. Ломоносова, руководители проф. Е. В. Захаров и проф. А. С. Ильинский (Россия, г. Москва, 2012);

• семинаре кафедры «Прикладная математика» Казанского (Приволжского) федерального университета, руководитель - проф. Н. Б. Пле-щинскнй (Россия, г. Казань, 2013);

• семинаре кафедры Elcctrical, Electronic, and Communication Engineering университета Chuo, руководитель - проф. К. Kobayashi (Япония, г. Токио, 2013);

• семинаре кафедры Elcctrical Engineering университета Nilion, руководитель - проф. Т. Yamasaki (Япония, г. Токио, 2013);

• семинаре Computational and Applied Mathematics университета Chalmers, руководитель - проф. S. Larsson (Швеция, г. Гстеборг, 2013);

• семинаре «Вычислительная математика н приложения», Институт вычислительной математики РАН, руководитель - чл.-корр. РАН, проф. Е. Е. Тыртышииков (Россия, г. Москва, 2013);

• научно-методологическом семинаре НИВЦ МГУ им. М. В. Ломоносова, руководитель - д.ф.-м.н., проф. А. В. Тихонравов (Россия, г. Москва, 2013);

• семинаре, руководимом акад. РАН, проф. В. А. Ильиным и акад. РАН, нроф. Е. И. Моисеевым (факультет ВМК, МГУ им. М. В. Ломоносова, Россия, г. Москва, 2013);

• международных конференциях «Days on Diffraction» (Россия, 1'. Санкт-Петербург, 2007, 2011, 2013);

• 13-й и 14-й международных конференциях «Mathematical Methods in Electromagnetic Theory» (Украина, г. Киев, 2010; г. Харьков, 2012);

• международных конференциях «Progress in Electromagnetic Research Symposium» (Китай, г. Сучжоу, 2011; Малайзия, г. Куала Лум-пур, 2012; Россия, г. Москва, 2012);

• международном семинаре «Workshop on Largc-Scale Modeling» (Швеция, г. Сунне, 2012);

• международной конференции «International Symposium on Electromagnetic Theory» (Япония, г. Хнрошнма, 2013).

Работа была поддержана грантами РФФИ (.V 06-07-89063а, 20082009; № 12-07-97010-р_А, 2012-2013; № 11-07-00330-А, 2011-2012), ФЦП

7

(Л* 2.1.1/1647, 2009-2011; № 14.В37.21.1950, 2012-2013), программы ^эЬу (2012-2013, Швеция), Грантами Президента РФ (МК-2074.2011.1, 20112012; МК-90.2014.1, 2014-2015).

Публикации

Основные результаты диссертации отражены в 29 научных работах (2 монографии; 27 статей, из них 7 публикаций без соавторов). Все указанные статьи опубликованы в реферируемых журналах (статьи [3-25] опубликованы в журналах, входящих в перечень ВАК РФ, в которых рекомендуется публиковать основные научные результаты диссертаций на соискание ученых степеней доктора и кандидата наук), обе монографии также прошли рецензирование. В совместных работах профессору Ю. Г. Смирнову принадлежит первоначальная постановка задач, аспирантам Е. В. Зарембо и Е. Ю. Смолькину - программная реализация некоторых численных методов, Д. В. Валовику - получение конкретных результатов и их доказательства.

Структура и объем работы

Диссертация изложена на 155 страницах, включающих 11 рисунков, и состоит из введения, трех глав, списка литературы, содержащего 127 наименований, и четырех приложений.

Краткое содержание работы

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

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

Рассматриваются монохроматические ТЕ-волны Ее~гш(, Не~,ш<, распространяющиеся вдоль границы однородного, изотропного, немагнитного диэлектрического слоя Е := {(ж, у, г) е М3 : 0 ^ х ^ /¿}, здесь

Е = (0, Еу, 0)т, Н = (Нх, О, Н,)т— (1.1)

комплексные амплитуды; ш - круговая частота.

Слои £ расположен в декартовой системе координат Охуг между двумя полупространствами х < 0 и х > Н. Полупространства заполнены изотропной немагнитной средой без источников и имеют постоянные

диэлектрические проницаемости £i J: £q и £3 ^ £0 соответственно, где £о > 0 - диэлектрическая проницаемость вакуума. Считаем, что всюду fi = Цо, где /хо >0 — магнитная проницаемость вакуума.

Диэлектрическая проницаемость £ волновода И описывается следующим выражением: е = е2 + Д|Е|2), где £г > max(ei, £3), /бС1 [0, +ос), /(s2) ^ 0 и /(0) = 0.

Комплексные амплитуды (1.1) удовлетворяют стационарным уравнениям Максвелла

rot H = —iweЕ,

(1.2)

rotE = гш/iH, ;

условию непрерывности касательных составляющих компонент поля на границах раздела сред х = 0, х = h и условию излучения на бесконечности: электромагнитное поле затухает как 0(|ж|-1) при |а:| —> оо.

Доказано, что компоненты ТЕ-волн, распространяющихся вдоль границы слоя имеют представление

Еу = E„(a:)ei72, Нх = Y{x(x)eilz, Я, = Hz{x)élz, (1-3)

где 7 - неизвестный (действительный) спектральный параметр (постоянная распространения электромагнитной волны).

Обозначим А-ц := ur/4)£о- Подставляя компоненты (1.3) в систему (1.2), выполняя нормировку в соответствии с формулами х = кцх, S = 7 = ëj = = 1.2,3), обозначая Y(x) := Еу(х) и

опуская значок тильды, получаем

Y"(x) = (j2-s)Y(x), (1.4)

где

£i, х < 0;

£ = £2 + f(Y2), О^х^ h; (1.5)

гз, х > h.

Функция Y дифференцируема так, что

Y{x) б С\-оо, +00) п С2(—оо, 0) п С2(0, h) П C2(h, +00). (1.6)

Условия сопряжения для функции Y следуют из условий непрерывности касательных компонент электромагнитного поля и имеют вид

MUo = 0, MU = 0, [y']U = 0, [Y']\x=h = о, (1.7) .'Ac[/lU0 = xUm_o/(x-)- Jm+o/(x).

Определение 1.1. Число 7 = 7 такое, что при фиксированном значении У(0) ^ 0 (без потери общности можно считать Y(0) > 0) существует не равная тождественно нулю функция Y(x; 7), которая удовлетворяет уравнению (1.4), условиям (1.6), (1.7) и затухает как 0(|а.'|-1) при |ж| —» оо, будем называть собственным значением, а функцию У (ж; 7), соответствующую этому собственному значению, собственной функцией.

Задача Ре'- доказать существование собственных значений 7, удовлетворяющих определению 1.1.

Разыскиваются такие положительные значения 7, что справедливо max(ei,£3) < 72 < е2- Это условие соответствует классической задаче распространения ТЕ-волн в слое с постоянной диэлектрической проницаемостью1, поэтому мы придерживаемся его при изучении задачи Ре-Пусть к\ := Т -£i> 0, fcf := г2 - 72 > 0, к\ := -у2 - £3 > 0. При i < О II х > h уравнение (1.4) является линейным. Учитывая условие на бесконечности, получаем, что решения уравнения (1.4) в указанных областях имеют вид

Y(x) = {Х<°; (18)

Постоянная А в (1.8) отвечает значению F(0) (см. определение 1.1) и предполагается фиксированной (известной).

Постоянная В определяется условиями сопряжения (1.7). Внутри слоя £ уравнение (1.4) принимает вид

У" = - (k¡ + f(Y2)) Y. (1.9)

Первый интеграл рассматриваемого уравнения имеет вид

Y'2 + k¡Y2 + ifi{Y2) = С, (1.10)

где íp(Y2) = /0Г f{u)du\ С - постоянная.

Постоянная С вычисляется из первого интеграла (1.10) с использованием решения (1.8) и условий сопряжения (1.7) и равна

С= (£2-£l)A2 + v(A2). (1.11)

Неизвестная постоянная В определяется из уравнения

(£2 - £3)В2 + р(В2) = (£2 - £1)А2 + <р(А2), (1.12)

которое имеет два действительных решениях ±В'.

'Вайнштейн J1. А. Электромагнитные волны. - М.: Радио и связь, 1988.

10

Постоянные В и С не зависят ни от точки х = h, ни от спектрального параметра 7; кроме того, С > 0.

Доказано, что дисперсионное уравнение задачи Ре имеет вид

Ф(г,п) :=Tl + nT2 = h, (1.13)

где

к\ +СО

Ti= wdr], т2= wdr¡ и w = —-;

J J Щ + /(т) + rf

-ос

функция г s т(г)) определяется из уравнения (if + fc|)r + (р(т) = С; С определена формулой (1.11); п = 0, 1, 2,...

Важно отметить, что дисперсионное уравнение (1.13) не зависит от решений уравнения (1.12).

Фактически уравнение (1.13) является семейством (но не системой) уравнений для различных п. Необходимо решать относительно 7 каждое из получающихся уравнений. Другими словами, пусть а^ множество решений дисперсионного уравнения (1.13). Тогда ст^ = (Jío0"^ где (Tj содержит все положительные решения (н только их) уравнения Ф(7; j) — h = 0. Кроме того, cr¿ П crj = 0 для всех возможных i / j.

Пусть аЕ - совокупность всех собственных значений задачи Ре- Справедлива следующая теорема об эквивалентности, которая, в частности, утверждает, что cíe = <7е ■

Теорема 1.1 (об эквивалентности). Значение 7 = 7 является собственным значением задачи Ре тогда и только тогда, когда существует целое число п = п ^ 0 такое, что 7 = 7 удовлетворяет уравнению Ф(7;п) - h = 0.

Кроме того, собственная функция Y(x; 7) имеет в точности ri пулей при х G (0, К); если x¡ является i-м нулем функции Y(x;7), то

Теорема 1.2 (о периодичности). Пусть 7 собственное значение задачи Ре- Если собственная функция У(х; у) имеет более двух пулей при х 6 (0,/¿), тогда функция У(х;у) периодическая с периодом 2То.

Пусть Г = (max(el5 £3), £2) « '4f = illf h* = зирФ(7;А;).

72е Г 72еГ

Величина h¡ní всегда существует. Указанная sup существует не всегда,

например, sup не существует, если / = 0. Когда мы пишем /j£u мы

предполагаем, что указанный (конечный) sup существует.

Теорема 1.3. Пусть для некоторого р выполняется /ifnf < hgup, f 6 С1 [0,+оо) и h таково, что h?nf < h < h%up, тогда задача Ре имеет по крайней мере одно решение (собственное значение).

Теорема 1.4. Пусть Ф(7;/с) неограничепа при 72 6 Г, / Е С1 [0,+оо) и h таково, что для некоторого р выполняется hfn{ < h, тогда задача Ре имеет по крайней мере одно решение (собственное значение).

В теоремах 1.3, 1.4 утверждается существование собственных значений; в теоремах 1.5, 1.6 - их изолированность.

Теорема 1.5. Пусть для некоторого р выполняется h\ní < hgup, функция f(u) является аналитической функцией в С (как функция комплексного переменного и = Y2) и h таково, что h\ní < h < h^up, тогда задача Ре имеет по крайней мере одно решение (собственное значение).

Кроме того, если условие /ífnf < h^np справедливо для всехр, то множество собственных значений Ое задачи Ре является дискретным на D = {7 G R : 72 € Г}, т.е. на каждом отрезке I С D содерэюится не более конечного числа (изолированных) собственных значений.

Теорема 1.6. Пусть функция Ф(7;А;) неограничепа при 72 € Г, функция f(u) является аналитической функцией в С (как функция комплексного переменного и = Y2) и h таково, что для некоторого р выполняется h?n{ < h, тогда задача Ре имеет по крайней мере одно решение (собственное значение).

Кроме того, если условие h-mí < h%ap справедливо для всехр, то множество собственных значений ag задачи Ре является дискретным на D = {7 б R : 72 £ Г}, т.е. на каэ/едом отрезке I С D содержится не более конечного числа (изолированных) собственных значений.

Величины h¡ní и h^up можно находить численно.

Условие 72 < £о является точным. Действительно, для функции / = О получаем линейную задачу. В такой линейной задаче необходимо 72 < е2.

Также в первой главе исследованы две задачи для конкретных функций нелинейности: задача Ргл / = «К2 (закон Ксрра); задача РЕ2: / = (нелинейность с насыщением). Эти нелинейности находят ши-

рокое применение в нелинейной оптике.

В теоремах 1.7, 1.8 используются следующие обозначения для собственных значений 7 задач Peí, Ре2'■ 7i значит, что все собственные значения упорядочены по возрастанию; 7(771) значит, что это собственное значение является решением дисперсионного уравнения (1.13) при

п = т; значком P¿n обозначается линейная задача с постоянным е = е2, такая задача имеет не более конечного числа собственных значении, обозначаемых 7

Теорема 1.7 (Задача Peí)- Пусть rriiii(er1, > е0, тах(е1,ез) < г2, а > 0 и У(0) 0. D этом случае для любого h > 0 задача Ре\ имеет бесконечное число положительных собственных значений 7,- (с точкой накопления на бесконечности).

Собственные значения 7,- имеют следующие свойства:

1) сели 7i,72, ■■. все решения задачи Ре\, то

7и 7г. • • ■ е (шах(£Ь е3), +оо) и lim 7? -> со;

J—>оо J

2) если задача Pl¿" имеет р решений 7i < 72 < ... < 7Р, ото существует Qo > 0 такое, что для любого а = а' < а0 справедливо

if 6 (тах(еьез),£2) и Нш 7; = 7;, i = Т~р,

a'—>Q

где 7ь ... ,7j, - первые р решений задачи Peí при а = а'; 2 У если q > р, то lim 7^ = +00;

а'->+0

3) для больших значений 7 и произвольно малого Д > 0 енухшедливо следующее асимптотическое двойное неравенство:

(1 - Д)7.(т) ^ 7(т) < >/2(1 + А)у.(т +1),

гдеу1{т) = £2+(/_1 (j^))2, /_1 есть обращение функции f (t) = Iní; 3') если \/'2aC < 1, где С = (е2 — £i)A2 + 0,5аЛ4, то для больших значений 7 справедливо более простое асимптотическое неравенство 7М ^ 7о(ш), где 7о2(т) = е2 + (f Ih(2qC))2;

4) если собственное значение 7,- —> оо, то max |У(а;;7,)| —¥ оо.

хе(0М)

Теорема 1.8 (Задача Рк)- Пусть miii(ei,£r3) ^ е0, max(ei,6.¡) < е2, а, ß > 0 и К(0) ф 0. Тогда существует hm\n > 0 такое, что для любого h > /¿mill задаем Ре2 имеет конечное число (и не менее одного) полоо/си-тельных собственных значений 7, .

Для всякого собственного значения 7; задачи Ре2 справедливо, что

тах(£Ь e-¿) < 7f <е2 + aß~\

кроме того, для всяких допустимых т и т + 1 справедливо

max 7(771 + 1) < шах 7(т),

где max берется среди всех положительных решений уравнения (1.13) с заданным п = т.

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

Рассматриваются монохроматические ТМ-волны Ее"'"', Не~!Ш<, распространяющиеся вдоль границы однородного, изотропного, немагнитного диэлектрического слоя Е := {(ж, у, z) R3 : 0 ^ х ^ /г}, здесь

Е = (Ех, 0, Ez)'r, Н = (0, Ну, 0)т— (2.1)

комплексные амплитуды; ш - круговая частота.

Слой Е расположен в декартовой системе координат Oxyz между двумя полупространствами х < 0 и х > h. Полупространства заполнены изотропной немагнитной средой без источников и имеют постоянные диэлектрические проницаемости е = Е\ ^ £о и £ = £3 ^ £q соответственно, где га > 0 - диэлектрическая проницаемость вакуума. Считаем, что всюду ц = //о, где /хо > 0 - магнитная проницаемость вакуума.

Диэлектрическая проницаемость е волновода Е описывается следующим диагональным тензором 3x3:

е = diag (ехх, £уу, ezz), (2.2)

элементы £хх, £,, имеют вид

£хх = £/ + £of (а\Ех\2 + b\Ez\2) , е„ = £д + £09 (с\Ех\2 + d\Ez\2) , (2.3)

где £; > max(£i,£3), £q ^ 0 - постоянные (вещественные) составляющие диэлектрических проницаемостей ехх, £zz; a, b,c,d- неотрицательные постоянные, не все равные нулю;

/еС1 [0,+оо), рбС^О.+оо), f(s2)> 0, g(s2)> О

и

/(0) = ,9(0) = 0.

Кроме того, функции / и д являются таковыми, что выполняется соотношение щщ = К0Т0Р0МУ удовлетворяют многие типы нели-нейностей1. В главе 2 доказано, что без этого условия можно обойтись. Элемент еуу не оказывает влияния на распространение ТМ-волн.

"Joseph R. I., Cliristodoulides D. N. // Optics Letters. - 1987. - V. 12, .V 10. - P. 826-828.

Комплексные амплитуды (2.1) удовлетворяют стационарным уравнениям Максвелла

rot Н = —г'шеЕ,

2.4

rot Е = wj/iH,

условию непрерывности касательных составляющих компонент поля на границе раздела сред х = 0, х = h и условию излучения на бесконечности: электромагнитное ноле затухает как 0(\х\ при —> оо.

Доказано, что компоненты ТМ-волн, распространяющихся вдоль границы слоя Е, имеют представление

Ех = Ех{х)е^, Ег = Vz{x)eilz, Ну = Н^К^, (2.5)

где 7 - неизвестный (действительный) спектральный параметр (постоянная распространения электромагнитной волны).

Обозначим кц := иР^фо. Подставляя компоненты (2.5) в систему (2.4), выполняя нормировку в соответствии с формулами х = кдх,

~JZ — h

7 = Ь ёг = * (г = 1,2), = = *

-f-, обозначая

(2.0)

(2.7)

(1х

Z(x) := Ег, Х{х) := гЕг и опуская значок тильды, получаем

-г" + =

где с учетом предыдущего

£1, х < 0;

£; + /{аХ2 + Ьг2),

^£3, X > /г;

(£ь х < 0;

£д + д{сХ2 + аг2), 0 < X ^ к; £з, х > к.

Функции X, Z дифференцируемы так, что Х(х) £ С1 ( оо, 0] П С1 [0, /г] П С1 [/г, +оо),

г{х) € С(—оо, +оо) П С1 ( оо, 0] Г) С^О, Л] п С1 [/г, +оо) П (2.8)

П С2(—оо, 0) П С2(0, к) П С2(Н, +оо).

Условия сопряжения для функций X, Z следуют из условий непрерывности касательных компонент электромагнитного поля и имеют вид

VfX — Z']\x=i) = 0, [Z]\r=Q = 0, bx-z'\\x=h = о, [Z]u = o.

(2.9)

Определение 2.1. Число 7 = 7 такое, что при фиксированном значении Х(0) ф 0 (без потери общности можно считать АГ(0) > 0) существуют не равные тождественно нулю функции Х(х;7), ^(а:;7), которые удовлетворяют системе уравнений (2.6), условиям (2.8), (2.9) и затухают как О(|д:| !) при \х\ оо, будем называть собственным значением, а функции Х(х;у), ^(ж;7), соответствующие этому собственному значению, -собственными функциями.

Задача Рм'- доказать существование собственных значений 7, удовлетворяющих определению 2.1.

Разыскиваются такие положительные значения 7, что справедливо неравенство тах(е1,е3) < 72 < £¡. Это условие соответствует классической задаче распространения ТМ-волн в слое с постоянной диэлектрической проницаемостью (см. Вайнштспн Л. А.), поэтому мы придерживаемся его при изучении задачи Рд/.

Пусть Щ := 72 - £1 > 0, к'1 := -у2 - е3 > 0.

При х < 0 и х > /г система (2.6) является линейной. Учитывая условие на бесконечности, получаем, что решения системы (2.6) в указанных областях имеют вид

уы = \Аекхх' * < 2(х) = х <

(Дс"^«"*), х > /г; ^ {-Ву^кзе-^1^, х>к.

(2.10)

Постоянная А в (2.10) отвечает значению А'(0) (см. определение 2.1) и предполагается фиксированной (известной), а постоянная В определяется условиями сопряжения (2.9).

Внутри слоя Е система (2.6) принимает вид

(2.11)

-г" + 1х' = {еа + д)г,

Система (2.11) может быть переписана в нормальной форме

х, = 72(е,„ + .9) + 2(£/-72 + /)*2/^

7 (2*2/^+ £/ + /) ' (2.12)

Первый интеграл системы (2.11) имеет вид

Х\е; - 72 + I? + 72 (От/ - 72)А'2 + ед2Г2) + 72С = С, (2.13)

где с ^ с(х2, ^2) ЕЕ £ ¡с%2+'122 д{з)аа; С ■

постоянная.

Обозначим через

Л'о := Л"(0 + 0), Хк := Х{к - 0), := Я(0 + 0), Я,, := г {к - 0)

предельные значения функции Л' и Z на границах слоя £ изнутри. Справедливы формулы

е1А = (е/ + ПаХ?}+ЪкЪ-2А°-))Х( ,, = £3В=(е1 + 1(аХ1 + ЪкЬ-1В2))Х,„ гк = -У~1к3В,

которые получены из условии сопряжения (2.9) и решений (2.10).

Поскольку постоянная А предполагается известной, то Л'о определяется из первого уравнения (2.14).

Обозначим /о := / (аХЦ + Ьг^) и (Зо := (? (Х$, Тогда, используя первый интеграл (2.13), подставляя х = 0, найдем

С = Х2(£/ - 72 + /„)2 + 72 ((£/ - 72)Х2 + еаг%) + С0. (2.15)

Неизвестные значения X/, и Еь определяются из системы уравнений

- -уезк^г,, = (е/ + //,)А'„;

(е} - 72 + ЬШ + Г ((*/ - Т)Х1 + ея2%) + С„ = С,

где /„, := ДаХ2 + и С„. := С(Х2, Я2).

При известном значении постоянная В определяется но формуле В =

Доказано, что дпснсрспонное уравнение задачи Рд/ имеет вид

Ф(7; п) := -Т1 + пТ2 = /г, (2.17)

где

4(0) +со

П = I Т2 = ^ -шс1т],

Тх =

V СО

величины ч(0), 77(/г) определяются формулами

Ч(0) = ^1>0 и ч(Л) = -^<0; л-Х п-'з

а функция ы; имеет следующий вид:

ш =_+ я_

72(е/ + /)(£я + .д) + (£/-72 + /)г72'

где

Г- f< \ - ( cr]2 + d(sf + f)2\ f = f(r), (J = g [т , J );

функция r = r(7j) определяется из уравнения

- Г + /)2 + 72(£/ - 7=У + 7Ч(Е/ + /)2 , , . _ г Т-aí?2 + 6(£/ + /)2-+ 7 <?(М) = С,

ГДС s = an2+b(e1+j)l' í = ar?+£j+f)*> а функции / И 0 имеют ТОТ ЖС СМЫСЛ, что и выше; постоянная С определяется формулой (2.15); тг = 1,2,3,...

Важно отметить, что дисперсионное уравнение (2.17) не зависит от решений системы (2.16).

Фактически уравнение (2.17) является семейством (но не системой) уравнений для различных п. Необходимо решать относительно 7 каждое из получающихся уравнений. Другими словами, пусть af¡ -- множество решений дисперсионного уравнения (2.17). Тогда = где <jj содержит все положительные решения (и только их) уравнения Ф(7; j) — h = 0. Кроме того, cr¿ П <jj = 0 для всех возможных i ф j.

Пусть - совокупность всех собственных значений задачи P\¡. Справедлива следующая теорема об эквивалентности, которая, в частности, утверждает, что ом = af¡.

Теорема 2.1 (об эквивалентности). Значение 7 = 7 является собственным значением задачи Рм тогда и только тогда, когда существует целое число п = ñ ^ 0 такое, что 7 = 7 удовлетворяет уравнению Ф(7;п) - h = 0.

Кроме того, собственная функция Z(x; 7) имеет в точности п нулей при х е (0 ,h); если x¡ является i-м нулем функции Z{x\ 7), то

+оо

xi = / wdri+ (г - 1)Т2.

ч(о)

Теорема 2.2 (о периодичности). Пусть 7 - собственное значение задает Рм. Если собственная функция Z(x; 7) имеет более двух нулей при х G (0, h), тогда функция Z(x; 7) периодическая с периодом 2

Пусть Г - (max(ei, £3),£f) и /¿fnf = iiif Ф(7,k), h* = sup Ф(у,к).

гег 7зеГ

Be

личина всегда, существует. Указанная sup существует не всегда, например, sup не существует, если / = 0, g = 0. Когда мы пишем /4ир, мы предполагаем, что указанный (конечный) sup существует.

Теорема 2.3. Пусть для некоторого р выполняется к'(п! < к?ир, / € С1 [0, +оо), д 6 С1 [0,+оо), система (2.16) имеет действительное решение (Х/,,^/,) и /г таково, что < к < кЦир, тогда задача Рм имеет по крайней мерс одно решение (собственное значение.).

Теорема 2.4. Пусть Ф(7, к) неог]хтиченна при 72 £ Г, / £ С1 [0, +ос), д Е С1 [0,+оо), система (2.16) имеет действительное решение (X/,, И/,) и к таково, что для некоторого р выполняется < к, тогда задача Р\1 имеет по крайней мере одно решение (собственное значение).

В теоремах 2.3, 2.4 утверждается существование собственных значений; в теоремах 2.5, 2.6 - их изолированность.

Пусть правые части системы (2.12) являются аналитическими функциями (как функции комплексных переменных (и, и) = (А2,.^2)) в некотором шаре 13(0, П) с центром в нуле и радиусом Л > 0. Шар выбран таким образом, что |2Х2/у2 + /| < е/ — Д, где Д > 0 - фиксированное достаточно малое число. В силу непрерывности функции / такое число Л всегда существует, поскольку /(0) = 0 и /' ограничена на В(О, Л). Теорема 2.5. Пусть для некоторого р выполняется к?п1 < /г£пр, система (2.16) имеет действительное решение (А/,, .£/,), щмвыс части системы (2.12) являются аналитичеекгши функциями (как функции комплексных переменных (и, и) = (А-2, 22)) в некотором шире В(0,П), где Я > 0, число М является наибольшим из модулей максимумов правых частей системы (2.12) и к ^ т^ таково, что < к < /гЦ?1ф, тогда задача Рм имеет по крайней .мере одно решение.

Кроме того, если условие /¡?пГ < справедливо для всехр, то множество собственных значений ам задачи Рм является дискретным на И = {7 € II '. 72 € Г}, т.е. па каждом отрезке I С. И содержмтся не более конечного числа (изолированных) собственных значений.

Теорема 2.6. Пусть Ф(-у,к) пеог]мничснпа при 72 6 Г, система (2.16) имеет действительное решение (А/,,.!?/,), пщвые части системы (2.12) являются аналитическими ф/ункциями (как функции комплексных переменных (и, и) = (X2, Я2)) в некотором шаре. Л(0, Л), где Л > 0, число М является наибольшим из модулей максимумов правых ■частей системы (2.12) и /1 ^ тщ таково, что для некоторого р выполняется к?п{ < к, тогда задача Рм имеет по крайней мере одно решение.

Кроме того, если условие к?п1 < к? справедливо для всехр, то множество собственных значений &м задает Рм является дискретным па В = {7 £ К : 72 £ Г}, т.е. па каждом отрезке / С О содерлсится не более конечного ''тела (изолированных) собственных значашй.

Величины hinl и h¡fup можно находить численно.

Заметим, что условие -у- < £f является точным. Действительно, для функций / = 0, д = 0 получаем линейную задачу. В такой линейной задаче необходимо 72 < e¡.

В главе 3 исследуется нелинейная двухнараметрпческая задача сопряжения на собственные значения для системы уравнений Максвелла, описывающая распространение связанных электромагнитных ТЕ-ТМ-волн в слое с ксрровской нелинейностью.

Рассмотрим электромагнитные волны

ew-»r = + Ел/е-;ы»(,

= НЕе-пк* + Нл/е--»(, 1 ' >

распространяющиеся вдоль границы однородного, изотропного, немагнитного диэлектрического слоя £ {{x,y,z) € R3 : —h íj х ^ h}, здесь

Ея= (0,ЕУ,0)>\ Елг= (Ех,0,Егу\

П,.: = (НХ,(1 Hz)'\ НЛ/= (0, Ну, 0)т - (3.2)

ТЕ-волны ТМ-волны

комплексные амплитуды; u¡e,ojm - круговые частоты; ( • )т - операция транспонирования.

Слой Е расположен в декартовой системе координат Oxyz между двумя полупространствами х < 0 и х > h. Полупространства заполнены изотропной немагнитной средой без источников и имеют постоянные диэлектрические проницаемости £ = £\ £q п £ = е^ ^ £ц соответственно, где £о > О - диэлектрическая проницаемость вакуума. Считаем, что всюду ц = fia, где ¡iо > 0 - магнитная проницаемость вакуума.

Диэлектрическая проницаемость £ волновода £ описывается законом Ксрра1, который имеет вид £ = е2 + а|ЕШьШЛ,|2, где £2 > max(ei,e3), а > 0.

Поля (3.1) удовлетворяют сйстеме уравнений rot (Еве-'"*' + Еме-^>') - í/íwbHec-^-' +

rot (НЕеЧи»* + = -ieus^Ee4""1 - i£UJhlЕд/е^'4, ( ' '

условию непрерывности касательных составляющих компонент поля на границах раздела сред х = —h, х = h и условию излучения на бесконечности: электромагнитное поле затухает как 0{\х\~1) при \х\ —> ос.

'Делоне II. Б. Взаимодействие лазерного излучения с веществом. — М.: Наука, 1989 ; Ландау JI. Д., Лифшиц Е. М. Теоретическая физика: учеб. пособ. для вузов: 10 т. Т. VIII. Электродинамика сплошных сред. - М.: Физматлит, 2001 ; Шен И. Р. Принципы нелинейной оптики. -М.: Наука, 198!).

Доказано, что компоненты ТЕ-ТМ-волн, распространяющихся вдоль границы слоя имеют представление

Ех = Е1(ж)е*7л'-, Еу = Ev(a:)e'7';", Ez = Ег(з;)е'^,

Нх = Вх(х)е^=, Ну = Яу(х)с^'=, Н: - [6Л)

где 7е, 7л/ ~ пара неизвестных (действительных) спектральных параметров (пара постоянных распространения связанной ТЕ-ТМ-волны).

Пусть А'5 = Подставляя ноля (3.2) с компонентами (3.4) в

систему (3.3), нормируя в соответствии с формулами х = кцх, ^ = fcojj, IE = 1м = £j = %,a = f0, обозначая X := iEx, У := Е„', Z := е!,

9 —2

т := oj^uj^ и опуская значок тильды, получаем

7л/- Z') = еХ,

ГЕУ - У" = теУ, (3.5)

%мХ' - Z" = eZ, где

х < —к;

s= <£2 + a(X2 + Y2 + Z2), -к^х^к;

x > к.

Функции X, Y, Z дифференцируемы так, что

Х(х) е С1 (-оо, 0] П С^О, к] П С1 [к, +ос), У (ж) 6 С^—оо, +оо) П С2(-оо, 0) П С2( 0, к) П С2{к, +оо), Z{x) е С(—оо, +оо) П С1 ( оо, 0] П С1 [0, к] П С1 [к, +оо) П

П С2(—оо, 0) П С2(0, /г) П С2(к, +оо).

Условия сопряжения для функций Л', У, Z следуют из условии непрерывности касательных компонент электромагнитного ноля и имеют вид

[Z' - 7.UA1U±/, = 0, [У]|Х=±Л = 0, [У']|Х==Л. = О, [Z\\x=±h = 0. (3.7)

Определение 3.1. Пару чисел (7е,7м) такую, что при заданных значениях Л'(0), У(0) существует не равный тождественно нулю вектор W = (Л', У, Z)T такой, что функции X, У, Z удовлетворяют системе уравнений (3.5), условиям (3.6), (3.7) и затухают как 0(|ж|-1) при |ж| —» оо, будем называть парным собственным значением. Вектор IV, который соответствует парному собственному значению, будем называть собственным вектором, а компоненты X(x;%,%j), Y(x;je,1m), Z(x;je,7m) вектора W - собственными функциями.

Задача Р: доказать существование парных собственных значений (7е, 7л/), удовлетворяющих определению 3.1.

Пусть кгЕ1 := 71-теи к2Е-л := Те-те3, к\п := 7л/к2ш := 7л/ Решения (линейной) системы (3.5) при х < —Н и х > /г в соответствии с условиями излучения имеют вид

' Х(х) = С^е^-^'»», У (я) = С^е-^-'1^, (3.8)

Х(х) = cj^c^*»», У (я) =

где c[jh\ CE h\ off и CE^ - постоянные интегрирования. Постоянные СЕ' предполагаются заданными (они отвечают значениям Х(0), Y{0) в определении 3.1), а постоянные С\,Н\ Сопределяются пз условий сопряжения.

Пусть fc| = Г£2 - 7|, к2, = £2 - и fx = (X2 + Y2 + Z2)X, fY = (X2 + У2 + Z2)Y, fz = (A'2 + Y2 + Z2)Z, тогда внутри слоя К систему (3.5) можно записать в виде

х = --умк^г' - akrffx, Y" + &|У = -ат/у, (3.9)

Z" + A^Z = -ae^khfz - ае^Ум/'х-

Обращая линейные части второго и третьего уравнений в системе (3.9), перейдем к интегральным уравнениям.

Пусть Le = ¿2 + к%, Lm = + k\j. Тогда функции Грина для следующих краевых задач

Г LeGe = —6{х — s), Г

I OxGe|х=_„ = OxGEU = 0; 11 \

LmGm = —5(х — s), gm\x=-h = GA,\x=h = 0

(3.10)

имеют вид

COskp;(x+fl) COS kf.;(s~h) ^ h

Ge(x,s)={ ^t^t,(3.11)

kKsu\2kKh ' s < x n,

sin км(x+li) sin км(s-h) ^ ,

n,,(rc.)-J kMSm2kMh ' ' /о 19")

UM{X,b) — ^ sinki,(x-h) sinkM(s+h) „__.

км sin 2/,'д//(

s < х < h.

Переходя от системы дифференциальных уравнений (3.9) к интегральным (используя найденные функции Грина и вторую формулу Грина) и затем используя условия сопряжения (3.7) в полученной системе

интегральных уравнений, получаем систему дисперсионных уравнений

ФвЫти) := = О,

ФЕЫЛМ) ■= С^кши^ги) - a^tolM) = о, ^ ) где

SE{h, 7е) ■= (к% - kEikE3) sin 2kEh - кЕ (кЕ1 + кЕЯ) cos 2kEh\ (3.14) f)Ai{h, 7л/) := (е^зА-л/ - е\кткм3) sin 2kÁIh-

-е2км [£\кш + е3кмi) cos 2kMh,

Qe{K 7e, 7л/) := т{кЕI cos 2/j£/í - kE sin 2kEh) cos kE(x + h) fYdx—

~ ткЕi /Д cos kE(x - h)fYdx, Qm{Kie,1m) '•= —2r¡l¡(e\k\i sin iksih - e2A\m cos 2kMh)fx(h -0)x x sin 2кмН + км{£\км sin 2k±\¡h — £2&лл cos 2kÁ\¡h) x x [iMkirfz sin км(х + h) - 7л//л- cos км{х + /¿)] dx~ - 2-yl¡eo км ifx{—h + 0) sin 2 kMh+ + еЖшкм /Д [lMkMfz sin k.,u(x - /г) - 7j,fx cos kM(x - /г)] dx. Доказано, что в интегральной форме система дифференциальных уравнений (3.9) может быть записана в виде

u = qNj (|u|2u) + h, (3.16)

где u = (X, Y, ZY и |u|2 = X2 + Y2 + Z2, \\ufc = \\Xfc + ||Y||2; + ||Z||2, и IHIc = max |u(x)|, а С[-/г,/г] = C[-h,h] x C[-h,hJ x C[-h,h\-

линейный оператор Ni имеет вид

N := qN] H Ni := К + J + К + J,

где

KS = f-h K(a;, s)g(x)dx, Kg = fhh K(x, s)g(x)dx, g = (.91, .92, <73 )T, матрицы ядер K(a;, s) и К(ж, s) имеют вид

(qnD2sGu 0 qi3daGni\ 0 q22GE 0

Q-JAGM 0 Q-MG.M J

9u = 4l"ñ¡ku, 9i3 = Чм, 922 = т, 9з1 = -e^Ju, Ф»з =

К (x,s) =

((/и lim dz,G \¡ 0 r/13 lim OsGm\

s-i-h+0 s->-/i+0

0 922 G£|s^_/i+0 0

<731 lim dxGu 0 (}i 3 lim Gm

Я-+-Л+0 s—»—/j+0

911 = 1mPi(s), Qi3 = -kl,pi{s), (¡22 = q(s), q31 = ~^P2(s), № = k\iP2{s)

.. „ /.Л _ ^..klLL_cosk-M(s-h)_ , ч _ А-д/i amkst(s-ll)

" Viy*) k'i, E-ifrA/sm2A-,uft-E-2fcjíi<-os2A-4,/i> W) — ьтЗАгд/Л-Ыдп cos2/.A/ft>

( \ A'l-i coakf. (s - h) T T

q(s) = tk¡:Án2kt:h_kmC0¡2kEh, а операторы J, J предсташшы следующим образом:

«,2 . _ /1 0 0\ fAipiis) 0 0\

j = _2íl+Ü о oo,j= о оо

\0 О 0/ \27л/^-/Р2(б') О О/

и h = (/ii, /¿2, /¿з) ', где

. _ СЙ'^Ч/з Г E2kmcosk„(a-h)__, , , ,ч"|

— kuam2kuh [^км sin2kMh-e2kM1 cos2kK,h ^ Lub ^ ^ '^J '

I _ fк/:?, Г COS kf;(s — h) I 7 / I ; ч"|

rt2 ~~ л [i,.;sin2i-t;a-te1 cos2k,.jh + cos (s + ft)j ,

. -cffí-ш Г_sinfcAf(s-/!)__, 4inI. /„ , /л]

"3 lMsm2kMh [скм^кмк- £2кшав2кы1г + bill ЛДДб -f fljj .

Операторы К, К, J, J являются линейными.

Теорема 3.1. Пусть В,ъ = {и : ||и|| ^ го} шар радиуса г о с центром в нуле. Также пусть выполняются два условия:

q := 3arg||Ni|| < 1, m-gHN^ + ||h|| ^ г0. (3.17)

Тогда существует единственное решение и 6 ВТа уравнения (3.16). Справедлива следующая теорема единственности. Теорема 3.2. Если а ^ А2, где А = п10г^а уравнение (3.16)

имеет единственное решение и в шаре Вг, = {и : ||и|| ^ г»} и и € 6 С [-h,h], ||и]| ^ г„ где г. = "^=cos (l arceos (§ ||h|| уЩЩ) - f).

Отметим, что Л > 0 не зависит от а.

Также справедлива теорема о непрерывной зависимости решения операторного уравнения (3.16) от спектрального параметра.

Теорема 3.3. Пусть ядра матричного оператора N и правые части h уравнения (3.16) непрерывно зависят от параметра 7 € Го, где Гц некоторый действительный отрезок. Пусть также ||h|| ^ Тогда

решение 11(7) уравнения (3.16) для 7 е Го существует, единственно и непрерывно зависит от параметра 7 G Го.

Переходя к пределу при а —» 0 в системе дисперсионных уравнении (3.13) для нелинейного слоя, получаем известные уравнения для линейных задач о распространении ТЕ- и TM-волн в слое с постоянной диэлектрической проницаемостью, а именно: .9я(/г, 7е) = 0, <ш(Л,7.и) = 0, где функции (/в(/г,7£), <ш(Л,7д/) определены формулами (3.14), (3.15) (эти дисперсионные уравнения, определяющие собственные значения, хорошо известны и полностью изучены).

Для собственных значений je и 7м указанных линейных задач справедливы неравенства:

inax(TE:i, тез) <j%< те2, max(¡rb е3) < jj¡ < е2. (3.18)

В задаче разыскиваются решения, удовлетворяющие этим же неравенствам. Обозначим

■2 /2 -2 -2

\Е .__Js,m ХЕ _ Jc.tn .Ц _ Js.tn .д/ Je,til

1 — T£°- 4/г2 ' Ас,"> '— Т£2 - . Km "— £2_ 4/¡2' Ac,m '— 62 ~ '

Ще js¡m = 7Гm - тп-й неотрицательный корень уравнения sin ж = О,

, • 7r(2m+l) „ _

ft Je,ш = —2— _ m'n положительны!! корень уравнения cos ж = 0; m = 0,1,2,...

Доказано, что если уравнение c¡E[h,jE) = 0 имеет 1е корней 7^', . ■ .,jgh\ то существует целое число 0 < тле < 1е такое, что

7Г} €

1е+1) 6 ' ¿ = rnE,lE-l.

Доказано, что если уравнение дм(И,7.1/) = 0 имеет 1М корней 7^', ) су] II

(¡+i)

... то существует целое число 0 < < l4\¡ такое, что

= пгм,1м ~ 1-

Мы предполагаем, что уравнения яе(Л.,7е) = 0 и д.и(/1,7м) = 0 имеют 1е и 1м корней соответственно (выбором /г всегда можно добиться того, что оба уравнения будут иметь корни, а одно из них будет иметь заданное число корней).

Точки ^ ^, г = 0,1е, ] = 0являются простыми полюсами функций Грина (3.11), (3.12). Пусть 5f > 0 и 6-' > 0, г = 0,1е,

j = 0, /д/, - достаточно малые числа такие, что локализация собственных значений 7^', г = 1,1е и 7^', j = 1,1м, линейных задач, указанная выше, не нарушается.

Построим отрезки

г = 0, ше — 1; I = гпЕ, 1е - 1; г = 0,тм - 1; г = гпм, 1м - 1-

При наших предположениях функция дЕ (7е) имеет различные знаки на концах отрезков Tf и обращается в нуль в точках 7^'. Тот же самый вывод имеет место для функции дм ('уд/)• Обозначим Г : — Гь х Гл/, где

:= 11!=!^, г" := и!=1 г£'.

Основным результатом третьей главы является следующая

Теорема 3.4. Пусть е2 > гпах(£1,£3) > О и уравнения <?£(/1,7в) = О, 5л/(Л, 7л/) = 0 имеют 1е и 1м решений соответственно. Тогда найдется значение а0 > 0 такое, что для всякого 0 < а ^ «о существует по крайней мере 1е-1м парных собственных значений Е^хГ^,

г = 1, 1е, j = 1,1м, задает Р.

Из теоремы 3.4 следует, что при указанных предположениях существуют поверхностные связанные ТЕ-ТМ-волны, распространяющиеся вдоль границы изотропного слоя с керровской нелинейностью. Так как в линейном слое связанных волн не существует, то из утверждения теоремы следует, что в нелинейном слое существует новый режим распространения волн. При доказательстве теоремы 3.4 получена теоретическая оценка величины »о через нормы соответствующих операторов.

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

Публикации по теме диссертации

Монографии

[1] Валовик Д. В., Смирнов Ю. Г. Распространение электромагнитных волн в нелинейных слоистых средах. Пенза : Изд-no ПГУ, 2010. 264 с.

[2] Smirnov Yu. G., Valovik D. V. Electromagnetic wave propagation in nonlinear layered waveguide structures. - Penza : PSU Press, 2011. — 248 p.

Публикации в журналах, рекомендованных ВАК

[3] Валовик Д. В., Смирнов Ю. Г. К задаче о распространении нелинейных связанных электромагнитных ТЕ-ТМ-волн в слое // Жур. вычнел. матем. и матсм. физ. - 2011. - Т. 54, № 3. - С. 504-518.

[4] Smirnov Yu. G., Valovik D. V. Problem of nonlinear coupled electromagnetic TE-TE wave propagation // J. Math. Phys. - 2013. - V. 54, 8. - P. 083502-1-13.

[5| Валовик Д. В., Смолькин Е. Ю. Расчет постоянных распространения неоднородных нелинейных двухслойных круглых цилиндрических волноводов методом задачи Коши // Радиотехн. и электроника. - 2013. - Т. 58, № 8. - С. 759-707.

[6] Валовик Д. В., Смирнов Ю. Г., Смолькин Е. Ю. Нелинейная задача сопряжения на собственные значения, описывающая распространение ТЕ-волн в двухслойных цилиндрических диэлектрических волноводах // Жур. вычнел. магем. и матсм. физ. - 2013. - Т. 53, № 7. - С. 1150-1161.

[7] Smirnov Yu. G., Valovik D. V. Coupled electromagnetic transvcrsc-elcctric-transverse magnetic wave propagation in a cylindrical waveguide with Kerr nonlinearity // J. Math. Phys. - 2013. - V. 54, .V 4. - P. 043500-1-22.

[8] Valovik D. V. On the problem of nonlinear coupled clcctromagnctic transverse electric-transverse magnetic wave propagation // Л. Math. Phys. - 2013. - V. 54, № 4. -P. 042902-1-14.

[9] Валовик Д. В., Зарембо Е. В. Решение нелинейной краевой задачи на собственные значения для электромагнитных ТМ-волн, распространяющихся в слое с керровской нелинейностью, методом задачи Коши // Радиотехн. и электроника. -2013. - Т. 58, Л'' 1. - С. 69-72.

[10| Валовик Д. В., Зарембо Е. В. Метод задачи Коши решения нелинейной задачи сопряжения на собственные значения для ТМ-волн, распространяющихся в слое с произвольной нелинейностью /,/ Жур. вычнел. матем. и матсм. физ. - 2013. -Т. 53, № 1. - С. 74 89.

[11] Smirnov Yu. G., Valovik D. V. Coupled electromagnetic TE-TM wave propagation in a layer with Kerr nonlinearity //J. Math. Phys. - 2012. - V. 53, S' 12. - P. 123530-1-24.

[12) Валовик Д. В., Смирнов Ю. Г. О распространении связанных электромагнитных ТЕ- и ТМ-волн в плоском слое с керровской нелинейностью // Известия высших учебных заведений. Поволжский регион. Физ.-матем. науки. - 2012. - № 4(24). -С. 21-48.

[13] Валовик Д. В., Смолькин Е. Ю. Численное решение задачи о распространении электромагнитных ТМ-волн в круглом диэлектрическом волноводе, заполненном нелинейной средой // Известия высших учебных заведений. Поволжский регион. Физ.-матем. пауки. - 2012. - № 3(23)'. - С. 29-37.

|14] Валовик Д. В. Распространение электромагнитных ТЕ-волн в нелинейной среде с насыщением // Радиотехп. и электроника. - 2011. - Т. 56, X® 11. - С. 1329-1335.

[15] Валовик Д. В. Задача о распространении электромагнитных ТМ-волн в слое с произвольной нелинейностью // Жур. вычисл. матсм. и матем. физ. -2011. - Т. 51, 9. - С. 1729-1739.

[1С] Валовик Д. В., Смирнов Ю. Г. Нелинейные эффекты в '¡адате о распространении электромагнитных ТМ-волн в слое с ксрровской нелинейностью // Радиотехн. и электроника. - 2011. - Т. 56, № 3. - С. 309-314.

[17) Валовик Д. В., Смирнов Ю. Г. Распространение ТМ-ноляризованных электромагнитных волн в диэлектрическом слое из нелинейного метаматериала // Известия высших учебных заведений. Поволжский регион. Физ.-ыатсм. науки. - 2010. -■V 3(15). - С. 71-87.

[18) Валовик Д. В. Задача о распространении электромагнитных волн в слое с произвольной нелинейностью (II. ТМ-волны) // Известия высших учебных заведений. Поволжский регион. Физ.-матем. науки. - 2010. - № 2(14). - С. 54-65.

[19) Валовик Д. В. Задача о распространении электромагнитных волн в слое с произвольной нелинейностью (I. ТЕ-волны) // Известия высших учебных заведений. Пополжский регион. Физ.-матем. науки. - 2010. - № 1(13). - С. 18-27.

[20) Валовик Д. В., Смирнов Ю. Г. Расчет постоянных распространения и нолей для поляризованных электромагнитных ТМ-волн в нелинейном анизотропном слое // Радиотехн. и электроника. - 2009. - Т. 54, Л'! 4. - С. 411-417.

[21) Валовик Д. В., Смирнов Ю. Г. О распространении ТМ-иоляризованных электромагнитных волн в нелинейном слое с нелинейностью, выраженной законом Керра // Жур. вычисл. матсм. и матсм. физ. - 2008. - Т. 48, № 12. - С. 2185-2194.

)22) Валовик Д. В., Смирнов Ю. Г. Нелинейная задача на собственные значения для ТМ-иоляризованных электромагнитных волн в нелинейном слое // Изв. вузов. Математика. - 2008. - № 10. - С. 70-74.

[23] Валовик Д. В., Смирнов Ю. Г. Расчет постоянных распространения ТМ-иоля-ризованных электромагнитных волн в нелинейном слое // Радиотехн. и электроника. - 2008. - Т. 53, .V 8. - С. 934-940.

[24) Валовик Д. В. О существовании решений нелинейной краевой задачи на собственные значения для ТМ-иоляризованных электромагнитных волн // Известия высших учебных заведений. Поволжский регион. Физ.-матем. науки. - 2008. - № 2. -С. 86-94.

[25] Валовик Д. В., Смирнов Ю. Г. Распространенно ТМ-ноляризованных электромагнитных волн в нелинейном анизотропном слое // Известия высших учебных заведений. Поволжский регион. Физ.-матем. пауки. - 2007. - № 4. - С. 51-59.

Публикации в других журналах

[26| Smirnov Yu. G., Valovik D. V. Boundary eigenvalue problem for Maxwell equations in a nonlinear dielectric layer // Appl. Math. - 2010. - № 1. - P. 29-36.

[27] Smirnov Yu. G., Valovik D. V. Nonlinear effects of electromagnetic TM wave propagation in anisotropic layer with Kerr nonlincaritv // Adv. Math. Pliys. - 2012. -V. 2012. - P. 1-21.

[28] Valovik D. V. Electromagnetic wave propagation in nonlinear layered waveguide structures. Computational approach to determine propagation constants // Springer Proceedings in Mathematics & Statistics, Volume 52: Inverse Problems and Large Scale computations. - 2013. - P. 71-91.

[29] Smirnov Yu. G., Smol'kin E. Yu., Valovik D. V. Nonlinear Double-Layer Bragg Waveguide: Analytical and Numerical Approaches to Investigate Wavcguiding Problem // Adv. in Numer. Anal. 2014. V. 2014. P. 1-11.

Научное издание

Нелинейные одно-и двухпараметрическне

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

Специальность 01.01.02 Диффереициальиые уравнения, динамические системы и оптимальное управление

Редактор А. Г. Тсыникова Технический редактор А. Г. Тсыникова Компьютерная верстка Д. В. Валовика

Подписано в печать 04.07.2014. Формат 60x84-^. Усл. нсч. л. 1,03. Заказ № 008474. Тираж 100.

Пенза, Красная, 40, Издательство ПГУ Тел./факс: (8412) 56-47-33; e-mail: iic@pnzgu.ru