$Cast AIv
Электростатическая задача, для решения уравнений движения частиц использую leapfrog, и с ним всё нормально, но в сильно неоднородных полях, при переходе частицы из одной ячейки в другую, постепенно накапливаются ошибки.
попробовал, ради интереса, интегрировать не leapfrog'ом, а полуаналитическим методом, когда полная энергия сохраняется с машинной точностью, а момент перехода чатицы из одного узла в другой фиксируется и, соответсвенно применяются уже новые поля. Проверил, как оно решает — всё отлично, какое бы неоднородное поле не было, фазовая траектория всегда представляет собой замкнутую линию. Вставил этот интегратор в программу, но получил лишь нефизичные усиливающиеся ленгмюровские колебания, при том что шаг времени и сетки выбран такой же как раньше. Время дискретизируется с удвоенной плазменной частотой, а сетка с шагом дебая. В чём дело и почему нельзя применять полуаналитическое описание? В книжках и статьях ничего на эту тему не нашёл.
