Методы Галеркина и коллокации для решения объемного сингулярного интегро-дифференциального уравнения в задачах дифракции на диэлектрических телах тема диссертации и автореферата по ВАК РФ 01.01.07, кандидат физико-математических наук Миронов, Денис Алексеевич

• Миронов, Денис Алексеевич
• кандидат физико-математических наук
• 2010

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

Введение.

1 Дифференциальная формулировка задач дифракции.

1.1 Постановка задачи дифракции на диэлектрическом теле в свободном пространстве.

1.2 Постановка задачи дифракции на диэлектрическом теле в прямоугольном волноводе.

2 Тензорные функции Грина.

2.1 Тензорная функция Грина для свободного пространства.

2.2 Тензорная функция Грина прямоугольного волновода.

3 Метод интегро-дифференциальных уравнений для решения задач дифракции на диэлектрическом теле.

3.1 Интегро-дифференциальное уравнение электрического поля.

3.2 Объемное сингулярное интегральное уравнение.

3.3 Метод псевдодифференциальных операторов для исследования объемного сингулярного интегрального уравнения электрического поля.

3.3.1 Задача дифракции и уравнение электрического поля.

3.3.2 Теорема о композиции, эллиптичность и фредгольмовость.

3.3.3 Уравнение электрического поля как псевдодифференциальное уравнение.

3.3.4 Эллиптичность и фредгольмовость ПДО задачи.

3.3.5 О компактности оператора Я.

4 Численные методы решения интегро-дифференциальных уравнений. 47 4.1 Проекционные методы.

4.2 Метод Галеркина для решения интегро-дифференциального уравнения задачи дифракции на диэлектрическом теле в свободном пространстве.

4.2.1 Метод Галеркина.

4.2.2 Метод Галеркина для решения интегро-дифференциального уравнения электрического поля.

4.2.3 Формирование матрицы коэффициентов в методе Галеркина.

4.3 Метод коллокации для решения интегро-дифференциального уравнения задачи дифракции на диэлектрическом теле в прямоугольном волноводе.

4.3.1 Метод коллокации.

4.3.2 Метод коллокации для решения интегро-дифференциального уравнения электрического поля.

4.3.3 Формирование матрицы коэффициентов в методе коллокации.

4.4 Учет симметрии матриц.

4.4.1 Симметрия матрицы, получаемой при решении задачи с использованием метода Галеркина.

4.4.2 Симметрия матрицы, получаемой при решении задачи с использованием метода коллокации.

4.5 Параллельный алгоритм формирования матрицы.

5 Результаты расчетов.

5.1 Результаты по вычислению коэффициентов матрицы, получаемой при использовании метода Галеркина.

5.2 Результаты по вычислению коэффициентов матрицы, получаемой при использовании метода коллокации.

5.2.1 Учет симметрии матрицы при использовании метода коллокации

5.3 Решение системы линейных алгебраических уравнений.

5.3.1 Результаты по задаче дифракции на диэлектрическом теле в свободном пространстве.

5.3.2 Результаты по задачи дифракции на диэлектрическом теле в прямоугольном волноводе.

5.4 Замечания по отладке программ при использовании многопроцессорных вычислительных комплексов.

Введение диссертации (часть автореферата) на тему «Методы Галеркина и коллокации для решения объемного сингулярного интегро-дифференциального уравнения в задачах дифракции на диэлектрических телах»

Задачи дифракции электромагнитных волн на рассеивателях с конечной проводимостью (на диэлектрических телах) имеют широкий диапазон применений в микроволновой технике, медицине, физике плазмы и в других областях. Задача дифракции на неоднородном теле произвольной формы имеет тесную связь с вопросами о влиянии электромагнитного излучения на биологические объекты.

Исследование задач дифракции электромагнитных волн на диэлектрических телах является классической проблемой электродинамики. Традиционная (физическая) теория дифракции создавалась на протяжении нескольких столетий X. Гюйгенсом, О. Френелем, Г. Гельмгольцем, Г.Р. Кирхгофом, Д. Лармором и другими авторами. Для понимания волновых процессов и расчета дифракционных полей большое значение имеет принцип Гюйгенса, согласно которому распространение волн обусловлено действием вторичных источников. Френель уточнил принцип Гюйгенса, приняв во внимание интерференцию сферических волн, излучаемых вторичными источниками. Дальнейшее уточнение принципа Гюйгенса - Френеля принадлежит Кирхгофу, который дал его строгую формулировку, основываясь на уравнении Гельмгольца. В современной теоретической оптике приближенное решение дифракционных задач производится почти исключительно с помощью принципа Гюйгенса - Кирхгофа. Электродинамическая (векторная) формулировка принципа Гюйгенса была дана Котлером.

Благодаря работам А. Пуанкаре и А. Зоммерфельда стало ясно, что в задачах дифракции электромагнитных волн речь идет о некоторой краевой задаче математической физики. В общей постановке задача состоит в нахождении решений уравнений Максвелла, удовлетворяющих определенным краевым условиям. К этому надо добавить «условия излучения» (Зоммерфельд, 1912), состоящие в том. что вся энергия, излучаемая источником, должна уходить в бесконечность. Кроме того, следует учитывать особое поведение полей в окрестности края поверхности гонкого экрана. Первое аналитическое решение задачи дифракции на идеально проводящей полуплоскости было дано Зоммерфельдом [1,2]. Уже это решение позволило сделать ряд важных выводов о поведении электромагнитного поля в ближней и дальней зоне, об особенности полей в окрестности края тонкого экрана, о поведении полей на бесконечности и т.д.

Релеем [3] рассмотрена двумерная задача возбуждения диэлектрического цилиндра кругового сечения плоской волной. Решение находится методом разделения переменных и представляет собой ряд по тригонометрическим и цилиндрическим функциям (ряд Рэлея), который сходится при всех значениях к0К, где = 2л IЛ - волновое число, Я- длина падающей волны, Я- радиус цилиндра. Если к0Я<< 1, то ряд сходится достаточно быстро. При к0Я » 1 сходимость ряда может быть улучшена при помощи преобразования Ватсона [4]. В резонансном случае (т.е. для цилиндров с диаметром порядка длины волны) ряд Рэлея сходится плохо, и для его суммирования необходимо применять специальные методы [5]. Методом разделения переменных получено решение для диэлектрического шара [6] и для некоторых других тел простой формы.

Решение задач дифракции на телах сложной формы невозможно аналитически, поэтому необходимо разрабатывать численные (приближенные) методы решения таких задач. Для тел, размеры которых значительно меньше или значительно больше длины волны, возможно применение асимптотических методов [7]. Однако в «резонансном» случае, когда размеры тела сравнимы с длиной волны, применение асимптотических методов невозможно, и приходится решать задачу численно с помощью современных компьютеров.

При этом необходимо решать трехмерные векторные краевые задачи для системы уравнений Максвелла в полной электродинамической постановке. Решение таких задач является в настоящее время одной из самых актуальных проблем в электродинамике и с приемлемой для практики точностью на электродинамическом уровне строгости математическими методами требует очень большого объема вычислении, что зачастую невозможно даже на самых современных суперкомпьютерах. Многочисленные дорогостоящие пакеты прикладных программ для решения задач электродинамики (Апз!б, Ошкшауе и т.д.), имеющиеся на рынке программных продуктов, не используют кластерных технологий, решают задачу традиционными конечно-разностными методами или методами конечных элементов и не дают удовлетворительных по точности результатов. Несмотря на достаточно высокую цену предлагаемых продуктов, результаты, анонсируемые их производителями, сильно разнятся и не могут считаться надежными [8]. Последнее объясняется, прежде всего, отсутствием интереса у разработчиков ПО к данному кругу задач. Кроме того, указанные пакеты являются лишь модификациями программ для решения задач, связанных с техникой связи, и не имеют хорошего теоретического обоснования.

При решении краевых задач в неограниченных областях конечно-разностные методы и методы конечных элементов встречают принципиальные трудности: область, в которой решается задача, должна быть сделана конечной. Такая редукция области приводит к появлению неконтролируемой погрешности, причем размеры области для уменьшения погрешности должны быть достаточно велики. Конечно-разностные методы и методы конечных элементов в такой ситуации обычно приводят к очень большим, но разреженным матрицам (порядка 109 и более).

В качестве альтернативного метода решения задач в неограниченных областях применяется метод интегральных или интегродифференциальных уравнений. В этом случае задача сводится к интегральному или интегродифференциальному уравнению в области неоднородности, которая по размерам существенно (на порядки) меньше области решения задачи в случае применения конечно-разностных методов и методов конечных элементов. После дискретизации получается конечномерная система уравнений с матрицей существенно меньшего порядка (103 -104). До недавнего времени для задач дифракции резонансного диапазона не существовало универсального метода решения всех перечисленных выше задач, приемлемого с точки зрения вычислительных ресурсов, и для решения более общих задач приходилось вносить в их постановку различные упрощения, например, считать диэлектрическое тело «телом вращения» или заменять непрерывное распределение параметров среды кусочно-постоянным. Впервые объемные интегральные уравнения были выведены в [9]. В работах [10, 11] были приведены объемные сингулярные интегральные уравнения относительно вектора электрического поля, описывающие трехмерные задачи дифракции на неоднородных диэлектрических рассеивателях. К достоинствам этих уравнений следует отнести простоту и универсальность учета неоднородности любого типа, т.е. анизотропии. Опишем кратко результаты, полученные в этих работах.

-г

Рассмотрен следующий класс задач. Пусть область Qc: К , ограниченная поверхностью 5, характеризуется произвольной тензор-функцией ¿(х), обращающейся вне <2 в константу - диэлектрическую проницаемость £0 свободного пространства. Магнитная проницаемость во всем пространстве предполагается постоянной. Требуется определить электромагнитное поле, возбуждаемое в данной среде внешним монохроматическим полем, представляющим собой либо плоскую волну, либо поле стороннего тока. Вводя электрический ток поляризации, применяя известные формулы векторного потенциала и теорему о дифференцировании интегралов со слабой особенностью [12], для рассматриваемой задачи получено [13] следующее векторное сингулярное интегральное уравнение относительно электрического поля:

E(*) = E°(*)-i j

Q\ £0 J g(*) j

Е(*) +

0 J

Е (y)G(r)dy+ (0.1) v

•P-i

E (y),grad X J gradxG(r)dy,

G(r) - exP(;/H)r) фундаментальное где и{г)-——-— - фундаментальное решение трехмерного уравнения 4ят

Гельмгольца в свободном пространстве, соответствующее временной зависимости ехр(—¿Ш), а Е°(.т) - известная гладкая функция. При достаточно общих предположениях о компонентах тензора проницаемости, оператор уравнения является ограниченным в пространстве и в П^] получена оценка его нормы.

Уравнение (0.1) относится к классу многомерных сингулярных интегральных уравнений, общая теория которых построена в [12]. Основные результаты этой теории получены для уравнений на многообразиях без края, поэтому для того, чтобы применить их к (0.1), необходимо предварительно определить продолжение оператора уравнения в ). В [13] такое продолжение построено, исследована его взаимосвязь с исходным оператором, и в явном виде выписан символ продолженного оператора (0.1).

Утверждение 0.1. Пусть компоненты тензора диэлектрической проницаемости непрерывны во всем пространстве. Тогда для того, чтобы оператор уравнения был нетеров в необходимо и достаточно выполнения условия з п!" X аг (в, ф)а/(в, ф)Еу {х) > 0,

0.2) где а}{в,ф) - декартовы координаты точек единичной сферы, е^ — компоненты тензора ё в декартовой системе координат.

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

Утверждение 0.2. Пусть декартовые составляющие исходного поля и компоненты тензора диэлектрической проницаемости непрерывны по Гельдеру во всем пространстве, выполняется условие (0.2), а тензор системы (0.1) существует, единственно в ¿2(<2) и является классическим решением задачи дифракции, т.е. удовлетворяет поточечно уравнениям Максвелла и условиям излучения.

Для изотропной среды требование непрерывности диэлектрической проницаемости во всем пространстве удается снять. В [14] для системы (0.1) в явном виде построен регуляризирующий оператор, и теорема существования и единственности решения задачи дифракции доказана при выполнении неравенства л е(х)-£ (х) 21 положительно определен почти всюду в О. Тогда решение

1тг:(л:) > 1т .

0.3)

В [13] также проведено исследование системы (0.1) при минимальном требовании ограниченности компонент тензора диэлектрической проницаемости. При условии

Нх)-ё\х)-2/W

2 i обобщающем (0.3) на анизотропный случай, доказана теорема существования и единственности решения (0.1) в L2(Q).

Подробное исследование задачи дифракции в свободном пространстве на сложных диэлектриках с применением метода интегральных уравнений проведено A.C. Ильинским, А.Б. Самохиным и Ю.Ю. Капустиным в работах [15-17].

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

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

В работе Ю.А. Еремина, В.И. Ивахненко [19] рассматривается проблема моделирования царапины на стенке волновода. Исходная дифференциальная задача сводится с использованием тензорной функции Грина к интегральному уравнению, для приближенного решения которого построен и обоснован численный метод.

В работе Е.М. Карчевского [20] метод сингулярных интегральных уравнений успешно применен при исследовании проблемы о собственных модах волновода с размытой границей.

До недавнего времени главной областью применения волноводов и резонаторов была техника СВЧ. К основным вопросам, которые рассматривались исследователями, можно отнести следующие: определение собственных частот и волн для волноводов и резонаторов разной формы; исследование характера волновых процессов в волноводах и резонаторах; изучение влияния на поля потерь в волноводах и резонаторах и определение добротности; исследование диэлектрических, оптических и квазиоптических волноводов и резонаторов. Эти и ряд других вопросов исследованы в монографиях [21-23]. Наиболее полно математический аппарат для решения задач электродинамики применительно к волноводам и резонаторам исследован в работах В.В. Никольского [24 - 27].

В данной работе рассматриваются две задачи дифракции на диэлектрическом теле в свободном пространстве и в бесконечном прямоугольном волноводе. Несмотря на различные краевые условия, получаем однотипные интегро-дифференциальные уравнения, используя функции (тензоры) Грина.

Для численного решения задачи дифракции в свободном пространстве ' используется метод Галеркина. Выбор метода Галеркина обусловлен тем, что этот метод обладает большой устойчивостью к погрешностям вычислений. При его использовании получается плотная матрица с блочно-теплицевой структурой. При заполнении получаемой матрицы используются средства мощных вычислительных комплексов - кластеров (суперкомпьютеров), так как процесс заполнения при использовании ПЭВМ занимает неприемлемо большое время. Для запуска программ на кластерах широко используется реализация стандарта MPI [28] - OpenMPICH [29]. Показаны проблемы, встречающиеся при отладке и запуске программы, реализованной под

OpenMPICH и пути их решения. Показаны методы распределения данных для одновременного использования нескольких узлов. Получены первые результаты. Также для запуска программы заполнения матрицы была использована система «х-сош» [30] позволившая заполнить матрицу для исследований еще большего размера, чем реализация MPI.

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

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

Во второй главе описаны функции Грина, использованные для перехода к интегро-дифференциальным и интегральным уравнениям.

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

В четвертой главе описаны численные методы решения задач. В первой части исследуются проекционные методы, показаны основные утверждения о сходимости методов. Во второй части исследуется метод Галеркина, использованного для случая свободного пространства, в третьей части — метод коллокации для случая прямоугольного волновода. Доказываются теоремы о сходимости методов для решения соответствующих интегро-дифференциальных уравнений. В четвертой части описываются методы оптимизации на основе свойств симметрии получаемых матриц в использованных численных методах для сокращения времени вычислений и сокращения выделяемой памяти. В пятой части, в терминах стандарта MPI, описываются методы распределения вычислений для обоих случаев, а также, в терминах системы «х-сош», метод раздачи заданий узлам (только случай свободного пространства) для одновременного использования произвольного числа узлов суперкомпьютерного комплекса.

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

Заключение диссертации по теме «Вычислительная математика», Миронов, Денис Алексеевич

Результаты работы были использованы для решения обратной задачи дифракции на диэлектрическом теле, размещенном в прямоугольном волноводе - задачи определения диэлектрической проницаемости материалов [65], [66].

Заключение

Данная работа содержит следующие основные результаты:

1. Краевые задачи дифракции на диэлектрическом теле в свободном пространстве и в прямоугольном волноводе сведены к решению векторных сингулярных интегро-дифференциальных уравнений.

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

3. Разработан метод Галеркина для численного решения векторного сингулярного интегро-дифференциального уравнения для задачи дифракции на диэлектрическом теле в свободном пространстве. Доказана теорема о сходимости метода Галеркина.

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

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

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

1. Зоммерфельд А. Дифференциальные уравнения в частных производных в физике. М.: Иностр. литература, 1950.

2. Зоммерфельд А. Оптика. М.: Иностр. Литература, 1973.

3. Lord Rayleigh // Ph.Mag., Sc. Papers. 1918. V.36.

4. Хенл, Мауэ, Вестпфаль. Теория дифракции. М.: Мир, 1964.

5. Воронцов А.А., Мировицкая С.Д. Специальные функции теории задачи рассеяния. ~М.: Радио и связь, 1991.

6. Ван де Хюлст Г. Рассеяние света малыми частицами. М.: И.Л., 1961.

7. Кравцов Ю.А., Орлов Ю.И. Геометрическая оптика неоднородных сред. -М.: Наука, 1980.

8. Yakovlev V.V. Commercial ЕМ codes suitable for modelling of microwave Heating a comparative review // Proceedings of the 3rd international workshop. - Berlin: Springer, 2001. P. 87-95.

9. Хижняк H.A. Функция Грина уравнений Максвелла для неоднородных сред // Ж. тех. физ. 1958. Т.28, №7. С. 1592-1609.

10. Livesay Р.Е., Chen К. Electromagnetic fields included indide arbitrary biological bodies // IEEE Trans., 1974. V.MTT-22, №12.

11. Yaghjian A.D. Electric dyadic Green's functions in the source region // Proc. IEEE. 1980. V.68, №2,

12. Михлин С.Г. Многомерные сингулярные интегралы и интегральные уравнения. М.: Физматгиз, 1962.

13. Самохин А.Б. Исследование задач дифракции электромагнитных волн в локально-неоднородных средах // ЖВМиМФ. 1990. Т.ЗО, №1.

14. Самохин А.Б. Дифракция электромагнитных волн на локально-неоднородном теле и сингулярные интегральные уравнения // ЖВМиМФ. 1992. Т.32, №5.

15. Ильинский А.С., Капустин Ю.Ю., Самохин А.Б. Математическая модель задачи дифракции на неоднородном цилиндрическом теле // ЖВМиМФ. 1998. Т.28, №9.

16. Ильинский А.С., Капустин Ю.Ю., Самохин А.Б. Метод сингулярного интегрального уравнения для решения задачи дифракции на неоднородном теле // Межв. сб. «Дифракция и распространение радиоволн». М.: МФТИ, 1998.

17. Савостьянов Д.В., Тыртышников E.E. Применение многоуровневых матриц специального вида для решения прямых и обратных задач электродинамики // Вычислительные методы и программирование: новые вычислительные технологии. 2006. Т. 7, №1. С. 1-16.

18. Еремин Ю.А., Ивахненко В.И. Строгие и приближенные модели царапины на основе метода интегральных уравнений // Дифф. уравнения. 2001. Т.37, №10. С. 1386-1394.

19. Гуревич А.Г. Полые резонаторы и волноводы. М.: Мир, 1974.

20. Вайнштейн Л.А. Открытые резонаторы и открытые волноводы. М.: Сов. радио, 1966.

21. Григорьев А.Д., Янкевич В.Б. Резонаторы и резонаторные замедляющие системы СВЧ. М.: Радио и связь, 1984.

22. Никольский B.B. Вариационные методы для внутренних задач электродинамики. М.: Наука, 1967.

23. Никольский В.В. Математический аппарат электродинамики. М.: изд-воМИРЭА, 1973.

24. Никольский В.В., Никольская Т.И. Декомпозиционный подход к задачам электродинамики. М.: Наука, 1983.

25. Никольский В.В. Никольская Т.И. Электродинамика и распространение радиоволн. М.: Наука, 1989.

26. MPI: А Message Passing Interface Standart. Version 1.0. - University of Tennessee, May 5, 1994.

27. W.Gropp, E.Lusk. Technical Report ANL-96/5, Argonne National Laboratory, 1996.

28. Филамофитский М.П. Система поддержки метакомпьютерных расчетов x-com: Архитектура и технология работы // Вычислительные методы и программирование. НИВЦ МГУ им. М.В. Ломоносова, 2004. Т.5. С. 1-9.

29. Ильинский A.C., Смирнов Ю.Г. Дифракция электромагнитных волн на проводящих тонких экранах. М.: ИПРЖР, 1996.

30. Adams R. Sobolev spaces. New York.: Academic press, 1975.

31. Самохин А.Б. Интегральные уравнения и итерационные методы в электромагнитном рассеянии. М.: Радио и Связь, 1998.

32. Цупак A.A. Метод Галеркина для решения интегрального уравнения в задаче дифракции на локально неоднородном теле в случае Н-поляризации // Труды математического центра имени Н.И.

33. Лобачевского. Казань: НИИММ им. Чеботарева, 2000. Т. 6. С. 240248.

34. Tsupak А.А. Vector integral equation method for diffraction problem in a cavity resonator // Processing and abstracts of 2001 Far-Eastern school-seminar on mathematical modeling and Numerical Analysis. Nakhodka, Russia: August 22-28, 2001. P. 200

35. Михлин С.Г. Сингулярные интегральные уравнения // Успехи математических наук. 1948. Т. 3, №3. С. 29-112.

36. Колтон Д., Кресс Р. Методы интегральных уравнений в теории рассеяния. М.: Мир, 1987.

37. Канторович Л.В., Акилов Г.П. Функциональный анализ. М.: Наука, 1984.

38. Егоров Ю. В., Шубин М. А. Линейные дифференциальные уравнения с частными производными. Элементы современной теории. / Итоги науки и техники. Сер. Современные проблемы математики. Фундаментальные направления. -М.: ВИНИТИ, 1988. Т. 31. С. 5-125.

39. Шубин М. А. Псевдодифференциальные операторы и спектральная теория. М.: Наука, 1978 (переиздание: М.: Добросвет, 2005).

40. Ремпель LLL, Шульце Б.-В. Введение в общую теорию индекса эллиптических краевых задач. М.: Мир, 1986.

41. Като, Т. Теория возмущений линейных операторов М.: Мир, 1972.

42. Егоров, Ю. В. Лекции по уравнениям с частными производными. Дополнительные главы — М.: Изд-во Московского университета, 1985.

43. Панич, О. И. Введение в общую теорию эллиптических краевых задач -Киев: ВищаШкола, 1986.

44. Мазья, В. Г. Пространства Соболева Л.: Изд-во Ленинградского университета, 1985.

45. Kress R. Linear Integral Equations. Applied Mathematical Sciences. V.82. -Springer-Verlag. New-York Inc., 1989

46. Вайнштейн Л.А. Электромагнитные волны M:. Радио и связь. 1988.

47. Воеводин В.В., Воеводин Вл. В. Параллельные вычисления. — СПб.: БХВ-Петербург, 2002.

48. Жуматий С.А. Система анализа производительности параллельных программ на кластерных установках // Вычислительные методы и программирование. НИВЦ МГУ им. М.В. Ломоносова, 2005. Т.6 С. 57-64.

49. Шилдт Г. Полный справочник по С++. 4-е изд. М.: Вильяме, 2003.

50. Дацюк В.Н., Букатов A.A., Жегуло А.И. Методическое пособие по курсу "Многопроцессорные системы и параллельное программирование" Введение в организацию и методы программирования многопроцессорных вычислительных систем. -Ростов-на-Дону, 2000.

51. Медведик М.Ю., Смирнов Ю.Г. Численное решение объемного сингулярного интегрального уравнения методом коллокации. // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. 2009. № 4. С. 55-71.

52. Миронов Д.А. Математическое моделирование с использованием мощных вычислительных комплексов. // Надежность и качество. Труды международного симпозиума. Россия, Пенза, 23-31 мая, 2005 г. -Пенза: ПГУ, 2005. С. 141-145.

53. Немет Э., Снайдер Г., Хейн Т. Руководство администратора Linux. -М.: Вильяме, 2004.

54. Немнюгин С., Чаунин М. Эффективная работа: Unix. СПб.: Питер, 2001.

55. Колисниченко Д.Н. Linux-сервер своими руками. 3-е изд. СПб.: Наука и Техника, 2005.

56. Smirnov Y., Shestopalov Y., Mironov D. Analysis of Inverse Scattering in a Waveguide using the Method of Volume Singular Integral Equation // URSI International Symposium on Electromagnetic Theory (EMTS 2010). Berlin: August 16-19, 2010. P. 532-534.

57. Миронов Д.А., Смирнов Ю.Г. О существовании и единственности решений обратной краевой задачи определения диэлектрической проницаемости материалов. // ЖВМиМФ. 2010. Т.50, №9. С. 1587-1597.