Физический решатель оптики

Решатель физической оптики (PO) в Antenna Toolbox™ позволяет вам решать для RCS объекта.. В физической оптике зарегистрированный инцидент используется, чтобы вычислить токи на поверхность структуры в ответ на посягающую плоскую волну. С доступными токами можно получить рассеянное поле в желаемых точках в далеком поле.

Субдомен основные функции RWG и дополнительные размерности

Знакомые основные функции Вильтона Рао Глиссона (RWG) на треугольниках основаны [2].

В изображении для двух произвольных треугольных закрашенных фигур _trn+ и _trn-, имеющий области _An+ и _An- и совместно использующий общее ребро_{, ln} основные функции имеет форму

где вектор, чертивший от свободной вершины треугольного _trn+ к наблюдательному посту; вектор, чертивший от наблюдательного поста до свободной вершины треугольного _trn-. Основная функция является нулем вне двух смежных треугольников. Векторная основная функция RWG линейна и не имеет никакого потока (то есть, не имеет никакого нормального компонента) через его контур.

От [1], наряду со стандартным определением, этот метод требует двух модульных векторов нормали и 2D единичных векторов, также показанных на рисунке. Вектор является плоскостью треугольного _trn+; оба вектора перпендикулярны ребру _ln. Они заданы в центре ребра _ln, которые обозначаются. Направления

также показаны на рисунке. Этот метод принимает, что векторы нормали правильно (угол между смежным должен быть меньше 180 градусов), и исключительно заданный. Определенная векторная ориентация (например, внешние или внутренние векторы нормали) не имеет значения. Мы затем формируем два вектора векторного произведения,

и установите, что оба таких единичных вектора, направленные вдоль ребра, идентичны,

Только вектор в конечном счете необходим.

Поверхностная плотность тока, $\vec{J} (\vec{r})$ , на целой металлической поверхности расширен в _NPO основные функции RWG. эти основные функции (или области) могут наложиться и быть произвольно распределены на пробеле.

$\vec{J} (\vec{r}) = \sum_{n = 1}^{N_{P O}} I_{n}^{P O} {\vec{f}}_{n} (\vec{r}) (4)$

Нахождение _IPO

Подходящее Приближение ФО имеет форму:

$\vec{J} (\vec{r}) = 2 δ (\vec{r}) [\vec{n} (\vec{r}) \times \vec{H} (\vec{r})] (5)$

где δ составляет теневые эффекты. Если наблюдательный пост находится в затененной области, δ должен быть нулем. В противном случае это равняется ±1 в зависимости от направления падения относительно ориентации вектора нормали $\vec{n} (\vec{r})$ . Используя Eq. (4) урожаи:

$\sum_{n = 1}^{N_{P O}} I_{n}^{P O} {\vec{f}}_{n} (\vec{r}) = 2 δ (\vec{r}) [\vec{n} (\vec{r}) \times \vec{H} (\vec{r})] (6)$

Касательно [3] основы изящный способ выразить неизвестные $I_{n}^{P O}$ явным образом, с помощью интересного изменения метода коллокации. Во-первых, рассмотрите узел коллокации, который стремится к центру ребра ${\vec{r}}_{n}$ из центральной основной функции ${\vec{f}}_{n} (\vec{r})$ и расположен в плюс треугольник. Умножьте Eq. (6) вектором ${\vec{t}}_{n}^{+}$ . Поскольку нормальный компонент основной функции в ребре - все до одного, другие основные функции, совместно использующие тот же треугольник, не имеют никакого нормального компонента в ребре, результат становится

$I_{n}^{P O} = 2 δ ({\vec{r}}_{n}) {\vec{t}}_{n}^{+} \cdot [{\vec{n}}_{n}^{+} \times \vec{H} ({\vec{r}}_{n})] (7 a)$

Повторите ту же операцию с минус треугольник:

$I_{n}^{P O} = 2 δ ({\vec{r}}_{n}) {\vec{t}}_{n}^{-} \cdot [{\vec{n}}_{n}^{-} \times \vec{H} ({\vec{r}}_{n})] (7 b)$

Добавьте уравнения 7 (a) и 7 (b), разделите результат на два и преобразуйте тройное векторное произведение, чтобы получить:

$I_{n}^{P O} = \frac{2 δ ({\vec{r}}_{n}) \vec{H} ({\vec{r}}_{n}) \cdot ([{\vec{t}}_{n}^{+} \times {\vec{n}}_{n}^{+}] + [{\vec{t}}_{n}^{-} \times {\vec{n}}_{n}^{-}])}{} 2 (8)$

Поэтому согласно уравнениям (2) и (3),

$I_{n}^{P O} = 2 δ ({\vec{r}}_{n}) \vec{H} ({\vec{r}}_{n}) \cdot {\vec{l}}_{n} (9)$

Определение освещенных или теневых областей

Вычисление $δ (\bar{r})$ потребности составлять эффект затенения. Для простых выпуклых структур использование нормального, чтобы протестировать против направления излучения указало бы на освещенную или теневую область. Если нормальный из треугольника указывает в противоположном направлении излучения, то поверхность освещается. Если нормальный из треугольника находится в том же направлении, то поверхность затенена. Но этот простой тест перестал работать, когда объект невыпукл, как имеет место в более комплексных структурах. Чтобы обработать это, выполните треугольный сегментом тест на пересечение, чтобы строго определить значение $δ (\bar{r})$ . Значение $δ (\bar{r})$ 0 для теневых поверхностей или ±1 в зависимости от направления падения относительно ориентации вектора нормали. Чтобы реализовать это относительно основных функций RWG, которые формируются о поверхности Области ФО, проверяйте на обе произвольных треугольных закрашенных фигуры $t r_{n}^{+}$ и $t r_{n}^{-}$ быть в освещенной области и только затем считать вклад сделанным ребром к вычислению ПО текущий. Если любой треугольник находится в теневой области, значение дельты оценено, чтобы обнулить, и поэтому ребро не способствует.

Ссылки

[1] У. Джейкобус и Ф. М. Лэндсторфер, “Улучшенная Формулировка MM По для Рассеивания от 3D Отлично Проводящих Тел Произвольной Формы”, Сделка IEEE. Антенны и Распространение, издание AP-43, № 2, стр 162-169, февраль 1995.

[2] С. М. Рао, Д. Р. Вилтон и А. В. Глиссон, “Электромагнитное рассеивание поверхностями произвольной формы”, Сделка IEEE. Антенны и Распространение, издание AP-30, № 3, стр 409-418, май 1982.

[3] С. Макаров, антенна и моделирование EM в MATLAB, Вайли, Нью-Йорк, 2002.

Документация