An Efficient Parallel Implementation of the SIMPLE Algorithm Based on a Multigrid Method

Abstract

This paper deals with a parallel implementation of the SIMPLE (Semi-Implicit Method for Pressure-Linked Equations) algorithm to numerically solve the Navier–Stokes system of equations for viscous incompressible flows. A mechanism of interprocess communication using a mesh decomposition with virtual cells and an algebraic multigrid method is proposed. A method of distributed matrix storage and an algorithm for matrix-vector operations reducing the number of interprocess communications are presented. The results of a series of numerical experiments on structured and unstructured grids (including a problem of external aerodynamics) are presented. Based on the results obtained, an analysis of the influence of multigrid solver settings on the overall efficiency of the algorithm is made. It is shown that the parallel algorithm based on an algebraic multigrid technique proposed for the SIMPLE method makes it possible to efficiently calculate problems on hundreds of processors.