Доброго времени суток!
Как известно, чтобы моделировать не резонаторы в FDTD используются специального рода граничные условия, так называемые поглощающие. То есть мы не можем посчитать граничное значение магнитного или электрического поля т.к. необходимый для подсчета сосед в сетку не входит, в этом случае мы просто присваиваем сл элементу значение предыдущего. Например, если у нас размер 200, то для электрического поля пишем
![$Ez[0]=Ez[1]$ $Ez[0]=Ez[1]$](https://dxdy-03.korotkov.co.uk/f/e/f/2/ef290daabe78d98bc88c83220fd4d16882.png)
, аналогично для магнитного поля
![$Hy[199]=Hy[198]$ $Hy[199]=Hy[198]$](https://dxdy-02.korotkov.co.uk/f/9/a/0/9a0168caac8fe69ddd494c0754d5234182.png)
. Однако такой простой способ работает только для случая, когда за один временной шаг волна проходит один пространственный шаг, то есть

. Если мы добавляем проницаемость, то вместо поглощения мы наблюдаем отражение. Не очень ясно в чем природа этого отражения и как с ним бороться (кроме применения других условий)?