An image of multiple clusters forming in Hagen-Poiseuille flow is shown below. We have shown the ability to simulate thousands of deformable particles and capture non-Newtonian characteristics that agree well with experimental observations. Our approach has shown the ability to scale on as many as 65,536 computational cores of Argonne National Laboratory’s IBM Blue Gene/P. When computing the dynamics of 1,000’s to 100,000’s of deformable particles suspended in fluid, datasets can easily reach 1 terabyte or greater for a single simulation. As an example, one such simulation with 100,000 deformable spheres (each with 254 triangular surfaces) requires a LB fluid domain in excess of 2 billion lattice grid points. When fluid and solid data sets from CFD codes are O(100 GB), the use of distributed memory architectures is required for post-processing.