Computational efficiency is often best served by maintaining the spatial proximity of a particle’s data to its location on the physical domain because the time advancement of particles is dependent upon the Eulerian velocity values in the overlapping mesh cells [4,5,14,1].