Passive particle transport of 200 particles

This video shows a 2D passive particle transport of 200 particles in a channel with an obstacle. The motion is only given by the update position differential equation X^n+1 = X^n + v * dt, where n+1 is the future postion, n is the current position and v is the flow velocity, which does not depend on the particle's motion. As can be seen, particles in the middle of the channel have slow speed, while the others are faster. The numerical simulation was developed using the Finite Element method in an unstructured triangular mesh generated using the Gmsh software.