CC BY
14
_RUDN Journal of MlPh_ 2018 Vol. 26 No. 4 321-330
Вестник РУДН. Серия МИФ http://journals.rudn.ru/miph
Математическое моделирование
УДК 519.633.2
DOI: 10.22363/2312-9735-2018-26-4-321-330
О нормальных модах закрытого волновода с разрывным
заполнением М. Д. Малых
Кафедра прикладной информатики и теории вероятностей Российский университет дружбы народов ул. Миклухо-Маклая, д. 6, Москва, Россия, 117198
Рассматривается волновод постоянного поперечного сечения S с идеальным проведением стенками. Предполагается, что заполнение волновода не изменяется вдоль его оси и описывается кусочными непрерывными функциями е и /л на поперечном сечении волновода. Показано, что возможно сделать замену переменных, которая позволяет работать только с непрерывными функциями.
Вместо разрывных поперечных компонент электромагнитного поля Е и Н мы предлагаем использовать четыре потенциала ие,ин и ve,vh. Мы можем доказать как обобщение теоремы Тихонова-Самарского, что любое поле в волноводе допускает представление в такой форме,
o 1
если мы рассматриваем потенциалы ие,иь как элементы пространства Соболева W¿(S), а потенциалы ve,vh, как элементы пространства Соболева W}(S).
Если е и /л — кусочные постоянные функции, то уравнения Максвелла, записанные в четырёх потенциалах, сводятся к двум независимым системам. Это обстоятельство даёт нам новый подход к исследованию спектральных свойств волноводов. Во-первых, мы можем доказать полноту системы нормальных волн в закрытых волноводах, используя стандартные функциональные пространства. Во-вторых, мы можем предложить новую технику для вычисления нормальных волн, используя стандартные конечные элементы. В конце статьи представлена программа, написанная на языке FreeFem+—+, для вычисления дисперсионных линий волновода. Также рассмотрен вопрос о вычислении мод при больших значениях к = ш/с.
Ключевые слова: волновод, уравнения Максвелла, пространства Соболева, метод конечных элементов, нормальные моды
1. Введение
Закрытые волноводы являются наиболее изученным объектом математической теории волноведущих систем. Ещё в середине прошлого века были предложены корректные модели распространения волн по таким системам, в частности была рассмотрена весьма сложная как для теоретического анализа, так и для применения численных методов векторная модель, в которой рассматриваются поля, удовлетворяющие уравнениям Максвелла. Наиболее изученными оказались волноводы, заполненные однородным веществом, исчерпывающая теория которых была построена в ставших уже классическими работах А.Н. Тихонова и А. А. Самарского [1].
Тем не менее математическая сторона моделирования волноводов с сложным заполнением на основе уравнений Максвелла не развита в полной мере. Дело в том, что ключевым техническим моментом в теории полого волновода является использование потенциалов Герца или функции Боргниса [2, 3]. Этот приём не удалось обобщить на случай волноводов с переменным заполнением, и это заставляет работать с компонентами электромагнитных полей в специально конструируемых пространствах Соболева, для которых приходится в теории доказывать свои теоремы вложения [4-7], а на практике — использовать смешанные конечные элементы [8—
Статья поступила в редакцию 14 сентября 2018 г.
10]. Это, в частности, существенно осложняет использование свободных и хорошо отлаженных конечно элементных программ (FEA software) для расчёта волноводов со сложным заполнением.
Современное развитие волоконной оптики и создание мета-материалов вновь возвращает нас к модели закрытого волновода, заполнение которого описывается кусочно-постоянной функцией. Вероятно, наиболее очевидным примером являются многожильные оптические волокна (multicore fiber), которые активно исследуются в последнее время, поскольку пропускная способность обычных одноканальных волокон использована уже почти полностью [11] и, очевидно, недостаточна для перемещения огромных массивов данных, что является ключевым моментом на пути к мобильным сетям 5G. Следует также заметить, что распространение волноводных мод по открытым оптическим волноводам, равно как и взаимодействие нескольких таких волноводов, можно корректно моделировать, поместив всю волноведующую систему в ящик, на границе которого поставлены условия идеальной проводимости [12]. Не трудно заметить, что система «волокно + ящик» представляет собой ни что иное, как закрытый волновод со сложным заполнением.
Конечной целью наших исследований [13-15] будет разработка численных методов исследования волноводов в стандартных функциональных пространствах и их реализация на высокоуровневых языках программирования, таких как FreeFemH—H [16,17], или в системах компьютерной алгебры. При этом мы стремимся использовать как можно более простые, а следовательно, и надёжные конструкции. На наш взгляд, ключевая проблема в исследовании волноводов с разрывным заполнением — необходимость работать с разрывными компонентами полей. Обычно потенциалы рассматривают как способ понижения системы уравнений в частных производных, то есть как способ интегрирования системы. Мы же хотим взглянуть на введение потенциалов как на замену переменных, при которой мы переходим от разрывных функций к гладким потенциалам. В этой статье мы опишем такую замену и покажем, как при её помощи вычислить дисперсионные кривые во FreeFemH—H.
2. Декомпозиция Гельмгольца
Рассмотрим закрытый волновод постоянного сечения 5 с кусочно непрерывным заполнением е и не меняющимся вдоль оси волновода. Линию разрыва заполнения обозначим как Г. Направим ось Ох декартовой системы координат по оси волновода.
Под электромагнитным полем в закрытом волноводе 5 х 2 х Т с заполнением е,^ будем понимать векторные поля Е,Н, компоненты которых определены на — Г) х 2 х Т, при условии, что сужение Е, Н и их частных производных по г и £ на сечение 5 при любых значениях г и £ являются кусочно гладкими функциями, удовлетворяющими
1. уравнениям Максвелла
rotE = dt Н, с
rotH = + -dtE, с
V • еЕ = 0,
V • ¡JÉ = 0
(1)
внутри волновода 5 х 2 х Т; 2. условиям идеальной проводимости стенок волновода
Ё х п = 0, Н • п = 0 в регулярных точках границы дБ х 2 х Т;
(3)
3. условиями сопряжения
f[E х п] = 0, [еЁ ■ п] = О \[Й х п] = 0, [рН ■ п] = 0
в регулярных точках границы разрыва заполнения Г х Z х Т. Наконец, примем для краткости, что
Л± = (АХ,АУ, Of and V = (дх,ду, 0)т, V' = (-ду,дх, 0)т. Связь между полями и потенциалами зададим следующим образом:
Е± = Vue + 1 V've, Н± = Vvh + 1 V'uh. (4)
e ^
Каждая из этих формул представляет собой двумерный аналог декомпозиции Гельмгольца, хорошо известной в теории упругости [18].
Замечание 1. В электродинамике для поля Н± такие потенциалы возникали при доказательстве полноты системы нормальных мод в качестве вспомогательной конструкции [5]. Все четыре потенциала были введены в нашей работе [14] для гладкого заполнения без коэффициентов 1/е и 1/^, важных только для разрывного случая.
Теорема 1. Для любого электромагнитного поля Е,Н в волноводе найдутся
о
такие функции ие, ии переменных х, £ со значениями в пространстве Соболева Ш2 (5) и такие функции уе, переменных х, £ со значениями в пространстве Соболева ), что справедливо равенство (4). Указанное представление единственно с точностью до аддитивных констант.
Теорема 1 означает, что при переходе от переменных Е, Н к четырём потенциалам и двум компонентам Ez, Hz по формулам (4) не теряются решения уравнений Максвелла. При этом условия
o
Ue,Uh,Ez e W\(S) and Ve,Vh,Hz e W^(S)
заменяют нам условия на разрывах заполнения, равно как и граничные условия. Поскольку потенциалы являются элементами пространств Соболева, далее естественно рассматривать уравнения Максвелла в слабой форме [19].
Обратимся теперь к случаю, когда заполнение волновода является кусочно постоянным. Уравнения Максвелла дают, что потенциалы ие, Uh и Ez — элементы
o
W^(S), связанные уравнениями
e(Vu, Vue)dxdy = dz J J euEz dxdy,
s s
с
H
s
— (Vu, Vuh)dxdy = -dt JJ еиЕг
s
dxdy
для любой и из Со°(Б), где Ег = дгие + дги}, потенциалы уе, V} и Нг — элементы ^¡¡(^), связанные уравнениями
-(Уу, Ууе)Лхйу = дг
(6)
^(Уь, Уv}l)dxdy = дг / йхйу,
для любой V из С^(5), где Нг = дг— дгье.
Уравнения (5) и (6) можно использовать и для конструирования полей в волново-
о
де. Если ие, и} и Ех из Ш¡(5) удовлетворяют уравнениям (5), а уе, V} и Нх из ^¡(^О удовлетворяют уравнениям (6), то поле Е,Н, вычисленное по формулам (4), удовлетворяет уравнениям Максвелла в обобщённом смысле. Более того, если это поле имеет вне разрывов заполнения непрерывные частные производные 1-го порядка по всем переменным, а на разрывах заполнения — разрывы 1-го рода, то это поле вне разрывов заполнения соответствует уравнениям Максвелла (1), условиям сопряжения (3) на разрывах заполнения и краевым условиям идеальной проводимости (2).
Поскольку система уравнений Максвелла распалась на две независимые системы,
электромагнитное поле Е,Н в волноводе, заполнение которого описывается кусочно постоянными функциями е и представляет собой суперпозицию ТЕ- и ТМ- полей.
3. Нормальные ТМ-моды
Покажем, как развитую теорию можно применить к отысканию нормальных мод волновода. Напомним, что под нормальной модой волновода будем понимать нетривиальное поле вида
Е = Е(х, Н = Н(х,
где ш — положительное число (круговая частота моды), а 7, вообще говоря, число комплексное (волновое число). Моды с вещественным волновым числом называют нормальными волнами, распространяющимися вдоль или против оси волновода.
Нормальная ТМ-мода волновода описывается потенциалами
ие = йев^"-™*, ин = ,
которые удовлетворяют уравнениям
/ г г г г г г
„2
е(Уи, Уйе)ёхёу = —7 JJ еийес!хс1у + к^ JJ euйhdxdy,
в в
— (Уи, Vйh)dxdy = —kl JJ ^и^^еЛхйу + к2 JJ еий}dxdy,
(7)
5 в
о
для любой и Е Ш2(5); здесь для краткости используется обозначение к = ш/с. Перепишем эту систему уравнений в операторном виде, используя стандартную
технику теории пространств Соболева [20],
(Аейе = -7 2Вейе +
н, , :
2 (8)
А1 йн = -к^Вейе + к2Вейн
где Ае,А 1 ,Ве — ограниченные самосопряжённые операторы, а Ве ко всему прочему ещё и вполне непрерывный.
Для теоретического анализа удобно исключить йн и оставить одно уравнение
1
Аейе = -72 + В^^Л1 - Ве^ В^ йе. (9)
Поэтому задачу об отыскании всех нормальных ТМ-мод при заданной частоте ш можно рассматривать как задачу на собственные значения для операторного пучка, в которой Л = 72 рассматривается как спектральный параметр. В силу самосопряжённости операторов и полной непрерывности оператора Ве собственные
о
функции этой задачи образуют базис пространства Ш2(5). Отсюда нетрудно вывести, что любое ТМ-поле в волноводе можно представить в виде суперпозиции нормальных ТМ-полей.
Замечание 2. Следует подчеркнуть, что мы установили именно базисность системы нормальных мод, полнота системы нормальных мод для волноводов со сложным заполнением была установлена в работах А.Л. Делицына [6,7]. Этот результат является естественным обобщением теоремы о базисности системы нормальных мод полого волновода, установленной ещё в 1940-х годах А. Н. Тихоновым и А. А. Самарским [1,3,21].
При фиксированной частоте ш имеется бесконечное число нормальных мод, из них лишь конечное число представляют собой бегущие волны. Прочие моды имеют чисто мнимые значения 7. Наши вычисления мод волновода, заполненного веществом, основанные на неполном методе Галёркина [15], приводили в высших модах к волновым числам, имеющим как вещественную, так и мнимую части, одна из которых была очень мала. Теперь ясно, что это был чисто численный эффект.
Чтобы описать параметры распространяющихся нормальных волн, удобно использовать дисперсионную кривую. Все точки ¿^-плоскости, где эта задача на собственные значения имеет нетривиальное решение, образуют кривую, которую будем называть дисперсионной кривой волновода. Для её построения естественно использовать метод усечения: мы будем использовать пространство конечных элементов вместо пространства Соболева и изменим операторами разреженными матрицам, порождёнными теми же самыми билинейными формами.
Для вычисления этих матриц и дальнейших манипуляций с блочными разреженными матрицами мы использовали свободное программное обеспечение РгееРеш +—+ [17]. На этом языке, не опускаясь на уровень матриц и чисел, можно задать пространство конечных элементов и матрицы операторов Ае, А1, Ве. Однако вычисление
обратной матрицы приведёт к задаче на собственные значения с неразреженными матрицами, что существенно усложнит её решение. Поэтому для вычислений удобнее пользоваться исходной системой (8). Подставляя 7 = к@, мы можем переписать нашу задачу на собственные значения (8) в блочно-разреженном виде:
(Ае 0 \ /йе\ 2 (-Р2В< №е\ (йе\
\ 0 \йн) = к у-рв£ ве) (й^ .
Эту задачу можно решить стандартными средствами, если рассматривать к2 как собственное значение, а /3 — как параметр. Для построения дисперсионной кривой волновода теперь достаточно решать эту стандартную задачу, меняя значения /3 с некоторым шагом. Это позволяет нам работать во ГгееГетН—Ь со всеми волноводами, границы которых могут быть описаны параметрически с помощью элементарных функций, а заполнение — с помощью алгебраических неравенств.
Вычислительные возможности этой программы проиллюстрируем примером.
Пример 1. На рис. 1 представлена дисперсионная кривая для волновода с сечением
5" = {0 < ж < 1} х {0 < у < 1} и кусочно постоянным заполнением
1,2, (Ж-о, 5)2 +(у-0, З)2 <0,5,
1,
иначе
ц
(10)
7
12
10
б 8
10
12
14
Рис. 1. Дисперсионная кривая для примера 1. Сплошной линией отмечена ветвь 71 (к), пунктиром — ветвь 72(к), которая с графической точностью
совпадает с третьей ветвью
Мы использовали 2120 треугольников для построения пространства конечных элементов и меняли значения /3 с шагом А/3 = 0.1. Чтобы проверить сходимость, мы сделали ряд числовых экспериментов для полого волновода.
4. Высокочастотный предел
В оптических задачах величина к = ^ весьма велика, поэтому весьма полезно взглянуть на высокочастотный предел [22]. Мы не можем перейти к пределу к —)- оо непосредственно в (9), но мы можем сделать это после усечения. Дело в том, что оператор Ве является вполне непрерывным и поэтому необратим. Однако МКЭ вносит некоторую регуляризацию: матрица Ве становится обратимой, причём норму обратной матрицы можно оценить через характерный линейный размер Ь, используемой сетки как
цб^1!! ~ /г2.
Это обстоятельство позволяет применить формулу Неймана к задаче (9) после её дискретизации. В результате получится, что
В А т1 А1 - вЛ В, = -В, {в-1 + 4 В-1 А1 В-1 + о( Л I I В
(f*i- в-)- в- = -В- (V + FВ-1Аiв+° (£)
*+° (£))
йе = ^ (А1 + i
Если обозначить собственные значения пучка Ае - ¡32 А1 как ¡Зп, п = 1, 2,..., то собственные значения задачи (9) асимптотически равны
Тп = Рпк + ....
Таким образом, нормальные моды имеют вид
¿(ж, , н(х, уУ*^"-™*.
Условие применимости формулы Неймана состоят в том, что к~2к~2 ^ 1, т.е. шаг сети И, должен быть меньше длины волны в вакууме. Мы полагаем, что найденная формула может быть полезна для управления величиной отношения 7/к волноводных мод путём изменения заполнения волновода.
и
5. Заключение
В настоящей статье мы хотели показать, что теоретические и численные исследования электромагнитных полей в закрытых волноводах, заполнения которых описываются кусочно-постоянными функциями е и возможно вести в стандартных пространствах Соболева и при помощи обычных конечных элементов. От разрывных компонент полей мы предлагаем перейти к четырём непрерывным скалярным функциям — потенциалам. Этот приём позволяет в теории легко обосновать ба-зисность системы нормальных мод такого волновода, а на практике предложить способ приближенного вычисления нормальных мод, использующий стандартные средства численного анализа, разработанные для скалярных краевых задач математической физики. Для подтверждения этого проделаны численные эксперименты в среде ЕгееЕеш++.
Благодарности
Автор признателен профессору А.Н. Боголюбову (МГУ), профессору Л.А. Севастьянову (РУДН) и доценту К. П. Ловецкому (РУДН) за обсуждение работы и весьма интересные замечания.
Публикация подготовлена при поддержке Программы РУДН «5-100» и при частичной поддержке грантов РФФИ №№ 18-07-00567 и 18-51-18005.
Литература
1. Самарский А. A., Тихонов А. Н. О представлении поля в волноводе в виде суммы полей ТЕ и ТМ // Журнал технической физики. — 1948. — Т. 18, № 7. — С. 959-970.
2. Zhang K., Li D. Electromagnetic Theory for Microwaves and Optoelectronics. — 2 edition. — Berlin: Springer, 2008.
3. Могилевский И. Е., Свешников А. Г. Математические задачи теории дифракции. — Москва: Физический факультет МГУ, 2010.
4. Боголюбов А. Н., Делицын А. Л., Свешников А. Г. О полноте системы собственных и присоединённых функций волновода // Журнал вычислительной математики и математической физики. — 1999. — Т. 38, № 11. — С. 1891-1899.
5. Делицын А. Л. Об одном подходе к вопросу о полноте нормальных волн волновода с магнитодиэлектрическим заполнением // Дифференциальные уравнения. — 2000. — Т. 36, № 5. — С. 629-633.
6. Боголюбов А. Н, Делицын А. Л., Малых М. Д. О корневых векторах цилиндрического волновода // Журнал вычислительной математики и математической физики. — 2001. — Т. 41, № 1. — С. 126-129.
7. Делицын А. Л. О полноте системы собственных векторов электромагнитных волноводов // Журнал вычислительной математики и математической физики. — 2011. — Т. 51, № 10. — С. 1883-1888.
8. Делицын А. Л. О проблеме применения метода конечных элементов к задаче вычисления мод волноводов // Журнал вычислительной математики и математической физики. — 1999. — Т. 39, № 2. — С. 315-322.
9. Делицын А. Л., Круглов С. И. Применение метода смешанных конечных элементов для вычисления мод цилиндрических волноводов с переменным показателем преломления // Журнал радиоэлектроники. — 2012. — № 4.
10. Lezar E., Davidson D. B. Electromagnetic Waveguide Analysis // Automated Solution of Differential Equations by the Finite Element Method. — The FEniCS Project, 2011. — Pp. 629-643.
11. Coffey V. C. Novel Fibers Use Space to Extend Capacity Limits // Photonics Spectra. — 2013. — Vol. 4, No 7.
12. Моделирование распространения поляризованного света в тонкоплёночной волно-водной линзе / Д. В. Диваков, М. Д. Малых, А. Л. Севастьянов и др. // Вестник РУДН. Серия: Математика. Информатика. Физика. — 2017. — Т. 25, № 1. — С. 56-68.
13. On the Representation of Electromagnetic Fields in Closed Waveguides Using Four Scalar Potentials / M. D. Malykh, L. A. Sevastianov, A. A. Tiutiunnik, N. E. Niko-laev // Journal of Electromagnetic Waves and Applications. — 2018. — Vol. 32, No 7. — Pp. 886-898.
14. О сведении уравнений Максвелла в волноводах к системе связанных уравнений Гельмгольца / М. Д. Малых, А. Л. Севастьянов, Л. А. Севастьянов, А. А. Тю-тюнник // Вестник РУДН. Серия: Математика. Информатика. Физика. — 2018. — Т. 26, № 1. — С. 39-48.
15. Тютюнник А. А. О вычислении электромагнитных полей в закрытых волноводах с неоднородным заполнением // Вестник РУДН. Серия: Математика. Информатика. Физика. — 2018. — Т. 26, № 2. — С. 129-139.
16. Hecht F. New Development in FreeFem++ // J. Numer. Math. — 2012. — Vol. 20, No 3-4. — Pp. 251-265. — ISSN 1570-2820.
17. Hecht F., 2018. — Freefem++. — Laboratoire Jacques-Louis Lions, Universite Pierre et Marie Curie, Paris, 3 edition. — www.freefem.org.
18. Ляв Д. Теория упругости. — Москва, Ленинград: ГТТИ, 1939.
19. Дюво Г., Лионс Ж.-Л. Неравенства в механике и физике. — Москва: Наука, 1980.
20. Stummel F. Rand- und Eigenwertaufgaben in Sobolewschen Raumen. — Berlin-Heidelberg-New York: Springer, 1969.
21. Chew W. C. Lectures on Theory of Microwave and Optical Waveguides. — 2012. — http://wcchew.ece.illinois.edu/chew/course/tgwAll20121211.pdf.
22. Бабич В. М., Булдырев В. С. Асимптотические методы в задачах дифракции коротких волн. — Москва: Наука, 1972.
UDC 519.633.2
DOI: 10.22363/2312-9735-2018-26-4-321-330
On Normal Modes of the Closed Waveguide with Discontinuous
Filling Mikhail D. Malykh
Department of Applied Probability and Informatics Peoples' Friendship University of Russia (RUDN University) 6, Miklukho-Maklaya str., Moscow, 117198, Russian Federation
We consider a waveguide of a constant cross-section S with ideally conducting walls. We assume that the filling of waveguide doesn't change along its axis and is described by the piecewise continuous functions £ and ^ defined on waveguide cross-section. We show that it is possible to make substitution which allows to work only with continuous functions.
Instead of noncontinuous cross-components of an electromagnetic field E and H we offer to use four potentials ue,uh and ve,vh. We can prove as the generalization of Tikhonov-Samarskii theorem that any field in the waveguide allows representation in such form if we consider the
o 1
potentials ue,uh as elements of Sobolev space W2(S) and the potentials ve,vu as elements of Sobolev space W2l(S).
If e and ^ are the piecewise constant functions then Maxwell's equations written in four potentials reduce to a pair of independent systems. This statement give us new approach to the investigation of spectral properties of waveguides. First, we can prove the completeness of the system of the normal waves in closed waveguides using standard functional spaces. Secondly, we can offer new technique for calculation of the normal waves using standard finite elements. FreeFem++ program for calculation of disperse lines of waveguides is presented. The question of calculation of modes at great values of k = is also considered.
Key words and phrases: waveguide, Maxwell's equations, Sobolev's spaces, finite element method, normal modes
References
1. A. A. Samarskiy, A. N. Tikhonov, On the Representation of a Field in a Waveguide in the Form of a Sum of Fields TE and TM, Technical Physics. The Russian Journal of Applied Physics [Zhurnal tekhnicheskoy fiziki] 18 (7) (1948) 959-970, in Russian.
2. K. Zhang, D. Li, Electromagnetic Theory for Microwaves and Optoelectronics. 2nd ed., Springer, Berlin, 2008.
3. I. E. Mogilevskii, A. G. Sveshnikov, Mathematical Problems of the Theory of Diffraction, Faculty of Physics MSU, Moscow, 2010, in Russian.
4. A. N. Bogolyubov, A. L. Delicyn, A. G. Sveshnikov, On the problem of the Excitation of a Waveguide with an Inhomogeneous Medium, Computational Mathematics and Mathematical Physics 38 (11) (1999) 1815-1823.
5. A. L. Delicyn, On One Approach to the Question of the Completeness of Normal Waves of a Waveguide with a Magnetodielectric Filling, Differentsialnye Uravneniya 36 (5) (2000) 629-633, in Russian.
6. A. N. Bogolyubov, A. L. Delicyn, M. D. Malykh, On the Root Vectors of a Cylindrical Waveguide, Computational Mathematics and Mathematical Physics 41 (1) (2001) 121-124, in Russian.
7. A. L. Delicyn, On the Completeness of the System of Eigenvectors of Electromagnetic Waveguides, Computational Mathematics and Mathematical Physics 51 (10) (2011) 1771-1776.
8. A. L. Delicyn, Application of the Finite Element Method to the Calculation of Modes of Dielectric Waveguides, Computational Mathematics and Mathematical Physics 39 (2) (1999) 298-304, in Russian.
9. A. L. Delicyn, S. I. Kruglov, Application of a Method of the Mixed Finite Elements for Calculation of Modes of Cylindrical Waveguides with Variable Index of Refraction, Journal of radio electronics (4), in Russian.
10. E. Lezar, D. B. Davidson, Electromagnetic Waveguide Analysis, in: Automated solution of differential equations by the finite element method, The FEniCS Project, 2011, pp. 629-643, in Russian.
11. V. C. Coffey, Novel Fibers Use Space to Extend Capacity Limits, Photonics Spectra 4 (7), in Russian.
12. D. V. Divakov, M. D. Malykh, A. L. Sevastianov, L. A. Sevastianov, Simulation of Polarized Light Propagation in the Thin-Film Waveguide Lens, RUDN Journal of Mathematics, Information Sciences and Physics 25 (1) (2017) 56-68, in Russian.
13. M. D. Malykh, L. A. Sevastianov, A. A. Tiutiunnik, N. E. Nikolaev, On the Representation of Electromagnetic Fields in Closed Waveguides Using Four Scalar Potentials, Journal of Electromagnetic Waves and Applications 32 (7) (2018) 886-898.
14. M. D. Malykh, A. L. Sevastianov, L. A. Sevastianov, A. A. Tyutyunnik, On the Reduction of Maxwell's Equations in Waveguidesto the System of Coupled Helmholtz Equations, RUDN Journal of Mathematics, Information Sciences and Physics 26 (1) (2018) 39-48, in Russian.
15. A. A. Tyutyunnik, On the Calculation of Electromagnetic Fields in Closed Waveguides with Inhomogeneous Filling, RUDN Journal of Mathematics, Information Sciences and Physics 26 (2) (2018) 129-139, in Russian.
16. F. Hecht, New Development in FreeFem++, J. Numer. Math. 20 (3-4) (2012) 251-265.
17. J. Love, Theory of Elasticity, GTTI, 1939, in Russian.
18. G. Duvaut, J.-L. Lions, Les inequations en mecanique et en physique, Dunod, Paris, 1972.
19. F. Stummel, Rand- und Eigenwertaufgaben in Sobolewschen Raumen, Springer, Berlin-Heidelberg-New York, 1969.
20. W. C. Chew, Lectures on Theory of Microwave and Optical Waveguides (2012). URL http://wcchew.ece.illinois.edu/chew/course/tgwAll20121211.pdf
21. V. M. Babich, V. S. Buldyrev, Short-Wavelength Diffraction Theory: Asymptotic Methods, Springer, Berlin, 1991.
© Малых М.Д., 2018
This work is licensed under a Creative Commons Attribution 4.0 International License Для цитирования:
Малых М.Д. О нормальных модах закрытого волновода с разрывным заполнением // Вестник Российского университета дружбы народов. Серия: Математика. Информатика. Физика. — 2018. — Т. 26, № 4. — С. 321-330. — DOI: 10.22363/23129735-2018-26-4-321-330.
For citation:
Malykh M. D. On Normal Modes of the Closed Waveguide with Discontinuous Filling, RUDN Journal of Mathematics, Information Sciences and Physics 26 (4) (2018) 321-330. DOI: 10.22363/2312-9735-2018-26-4-321-330. In Russian.
Сведения об авторах:
Малых Михаил Дмитриевич — кандидат физико-математических наук, доцент кафедры прикладной информатики и теории вероятностей РУДН (e-mail: malykh-md@ rudn.ru, тел.: +7(495)9522572)
Information about the authors:
Malykh Mikhail D. — Candidate of Physical and Mathematical Sciences, assistant professor of Department of Applied Probability and Informatics of Peoples' Friendship University of Russia (RUDN University) (e-mail: [email protected], phone: +7(495)9522572)
