Применение проекционных методов к исследованию волноведущих и резонансных систем с особенностями тема диссертации и автореферата по ВАК РФ 01.01.03, кандидат физико-математических наук Ерохин, Александр Игоревич

• Ерохин, Александр Игоревич
• кандидат физико-математических наук
• 2012

Оглавление диссертации кандидат физико-математических наук Ерохин, Александр Игоревич

ВВЕДЕНИЕ

ГЛАВА I АСИМПТОТИКА СКАЛЯРНЫХ СПЕКТРАЛЬНЫХ ЗАДАЧ С ГРАНИЧНЫМИ УСЛОВИЯМИ ДИРИХЛЕ И НЕЙМАНА ДЛЯ ОПЕРАТОРА ЛАПЛАСА В ОКРЕСТНОСТИ УГЛОВЫХ ТОЧЕК

§1.1 Выделение особенности решения уравнения Пуассона с граничными условиями Дирихле при наличии в области входящего угла

§1.2 Асимптотика решения спектральной задачи Дирихле в окрестности угловой точки

§1.3 Выделение особенности решения уравнения Пуассона с граничными условиями Неймана при наличии в области входящего угла

§1.4 Асимптотика решения спектральной задачи Неймана в окрестности угловой точки

ГЛАВА II СКАЛЯРНАЯ ТРЕХМЕРНАЯ ЗАДАЧА ВОЗБУЖДЕНИЯ ЭЛЕКТРОМАГНИТНЫХ ВОЛН В ВОЛНОВОДЕ ПРИ НАЛИЧИИ ВХОДЯЩИХ РЕБЕР

§2.1 Постановка задачи, доказательство существования и единственности решения

§2.2 Решение спектральной задачи с граничными условиями Дирихле методом конечных элементов с учетом особенности решения в окрестности угловой точки

§2.3 Решение задачи возбуждения электромагнитных волн неполным методом Галеркина, доказательство сходимости

§2.4 Результаты численного расчета задачи

ГЛАВА III ВЕКТОРНАЯ ТРЕХМЕРНАЯ ЗАДАЧА ВОЗБУЖДЕНИЯ ЭЛЕКТРОМАГНИТНЫХ ВОЛН В ВОЛНОВОДЕ ПРИ НАЛИЧИИ ВХОДЯЩИХ РЕБЕР

§3.1 Постановка задачи, доказательство существования и единственности решения

§3.2 Решение спектральной задачи с граничными условиями Неймана методом конечных элементов с учетом особенности решения в окрестности угловой точки

§3.3 Решение векторной задачи возбуждения электромагнитных волн неполным методом Галеркина, доказательство сходимости

§3.4 Результаты численного расчета задачи

ГЛАВА IV РАСЧЕТ РЕЗОНАТОРА С АКСИАЛЬНО-СИММЕТРИЧНЫМ ДИЭЛЕКТРИЧЕСКИМ ЗАПОЛНЕНИЕМ КУСОЧНО-ПОСТОЯННОГО РАДИУСА

§4.1 Математическая постановка задачи, доказательство существования решения, сходимость приближенного решения

§4.2 Результаты численного расчета задачи методом конечных элементов

Введение диссертации (часть автореферата) на тему «Применение проекционных методов к исследованию волноведущих и резонансных систем с особенностями»

Актуальность и практическая значимость.

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

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

Геометрические особенности могут являться следствием создания структур сложной геометрической формы или просто имеющих неоднородность поверхности [1]. Примером такой электродинамической конструкции служит волновод с входящими ребрами. Подобные системы широко используются в микроволновых устройствах и цепях [2].

Входящие углы в волноводе могут быть по многим причинам, как по техническим, например, вследствие состыковки нескольких волноводов

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

Гребенчатые волноводы обладают важными для практического применения свойствами, одним из которых является эффект аномально малого затухания некоторых типов нормальных волн при определенных соотношениях между геометрическими параметрами и длиной волны [3]. Данный эффект обнаружен экспериментально [4] и численно исследован в работах [5-7] для случая круглого волновода и в работе [8] для случая прямоугольного волновода. В работе [9] проведено аналитическое исследование эффекта аномально малого затухания некоторых типов нормальных волн в плоском гребенчатом и азимутально-гофрированном волноводах. Для плоского гребенчатого волновода наличие эффекта аномально малого затухания было подтверждено численным экспериментом [10].

Гребенчатые волноводы находят применение в качестве фильтра для пассивного микроволнового устройства или настроечного элемента в нем [11,12]. Гребенчатый прямоугольный [13] и гребенчатый цилиндрический [14] волноводы часто используются в двухмодовых полых или диэлектрических резонаторных фильтрах. Многогребенчатые прямоугольные волноводы применяются как настроечные элементы во входящих коаксиальных фильтрах [15].

С помощью входящих углов в волноводе можно также моделировать наличие дефектов в поверхности системы (царапины, изломы и др.) или щупы, возбуждающие волновод.

Очень часто математические модели, описывающие системы с физическими особенностями, имеют особенности математического характера. Так в случае волновода с входящим ребром известно, что решение краевой задачи, описывающей данную систему, имеет сингулярность, которая значительно осложняет численный эксперимент [16-18]. Решить эту проблему можно с помощью выделения этой особенности инструментами математической физики и корректировки численного алгоритма с учетом полученной особенности. Теоретическое исследование однородных по длине вол-новедущих систем, содержащих входящие ребра, проведено в работах [1921]

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

Классическая теория [22-24] дает аналитическое описание резонаторов с идеально проводящей поверхностью наиболее простых форм.

В случае наличия потерь в поверхности волновода можно использовать модель поверхности системы с импедансными граничными условиями [25]. В работах [25-27] решена задача расчета резонансных волновых чисел и добротности собственных колебаний сферического резонатора в приближении теории возмущений для случая малого по модулю поверхностного импеданса. Численный алгоритм решения данной задачи для случая произвольного по модулю комплексного поверхностного импеданса рассмотрен в работе [28]. Использованная в последней работе техника применена для расчета многослойного сферически-слоистого резонатора в работах [29-31].

С помощью выбора материала заполнения резонатора и его структуры можно существенно менять спектр частот, распределение полей, а также значительно увеличить добротность [30,31]. Свойство высокодобротности сферического резонатора может быть использовано при его применении в метрологии [32].

К резонаторам сравнительно простых форм также относятся цилиндрические резонаторы [22-24]. Среди цилиндрических резонаторов большой интерес представляют резонаторы аксиально-симметричных форм кусочно-постоянного радиуса. Помимо перечисленных свойств резонаторов, предполагается, что такие элементы могут найти применение при конструировании проводящих систем с малыми потерями, а также для генерации СВЧ-энергии [33]. В частности, такие резонаторы экспериментально исследовались на кафедре колебаний физического факультета МГУ [34].

Рассматривая особенности резонаторов и волноводов в виде заполнения материалами с различными электромагнитными свойствами, следует отметить приобретающие все большую популярность метаматериа-лы, искусственно созданные вещества, проявляющие при взаимодействии с электромагнитным полем свойства, существенно отличные от свойств природных материалов. Примером новых метаматериалов служат киральные среды - среды, содержащие киральные объекты, т.е. объекты, которые не совмещаются со своим зеркальным отображением (лист Мебиуса, неправильный тетраэдр, различные спирали [35-37]). Киральные среды исследуются в работах [38-41]. В работах [42-45] рассмотрены важные для практического применения свойства, которые могут проявляться в различных киральных средах.

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

Некоторые методы численного исследования.

Метод поперечных сечений.

В 1955 году в работе С.А. Щелкунова и П.Е. Краснушкина для расчета волноводов с плавными нерегулярностями был предложен метод поперечных сечений [47-48], впоследствии развитый в работах А.Г. Свешникова, A.C. Ильинского, Б.З. Каценеленбаума, В.П. Моденова и других авторов. В данном методе поле в любом сечении волновода представляется в виде бесконечной суммы полей нормальных волн [49-51], распространяющихся в обоих направлениях. Коэффициенты разложения поля являются функциями продольной координаты и удовлетворяют бесконечной системе обыкновенных дифференциальных уравнений первого порядка. Сходимость данного метода была доказана для слабо нерегулярных волноводов на физическом уровне строгости.

Метод частичных областей.

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

Неполный метод Галеркина.

В 1915 году для приближенного решения краевых задач для дифференциальных уравнений Б.Г. Галеркиным был предложен метод, который заключается в поиске решения в виде суммы непрерывной функции, удовлетворяющей краевым условиям задачи, и конечного ряда, члены которого берутся из полной в классе непрерывных функций системы линейно-независимых функций, удовлетворяющих однородным краевым условиям. Коэффициенты разложения находятся из условия ортогональности невязки к базисным функциям, которое приводит к системе алгебраических уравнений. Данный метод [54,55] относится к группе проекционных методов [50,51], так как приближенное решение проектируется на подпространство базисных функций. Поскольку И.Г. Бубнов предложил аналогичный метод для построения приближенного решения вариационных задач, то данный метод также носит название метода Бубнова-Галеркина.

В начале 60-х годов А.Г. Свешниковым был предложен неполный метод Галеркина [56-58], который приводит к конечной системе обыкновенных дифференциальных уравнений. Благодаря специальному выбору проекционных соотношений в неполном методе Галеркина построенное приближенное решение удовлетворяет тому же энергетическому равенству, что и точное решение, что позволяет провести строгое математическое доказательство сходимости метода. Зачастую для решения жесткой системы дифференциальных уравнений, получаемых при численной реализации неполного метода Галеркина, приходится использовать специальные трудоемкие методы, в частности метод направленной ортогонализации [59].

Метод конечных разностей.

Одним из наиболее мощных и универсальных методов приближенного решения краевых задач для дифференциальных уравнений в частных производных является метод конечных разностей [60], с успехом применяемый для моделирования волноведущих систем с физическими особенностями [61-70].

При решении нестационарных задач расчета волноведущих систем методом конечных разностей используются различные алгоритмы дискретизации уравнений Максвелла, основой которых является идея, впервые предложенная Кейном Йе в 1966 году [71,72], которая заключается в специальном выборе пространственного шаблона для представления электромагнитных полей, который наиболее полно отражает уравнения Максвелла.

Стационарные задачи расчета волноводов разделяются на два класса: спектральные задачи и задачи возбуждения. При решении спектральных задач, т.е. задач расчета однородной по длине волноведущей системы со сложной геометрией поперечного сечения и с неоднородными в поперечном направлении электромагнитными свойствами ее заполнения, строятся консервативные конечно-разностные схемы, либо применяются проекци-онно-сеточные методы. При решении задач возбуждения, т.е. задач расчета неоднородных по длине волноводов, используются различные подходы, в частности применяются методы ограничения области для случая локального характера неоднородности. Для решения полученной конечно-разностной задачи применяются прямые и модифицированные итерационные методы с комплексными итерационными параметрами, а также метод немонотонной прогонки [73,74]. Для расчета стационарной задачи также применяются различные подходы, например, основанные на замене данной задачи на близкую к ней нестационарную, для которой имеются устойчивые алгоритмы численного решения (метод параболического уравнения, метод опорной волны). Стоит также отметить подход, который заключается в регуляризации уравнения Гельмгольца, в использовании алгоритма ограничения спектра и выделения спектральной полосы и алгоритмов коррекции спектра параболического уравнения [52,75,76].

Метод конечных элементов.

Одним из наиболее популярных в последнее время методов расчета сложных волноведущих систем является метод конечных элементов (МКЭ) [40,77-81]. Математические основы этого метода были сформулированы в 1943 году известным математиком Р. Курантом [82-84].

В МКЭ приближенное решение ищется в виде разложения по базисным функциям, имеющим конечный носитель. В зависимости от геометрии рассматриваемой области происходит ее разбиение на подобласти, на которых искомая функция локально аппроксимируется кусочно-полиномиальными функциями. В случае простых областей в качестве стандартных конечных носителей используются треугольники или четырехугольники. Если априорно известно, что решение в некоторых точках может иметь особенность, в МКЭ используется нерегулярная сетка со сгущением в окрестностях данных точек, либо, если известно поведение предполагаемой сингулярности, в пространство базисных функций добавляется новая функция, хорошо аппроксимирующая особенность. Для улучшения точности приближенного решения применяются Ь-версия, р-версия или их комбинация - Ь-р версия метода конечных элементов. В случае применения Ь-версии МКЭ для повышения точности найденного результата производится сгущение сетки, в том числе и нерегулярное. При использовании р-версии для уточнения решения повышается степень используемых полиномов, но на практике оказывается, что эффективнее применять либо Ь-версию метода, либо Ь-р версию метода с невысокой степенью полиномов, обычно не выше третьей. Коэффициенты разложения приближенного решения находятся из проекционных соотношений или минимизации функционала, что приводит к системе линейных алгебраических уравнений с ленточной матрицей, т.е. к разностной схеме. Таким образом метод конечных элементов относится к группе проекционно-сеточных методов, которые являются комбинацией проекционных и конечно-разностных методов. Так как нередко при решении задач электродинамики получаемые матрицы являются разреженными симметричными, но незнакоопределенными, для их обработки приходится применять специальные устойчивые алгоритмы [85-87].

При использовании метода конечных элементов стоит отметить весьма мощный и перспективный метод атомарных функций [88-90].

При применении метода конечных элементов к задачам расчета волноведущих систем возникает проблема появления фиктивных решений [91,92], так называемых "духов", которые не соответствуют реально распространяющимся модам в волноводе. Существует несколько способов устранения несуществующих решений, которые можно разделить на два типа: апостериорный способ, когда вычисляются все решения, а потом из них удаляются "духи", и априорный, когда используется такая математическая постановка задачи или строятся такие алгоритмы ее решения, что фиктивные решения не появляются. Одним из априорных методов борьбы с нефизическими решениями является метод штрафов [93], однако показано [94], что и при применении данного метода не всегда удается решить проблему фиктивных решений. Другим априорным методом борьбы с "духами" является метод смешанных конечных элементов, который зарекомендовал себя как весьма эффективный и надежный метод расчета волноведущих систем со сложной геометрией и анизотропным заполнением [95-97].

Методы ограничения области.

При решении задач возбуждения волноведущих систем зачастую применяются методы ограничения области. Среди таких методов следует отметить метод эффективных граничных условий [52], использование бесконечных конечных элементов [77,98,99] при решении задачи МКЭ, метод эффективных интегральных граничных условий [100] и др. Отдельно следует отметить метод ограничения области, связанный с постановкой парциальных условий излучения, предложенных А.Г. Свешниковым [50,51,101], который оказался весьма эффективным при численном моделировании волноведущих систем. В случае нестационарных задач для ограничения области удобно использовать "идеально согласующиеся граничные условия" [102].

Методы исследования в данной работе.

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

Для численного исследования резонансной системы в работе применяется проекционно-сеточный метод - МКЭ. Вследствие аксиальной симметрии резонатора трехмерная задача расчета спектра сводится к двумерной, а наличие прямых углов в сечении определяет выбор конечных элементов.

Для исследования волновода применяется комбинация проекционных методов, а именно проекционного метода в виде неполного метода Галер-кина и проекционно-сеточного метода конечных элементов.

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

Содержание работы.

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

Заключение диссертации по теме «Математическая физика», Ерохин, Александр Игоревич

Заключение

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

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

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

3. На основе доказанных теорем делается вывод о целесообразности применения предложенного алгоритма к расчету рассматриваемых систем.

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

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

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

Список литературы диссертационного исследования кандидат физико-математических наук Ерохин, Александр Игоревич, 2012 год

1. Шестопалое В.П., Кириленко А.А., Рудь Л. А. Резонансное рассеяние волн. Волноводные неоднородности. Т. 2. Киев: Наукова думка, 1986.

2. Schiff В., Yosibash Z. Eigenvalues for Waveguides Containing Re-Entrant Corners by a Finite-Element Method with Superelements. IEEE Transactions on Microwave Theory and Techniques. 2000. 48, № 2. 214220.

3. Моденов В.П., Слепян Г.Я. К расчету собственных волн круглого азимутального гребенчатого волновода с учетом потерь в стенках. Радиотехника и электроника. 1984. 24, № 9. 1829-1832.

4. OlverA.D., Clarricoats P. J.B., Chong S.L. Experimental determination of attenuation in corrugated circular waveguides. Electronics Letters. 1973. 9, № 18. 424-426.

5. Clarricoats P.J.B., Saha P.K. Attenuation in corrugated circular waveguide. Electronics Letters. 1970. 6, № 12. 370-372.

6. Clarricoats P.J.B., Olver A.D., Chong S.L. Attenuation in corrugated circular waveguides. Part 1: Theory. Electrical Engineers. Proceedings of the Institution of Electrical Engineers. 1975. 122, № 11. 1173 1179.

7. Альховский Э.А., Ильинский A.C. О влиянии гофра на коэффициент затухания основной волны в гофрированном волноводе. В Прямые и обратные задачи теории антенн. 176-179. М: МГУ, 1976.

8. Родионов А.А., Раевский С.Б. Расчет дисперсионных характеристик и коэффициента затухания прямоугольного гофрированного волновода. Известия вузов. Радиоэлектроника 1977. 20, № 9. 69-73.

9. Слепян Г.Я. Нормальные волны гребенчатых волноводов с конечной проводомостью стенок. Радиотехника и электроника 1981. 26, № 9. 1833-1839.

10. Ю.Моденов В.П., Магатаев А.В. Численное исследование эффекта аномально малого затухания электрических волн плоского гребенчатого волновода. Вестн. Моск. Ун-та. Сер. 3. Физика. Астрономия. 1983. 24, № 4.

11. M Guglielmi, R.C. Molina, and A. Alvarez. Dual-mode circular waveguide filters without tuning screws. IEEE Microwave Guided Wave Lett. 1992. 2. 457^58.

12. К E. Boria, G. Gerini, and M. Guglielmi. Computer aided design of reentrant coaxial filter including coaxial excitation. IEEE MTT-S Int. Microwave Symp. Dig. 3. 1999. 1131-1134.

13. Назаров С.А., Пламеневский Б.А. Эллиптические задачи в областях с кусочно-гладкой границей. М: Наука, 1991.

14. Боголюбов А.Н., Делицын А.Л., Могилевский И.Е., Свешников А.Г. Проблема вычисления мод волноводов при наличии входящих ребер. Журнал радиоэлектроники (электронный журнал), http://ire.cplire.ru. 2000. №8. 1-33.

15. Луи de-Бройль. Электромагнитые волны в волноводах и полых резонаторах. М: иностр. лит., 1948.

16. Гуревич А.Г. Полые резонаторы и волноводы. М.: Сов. радио, 1952. 2А.Вайнштейн Л.А. Электромагнитные волны. М.: Сов. радио, 1957.

17. Смайт B.C. Электростатика и электродинамика. М.: иностр. лит., 1954.

18. Фел С.С., Фридберг П.Ш., Левинсон И.Б. Радиофизика и электроника. 1962. 6, № 11. 1125-1130.21 .Ильинский A.C., Слепян Г.Я. Колебания и волны в электродинамических системах с потерями. М.: МГУ, 1983.

19. Моденов В.П. Собственные колебания электромагнитного резонатора с импедансной сферической поверхностью. Радиотехника и электроника. 2000. 45, № 10. 1198-1201.

20. Моденов В.П., Чулков Ф.М. Электродинамика и техника СВЧ и КВЧ. 2000. 8, № 3-4. 89.

21. Моденов В.П., Гущин P.A., Ерохин А.И., Шапкина Н.Е. Математическое моделирование сферически-слоистого диэлектрического резонатора. Журнал радиоэлектроники (электронный журнал), http://jre.cplire.ru. 2006. № 4.

22. Моденов В.П., Гущин P.A., Ерохин A.M., Шапкина Н.Е. Исследование добротности многослойного сферического диэлектрического резонатора. Вестн. Моск. Ун-та. Сер. 3. Физика. Астрономия. 2007. № 2. 36.

23. Ъ2.Кравченко В.Ф. Методы определения скорости света, основанные на импедансных измерениях сверхпроводников. Радиотехника. 1995. № 10. 108-117.

24. ЪЪ.Арзин А.П., Винтизенко И.И., Канавец В.И., Сулакшин A.C., Штейн Ю.Г. Экспериментальное исследование взаимодействия сильноточного релятивистского пучка с сверхразмерной замедляющей структурой. Письма в ЖТФ. 1986. 12, № 16. 970-973.

25. ЪА. Афонин Д.Г., Малышкин А.К. Электродинамическая система с периодической структурой. 12 Международная конференция "СВЧ техника и телекоммуникационные технологии". Сентябрь 2002. 379380.

26. Lindman K.F. Annalen der Physik. 1920. 63, № 4. 621-644.

27. Кацеленбаум Б.З., Коршунова Е.Н., Сивов А.Н., Шатров АД. Киральные электродинамические объекты. Успехи физических наук. 1997. 167, № 11. 1201-1212.

28. Ъ9.Моденов В.П., Ромашин А.В. Задача дифракции электромагнитных волн на биизотропном включении в цилиндрическом волноводе. Электромагнитные волны и электронные системы. 2005. 10, № 8. 2328.

29. V.K. Varadan, V.V. Varadan, A. Lakhtakia. On the possibility of designing anti-reflection coatings using chiral composites. J. Wave-Material Interaction. 1987. 2, № 1. 71-81.

30. P. Pelet and N. Engheta. The theory of chirovaveguides. IEEE Trans. Antennas Propogat. 1990. 38. 90-98.

31. Al.S.A. Schelkunoff. Conversion of Maxwell's equations into generalized telegraphist's equation. Bell System Technical Journal. 1955. 34, № 5. 995-1043.

32. АЪ.Кацелененбаум Б.З. Теория нерегулярных волноводов с медленно меняющимися параметрами. М : Наука, 1961.

33. Краснушкин П.Е. Метод нормальных волн в применении к проблеме дальних радиосвязей. М: Изд-во МГУ, 1947.5Q.Ильинский А.С., Кравцов В.В., Свешников А.Г. Математический модели электродинамики. М: Высшая школа, 1991.

34. Свешников А.Г., Могилевский И.Е. Математические задачи теории дифракции. М: Физический факультет МГУ, 2010.

35. Ы.Завадский В.Ю. Вычисление волновых полей в открытых областях и волноводах. М: Наука, 1972.5Ъ.Григорьев А.Д., Янкевич В.Б. Резонаторы и резонаторные замедляющие системы СВЧ. Численные методы расчета и проектирования. М: Радио и связь, 1984.

36. Калиткин Н.Н. Численные методы. М: Наука, 1978.

37. Марчук Г.И. Методы вычислительной математики. М: Наука, 1980.

38. Быков A.A., Ильинский A.C. Прямой численный метод исследования электродинамических свойств полого диэлектрического трансформатора. Радиотехника и электроника. 1982. 27, № 9. 17061710.

39. Самарский A.A. Теория разностных схем. М.: Наука, 1983.61 .Свешников А.Г., Боголюбов А.Н., Митина И.В. Расчет газово-диэлектрического световода конечно-разностным методом. Радиотехника и электроника. 1982. 27, № 3. 401-408.

40. Боголюбов А.Н., Едакина Т.В. Применение вариационно-разностных методов для расчета диэлектрических волноводов. Вестн. Моск. Унта. Сер. 3. Физика. Астрономия. 1991. 32, № 2. 6-14.

41. Боголюбов А.Н., Едакина Т.В. Расчет диэлектрических волноводов со сложной формой поперечного сечения вариационно-разностным методом. Вестн. Моск. Ун-та. Сер. 3. Физика. Астрономия. 1993. 34, № 3. 72-74.

42. Свешников А.Г., Боголюбов А.Н., Минаев Д.В., Сычкова A.B. Расчет диэлектрических волноведущих систем конечно-разностным методом. Радиотехника и электроника. 1993. 38, № 5. 804-809.

43. Боголюбов А.Н., Делицын A.JJ., Красилъникова А.В., Минаев Д.В., Свешников А.Г. Математическое моделирование волноведущих систем на основе метода конечных разностей. Зарубежная радиоэлектроника. Успехи современной радиоэлектроники. 1998. № 5. 39-54.

44. Полянский Э.А. Метод коррекции решения параболических уравнений в неоднородном волноводе. М: Наука, 1985.

45. Боголюбов А.Н., Буткарев И.А. Применение метода конечных элементов к исследованию волноводного перехода. Вестн. Моск. Унта. Сер. 3. Физика. Астрономия. 2003. № 4. 6-9.

46. Стренг Г., Фикс Дж. Теория метода конечных элементов. М: Мир, 1977.

47. ЪЪ.Норри Д., de Фриз Ж. Введение в метод конечных элементов. М: Мир, 1981.

48. Марчук Г.И., Агошков В.И. Введение в проекционно-сеточные методы. М: Наука, 1981.

49. J.R. Bunch, D.J. Rose. Sparse matrix computations. New York: Academic Press, 1976.

50. J.R .Bunch, L. Kaufman. Some stable methods for calculating inertia and solving symmetric linear systems. Mathematics of Computation. 1977. 31. 163-179.

51. Райе Дж. Матричные вычисления и математическое обеспечение. М: Мир, 1984.

52. Ш.Кравченко В. Ф. Лекции по теории атомарных функций и некоторым их приложениям. М: Радиотехника, 2003.

53. Ю.Кравченко В. Ф., Рвачев В. Л. Алгебра логики, атомарные функции и вейвлеты в физических приложениях. М: Физматлит, 2006.

54. Кравченко В.Ф., Лабунько О.С., Лерер A.M., Синявский Г.П. Вычислительные методы в современной радиофизике. М: Физматлит, 2009.

55. Г. Angkaew, М. Matsuhara, and N. Kumagai. Finite-element analysis of waveguide modes: A novel approach that eliminates spurious modes. IEEE Trans. Microwave Theory Tech. 1987. MTT-35. 117-123.

56. Боголюбов A.H., Делш(ын А.Л. Расчет диэлектрических волноводов методом конечных элементов, исключающий появление нефизических решений. Вестн. Моск. Ун-та. Сер. 3. Физика. Астрономия. 1996. 37, № 1. 9-13.

57. Rahman В.М.А., Davies J.B. Penalty function improvement of waveguide solution of finite element. IEEE Trans. Microwave Theory Tech. 1984. MTT-32, №8. 922-928.

58. Bermudes A., Pedreira D.G. A finite element method for computation of waveguides. Numer. Math. 1992. 61, № 2. 39-57.

59. Nedelec J.C. Mixed finite elements in R3. Numer. Math. 1980. 35. 315341.

60. Raviart P.A., Thomas J.M. A mixed finite element method for 2-nd order elliptic problems. In. A.Dold, B. Eckman ets. Mathematical Aspects of121the Finite Element Methods. Lecture Notes in Math. (Springer Berlin Heidelberg). 1983. 606. 292-315.

61. Kikuchi F. Mixed and penalty formulations for finite element analysis of an eigenvalue problem in electromagnetism. Computer Methods in Applied Mechanics and Engineering. 1987. 64. 509-521.

62. Боголюбов А.Н., Свешников А.Г., Лопушенко В.В. Расчет градиентных оптических волокон конечно-разностным методом с использованием эффективных граничных условий. Вестн. Моск. Унта. Сер. 3. Физика. Астрономия. 1991. 32, № 2. 6-14.

63. Свешников А.Г. Принципы излучения. ДАН СССР. 1950. 3, № 5. 517-520.

64. Johnson S.G. Notes on perfectly matched layers (PMLs). math.mit.edu (Массачусетский технологический институт, отделение математики), 2010.

65. Винер Н., Пэли Р. Преобразование Фурье в комплексной области. М: Наука, 1964.

66. Рихтмайер Р. Принципы современной математической физики. М: Мир, 1982.

67. Колмогоров А.Н., Фомин С.В. Элементы теории функций и функционального анализа. М: Наука, 1989.

68. Морс Ф.М., Фешбах Г. Методы теоретической физики. Т. 1. М: иностр. лит., 1958.

69. Владимиров B.C. Уравнения математической физики. М: Наука, 1981.

70. Свешников А.Г., Боголюбов А.Н., Кравцов В.В. Лекции по математической физике. М: МГУ, 1993.

71. ИосидаК. Функциональный анализ. М: Мир, 1967.

72. Свешников А.Г., Тихонов А.Н. Теория функций комплексной переменной. М: Наука, 1967.

73. Лаврентьев М.А., Шабат Б.В. Методы теории функций комплексного переменного. М: Наука, 1987.

74. Ладыженская О.А. Краевые задачи математической физики. М: Наука, 1973.

75. Roberts J.E., Thomas J.-М. Mixed and Hybrid Methods. In Handbook of Numerical Analysis. V.2. ed. Ciarlet P.G., Lions J.L., North-Holland, 1991. 523-639.

76. DeclouxJ., NassifN., RappazJ. On spectral approximation. Part 1. The problem of convergence. RAIRO Numerical Analysis. 1978. 12, №2. 97-112.

77. Делицын А.Л. О проблеме применения метода конечных элементов к задаче вычисления мод диэлектрических волноводов. Журнал выч. мат. и мат. физ. 1999. 39, № 2. 315-322.

78. Bramble J., Osborn J. Rate of Convergence Estimates for Nonselfajoint Eigenvalue Approximation. Mathematics of Computation. 1973. 27, № 123. 525-549.

79. Банков C.E., Куру шин А. А., Разевиг В Д. Анализ и оптимизация трехмерных СВЧ-структур с помощью HFSS. М: COJIOH-Пресс, 2005.