О точности моделирования интеграла кулоновских столкновений методом Монте-Карло

В предыдущих работах метод прямого моделирования типа Монте-Карло для кинетического уравнения с нелинейным оператором кулоновских столкновений подробно объясняется для случая многокомпонентной плазмы. Дана строгая оценка порядка аппроксимации интеграла столкновений Ландау–Фоккера–Планка (ЛФП) интегралом Больцмана. В данной публикации точность аппроксимации исследуется численно. Показано, что порядок аппроксимации выше, чем это удаëтся доказать аналитически. Анализируется возможность выбора оптимальных шагов по времени. Приводится сравнение с аналитическими формулами в задаче электрон-ионной температурной релаксации.