Совместное разностное решение уравнений Даламбера и Максвелла. Одномерный случай

СОВМЕСТНОЕ РАЗНОСТНОЕ РЕШЕНИЕ УРАВНЕНИИ ДАЛАМБЕРА И МАКСВЕЛЛА.

ОДНОМЕРНЫЙ СЛУЧАЙ

Головашкин Д.Л.1, Яблокова Л.В.2

1 Институт систем обработки изображений РАН,

2 Самарский государственный аэрокосмический университет имени академика С.П. Королёва

Аннотация

Предложена методика отыскания совместного разностного решения волнового уравнения и системы уравнений Максвелла, позволяющая совместить достоинства и избежать недостатков обоих упомянутых численных методов нанофотоники. В одномерном случае на тестовых примерах продемонстрированы сходимость такого решения, возможность наложения РМЬ-слоёв и задания падающей волны по технологии TF/SF.

Ключевые слова: уравнение Даламбера, уравнения Максвелла, разностная схема, РМЬ-слой, методика TF/SF.

Введение

Разностный метод решения уравнений Максвелла (FDTD) на настоящий момент является наиболее популярным инструментом математического моделирования в нанофотонике и нанооптике, о чём свидетельствует многочисленность реализаций метода, адаптированных к решению широкого класса задач на различных вычислительных архитектурах. Вместе с тем высокие требования к системным ресурсам (в частности, к объёму оперативной памяти) сдерживают применение метода в случаях, когда аппаратная база (например, бюджетные видеопроцессоры) характеризуется высоким быстродействием, но весьма ограниченным объёмом ОЗУ.

Указанный недостаток принято преодолевать учётом природы электромагнитного излучения (введением подвижных сеточных областей [1]), особенностей строения элементов нанофотоники (декомпозицией сеточной области [2]), разработкой новых приёмов программирования (метод пирамид [3]) или обращением к разностному решению волнового уравнения [4]. Комбинированию последнего подхода с FDTD-методом и посвящена предлагаемая работа.

По сравнению с разностным решением уравнений Максвелла, аналогичное решение уравнения Даламбера позволяет в двумерном случае на треть сократить объём выделяемой оперативной памяти (сохраняются два временных сечения электрического поля вместо трёх проекций электромагнитного) и на 10% снизить число арифметических операций (10 при вычислении поля в одном узле сеточной области вместо 11 при вычислении трёх сеточных функций в ячейке Yee). Впрочем, при рассмотрении трёхмерной задачи метод FDTD становится более предпочтительным, что ограничивает область применения развиваемого подхода двумерным случаем.

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

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

Перенесение современных методик, разработанных для решения уравнений Максвелла, на разностное решение волнового уравнения не всегда приводит к упрощению задачи. Так, в работе [5] авторы обращаются к заданию более сложного «прозрачного» источника вместо хорошо зарекомендовавшей себя технологии разделения полей TF/SF в силу неработоспособности последней при решении уравнения Даламбера.

Поэтому представляется целесообразным совместное разностное решение уравнений Даламбера и Максвелла, где первое привлекается для снижения требований к системным ресурсам, а участие второго позволяет использовать без модификации современные эффективные методики. В настоящей статье это методики задания поглощающих слоёв (РМЬ) и формирования падающей волны (TF/SF), рассмотренные на примере одномерного случая с целью верификации работоспособности заявленного подхода.

1. Согласование разностных решений уравнений Даламбера и Максвелла на границе сеточных подобластей Для одномерных уравнений Максвелла

Н'О

дИу

д1

д_Ех_ дz

дt

и

дz

запишем разностную схему Yee [б]

Яп+о.з гтп-о,э

V - И V

ук+О,5 ук+О,5

En - En

E"+1 - E” x

h

ЯП+О,Э 7ГИ+и,З

V -

ук+О,5 ук-О,5

h

где сеточная проекция электрического поля на ось X - E"t определена в узлах

{(4, zk): tn = nht, n = 0, 1, N = T/ht, zk=LW + khz, k = 0, .., K = LM/hz},

Б

О

О

h

h

z

Б

О

а сеточная проекция магнитного поля на ось Y -

Н„ +0,5 в узлах

Л+0,5 ^

{(/„+0,5, ^+0,5): 4+0,5 = („ + 0,5)к, п = 0, 1, N-1,

2к+0,5 = ь№ + (к + 0,5)кг, к = 0, .., А-1}

сеточной области Вм, наложенной (рис. 1) на область вычислительного эксперимента Бм (0 < / < т, Ья<г <Ья+Ьм), для простоты заданной в свободном пространстве.

Аналогично для одномерного волнового уравнения

д2 E,

2 д2E -2 x- _ 0

дt2 дx2

запишем разностную схему [7]

ЕП+1 _ 2E” + E

^n-1

h2

2 E” _ 2E” + E”

m+1_________________Am-1

h

где функция Е”х определена в узлах

{(/„, 2т): /„ = пки п = 0, 1, .., N = Шь 2т = ткг, т = 0, .., М= Ья/к}

сеточной области ВМ, наложенной (рис. 1) на Вя

(0 < / < т, 0 < г < Ья).

Г

0

-DL

-Dff-

~л0 0,5 1 1,5 2 ... К-0,5 К

1 ... М-1 М

—Q-O—0^гСЬг0-----л--П »

t I О--О.....О, п-0,5

о

D

w

г W

D

м

Lw+jLM z

Рис. 1. Объединение сеточных областей.

Квадратами изображены узлы для проекции электрического поля, окружностями - для магнитного

Задавшись целью совместного отыскания разностного решения уравнений (1) и (3), согласуем вычисления на и 0м, полагая значения Е" , опре-

у хм ' 1

делённые на Б^ при т = М (крайний правый узел), и Е„ , определённые на ОМ при к = 0 (крайний левый узел), равными (рис. 1). На границах объединённой области В = DWuDM установим электрическую стенку. За начальное условие примем отсутствие поля в момент времени / = 0.

Алгоритм перехода на следующий временной слой при реализации совместного решения состоит в определении сеточных функций в узлах Б^ и Вм в соответствии с табл. 1, где (5) имеет вид:

E - 2En + En

En _ JFn + Fn

Ex 2Exu ^ Ex„ hi

В (5) Е^ определено на области Б^, Е"Хм ] - на

Б№, остальные значения относятся к обеим областям.

Для подтверждения правомерности предложенного подхода к совместному решению на одной вычислительной области уравнений Максвелла и вол-

нового авторы провели четыре серии вычислительных экспериментов. При этом использовались операционная система Windows XP Professional SP3, вычислительный комплекс Octave 2.9.15 и процессор AMD Opteron 246.

Таблица І. Расчёт совместного решения (І) и (3) на временном слое n+l

Узлы сеточных областей Подобласти вычислительных областей Решаемые уравнения

m = 0 2 = 0, левая граница Б и Б№ E”+1 = 0 x0

0 < m < M-1 0 < 2 < LW, Бш без границ

m = M и k = 0 2 = LW, общая граница В№ и Б^

0 < k < K ^¥< 2 < LW+LM, Б^ без границ

k II 2 = + ^^, правая граница Б и Б^ E”+1 = 0 xK

Эксперименты проводились при различных значениях дискретизации сеточной области Q, Qt и QT, где первый параметр характеризовал число узлов сеточной области по пространству (приходящееся на одну длину волны); второй - количество узлов по времени (приходящееся на временной интервал, за который плоский волновой фронт в вакууме пройдёт расстояние в одну длину волны); третий - «длительность» запускаемого цуга в длинах волн. При этом они менялись от (10, 20, 5) - минимальных значений, удовлетворительно описывающих распространение плоской однородной волны в свободном пространстве, до (100, 200, 15) - соответствующих весьма низким величинам погрешности. Протяжённость области В составляла 20Х, что позволяло не учитывать влияние отражённой от границ волны в области регистрации.

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

Излучение вводилось в область посредством «жёсткого» источника [7], расположенного в центре.

Равномерная погрешность

є _max 1 - A,\\x100% k \ 1 k|l

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

M

M

2

h

ве плоской однородной монохроматической волны равен 1) на отрезке [10,5Х, 11,5Х].

Таблица 2. Погрешности е первой и второй серий экспериментов

Q, Qt QT

5 10 15

10, 20 0,7268 0,2434 0,1275

20, 40 0,2217 0,1080 0,0371

50, 100 0,0716 0,0234 0,0113

100, 200 0,0234 0,0083 0,0046

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

Во второй серии решалось волновое уравнение при 1^ = 0 и D = Бш. Дискретизация сеточной области для электрической компоненты поля, расположение и тип источника, размеры вычислительной области и области регистрации погрешности совпадали с предыдущим случаем.

Погрешности всех вычислительных экспериментов второй серии в точности совпали с результатами, представленными в табл. 2, что неудивительно. Известна процедура получения волнового уравнения из уравнений Максвелла [8].

Третья серия относится к совместному решению уравнений (1) и (3). Совокупная протяжённость объединённых вычислительных областей и Б составила 20Х (по 10Х на каждую). Источник располагался в области решения волнового уравнения на расстоянии половины длины волны от границы раздела (рис. 2а), область регистрации [10Х, 11Х] помещалась в Бм. Таким образом, расстояние от источника до этой области не менялось по сравнению с предыдущими экспериментами.

Рис. 2. Вычислительные области третьей (а) и четвёртой (б) серий экспериментов. Звёздочкой обозначен «жёсткий» источник, в обведённой пунктиром области регистрируется погрешность

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

Причиной тому является различие в дифференциальных шаблонах схем (2) и (4), несущественное при раздельном решении (оба шаблона центральные

и характеризуются одинаковыми порядками аппроксимации) и важное при совместном (в них участвуют разные сеточные функции).

Таблица 3. Погрешности е третьей серии экспериментов

Q, Qt QT

5 10 15

10, 20 0,5533 0,1907 0,0926

20, 40 0,3531 0,0562 0,0313

50, 100 0,0687 0,0217 0,0112

100, 200 0,0235 0,0081 0,0047

Для определения значения напряжённости электрического поля Е"х+ на стыке DW и Df по (5) используется центральная по пространству разность, в которую входят: общее значение E" , E" из DW и

А ^ ХМ хм-1 п

е; из DM . При этом Е"м 1 определяется из (4) по тому же шаблону, а Е; - уже через сеточные значения

магнитного поля (2), взятые на другом расстоянии от узла k = 1 (хотя шаблон при этом тоже центральный).

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

При расположении источника в области D на расстоянии половины длины волны от границы раздела (четвёртая серия экспериментов, рис. 2б), а области регистрации [9Х, 10Х] в DW результаты отличаются от представленных в таблицах 2 и 3 весьма незначительно - тоже на одну десятитысячную процента для густой сеточной области.

2. Наложение поглощающего слоя при совместном решении Важной задачей при моделировании дифракции электромагнитного излучения разностным методом является помещение поглощающих слоёв (PML-perfectly matched layer [9]) у границ вычислительной области для исключения отражения от последних. По замыслу авторов настоящей работы, рассмотрение задачи дифракции на оптическом элементе уместно связывать с решением волнового уравнения в области непосредственного нахождения этого элемента, а организацию поглощения у границ вычислительной области (вне элемента) - с решением уравнений Максвелла, моделирующим PML - слой. При этом сокращаются требования к системным ресурсам ЭВМ по сравнению с FDTD-методом (в двумерном случае) и используется современная методика наложения PML-слоя вместо устаревшего подхода Мура [10], традиционно сопровождающего решение волнового уравнения.

Представляя уравнения Максвелла в поглощающем слое вместо (1) записывают:

дНу дЕх

ц0-------у + а Ну =-------------х

0 дt у дz

дЕх г, дНу

є0— +аЕх =---------------

0 дt х дz

где второе слагаемое в левой части определяет поверхностную плотность электрического и магнитного (воображаемого) токов. Следующее соотношение электрической и магнитной проводимостей ст/е0 = СТ /ц и плавное возрастание значения ст от начала слоя к концу (при правильно подобранных параметрах) обеспечивает поглощение излучения без отражения [9].

Простейшая разностная аппроксимация уравнений Максвелла в поглощающем слое имеет вид:

/-т+0,5 ц-т-0,5

Т т-0,5

^0

Я:+0,5 ттп.

- — Н -

yk+0,5 yk+0,5

h

Em+l — Em

і tt:

t0k+0,5 H-k+0,5

e: — e:

h

+ akE:t = —

Н:+0,5 — h:+0,5 yk+0,5 yk—0,5

h

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

При постановке вычислительных экспериментов расположим поглощающий слой длины X у правого края Б.

Тогда с задаётся на отрезке [Ь№+Ьм - X, ЬТ+Ьм] по формуле:

а = ат

X

где сттах и q зависят от дискретизации сеточной области согласно табл. 4 из [11].

Верифицируя предположение о возможности совместного решения волнового уравнения и системы уравнений Максвелла с поглощающим слоем, сначала (для получения эталонного решения) рассмотрим случай Б = Бм (пятая серия экспериментов, рис. 3а); = 0, Ьм = 12Х. Для исследования отра-

жённой волны зададим «прозрачный» источник [2] (в точке г = 9,5Х) и будем регистрировать в подобласти [10Х, 11Х] величину 8= тах| Л°ктр\ (Лотр- ком-

к I I

плексная амплитуда отражённой волны), являющуюся абсолютной погрешностью вычислительного эксперимента (отражённой от PML-слоя волны в идеальном случае быть не должно).

2>м____________ ._____________,

а)------@-------г - -| |_

Х/2

6)-----------------------@------------------------|=

з-

шрм

Рис. 3. Вычислительные области пятой (а) и шестой (б) серий экспериментов. Символом «@» обозначен «прозрачный» источник, в обведённой пунктиром области регистрируется погрешность, в закрашенной расположен поглощающий слой

Отметим удовлетворительное функционирование поглощающего слоя; даже на самой грубой сеточной области (<2 = 10, Qt = 20) при работе с не вполне монохроматической волной ^Т = 5) модуль комплексной амплитуды отражённого поля составил треть

процента от модуля амплитуды (теоретического) падающей волны (табл. 4).

Таблица 4. Погрешности 8*100 пятой серии вычислительных экспериментов

Q, Qt QT

5 10 15

10, 20 0,3546 0,1032 0,1463

20, 40 0,1713 0,1238 0,1261

50, 100 0,1127 0,1068 0,1068

100, 200 0,0129 0,0115 0,0115

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

При отыскании совместного разностного решения (1) и (3) (шестая серия экспериментов) примем = 10Х, Ьм = 2Х, расположенный в том же месте вычислительной области «прозрачный» источник теперь находится в БТ (рис. 3б).

Таблица 5. Погрешности 8*100 шестой серии вычислительных экспериментов

Q, Qt QT

5 10 15

10, 20 0,1287 0,1312 0,1483

20, 40 0,1645 0,1232 0,1322

50, 100 0,1136 0,1068 0,1068

100, 200 0,0127 0,0115 0,0115

Сравнивая результаты вычислительных экспериментов (табл. 5) с данными пятой серии (табл. 4), отметим их совпадение в диапазоне Q = 50, Qt = 100, QT = 10 - Q = 100, Qt = 200, QT = 15 и незначительность отличий при остальных параметрах, что свидетельствует о возможности совместного разностного решения уравнений Даламбера и Максвелла в случае, когда последние используются для моделирования поглощающего слоя.

3. Задание падающей волны по технологии TF/SF при совместном решении

Формирование падающей волны через задание «прозрачного» источника (как это было сделано в предыдущем параграфе) не всегда удобно. Излучение от такого источника распространяется не только в желаемом исследователю направлении, но и в противоположном (этим, например, характеризуется пакет MEEP [12], реализующий FDTD-метод). Кроме того, поглощению у границ области в этом случае подлежит результирующее поле, состоящее из падающего и рассеянного. Всё это зачастую приводит к увеличению погрешности вычислительного эксперимента.

Задание падающей волны по классической технологии TF/SF (total field / scattered field [13]) в рамках FDTD-метода позволяет избежать распространения в противоположном направлении и появления падаю-

є

0

щего поля в области расположения РМЬ-слоёв (там окажется только рассеянная волна).

Основой технологии является разделение полей при решении (1) на результирующее (в подобласти кЬ < к < ^ области БМ) и рассеянное (к < кЬ и к > кр), при котором в соответствующих узлах вместо (2) записываются разностные уравнения

укі-0,5 Укі-0,5

(Еп -Еп )-Еп

\ Укі Укі } Укі -1

Еп+1 - Еп

Хкі Хкі

Нп+ 0,5

Укі+0,5

_(Н 0,5 + Н 0,5

\ Укі-0,5 Укі-0,5

)

НП+0,5 нп-0,5

Ук% +0,5 Ук% +0,5

Еп -

Х \ Хкк ХкЕ )

Еп+1 - Е

Хкк Хкк

п+1

Хк

(Нп+ 0,5 + Н п+0,5

Укд+0,5 укд+0,5 /

- нп

К

hx

где сеточные функции под тильдой соответствуют аналитически [13] или численно [14] заданному падающему полю.

В работе [10] высказывается предположение о возможности использования приведённой методики при решении волнового уравнения. Однако авторам настоящей статьи не известно ни одной успешной реализации (включая собственную попытку [5]) этого замысла. Выходом из сложившейся ситуации может стать совместное решение уравнений Даламбера и Максвелла, при котором область БТ содержит исключительно общее поле, а к Б относится граница раздела общего и рассеянного полей.

Предваряя экспериментальную проверку этого утверждения, проведём серию вычислительных экспериментов (седьмую) для определения погрешностей FDTD-метода при задании падающей волны по методике TF/SF. Примем для этого = 0, Ьм = 20Х, Б = Бм, задав результирующее поле на отрезке [9Х, 11Х] (рис. 4а).

Ом

а)------

рассеянное поле 6)-----------#.

результирующее поле

Вм

-#-----------

рассеянное поле -#-----------

'I О”

Рис. 4. Вычислительные области седьмой (а) и восьмой (б) серий экспериментов. Символом «#» обозначена граница раздела результирующего и рассеянного полей, в обведённой пунктиром области регистрируется погрешность

Регистрация погрешности є производится на отрезке [9,5Х, 10,5Х] в центре D. Поглощающие слои у границ области не используются. При численном задании падающей волны в (6) результаты совпадают с представленными в табл. 2, для аналитического способа задания - приведены в табл. 6.

Таблица 6. Погрешности е седьмой серии экспериментов при аналитическом задании падающей волны

Q, Qt QT

5 10 15

10, 20 4,0678 3,1662 2,2889

20, 40 1,2341 0,733 0,7927

50, 100 0,4902 0,2658 0,1752

100, 200 0,126 0,0899 0,0793

Существенное отличие (на порядок) результатов из обеих таблиц традиционно объясняется [14, 15] взаимной компенсацией погрешностей двух разностных решений (основного и вспомогательного) задачи (1), используемых в рамках численной методики TF/SF. При аналитическом задании волны в (6) такой компенсации не происходит.

Совместное решение (1) и (3) (восьмая серия экспериментов) сопроводим разбиением О на две подобласти 0м = ОМ и ОМ , где 0^ определена на отрезке [0, Iм], ОМ - на [Iм+Ь№, 21м+і\ а 0№ длиной Ь№ поместим между ними (рис. 4б). Примем при этом Ьм = 9,5Х, Ь№ = X и будем формировать значение є на том же отрезке, что и ранее. Теперь он

тл№

совпадает с О , а левая и правая границы результирующего поля по-прежнему принадлежат 0м: левая находится в , правая - в . Таким образом, в рамках одного подхода реализуется технология TF/SF через (2), (6) и возможность исследования дифракции на оптическом элементе (его можно поместить в посредством решения (4).

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

Таблица 7. Погрешности е восьмой серии экспериментов

Q, Qt Задание падающей волны в

Аналитическое

QT

5 10 15

10, 20 3,2787 2,1234 2,0421

20, 40 1,7547 1,0815 0,5867

50, 100 0,4257 0,2717 0,1958

100, 200 0,1527 0,0949 0,0709

Q, Qt Численное

QT

5 10 15

10, 20 0,6168 0,3302 0,1740

20, 40 0,3476 0,0662 0,0537

50, 100 0,0668 0,0189 0,0112

100, 200 0,024 0,0081 0,0046

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

0

К

К

Х

к

0

К

К

Х

Ц

0

К

К

Х

R

к

0

(табл. 6, 7), при численном - несколько десятитысячных (табл. 2, 7). Очевидно, большая разница в первом случае обусловлена неудачной (аналитической) методикой задания падающей волны, которая сама по себе характеризуется значительной погрешностью. Наличие упомянутой разницы (впрочем, весьма несущественной, что свидетельствует о правомерности совместного решения) объясняется причиной, выявленной при анализе третьей серии экспериментов.

Заключение

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

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

Благодарности Исследования выполнены при поддержке грантов Президента РФ МД-6809.2012.9, НШ-4128.2012.9 и грантов РФФИ №11-07-13164-офи-м-2011-РЖД, № 10-07-00553.

Литература

1. Fidel, B. Hybrid ray-FDTD moving window approach to pulse propagation / B. Fidel, E. Heyman, R. Kastner and R.W. Zioklowski // Journal of Computational Physics. -1997. - Vol. 138, Issue 2. - P. 480-500.

2. Головашкин, Д.Л. Декомпозиция сеточной области при разностном решении уравнений Максвелла / Д.Л. Головашкин, Н.Л. Казанский // Математическое моделирование. - 2007. - Т. 19, № 2. - C. 48-58.

3. Головашкин, Д.Л. Решение сеточных уравнений на графических вычислительных устройствах. Метод пирамид / Д.Л. Головашкин, А.В. Кочуров // Вычислительные технологии. - 2012. - Т. 17, № 3. - С. 55-69.

4. Козлова, Е.С. Моделирование распространения короткого двумерного импульса света / Е.С. Козлова,

В.В. Котляр // Компьютерная оптика. - 2012. - Т. 36, № 2. - C. 158-164.

5. Астраханцева, М.А. Методика формирования падающей волны при разностном решении двумерного волнового уравнения / М.А. Астраханцева, Д.Л. Головашкин, Л.В. Яблокова // Труды седьмой Всероссийской научной конференции с международным участием «Математическое моделирование и краевые задачи», Ч. 2. - Самара, 2010. - С. 298-301.

6. Yee, K.S. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media / K.S. Yee. // IEEE Trans. Antennas Propag. - 1966. -AP-14. - P. 302-307.

7. Taflove, A. Computational Electrodynamics: The Finite-Difference Time-Domain Method: 2nd. ed. / A. Taflove,

S. Hagness // Boston: Arthech House Publishers, 2000. -852 p.

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

9. Berenger, J.-P. A perfectly matched layer for the absorption of electromagnetic waves / Jean-Pierre Berenger // Journal of computational physics. - 1994. - Vol. 114. -P. 185-200.

10. Mur, G. Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations / G. Mur / IEEE Trans. Electromagnetic Compability. - 1981. - Vol. 23. - P. 377-382.

11. Головашкин, Д.Л. Методика формирования падающей волны при разностном решении уравнений Максвелла (двумерный случай) / Д.Л. Головашкин, Н.Л. Казанский // Автометрия. - 2007. - Т. 43, № 6. - С. 78-88.

12. http://ab-initio.mit.edu/wiki/index.php/Meep .

13. Umashankar, K. A novel method to analyze electromagnetic scattering of complex objects / K. Umashankar, A. Taflove // IEEE Trans. Electromagn. Compat. EMC-24, 1982. - P. 397-405.

14. Prather, D.W. Formulation and application of the finite-difference time-domain method for the analysis of axially symmetric diffractive optical elements / D.W. Prather, S. Shi // J. Opt. Soc. Am. A. - 1999. - Vol. 16, N 5. -P. 1131-1142.

15. Головашкин, Д.Л. Методика формирования падающей волны при разностном решении уравнений Максвелла. Одномерный случай / Д.Л. Головашкин, Н.Л. Казанский // Автометрия. - 2006. - Т. 42, № 6. -

C. 78-85.

References

1. Fidel, B. Hybrid ray-FDTD moving window approach to pulse propagation / B. Fidel, E. Heyman, R. Kastner and R.W. Zioklowski // Journal of Computational Physics. -1997. - Vol. 138, Issue 2. - P. 480-500.

2. Golovashkin, D.L. Mesh domain decomposition fdtd /

D.L. Golovashkin, N.L. Kazanskiy // Мathematic modeling. - 2007. - Vol. 19, N 2 - P. 48-58. - (In Russian).

3. Golovashkin, D.L. Solution of difference equations in graphical D.L. Golovashkin, A.V. Kochurov // Computing technologies. - 2012. - Vol 17, N 3. - P. 55-69. - (In Rus-sian).computing devices. Method of pyramids /

4. Kozlova, E.S. Model operation of distribution of a short twodimensional impulse of light / E.S. Kozlova, V.V. Kotlyar // Computer optics. - Vol. 36, N 2. - P. 158164. - (In Russian).

5. Astrakhantseva, M.A. Method of formation of the incident wave in a two-dimensional finite-difference solution of the wave equation / M.A. Astrakhantseva, D.L. Golovashkin, L.V. Yablokova // Proceedings of the Seventh All-Russian Scientific Conference with International Participation “Mathematical modeling and boundary problems”, Issue 2. - Samara, 2010. - P. 298-301.

6. Yee, K.S. Numerical solution of initial boundary value problems involving Maxwell’s equations in isotropic media / K.S. Yee. // IEEE Trans. Antennas Propag. - 1966. -AP-14. - P. 302-307.

7. Taflove, A. Computational Electrodynamics: The Finite-Difference Time-Domain Method: 2nd. ed. / A. Taflove,

S. Hagness // Boston: Arthech House Publishers, 2000. -852 p.

8. Nikolsky, V.V. Electrodynamics and radio waves propagation / V.V. Nikolsky, T.I. Nikolskaya. - Moscow: “Nauka” Publisher, 1999. - 544p. - (In Russian).

9. Berenger, J.-P. A perfectly matched layer for the absorption of electromagnetic waves / Jean-Pierre Berenger // Journal of computational physics. - 1994. - Vol. 114. - P. 185-200.

10. Mur, G. Absorbing boundary conditions for the finite-difference approximation of the time-domain electromagnetic field equations / G. Mur / IEEE Trans. Electromagnetic Compability. - 1981. - Vol. 23. - P. 377-382.

11. Golovashkin, D.L Technique of formation of the incident wave at the difference solution of Maxwell's equations (twodimensional case) / D.L. Golovashkin, N.L. Kazanskiy // Journal “Avtometriya”. - 2007. - Vol. 43, N 6. - P. 78-88. -(In Russian).

12. http://ab-initio.mit.edu/wiki/index.php/Meep .

13. Umashankar, K. A novel method to analyze electromagnetic scattering of complex objects / K. Umashankar, A. Taflove // IEEE Trans. Electromagn. Compat. EMC-24, 1982. - P. 397-405.

14. Prather, D.W. Formulation and application of the finite-difference time-domain method for the analysis of axially symmetric diffractive optical elements / D.W. Prather,

S. Shi // J. Opt. Soc. Am. A. - 1999. - Vol. 16, N 5. -P. 1131-1142.

15. Golovashkin, D.L. Technique of formation of the incident wave at the difference solution of Maxwell's equations. Dimensional case / D.L. Golovashkin, N.L. Kazanskiy // Journal “Avtometriya”. - 2006. - Vol. 42, N 6. - P. 78-85. - (In Russian).

JOINT FINITE-DIFFERENCE SOLUTION OF THE DALAMBER AND MAXWELL’S EQUATIONS.

ONE-DIMENSIONAL CASE

D.L. Golovashkin \ L. V. Yablokova 2

1 Image Processing Systems Institute of the RAS,

2 S.P. Korolyov Samara State Aerospace University

Abstract

The technique of searching of the collateral finite-difference solution of a wave equation and set of equations of Maxwell is offered, allowing to combine advantage and to avoid shortcomings of both made mention numerical methods of nanophotonics. In an one-dimensional case on test examples convergence of such decision, possibility of imposing of PML layers and a task of an incident wave on the TF/SF technology are shown.

Key words: Delamber equation, Maxwell equation, difference scheme, PML layer, TF/SF technique.

Сведения об авторах

Головашкин Димитрий Львович, доктор физико-математических наук, доцент, ведущий научный сотрудник Института систем обработки изображений РАН. Область научных интересов: разностное решение уравнений Максвелла (FDTD-метод); дифракционная оптика; векторные и параллельные матричные вычисления.

E-mail: dimitriy@smr.ru.

Dimitry Lvovich Golovashkin, Doctor of Physical and Mathematical Sciences, Associate Professor, Leader Researcher of Image Processing Systems Institute Russian Academy of Scie-nes. Scientific interests: FDTD-method, subwave optics, vector and parallel algorithms for matrix computation.

Яблокова Людмила Вениаминовна, 1978 года рождения. В 2000 году окончила Самарский государственный педагогический университет. Работает в СГ АУ (Самарский государственный аэрокосмический университет имени академика С.П. Королёва) старшим преподавателем кафедры прикладной математики. Область научных интересов: разностные схемы.

E-mail: lyablokova@yandex. ru .

Lyudmila Veniaminovna Yablokova, (b. 1978) graduated from Samara State Pedagogical University in 2000. She works as the senior teacher at the S.P. Korolyov Samara State Aerospace University Applied Mathematics sub-department. Research interests: finite-difference approximations.

Поступила в редакцию 26 сентября 20I2 г.