Loading [MathJax]/jax/output/CommonHTML/jax.js

Friday, September 26, 2014

Collisions And More Collisions

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(ρ)vn12{Nn(ρ)1}vn

v=12{Nn(ρ)+1}vn

vn=2Nn(ρ)+1v

and

vs=v2Nn(ρ)Nn(ρ)+1v={12Nn(ρ)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.