Распространение радиально-ограниченных вихревых пучков в ближней зоне. Часть I. Алгоритмы расчёта

На примере дифракции плоской волны на круглой апертуре в ближней зоне (порядка нескольких длин волн) проведено сравнение алгоритмов расчёта с использованием векторного интегрального преобразования Рэлея-Зоммерфельда (РЗ) и разложения по плоским волнам (ПВ) по точности и скорости вычислений. В скалярном случае, что соответствует вычислению поперечных компонент электрического поля, результаты различаются только в очень близкой к апертуре области. В векторном случае при расчёте продольной компоненты в методе ПВ возникает особенность в области спектральных частот, радиус которых близок к единице. Предложены различные варианты обхода этой особенности. На расстоянии нескольких длин волн результаты двух рассматриваемых алгоритмов совпадают и отличаются от полученных с помощью конечно-разностного временного метода (FDTD) только масштабно (среднеквадратичное отклонение с учетом масштаба составляет менее 2%). Таким образом, рассмотренные в данной работе алгоритмы позволяют получать за существенно меньшее время структурно верную (но несколько завышенную по амплитуде) картину дифракции в ближней зоне. Такая амплитудная "завышенность" может быть связана с тем, что рассмотренные методы РЗ и ПВ не подразумевают наличия y-компоненты в дифракционной картине изначально x-поляризованного поля. Во второй части статьи рассматривается модификация метода ПВ, позволяющая учесть наличие всех векторных компонент.