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

• Головашкин, Димитрий Львович
• кандидат физико-математических наук
• 2000

Оглавление диссертации кандидат физико-математических наук Головашкин, Димитрий Львович

ВВЕДЕНИЕ.

1. РАЗНОСТНЫЕ СХЕМЫ ДЛЯ УРАВНЕНИЙ МАКСВЕЛЛА.

1.1. Постановка смешанной задачи Дирихле и Неймана для уравнений Максвелла.

1.2. Явно-неявная разностная схема, записанная в трехмерной декартовой системе координат.

1.3. Явно-неявные разностные схемы, записанные в двумерной декартовой системе координат

1.4. Неявные разностные схемы, записанные в двумерной декартовой системе координат.

Выводы.

2. ИССЛЕДОВАНИЕ ПРОХОЖДЕНИЯ ЭЛЕКТРОМАГНИТНОГО ИЗЛУЧЕНИЯ ЧЕРЕЗ ЦИЛИНДРИЧЕСКИЕ ДИФРАКЦИОННЫЕ КОРОТКОФОКУСНЫЕ МИКРОЛИНЗЫ.

2.1. Постановка задачи.

2.2. Анализ распределения интенсивности оптического излучения на главной оптической оси рефракционной цилиндрической микролинзы и дифракционных цилиндрических микролинз.

Выводы.

3. РАСЧЕТ ЭНЕРГИИ, ОТРАЖЕННОЙ ОТ АНТИОТРАЖАЮЩЕЙ БИНАРНОЙ СТРУКТУРЫ.

3.1. Применение разностного решения уравнений Максвелла для расчета отраженной энергии.

3.2. Применение теории эффективных сред для расчета отраженной энергии.

Выводы.

4. РЕАЛИЗАЦИЯ ПАРАЛЛЕЛЬНЫХ ВЫЧИСЛЕНИЙ ПРИ РАЗНОСТНОМ РЕШЕНИИ УРАВНЕНИЙ МАКСВЕЛЛА.

4.1. Выбор явно-неявной схемы при решении уравнений Максвелла, записанных в двумерной декартовой системе координат.

4.2. Реализация поперечной прогонки.

4.3. Распараллеливание алгоритма решения разностных уравнений неявных схем, записанных в двумерной декартовой системе координат.

4.4. Распараллеливание алгоритма решения разностных уравнений явно-неявной разностной схемы, записанной в трехмерной декартовой системе координат.

Выводы.

Введение диссертации (часть автореферата) на тему «Моделирование распространения электромагнитного излучения в системах, содержащих дифракционные оптические элементы, методом разностного решения уравнений Максвелла»

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

Актуальность темы.

В последние годы особый интерес в оптике вызывают задачи, в которых необходимо учитывать поляризацию электромагнитной волны, френелевское отражение и распределение интенсивности, как в самом оптическом элементе, так и в его выходной плоскости. Это связано с появлением интегральных оптических схем, антиотражающих решеток, дифракционных микролинз для ввода излучения в оптическое волокно. Методы геометрической и скалярной оптики не позволяют исследовать прохождение излучения через данные системы. Так, геометрическое приближение не позволяет в полном объеме учитывать дифракцию на дифракционной микролинзе, что делает лучевой подход ограани-ченным для исследования микролинз. В рамках скалярного приближения известна работа Н. JI. Казанского [28], в которых рассматривается задача распространения излучения через микролинзу в приближении Френеля-Кирхофа. F. Wyrowski для этих целей использовал приближение Френеля [102]. Однако для исследования короткофокусных микролинз, используемых для ввода излучения в волокно, необходимо рассчитывать распределение интенсивности на расстоянии нескольких длин волн от микролинзы. Такую задачу рассмотрели M.S. Mirotznik, D.W. Prather и J.N. Mait решая уравнение Гельмгольца методом конечных элементов [87]. Однако при этом не исследовалось прохождение излучения через рефракционную линзу и не оценивалось распределение интенсивности внутри предложенных дифракционных микролинз. Подход, предложенный в [87] не позволяет исследовать оптический элемент в волноводе. Поэтому 5 актуально исследование прохождения электромагнитной волны через рефракционную и дифракционные микролинзы, заключенные в волновод.

При изготовлении дифракционных микролинз для мощных лазеров применяют материалы, обладающие высокой устойчивостью к термическому разрушению. Наиболее прочным материалом является алмаз, однако он обладает высоким показателем преломления, что обусловливает отражение большой доли энергии от границы раздела алмаз-воздух. Для снижения доли отраженной энергии алмазные оптические элементы просветляют, формируя на поверхности антиотражающий субволновый рельеф. Принципы работы антиотражаю-щих поверхностей глубоко исследовались в работах D.H. Raguin и G.M. Morris [92-96]. Методика расчета бинарного антиотражающего рельефа, согласно с теорией эффективных сред, представлена в работе М.Е. Motamedi, W.H. Southwell и W.J. Gunning [89]. Однако аналогия с теорией тонких пленок не является очевидной и требует проверки. Поэтому актуален поиск оптимальной формы рельефа антиотражающей поверхности с помощью численного решения уравнений Максвелла.

Решение уравнений Максвелла для задач оптики требует применения численных методов. Известен аналитический подход, рассматривающий прохождение электромагнитной волны через дифракционные решетки, однако для общего случая аналитические методы отсутствуют.

В оптике широко применялись разностные схемы. Так A. Renaut составил неявную разностную схему для решения волнового уравнения [98], а Т. Fla записал консервативную разностную схему для решения нелинейного уравнения Шредингера [83]. S.T. Chu, W.P. Huang и S.K. Chaudhuri представили явную разностную схему для уравнений Максвелла [82], которая не является устойчивой. При решении задачи построения разностных схем для уравнений Максвелла вполне очевиден отказ от явных схем. Даже будучи устойчивыми, они, как правило, налагают жесткие ограничения на шаги сетки, что делает неявные схемы более предпочтительными. 6

В классических работах Самарского A.A. и Марчука Г.И. представлен большой выбор неявных разностных схем [44,45,59-62], однако для уравнений Максвелла схемы не приведены, поэтому их разработка является актуальной.

Разностное решение задач оптики требует существенных временных затрат. Следовательно, актуальна разработка алгоритмов, реализующих параллельные вычисления. Особенности организации параллельных вычислений для решения уравнений математической физики изложены И.Н Молчановым [49]. Н.Н Миренковым [47] представлен алгоритм распараллеливания продольно-поперечных прогонок, а Б.А. Головин [27] описал методику расчета характеристик параллельных вычислительных процессов. Однако перечисленные авторы не исследовали проблемы построения алгоритмов распараллеливания для решения разностных уравнений явно-неявных схем и повышения эффективности алгоритмов распараллеливания продольно-поперечных прогонок за счет асин-хронизации и использования как правой, так и левой прогонок.

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

A. Renaut в работе [98] предложил неявную разностную схему для волнового уравнения, записанного в двумерной декартовой системе координат. Схема имеет второй порядок аппроксимации по пространству и времени, требуя задания начального условия на двух временных слоях сетки. Однако невозможно развить предложенный в [98] подход для изучения волновых процессов в трехмерном случае. Решая волновое уравнение можно изучить поведение электромагнитной волны только одного типа (Н или Е). Те же ограничения верны для работы [87], в которой M.S. Mirotznik, D.W. Prather и J.N. Mait решают уравнение Гельмгольца методом конечных элементов. Т. Fla [83] использовал консервативную разностную схему для решения нелинейного уравнения Шре-дингера, которое учитывает нелинейные эффекты при прохождении излучения 7 через оптические элементы, но верно только для параксиального случая, а значит не может применяться для исследования микролинз. S.T. Chu, W.P. Huang и S.K. Chaudhuri [82] представили явную разностную схему для уравнений Максвелла, записанную на равномерной сетке. Показано, что схема имеет первый порядок аппроксимации по пространству и времени и выписаны условия устойчивости для одномерного, двумерного и трехмерного случаев. Однако в ходе вычислительных экспериментов не обнаружены параметры сетки, при которых приведенная схема являлась бы устойчивой.

В работе М.А. Голуба, H.JI. Казанского, В.А. Сойфера [28] подробно рассматривается задача распространения излучения через микролинзу в приближении Френеля-Кирхофа. Получена зависимость эффективности дифракционных микролинз от числа уровней квантования фазовой функции, представлены величины отклонений фокусов дифракционных микролинз от геометрического фокуса. M.S. Mirotznik, D.W. Prather и J.N. Mait [87] исследовали распределение интенсивности света при прохождении излучения через дифракционные микролинзы в непосредственной близи от линзы. Однако системного изучения влияния числа уровней квантования микролинз на эффективность и смещение фокуса не проводилось. Также в [87] не представлена интенсивность поля внутри дифракционных микролинз.

Теория эффективных сред применялась для определения антиотражающей способности субволнового рельефа с конца прошлого века. Тогда была сформулирована теория эффективных сред первого порядка [92]. В середине нынешнего века С. М. Рытов представил теорию эффективных сред второго порядка [99]. Для разнообразных приложений теория эффективных сред развивается в работах D.H. Raguin и G.M. Morris [92-96].

Для реализации параллельных вычислений с целью сокращения времени решения разностных уравнений, в работах Вальковского В.А. [8,9], Миренкова H.H. [8, 47] и Молчанова И.Н. [49] представлены разнообразные алгоритмы распараллеливания. Приводятся алгоритмы для явных и неявных схем, обсуж8 даются вопросы разделения памяти между процессами. Однако отсутствуют алгоритмы распараллеливания для явно-неявных схем, а приведенные алгоритмы реализации прогонки могут быть улучшены.

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

Для достижения указанной цели определены основные задачи диссертации:

1. Построение корректных и сходящихся разностных схем для решения системы уравнений Максвелла.

2. Изучение влияния числа уровней квантования дифракционной микролинзы на положение ее фокуса и эффективность.

3. Определение параметров бинарной субволновой структуры на алмазных пленках, обеспечивающей наименьшие потери на отражение инфракрасного (ИК) излучения.

4. Реализация параллельных вычислений с целью сокращения времени вычислительного эксперимента.

Структура и краткое содержание диссертации. Диссертация состоит из Введения, четырех Глав, Заключения и двух Приложений. В первой Главе предложены явно-неявные и неявные разностные схемы для системы уравнений Максвелла, записанной в трехмерной и двумерной декартовой системе координат, проводится исследование этих схем.

Заключение диссертации по теме «Применение вычислительной техники, математического моделирования и математических методов в научных исследованиях (по отраслям наук)», Головашкин, Димитрий Львович

ВЫВОДЫ

1. Расчет с помощью решения уравнений Максвелла позволяет определить положение фокусов микролинз и их эффективность.

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

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

95

3 РАСЧЕТ ЭНЕРГИИ, ОТРАЖЕННОЙ ОТ АНТИОТРАЖАЮЩЕЙ

БИНАРНОЙ СТРУКТУРЫ

3.1. Применение разностного решения уравнений Максвелла для расчета отраженной энергии

При использовании сред с большим показателем преломления заметное влияние на энергетическую эффективность оптического элемента начинает оказывать высокое отражение от границы среда - воздух. Например, от границы германий (ве) - воздух для длины волны 10,6. мкм отразится 36% энергии [38].

Традиционно [4,41,16,76] эта проблема решается с помощью наложения на границу тонкой пленки с показателем преломления пШт - »

3.1) и толщинои 11 =

3.2) где п; и п8:— соответственно показатели преломления среды из которой падает излучение, и среды, на которую падает излучение.

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

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

Другой подход к уменьшению отраженной энергии - формирование субволнового профиля на границе двух сред, уменьшающего отражение [99]. Для простоты в данной работе будет рассмотрен бинарный профиль, представленный в [89] (рис. 3.1)

96 Л

Рис. 3.1

Период бинарной антиотражающей структуры.

На рис. 3.1 период структуры обозначен Л, ширина выступа Ь, высота берется равной Ь. Ограничимся рассмотрением случая, когда на структуру падает Н волна и вектор электрического поля направлен вдоль штрихов рельефа.

В диссертации ставилась задача найти отраженную энергию от бинарной структуры при различных значениях заполнения я=Ь/Л для Л,=10,6 мкм. ♦ *

18 ,90 ] на границе раздела среда - воздух. Показатель преломления среды (поликристаллической алмазной пленки) составляет п5=2,4. [23*]. Актуальность задачи обусловлена возможностью применения поликристаллической алмазной пленки для изготовления дифракционных элементов, способных выдержать излучение СОг лазеров мощностью 10-20 кВт [11]. Традиционные материалы РЖ оптики, такие как селенид цинка или хлорид калия таких мощностей не выдерживают. $ $ $

Поставим вычислительный эксперимент [25 ,26 ,91 ], в котором волна Н01 будет падать на бинарный антиотражающий рельеф с Л

10,6 10,6 а = —из

3,4 ' 4Д4

3.2). Для постановки эксперимента использовалась схема 2 (1.81)-(1.83) со еле дующими параметрами дискретизации: Ь

7=—, Ь=—, N. — = 2000. При 2 100 у 100 Ь2 этом погрешность х<0,01. В отличии от экспериментов предыдущей главы для хорошего описания волны длиной 10,6 мкм достаточно использования более

97 грубых шагов дискретизации. Размеры исследуемой области составляли Ьу=50 мкм, Ь2=330 мкм.

Искалась усредненная по двадцати длинам волн прошедшей волны велитЦ чина 11Е^ последиэлектрика^У^ в воздухе после диэлектрика с антиотражающим о о бинарным рельефом и рассчитывалось значение доли отраженной энергии тЦ Е х после диэлектрика^У^ р = ---(3.3)

1 0 0 ТЬ /^х без диэлектрика ^У^^

0 0 тЦ где / без диэлектрика^1 " усредненный по двадцати длинам волн квадрат о о напряженности электрического поля в воздухе при отсутствии диэлектрика.

Из таблицы 3.1 видно, что менее всего энергии отражается от бинарного рельефа с заполнением q= 0,25, тогда доля отраженной энергии р=1,2471 %.

На рис. 3.2 представлен усредненный модуль амплитуды падающей и отраженной волн в эксперименте с q=0,25. Из рис. 3.2 видно, что на рельефе бинарной антиотражающей структуры возникают волны Н высоких порядков, которые порождены скачками профиля антиотражающей структуры.

99 ширина волновода (мкм) Рис. 3.2

Усредненный модуль амплитуды падающей и отраженной волн.

3.2. Применение теории эффективных сред для расчета отраженной энергии

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

Теория эффективных сред ставит в соответствие бинарному профилю тонкую диэлектрическую пленку с проницаемостью вЭфф. (рис. 3.3).

100 ж1™*"

- -""Л

-о :?т

Рис. 3.3

Сопоставление бинарному профилю диэлектрической пленки.

В теории эффективных сред нулевого порядка [93,97] £Эфф=б(0) . При ориентации вектора электрического поля вдоль штрихов рельефа е(0Че5+(1-Я)^. (3.4)

В теории эффективных сред второго порядка [94] вЭфф=е(2) . При ориентации вектора электрического поля вдоль штрихов рельефа

8(2)=£(0)(1 + Ав), (3.5)

К2 2 Л \2 (а-1)2 АБ = —q¿{l-q)

1 + я1сх -1

3.6) где а = п0

П:

При нормальном падении антиотражающий эффект должен наблюдаться при [88] X Л тах(п;,п5).

3.7)

При изготовлении антиотражающих профилей период профиля находят так [96] X А

3.8)

101

Расчет отраженной энергии происходит по известной формуле из теории тонких пленок [14]: где \У1 и лу2 - энергия, отразившаяся на границах 1 и 2 , найденная по формуле Френеля [14], р - энергия, отразившаяся от тонкопленочного покрытия, 8 - набег фазы в тонкопленочном слое толщиной с! при нормальном падении.

В таблице 3.2 представлена зависимость доли отраженной энергии, рассчитанной по формулам (3.4)-(3.10) теории эффективных сред, от значений параметра q.

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

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

Р = ехр(21б)

3.9)

1 + ехр(21б)

ЗЛО)

ЗАКЛЮЧЕНИЕ

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

Получены следующие основные результаты:

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

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

3. На основе разработанных разностных схем проведено моделирование прохождения электромагнитного излучения через рефракционную и дифракционные микролинзы и субволновые антиотражающие структуры. Моделирование позволило обнаружить эффект смещения фокусов дифракционных микролинз относительно фокуса рефракционной микролинзы, зависящее от числа уровней квантования. Методом вычислительного эксперимента найдены параметры субволновой структуры, обеспечивающей наименьшее отражение ИК излучения от алмазной пленки.

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

Список литературы диссертационного исследования кандидат физико-математических наук Головашкин, Димитрий Львович, 2000 год

1. Барский А. Б. Параллельные процессы в вычислительных системах. М.: Радио и связь, 1990.-255 с.

2. Березин И.С., Жидков Н.П. Методы вычислений. М.: Физматгиз, 1962.-Т.2.

3. Бицадзе A.B. Уравнения математической физики. М.: Наука, 1982.-146 с.

4. Борн М., Вольф Э. Основы оптики: Пер. с англ.- М.:Наука, 1973.- 720 с.

5. Боровиков В.А., Кинбер Б.Е. Геометрическая теория дифракции. -М.: Связь,1978.-248 с.

6. Вабищев П.Н. Численное моделирование. М.:МГУ, 1993.-152 с.

7. Ваганов Р.Б., Каценеленбаум Б.З. Основы теории дифракции . М.: Наука,1982.-272 с.

8. Вальковский В.А., Котов В.Е., Марчук А.Г., Миренков H.H. Элементы параллельного программирования. М.: Радио и связь, 1983.-239 с.

9. Вальковский В.А., Малышкин В.И. Синтез параллельных программ и систем на вычислительных моделях. Новосибирск: Наука, 1988.-126 с.

10. Ю.Валях Е. Последовательно-параллельные вычисления: Пер. с англ. М.: Мир, 1985.-456 с.

11. П.Велихов Е.П., Голубев B.C., Григорянец А.Г., Лебедев Ф.П., Николаев Г.А. Мощные газоразрядные CO2 лазеры и их применение в технологии. - М.: Наука, 1984.-105 с.

12. Взятышев В.Ф. Диэлектрические волноводы. М.: Советское радио, 1970.-206с.

13. И.Виноградова М.Б., Руденко О.В., Сухоруков А.П. Теория волн. М.: Наука,1979.-384 с.

14. Вычислительная оптика. Справочник под редакцией М.М. Русинова. JL: Машиностроение, 1984.-423с.

15. Годунов С.К., Рябенький B.C. Разностные схемы. -М.: Наука, 1973.-400 с.

16. Гончаренко A.M., Редько В.П. Введение в интегральную оптику. Минск:

17. Наука и техника, 1975.-152 с.148

18. Головашкин Д.JI. Моделирование распространения электромагнитного излучения методом конечных разностей// Методы компьютерной оптики/ под ред. В.А. Сойфера. М.:Физматлит, 2000.-С.224-231.

19. Головашкин Д.Л., Павельев B.C. Анализ прохождения электромагнитного импульса через антиотражающую структуру// Методы компьютерной оптики/ под ред. В.А. Сойфера. М. :Физматлит, 2000. -с.232-236.

20. Головашкин Д.Л. Разностная схема для уравнений Максвелла// Труды девятой межвузовской конференции. Самара, 1999.- с. 43-45.

21. Головашкин Д.Л., Сойфер В.А. Моделирование распространения электромагнитной волны посредством разностного решения уравнений Максвел-ла//Материалы международной молодежной школы по оптике. Саратов, 2000 (в печати).

22. Головашкин Д.Л., Дегтярев A.A., Сойфер В.А. Моделирование волноводно-го распространения оптического излучения в рамках электромагнитной теории//Компьютерная оптика. 1997.-N.17.-c.5-9.

23. Головашкин Д.Л., Дегтярев A.A. Алгоритм второго порядка точности по времени для решения уравнений Максвелла// Компьютерная оптика.- 1998.-N.18.-c.39-41.

24. Головашкин Д.Л., Павельев B.C., Сойфер В.А. Моделирование прохождения электромагнитной волны через алмазную антиотражающую структуру// Известия СНЦ РАН.- 1999.-т.1.-с.95-98.

25. Головашкин Д.Л., Сойфер В.А., Шустов В.А. Реализация параллельных вычислений при разностном решении уравнений математической физики// Известия СНЦ РАН.- 2000.-т.З.-(в печати)

26. Головашкин Д.Л., Павельев B.C., Сойфер В.А. Численный анализ прохождения света через антиотражающую алмазную структуру в рамках электромагнитной теории// Компьютерная оптика. 1999.-N.19.-c.44-46.

27. Головашкин Д.Л., Сойфер В.А. Анализ прохождения электромагнитного излучения через дифракционную линзу// Автометрия.-1999.-N6.-c. 119-121.149

28. Головин Б.А. Расчет характеристик и планирование параллельных вычислительных процессов.-М.: Радио и связь, 1983.- 272 с.

29. Голуб М.А., Казанский H.JL, Сойфер В.А. Математическая модель фокусировки лазерного излучения элементами компьютерной оптики // Научное приборостроение.-1993.-т.З.-Ы1.-с.8-28.

30. Горелик Г.С. Колебания и волны.-М.: Гос.изд.физ-мат. лит, 1959.-572 с.

31. Гантмахер Ф.Р. Теория матриц.-М.: Наука, 1988.-552 с.

32. Грицык В.В. Распараллеливание алгоритмов обработки информации в системах реального времени. Киев: Наук. Думка, 1981.- 215 с.32.3оммерфельд А. Оптика: Пер. с нем. М.: Иностранная литература, 1953г.-486с.

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

34. Каули Д. Физика дифракции. М.: Мир, 1979.-432 с.

35. Ключников A.C. Теория волновых процессов. Минск: БГУ, 1977.-176с.

36. Корн Г., Корн Т. Справочник по математике для научных работников и инженеров. -М.: Наука, 1977.-832 с.

37. Кравцов Ю.А., Орлов Ю.И. Границы применимости метода геометрической оптики и смежные вопросы// Успехи физических наук, 1980.-t.132.-N3.-с.475-496.

38. Кэй Д., Лэби Т. Таблицы физических и химических констант: Пер с англ. -М.: Физматгиз, 1962.-780с.

39. Ландау Л.Д., Лифшиц Е.М. Теория поля. М.: Наука, 1988.-Т.2.

40. Ландау Л.Д., Лифшиц Е.М. Электродинамика сплошных сред. М.: Физматгиз, 1959.-532С.

41. Ландсберг Г.С. Оптика. -М.: Наука, 1976.- 926 с.

42. Марешаль А., Франсон М. Структура оптического изображения. М.: Мир, 1964.-295 с.

43. Марков Г.Т., Петров Б.М., Грудинская Г.П. Электродинамика и распространение радиоволн. М.: Сов. радио, 1973.-376 с.150

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

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

46. Микеладзе Ш.Е. Численные методы математического анализа. М.: Гос. изд. Технико-теоретической литературы, 1953.-526 с.

47. Миренков H.H. Реализация продольно-поперечных прогонок на ВС "Минск-222"// Вычислительные системы. Новосибирск: ИМ СО АН СССР, 1968.-вып.30.-с.26-33.

48. Михлин С.Г., Смолицкий X.JI. Приближенные методы решения дифференциальных и интегральных уравнений .- М.: Наука, 1965.-383 с.

49. Молчанов И.Н. Введение в алгоритмы параллельных вычислений. Киев.:1. Наук. Думка, 1990.-127 с.

50. Морс Ф.М., Фешбах Г. Методы теоретической физики : Пер. с англ.-М.: Издательство иностранной литературы, 1958, 929 с.

51. Нестеренко Б.Б., Марчук В.А. Основы асинхронных методов параллельныхвычислений. Киев: Наук. Думка, 1989.- 171 с.

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

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

54. Пановский В., Филипс М. Классическая электродинамика: Пер. с англ. -М.: Физматгиз, 1963.-460 с.

55. Пярнпуу A.A., Хохлюк В.И. Реализация параллельных алгоритмов в прикладных задачах. М.: ВЦ АН СССР, 1988.-43 с.

56. Русинов М.М. Техническая оптика. М.: Физматгиз, 1961.-328 с.

57. Савельев И.В. Основы теоретической физики. -М.: Наука, 1975.-Т.1.

58. Садовничий В.А. Теория операторов. М.: Высшая школа, 1999.-367 с.

59. Самарский A.A. Теория разностных схем.- М.: Наука, 1989.-614 с.

60. Самарский A.A., Николаев Е.С. Методы решения сеточных уравнений. М.: Наука, 1978.-561 с.151

61. Самарский A.A. Введение в теорию разностных схем. -М.: Наука, 1971.-320 с.

62. Самарский A.A., Гулин A.B. Устойчивость разностных схем. -М.: Наука, 1973.-430 с.

63. Сена J1.A. Единицы физических величин и их размерности. М.: Наука, 1988.-430 с.

64. Сисакян И.Н., Сойфер В.А. Тонкая оптика, синтезируемая на ЭВМ// Физические основы прикладной голографии.-JI.: Физико-технический институт им. А.Ф. Иоффе АН СССР, 1984.-е. 142-164.

65. Снайдер А., Лав Дж. Теория оптических волноводов М.: Радио и связь, 1987.- 654 с.

66. Сойфер В.А. Введение в дифракционную микрооптику. Самара, 1996.-94 с.

67. Солимено С., Крозиньяни Б., Порто П.Ди. Дифракция и волноводное распространение оптического излучения. М.: Мир, 1989.-662 с.

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

69. Стреттон Д. Теория электромагнетизма :Пер. с англ. М.: Гостехиздат, 1948.-780 с.

70. Тамм И.Е. Основы теории электричества. М.: Наука, 1976.-614 с.

71. Тихонов А.Н., Самарский A.A. Уравнения математической физики. -М.: Наука, 1966.-724с.

72. Туров Е.А. Материальные уравнения электродинамики. -М.: Наука, 1983.-158с.

73. Фелсен Л., Маркувиц Н. Излучение и рассеяние волн. М.: Мир, 1978.-548с.

74. Физические величины. Справочник под редакцией И.С. Григорьева, Е.З. Мейлихова. М.: Энергоатомиздат, 1991.- 1232с.

75. Фок В.А. Проблемы дифракции и распространения электромагнитных волн. -М.: Советское радио, 1970.-316 с.

76. Хансперджер Р. Интегральная оптика. -М.: Мир, 1985.-384 с.

77. Хаусс X. Волны и поля в оптоэлектронике. М.: Мир, 1988.-430 с.152

78. Хиппель А.Р. Диэлектрики и волны :пер. с англ. -М.: Иностранная литература, 1960.-43 8с.

79. Хоар Ч.Р. Взаимодействующие последовательные процессы: пер. с англ. -М.: Мир, 1998.-264 с.

80. Янке Е., Эмде Ф., Лёш Ф. Специальные функции. М.: Наука, 1968.-344 с.

81. Ярив А. Введение в оптическую электронику. М.: Высшая школа, 1983.320 с.

82. Chu S.T., Huang W.P., Chaudhuri S.K. Simulation and analysis of waveguide based optical integrated circuits// Computer Physics Communications.- 1991.-Vol. 68.-P. 451-484.

83. Fia T. A numerical energy conserving method for the DNLS equation// Journal of computation physics.-1992.-N. 102.-P. 71-79.

84. Hatakoshi G., Fujima H.,Goto K. Waveguide grating lenses for optical couplers// Applied Optics.- 1984.- Vol. 23.-N. 1 l.-P. 1749-1753.

85. Kusuda Y., Saton A., Tanaka I. Application of high NA planar microlens (PML)//1.ternational Optical Design Conference. Rochester.-1994.- P. 280-284.

86. Li Y. Establishment of the maximum encircled energy in the geometrical focal plane// Opt Acta.-1984.- Vol. 31.-N. 10.-P. 1107-1118.

87. Mirotznik M., Prather D., Mait J. A hybrid finite element-boundary element method for the analysis of diffractive elements// Journal of modern optics.-1996.-Vol. 43.-N. 7.-P. 1309-1321.

88. Moharam M.G. and Gaylord Т.К. Diffraction analysis of dielectric surface-relief grating//JOSA.-1982.-N. 72.-P. 1385-1392.

89. Motamedi M.E., Southwell W.H., Gunning W.J. Antireflection surfaces in silicon using binary optics technology// Applied Optics.- 1993.- Vol. 31.-N. 22.-P. 43714373.

90. Pavelyev V.,Golovashkin D., Soifer V., Konov V., Kononenko V., Pimenov S., Prokhorov A. CVD diamond transmissive diffractive optics for C02 lasers// Diffractive Optics Jena 99.- P. 32-35.153

91. Raguin D.H., Morris G.M. Antireflection structured surfaces for the infrared spectral region// Applied Optics.- 1993.- Vol. 32.-N. 7.-P. 1154-1167.

92. Raguin D.H., Morris G.M. Analysis of antireflection-structured surfaces with continuous one-dimensional surface profiles// Applied Optics.- 1993.- Vol. 32.-N. 14.-P. 2582-2598.

93. Raguin D.H., Morris G.M. Design of 1-D anti-reflection structured surface using second-order effective medium theory// JOS A technical Digest.- 1992.-N. 9.-P. 44-46.

94. Raguin D.H., Morris G.M. Diffraction analysis of antireflection. surface-relief gratings on losslees dielectric surfaces// JOS A technical Digest.- 1990.-N. 15.-P. 122-123.

95. Raguin D.H., Morris G.M. Analysis of 1-D antireflection structured surfaces //JOSA technical Digest- 1991.-N. 17.-P. 153.

96. Rayleigh On the influence of obstacles arranged in rectangular order upon the propertis of a medium// Philos. Mag.- 1892.-N. 34.-P. 481-502.

97. Renaut R. Absorbing boundary conditions, difference operators, and stability// Journal of computation physics.- 1992.-N. 102.-P. 236-251.

98. Rytov S.M. The electromagnetic properties of finely layered medium// Soviet Phis.-1956.-N. 2.-P. 466-475.

99. Schweicher E. Introductory review to diffraction optics// Rev.X, 1985.-N. 116.

100. Swanson G., Veldkamp W. Binary lenses for use at 10.6 micrometers// Opt. Eng.-1985.- Vol. 24.-N. 5.-P. 791-795.

101. Wyrowski F., Bryngdahl O. Digital holography as part of diffractive optics// Report on Progress in Physics.- 1991.-Vol. 54, N. 12.-P. 1481-1571.