Градиент FDTD ошибается, хотя S11 уже выглядит сошедшимся
У резонатора есть длинная память. После короткого возбуждающего импульса энергия в нём уходит экспоненциально, а FDTD-решатель обрывают в момент, который кажется безопасным: спектр S11 почти застыл, портовый сигнал затих, достигнут заданный порог энергии. Для обычной частотной характеристики этого часто достаточно. Для вопроса «на сколько сдвинется S11, если чуть передвинуть стенку, патч или диэлектрик?» тот же расчёт может быть ещё далёк от ответа.
У производной есть лишний временной хвост
В arXiv:2610.28496 автор разбирает этот эффект для случая одной доминирующей затухающей моды, где изменение геометрии в основном сдвигает полюс резонанса. Ошибка резонансной части S-параметра после моделирования длительностью T спадает как exp(−T/τ). У ошибки его чувствительности появляется множитель (1+T/τ): производная отклика содержит временной фактор, поскольку небольшая перестройка резонансной частоты всё сильнее набегает по фазе на длинном хвосте колебаний.
Из этой разницы получается неприятная арифметика. Для 1% ошибки самого отклика в простой модели нужно около 4.6 постоянной затухания τ. Градиенту требуется около 6.6τ. Разница выглядит умеренной, пока симуляция не упирается в высокодобротную моду, для которой одна τ уже означает много тысяч временных шагов.
Гладкий S11 скрывает большую ошибку
В PEC-полости с модой TM110 и нагруженной добротностью QL=5025 автор сравнил расчёт с практически сходящимся длительным прогоном. При остановке на уровнях энергии от −30 до −60 dB ошибка чувствительности составляла 72.3%, 42.9%, 22.1% и 10.3%. Отклонение |S11| в этих же точках лежало между 0.50 и 0.03 dB. Это не доказательство полной сходимости комплексного S11: в данной структуре слабосвязанная резонансная часть была невелика относительно полного сигнала. Именно поэтому знакомый гладкий график может особенно убедительно скрывать неверный градиент.
Проверка градиента совпадением нескольких методов тоже не спасает автоматически. Automatic differentiation, дискретный adjoint и конечные разности, применённые к одному усечённому временному окну, дифференцируют один и тот же неполный DFT. Они способны согласиться между собой и коллективно ошибаться. Автор демонстрирует близкое к модели поведение в rfx, Meep и openEMS, однако это пока результаты одного препринта без независимой репликации и открытого набора моделей.
Короткий прогон можно продолжить аналитически
Работа предлагает вытащить хвост резонанса из короткого расчёта методом matrix pencil. Несколько зондов поля внутри объёма записывают затухающие колебания; алгоритм оценивает их комплексные полюса и продолжает отклик аналитически. В сообщённом примере той же полости трёх пробников и 26.3 ns моделирования хватило для ошибки чувствительности менее 1% вблизи TM110, тогда как прямому DFT потребовалось 1.6–1.8 μs. Заявленное 66-кратное сокращение относится именно к этому случаю, а не ко всем FDTD-моделям.
Формула перестаёт быть точной при нескольких близких модах, заметной зависимости residue от параметра, расстройке от резонанса, низкой добротности и длинном импульсе. Но инженерный вывод остаётся жёстким: длительность FDTD нужно выбирать отдельно для целевой величины. Если симуляция должна двигать геометрию по градиенту, критерием окончания обязан стать сам градиент, его экстраполяция или оценка модового хвоста.