Количество частиц: 900 ( x 30 рядов).

Начальные перемещения частиц – нулевые.

Начальное распределение вектора скорости частиц – случайное, с плотностью:

Равномерная плотность внутри эллипса Vx2 / Vx,max2 + Vy2 / Vy,max2 = 1; 0 вне эллипса,
Нормальное распределение плотности.

Задаваемые начальные среднеквадратические значения σ (Vx(0)) и σ (Vy(0)) указаны под графиками.

Распределение частиц по компонентам скорости Vx (слева) и Vy (справа).

Оси абсцисс – |Vx| и |Vy| (масштаб одинаковый).
Оси абсцисс – sqrt|Vx| и sqrt|Vy| (масштаб одинаковый).

Оси ординат – плотности распределения частиц.

Изобразить нормальное распределение с дисперсией, равной дисперсии вычисляемой плотности.

P Чувствительность
к мгновенным изменениям
|Vx| |Vy|
σ (Vx(0)) = σ (Vy(0)) =


Средние значения Vx2 и Vy2 по всем частицам.
(Усреднение производится также по времени с весом, пропорциональным времени):

TxxS = < Vx2 > =
TyyS = < Vy2 > =
< Vx2 > + < Vy2 > =
< VxVy > =

Отношение моментов распределения компонент скорости:
Для нормального распределения < V4 > / < V2 >2 / 3 = 1

< Vx4 > / < Vx2 >2 / 3 =
< Vy4 > / < Vy2 >2 / 3 =

Отношение стационарных и начальных квадратов компонент скорости:

TxxS + TyyS = k1 (Txx0 + Tyy0)
TxxSTyyS = k2 (Txx0Tyy0)
k1 =
k2 =
(1/k2 = )

Изменение кинетической энергии во времени. Горизонтальная линия – половина полной энергии.


График целиком. Начало графика. Окончание графика.

Полная энергия, приходящаяся на одну частицу: E = .
Время t = ( шагов)

Скорости частиц Перемещения частиц
Vy Uy
Vx Ux