CC BY
76
УДК 517.6
Ю. Г. Смирнов, М. Ю. Медведик, М. А. Максимова
РЕШЕНИЕ ЗАДАЧИ ДИФРАКЦИИ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ НА ЭКРАНАХ СЛОЖНОЙ ФОРМЫ1
Аннотация. Рассмотрена задача дифракции электромагнитной волны на бесконечно тонком идеально проводящем экране. Задача сведена к интегродиф-ференциальному уравнению. Предложен метод Галеркина. Получены численные результаты.
Ключевые слова: электромагнитная задача дифракции, краевая задача, интегральное уравнение, численный метод.
Abstract. The article considers a problem of diffraction of electromagnetic wave on the infinitely thin and perfectly conducting screen. The problem is reduced to an in-tegro-differential equation. The authors suggest the Galerkin method and present the numerical results.
Key words: electromagnetic problem of diffraction, boundary value problem, integral equation, numerical method.
Введение
Решение трехмерных векторных задач дифракции является в настоящее время одной из самых актуальных проблем в электродинамике. Решение этих задач с приемлемой для практики точностью на электродинамическом уровне строгости математическими методами требует очень большого объема вычислений и часто невозможно даже на самых современных компьютерах. Существующие в настоящее время многочисленные пакеты прикладных программ (Ansis, Quickwave и т.д.) не дают требуемой точности, поскольку используют методы конечных элементов. И особенно остро стоит проблема решения задач дифракции в резонансном диапазоне частот. При решении краевых задач в неограниченных областях конечно-разностные методы и методы конечных элементов встречают принципиальные трудности: область, в которой решается задача, должна быть сделана конечной. Такая редукция области приводит к появлению неконтролируемой ошибки, причем размеры области для ее уменьшения должны быть достаточно велики. Конечно-разностные методы и методы конечных элементов в такой ситуации обычно приводят к очень большим, но разреженным матрицам (порядка 109 и более). Одним из перспективных методов решения задач в ограниченных областях является метод поверхностных интегральных уравнений [1]. В этом случае задача сводится к интегральному или инте-гродифференциальному уравнению в области неоднородности, которая по размерам существенно (на порядки) меньше области решения задачи, в случае применения конечно-разностных методов и методов конечных элементов. После дискретизации получается конечномерная система уравнений с плотной (заполненной) матрицей. Таким образом, второй путь приводит
1 Работа выполнена при поддержке ФЦП «Научные и педагогические кадры
инновационной России» на 2009-2013 гг., соглашение № 14.В37.21.1950.
к необходимости решать системы уравнений с плотными матрицами, но существенно меньших порядков.
1. Постановка задачи
Пусть AcR2 ={хз =0}сR3 - ограниченная область с кусочно-
и и тл и
гладкой границей Г , состоящей из конечного числа простых дуг класса С , сходящихся под углами, отличными от нулевого. Задача дифракции стороннего монохроматического электромагнитного поля E0, H0 на бесконечно
тонком идеально проводящем экране A, расположенном в свободном про-
2 2 —1
странстве с волновым числом к, к =ю ц(е + /аю ), Im к > 0 (к Ф 0), состоит в определении рассеянного электромагнитного поля
E,H е С2 (R3 \ A) р|С(R \ Г§) f|C(R—\ Г§), (1)
S>0 S>0
удовлетворяющего:
- однородным уравнениям Максвелла:
Rot H = -ikE,
Rot E = ikH, x е R3\ A ; (2)
- краевым условиям для касательных составляющих электрического
поля на поверхности экрана:
ex|q=-et0|a ; (3)
- условиям конечности энергии в любом ограниченном объеме пространства:
E, H е 4c (R3) (4)
- условиям на бесконечности (условиям Зоммерфельда):
d ( E Л ( E Л —1 ( E Л
dr
H
— ik
v H
= o(r 1), = O(r 1), r: |x| ^<~. (5)
тт 77ПОЛН 77О , гр ттПОЛН тт0 . тт
Для полного поля Е = Е + Е, Н = Н + Н .
Будем предполагать, что все источники падающего поля находятся вне экрана П так, что для некоторого 8> 0
Е0 е С“(П8), П = {х :| х - у |< 8,у еП}, (6)
откуда следует, что
Ет°Пе С - (П). (7)
Обычно падающее поле - это либо плоская волна, либо электрический или магнитный диполь, расположенный вне П. В этих случаях условия (6), (7) выполнены. Поле Е°,Н0 является решением системы уравнений Максвелла в свободном пространстве без экрана.
Имеет место утверждение единственности для поставленной задачи. Утверждение 1 [2]. Задача (1)-(5) при 1т к > 0, к Ф 0, имеет не более одного решения.
Заметим, что утверждение сохраняет силу и для системы из конечного числа непересекающихся экранов.
2. Свойства решений уравнений электрического поля
Ниже кратко приведем основные результаты о разрешимости задачи (1)-(5). Для изучения задачи дифракции на экране П введем векторное пространство распределений Ж . Положим для любого вещественного 5 [3]
Н (П):={и|п :и е Н5 (2)},
Н5 (П) :={и е Н5 (2): 8ирр и сП}.
Скалярное произведение и норма в Н5 (2) определяются обычным образом:
(и,1)=|И2 =(и,и)5; ®:=(1+Й2) .
Здесь и всюду ниже, где не указана область интегрирования, подразумевается интеграл по Д2 . Н5 (п) является (замкнутым) подпространством
Н5 (2) с индуцированными скалярным произведением и нормой. Далее
Н5 (П) = Н5 (2)/Н(п); в Н5 (П) вводится скалярное произведение и
норма факторпространства. Пространства Н 5 (П) и Н5 (п) антидвойствен-ны друг к другу при всех 5 е Д; Н5 (п) можно получить замыканием С0° (П) в пространстве Н5 (2) [3].
В дальнейшем нас будут интересовать главным образом пространства вектор-функций, поэтому через и, V будем обозначать векторы и = (1,и2) ,
V = (11,V2) и т.д. При этом Н в записи и е Н уже понимается как декартово произведение двух экземпляров пространства Ня со скалярным произведением и нормой
(u,1 )5 = (иь 11 )5 + (u2,12 )5 = {(I)25 и• ОД ё^
И2=1Ы2+1Ы2 = {( ^>251и (^)|2 ё
Сохраним те же обозначения для пространств в векторном случае, так как во всех ситуациях из контекста ясно, о каком пространстве идет речь.
Определим гильбертово пространство Ж = Ж(п) как пополнение С0° (П) по норме
фпЛ 1 ^1
со скалярным произведением
И2 = /ж 1“'®l d^ + \w$-U©I di
( ( W
d 5,
JJ
где и обозначает преобразование Фурье распределения и .
Утверждение 2 [2]. Ж = {и е Н-1/2(П): и е Н-1/2(П)}.
Норма в Ж может быть записана в виде
IMlW =1 “II—1/2 +l ldlv “II—1/2.
Определим пространства w1 и w2 как подпространства w с условиями:
W := {“е w: V5 51“1 (5) + ^ (5) = 0};
W, := {и е W: V5 (5) — 5U (5) = 0},
“/ (5)е С~(2) , и определение имеет смысл.
Для пространства W':={jW(a))' антидвойственного к W , имеем
Утверждение 3 [2]. Пространство W' определяется следующим образом:
W' = { f A : f е H—1/2 (R2), f —-dXLе H—1/2 (r2 ).
3. Векторные потенциалы и представление решений
Наиболее естественным и эффективным способом решения задач (1)-(5) является представление полей E, H с помощью векторного потенциала [2]
E = ik—1 ^Graddiv ((1“) + к2 A1U), (8)
H = Rot (A1u); (9)
1 гегк1х—y| , . 3 -
A1U = — \-j----j-u (y )dy; x = (xb X2, X3 )е R \ A . (10)
4rcJ x — y A 1 Л
/1 2 \T
Вектор и = I u , и ) имеет смысл плотности тока на поверхности экрана A .
Будем предполагать, что
и еW(A) (11)
и
u,divuе С1 (A). (12)
Известно, что [2]
div (A1u) = A1 (div u), xg A; uеW.
Тогда (8) эквивалентно
E = ik—1 (rad A1 (divu) + к2 Лц. ), xе R3 \ A. (13)
Краевое условие (3) приводит к интегродифференциальному уравнению для u . Опуская точку x на A , из (13) получаем
gradA(divu) + k2Au = f; x = (xbx2)еА; (14)
Au =
k|x—y|
I \x — y u (y/dy ; (i5)
й
f = 4^ikET0A; f е С “(A). (16)
Здесь операция grad рассматривается в R .
4. Базисные функции
В качестве базисных функций будем использовать прямоугольные функции-крышки, заданные ниже. Функции-крышки определяются для пар смежных прямоугольных ячеек сетки, разбитой на прямоугольники, как изложено ниже [4]. Для пар прямоугольников Р+ и Р , которые имеют общую грань In, функция формы fn (x) определяется как
Фп (x) = <
L + л. + +
------(x — p )u , если x е Р ,
A+
L — - — _
------(x — p )u , если x е Р
A—
где Ь является длиной грани 1п; Л+ и А есть площади Р+ и Р соответ-
. — Л + Л —
ственно, а р , р и единичные векторы и , и определены, как показано на рис. 1.
5. Метод Галеркина
Для операторных уравнений в гильбертовых пространствах проекционный метод, строящийся с помощью ортопроекторов на конечномерные под-
пространства, приводит к методу Галеркина. Пусть X и Y - гильбертовы подпространства и A: X ^ Y - инъективный линейный ограниченный оператор. Пусть Xn с X и Yn с Y - последовательность подпространств таких, что dim Xn = dim Yn = n , и пусть Pn : Y ^ Yn - последовательность ортопроекторов. Тогда фп е Xn будет приближенным решением уравнения Аф = f с помощью проекционного метода, образованного посредством выбора пространств Xn и проекторов Pn , тогда и только тогда, когда
{АФп,g) = (f,g), для любого g е Yn . (17)
h
Рис. 1. Базисные функции «rooftop»
Действительно, выражение (17) эквивалентно уравнению Pn (АФп - f). Уравнение (17) называют уравнением Галеркина.
Допустим, что подпространства Xn и Yn являются линейными оболочками базисных и тестовых функций Xn = span[z1,..., zn }, Yn = span{v1,...,vn } . Представим фп в виде линейной комбинации базисных функций
n
ФП = £ ' Ynzn k=1
и подставим это выражение в формулу (17).
В результате получим эквивалентную систему линейных алгебраических уравнений порядка n
£Yk(Azk,vj)=(fvj) )=1,...,n, (18)
k=1
относительно неизвестных коэффициентов Yk .
6. Свойство аппроксимации подпространств базисных функций «rooftop»
Рассмотрим вопрос об аппроксимации непрерывно-дифференцируемой (векторной) функции f (x,y) = (f (x,y),f (x,y)), f е cl (П)
f е Cq (П), f е C0 (П)), в прямоугольнике П = [ 0, a ] х [ 0, b ] «rooftop» базисными функциями ф j (x, y) по методу, предложенному в статье [4]. Рассмотрим в П равномерную прямоугольную сетку с шагами Hq и Н по осям x и y с узлами My = f,yj), xt = ihQ, yj = ,/% (i = 0,...,Nq, j = 0,...,N2),
hi =
N1
h2 =
n2
■. Базисную функцию ф j (x, y), отвечающую ребру j, опре-
делим по правилу (рис. 2):
фj (^y/ =
(j /--
Г.+ j ''
sj
x,G|-— в П-;
sj- j
фj (x,y/ =
G y-yi(j -)— в п+.
sj
G, y2j'-- y / в П-, 2 ’sj- j
(19)
Ф j = 0 вне прямоугольников Пу, П-. Здесь lj - длина j -го ребра; S j -
площадь прямоугольника П .
Рис. 2. Вид базисных функций на плоскости
Рассмотрим конечномерное подпространство XN = span{^l,..,фN} , являющееся линейной оболочкой базисных функций фj , 1,...,N , где N - количество внутренних ребер сетки. Нетрудно проверить, что ф j ЕЖ (п), XN с Ж . Имеет место следующий результат.
b
a
Теорема 1 [5]. Для любого фєW верна оценка inf ||у-ф|^ ^0,
rnf ||¥-ф^ ^C0 ( + h2 —ф|c2(п/,
(20)
где С0 не зависит от Иу и И2, если фе Ж п С2 (п).
7. Расчет матричных элементов для экрана канонической формы
Для решения реальных физических задач требуется вычислять матрицы как можно большего размера, поэтому остро стоит вопрос об оптимизации расчета базисных элементов. Наиболее естественным способом сокращения вычислений является учет матричной симметрии элементов, также удобно применить параллельный алгоритм. Рассматриваемые матрицы являются ган-келевыми. Рассмотрим процедуру заполнения матричных элементов для прямоугольного экрана канонической формы.
Произведем заполнение матрицы, используя матричное уравнение (Афп, g) = (/, g). Каждый элемент матрицы получается путем вычисления четырехкратного интеграла
имеющего особенность в области интегрирования. Здесь х = (ху, Х2),
У = (Уі,У2), а G(х,у) - известная функция G(х,у) =--------------------------------:-:-. В каче-
стве базисных функций фу (х) и тестовых фj (у) выберем функции, предложенные в статье [4].
Пример расчетной сетки, на которой решалась матрица, представлен на рис. 3. Сначала нумеруются вертикальные ребра, затем горизонтальные.
Для того чтобы избежать возникновения ситуации, приводящей к делению на ноль, возможной в связи с особенностью в функции О (х, у), предлагается разнести внутренние узлы в квадратурной формуле. Например, для горизонтальных носителей возьмем точность интегрирования на носитель, равную п, а для вертикальных - (п + 2).
Первыми были получены результаты на квадратной пластине. Длина ребра квадратной пластины была равна длине волны X . Волновое число к равняется 271.
Для рис. 4 и 5 была рассмотрена задача с размером сетки 34 на 33 и точностью расчета интеграла 16 точки на носитель. Рисунок 4 отображает поведение поверхностных токов, направленных вдоль оси 0х. Рисунок 5 отображает поведение поверхностных токов, направленных вдоль оси 0У.
lij =- j G( x, y)div фі (x) • div ф j (y)ds + к2 J G( x, y )ф,- (x) • ф j (y)ds,
x - y|
Численные результаты
Рис. 3. Пример расчетной сетки
Решение задачи дифракции на экранах сложной геометрической формы субиерархическим методом
Описанный выше алгоритм позволяет получить решение задачи на экране канонической формы. Рассмотрим алгоритм построения решения задачи дифракции на экранах сложной геометрической формы [9-19]. Воспользуемся матрицей А, полученной методом Галеркина для прямоугольного экрана П. Для решения задачи на экране сложной геометрической формы О необходимо, чтобы новый экран полностью помещалось в экран п.
Введем вектор Щ, описывающий геометрию фигуры. Если носитель Пг- ( = 1,...,N) принадлежит экрану О, то Щ = 1, иначе Щ = 0. Решая
СЛАУ итерационным методом, будем поэлементно перемножать вектор геометрии на вектор решений. Таким образом, система будет решаться на некоторой подматрице Аг матрицы А. Подробное обоснование субиерархиче-ского метода приведено в [9].
Из уравнения А^ = / находим решение рассматриваемой задачи для
экрана сложной геометрической формы. На рис. 6-8 представлен пример расчета поверхностных токов для фигуры сложной геометрической формы. Примеры решения рассматриваемой задачи на различных экранах различными методами представлены в [6, 8, 10].
300
250
200
150
100
50
Рис. 4. Поверхностные токи их /
15 10
5 О
Рис. 5. Поверхностные токи /,/у /
Рис. 6. Геометрия экрана
Для экрана, представленного на рис. 9, были построены графики (рис. 10, 11). Фигура вырезалась из квадратной пластины размером ХхХ .
Было произведено теоретическое исследование задачи дифракции электромагнитного поля на экранах сложной формы. Получены численные результаты для решения уравнения электрического поля.
Рис. 7. Поверхностные токи
Рис. 8. Поверхностные токи
Рис. 9. Геометрия экрана
Рис. 10. Поверхностные токи
Рис. 11
Разработанные алгоритмы и программы могут быть использованы при
решении векторных задач электродинамики и при математическом моделировании сложных электродинамических процессов и объектов.
Список литературы
1. Самохин, A. Б. Интегральные уравнения и итерационные методы в электромагнитном рассеянии / A. Б. Самохин. - М. : Радио и связь, 1998. - 160 с.
2. Ильинский, А. С. Дифракция электромагнитных волн на проводящих тонких экранах / А. С. Ильинский, Ю. Г. Смирнов. - М. : Радиотехника, 1996. - 176 с.
3. Трибель, Х. Теория интерполяции, функциональные пространства, дифференциальные операторы / Х. Трибель. - М. : Мир, 1980. - 664 с.
4. Hanninen, I. Singularity subtraction integral formulae for surface integral equations with RWG, rooftop and hybrid basis functions / I. Hanninen, M. Taskinen and J. Sarvas // Prog. Electromagn. Res. PIER. - 2006. - V. 63. - P. 243-278.
5. Медведик, М. Ю. Применение функций крышек для решения задачи дифракции электромагнитных волн для решения задачи дифракции на экранах сложной формы / М. Ю. Медведик // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. 2012. - № 3 (23). - С. 84-98.
6. Rao, S. M. Electromagnetic Scattering by Surface of Arbitrary Share / S. M. Rao, D. R. Wilton and A. W. Glisson // IEEE Transactions on antennas and propagation. -1982. - V. Ap-30. - P. 409-417.
7. Смирнов, Ю. Г. О сходимости методов Галеркина для уравнений с операторами, эллиптическими на подпространствах, и о решении уравнения электриче-
ского поля I Ю. Г. Смирнов II Журнал вычислительной математики и математической физики. - 2007. - Т. 4, № 1. - С. 129-139.
S. Zwamborn, P. M. A Weak Form of the Conjugate Gradient FFT Method for Plate Problems I Peter M. Zwamborn and Peter M. van den Berg II IEEE Transactions on antennas and propagation. - 1991. - V. 39. - P. 224-22S.
9. Медведик, М. Ю. Применение субиерархического метода в задачах электродинамики I М. Ю. Медведик II Вычислительные методы и программирование. -2012. - Т. 13. - С. S7-97.
10. Медведик, М. Ю. Параллельный алгоритм расчета поверхностных токов в электромагнитной задаче дифракции на экране I М. Ю. Медведик, Ю. Г. Смирнов, С. И. Соболев II Вычислительные методы и программирование. - 2005. -Т. б. - С. 99-10S.
11. Медведик, М. Ю. Субиерархический параллельный вычислительный алгоритм и сходимость метода Галеркина в задачах дифракции электромагнитного поля на плоском экране I М. Ю. Медведик, Ю. Г. Смирнов II Известия высших учебных заведений. Поволжский регион. Сер. Естественные науки. - 2004. -№ 5. - С. 5-19.
12. Антонов, А. В. Разработка Web-ориентированного вычислительного комплекса для решения трехмерных векторных задач дифракции электромагнитных волн на основе субиерархических параллельных алгоритмов и ГРИД-технологий I А. В. Антонов, М. Ю. Медведик, Ю. Г. Смирнов II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2007. - № 4. -С. б0-б7.
13. Медведик, М. Ю. Субиерархический параллельный вычислительный алгоритм для решения задач дифракции электромагнитных волн на плоских экранах I М. Ю. Медведик, Ю. Г. Смирнов II Радиотехника и электроника. - 200S. - Т. 53, № 4. - С. 441-44б.
14. Медведик, М. Ю. Применение ГРИД-технологий для решения объемного сингулярного интегрального уравнения для задачи дифракции на диэлектрическом теле субиерархическим методом I М. Ю. Медведик, Ю. Г. Смирнов II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 200S. - № 2. - С. 2-14.
15. Медведик, М. Ю . Субиерархический метод для решения псевдодифференци-ального уравнения в задаче дифракции в слоях, связанных через отверстие I М. Ю. Медведик, И. А. Родионова II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2009. - № 3. - С. 59-70.
16. Медведик, М. Ю. Субиерархический подход для решения объемного сингулярного интегрального уравнения задачи дифракции на диэлектрическом теле в волноводе методом коллокации I М. Ю. Медведик, Д. А. Миронов, Ю. Г. Смирнов II Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2010. - № 2. - С. 32-43.
17. Медведик, М. Ю. Субиерархический метод решения задачи дифракции электромагнитных вол на диэлектрическом теле в прямоугольном волноводе I М. Ю. Медведик, Ю. Г. Смирнов II Радиотехника и электроника. - 2011. - Т. 5б, № S. - С. 940-945.
1S. Медведик, М. Ю. Субиерархический метод решения интегрального уравнения Липпмана - Швингера на телах сложной формы I М. Ю. Медведик II Радиотехника и электроника. - 20i2. - Т. 57, № 2. - С. i75-iS0.
19. Медведик, М. Ю. Субиерархический метод для решения псевдодифференци-ального уравнения в задаче дифракции в слоях, связанных через отверстие I М. Ю. Медведик, И. А. Родионова, Ю. Г. Смирнов II Радиотехника и электроника. -2012. - Т. 57, № 3. - С. 2S1-290.
Смирнов Юрий Геннадьевич
доктор физико-математических наук, профессор, заведующий кафедрой математики и суперкомпьютерного моделирования, Пензенский государственный университет
E-mail: [email protected]
Медведик Михаил Юрьевич
кандидат физико-математических наук, доцент, кафедра математики и суперкомпьютерного моделирования, Пензенский государственный университет
E-mail: [email protected]
Максимова Марина Александровна аспирант, Пензенский государственный университет
E-mail: [email protected]
Smirnov Yury Gennadyevich Doctor of physical and mathematical sciences, professor, head of sub-department of mathematics and supercomputer modeling, Penza State University
Medvedik Mikhail Yuryevich Candidate of physical and mathematical sciences, associate professor, sub-department of mathematics and supercomputer modeling,
Penza State University
Maximova Marina Alexandrovna
Postgraduate student,
Penza State University
УДК 517.6 Смирнов, Ю. Г.
Решение задачи дифракции электромагнитной волны на экранах сложной формы / Ю. Г. Смирнов, М. Ю. Медведик, М. А. Максимова // Известия высших учебных заведений. Поволжский регион. Физикоматематические науки. - 2012. - № 4 (24). - С. 59-72.
