Содержание к диссертации
Введение
Глава 1. Распространение электромагнитных полей в неограниченных средах с периодически изменяющимся по одной из декартовых координат показателем преломления 27
1.1 Введение 27
1.2 Н - и Е - постановки дисперсионных задач 28
1.3 Расчет дисперсии плоских волн в Я- постановке 33
1.4 Расчет дисперсии плоских волн в Е- постановке 53
1.5 Расчет дисперсии Н-иЕ- волн 64
1.6 Обсуждение результатов 82
1.7 Выводы 88
Глава 2. Симметричные волны круглого диэлектрического волновода с периодически изменяющимся вдоль оси показателем преломления 90
2.1 Введение 90
2.2. Постановка краевой задачи 92
2.3. Составление дисперсионных уравнений 95
2.4. Результаты расчета дисперсионных зависимостей 101
2.5. Обсуждение результатов 110
2.6. Выводы 128
Глава 3. Электромагнитные волны слабонаправляющего градиентного световода с периодически изменяющимся вдоль оси показателем преломления 131
3.1 Введение 131
3.2 Постановка краевой задачи 132
3.3 Дисперсионное уравнение 140
3.4 Расчет дисперсионных зависимостей, обсуждение результатов 142
3.5 Выводы 172
Глава 4. Диэлектрические волноводы с импедансными поверхностями 174
4.1 Введение 174
4.2 ДВ с периодически изменяющимся поверхностным импедансом 175
4.2.1 Постановка краевой задачи 175
4.2.2 Составление дисперсионного уравнения 178
4.2.3 Численные результаты 182
4.3 Круглый ДВ со спирально проводящей поверхностью 187
4.3.1 Постановка краевой задачи 187
4.3.2 Составление дисперсионного уравнения 189
4.3.3 Численные результаты 192
4.4 Круглый ДВ, покрытый резистивной пленкой 202
4.4.1 Постановка краевой задачи 202
4.4.2 Составление дисперсионного уравнения 204
4.4.3 Численные результаты 205
214 217 226
4.5 Выводы 213
Заключение
Список литературы
Приложение
- Расчет дисперсии плоских волн в Е- постановке
- Результаты расчета дисперсионных зависимостей
- Дисперсионное уравнение
- ДВ с периодически изменяющимся поверхностным импедансом
Введение к работе
Открытые направляющие структуры, в частности, диэлектрические волноводы (ДВ) находят широкое применение как линии передачи, а их отрезки как базовые элементы различных устройств во всех участках высокочастотного диапазона. В диапазонах СВЧ и КВЧ на основе ДВ строятся такие функциональные узлы, как линии задержки, антенны бегущей волны, открытые диэлектрические резонаторы, антенные облучатели. На основе диэлектрических волноводов с резистивными слоями строятся фиксированные и поляризационные аттенюаторы, согласующие устройства. Периодически -нерегулярные ДВ используются при создании различных частотно избирательных устройств. Слоистый круглый ДВ является строгой математической моделью оптического волокна и декомпозиционной базой различных функциональных узлов оптического диапазона.
Диссертация посвящена исследованию особенностей распространения электромагнитного поля в немагнитных неограниченных и ограниченных средах, диэлектрическая проницаемость которых периодически изменяется в направлении распространения поля (вдоль оси z). Примером открытой направляющей структуры с периодической зависимостью є от продольной координаты является оптическое волокно, на базе которого выполняются решетки Брэгга, используемые при создании частотно-избирательных устройств оптического диапазона. Кроме того, в диссертации исследуются особенности распространения электромагнитного поля в периодически-нерегулярных диэлектрических волноводах и в круглом однородном ДВ, покрытом тонкой резистивной пленкой.
Актуальность темы: Успехи в области оптических технологий приводят к тому, что телекоммуникационные системы переходят к передаче данных в магистральных сетях по оптоволокну. Целесообразность подобного перехода стимулируется высокой пропускной способностью таких сетей [1 - 8] и большой скоростью передачи информации.
Сегодня оптическое волокно достаточно широко используется и в телекоммуникационных, и в компьютерных сетях любого масштаба. Об оптических сетях, как о самостоятельной технологии, заговорили в связи с освоением метода плотного мультиплексирования с разделением по длине волны, а также с разработкой функциональных узлов оптического диапазона, таких как оптоволоконные широкополосные оптические усилители, оптические фильтры, достаточно точные волновые демультиплексоры, оптические мультиплексоры ввода/вывода каналов, узкополосные лазеры и ряда других [1 -7].
Брэгговские решетки широко используются в системах связи. В лазерной технике они применяются для достижения одномодовости. Волоконные Брэгговские решетки используются в различных устройствах оптического диапазона таких, как фильтры, мультиплексоры, компенсаторы дисперсии. Их основные достоинства - низкие потери, легкость соединения с другими участками волоконного тракта, низкий температурный коэффициент длины, простая конструкция, дешевизна.
Актуальность проводимых исследований определяется отсутствием методик, позволяющих производить теоретические расчеты характеристик распространения волн в волоконных структурах с периодически изменяющимся в продольном направлении показателем преломления. Создание таких методик позволит разрабатывать новые устройства и совершенствовать имеющиеся.
Оптоволоконная решетка Брэгга [3, 5, 6] представляет собой участок волокна, в сердцевине которого коэффициент преломления периодически изменяется вдоль направления распространения волны. Эти изменения можно вызвать воздействием ультрафиолетового излучения, прикладываемого с помощью интерферометра или фазовой маски. Таким способом получают пространственную дифракционную решетку, позволяющую разрешить главные максимумы дифракционной картины для различных длин волн [3, 5, 6]. Схематически волоконная решетка Брэгга изображена на рис. В.1.
Принцип ее действия заключается в следующем: две волны, распространяющиеся в одном направлении с фазовыми постоянными /?, и /?2, будут обмениваться энергией, если выполняется условие Брэгга: где d - период волоконной решетки.
В отражающем фильтре (рис. В.1) энергия прямой волны переходит в энергию рассеянной волны, распространяющейся в обратном направлении и имеющей ту же длину Л, при условии: |/?-(-/?)| = 2./? = ^ _ 2-л-п.д-
Полагая /? = —, где Я0 - длина волны падающего света в вакууме,
К neff - эффективный показатель преломления волновода или волокна, получаем, что волна отражается при условии: Л =2-neff-d.
Л0 называется Брэгговской длиной волны. Эффективность отражения уменьшается по мере расхождения длины волны падающего света с Брэгговской длиной волны. При прохождении через решетку света сложного спектрального состава отражается только Брэгговская длина волны, в то время как остальные длины волн проходят, не отражаясь [7].
Волоконные решетки классифицируются на короткопериодные и длиннопериодные [7]. Короткопериодные решетки также называются Брэгговскими решетками и имеют период порядка 0,5 мкм. Период короткопериодной решетки сравним с длиной волны. Длиннопериодные решетки имеют период от нескольких сотен мкм до нескольких мм. Период
Решетка Брэгга «олочка
Сердцевина
Модуляция показателя прел обілення
Рис. В.1 Однородная Брэгговская решетка с постоянными амплитудой изменения показателя преломления и периодом. длиннопериодной решетки намного больше длины волны. В настоящее время созданы волоконные Брэгговские решетки с чрезвычайно малыми потерями (0,1 дБ), высоким соответствием геометрических размеров решетки заданным (отклонение от заданных размеров не превышает ±0,05 нм) и высоким подавлением помех соседнего канала (40 дБ) [6]. В диссертационной работе предлагается методика описания физических процессов, происходящих в волоконной решетке Брэгга, не зависящая от геометрических параметров исследуемой решетки.
Такой параметр, как температурный коэффициент длины, характеризует изменение длины волокна в зависимости от температуры. У Брэгговских решеток он обычно составляет 1,25-10"2 нм/С. Однако, это изменение можно компенсировать, добавляя в среду решетки материалы, которые имеют отрицательный температурный коэффициент увеличения длины. Такие пассивные температурно-компенсированные решетки имеют температурные коэффициенты длины около 0,07 10"2 нм/С. Это обеспечивает очень малый сдвиг центральной длины волны в эксплуатационном диапазоне температур (100 С) и означает, что Брэгговские решетки могут работать фактически без контроля температуры [6]. Такие свойства решеток делают их очень привлекательными для выполнения на их основе различных функциональных узлов оптического диапазона.
Несмотря на большой интерес к волоконным решеткам, вопрос получения (записи) закона изменения их показателя преломления до сих пор не закрыт. В этом направлении постоянно ведутся исследования [8 - 12]. Запись волоконных решеток основана на использовании фоточувствительной сердцевины оптического волокна. Обычное кремниевое волокно при добавлении примеси германия становится чрезвычайно фоточувствительным. Подвергая это волокно воздействию ультрафиолетового излучения (УФ), можно вызвать изменение показателя преломления материала сердцевины. В таком волокне решетка может быть создана с помощью облучения волокна двумя интерферирующими ультрафиолетовыми пучками. Наведенная на сердцевину интерференционная картина вызывает периодические изменения показателя преломления вдоль волокна. Там, где интенсивность УФ излучения высокая, показатель преломления увеличивается, а где она мала - показатель остается без изменений. Требуемое для получения решеток изменение показателя преломления достаточно мало - около 10"4. Для производства решеток также могут быть использованы другие методы, например, применение фазовых масок [6]. Фазовая маска является дифракционным оптическим элементом. С ее помощью падающее на волокно УФ излучение "расщепляется" на лучи, которые, интерферируя между собой, прочерчивают решетку внутри волокна [6].
Недавно впервые была осуществлена запись решеток показателя преломления в германосиликатном световоде непрерывным излучением Аґ-лазера УФ диапазона (333-364 нм) [12]. Индуцированный показатель преломления в световоде с молекулярным 10% содержанием Ge02 в сердцевине составил 1.9-10-4 при плотности мощности излучения 1.7-105 Вт/см2. Полученные решетки обладают примерно такой же термостабильностью, что и решетки, записанные с помощью KrF - лазера (248 нм). Можно отметить, что излучение ближнего УФ диапазона может быть использовано для записи длиннопериодных решеток показателя преломления в германосиликатных световодах. Так как максимально индуцированный в сердцевине показатель преломления составил 1.9-10"4, то широко распространенный Аґ - лазер без удвоения частоты может быть использован для изготовления фотоиндуцированных решеток различных типов с большим периодом. Эти исследования показали, что данный способ фотовозбуждения может быть использован и для записи Брэгговских решеток, однако тогда для повышения когерентности УФ излучения необходимо выделение одной лазерной линии. При этом значительно снижается интенсивность УФ излучения и, по-видимому, эффективность наведения показателя преломления. Из проведенных экспериментов следует, что даже значительное увеличение времени облучения не позволяет компенсировать понижение интенсивности записывающего излучения.
Важным вопросом, возникающим при решении задач, связанных с передачей информации через волоконные решетки, является вопрос о влиянии отклонений геометрических параметров решетки, а также параметров среды от заданных на характеристики сигнала, передаваемого через данные решетки [13 - 16]. В частности, отклонения параметров решетки от расчетных могут привести к фазовым и амплитудным погрешностям передаваемого сигнала. Характеристики оптических узлов обычно деградируют за счет случайных флуктуации геометрических размеров решетки, а также за счет отклонения показателя преломления от заданного значения. Рэлеевское рассеяние, например, является типичным примером рассеяния, обусловленного случайными флуктуациями показателя преломления. Случайные колебания геометрии устройства, шероховатость поверхностей и хаотичность в периодичности решеток приводят к увеличению рассеяния во многих оптических узлах.
Таким образом, на основании вышеизложенного можно утверждать, что волоконные Брэгговские решетки, широко применяемые в волоконно-оптической связи, требуют высокой точности изготовления. Кроме того, они должны быть невосприимчивы к поляризационному воздействию, в них нежелательно отклонение показателя преломления от заданного значения. Решетка также должна быть лишена двойного лучепреломления и дихроизма. Поэтому процесс изготовления решетки должен сопровождаться теоретическим анализом ее характеристик с целью недопущения перечисленных выше нежелательных явлений.
Для анализа характеристик волоконных решеток Брэгга в процессе и после их изготовления используется теория оптических цепей [17]. Данная теория является аналогичной теории электрических цепей. Теория оптических цепей основана на неинтерферометрическом и интерферометрическом методах измерения характеристик оптических устройств.
Неинтерферометрические методы построены на прямом эксперименте. У передатчика модулируется интенсивность выходного луча, а приемник фиксирует прошедший через оптическое устройство сигнал. В данном эксперименте по связи мощности на выходе передатчика и мощности на входе приемника получают электрическую передаточную функцию, по которой вычисляется величина оптической передаточной функции и групповое время задержки. Данная схема измерений может быть дополнена устройствами, позволяющими анализировать поляризацию света, проходящего через оптическую систему. Однако чувствительность неинтерферометрического метода ограничена погрешностями прямых измерений.
Интерферометрическое измерение с фиксированной поляризацией в схеме с одним передатчиком является самой простой формой измерения в теории оптических цепей. Интерферометрический метод измерения с низкой когерентностью света позволяет получить импульсную характеристику устройства с высоким пространственным разрешением. Напротив, интерферометрический метод измерения с высокой степенью когерентности света позволяет определить коэффициент отражения, по которому может быть рассчитана импульсная характеристика. Пространственное разрешение ограничено диапазоном перестройки лазера по длине волны. Несмотря на ограниченность перестройки лазера, решетки с большой длиной могут быть исследованы благодаря высокой когерентности одномодового лазера. В любом случае импульсная характеристика позволяет определять структуру изотропной решетки. Полученные результаты измерений подвергаются анализу, по результатам которого можно исправлять ошибки изготовления оптического устройства. Для исправления показателя преломления оптического устройства применяют локальное ультрафиолетовое излучение.
Для расчета некоторых типов профилей Брэгговских решеток в оптическом диапазоне может быть применен так называемый транспортный метод (Transport method), широко использующийся в общей теории цепей при расчетах характеристик электронных приборов и устройств, работающих в СВЧ - диапазоне [6, 18, 19]. Данный метод обеспечивает хорошую точность практических расчетов Брэгговских волоконных решеточных структур. К сожалению, подробное его описание фактически невозможно встретить в широко доступной научной литературе, т.к. практические результаты, которые достигаются с его помощью, имеют большую коммерческую стоимость, поэтому заинтересованные в подобных исследованиях организации стараются не предавать их широкой огласке, а используют их исключительно для служебного пользования.
В настоящее время разрабатывается достаточно много оптических устройств, у которых показатель преломления гармонически изменяется вдоль направления распространения волн. Постоянно совершенствуются методы изготовления устройств такого типа. Создаются методы контроля их геометрических и физических параметров. Однако вопросы, связанные с математическим, описанием процессов прохождения волн в таких структурах, находятся в начальной стадии развития. Очевидно, что разработка методов расчета ДВ с периодически изменяющимися вдоль оси параметрами является сегодня одним из самых перспективных направлений прикладной электродинамики, которые способны продвинуть развитие и совершенствование элементной базы СВЧ, КВЧ и оптического диапазонов волн.
Целью диссертации является:
Разработка общей методики расчета характеристик распространения электромагнитных полей в неограниченных средах с периодически изменяющимся по одной из декартовых координат показателем преломления, в круглых открытых ДВ с периодически изменяющимся вдоль оси показателем преломления (в том числе градиентных); исследование особенностей распространения электромагнитных волн в открытых направляющих диэлектрических структурах с периодически изменяющимся вдоль их оси показателем преломления; исследование особенностей распространения волн в периодически-нерегулярных круглых ДВ и в диэлектрическом волноводе, покрытом тонкой резистивной пленкой; создание эффективных алгоритмов и программ, позволяющих проводить электродинамический расчет дисперсионных характеристик и характеристик затухания волн указанных направляющих структур.
Методы исследования.
Представленные в диссертационной работе теоретические результаты получены на основе методов Бубнова-Галеркина [20 - 25], метода частичных областей (МЧО) [26 - 37] и метода поверхностного тока [38 - 41] (МПТ).
Алгоритмы, созданные на основе методов Бубнова-Галеркина, МЧО, и МПТ, удобны для использования в системах автоматизированного проектирования (САПР) функциональных узлов СВЧ, КВЧ и оптического диапазонов волн, ввиду их универсальности и простоты алгебраизации функциональных уравнений, получаемых в результате реализации граничных условий.
Научная новизна. В диссертационной работе: предложен общий подход к исследованию характеристик распространения волн в ограниченных и неограниченных диэлектрических средах с периодически изменяющимся в направлении распространения показателем преломления; исследованы характеристики плоских волн в неограниченных средах с периодически изменяющимся в направлении распространения показателем преломления; исследованы характеристики Н - и Е - волн, распространяющихся в неограниченных средах с периодически изменяющимся показателем преломления; исследованы характеристики Н - волн, в волоконных световодах с периодически изменяющимся в направлении распространения показателем преломления; исследованы характеристики волн, распространяющихся в слабонаправляющих градиентных световодах с периодически изменяющимся вдоль оси показателем преломления. предложена методика расчета характеристик распространения в круглых периодически-нерегулярных ДВ; исследованы особенности распространения волн в круглом ДВ со спирально проводящей поверхностью и в диэлектрическом волноводе, покрытом резистивной пленкой.
Обоснованность и достоверность научных положений, выводов и рекомендаций, сформулированных в диссертации, подтверждается: использованием при расчете направляющих структур теоретически обоснованных методов; сравнением численных результатов, полученных различными методами; численной проверкой выполнения предельных переходов от рассматриваемых структур к структурам, решения краевых задач для которых достоверно известны; проверкой полученных результатов на сходимость.
Практическая ценность работы заключается: в создании алгоритмов и программ, позволяющих производить расчет дисперсионных характеристик волн, распространяющихся в волоконных световодах с периодически изменяющимся вдоль направления распространения волн показателем преломления; в создании алгоритмов и программ, позволяющих производить расчет дисперсионных характеристик волн, распространяющихся в слабонаправляющих градиентных световодах с периодически изменяющимся вдоль оси показателем преломления; в создании алгоритмов и программ, позволяющих производить расчет дисперсионных характеристик волн, распространяющихся в периодически-нерегулярных ДВ; в создании алгоритмов и программ, позволяющих производить расчет дисперсионных характеристик волн, распространяющихся в ДВ, покрытых резистивной пленкой.
Указанные алгоритмы и программы являются основой для создания системы компьютерного проектирования функциональных узлов СВЧ, КВЧ и оптического диапазонов волн.
Реализация и внедрение результатов.
Пакеты программ переданы в ННИПИ "Кварц", ФГУП НИИИС им. Седакова, Институт химии высокочистых веществ РАН.
Положения, выносимые на защиту:
Обоснование применения метода Галеркина для расчета характеристик волн, открытых направляющих диэлектрических структур с периодически изменяющимися в направлении распространения волн параметрами.
Модификация формулы Брэгга для открытых направляющих диэлектрических структур СВЧ, КВЧ и оптического диапазона с периодически изменяющимся в направлении распространения волн показателем преломления.
Приближенный метод расчета характеристик волн круглого открытого ДВ с периодически изменяющимся вдоль оси показателем преломления.
Результаты расчета дисперсии волн ДВ с периодически изменяющимся вдоль направления их распространения показателем преломления.
Постановка и решение дисперсионной задачи для волн в слабонаправляющем градиентном ДВ с периодически изменяющимся вдоль оси показателем преломления.
Результаты исследования дисперсионных свойств волн слабонаправляющего градиентного ДВ с периодически изменяющимся вдоль оси показателем преломления.
Формулировки краевых задач и алгоритмы расчета характеристик распространения волн периодически-нерегулярных ДВ и импедансных открытых направляющих структур СВЧ и КВЧ диапазонов.
Апробация работы.
Результаты диссертационной работы докладывались и обсуждались на:
Всероссийской научно-технической конференции "Информационные системы и технологии. ИСТ - 2002", Н.Новгород, 2002;
II Международной научно-технической конференция "Физика и технические приложения волновых процессов", Самара, 2003;
Всероссийской научно-технической конференции "Информационные системы и технологии. ИСТ-2004", Н.Новгорд, 2004;
III Международной научно-технической конференция "Физика и технические приложения волновых процессов", Волгоград, 2004.
Всероссийской научно-технической конференции "Информационные системы и технологии. ИСТ-2005", Н.Новгорд, 2005;
IV Международной научно-технической конференция "Физика и технические приложения волновых процессов", Н.Новгород, 2005.
Краткое содержание работы
Во введении проводится анализ современного состояния вопроса, ставится цель диссертационной работы, обосновывается ее актуальность, формулируются задачи исследований, определяется новизна полученных результатов и их практическая ценность, формулируются основные положения, выносимые на защиту, кратко излагается содержание диссертации.
В первой главе диссертации: приводятся результаты исследования дисперсионных свойств волн в неограниченных периодически неоднородных в направлении распространения волн диэлектрических средах в Н- и Е— постановках, полученные с использованием метода Бубнова-Галеркина. Методика составления дисперсионных уравнений сводится к следующим операциям. Записав уравнения Максвелла для случая распространения гармонической волны в изотропной электронейтральной среде без потерь и предположив отсутствие зависимости компонент поля от одной из поперечных координат, приходим к волновому уравнению относительно электрической (магнитной) компоненты волны. Используя метод разделения переменных, получаем дифференциальные уравнения относительно функции Z(z), описывающей продольную зависимость компоненты электрического (магнитного) поля. Разложив функции Z(z) и e(z) в ряды: mss-oo m-—co где у — постоянная распространения волны в бесконечной периодической среде, подставив эти ряды в полученное дифференциальное уравнение, с использованием процедуры Галеркина получаем бесконечную систему однородных алгебраических уравнений относительно неизвестных коэффициентов Ът. Записывая условие нетривиальности решений этой системы, получаем характеристическое уравнение относительно у.
Результаты решения этого уравнения показывают, что нормированная постоянная распространения линейно зависит от частоты не во всем частотном диапазоне. На дисперсионных характеристиках имеются горизонтальные участки, которым соответствует резкое (резонансное) возрастание нормированной постоянной затухания. Расположение т -го горизонтального участка на частотной зависимости фазовой постоянной и соответствующего ему участка резкого возрастания коэффициента затухания определяется соотношением: _4-neff.d/Ат~ /m + V где Лт - длина волны падающего света в вакууме (брэгговская длина волны), neff = 1 - эффективный показатель преломления среды, 2 d - период структуры, т = 0,1,2,.... Это условие брэгговского отражения волны от периодической структуры [5].
В главе рассматривается случай распространения плоских волн двух поляризаций под углом к оси периодичности среды. При этом плоские волны образуют поля с зависимостью от поперечной координаты, соответствующие волнам типа Ни Е, если связывать классификацию с осью z, вдоль которой периодически изменяется показатель преломления среды. Приводятся результаты расчета дисперсии и затухания.
Корректность постановки дисперсионной задачи и правильность расчетов проверяются по выполнению предельных переходов от периодически неоднородных сред к практически однородным средам (є2«є1) и к средам с малыми изменениями параметров на расстояниях порядка длины волны {d » Л). Достоверность результатов проверяется также сходимостью решений по волновым числам и по числу собственных функций, учитываемых в представлениях полей.
Во второй главе диссертации: приводятся результаты исследования дисперсионных свойств симметричных волн типа Н в круглом диэлектрическом волноводе с периодически изменяющимся вдоль оси показателем преломления. Дисперсионные уравнения составляются с использованием методов Бубнова-Галеркина и частичных областей (МЧО). При их решении применяется итерационный процесс.
В данном случае решается система двух трансцендентных уравнении. Первое из них составляется на основе дифференциального уравнения,' описывающего продольную зависимость поля, с использованием описанной выше процедуры Галеркина. Второе получается из граничных условий на поверхности ДВ.
Записав условия неразрывности тангенциальных компонент поля на границе раздела двух диэлектрических сред, подставив в эти условия выражения для электрической и магнитной компонент поля и приравняв в полученных уравнениях коэффициенты при одинаковых гармониках, получаем характеристические уравнения относительно а и хт-> гДе
Хт=лієц-М'<»2- -і-у 1 Поскольку все Хт связаны с продольным волновым числом р, возникает неоднозначность решения дисперсионной задачи, то есть мы получаем множество характеристических уравнений относительно а и р при различных т. Полагая, что продольная периодичность внутренней области слабо влияет на характер продольной зависимости поля во внешней области, то есть в этой области доминирует нулевая пространственная гармоника, из всей системы указанных характеристических уравнений оставляем лишь одно - с т = О. Таким образом, получаем систему двух трансцендентных уравнений относительно поперечного для внутренней области волнового числа а и постоянной распространения /. Для ее решения используется итерационный процесс. На основе составленного алгоритма были получены численные результаты.
Анализ численных результатов показал, что у частотных зависимостей фазовых постоянных волн в нижней части частотного диапазона имеются горизонтальные участки. В области частот, соответствующей этим "полочкам", у постоянной затухания наблюдаются всплески. Качественный вид частотных зависимостей для фазовой постоянной и постоянной затухания для случая ограниченной среды не изменился по сравнению со случаем неограниченной среды. Сохранилось и месторасположение "полочек" на частотной оси.
Анализ дисперсионных зависимостей показал, что центральные точки * . горизонтальных участков дисперсионных характеристик различных волн несколько смещены друг относительно друга, отсюда можно утверждать, что для различных волн эффективные показатели преломления в принципе разные, (neff=n$), но близкие по величине. Значение эффективного показателя преломления приближенно можно определить как: n,eff = c/v , где с - ' / phase скорость света в окружающей среде, Vphase - фазовая скорость волны, направляемой волокном.
Критические частоты волн, распространяющихся в волоконном световоде с периодическим изменением показателя преломления вдоль оси распространения, близки к критическим частотам волн H0q обычного h (регулярного) волокна, определяемым как: со"^" = — ч =, где h0 - q -й а-уІМєі-єіі) корень уравнения: J0(a-a) = 0. Благодаря этому при нахождении критической частоты волны Я01, на первом шаге итерационного процесса берем hn =3.832, для волны Н02 - /712=7.02, Я03 - hn =10.17, и так далее. Численные отличия значений критических частот, определяемые по формуле, справедливой для обычного (регулярного) волокна, объясняются периодичностью направляющей структуры.
Были получены частотные 'зависимости коэффициентов замедления основной т = 0 и высших т = ±1,±2 гармоник. Из этих зависимостей следует, что с ростом частоты коэффициент замедления Ру, всех гармоник
є асимптотически стремится к значению пп = у . Таким образом, на высоких частотах распространение волн в волокне с периодически изменяющимся вдоль оси показателем преломления происходит в целом по тем же канонам, что и в регулярном диэлектрическом волноводе. При этом, поскольку фазовые скорости гармоник с ростом их номера уменьшаются, коэффициенты замедления высших гармоник стремятся к асимптоте сверху.
На критической частоте коэффициент замедления основной (я?=0) гармоники, соответствующей волне H0q, как и должно быть, равен единице.
Были рассчитаны частотные зависимости действительной части л,л. \а2+Р2эффективного показателя преломления, определяемого как: neff = —=-, для у є j ju-co трех вышерассмотренных пространственных гармоник. Формула n,eff = c/v / phase определят физический смысл эффективного показателя преломления. При расчетах использовалась первая из вышеприведенных формул. Результаты, получаемые при использовании и той, и другой формул, одинаковы. Из полученных дисперсионных зависимостей видно, что эффективный показатель преломления волокна для основной гармоники (т = 0) не зависит от частоты, для других гармоник имеет место явно выраженная частотная зависимость. Причем для гармоник с высокими номерами она во всем частотном диапазоне имеет аномальный характер.
Достоверность полученных результатов проверена сходимостью решений по волновым числам и по числу собственных функций, учитываемых в представлениях полей.
Предложенный метод расчета дисперсионных характеристик волн волокна с периодически изменяющимся вдоль оси показателем преломления строго применим лишь для симметричных волн магнитного типа. Однако с помощью него можно установить общие физические закономерности распространения волн в периодически-нерегулярных волоконных световодах, которые характерны и для симметричных волн типа Е, и для гибридных волн.
В третьей главе диссертации: описываются постановка и решение краевой задачи для слабонаправляющего периодически-нерегулярного градиентного волокна с усеченным параболическим профилем. Дисперсионное уравнение составляется на основе процедуры Бубнова-Галеркина с представлением поля в направляющей структуре набором пространственных . гармоник. Приводятся дисперсионные зависимости основной LP00 волны для различных типов волокон (со слабой, средней и сильной градиентностью). Показывается, что, когда волокно обладает слабой градиентностью, центры частотных интервалов всплеска затухания строго соответствуют условию Брэгга. С увеличением градиентности центры указанных интервалов смещаются в область более низких частот. Объясняется это следующим образом: согласно формуле Брэгга значения длины волны, соответствующие максимальному отражению от периодической направляющей структуры, определяются [5] как К =****/*** гДе %~ эффективный показатель преломления сердцевины волокна, т - номер гармоники. С увеличением градиентности эффективный показатель преломления возрастает, что приводит к увеличению длины волны, соответствующей максимальному отражению (затуханию в периодической направляющей структуре). Приводятся дисперсионные зависимости для LPlp волн при различных значениях / и q для волокна со слабой градиентностью.
Представляемые результаты демонстрируют хорошую сходимость решений: амплитуды гармоник с ростом их номера быстро и равномерно убывают. Это подтверждает корректность выбранного представления продольной зависимости поля и алгебраизации дисперсионной задачи.
В четвертой главе диссертации: описываются постановка и решение краевых задач для диэлектрических волноводов с импедансными поверхностями.
В главе рассматривается круглый открытый ДВ с периодически изменяющимся поверхностным импедансом, круглый ДВ со спирально проводящей поверхностью и круглый ДВ, покрытый резистивной пленкой.
Процедуры составления дисперсионных уравнений в рассматриваемых задачах однотипны и сводятся к следующим операциям. Уравнения Гельмгольца в общем случае решаются относительно продольных компонент обоих векторов Герца. По найденным значениям векторов Герца (электрического и магнитного) для первой и второй областей находятся выражения для векторов напряженности электрического и магнитного полей. Найденные значения векторов напряженностей электрического и магнитных полей подставляются в граничные условия. Запись граничных условий: - для круглого открытого ДВ с периодически изменяющимся поверхностным импедансом имеет вид: Hzl (г = a) = Hz2 (г = а); Н^ (г = а) = Н^ (г = а); Ez2(r = a) = W(z)-H92(r = a); E^(r = a) = -W(z)-Hz2(r = a), где Щ^) = ^/є(г) - (2-я }импеданс внешней поверхности, є(г) = є1 +A-cos z ; - для круглого ДВ со спирально проводящей поверхностью: Esl(r = a) = Es2(r = a) = 0; Hx](r = а) = Hs2(r = a); Ezl(r = a) = Ez2(r = a);
Е„1(г = а) = Е<р2(г = а), где Еж = Ег -sin^ + Я,, -cosy/; Hs= H2-smy/ +H^-cosy/, sin ц/ = , d - шаг спирали, a - радиус ДВ;
2-я-а - для круглого ДВ, покрытого резистивной пленкой (параметры пленки а и єтаковы, что а»є-а>): Ezl(r = a) = Ez2(r = a); Epl(r = a) = Ep2(r = а); Hzl(r = а)-Hz2(r = а) = Д-сг-Е^(г = а); Н^г = а)-Н^(г = а) = -Д-а-Ел(г = а), где д-сг - поверхностная проводимость резистивной пленки, приводит к системам линейных однородных алгебраических уравнений. Приравнивая к нулю их главные определители, получаем дисперсионные уравнения волн, распространяющихся в рассматриваемых структурах.
На основе составленных алгоритмов были получены численные результаты в виде дисперсионных зависимостей волн, распространяющихся в рассматриваемых направляющих структурах.
Для ДВ с периодически изменяющимся поверхностным импедансом расчеты проводились в нулевом приближении. В связи с этим полученные, частотные зависимости соответствуют частотным зависимостям постоянных распространения нулевых гармоник волн. Из них следует, что постоянные распространения основных (т = 0) гармоник во всем частотном диапазоне являются чисто действительными величинами (Д, (й)) = 0), распространение основных гармоник начинается с нулевых частот. Таким образом, основные > гармоники не имеют критических частот. На низких частотах поле "прилипает". к импедансной поверхности, оказываясь слабо замедленным. При этом внешнее поле создается поверхностными магнитными токами, созданными скачком тангенциальных компонент электрического поля, имеющим место в предлагаемой математической модели. Из частотной зависимости коэффициентов замедления основных (т = 0) гармоник первых шести волн видно, что с ростом частоты коэффициенты замедления Р/ основных гармоник всех волн асимптотически стремится к значению щ = Аєу , при этом стремление к указанной величине происходит снизу. Из этого следует, что на высоких частотах в волокне с периодически изменяющимся на поверхности волокна показателем преломления имеет место диэлектрический эффект (как и в регулярном диэлектрическом волноводе), заключающийся во втягивании с ростом частоты поля в ДВ. Поперечные волновые числа II области для всех рассматриваемых волн имеют чисто мнимое значение во всем диапазоне .' частот. Этот факт подтверждает отсутствие критических частот у основных (т = 0) гармоник и говорит о том, что поле основных гармоник рассматриваемых волн имеет поверхностный характер.
Из дисперсионных зависимостей, полученных для круглого ДВ со спирально проводящей поверхностью, следует, что лишь одна волна не имеет критической частоты. Будем называть ее НЕ00. У всех остальных волн есть частоты, на которых постоянная распространения /Г обращается в нуль. Таким образом, на низких частотах волны рассматриваемой структуры уподобляются волнам экранированного волновода, образуемого спиральной линией. Из анализа частотных зависимостей постоянных затухания первых шести (по порядку следования корней дисперсионного уравнения) волн," распространяющихся в круглом ДВ со спирально проводящей поверхностью, следует, что на низких частотах все волны, за исключением первой, имеют всплеск затухания (/Г*0). Можно отметить, что с увеличением номера волны увеличивается и амплитуда этого всплеска. В области частот, где волны становятся затухающими, поперечные волновые числа становятся комплексными. Это говорит о том, что затухание связано с излучением. Коэффициент замедления волн с ростом частоты асимптотически стремится к значению Ту , при этом стремление к указанной величине происходит снизу.
Из анализа частотных зависимостей поперечных волновых чисел для внешней области II направляющей структуры следует, что на частотах ниже критических частот поверхностных волн мнимые части поперечных волновых чисел всех волн (за исключением первой) в области II больше нуля (lm(ar2)>0). При достижении критической частоты lm(or2) = 0. При дальнейшем увеличении частоты Im(ar2) < 0, что соответствует поверхностным волнам, у которых поле экспоненциально убывает по радиальной координате при удалении от поверхности ДВ. Таким образом, поле волн на частотах ниже критической нарастает при удалении от ДВ. Это говорит о том, что на частотах ниже критических мы имеем дело с вытекающими волнами, при достижении критической частоты вытекающие волны переходят в поверхностные.
В случае ДВ покрытого резистивной пленкой, численно были исследованы волны Е01 и ЕНи. Из частотных зависимостей видно, что при уменьшении частоты поверхностные волны переходят в быстрые волны.
Поведение решений дисперсионного уравнения для волны ЕНи имеет свои принципиальные особенности. Оставаясь в Ш-м квадранте, решения при уменьшении частоты переходят через линию fl' = k0 в область быстрых волн (^>hase >с)- При дальнейшем уменьшении частоты решения переходят в область II, где соответствуют вытекающим волнам. Частотная область у вытекающей волны ЕНп значительно шире, чем у волны Е01. На низких частотах вытекающая волна ЕНи так же, как и 01, переходит в медленную волну.
Существование быстрых поверхностных волн объясняется влиянием резистивной пленки. В обычном диэлектрическом волноводе медленные . поверхностные волны ЕНХт непосредственно переходят в вытекающие.
Численные исследования показывают, что при увеличении проводимости пленки частотные области существования вытекающих волн Е0т, ЕНЫ сужаются, минимальные значения относительной фазовой постоянной увеличиваются. Из частотных зависимостей видно, что у волны ЕНп затухание больше, чем у Е0]. Это объясняется тем, что распространение несимметричной ' волны сопровождается протеканием в пленке как продольных, так и азимутальных токов. При распространении симметричной волны в пленке протекают только продольные токи.
Показано, что нанесение на поверхность диэлектрических волноводов тонких резистивных слоев приводит к существенному изменению особенностей распространяющихся в ней волн, а в некоторых случаях и к изменению их спектра.
Результаты проведенных исследований опубликованы в работах [42 - 56].
Расчет дисперсии плоских волн в Е- постановке
Рассматриваем продольно-периодическую среду с диэлектрической проницаемостью, описываемой функцией (1.26).
Подставив выражение для диэлектрической проницаемости (1.26) в уравнение (1.23), имеем:Вывод уравнения (1.40) приведен в приложении П. 1.10.В системе уравнений, образуемой из (1.40) интегралы не берутся аналитически. Поэтому интегралы находятся численно. Система уравнений (1.40) с использованием метода редукции записывается в том или ином приближении.
Первое приближение решения системы. В этом случае т = 0,±1; к = 0,±\.На рис. 1.12 и 1.13 приведены частотные зависимости нормированных фазовых постоянных p d и нормированных постоянных затухания p"d для основной (т = 0) и высших гармоник (т = ±\). Данные зависимости рассчитаны в первом приближении, когда т = 0,±\; к = 0,±\. На рис. 1.14 и 1.15 приведены аналогичные частотные зависимости, рассчитанные во втором приближении. Все зависимости получены для среды с теми же параметрами, что и в случае Н- волн, рассмотренных выше.На рис. 1.12 и 1.14 зависимость /?0 0 и параллельные ейхарактеристики соответствуют электромагнитному полю,распространяющемуся вдоль оси z, зависимость р0 0 и параллельные ей характеристики - полю, распространяющемуся навстречу оси z. Нормированные фазовые постоянные Pmd линейно зависят от частоты не вовсем частотном диапазоне. На частотных характеристиках p d имеютсяРис. 1.13 Частотные зависимости нормированных постоянных затухания основной (т = 0) и высших пространственных гармоник (т = ±1) волны типа Е, рассчитанные в первом приближении, когда т = 0,±1; к = 0,±1, для среды с параметрами:Рис. 1.14 Частотные зависимости нормированных фазовых постоянных основной (т = 0) и высших пространственных гармоник (т = ±\,+2) волны типа , рассчитанные во втором приближении, когда m = 0,±1,±2; = 0,+1,+2, для среды спараметрами:Рис. 1.15 Частотные зависимости нормированных постоянных затухания основной (т = 0) и высших пространственных гармоник (w = ±l,±2) волны типа Е, рассчитанные во втором приближении, когда т = 0,±1,±2; = 0,±1,±2, для среды с параметрами: еу =3, єу =1, /горизонтальные участки, которым соответствует резкое (резонансное) возрастание нормированной постоянной затухания /? /,рис. 1.13 и 1.15.
Расположение горизонтальных участков частотных зависимостей фазовой постоянной и участков резкого возрастания затухания, как и для случая Я волн, определяется соотношениями (1.35), (1.36).
Анализ рис. 1.12 и 1.14, 1.13 и 1.15 показывает, что с увеличением частоты протяженность горизонтальных участков частотных зависимостей фазовой постоянной и величина затухания уменьшаются. Эффективность взаимодействия волн с периодической структурой уменьшается по мере расхождения длины волны падающего света с брэгговской длиной волны.На рис. 1.16 приведены частотные зависимости коэффициентов замедления для основной (т - 0) и высших гармоник (т = ±1,±2). Из рисунка видно, что с повышением частоты коэффициент замедлениярт = "/ - -уЩ = п,. Из этого следует, что на высоких частотахраспространение волн в среде с периодически изменяющимся вдоль направления распространения волны показателем преломления происходит по тем же правилам, что и в случае однородной диэлектрической среды.Поскольку задача о расчете характеристик распространения волн типа Е решалась в конечном приближении, т.е. с учетом конечного числа членов ряда (1.21), встает вопрос о необходимости исследования сходимости результатов решения дисперсионного уравнения, а также сходимости ряда (1.21). В таблице 1.3 приведены результаты проверки сходимости решений дисперсионных уравнений, получаемых в различных приближениях.
Рассматриваются постоянные распространения трех волн в среде соследующими параметрами: єу =3, єу =1, у =1, k0d = 0.5, d=\.Рис. 1.16 Частотные зависимости коэффициентов замедления основной (т-0) и высших пространственных гармоник (т = ±1,±2) волны типа Е, рассчитанные во втором приближении, когда т = 0,±1,±2; к = 0,±1,+2, для среды с параметрами:Из таблицы 1.3 видно, что сходимость по волновым числам для волн типа Е достигается в третьем приближении.
В таблице 1.4 приведены результаты исследования сходимости ряда (1.21). Для получения численных значений членов ряда (1.21) корень дисперсионного уравнения у подставляется в систему линейных однородных алгебраических уравнений относительно коэффициентов Ьт. Полагая b0 = 1 ипонижая порядок системы на 1, находим значения всех остальных коэффициентов Ът. Найденные значения коэффициентов Ът подставляем в (1.21), откуда находим члены ряда.
Из таблицы 1.4 видно, что последующие члены ряда гораздо меньше предыдущих, а члены ряда, учитываемые при решении задачи во втором приближении (т = ±2; к = ±2) вносят незначительную добавку в сумму ряда.
Для проверки физичности полученных результатов были рассмотрены предельные переходы к средам с є2 «є{ и к средам с малым изменением показателя преломления на расстоянии порядка длины волны (d » Л). На рис. 1.17 представлены частотные зависимости коэффициентов замедленияволн в неограниченной среде с параметрами: еу =3, єу =0,1, у = 1, d=\(случай є2«єх). На рис. 1.18 представлены частотные зависимости коэффициентов замедления волн для неограниченной среде с параметрами:представленных частотных зависимостей (рис. 1.17 и 1.18), коэффициент замедления во всех предельных случаях стремится к Je x , что соответствует однородным средам.1.5 Расчет дисперсии Н- и Е- волнВ данном разделе диссертации рассматривается распространение плоских волн двух поляризаций под углом к оси периодичности среды. В этом случае плоские волны образуют поля с зависимостью от поперечной координаты, соответствующие волнам типа Н и Е, если связывать классификацию с осью z, вдоль которой периодически изменяется показатель преломления среды.Решения уравнений (1.10), (1.14) с учетом периодичности структуры представляем в виде:
Результаты расчета дисперсионных зависимостей
На основе составленного алгоритма были получены численные результаты в виде частотных зависимостей нормированной постоянной распространения p0d, соответствующей нулевой гармонике, инормированных постоянных распространения pmd, Рт = Р0 +-—,соответствующих m-ным гармоникам. Также были рассчитаны частотные зависимости нормированного поперечного волнового числа II области хт и эффективного показателя преломления п. Смысл эффективного показателя преломления заключается в следующем: считаем, что показатель преломления в ограниченной области не изменяется по гармоническому закону, а имеет фиксированное значение. При этом на фиксированной частоте постоянная распространения волны в среде с п равна соответствующей постоянной распространения в среде с периодически изменяющимся показателем преломления.
На представленных рисунках использованы следующие обозначения:На рис. 2.2 и 2.3 представлены частотные зависимости нормированных постоянных распространения и постоянных затухания для основной (т = 0) и высших пространственных гармоник (т = ±\) волны Я01. Все зависимости Рис. 2.2 Частотная зависимость нормированной фазовой постоянной основной (/« = 0) и высших пространственных гармоник (т = ±1) волны Я01.Рис. 2.3 Частотная зависимость нормированной постоянной затухания основной (т = 0) и высших пространственных гармоник (т = ±1) волны Я0]. рассчитаны в первом приближении, когда т = 0,±1; к = 0,±1 для ДВ спараметрами: єу =3, єу =1, у =1, а = 10, d=l,. Значение, а на первомшаге итерации вычислялось как а = у . Из рис. 2.2 видно, что частотныезависимость фазовых постоянных волн имеют линейную зависимость не во всем диапазоне частот. Частотные зависимости фазовых постоянных имеют горизонтальные участки, которым соответствуют всплески постоянной затухания (рис. 2.3).
На рис. 2.4, 2.5 представлены действительная и мнимая составляющие поперечных волновых чисел, для внешней области II направляющей структуры, соответствующие дисперсионным характеристикам, приведенным на рис. 2.2 и 2.3. Они рассчитывались по найденному из дисперсионных уравнений с использованием итерационного процесса значению р, по формуле:На рис. 2.6, 2.7 представлены действительная и мнимая составляющие эффективного показателя преломления, значения которого определяются по формуле:Из рисунков видно, что с повышением частоты эффективного показатель преломления стремится к значению уЩ.
На рис. 2.8 представлены зависимости коэффициентов замедления основной и высших пространственных гармоник от частоты. Из рисунка следует, что с повышением частоты коэффициенты замедления стремятся к yj\ =nr Из этого следует, что на высоких частотах распространение волн вволокне с периодически изменяющимся вдоль оси показателем преломления происходит по тем же канонам, что и в регулярном ДВ. При этом, поскольку фазовые скорости гармоник с ростом их номера уменьшаются, коэффициенты замедления высоких гармоник стремятся к асимптоте сверху.
На рис. 2.9-2.15 изображены частотные зависимости, аналогичные рис. 2.2-2.8, но рассчитанные во втором приближении. Анализ данных рисунков показывает, что с увеличением номера приближения не происходит качественного изменения дисперсионных зависимостей, хотя имеются некоторые незначительные численные различия. Это подтверждает устойчивость постановки дисперсионной задачи.Из сравнения рис. 2.2, 2.3 и 2.9, 2.10 видно, что с увеличением номера приближения на единицу, количество пространственных гармоник возрастает на четыре.
Так как постоянные распространения волн определяются с помощью итерационного алгоритма, возникает необходимость проверки сходимости найденных решений. Такая проверка была произведена и ее результаты дляструктуры с параметрами: єу =Ъ, єу =1, у =1, 1 n dA = 1.46, а =10,Как видно из рис. 2.16 и 2.17, сходимость решений дисперсионного уравнения наступает на третьем шаге итерации, что говорит об устойчивости итерационного алгоритма.Из рис. 2.2 и 2.9 видно, что у зависимостей фазовых постоянных волн вблизи нулевых частот имеются горизонтальные участки. В области частот, соответствующей этим "полочкам", у постоянной затухания, рис. 2.3 и 2.10, наблюдаются всплески. Как видно из рис. 2.2 и 2.9, дисперсионные характеристики всех гармоник имеют горизонтальные участки, положение центральных точек которых определяется [6] формулой Брэгга:
Дисперсионное уравнение
Вывод уравнения (3.44) приведен в приложении П.3.2.Приравнивая нулю главный определитель системы (3.44), получаем дисперсионное уравнение волн рассматриваемой направляющей структуры, которое, будучи записанным при различных / и р, дает дисперсионные характеристики всех LP. волн, представляемых (в силу периодичностиструктуры) наборами пространственных гармоник. Порядок системы алгебраических уравнений (3.44), т.е. приближение, в котором решается дисперсионная задача, определяется числом учитываемых пространственных гармоник.
Раскрыв в (3.46) определитель, получаем дисперсионное уравнение относительно / в первом приближении. Аналогичным образом получаются дисперсионные уравнения относительно у в более высоких приближениях.
Постоянные распространения рт связаны с решением характеристического уравнения (3.46) у соотношениями:7ЇЇП
На основе разработанного алгоритма были рассчитаны дисперсионные зависимости основной волны LP волокна со слабой А = 0,005, среднейА. =0,05 и сильной Д. =0,5 градиентностью (правомерность расчетов для сильноградиентного волокна будет показана ниже) при Д_ =0,15, я = 10 мкм, d = 1 мкм. Расчет производился в первом приближении, когда учитывалисьгармоники т = 0,±1. Результаты расчетов приведены на рис. 3.2 - 3.10. Нарис. 3.2, 3.3 приведены нормированные значения продольных волновыхі а і ичисел основной /?0 = /L + //L и первых р±х = /?±1 + і/ЗіХ двухпространственных гармоник, представляющих поле волны P в волокне сослабой градиентностью; на рис. 3.4 - частотные зависимости коэффициентов замедления казанных гармоник. На рис. 3.5 - 3.7 и рис. 3.8-3.10 приведены аналогичные дисперсионные зависимости для волокон со средней и сильной градиентностью.
Из представленных зависимостей (рис. 3.3, 3.6, 3.9) видно, что, когда волокно обладает слабой градиентностью, центры частотных интервалов всплеска затухания строго соответствуют условию Брэгга [6]. С увеличением градиентности они (центры указанных интервалов) смещаются в область более низких частот. Объясняется это следующим образом: согласно формуле Брэгга значения длин волн, соответствующие максимальному отражению от периодической направляющей структуры, определяются [44] как:соответствующей максимальному отражению (затуханию в периодической направляющей структуре). Отметим, что для всех случаев волокон с той или иной степенью градиентности (слабой, средней или сильной) для всех значений р , найденных из дисперсионного уравнения, выполняетсяусловие:
При увеличении градиентности волокна происходит некоторое увеличение численных значений фазовых постоянных основной (т = 0) и высших (т = ±1) пространственных гармоник. Аналогичная зависимость характерна и для амплитуд всплесков затухания соответствующих гармоник.
Результаты расчетов демонстрируют общую закономерность: при со -»оо коэффициент замедления волны, образуемой всей совокупностью пространственных гармоник, стремится к значению 1 1/ (рис. 3.4, 3.7 V /є2 и 3.10). Коэффициент замедления основной (т = 0) гармоники для всех случаев градиентности волокна (слабой, средней и сильной) стремится к указанному значению снизу. На всех рисунках зависимости, расположенные над осью (k d), соответствуют прямым волнам, под осью (k d) - тем же самым волнам (направляющая структура взаимная), но распространяющимся в противоположную сторону. Для проверки корректности поставленной задачи были проведены исследования сходимости полученных решений в зависимости от числа членов ряда, представляющего зависимость поля от продольной координаты. Результаты этих исследований представлены в таблице 3.1. В таблице приведены члены ряда (3.36), полученные из системы уравнений (3.44) во 2-ом приближении при kQd = 7. Численные результаты демонстрируют хорошую сходимость: амплитуды гармоник с ростом их номера быстро и равномерно убывают. Это подтверждает корректность выбранного представления (3.36) продольной зависимости поля и методики алгебраизации дисперсионной задачи. На рис. 3.11, 3.12 представлены действительная и мнимая составляющие поперечных волновых чисел направляющей структуры при г а, соответствующие дисперсионным характеристикам, приведенным на рис. 3.2 и 3.3. Данные величины рассчитывались по найденному значению р, которое было определено из дисперсионного уравнения. Соотношение для расчета х имеет следующий вид: хт =4єігм-о)2-РІ. На рис. 3.11,3.12 введены следующие обозначения: %т = хт+іХт Из данных рисунков видно, что мнимая часть поперечного волнового числа во всем частотном диапазоне имеет отрицательное значение, то есть поле имеет поверхностный характер. Как следует из (3.34), при выполнении условия (3.2) дисперсионные зависимости волн с / 0, р 0, радиальные распределения полей которых описываются функцией (3.35), качественно должны мало отличаться от приведенных выше. Численные исследования подтверждают это. Так же это видно из рис3.2 - 3.3 и 3.13 - 3.20. При изменении значений / и р
ДВ с периодически изменяющимся поверхностным импедансом
Открытые диэлектрические структуры с различными импедансными поверхностями находят широкое применение при построении функциональных узлов СВЧ, КВЧ и оптического диапазонов волн. Гофрированные ДВ находят применение при построении антенных облучателей, частотно-избирательных, фазосдвигающих и согласующих устройств. Диэлектрические волноводы, покрытые тонкими резистивными пленками, используются при построении модовых фильтров, фиксированных поглощающих аттенюаторов, поляризационных аттенюаторов, согласованных нагрузок. Краевые задачи для ДВ с импедансными поверхностями могут иметь различные формулировки. В частности, при рассмотрении гофрированных диэлектрических волноводов могут использоваться метод частичных областей (МЧО) и импедансный метод [26 - 37]. В первом случае реализация граничных условий на поверхности ДВ осуществляется с использованием условий ортогональности собственных функций краевых задач для области гофра и прилегающих к нему слоев. При этом с учетом периодичности направляющей структуры поле в согласуемых областях представляется бесконечными суммами пространственных гармоник. Во втором случае вводится поверхностный импеданс как функция продольной координаты. В результате реализации граничных условий (в данном случае импедансных) образуется система функциональных уравнений, алгебраизация которых производится с использованием процедуры Галеркина [20 - 25].
Диэлектрические волноводы, покрытые тонкими резистивными пленками, обычно исследуются с использованием метода поверхностного тока (МПТ) [38 - 41], в котором для тангенциальных компонент магнитного поля в месте расположения пленки вводятся разрывные граничные условия.
В главе рассматриваются три задачи: ДВ с периодически изменяющимся поверхностным импедансом; круглый ДВ со спирально проводящей поверхностью; круглый диэлектрический волновод, на поверхность которого нанесена тонкая резистивная пленка. Первая задача решается импедансным методом, вторая методом частичный областей, третья -МПТ.
На рис. 4.1 показан круглый ДВ с периодически изменяющимся поверхностным импедансом. На поверхности ДВ диэлектрическая проницаемость периодически изменяется вдоль оси z: є = e(z).Запись электрического и магнитного векторов Герца, для областей I и II, соответственно, имеет вид:Hi2)(aim -f) — функции Бесселя и Ханкеля 2-го рода, соответственно [61 - 66]. Граничные условия на поверхности направляющей структуры записываем в виде:
Представление полей в виде (4.11 -4.18) предполагает возможностьвыполнения нулевого граничного условия на бесконечности(экспоненциального убывания поля во внешней области по радиальнойкоординате).
Таким образом, имеем систему 4-х линейных однородных алгебраических уравнений относительно неизвестных: А10, В10, А20, В20.
Приравнивая нулю главный определитель системы, получаем дисперсионное уравнение, которое решается совместно с уравнениями:В этом случае т = 0;±1; к = 0;±1.К уравнениям (4.28) - (4.29) добавляются аналогичные уравнения, записанные при к = ±\. Таким образом, вместо 2-х уравнений образуется 6 уравнений с неизвестными: Ах_х\ А10; Ап; А1Л\ А20; A2i; 5М; Bia\ Вп; B2_t; В20;В21. Уравнения (4.30) - (4.31) записываем трижды при к = 0;±\. При этом суммы будут содержать по три члена с номерами т = 0;±\. В результате добавляется еще 6 уравнений. Таким образом, система будет содержать 12 уравнений относительно перечисленных выше 12-ти неизвестных.
Приравнивая нулю главный определитель системы, получаем дисперсионное уравнение в 1-ом приближении, которое решается совместно с уравнениями:относительно неизвестных: а( 2), Д,, через которые выражаются остальныеПри расчетах диэлектрическую проницаемость в записи импеданса представляли в виде:На основе составленного алгоритма были получены численные результаты в виде частотных зависимостей постоянных распространения, коэффициентов замедления и нормированных поперечных волновых чисел а2 во второй области для первых шести (по порядку следования корней дисперсионного уравнения) волн ДВ со спирально проводящей поверхностью.
Все зависимости рассчитаны для ДВ с параметрами: у =3, еу =1,На рис. 4.6 приведены частотные зависимости фазовых постоянных первых шести волн, распространяющихся в рассматриваемой структуре. Из рисунка видно, что лишь одна волна не имеет критической частоты. Будем называть ее НЕЮ. У всех остальных волн есть частоты, на которыхпостоянная распространения р обращается в нуль. Таким образом, на низких частотах волны рассматриваемой структуры уподобляются волнам экранированного волновода, образуемого спиральной линией. На рис. 4.7 представлены частотные зависимости постоянных затухания первых шести волн, распространяющихся в круглом ДВ со спирально проводящей поверхностью, соответствующие фазовым постоянным, приведенным на рис. 4.6. Из рисунка видно, что на низких частотах все волны, за исключением первой, имеют всплеск затухания (р" 0). Можно отметить, что с увеличением номера волны увеличивается и амплитуда этого всплеска. В области частот, где волны становятся затухающими, поперечные волновые числа, становятся комплексными. Это говорит о том, что затухание связано с излучением.