Задачи распространения электромагнитных волн в различных средах были и остаются актуальными в связи с их широким практическим применением. Необходимость теоретического исследования существования и свойств собственных волн диктуется практической потребностью передачи энергии поля на большие расстояния с минимальными потерями. Успехи в разработке данного направления электродинамики привели к построению различных классов волноведущих структур.
Распространение электромагнитных волн в волноводах с заполнением линейной средой (то есть когда диэлектрическая и магнитная проницаемости не зависят от электромагнитного поля) -тема классической электродинамики [9], [28], [29], [30].
В случае волновода кругового сечения и постоянных электрической и магнитной проницаемостей уравнения Максвелла решаются в цилиндрических координатах, при использовании метода разделения переменных появляется линейное обыкновенное дифференциальное уравнение, называемое уравнением цилиндрических функций или уравнением Бесселя, решение которого является комбинацией цилиндрических функций. Собственные функции и собственные значения определяются как решения краевых задач с дополнительными условиями на контуре для решений и их первых производных [9]. С появлением нелинейной оптики предметом изучения в электродинамике стали сильные волновые поля, в которых начинает проявляться нелинейность сред. Качественно новыми эффектами нелинейной оптики стали порождение средой высших гармоник, а также самовоздействие волнового процесса, распространяющегося в нелинейной среде [9]. При распространении через жидкость или газ волны, создаваемой лазером, учитывают нелинейность поляризации среды, вызываемую целым рядом факторов [9]. Помимо поведения электронов в сильном электромагнитном поле существенно механическое воздействие поля на вещество: возникает давление, пропорциональное средней мощности волны, в результате чего в областях сгущения увеличивается диэлектрическая проницаемость. При распространении резко неоднородной волны — луча лазераможно сказать, увеличивается оптическая плотность среды в области сильного поля. Иными словами, в определенных условиях волновому процессу сопутствует образование канала, направляющего его энергию — нечто вроде диэлектрического волновода. Это называется самоканализацией, если канал сужается, наступает так называемая самофокусировка [9].
Исследования данной тематики претерпели значительную эволюцию в течение последних 50 лет, с того времени, как было обнаружено явление самофокусировки электромагнитного поля в нелинейной среде. Впоследствии основным предметом исследователей стали волноводы различных конфигураций и с различным заполнением среды. Эффекты самофокусировки используются в лазерах и оптоэлектронных приборах, для построения устройств обработки сигнала или блокировки моды в волоконных лазерах.
Квантовая теория самофокусировки была разработана Талановым для плазмы и группой американских ученых R.Y.Chiao, Е. Garmire, С.Н. Townes [22] для твердых тел в 1964 году, где исследуются условия, при которых электромагнитный луч может распространяться без затухания в средах с диэлектрической проницаемостью, возрастающей пропорционально интенсивности электромагнитного поля, и постоянной в его отсутствие. В работе [22] представлено теоретическое объяснение самофокусировки. Если при воздействии электромагнитного поля диэлектрическая постоянная обуславливает критический угол внутреннего отражения, превышающий дифракционный угол отклонения луча, то рассеяние путем дифракции не будет иметь места. Таким образом, при мощности луча (уменьшающейся пропорционально квадрату длины волны), выше некоторого критического уровня этот луч может быть сфокусирован при любом произвольном радиусе. В обычных диэлектрических материалах для достижения самофокусировки достаточно мощности обычного лазерного светового луча, т. е. порядка 108 Вт.
Авторами [22] получено точное решение волнового уравнения в форме линейно поляризованного плоского светового луча и в форме цилиндрически поляризованного луча кругового сечения с поперечной составляющей электромагнитного поля .
Исследование [22] вызвало большой интерес и послужило толчком к началу многочисленных исследований. Так, например, в работе [26] советских авторов Д. И. Абакарова, А. А. Акопяна и С. И. Пекара рассматривается новый вид самофокусирующегося луча, на этот раз с поперечным сечением произвольной формы и поперечными размерами, значительно превосходящими длину световой волны. Применение численных методов в этом случае в нулевом приближении дает результат [22], а в более высоком приближении констатируется появление продольного электромагнитного поля.
Расслоение среды на области, в одних из которых электромагнитное поле поперечно, а в других появляется также и продольная составляющая, было обнаружено также В. И. Елеонским и В. П. Силиным [20] при исследовании сходной задачи отражения наклонно падающих р — поляризованных волн от среды с нелинейной диэлектрической проницаемостью. Рассматривается отражение плоско-поляризованных волн (вектор электрической составляющей которых лежит в плоскости падения) от полупространства, заполненного средой с нелинейной диэлектрической проницаемостью. Чередование областей с поперечной и продольной составляющими поля объясняется чередованием продольной и поперечной степеней свободы электромагнитного поля в нелинейной среде.
Наряду с рассматриваемыми проблемами плоской геометрии, В. И. Елеонский представил теорию цилиндрических самофокусирующихся волноводов, в основном базирующуюся на качественном анализе фазовых траекторий [25], а также численные расчеты основной и низшей неосновной ТЕ-мод.
Y. Chen, основываясь на вариационной технике, получил аналитические решения для фундаментальной моды в свободном пространстве в нелинейных средах, используя аппроксимацию Гауссовой функцией [17], [24].
R.A.Sammut, C.Pask. [15] использовали вариационную формулировку волнового уравнения для волноводов с произвольной нелинейностью оптического волокна и представили аналитическую аппроксимацию для поля.
Позже D. Sjoberg [23] анализировал распространение электромагнитных волн в нелинейных цилиндрических волноводах методом возмущений, используя коэффициент нелинейности как параметр возмущения.
Для моделирования распространения волн в средах с нелинейностью Керра также используются методы распространяющегося луча [37], метод линий [38], конечно — разностные методы [8].
Наряду с вышеописанными разработками существуют также направления, посвященные проблемам устойчивости. Изучение устойчивости распространения волн необходимо при исследовании возможности практического применения эффекта самофокусировки. Важными в этой области являются исследования Н. Н. Ахмедиева [5], [6], [18], которые посвящены нестационарной проблеме нестабильности модуляции основной моды цилиндрического волновода с нелинейным заполнением среды по закону Керра. Эффект нестабильности модуляции в одномодовом оптическом волокне к этому времени был уже предсказан, исследован теоретически и практически наблюдаем [31] - [36]. Рассматривается распространение высокои низкочастотного поля в волноводе круглого сечения, диэлектрическая проницаемость среды внутри и вне волновода описывается законом Керра. В одномодовом оптическом волокне при нелинейных эффектах поперечная структура пучка определяется самим волноводом в линейном режиме, а продольная функция для относительно низких частот определена сочетанием нелинейности и дисперсии волокна. Для высоких частот или больших радиусов волновода должна также быть учтена дифракция, и ее влияние при высоких частотах на нестабильность модуляции исследовано автором. Эффект нестабильности модуляции вызывает проявление отдельных световых групп, однако, дальнейшее поведение этих групп является открытым вопросом. Наиболее их вероятное поведение в среде с нелинейностью по закону Керра — это пространственно — временное затухание или коллапс. Таким образом, в результате нестабильности модуляции на выходе оптического волокна имеем затухание в пространстве и времени вибраций. В среде с насыщенной нелинейностью, где также существуют подобные эффекты, световые группы при определенных условиях могут распространяться на достаточно длинные дистанции в виде трехмерных стационарных солитонов.
Вопросу нелинейных волновых взаимодействий посвящена работа А. П. Сухорукова [60]. Распространение волн в нелинейной среде в слое изучала группа исследователей под руководством В. П. Шестопалова [54] — [59], сходной проблемой о плоском трехслойном волноводе занимались H.W. Shurmann, B.C. Серов и Ю. В. Шестопалов [1] — [4], в работах которых описываются отражение и распространение плоской ТЕ-поляризованной волны. Аналогично случаю, рассматриваемому в настоящей диссертационной работе, ищутся решения уравнений Максвелла, удовлетворяющие условиям непрерывности на границе раздела сред. При этом также используется метод разделения переменных.
Представлено общее аналитическое решение уравнения Гельмгольца, описывающее рассеяние плоской монохроматичной волны в слое с нелинейностью по закону Керра. Все среды предполагаются непоглощающими, немагнитными, изотропными и однородными. Полученные результаты содержат условия существования физически допустимых решений, исследование зависимости решений от различных параметров задачи, проведено сравнение с линейным случаем.
Изучая распространение волн в слое, разные авторы использовали различные подходы, вводя специальные условия по отношению к нелинейной системе Фабри — Перота: J.H.Marburger и.
F.S.Felber [39] упростили анализ путем введения граничных условий, предполагающих, что нелинейный слой отделен от линейных сред идеальными зеркаламиJ. Danikaert [40] рассматривает реакцию твердого резонатора Фабри — Перота, включающую нелинейное поглощение и наклонное падение поперечно — электрического и поперечно — магнитного полейM. Haeltermann [41] и G. Vitrant [42] представили объединенную нелинейную теорию для поперечных эффектов резонаторов ФабриПерота, упростив численные расчеты и хорошо объяснив понятие оптической бистабильности.
Решения нелинейных уравнений Гельмгольца были представлены частично, для специальных случаев) в виде эллиптических функций Якоби [43] — [47]. Как показано в работе одного из авторов [48], общее решение более эффективно может быть представлено в форме эллиптических функций Вейерштрасса. Хорошо известно [14], [27], что эллиптические функции, в общем случае, обладают полюсами, однако, в предыдущих исследованиях е общее решение может быть использовано непосредственно для определения критических значений коэффициента отражения в различных случаях. Также отсутствовал анализ неограниченных решений в случае отрицательной константы Керра, в частности, никаких условий для существования таких решений.
В исследованиях [1] - [4] представлено аналитическое решение нелинейного уравнения Гельмгольца в слое, которое используется при анализе физической проблемы. Показано, как линейный случай появляется из нелинейного.
В работах [1] - [4] также проводится анализ так называемой «отсечки». Это случай неограниченной напряженности поля, и который особо выделяется с физической точки зрения, так как ему нет аналогов в линейной оптике. Необходимым условием для существования неограниченной напряженности поля является отрицательное значение коэффициента нелинейности. Как следствие, диэлектрическая проницаемость может равняться нулю внутри слоя. В постановке других физических задач явление неограниченной напряженности поля хорошо известно. По отношению к рассматриваемой ситуации случай неограниченного поля является, скорее всего, артефактом [4]. Существует вероятность, что эта особенность пропадет в случае введения в диэлектрическую проницаемость условий поглощения, но, с другой стороны, нельзя назвать такую диэлектрическую проницаемость нефизичной или искусственной, так как она дает начало многим хорошо известным явлениям нелинейной оптики. С математической точки зрения случай неограниченности поля не нуждается в особом пояснении, так как нелинейные обыкновенные дифференциальные уравнения, к которым сводится задача, обладают особенностью в полюсах. Тем не менее, физическая интерпретация «почти неограниченного» поля, например, для слабо поглощающих сред, пока не изучена [4].
Проводя сравнительный анализ этой работы с предыдущими исследованиями [50], [45], [46], можно отметить, что построение диэлектрической проницаемости 8, процедура решения нелинейного уравнения Гельмгольца и формулировка граничных условий являются традиционными, в то время как подход представления решений в виде эллиптических функций Вейерштрасса вместо использующихся ранее Якобиановых является более эффективным и имеет явное преимущество, описывая более общий физический случай.
Тем не менее, и данный подход оказался не универсальным, аналитического решения для общей задачи не было получено. Применение этого метода к изучению задачи о распространении волн в структурах кругового сечения приводит к неразрешимым дифференциальным уравнениям. Изучение задачи в такой постановке потребовало разработки нового математического аппарата. Кроме того, при всем многообразии проведенных исследований, дисперсионные соотношения теоретически не исследованы, не выяснены до конца вопросы, связанные со сходимостью методов, нет базовых теоретических результатов о существовании и единственности решений, поэтому данная диссертационная работа, в которой проведены такого рода исследования, является актуальной и практически важной.
Работа посвящена решению нелинейной краевой задачи на собственные значения для нелинейного уравнения Гельмгольца, к постановке которой приводит вопрос об изучении распространения собственных волн в цилиндрических диэлектрических волноводах кругового сечения, заполненных средой с нелинейностью по закону Керра. Такого рода задачи так или иначе сводятся к изучению сложной оператор — функции, нелинейно зависящей от спектрального параметра, которое весьма затруднительно традиционными методами теории дифракции. В связи с этим перспективным является метод оператор — функций, впервые примененный к задачам электродинамики в работах А. С. Ильинского и Ю. В. Шестопалова [51] - [53]. Для сведения задачи о собственных волнах к задаче на собственные значения с нелинейным вхождением спектрального параметра для интегральных уравнений используется метод функции Грина. Впервые для доказательства существования решений краевой задачи на собственные значения используются методы, базирующиеся на теореме Шаудера и принципе сжимающих отображений. Следует отметить, что такие методы использовались ранее для решения интегральных уравнений [12], но не были использованы в целом при решении задачи на собственные значения. При использовании предлагаемого метода спетральный параметр X изначально не фиксируется, а для каждого его значения решается уравнение краевой задачи, а затем полученное решение подставляется в дисперсионное соотношение. Принцип сжимающих отображений применяется для доказательства существования решений интегральных уравнений. В качестве численного метода предлагается итерационый метод, сходимость которого теоретически обоснована в работе.
Итак, основными целями настоящей работы являются:
Строгая постановка задачи о распространении собственных волн в цилиндрических диэлектрических волноводах круглого сечения с нелинейным заполнением среды по закону Керра как краевой задачи на собственные значения для системы уравнений Максвелла.
Разработка математического аппарата для исследования задачи о собственных волнахдоказательство базовых теоретических результатов о существовании и единственности решений дисперсионных уравнений (относительно спектрального параметра) и интегральных уравнений, отвечающих краевой задаче (относительно собственных функций).
Построение, обоснование и реализация эффективных численных методов для расчета собственных значений и соответствующих им собственных функций для поставленной задачи.
Работа состоит из трех глав и двух приложений.
Первая глава посвящена постановке задачи. Описывается класс волноведущих структур, которые будут рассмотрены в работе, формулируется задача о собственных волнах для однородной системы уравнений Максвелла. Для рассматриваемых двух случаев ТЕи ТНполяризации распространяющихся волн используется единый алгоритм сведения поставленной краевой задачи для уравнений Гельмгольца к нелинейной краевой задаче на собственные значения для обыкновенных дифференциальных уравнений с помощью метода разделения переменных.
Вопрос о точных решениях полученных нелинейных обыкновенных дифференциальных уравнений вынесен на отдельное рассмотрение в Приложение 1.
Вторая глава, посвященная исследованию разрешимости задач на собственные значения, содержит основные теоретические результаты. В первом пункте главы с использованием метода нелинейной оператор — функции строится функция Грина для краевой задачи и исследуются ее свойства. Во втором пункте выводятся интегральные уравнения, отвечающие краевой задаче, с помощью построенной функции Грина. Третий пункт содержит утверждения о существовании и единственности решений интегральных уравнений краевой задачи, опирающиеся на принцип Шаудера. Приведена и доказана соответствующая теорема, в которой также дается оценка коэффициента нелинейности, показывающая границы применимости предложенного метода. Четвертый пункт решает вопрос о непрерывной зависимости решения краевой задачи от спектрального параметра. Это утверждение будет одним из ключевых при доказательстве существования нетривиального решения поставленной задачи.
Основной теоретический результат содержится в пятом пункте второй главы, где, после приведения дисперсионного соотношения к нормализованному виду, даются достаточные условия существования ненулевых решений дисперсионного уравнения относительно спектрального параметра. Таким образом, при определенных теоремой условиях, существуют осесимметричные ТЕполяризованные волны, распространяющиеся без затухания в цилиндрических диэлектрических волноводах кругового сечения, заполненных средой с нелинейностью, выраженной законом Керра.
Это утверждение является обобщением известного соответствующего результата в линейной оптике. В этом же пункте установлено, что спектр рассматриваемой краевой задачи дискретен, границы каждого из m интервалов, внутри которого существует, по крайней мере, одно решение, определены с использованием ш-ых положительных корней функций Бесселя J0, Ji и Jb, а условия распространения волн без затухания зависят не только от малости параметра нелинейности, но также и от радиуса и вещественного параметра волновода.
В шестом и седьмом пунктах главы, после того, как установлена непустота спектра, доказывается существование приближенных решений параметра на каждом из рассматриваемых промежутков и корректность определения краевой задачи в п-ом приближении как функции от соответствующего решения спектрального параметра в том же приближении. Далее приводится формулировка численного метода решения интегрального уравнения краевой задачи, в качестве которого используется итерационный метод, и доказывается его сходимость.
При реализации численных расчетов необходимо знать границы выбора параметров, при которых выполняются условия полученных выше теоретических результатов, в первую очередь, теоремы о существовании решения решений дисперсионного уравнения. Заключительный восьмой пункт второй главы содержит некоторые оценки операторов задачи на собственные значения.
Третья глава посвящена разработке и реализации предложенного итерационного метода, графические результаты расчетов вынесены в Приложение 2.
В первом пункте третьей главы описывается алгоритм решения задачи на собственные значения в нулевом и первом приближениях, причем каждое действие алгоритма обосновано в силу доказанности базовых теоретических результатов о существовании точных и приближенных решений задачи, а также сходимости приближенных решений к точным. Начиная с нулевой итерации собственной функции Uo в формуле итерационного процесса, используется последовательность дисперсионных соотношений (решения которых, являющиеся приближениями собственных значений, сходятся к точным решениям) для получения соответствующего нулевого приближения спектрального параметра. Далее полученные нулевые приближения собственного значения и собственной функции применяются для следующего шага итерационного процесса, то есть первого приближения собственной функции иь и процесс повторяется.
Достоинством этого метода является хорошее приближение собственной функции и собственного значения уже в нулевом приближении (сравнение производится с решением задачи полным методом). Недостаток такого подхода — чрезмерное усложнение дисперсионного соотношения уже на первой итерации.
Во втором пункте третьей главы предложен альтернативный метод полного решения краевой задачи на собственные значения с заданной точностью. Сначала выбирается интервал поиска собственных значений с учетом условий соответствующих теорем, потом определяется шаг деления рассматриваемого интервала и вычисляются промежуточные собственные значения в каждом узле. Для каждого такого значения параметра решается итерационное уравнение, пока собственная функция не будет определена с требующейся точностью. После этого каждая пара собственного значения и соответствующей собственной функции подставляется в дисперсионное уравнение на предмет проверки знака последнего в узле. Заключительным процессом является нахождение интервала разбиения, на концах которого дисперсионное уравнение меняет знак. Руководствуясь теоремой о непрерывной зависимости от параметра, полагаем корнем среднюю точку отрезка.
Заключение
подытоживает основные результаты работы и указывает на возможности применения предложенного подхода в смежных областях.
В заключение приводится список литературы, состоящий из 61 источника, в том числе 8 работ автора работы, включенных в труды конференций и семинаров, в том числе и международных.
4 Заключение.
Из теоремы 2 следует, что при условиях, сформулированных в пункте 3.3, существует осесимметричное распространение ТЕ-поляризованных волн без затухания в цилиндрических диэлектрических волноводах кругового сечения, заполненных идеально проводящей, немагнитной, изотропной средой с диэлектрической нелинейностью, выраженной законом Керра. Этот результат обобщает известное соответствующее утверждение для диэлектрических волноводов круглого сечения с заполнением линейной средой. В ходе исследования были достигнуты следующие основные результаты:
1. Представлена строгая постановка задачи распространения собственных волн в цилиндрических диэлектрических волноводах с нелинейным заполнением средой по закону Керра, задача сведена к краевой задаче на собственные значения для уравнений Максвелла.
2. Получены все необходимые теоретические результаты для исследования распространения собственных волн: а). Доказаны теоремы о существовании точных решений задачи на собственные значения. б). Доказаны теоремы о существовании приближенных решений задачи на собственные значения. в). Доказаны теоремы о сходимости приближенных решений к точным.
3. Предложен и реализован численный метод расчета спектрального параметра и соответствующей собственной функции.
4. На основании анализа дисперсионного соотношения возможна продуктивная разработка анализа стабильности различных мод. Известно, что интеграл потока энергии Р = / dpp3z (р, 7) (где обозначает усрененную по времени z-компоненту ветора Пойтин-га) может быть использован в некоторых случаях для анализа стабильности. z может быть вычислено посредством итерационных решений un (7,/?).
5. Посредством предложенного подхода, который является достаточно обобщенным, могут быть исследованы другие виды нелинейности (например, насыщение, нелинейность более высокого порядка и т. д.).