Implementations of multi-fluid hydrodynamic simulations on distributed memory computer with a fully parallelizable preconditioned Bi-CGSTAB method

This paper reports recent work in implementing a multi-phase fluid dynamic code on a super computer system with distributed memory architecture. In th e numerical model. the dynamics of different fluids is formulated based on a unified framework for general continuum, The computational code was const ructed so that it is easily implemented on a distributed memory machine. A projection procedure for solving the flows of various compressibilities res ults in a pressure Poisson equation which has large discontinuities and sin gularities in the coefficients and the right-hand side terms, A preconditio ned Bi-CGSTAB method is used to solve the Poisson equation. A tridiagonal f actorization preconditioner is adopted to utilize the large scale pipeline structure of the hardware, To get perfect parallelism, we made use of a Jac obi splitting to invert the tri-diagonal system in the partitioning directi on. The resulting preconditioner is tested with a problem that has large ju mps and singularities. A satisfying vector processor utilization and a para llel speed-up of the fluid code are achieved. (C) 2001 Elsevier Science B,V , All rights reserved.