Куранта–Фридрихса–Леви [2]:
Δ
t
6
1
s
1
Δ
x
2
+
1
Δ
y
2
=
Δ
x
c
√
2
.
(6)
На границах расчетной области вводятся граничные поглотители,
которые предотвращают нефизическое отражение электромагнитно-
го поля, выходящего из расчетной области, обратно в моделируемое
пространство. В данной работе используются граничные поглотители
Мура второго рода [4].
Для введения поля в расчетную область, а также для нормирования
параметров рассеянной электромагнитной волны воспользуемся мето-
дом полного и рассеянного полей (total field / scatter field (TF/SF)) [2].
В соответствии с указанным методом в расчетной области выделяется
замкнутый контур, охватывающий рассеиватель (рис. 2).
Внутри TF/SF контура существует как падающее на частицу поле,
так и поле, рассеянное частицей:
~E
total
=
~E
inc
+
~E
scat
;
~H
total
=
~H
inc
+
~H
scat
,
(7)
где
~E
inc
и
~H
inc
—
падающее поле,
~E
scat
и
~H
scat
—
рассеянное поле,
~E
total
и
~H
total
—
полное поле. Вне TF/SF контура присутствует только рассе-
янное поле.
Расчет полей осуществляется в двух модельных областях про-
странства одинаковых размеров, в одной из которых рассеиватель от-
сутствует, а в другой — присутствует. В пустой расчетной области зада-
ется источник электромагнитного поля, причем плоская волна вводит-
ся с помощью метода мягкого источника [4]. Электромагнитная волна
из вспомогательной расчетной области вводится через TF/SF границу
в пространство, содержащее рассеиватель. Рассчитывается полное по-
ле после дифракции волны на частице. Параллельно с вычислением
Рис. 2. Иллюстрация к методу полного и рассеянного поля [2]
ISSN 1812-3368. Вестник МГТУ им. Н.Э. Баумана. Сер. “Естественные науки”. 2012
49