Consider the particle of space, mps with velocity v in collision with its neighbors. This collision is elastic and we have,
conservation of momentum,
mpsv=Nn(ρ)mpsvn+mpsvs
where Nn(ρ) is the number of neighbors the particle collides with, a number that depends on density, ρ.
v=Nn(ρ)vn+vs ---- (1)
and conservation of kinetic energy,
12mpsv2=Nn(ρ)12mpsv2n+12mpsv2s
v2=Nn(ρ)v2n+v2s --- (2)
Squaring equation (1),
v2=(Nn(ρ)vn+vs)2=N2n(ρ)v2n+v2s+2Nn(ρ)vn.vs
minus equation (2)
0={N2n(ρ)−Nn(ρ)}v2n+2Nn(ρ)vn.vs
0=vnNn(ρ){(Nn(ρ)−1)vn+2vs}
Since vn and Nn(ρ) are not zero,
(Nn(ρ)−1)vn+2vs=0
−vs=12{Nn(ρ)−1}vn
From (1),
v=Nn(ρ)vn−12{Nn(ρ)−1}vn
v=12{Nn(ρ)+1}vn
vn=2Nn(ρ)+1v
and
vs=v−2Nn(ρ)Nn(ρ)+1v={1−2Nn(ρ)Nn(ρ)+1}v
vs=−Nn(ρ)−1Nn(ρ)+1v
when Nn(ρ)=1
vn=v and vs=0
The moving particle stopped and the once stationary particle moves forward with velocity v.
when Nn(ρ)=2
vn1=23v, vn1=23v and vs=−13v
when Nn(ρ)=3
vn1=12v, vn2=12v, vn3=12v and vs=−12v
The following diagram shows the cascade of collisions as time progresses,
The envelop of this cascade will always lead the shape of this velocity profile. Particle collision within the envelop will not gain enough velocity to surpass the envelop. As such we have a velocity curve as show below for the case of Nn(ρ)=2
The peak of this curve is reduced by subsequent collisions, each time by the collision factor,
Cc=2Nn(ρ)+1
For particles in the reverse direction, they are first reduced by the reverse reduction factor,
Rc=−Nn(ρ)−1Nn(ρ)+1
and then by Cc for each subsequent collision they encounter.
This velocity factor vs collision count curve is not the shape profile of the packet of energy. For that we need to integrate over time.