Немонотонность схемы FDTD при моделировании границ раздела между диэлектриками

Показано, что схема Йе метода FDTD является немонотонной. Продемонстрировано, что при моделировании диэлектрической слоистой среды фактический порядок точности вычислений понижается до первого. Для повышения точности предложено проводить расчеты на сгущающихся сетках по методу Ричардсона-Калиткина. В качестве тестовой задачи рассмотрено распространение гауссова импульса через кремниевую пластинку, помещенную в воздухе. Проведены вычисления с повышением порядка точности до четвертого.