A High-order Method for Simulating Flows With Drag-coupled Particles

Jay Stotsky, University of Colorado

Fluid flows with embedded small particles often arise in the simulation of colloids, plasmas and biological fluids. The presence of small particles in a fluid flow induces singular point forces to act on the fluid. Without careful treatment, these forces are not accurately accounted for in standard numerical algorithms for solving fluid motion equations. In this poster, two efficient algorithms for obtaining solutions to fluid flows with large numbers of suspended particles are presented. One method is based on Particle-Particle-Particle-Mesh (PPPM) methods. The second is based on the Method of Local Corrections (MLC). We will analyze the error properties of these methods through analytical derivations and numerical simulation. Our results show that the PPPM-based method can be implemented efficiently, but is limited to second-order convergence, whereas the MLC-based algorithm, when combined with a special class of stencils, can yield higher-order accuracy.

Abstract Author(s): J. Stotsky