CC BY
35
УДК 517.958, 537.876.46
М. А. Москалева
ИССЛЕДОВАНИЕ ЗАДАЧИ ДИФРАКЦИИ ЭЛЕКТРОМАГНИТНОЙ ВОЛНЫ НА СИСТЕМЕ ПЕРЕСЕКАЮЩИХСЯ ТЕЛ И ЭКРАНОВ1
Аннотация.
Актуальность и цели. Математическое моделирование процесса дифракции электромагнитных волн на экранах и телах различной формы является важным аспектом в современной электродинамике. Целью данной работы является исследование задачи дифракции электромагнитных волн на системе пересекающихся тел и экранов сложной формы численным методом.
Материалы и методы. Рассматривается метод решения задачи дифракции электромагнитной волны на системе тел и экранов сложных форм. Поставленная задача дифракции представлена в виде системы интегродифференциаль-ных уравнений. Полученная система решена проекционным методом. Представлены численные результаты.
Результаты. Визуализировано поведение отраженного поля, полученное математическим моделированием, от систем пересекающихся тел и экранов. Разработаны программа и алгоритм, позволяющие определять модули решения системы интегродифференциальных уравнений, к которым сведена задача дифракции электромагнитной волны.
Выводы. Разработанные программы и алгоритмы могут быть использованы при решении векторных задач электродинамики и при математическом моделировании сложных электродинамических процессов и объектов, например, при решении задач дифракции в резонансных диапазонах частот.
Ключевые слова: задача дифракции, система интегродифференциальных уравнений, метод Галеркина, сеточные базисные функции, обобщенная расчетная сетка.
M. A. Moskaleva
A NUMERICAL METHOD FOR SOLVING ELECTROMAGNETIC WAVE DIFFRACTION PROBLEMS ON THE SYSTEM OF CROSSING BODIES AND SCREENS
Abstract.
Background. Mathematical modeling of diffraction of electromagnetic waves on screens and bodies of various forms is an important aspect in modern electrodynamics. The goal of this work is to research a problem of electromagnetic wave diffraction on a system of crossing bodies and screens of irregular shapes by the numerical method.
Materials and methods. The method for solving the electromagnetic wave diffraction problem on the system of bodies and screens of irregular shapes was considered. The set problem was represented as a system of integro-differential equations. This system was solved by the projection method. The numerical results are presented.
1 Работа выполнена при финансовой поддержке Минобрнауки РФ (грант в рамках Госзадания № 2.1102.2014/^).
Results. The authors visualized the behavior of the incident field, resulting in mathematical modeling of the system of crossing bodies and screens. The developed program and algorithm allow to define a system of integro-differential equation solution modulus to which the electromagnetic wave diffraction problem is reduced.
Conclusions. The developed program and algorithm can be used for solution of vector problems of electrodynamics and mathematical modeling of electrodynamic process and objects, for example, diffraction problem solution in resonant frequency bands.
Key words: diffraction task, system of integro-differential equations, Galerkin method, grid basis function, generalized computational grid.
1. Постановка задачи
3
Пусть в пространстве R расположена система тел Qi и экранов Qj (i = 1,...,I; j = !,•••,J). Пример показан на рис. 1.
Рассмотрим задачу дифракции электромагнитной волны на данной системе тел и экранов. Обозначим объединение конечного числа экранов Пj
следующим образом:
П = У П у.
)
При этом Пj представляют собой связные ориентируемые незамкнутые и непересекающиеся ограниченные поверхности класса Св R^3.
Граница дПj := Пj \ Пj поверхности П j есть кусочно-гладкая кривая,
состоящая из конечного числа простых дуг класса Сбез точек самопересечения, сходящихся под углами, отличных от нулевого; дП := У дП. Предпо-
}
лагаем, что экраны являются идеально проводящими.
№ 1 (37), 2016 Физико-математические науки. Математика
Определим также трубчатые окрестности ЭП5 края экрана:
ЭП5 := {х е Ж3 : dist(x,ЭП) < 5}.
Предполагаем, что Q^ - ограниченные области, границы которых дQi = Qi \ Qi - кусочно-гладкие замкнутые ориентируемые поверхности,
состоящие из конечного числа поверхностей класса С1. Определим Q := . Будем полагать, что тела Qi имеют диэлектрическую проницаемость £(x).
В исследуемой системе тела и экраны пересекаются, т.е. Q Рассмотрим задачу дифракции электромагнитной волны Eo, Ho
с гармонической зависимостью от времени вида e-lШ на системе тел и экранов. Источником падающего поля может быть, например, ток jo E .
Ниже P+, P-
- произвольные области, внешняя и внутренняя соответственно к П и такие, что ПсЭР ±.
Требуется определить полное электромагнитное поле Щ, удовлетворяющее:
- уравнениям Максвелла:
rotH = -iweE + j0_g, rotE = iw^oH
(Об uQ)c); (1)
- условиям непрерывности касательных компонент на границе области неоднородности:
= [HX =0; (2)
- краевым условиям на поверхности экрана П (за исключением точек края экрана):
Eт|п =0; (3)
- условиям конечности энергии в любом ограниченном объеме пространства:
^Hе Ь2М(Ж3) = И°ос(Ж3); (4)
- условиям Сильвера - Мюллера:
Es, H5 = о(1/ г), 1т к >0
^ хeг -Es = о(1/г), Es хeг -^ = о(1/г), г (5)
Es, H 5 = О (1/ г), 1т к = 0,
для рассеянного поля
Е5 = E - E0, И, = H - ^ (г = |Х, ег = х /\х\).
Постановка задачи дифракции электромагнитной волны для случая системы непересекающихся тел и экранов и описание свойств уравнений электрического поля подробно представлено в [1].
2. Интегродифференциальные уравнения задачи дифракции
Перейдем к системе интегродифференциальных уравнений электрического поля сформулированной задачи дифракции. Согласно формулам из [2] получаем следующую систему:
J-(ko + graddiv) J G (x, y )J (y )dy-
Q
-(ko + graddivx) JG(x,y)u(y)dsy = Eo (x), xe Q;
Q
-(kg + graddiv) J G (x, y )J (y )dy -
Q
■(ko + graddiv T) J G (x, y )u (y )ds}
Q
= Eg т (x), x e Q .
(6)
H
где G (x, y ) =
1 exp (iko\x - y|)
- известная функция Грина.
4n |x - y|
Перепишем систему интегродифференциальных уравнений в оператор-
ном виде:
LV = f
(7)
здесь L =
Г L11
L
12 >
L21 L22
V = (J,u), f =
Г Eo ^
-Ei
Ниже описан численный метод решения этой системы уравнений.
3. Дискретизация задачи
Будем решать задачу методом Галеркина. Для этого построим равномерную расчетную сетку на системе тел и экранов. Пусть рассматриваемая система состоит из плоского экрана Q, имеющего прямоугольную форму, и тела Q , являющегося прямоугольным параллелепипедом. Разобьем экран и тело на элементарные ячейки - конечные элементы. Для плоского экрана конечные элементы представляют собой прямоугольники, а для тела - прямоугольные параллелепипеды.
Определим сеточные базисные функции на экране и теле. На экране введем функции rooftop по методу, предложенному в статье [3]. Данная функция 9i (x1, x2, x3), отвечающая ребру i, определяется по правилу
фу (хъ Х2, Х3) = <
(x1 - x1,i-b х2 - x2,i-Ъ х3 - x3,i-1) S+ в Pi
(x1,i+1 - x1,x2,i+1 - x2,x3,i+1 - x3)в Pi
Si
(8)
где Ц является длиной ребра i; £ + и £ есть площади Р + и Р соответственно (рис. 2).
Рис. 2
На теле введем базисные функции «крышки» [4], пусть :=1 Xj,k - хj,k-1 | (где у = 1,.., 3 ), тогда получим формулы для у]к1т :
Vklm (x) =
1 - V^J - xJ'k I, x enklm, 0, x ¿Щт.
Пример базисных функций на теле изображен на рис. 3.
(9)
Рис. 3
Применяя метод Галеркина к системе (6), приходим к решению матричного уравнения
Lu = f,
здесь u - столбец неизвестных коэффициентов при базисных функциях; L - основная матрица СЛАУ; f - столбец правой части.
Основная матрица имеет блочный вид:
L =
< L11 т 21
г 12 >
22
Блоки матрицы, расположенные по главной диагонали, соответствуют решению задачи дифракции только на теле и только на экране.
4. Численное решение. Результаты расчетов
Рассмотрим систему 0 (см. рис. 1), в которой экран О имеет прямоугольную форму размером ХхХ, а тело Q является прямоугольным
XXX „
параллелепипедом размером ^х 2х 2 . Экран расположен в плоскости
ОХ1Х2, Х3 = 0 . Центр тела совпадает с центром системы координат. Таким образом,
° = |Х еД3: ХЬ Х2 е^Х']; Хз =0
Q = ^ х е Д3: х^ е ^-X X- | \, i = 1,2,3.
Пусть падающее поле распространяется гармонически и направлено вдоль оси ОХ2 , диэлектрическая проницаемость тела постоянна и определяется параметром ^ = 3п. В этом случае результат решения задачи представлен на рис. 4-9. Рисунок 4 иллюстрирует распространение модулей поверхностных токов на экране вдоль оси 0x1.
Рис. 4. Распространение поверхностных токов на экране вдоль оси 0x1
В соответствии с теорией нормальная компонента поверхностных токов, расположенная вдоль оси 0x1, стремится к нулю, касательная компонента поверхностных токов, расположенная вдоль оси 0x2, неограниченно возрастает.
Рис. 5. Распространение поверхностных токов на экране вдоль оси Ох2
Рис. 6. Распространение поля внутри тела на первом слое расчетной сетки
"D 7bo4y3 .bet1
Рис. 7. Распространение поля внутри тела на пятом слое расчетной сетки
Рисунок 5 иллюстрирует распространение модулей поверхностных токов на экране вдоль оси Ох2 .
Рис. 8. Распространение поля внутри тела на одиннадцатом слое расчетной сетки
Рис. 9. Распространение поля внутри тела на пятнадцатом слое расчетной сетки
В соответствии с теорией нормальная компонента поверхностных токов, расположенная вдоль оси 0x2, стремится к нулю, касательная компонента поверхностных токов, расположенная вдоль оси 0x1, неограниченно возрастает.
Рисунок 6 иллюстрирует распространение поля внутри тела на первом слое расчетной сетки, расположенном перпендикулярно оси 0х\, при X X п
Х1 = —4 + 32 . Рисунок 7 иллюстрирует распространение поля внутри тела на
пятом слое расчетной сетки, расположенном перпендикулярно оси 0х\, при X 5Х
Х1 =---\--.
1 4 32
Рисунок 8 иллюстрирует распространение поля внутри тела на одиннадцатом слое расчетной сетки, расположенном перпендикулярно оси
X 11Х _
0x1 при Х1 =— 4 + —. Рисунок 9 иллюстрирует распространение поля
внутри тела на пятнадцатом слое расчетной сетки, расположенном
X 15Х
перпендикулярно оси 0x1 при х = ——
Наблюдается симметрия в характере распространения поля внутри тела на слоях расчетной сетки - первый слой симметричен пятнадцатому (последнему) слою, пятый - одиннадцатому и т.д.
Рассмотрим систему 0, в которой экран О имеет прямоугольную форму размером ХхХ, а тело Q является прямоугольным параллелепипедом
размером Хх2хХ . Экран расположен в плоскости Ох^ , Х3 =-—. Центр
тела совпадает с центром системы координат. Таким образом,
° = \ х ей3: Х1, х2 е ^ —^ Х|; хз =■— \,
Q = \ х е Я3: х, е ( — 2 X , / = 1,2,3. Система 0 представлена на рис. 10.
Рис. 10. Система 0
Пусть падающее поле распространяется гармонически и направлено вдоль оси 0х2, диэлектрическая проницаемость тела постоянна и определяется параметром ^ = 3п.
В этом случае результат решения задачи представлен на рис. 11-16. Рисунок 11 иллюстрирует распространение модулей поверхностных токов на экране вдоль оси 0х1. В соответствии с теорией нормальная компонента поверхностных токов, расположенная вдоль оси 0х1, стремится к нулю, касательная компонента поверхностных токов, расположенная вдоль оси 0х2, неограниченно возрастает.
"D:/5UHW)l_lGP,tiT
Рис. 11. Распространение поверхностных токов на экране вдоль оси 0x1
"DVstrtiniJGP.M"
Рис. 12. Распространение поверхностных токов на экране вдоль оси 0x2
'Di/bmlvl-W
Рис. 13. Распространение поля внутри тела на первом слое расчетной сетки
Рисунок 12 иллюстрирует распространение модулей поверхностных токов на экране вдоль оси 0x2 .
Рис. 14. Распространение поля внутри тела на пятом слое расчетной сетки
Рис. 15. Распространение поля внутри тела на одиннадцатом слое расчетной сетки
Рис. 16. Распространение поля внутри тела на пятнадцатом слое расчетной сетки
В соответствии с теорией нормальная компонента поверхностных токов, расположенная вдоль оси 0х2, стремится к нулю, касательная компо-
нента поверхностных токов, расположенная вдоль оси Oxi, неограниченно возрастает.
Рисунок 13 иллюстрирует распространение поля внутри тела на первом слое расчетной сетки, расположенном перпендикулярно оси Oxi, при XX
xi = — 4 + — . Рисунок 14 иллюстрирует распространение поля внутри тела на пятом слое расчетной сетки, расположенном перпендикулярно оси Ox1 ,
X 5X _
при xi = --4+ "32 . Рисунок 15 иллюстрирует распространение поля внутри
тела на одиннадцатом слое расчетной сетки, расположенном
перпендикулярно оси Oxi, при xi = —4 + ~32 . Рисунок 16 иллюстрирует
распространение поля внутри тела на пятнадцатом слое расчетной сетки,
X 15X
расположенном перпендикулярно оси 0xi, при xi = — 4 +
При помощи субиерархического метода [5] из данных систем тел и экранов исключили поочередно тело и экран. Результаты, полученные в ходе решения задачи, совпадают с результатами, полученными только на теле [6] и на экране [2].
Список литературы
1. Smirnov, Yu. G. Integrodifferential Equations of the Vector Problem of Electromagnetic Wave Diffraction by a System of Nonintersecting Screens and Inhomogene-ous Bodies / Yu. G. Smirnov, A. A. Tsupak // Advances in Mathematical Physics. -2015. - Vol. 2015. - 6 p.
2. Максимова, М. А. Численное решение задачи дифракции электромагнитных волн на системе тел и экранов / М. А. Максимова, М. Ю. Медведик, Ю. Г. Смирнов, А. А. Цупак // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. - 2014. - № 3 (31). - С. 114-133.
3. Hanninen, I. Singularity subtraction integral formulae for surface integral equations with RWG, rooftop and hybrid basis functions / I. Hanninen, M. Taskinen, and J. Sar-vas // Prog. Electromagn. Res. PIER. - 2006. - Vol. 63. - P. 243-278.
4. Смирнов, Ю. Г. Существование и единственность решения объемного сингулярного интегрального уравнения в задаче дифракции / Ю. Г. Смирнов, А. А. Цупак // Журнал вычислительной математики и математической физики. -2004. - Т. 44, № 12. - С. 2264-2274.
5. Медведик, М. Ю. Субиерархический параллельный вычислительный алгоритм для решения задач дифракции электромагнитных волн на плоских экранах / М. Ю. Медведик, Ю. Г. Смирнов // Радиотехника и электроника. - 2008. - № 4. -С. 441-446.
6. Medvedik, M. Yu. Volume Singular Integral Equation Method for Solving Problems of Diffraction of Electromagnetic Waves by a Dielectric Inhomogeneous Body in a Rectangular Waveguide / M. Yu. Medvedik, Yu. G. Smirnov // Forum for Electromagnetic Research Methods and Application Technologies (FERMAT). - 2014. -Vol. 2, № 5. - P. 1-5.
References
1. Smirnov Yu. G., Tsupak A. A. Advances in Mathematical Physics. 2015, vol. 2015, 6 p.
2. Maksimova M. A., Medvedik M. Yu., Smirnov Yu. G., Tsupak A. A. Izvestiya vysshikh uchebnykh zavedeniy. Povolzhskiy region. Fiziko-matematicheskie nauki [ University proceedings. Volga region. Physical and mathematical sciences]. 2014, no. 3 (31), pp. 114-133.
3. Hanninen I., Taskinen M. and Sarvas J. Prog. Electromagn. Res. PIER. 2006, vol. 63, pp. 243-278.
4. Smirnov Yu. G., Tsupak A. A. Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki [Journal of calculus mathematics and mathematical physics]. 2004, vol. 44, no. 12, pp. 2264-2274.
5. Medvedik M. Yu., Smirnov Yu. G. Radiotekhnika i elektronika [Radio engineering and electronics]. 2008, no. 4, pp. 441-446.
6. Medvedik M. Yu., Smirnov Yu. G. Forum for Electromagnetic Research Methods and Application Technologies (FERMAT). 2014, vol. 2, no. 5, pp. 1-5.
Москалева Марина Александровна аспирант, Пензенский государственный университет (Россия, г. Пенза, ул. Красная, 40)
E-mail: m.a.moskaleva1@gmail.com
Moskaleva Marina Aleksandrovna Postgraduate student, Penza State University (40 Krasnaya street, Penza, Russia)
УДК 517.3 Москалева, М. А.
Исследование задачи дифракции электромагнитной волны на системе пересекающихся тел и экранов / М. А. Москалева // Известия высших учебных заведений. Поволжский регион. Физико-математические науки. -2016. - № 1 (37). - С. 37-49.
