Метод дискретных вспомогательных источников использован для моделирования в резонансной частотной области электромагнитного рассеяния на трех импедансных телах. Приведены некоторые результаты моделирования, касающиеся влияния величины и типа поверхностного импеданса на сечения рассеяния одной из возможных структур, состоящих из трех эллипсоидов.
Simulating of electromagnetic scattering from three impedance bodies.pdf Изучение электромагнитных полей, рассеянных структурами, состоящими из нескольких тел, размеры которых сравнимы с длиной волны падающего на структуру поля, имеет большое значение для решения ряда практически важных проблем, например таких, как радиолокационная заметность, идентификация объектов, электромагнитная совместимость и др. Особый интерес представляет случай, когда расстояние между телами структуры много меньше длины волны. Корректная (с учётом электромагнитного взаимодействия тел) постановка исследований подобного рода приводит к необходимости решения граничных задач теории рассеяния на системах тел. Для задач рассматриваемого класса речь идёт о нахождении решений уравнений Максвелла, удовлетворяющих заданным граничным условиям на поверхностях тел и условиям излучения на бесконечности. В подавляющем большинстве случаев получить аналитическое решение таких задач не удаётся, поэтому используются различные численные методы. Например, в [1] для решения задачи электромагнитного рассеяния на двух и трёх диэлектрических сферах использован метод граничных элементов, а в [2] для решения подобной задачи использован метод интегральных уравнений. В последние годы применительно к решению задач электромагнитного рассеяния (в том числе и на группах тел) существенно развит метод дискретных источников. В этом методе неизвестное поле в рассматриваемой области и на её границах представляют в виде конечной линейной комбинации полей некоторой системы источников, размещённых вне этой области. Такая конструкция удовлетворяет системе уравнений Максвелла и условиям излучения (где это необходимо). Коэффициенты линейной комбинации определяются путём удовлетворения граничным условиям на поверхностях рассеивателей. В силу своей идейной про- Математическая постановка задачи имеет вид Vx Ee = m^eHe, Vx He = -тгеЁе в De; (1) nq X Ee - Zq (P>q X ("q X He)) =-"q X E0 + Zq (P>q X ("q X H0)) На Sq , q = 1,2,3; (2) ; 4^еЙе} x R / R + {; -^Ee} = O(R~>), R ^ ю, (3) 1 где nq - единичные векторы нормалей к поверхностям Sq; R = (x2 + y2 + z2)2; a x b - векторное произведение, Zq - поверхностный импеданс соответствующего тела, Re Zq > 0 (q = 1,2,3). 2. Модель рассеянного поля Модель рассеянного поля строится следующим образом. Введем (см. рис. 1) внутри каждого из рассеивателей Dq вспомогательную поверхность Seq = KeqSq, подобную поверхности рассеивателя Sq в смысле гомотетии с центром в точке Oq. Если поверхность Sq является центральной, центр гомотетии выбираем так, чтобы он совпадал с центром поверхности. Коэффициенты гомотетии (подобия) Ke q характеризуют удаление вспомогательных поверхностей от поверхностей соответствующих тел, их значения лежат в интервале 0 < Ke < 1 (при Keq = 0 вспомогательная поверхность стягивается в точку, при Keq = 1 она совпадает с поверхностью соответствующего тела). Выберем на каждой из вспомогательных поверхностей Se конечную сово- N купность точек {Mn }n=q1 и в каждой точке Mn разместим пару независимых вспомогательных элементарных электрических диполей с моментами pn,q = P.n,qeтn,q, Р„q = P-n,qeTn,q, q = 1,2,3 , ориентированными вдоль единичных направлений eTn,q, eTn,q , выбранных в плоскости, касательной к Seq в точке Mn q, и излучающих в однородную среду с параметрами ee и |e. Представим неизвестное рассеянное поле {Ee, He} в De в виде суммы полей введенных вспомогательных диполей 3 N 3 N Ee (M) = (i / ^ )£ £ Vx (Vx П ^), He (M) =£ £ Vxfl n,, q=1 n=1 q=1 n=1 ППn q = Te (M,Mn q)pn,q, Te (M,Mnq) = e^RuM^ )/(4nRMM^), (4) pn,q = pn,qsn,q + pn,qgn,q, M e De. 1 Т T2 T2 e Здесь ke = ^se|e - волновое число в среде De; RMM - расстояние от точки на Se q до точки M в De; T e (M, Mnq) - фундаментальное решение уравнения Гельмгольца в области De; p^, ppq (q = 1,2,3, n = 1, Nq) - неизвестные комплексные постоянные (дипольные моменты); Nq - число точек размещения диполей на вспомогательной поверхности Se . Поле (4) удовлетворяет уравнениям Максвелла (1) и условиям излучения (3) в области De. Для того чтобы удовлетворить граничным условиям (2), необходимо соответствующим образом выбрать значения дипольных моментов pTn,q, pTn,q (q = 1,2,3, n = 1, Nq). Используем для этого метод коллокаций. Пусть Mj (j = 1, Lq) - точки коллокации на поверхности Sq; Lq - число точек коллокации на Sq. Тогда для определения неизвестных р^, p^ (q = 1,2,3, n = 1, Nq) получим следующую систему линейных алгебраических уравнений с комплексной матрицей размером 2(Lj + L2 + L3) x 2(Nx + N2 + N3): nq x El,q - Zq X X Hj,q x)) = = -nq x Ejq + Zqn x (nq x Я0%x)), q = 1,2,3, j = iTq, (5) где nq, Ejq, He q и E0 q, H0' q - значения векторов нормали и компонент рассеянного (4) и возбуждающего полей в точке j на поверхности тела с номером q. Решение системы (5) определяем путём минимизации функции Ф = Ц| nq x (Ejq + E0,q ) - Zq (nj x (Ц x (H jq + H'^ ))) |2 . (6) q=1 j=1 После решения задачи минимизации (определения неизвестных дипольных моментов pTn,q, pTn,q (q = 1,2,3, n = 1, Nq)) необходимые характеристики рассеянного поля определяем из (4). В частности, для компонент рассеянного поля в дальней зоне имеем Eefi(M) = (^ /6e)1/2Hе,ф(M) = (exp(ikeR)/keR)De(9,ф) + 0(R~2), E^ (M) = -(Це / Se )1/2 He,e (M) = (exp(ikeR)/ keR) Dф (9, ф) + 0(R~2), (7) где компоненты диаграммы рассеяния De (9, ф) и D^, (9, ф) определены выражениями 3 Nq De(9,ф) = (iroke^e /4n)^^Gnq(9,ф)[(ео89cosфcosa^ + cos9sinфcospf'q q=1 n=1 - sin9 cos yin,q)p^ + (cos 9 cos фcos a^ + cos 9 sinфcos p^ - sin9 cos y!)pp^ ], Dф (9, ф) = ^ ± N±Onq (9, ф) x 4П q=1 n=1 x[(cos ф cos p^ - sin фcos a^)p^ + (cos фcos P!^ - sinфcos apq)pp^ ], Gnq (9, ф) = exp{-ike (sin9 cos ф^ + sin9 sinф>^ + cos 9zn,,q )}' (8) в которых cos ajn,q ,cos Pjn,q ,cos yjn,q и cos a2,q ,cos P2q ,cos у 2,q - направляющие косинусы единичных векторов eTn,q и eTn,q ; xnq, ynq, znq - декартовы координаты точки Mn ; 6 и ф - общепринятые угловые сферические координаты точки наблюдения M. Контроль точности модели (4) осуществляем путём вычисления относительного значения функции (6) на сетке точек, промежуточных по отношению к точкам коллокации, выбираемых на поверхностях Sq всех тел, входящих в систему: 3 L'q Д = (Ф' / Ф0) , Ф0 =ХЖ; X E0,q - Zq n X (nj X H J, q ))|2, (9) q=1 j=1 где Ф' - значение функции (6) на указанной выше совокупности точек; Ф0 - значение соответствующей нормы падающего поля на этой же совокупности точек; Lq' - число промежуточных точек на поверхности рассеивателя с номером q . 3. Численные результаты На основании изложенной выше модели создана программа для расчёта компонент рассеянного поля и контроля точности полученного решения. Предполагается, что тела, образующие рассеивающую структуру, являются трёхосными эллипсоидами. Входными величинами программы являются координаты центров эллипсоидов и ориентация осей эллипсоидов в глобальной системе координат, величины полуосей (в длинах волн), значения нормированных импедансов эллипсоидов Z'q(Z' = Zq /(|e /ee)12), возбуждающее поле {E0,H0}, параметры подобия Ke , числа точек размещения диполей Nq и точек коллокации Lq для каждого из эллипсоидов, образующих исследуемую структуру. Координаты точек размещения вспомогательных диполей и точек коллокации, а также направляющие косинусы касательных направлений, вдоль которых ориентируются диполи и ставятся граничные условия, первоначально вычисляются в локальной системе координат, связанной с соответствующим эллипсоидом, а затем осуществляется пересчёт этих величин для глобальной системы отсчёта. Минимизацию функции (6) осуществляем методом сопряжённых градиентов; итерационный процесс останавливается при условии, если относительное изменение функции на каждой из десяти последних итераций не превышает 0,001. При помощи данной программы выполнена серия вычислительных экспериментов, направленных на изучение влияния величины и типа поверхностного импеданса на сечения рассеяния различных структур. Ниже представлены результаты для одной из исследованных структур. Структура состоит из трех эллипсоидов. Центр первого эллипсоида совмещён с началом декартовой системы координат, центр второго эллипсоида размещен в точке с координатами (3,628, 0, 0), а центр третьего - в точке с координатами (-3,628,0,0), т.е. центры второго и третьего эллипсоидов расположены на оси x симметрично относительно центра первого эллипсоида. Полуоси эллипсоидов kea, keb, kec ориентированы вдоль осей x,y, z соответственно и принимают значения keax = 1,0, kebj = 1,5, kecx = 2,0 для первого эллипсоида и kea2 = kea3 = 2,0, keb2 = keb3 = 1,5, kec2 = kec3 = 1,0 для второго и третьего. Легко видеть, что наименьшее расстояние Д/ между точками поверхностей эллипсоидов равно 0,1 X, где X - длина волны возбуждающего поля. Эллипсоид, с которым связана система координат, считается идеально проводящим (Z' = 0), второй и третий эллипсоиды предполагаются импедансными. В качестве возбуждающего поля выбрана линейно поляризованная плоская волна. Предполагается, что волна падает на структуру таким образом, что вектор ke лежит в плоскости xz и образует с осью z угол у . Рис. 2. Исследованная структура При решении задачи параметры метода для каждого эллипсоида выбраны одинаковыми: Ке1 = Ke,2 = Ke,3 = 0,6, N1 = N2 = N3 = 168, Ц = L2 = L3 = 336 . В локальных системах координат с центрами в центрах эллипсоидов точки размещения диполей и точки коллокации распределены следующим образом. В каждом из четырнадцати полусечений ф = const, отстоящих друг от друга на угловое расстояние Дф = 25,7°, равномерно по углу 9 выбраны двенадцать точек размещения диполей. Для точек коллокации алгоритм их расположения по углу 9 выбран таким же, как и для точек размещения диполей, но выбирают их как в полусечениях ф = const, определённых для точек размещения диполей, так и посередине между ними. На рис. 3, 4 представлены бистатические сечения рассеяния в E -плоскости описанной выше структуры, возбуждаемой плоской волной, падающей под углом у = 0°, при различных значениях нормированного поверхностного импеданса второго и третьего эллипсоидов. E -плоскость - это плоскость, в которой лежат векторы E0 и ke возбуждающей волны; в сферической системе координат эта плоскость состоит из двух полуплоскостей: ф = 0° и ф = 180°. По оси абсцисс отложено значение угла 9 в градусах, по оси ординат - значение сечения рассеяния ст(9, ф) = lim 4nR2[| Eefi (9,ф) |2 +1E^ (9, ф) |2]/1 E0 |2, нормированное на квадрат длины волны и выраженное в децибелах. В выбранной форме представления результатов направлению прямого рассеяния соответствует угол 6 = 0°, направлению обратного рассеяния - угол 6 = 180°. g/А2, дБ - -10 - -20 - -30 - Рис. 3. Бистатические сечения рассеяния в E -плоскости при различных значениях нормированного комплексного индуктивного импеданса второго и третьего эллипсоидов 150 100 50 0 50 100 6, град 0 - -40 " аА2, дБ 0 - -10 - -20 - -30 - Рис. 4. Бистатические сечения рассеяния в E -плоскости при различных значениях нормированного реактивного емкостного импеданса второго и третьего эллипсоидов Рис. 3 относится к случаю, когда второй и третий эллипсоиды характеризуются комплексным индуктивным импедансом (величина реальной части импеданса определяет поглощение электромагнитной энергии рассеивателем). Кривая 1 - это бистатическое сечение рассеяния структуры, в которой второй и третий эллипсоиды также являются идеально проводящими ( Z'2 = Z3 = 0 ), кривые 2-4 - биста-тические сечения рассеяния структуры, в которой поверхностные импедансы Z2 и Z^ второго и третьего эллипсоидов равны и соответственно принимают значения 0,1 - 0,1/, 0,3 - 0,3/ и 0,5 - 0,5/ . Рис. 4 относится к случаю, когда второй и третий эллипсоиды исследуемой структуры характеризуются реактивным ёмкостным поверхностным импедансом. Кривые 1-4 на этих рисунках - это бистатиче-ские сечения рассеяния структуры, в которой поверхностные импедансы Z2 и Z3 второго и третьего эллипсоидов соответственно принимают значения 0, 0,1/, 0,3/ и 0,5i. Рис. 5 и 6 характеризуют зависимость сечения обратного рассеяния собр рассматриваемой структуры от угла у падения плоской волны. Рис. 5 относится к структуре, в которой импедансы второго и третьего эллипсоидов являются комплексными и индуктивными, а рис. 6 - к структуре, в которой импедансы второго и третьего эллипсоидов являются чисто реактивными и емкостными. Кривая 1 на рис. 5 - это сечения обратного рассеяния структуры, в которой все три эллипсоида являются идеально проводящими; кривые 2 и 3 - сечения обратного рассеяния структуры, в которой нормированные поверхностные импедансы Z2 и Z3 второго и третьего эллипсоидов равны и соответственно принимают значения 0,1 - 0,1/ и 0,5 - 0,5/ . Кривая 1 на рис. 6 - это также сечения обратного рассеяния структуры, в которой все три эллипсоида являются идеально проводящими, а кривые 2 и 3 -сечения обратного рассеяния этой же структуры, в которой нормированные поверхностные импедансы Z2 и Z3 также равны, но соответственно принимают значения 0,1/ и 0,5/. с/Х2, дБ 0 -10 -20 -30 0 50 100 150 у, град Рис. 5. Зависимость сечения обратного рассеяния от угла падения плоской волны при различных значениях комплексного индуктивного импеданса второго и третьего эллипсоидов На рис. 5 и 6 по оси абсцисс отложено значение угла падения у плоской волны в градусах, по оси ординат - значение сечения обратного рассеяния, нормированное на квадрат длины волны и выраженное в децибеллах. Анализ результатов, представленных на рис. 3 — 6, позволяет сделать следующие выводы. Замена в рассматриваемой структуре двух идеально проводящих рассеивателей на рассеиватели с теми же геометрическими параметрами, но характеризуемые комплексными индуктивными импедансами приводит к уменьшению сечения обратного рассеяния при любых рассмотренных углах падения плоской волны. Например, при падении плоской волны нормально линии, на которой расположены центры тел структуры, сечение обратного рассеяния уменьшается на 10 дБ. Этого нельзя сказать для ситуации, когда те же рассеиватели характеризуются чисто реактивными емкостными импедансами. В этом случае (см. рис. 6) для ряда углов падения волны имеет место увеличение сечения обратного рассеяния, а при падении волны нормально линии, на которой расположены центры тел структуры, сечение обратного рассеяния изменяется не более чем на 1 дБ. Заключение Таким образом, в данной работе методом дискретных вспомогательных источников построен численный алгоритм решения задачи электромагнитного рассеяния на трех импедансных телах. Алгоритм реализован как компьютерная программа для расчета компонент рассеянного поля. Приведены некоторые результаты численных расчетов, касающиеся влияния величины и типа поверхностного импедаса на сечения рассеяния одной из структур, состоящей из трех эллипсоидов.
Hall W.S., Mao X.Q. Boundary element method calculation for coherent electromagnetic scattering from two and three dielectric spheres // Engineering Analysis with Boundary Elements. 1995. V. 15. P. 313—320
Ewe W.-B., Ы L.-W., Leong M.-S. Solving mixed dielectric/conducting scattering problem using adaptive Integral method // Progress in Electromagnetic Research, PIER 46. 2004. P. 143-163.
Дмитренко А.Г., Уринов Р.И. Рассеяние электромагнитной волны на структуре из конечного числа трехмерных импедансных тел // Известия вузов. Радиофизика. 2012. Т. 55. № 4. С. 299-308.