есть некоторый маленький объём, заполненный односортными частицами, в пределах объёма плотность считается постоянной. требуется расчитать диффузию в пространстве скоростей.
во время расчёта идёт итеративный процесс, с шагом по времени dt, характерное время между столкновениями dt_{st}. на каждой итерации по заданной вероятности парных столкновений разыгрывается столкнулась ли частица или нет. после этого номера столкнувшихся частиц объединяются в пары и рассчитываются сами столкновения — считаем, что энергия распределяется равновероятно между частицами, а углы отклонения определяются из законов сохранения.
проблема — если шаг моделирования сравним со временем столкновений, то очень быстро функция распределения частиц по энергиям приобретает максвелловскую форму (http://ompldr.org/vZDFiMA/screen-maxwell.png, по x — энергия, по y — кол-во частиц), но в реальной задаче шаг моделирования много меньше шага столкновений (т.е. столкновения редки) и при моделировании с тем же _полным_ временем вместо максвелловского распределения получается что то вроде гаммараспределения http://ompldr.org/vZDFiMQ/screen-nomaxwell.png (но средняя энергия и импульс системы нормально сохраняются)
в чём здесь ошибка?
