Nonoverlapping domain decomposition applied to a computational fluid mechanics code

Abstract

The purpose of this paper is the description of the development and implementation of the linear part of a numerical algorithm for the simulation of a Newtonian fluid flow and the parallelization of that code on several computer architectures. The test problem treated is the steady state, laminar, incompressible, isothermic, 2D fluid flow (ex-tendible to 3D case), the Navier-Stokes equations being discretized by a fully coupled finite volume method. For this problem, sparse data structures, nonstationary iterative methods and several preconditioners are applied. The numerical results allow the conclusion that the fully coupled version can compete with the decoupled classic SIMPLE method (Semi-Implicit Method for Pressure-Linked Equations), by using the Krylov subspace methods. Parallel versions of the coupled method based on nonoverlapping domain decomposition are discussed. © Springer-Verlag Berlin Heidelberg 1999.