Parallel computation of natural convection in trapezoidal porous enclosures

In this study, coupled non-linear partial differential equations governing the natural convection from an isothermal wall of a trapezoidal porous enclosure have been solved numerically by finite element method (FEM) in conjunction with GMRES, a Krylov subspace based solver. In view of the enormous amount of computation, a parallel numerical algorithm for incomplete LU-conjugate gradient (ILU-CG) solver on eight-noded ANUPAM cluster under MIMD paradigm based on ANULIB message passing library has been developed. Parallel computations have been carried out for various values of flow and geometric parameters both under Darcian and non-Darcian assumptions on the porous model. Cumulative heat fluxes and Nusselt number (Nu) associated with convection process are presented through computer generated plots. (C) 2004 IMACS. Published by Elsevier B.V. All rights reserved.