This method consists of a discretization of the momentum equations to obtain the velocities at the faces of a finite volume, in terms of the values of these variables at the grid points followed by the integration of the momentum and continuity equations in the finite volumes. This integration leads to equations where the values of the variables at the cell faces are to be replaced by the expressions obtained at the previous stage.

The linear system to be solved at each nonlinear iteration connects values of velocities and pressure at each grid point in each equation. The coefficient matrix is large, nonsymmetric, sparse, with non-null entries on the diagonal. The characteristics of these linear systems indicate the use of nonstationary iterative methods, for instance preconditioned GMRES, for their solution.

The strategy used to parallelize this method was nonoverlapping domain decomposition. The linear system to be solved at each outer iteration is then equivalent to the solution, on the interfaces separating subdomains, of the Schur complement reduced system followed by the update of the solution in the subdomains.

Numerical tests will be reported showing the number of outer iterations and elapsed times of this code on a Transputer based machine, compared to the corresponding results on a cluster of workstations.

purchase the full-text of this chapter (price £20)

go to the previous chapter
go to the next chapter
return to the table of contents
return to the book description
purchase this book (price £95 +P&P)