A parallel time-dependent incompressible flow solver and a
parallel multigrid elliptic kernel are described. The flow solver is
based on a second-order projection method applied to a staggered
finite-difference grid. The multigrid algorithms implemented in the
elliptic kernel, which is needed by the flow solver, are V-cycle and
full V-cycle schemes. A grid-partition strategy is used in the parallel
implementations of both the flow solver and the multigrid elliptic ker-
nel on all fine and coarse grids. Numerical experiments and parallel
performance tests show the parallel solver package is numerically
stable, physically robust and computationally efficient. Both the mul-
tigrid elliptic kernel and the flow solver scale very well to a large
number of processors on the Intel Paragon and the Cray T3D for
computations with moderate granularity. The solver package has
been carefully designed and coded so that it can be easily adapted
to solving a variety of interesting two and three-dimensional flow
problems. The solver package is portable to parallel systems that
support MPI, PVM and Intel NX for interprocessor communications.
