Appeared in FIFTH INTL. SYMP. ON SOLVING IRREGULAR STRUCTURED PROBLEMS IN PARALLEL: IRREGULAR '98, LNCS No. 1457, pp. 226-237, Springer, August 1998

Parallel Simulation of Particulate Flows

Matthew G. Knepley - Vivek Sarin - Ahmed H. Sameh
 
Department of Computer Sciences,
Purdue University,
West Lafayette, IN 47907
{knepley,sarin,sameh}@cs.purdue.edu
Phone: (765)494-7816, FAX: (765)494-0739

Abstract:

Simulation of particles in fluids requires the solution of nonlinear Navier-Stokes equations for fluids coupled with Newton's equations for particle dynamics, in which the most time consuming part is the solution of nonsymmetric and indefinite sparse linear systems. In this paper, we present a comprehensive algorithm for the simulation of particulate flows in two dimensional domains. A backward Euler method is used for time evolution, and a variant of Newton's method is used to solve the nonlinear systems. The linear systems are solved efficiently by a novel multilevel algorithm that generates discrete divergence-free space for the incompressible fluid. Unlike incomplete factorization preconditioners, our technique has the desirable properties of robust and effective preconditioning along with efficient implementation on parallel computers. We present experiments on the SGI Origin2000 that demonstrate the parallel performance of our algorithm, and discuss various aspects of the simulation package and the associated software design.


| AEM Home | Institute of Technology |

| Academics | Research | People | Information | Contact AEM |
| Major Projects | by Specialty | Other Programs |

Solid-Liquid Flows Home | Publications

Last updated October 16, 2000