CC BY
31
УДК 537.874.6
А.Г. Дмитренко, Р.И. Уринов
Решение задачи электромагнитного рассеяния на двух трехмерных импедансных телах методом дискретных источников
На основе метода дискретных источников построен алгоритм численного решения задачи электромагнитного рассеяния на двух трехмерных импедансных телах. Алгоритм реализован в виде компьютерной программы для расчета компонент рассеянного поля. Приведены некоторые результаты численных расчетов, иллюстрирующие влияние величины и типа поверхностного импеданса на бистатические сечения рассеяния.
Ключевые слова: электромагнитное рассеяние, импедансное тело, метод дискретных источников, сечение рассеяния.
Изучение электромагнитных полей, рассеянных структурами, состоящими из нескольких близко расположенных тел, размеры которых сравнимы с длиной волны, имеет большое значение для решения ряда практически важных проблем, например таких как влияние этих тел на характеристики антенн, идентификация объектов и др. Корректная (с учётом электромагнитного взаимодействия тел) постановка исследований подобного рода приводит к необходимости решения граничных задач теории рассеяния на системах тел. В подавляющем большинстве случаев получить аналитическое решение таких задач не удаётся, поэтому используются различные численные методы.
В данной работе предложенный ранее в [1] вариант метода дискретных источников для решения задач электромагнитного рассеяния на структурах из произвольного конечного числа трехмерных импедансных тел использован для численного решения задачи рассеяния на двух телах.
Геометрия структуры показана на рис. 1.
Будем рассматривать стационарную (зависимость от времени выбрана в виде ехр(-іюґ)) задачу дифракции электромагнитного поля {Ео,Ho} на структуре, состоящей из двух непересекающихся импедансных тел Dq, ограниченных гладкими поверхностями Sq ( =1,2). Эта структура размещена в однородной среде De с диэлектрической и магнитной проницаемостями 8е и Цє в декартовой системе координат Oxyz . Требуется найти
рассеянное поле Е,Не } в области De.
Математическая постановка задачи имеет следующий вид:
VxEe = №ЦеНе, УхНе = -№&еЕе в De, (1)
Пя х Ее - Zq (пч X (нч X Не )) = -Пч X Ео + 2Ч (^ х (пч х Но)) на Sq, (2)
{л/М^;Л/Й7Не}хR/R + } = O(R_1),^ ^ю, (3)
1
где ^ - единичные векторы нормалей к поверхностям Sq ; R = (x2 + у2 + z2)2; aх Ь - векторное произведение, Zq - поверхностный импеданс соответствующего тела, ЯеZq > 0(д = 1,2).
Ое (ве,ре)
{Ее, не}
{Ёо, Но} Рис. 1. Геометрия задачи
Суть метода решения задачи заключается в следующем. В соответствии с [1] введём (см. рис. 1) внутри каждого из рассеивателей Бд вспомогательную поверхность Бе,д = КедБд, подобную поверхности рассеивателя Бд в смысле гомотетии (подобия) с центром в точке Од . Коэффициенты гомотетии Ке д характеризуют удаление вспомогательных поверхностей от поверхностей соответствующих тел, их значения лежат в интервале 0 < Кед <1.
Выберем на каждой из вспомогательных поверхностей Бе,д конечную совокупность точек
[МП д }п=1 и в каждой точке МП д разместим пару независимых вспомогательных элементарных
^п,д п,д ^п,д ^п,д п,д ^п,д л Л электрических диполей с моментами р — р 1 е 1, р — р 1 е 1, д = 1,2, ориентированными
Т1 Т1 Т1 т2 т2 т2
V, -п,д -п,д г о
вдоль единичных направлений е^ , выбранных в плоскости, касательной к Бе,д в точке
Мп д, и излучающих в однородную среду с параметрами ее и це. Представим неизвестное рассеянное поле {Ее, Не } в Бе в виде суммы полей введенных вспомогательных диполей:
2 Кд 2
Ее (М) = (//Ш8е ) X ЕУХ(УхЙп,д ), Не (М) = £ ^УхППд ,
д=1п=1 д=1п=1
П пд = ^е (М Мп,д ) , ^ е (М Мп,д ) = ехр(/1е^ММп>д )/(4п&-ММ„а ), (4)
рпд = р”,де”,д + ртп’де”,д, М е Бе .
Здесь £е = ю^/ее|ае - волновое число в среде Ве ; Кммпд - расстояние от точки Мпд на 5(
е,д
до точки М в Бе; Те (М,МПд) - фундаментальное решение уравнения Гельмгольца в области
Бе; ртп,д,р^,д (д = 1,2,п = 1,Кд) - неизвестные комплексные постоянные (дипольные моменты); N - число точек размещения диполей на вспомогательной поверхности Бе д . Поле (4) удовлетво-
ряет уравнениям Максвелла (1) и условиям излучения (3) в области Бе. Для того чтобы удовлетво -рить граничным условиям (2), необходимо соответствующим образом выбрать значения дипольных
моментов рп,д,рп,д(д = 1,2,п =1,Ыд). Используем для этого метод коллокаций. Пусть М. (. = 1,Ьд) -
Т1 т2 4
точки коллокации на поверхности Бд; Ьд - число точек коллокации на Бд. Тогда для определения
п д п д
неизвестных р ,ч, р ,ч (д = 1,2, п = 1,Ыд) получим следующую систему линейных алгебраических
Т1 т2 4
уравнений с комплексной матрицей размером 2(Ь + Ь,2) х 2(#1 + N2):
пд х Ее,д -2д (пд х (пд х Не,д х)) = -Щ х Щ д + 1д (П1д х (П3д х Н0,д х)) д =1,2,. =1 Ьд, (5)
где пд,Ее,д,Не д и Е0д,Н0ц - значения векторов нормали и компонент рассеянного (4) и
возбуждающего полей в точке . на поверхности тела с номером д. Решение системы (5) определяем путём минимизации функции
2 Ьд . .....
Ф= X Щ х(Е.л +Е.д)-2д(п. х(п. х(Н.д +Н0,д)))|2 . (6)
д=1.=1
Задачу минимизации функции (6) решаем методом сопряжённых градиентов.
После решения задачи минимизации (определения неизвестных дипольных моментов
рп,д, рп,д (д = 1,2, п = 1,Ыд)) необходимые характеристики рассеянного поля определяем из (4).
Т1 т2 4
На основе изложенного выше подхода создана компьютерная программа, которая, в частности, использована для исследования влияния величины и типа поверхностного импеданса на бистатиче-ские сечения рассеяния. Некоторые результаты приведены на рис. 2.
ст / X2
Рис. 2. Бистатические сечения рассеяния в ^-плоскости структуры из двух эллипсоидов при различных значениях нормированного комплексного индуктивного импеданса второго эллипсоида. Кривая 1 - 22 = 0; кривая 2 - 22 = 0,1 -0,1/; кривая 3 - 22 = 0,3 - 0,3/; кривая 4 - 22 = 0,5 - 0,5/
Они относятся к структуре, состоящей из двух эллипсоидов с различными значениями полуосей. Центр первого эллипсоида совмещён с началом декартовой системы координат, центр второго эллипсоида размещен в точке с координатами (3,0628,0,0), т.е. на оси х. Полуоси эллипсоидов кеа, кеЪ, кес ориентированы вдоль осей х, у,2 соответственно и принимают значения кещ = 1,0, кеЪ = 1,5, кеС1 = 2,0 для первого эллипсоида и кеа2 = 2,0, кеЪ2 =1,5, кеС2 =1,0 для второго. Легко видеть, что наименьшее расстояние Д/ между точками поверхностей эллипсоидов равно 0,1Х, где X -длина волны возбуждающего поля. Эллипсоид, с которым связана система координат, считается
идеально проводящим (2' = 0), второй эллипсоид предполагается импедансным. В качестве возбуждающего поля выбрана линейно поляризованная плоская волна, распространяющаяся вдоль оси
2 , электрический вектор Ед которой ориентирован вдоль оси х. На рисунке представлены бистатические сечения рассеяния в Е -плоскости (плоскости векторов Е0 и ке ) при различных значениях
нормированного поверхностного импеданса 22 = 22/ второго эллипсоида. По оси абсцисс от-
V 8е
ложено значение угла 0 в градусах, по оси ординат - значение сечения рассеяния, нормированного на квадрат длины волны и выраженного в децибелах. В выбранной форме представления результатов направлению прямого рассеяния соответствует угол 0 = 0°, направлению обратного рассеяния -угол 0 = 180°. Параметры метода для каждого эллипсоида выбраны одинаковыми: *е,1 = Ке,2 = 0,6, N1 = N2 = 168, 11 = 12 = 336 .
Как видно из рисунка, замена в рассматриваемой структуре идеально проводящего рассеивателя на рассеиватель, характеризуемый комплексным индуктивным импедансом, приводит к заметному уменьшению сечения обратного рассеяния всей структуры. В рассмотренном диапазоне изменения поверхностного импеданса это уменьшение превышает 12 дБ. В то же время изменение комплексного индуктивного импеданса одного из рассеивателей почти не влияет на сечение прямого рассеяния. В рассмотренном диапазоне изменения импеданса сечение прямого рассеяния изменяется в пределах 0,5 дБ.
Литература
1. Дмитренко А.Г. Рассеяние электромагнитной волны на структуре из конечного числа трехмерных импедансных тел / А.Г. Дмитренко, Р.И. Уринов // Изв. вузов. Радиофизика. - 2012. - Т. 55, № 4. - С. 299-308.
Дмитренко Анатолий Григорьевич
Д-р физ.-мат. наук, профессор каф. исследования операций
Национального исследовательского Томского государственного университета (НИТГУ)
Тел.: (382-2) 41-31-69 Эл. почта: [email protected]
Уринов Радик Истамович
Аспирант каф. исследования операций НИТГУ
Тел.: +7-923-412-64-52
Эл. почта: [email protected]
Dmitrenko A.G., Urinov R.I.
Solution of electromagnetic scattering problem for two three-dimensional impedance bodies by discrete source method
We construct the algorithm of numerical solution of electromagnetic scattering for two three-dimensional impedance bodies on the basis of discrete source method. The algorithm is realized as a computer code. Some numerical results illustrating the influence of quantity and type of surface impedance are described.
Keywords: electromagnetic scattering, impedance body, discrete source method, scattering cross section.
