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

На правах рукописи
СОЛОВЬЕВ АНДРЕЙ СЕРГЕЕВИЧ
ПРИМЕНЕНИЕ МЕТОДОВ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ И
ПЛОСКИХ ВОЛН ДЛЯ РАСЧЕТА ДИФРАКЦИИ НА
ДИЭЛЕКТРИЧЕСКИХ СТЕРЖНЯХ И ПОИСКА
СОБСТВЕННЫХ ВОЛН УСИЛИВАЮЩИХ
МИКРОСТРУКТУРНЫХ ВОЛОКОН
01.04.21 – Лазерная физика
Автореферат
диссертации на соискание ученой степени
кандидата физико-математических наук
Саратов – 2008
Работа выполнена на кафедре лазерной и компьютерной физики Саратовского
государственного университета им. Н.Г. Чернышевского (г. Саратов, Россия)
Научный руководитель:
доктор физико-математических наук, профессор
Леонид Аркадьевич Мельников
Официальные оппоненты:
доктор физико-математических наук, профессор
Давидович Михаил Владимирович
кандидат физико-математических наук, старший
научный сотрудник ИРЭ РАН
Морозов Юрий Александрович
Ведущая организация:
Научный центр волоконной оптики при институте
общей физики РАН, г. Москва
Защита состоится «28» октября 2008 г. в 17.30 на заседании диссертационного совета
Д212.243.05 по специальности 01.04.21 – Лазерная физика при Саратовском государственном
университете им. Н.Г. Чернышевского по адресу: 410012, г. Саратов, ул. Астраханская, д.83,
корп. 3, ауд. 34.
С диссертацией можно ознакомиться в Зональной научной библиотеке им. В.А. Артисевич
Саратовского государственного университета
Автореферат разослан 26 сентября 2008 г.
Ученый секретарь диссертационного совета
В.Л. Дербов
Общая характеристика работы
Актуальность темы. Задачи дифракции электромагнитной волны на диэлектрической решетке актуальны для многих областей физики и техники, и являются стандартными задачами электродинамики, решение которых с использованием различных
методов можно найти в большом числе публикаций. Периодические решетки находят
широкое применение в различных областях науки и техники: спектроскопии и антенной технике, физике твердого тела, милли- и субмиллиметровой радиотехнике. Анализ
физических процессов, протекающих в дифракционных решетках, связан с решением
важных прикладных задач: созданием частотных и поляризационных фильтров, новых
источников излучения, дисперсионных открытых резонаторов и т.д.
Особую роль играют теоретические исследования, значительно сокращающие или
полностью исключающие экспериментальную проработку. С их помощью можно осуществить модельный синтез и оптимизацию различных устройств. Частично эта актуальная проблема решена в рамках математической теории дифракции плоских волн на
дифракционной решетке, позволившей эффективно проанализировать ряд практически
интересных ситуаций, накопить большой объем данных об аномальных и резонансных
режимах рассеяния, решить некоторые малопараметрические обратные задачи. Для решеток со сложной структурой решение задачи дифракции возможно лишь численно: при
этом эффективность алгоритмов и сходимость методов являются крайне важными.
Для решения проблемы дифракции на диэлектрических периодических структурах
используются различные численные и полуаналитические методы. Оптимальным представляется разработка и использование методов, наиболее удобных для решения целого
класса задач. В данной работе такими задачами являются задачи дифракции лазерного монохроматического поля на периодических диэлектрических структурах со сложной
формой поперечного сечения образующих ее элементов (решетки и фотонные кристаллы) или распространения волн в таких структурах (микроструктурные волокна).
Одним из наиболее широко применяемых методов расчета является метод, использующий ряды Фурье и разложения по плоским волнам, который в применении к анализу
дифракции на периодических структурах рассмотрен, например, в статье1 . В этом слу1
Lifeng Li. Use of Fourier series in the analysis of discontinuous periodic structure. // J. Opt. Soc. Am. A.
1996. V. 13, No. 9
3
чае необходимо разлагать в ряд все периодические функции, в том числе и диэлектрическую проницаемость. Полученная система уравнений решается относительно неизвестных амплитуд пространственных гармоник. Существуют также варианты матричных методов; так, например, был предложен метод матричной рекурсии для моделирования процесса дифракции на диэлектрических решетках2 . Решетка с неоднородностями
произвольной формы представляется в виде множества слоев, лежащих друг на друге,
каждый из которых представляет собой периодическую структуру с неоднородностями
прямоугольной формы. При переходе от одного слоя к другому используют матрицу передачи. Эти методы, при всех их достоинствах, имеют и ряд недостатков (проблемы со
сходимостью и трудности с реализацией для элементов со сложным профилем поперечного сечения).
В последнее время большое внимание уделяется также новому типу оптических материалов - фотонным кристаллам. Главной особенностью фотонно-кристаллических
структур является наличие в них запрещенных частотных зон, запрещающих распространение внутри кристалла электромагнитных волн с частотами лежащими в этой зоне.
Это свойство используется исследователями для проектирования и создания различных оптических устройств, таких как оптические волокна3 , одномодовые волноводы4
и лазеры5 .
В последнее время наряду со значительным прорывом в области технической реализации процесса создания фотонных кристаллов, был разработан целый ряд различных
математических инструментов и методов, позволяющих смоделировать процесс распространения электромагнитной волны в фотонных кристаллах, рассчитать и проанализировать их дисперсионные характеристики. Здесь также имеются трудности в реализации
известных методов для структур со сложным профилем поперечного сечения.
Для преодоления этих трудностей возможно использование метода интегральных уравнений. Основными преимуществами этого метода являются хорошая сходимость, возможность исследовать решетки с произвольной формой поперечного сечения
брусьев, высокая точность, а также возможность расчета различных сечений дисперси2
Lifeng Li. Formulation and comparison of two recursive matrix algorithms for modeling layered diffraction
gratings. // J. Opt. Soc. Am. A. 1996. V. 13, No. 5
3
Knight J.C., Broeng J., Birks T.A. Photonic band gap guidance in optical fibers. // Science. 1998. V. 282,
P. 1476
4
Cregan R.F., Mangan B.J., Knight J.C., Birks T.A. // Science. 1999. V. 285, P. 1537
5
Bullock D.L., Shih C., Margulies R.S. // J. Opt. Soc. Am. 1993. B. V. 10, P. 399.
4
онных поверхностей без каких-либо изменений алгоритма.
Микроструктурные волокна - волокна, имеющие очень сложную структуру в поперечном сечении, образованную регулярным или нерегулярным расположением воздушных отверстий, параллельных оси волокна, обладают множеством замечательных
свойств, таких как возможность регулировки дисперсии и нелинейности, возможность
получения больших значений этих величин и возможность создания одномодового режима для широкого интервала длин волн.
Особое внимание со стороны исследователей уделяется активным микроструктурным волокнам с сердцевиной, допированной ионами неодима, эрбия6 или иттербия.
С развитием волоконных лазеров и волоконных усилителей, появилась необходимость в точных алгоритмах расчета характеристик активных микроструктурных волокон, процессы в которых могут быть описаны с использованием приближения комплексного показателя преломления7 .
Таким образом, задача исследования периодических структур, таких как дифракционные решетки и фотонные кристаллы со сложной формой элементов их образующих,
а также микроструктурные волокна, является актуальной на данный момент, наряду с
задачей разработки эффективных численных методов, использующихся для моделирования процесса распространения электромагнитных волн в подобного рода структурах.
Целью
настоящей работы являются
- исследование процесса дифракции электромагнитной волны на дифракционной
решетке, состоящей из диэлектрических брусьев и разработка алгоритма расчета коэффициентов пропускания и отражения по амплитуде для решетки на основе
метода интегральных уравнений;
- исследование процесса дифракции электромагнитной волны на дифракционной
решетке, состоящей из диэлектрических брусьев, расположенных на подожке, и
разработка алгоритма расчета коэффициентов пропускания и отражения по амплитуде для решетки, находящейся на подложке, на основе метода интегральных
уравнений;
6
Cucinotta A., Poli F., Selleri S. Design of Erbiumdoped triangular photonic-crystal-fiber-based amplifiers.
// IEEE Photon. Tech. Lett. 2004. V. 16, P. 2027-2029.
7
Sader J.E. Method for analysis of complex refractive-index-profile fibers. // Opt. Lett. 1990. V. 15, P.
105–107.
5
- исследование процесса распространения электромагнитной волны в двумерном
фотонном кристалле и расчет различных сечений дисперсионных поверхностей
кристалла при помощи алгоритма, разработанного на основе метода интегральных
уравнений;
- исследование процесса распространения электромагнитной волны в активных
микроструктурных волокнах различных конфигураций и расчет их дисперсионных
характеристик, а также изучение влияния структуры волокна на его дисперсионные характеристики;
- исследование возможности выборочного усиления основной моды и мод высших
порядков, изучение влияния структуры волокна на коэффициент усиления мод с
различными порядковыми номерами.
Новые научные результаты. В диссертации предложены новые выражения для
функции Грина, использующиеся в разработанных алгоритмах решения задач на поиск
коэффициентов пропускания и отражения по амплитуде для решеток, состоящих из диэлектрических брусьев, находящихся в свободном пространстве или расположенных на
подложке.
Необходимость получения новых выражений функции Грина обусловлена наличием в
них логарифмических особенностей в неявном виде, что приводит к снижению точности
расчетов и ухудшению сходимости алгоритмов.
В результате анализа результатов расчетов, выполненных по предложенным алгоритмам, основанных на новых выражениях функции Грина, было показано, что форма
поперечного сечения брусьев решетки оказывает существенное влияние на спектральные и угловые зависимости модуля коэффициента пропускания по амплитуде. Также
было показано, что на эти зависимости оказывает влияние, в том числе, наличие подложки и величина ее диэлектрической проницаемости.
Представленный алгоритм позволяет рассчитать коэффициенты пропускания и отражения по амплитуде для решеток, состоящих из брусьев с любой формой поперечного
сечения.
Новое выражение функции Грина получено также для задачи на поиск собственных волн фотонного кристалла и нахождения его дисперсионных характеристик. Отличительными особенностями алгоритма, основанного на новом выражении для функции
6
Грина, также являются универсальность, то есть возможность его использования для
различных форм поперечного сечения брусьев и хорошая сходимость.
Особенностью, предложенного алгоритма, является также возможность его использования для получения результатов необходимых для построения различных видов дисперсионных характеристик фотонного кристалла: зависимости частоты волны от направления ее распространения и зависимости - и -составляющих при заданном значении частоты распространяющейся волны.
В диссертации также представлен детальный анализ процесса распространения
электромагнитной волны в активных микроструктурных волокнах различных конфигураций, основанный на результатах, полученных методом плоских волн.
Анализ показал, что структура активного волокна, оказывает существенное влияние
на его дисперсионные характеристики и на зависимость коэффициента усиления от длины волны.
Появление в микроструктурном волокне поглощения или усиления приводит к изменению значения эффективного показателя преломления, при этом коэффициент усиления моды зависит от структуры волокна и порядкового номера моды. На основании
расчетов найдены оптимальные значения размеров активной сердцевины для волокон с
заданными конфигурациями, позволяющие наиболее эффективно усиливать основную
моду по сравнению с модами высших порядков. Предложены конфигурации активных
волокон, позволяющие выборочно усиливать моды высших порядков.
Достоверность результатов, полученных в настоящей работе, подтверждается
сравнением с результатами, полученными другими методами, и совпадением результатов, полученных применением разработанных методов в частных случаях, с результатами, полученными методами, оптимизированными для более простых структур.
Научно-практическое значение работы. Предложенные в работе новые алгоритмы расчета коэффициентов пропускания и отражения по амплитуде могут быть использованы для исследования решеток, состоящих из диэлектрических брусьев с произвольной формой поперечного сечения, находящихся в пространстве или расположенных
на подложке.
Результаты этих исследований могут быть использованы для оценки формы поперечного сечения брусьев решетки, определения наличия подложки и определения зна-
7
чения ее диэлектрической проницаемости по амплитуде отраженной или прошедшей
сквозь решетку волны. Подобная методика может быть применима, в частности, для
расчетов концентрации веществ в среде, окружающей дифракционную решетку.
Алгоритм расчета дисперсионных характеристик фотонных кристаллов, основанный
на методе интегральных уравнений, может быть использован для нахождения запрещенных зон и построения изочастот. Результаты подобных расчетов важны при проектировании различных оптических устройств и приборов, таких как фильтры, резонаторы,
волноводы и т.п.
Результаты исследования активных микроструктурных волокон, представленного в
диссертации, расширяют общее представление о влиянии структуры активного волокна
на процесс распространения электромагнитной волны по нему. Предложенные в работе
конфигурации волокон могут быть использованы для выборочного усиления основной
моды и мод высших порядков. Выборочное усиление мод может найти свое применение
при проектировании частотных преобразователей, волоконных лазеров и волоконных
усилителей.
Исследования, проведенные в работе, вносят вклад в развитие метода интегральных
уравнений и представление о возможностях его использования для решения задач лазерной физики и оптики. Предложенный алгоритм может быть использован, в том числе,
и для решения ряда других аналогичных задач после внесения некоторых незначительных изменений.
Апробация работы. Результаты, полученные в ходе исследований, представленных в данной работе, были доложены на следующих конференциях:
1. Saratov Fall Meeting 2004, Россия, Саратов, 21-24 сентября, 2004.
2. Saratov Fall Meeting 2005, Россия, Саратов, 27-30 сентября, 2005.
3. Saratov Fall Meeting 2006, Россия, Саратов, 26-29 сентября, 2006.
4. Третья Международная конференция "Стеклопрогресс - XXI", Россия, Саратов,
22-25 мая, 2006.
5. III Russian-Finnish Meeting Photonics and Laser Symposium, Russia, Moscow,
June 14-17, 2007.
8
6. X Международные Чтения по квантовой оптике, Россия, Самара, 18-22 сентября,
2007.
7. Saratov Fall Meeting 2007, Россия, Саратов, 25-28 сентября, 2007.
8. Российский семинар по волоконным лазерам, Россия, Саратов, 1-4 апреля, 2008.
9. Всероссийская конференция по волоконной оптике, Россия, Пермь, 10-12 октября, 2007.
10. International Conference "Laser Optics 2008", Russia, St.Petersburg, June 23-28,
2008.
11. III Научная конференция для молодых ученых, Россия, Саратов, 25-27 июня,
2008.
Личный вклад автора состоит: в участии в постановках задач, в проведении аналитических выкладок, разработке алгоритмов решения задач и их реализации, в получении
численных результатов и их анализе, в поиске и анализе литературных источников.
Защищаемые положения и результаты.
1.
Новое выражение для ядра интегрального уравнения, отличающееся явно выде-
ленной логарифмической особенностью, для задач на нахождение коэффициентов пропускания и отражения по амплитуде дифракционных решеток, состоящих из диэлектрических брусьев, находящихся в свободном пространстве или расположенных на подложке. Использование этого выражения и соответствующей реализации алгоритма решения интегрального уравнения позволяет рассчитывать коэффициенты пропускания и
отражения решетки при любой форме поперечного сечения брусьев с хорошей сходимостью и точностью вычислений. Угловые и спектральные коэффициенты пропускания
существенно зависят от формы поперечного сечения брусьев и значения диэлектрической проницаемости подложки.
2.
Новое выражение для ядра интегрального уравнения, отличающееся явно вы-
деленной логарифмической особенностью, для расчета дисперсионных характеристик
9
двумерного фотонного кристалла. Алгоритм, основанный на использовании нового выражения для ядра, обладает хорошей сходимостью и позволяет рассчитывать различные
сечения дисперсионных поверхностей без каких-либо изменений алгоритма.
3.
Для активных микроструктурных волокон с заданной структурой существует оп-
тимальный размер активной сердцевины, позволяющий получить максимальную разницу между эффективными усилениями основной моды и мод высших порядков. Существует возможность создания активного микроструктурного волокна, выборочно усиливающего моды высших порядков.
Структура и объем работы. Диссертация состоит из введения, трех глав и заключения, содержит 118 страниц текста, включая иллюстрации. Список литературы на
10 страницах включает 94 наименования.
Краткое содержание работы
Во введении обосновывается актуальность, отмечается новизна и раскрывается
научно-практическая ценность работы, формулируются цели исследований, представляются защищаемые положения и результаты, а также краткое содержание диссертации.
Первая глава посвящена описанию метода интегральных уравнений в применении
к задачам на поиск коэффициентов пропускания или отражения по амплитуде для дифракционных решеток, состоящих из диэлектрических брусьев с произвольной формой
поперечного сечения брусьев, находящихся в свободном пространстве или расположенных на подложке см. Рис. 1.
В Главе предлагается методика выделения логарифмической особенности ядра интегрального уравнения в явной форме, что позволяет увеличить точность и улучшить
сходимость алгоритма. Она заключается в нахождении аналитического выражения для
суммы ряда, к членам которого стремятся члены суммы ядра интегрального уравнения с
ростом их порядкового номера. Аналитическое выражение для суммы ряда также имеет
10
y
θ
ε0
h
r
ε2
εh
x
0
εs
L
Рис. 1: Геометрия структуры.
логарифмическую особенность, но уже выделенную в явной форме, что позволяет избежать интегрирования в особенной точке путем перехода от интеграла по площади к
интегралу по контуру.
Также в Главе представляются спектральные и угловые зависимости модуля коэффициента пропускания дифракционных решеток, состоящих из диэлектрических
брусьев с формой поперечного сечения в виде правильного многоугольника, для случая
отсутствия подожки и для случая наличия подложки.
Демонстрируется влияние формы поперечного сечения брусьев и значения диэлектрической проницаемости подложки на зависимости модуля коэффициента пропускания
от угла падения и длины волны.
Анализируется сходимость алгоритма см. Рис. 2 при увеличении количества гармоник для задачи дифракции на решетке с квадратной формой поперечного сечения
брусьев, диэлектрической проницаемостью подложки рами: ¼
, ¾
, ¾
, и следующими парамет-
см. Рис. 1 и при увеличении количе-
ства азимутальных узлов для различных форм поперечного сечения брусьев, проводится
сравнение полученных результатов с результатами, полученными другими исследователями8 с использованием других методов, а также полученными на основе более простых
приближений см. Таблицу 1.
Вторая глава посвящена описанию метода интегральных уравнений в применении
к задаче на поиск дисперсионных характеристик двумерного фотонного кристалла, образованного круглыми воздушными отверстиями в диэлектрическом материале, упорядоченных так, что элементарная ячейка является квадратом.
8
Ф.Г. Богданов, Г.Ш. Кеванишвили. Дифракция волн на решетках и волноводных неоднородностях.
"Самшобло", Тбилиси, 1994.
11
Таблица 1: Сравнение результатов Ф.Г. Богданова (I) c результатами, полученными по
алгоритму для случая отсутствия подложки (II), и результатами, полученными по алго-
(III), отношение периода
ритму для случая решетки, находящейся на подложке с к длине волны
.
¾
2.1
5.4
¾
I
II
III
0
0.25 0.996579 0.9965788 0.996579
0
0.5
0
0.75 0.990497 0.9904977 0.990506
60
0.25 0.981983 0.9819877 0.981989
60
0.5
60
0.75 0.320571 0.3213180 0.337977
0
0.25 0.928648 0.9286475 0.928681
0
0.5
0
0.75 0.813697 0.8135704 0.816695
0.975028 0.9750275 0.975065
0.785866 0.7859282 0.786197
0.884717 0.8847234 0.884885
Также как и в первой главе, ядро интегрального уравнения было выражено в новом
виде, что в конечном итоге дало возможность избежать вычисления значения функции
Грина в особенной точке. Алгоритм, основанный на новом выражении для ядра интегрального уравнения, хорошо сходится см. Рис. 2 и позволяет получить высокую точность результатов.
По разработанному алгоритму была рассчитана зависимость квадрата волнового
числа от составляющей вектора распространения волны внутри фотонного кристалла
при заданном значении
составляющей. Расчеты проводились для следующих пара-
метров: диэлектрическая проницаемость вмещающей среды ¼
мкм, период структуры по обоим направления , радиус отверстий
мкм.
Достоверность результатов подтверждена их сравнением с аналогичными результатами, полученными при использовании метода плоских волн см. Рис. 3.
Кроме того, представленный метод используется в работе для построения изочастоты - зависимости пространственных составляющих вектора распространения волны внутри кристалла для заданной длины волны. Тем самым иллюстрируется гибкость
метода интегральных уравнений, позволяющего решать задачи разных типов.
12
Рис. 2: Кривые сходимости. Зависимость модуля амплитуды нулевой гармоники прошедшей волны от числа гармоник, использующихся при разложении ядра интегрального уравнения, для задачи дифракции электромагнитной волны на решетке с квадратной
формой поперечного сечения брусьев (слева) , зависимость квадрата волнового числа ¾
от числа гармоник, использующихся при разложении ядра интегрального уравнения, для
задачи на поиск дисперсионных характеристик двумерного фотонного кристалла (справа).
Третья глава
посвящена исследованию активных микроструктурных волокон и
влиянию их структуры см. Рис. 4 на зависимости эффективного показателя преломления
и эффективного усиления от частоты волны.
Задача на поиск собственных значений решается методом плоских волн, при этом
для формирования симметричной суперячейки см. Рис. 4 используется косоугольная система координат.
В Главе исследуются различные конфигурации волокон, для которых рассчитываются и строятся дисперсионные характеристики. Для получения численных результатов
используется метод плоских волн и приближение комплексного показателя преломления.
Иллюстрируется влияние наличия усиления или поглощения в волокне на вещественную и мнимую части эффективного комплексного показателя преломления.
Строятся зависимости мнимой части эффективного показателя преломления от порядкового номера моды (моды нумеруются в порядке убывания значения вещественной части эффективного комплексного показателя преломления) для различных радиусов активной сердцевины. Показывается, что для микроструктурных волокон с заданной конфигурацией, существуют оптимальные для выборочного усиления основной мо-
13
Рис. 3: Зависимость ¾ от при , рассчитанная методом интегральных уравнений
(слева) и методом плоских волн (справа).
Рис. 4: Структура исследуемого волокна (слева) и суперячейка, использующаяся при
анализе структуры (справа).
Рис. 5: Зависимость мнимой части эффективного комплексного показателя преломления (усиления) от номера моды для микроструктурных волокон с активными сердцевинами различных диаметров (слева), зависимость разницы усиления между 2-й и 4-й
модой от радиуса активной сердцевины (справа).
ды размеры активной сердцевины см. Рис. 5.
14
Предлагаются конфигурации активных микроструктурных волокон, позволяющие
существенно увеличить эффективный коэффициент усиления мод высших порядков и
выделить их тем самым по отношению к модам с другими порядковыми номерами.
В заключении подводится итог проведенных в работе исследований, формулируются основные заключения и выводы.
Основные результаты и выводы
1.
Предложен алгоритм решения задачи дифракции электромагнитной волны на
диэлектрической решетке из параллельных цилиндрических брусьев, находящихся в
свободном пространстве или расположенных на подложке.
Аналитические преобразования, представленные в работе, позволяют получить
функции Грина для этих задач в новом виде, в котором особенность выделена в явной
форме, что дает возможность заменить интегралы непосредственно на квадратурные
суммы.
Представленные алгоритмы решения двух вышеописанных типов задач, позволяют
рассчитать коэффициенты пропускания и отражения по амплитуде для дифракционных
решеток с высокой точностью, в отличие от алгоритмов, использующих функцию Грина
свободного пространства.
Этот факт подтверждается совпадением результатов, полученных предложенным
методом с результатами, полученными другими методами, и результатами, полученными методами, оптимизированными для более простых структур.
Были проведены оценки, показывающие высокую точность и хорошую сходимость
алгоритма, в том числе и для случаев сложной формы поперечного сечения брусьев.
2.
Построены угловые и спектральные зависимости коэффициента пропускания
для различных форм поперечного сечения диэлектрических брусьев, образующих решетку. Было продемонстрировано существенное влияние формы поперечного сечения
брусьев на кривые зависимостей коэффициентов пропускания решеток от угла падения и длины падающей волны. Этот факт дает возможность использования подобного
алгоритма для теоретических расчетов и оценок форм поперечного сечения микро- и
15
наноразмерных брусьев по спектрам отражения или пропускания исследуемых дифракционных решеток.
3.
При решении задачи на нахождение коэффициента пропускания дифракцион-
ной решетки, расположенной на подложке, были построены угловые и спектральные
зависимости коэффициента пропускания для случая отсутствия подложки и для случая
присутствия подложки с диэлектрической проницаемостью . Было продемон-
стрировано существенное влияние наличия подложки на кривые угловых и спектральных зависимостей, что позволяет использовать предложенный алгоритм для теоретического исследования спектров пропускания или отражения решетки, расположенной
на подложке, с целью определения наличия подложки или значения ее диэлектрической
проницаемости. Подобный метод, в частности, может быть использован для определения концентрации веществ в среде, окружающей дифракционную решетку.
4.
Предложен алгоритм нахождения дисперсионных характеристик двумерного
фотонного кристалла, основанный на методе интегральных уравнений.
Аналитические преобразования, описанные в работе, позволяют представить функцию Грина в новой форме с явно выделенной логарифмической особенностью, что позволяет найти точное значение интеграла от функции, ее содержащей.
Достоверность результатов, подтверждающаяся хорошим согласием значений, рассчитанных двумя разными методами: методом интегральных уравнений и методом плоских волн, доказывает правильность формулировки предлагаемого алгоритма и нового
выражения для функции Грина.
5.
Использование высокоточных квадратурных формул позволяет свести инте-
гральное уравнение к системе однородных линейных алгебраических уравнений. Условие равенства определителя этой системы нулю представляет собой дисперсионное
уравнение, которое позволяет найти различные дисперсионные характеристики фотонного кристалла, такие как зависимость частоты волны от направления ее распространения и зависимость - и -составляющих при заданном значении частоты распространяющейся волны без каких-либо изменений алгоритма.
Этот факт показывает гибкость метода интегральных уравнений, позволяющего рассчитывать различные дисперсионные характеристики без каких-либо изменений программной реализации алгоритма.
16
6.
Представлены результаты исследования различных конфигураций активных
микроструктурных волокон.
Для исследования усиления в микроструктурных волокнах использовалось приближение комплексного показателя преломления. В таком приближении коэффициент усиления моды пропорционален мнимой части эффективного показателя преломления соответствующей моды.
7.
Получены дисперсионные характеристики для волокна с воздушными отверсти-
ями целиком состоящего из активного стекла, то есть с однородным распределением
усиления. Показано, что зависимость вещественной части эффективного показателя
преломления активного волокна отличается от аналогичной зависимость для волокна,
изготовленного из обыкновенного стекла. Также показано, что мнимая части эффективного показателя усиления активного волокна зависит от длины распространяющейся
волны.
8.
Получены дисперсионные характеристики для волокна с активной сердцевиной.
Показано, что зависимость вещественной части эффективного показателя преломления
волокна с активной сердцевиной отличается от аналогичной зависимости для волокна
без активной сердцевины. Также показано, что мнимая часть эффективного показателя
усиления волокна с активной сердцевиной зависит от длины распространяющейся волны.
9.
Получены зависимости коэффициентов усиления мод от их порядкового номера
для различных размеров активных сердцевин в волокне с переменным радиусом воздушных включений и волокне с большим радиусом воздушных включений. Показано,
что существуют оптимальные значения радиусов активных сердцевин для этих двух типов волокон, при которых разница между усилением основной моды и усилением мод
высших порядков максимальна.
10. Были получены зависимости коэффициентов усиления мод от их порядкового
номера для активных волокон с несколькими активными сердцевинами, которые показывают, что определенное пространственное расположение активных сердцевин может
привести к выделению мод высших порядков, за счет наличия у них максимального эф-
17
фективного показателя распределения, превосходящего в том числе и показатель усиления основной моды.
Основные публикации по теме работы
1. A.S. Soloviev, I.S. Nefedov. Analysis of the wave diffraction on a polygon dielectric
grating by integral equation method. // Proceedings of SPIE 2004. 2005. V. 5773. P. 59-64.
2. А.С. Соловьев. Расчет дисперсии в стеклянном микроструктурном волокне с
поглощением или усилением. // Сборник докладов 3-й Международной конференции
"Стеклопрогресс - XXI". 2007. C. 228-233.
3. I.S. Nefedov, A.S. Soloviev. Analysis of the wave diffraction on a polygon dielectric
grating, placed on a dielectric substrate, by integral equation method. // Proceedings of
SPIE 2006. 2007. V. 6537.
4. A.S. Soloviev, A.I. Konyukhov, L.A. Melnikov, S.A. Akishin. Mode gain in
microstructure optical fibers. // Proceedings of SPIE 2006. 2007. V. 6537.
5. А.И. Конюхов, А.С. Соловьев, Л.А. Мельников, С.А. Акишин. Усиление направляемых мод в микроструктурных оптических волокнах. // Известия СГУ. 2007. Т. 7. №.
2. С. 30-36.
6. И.С. Нефедов, А.С. Соловьев. Дифракция на решетке диэлектрических брусьев с
сечением правильного многоугольника, расположенной на подложке. // Оптика и спектроскопия. 2008. Т. 104. № 3. С. 486-493.
7. А.С. Соловьев, И.С. Нефедов. Метод интегральных уравнений в задаче дифракции
волны на решетке из параллельных диэлектрических брусьев с сечением правильного
многоугольника. // Оптический журнал. 2008. Т. 75. № 5. С. 12-17.
8. А.И. Конюхов, А.С. Соловьев, Л.А. Мельников, С.А. Акишин. Усиление направляемых мод в микроструктурированных оптических волокнах. // II Российский семинар
по волоконным лазерам 2008. Материалы семинара. 2008. С. 18-19.
9. А.С. Соловьев, Л.А. Мельников, А.И. Конюхов. Усиление направляемых мод в
микроструктурных оптических волокнах. // Фотон-Экспресс. 2007. № 6 (62). С. 95-96.
18
СОЛОВЬЕВ Андрей Сергеевич
ПРИМЕНЕНИЕ МЕТОДОВ ИНТЕГРАЛЬНЫХ УРАВНЕНИЙ И
ПЛОСКИХ ВОЛН ДЛЯ РАСЧЕТА ДИФРАКЦИИ НА
ДИЭЛЕКТРИЧЕСКИХ СТЕРЖНЯХ И ПОИСКА СОБСТВЕННЫХ
ВОЛН УСИЛИВАЮЩИХ МИКРОСТРУКТУРНЫХ ВОЛОКОН
Автореферат
Подписано к печати 25.09.2008. Формат 60 × 84 1/16.
Бумага офсетная. Гарнитура «Times».
Усл. печ. л. 1,16 (1,25). Тираж 100 экз. Заказ 185-T.
Типография Саратовского государственного университета
410012, Саратов, ул. Большая Казачья, 112а, корпус 8.