Jerry Brackbill (Los Alamos National Laboratory) and I developed a numerical method to solve equations for a liquid containing solid suspended particles . Suspensions are common, and occur in blood flow, transport of slurries, movements of sediments in river beds, and advanced composites manufacturing. The numerical approach is a modification of Peskin's  immersed boundary method, combined with features from a particle-in-cell code [3-4]. Equations of motion for the fluid are solved using a finite difference method on a computational grid, with suspended particles replaced by a force density in the fluid equations. The force density is computed by discretizing stress-strain constitutive equations for the solid, using data provided by a collection of material points which mark regions in the fluid occupied by solids. The method retains the advantages of the immersed-boundary method in that one set of equations holds in the entire computational domain and the complication of internal boundaries is eliminated. Also, computational work is linear in the number of material points, so the method has the potential to handle a large number of particles. Our new addition is the manner in which the force density is constructed, with material properties of the suspended solids clearly specified by stress-strain constitutive equations. Since the material points are not physically or logically connected, particles of any size and shape are easily treated. Unlike most other methods for simulating suspensions, this method is not restricted to Stokes' flow and, in principle, the suspending fluid can be non-Newtonian.
A two-dimensional implementation of the method has been used to establish validity of this formulation by comparing numerical results for elastic vibrations, and particle settling in viscous fluids, to theory and analysis . Part of this work was supported by a Sandia-University Research Project grant and an NSF grant.
Last updated: September, 1998
Copyright © 1998, Deborah Sulsky