Lawrence Berkeley National Laboratory

We present a variation of an adaptive projection method for computing solutions to the incompressible Navier-Stokes equations with suspended particles. To compute the divergence-free component of the momentum forcing due to the particle drag, we employ an approach which exploits the locality and smoothness of the Laplacian of the projection operator applied to the discretized particle drag force. We present convergence and performance results to demonstrate the effectiveness of this approach.