Количество частиц: 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 по всем частицам.
(Усреднение производится также по времени с весом, пропорциональным времени):
Отношение моментов распределения компонент скорости:
Для нормального распределения
< V4 > / < V2 >2 / 3 = 1
Отношение стационарных и начальных квадратов компонент скорости:
TxxS + TyyS = k1 (Txx0 + Tyy0)Изменение кинетической энергии во времени. Горизонтальная линия – половина полной энергии.
Полная энергия, приходящаяся на одну частицу: E = .
Время t = ( шагов)