Highly accurate schemes for 3D Maxwell equations with Lorentz media on the basis of alternate direction implicit time stepping algorithms

Предложен алгоритм расчета трехмерной нестационарной задачи дифракции, описываемой уравнениями Максвелла в лоренцевых средах. Алгоритм основан на неявной схеме переменных направлений с псевдоспектральной аппроксимацией пространственных производных. Порядок интегрирования по времени может составлять 2, 4, 6. Лоренцевская дисперсия учитывается путем введения вспомогательных неизвестных в исходную систему уравнений первого порядка. Вычислительная сложность алгоритма оценивается величиной $O(N^3)\log N$ операций на шаг по времени, где $N$ — число точек сетки по одному направлению.