Software for Parallel Computation of 3D Thermal Convection

A computer program solves the differential equations of three-dimensional (3D) thermal convection of an incompressible fluid by use of a parallel-processing, finite volume numerical scheme. The equations of conservation of momentum, and energy are integrated over macroscopic control volumes on a normal, staggered grid. Upwind interpolation functions are used to prevent spurious numerical oscillations at high Rayleigh numbers. The resulting discretized equations, including a pressure equation that demands most of the computation time, are solved by a parallel-processing, multigrid method. The multigrid aspect of the method involves the use of a hierarchy of grids of different mesh sizes to obtain a solution on the finest grid. It has been proven, both theoretically and practically, that the multigrid aspect affords rapid convergence on a solution. The program could be used, for example, to predict oceanic convection currents. The effectiveness of the program has been demonstrated by applying it to test cases on several parallel-computing systems.