A Parallel Incompressible Flow Solver Package with a Parallel Multigrid Elliptic Kernel

The development and applications of a parallel, time-dependent incompressible Navier?Stokes 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 inplemented in the parallel 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 kernel 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 multigrid elliptic kernel and the flow solver scale very well to a large number of processors on Intel Paragon and 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 NX for interprocessor communications.