The aim of the present research and development work is to develop the computer program to
simulate the steady two-dimensional turbulent flow. The finite volume method is used to numerically
solve the flow governing equations. The Navier-Stokes equations are solved for the velocity field and the
SIMPLE algorithm is used to adjust the velocity field to satisfy the conservation law of mass. Since all the
variables are stored at the center of each control volume, the Rhie-Chow interpolation is used to avoid the
decoupling between the velocity and the pressure. The correct velocity field is then used to solve the k −
and ε − equations. The eddy viscosity, that represents the influence of turbulence on the mean flow field,
can be calculated from those values of k and ε obtained. The boundary layer on a flat plate is employed
as a test case because it is one of the standard benchmark problems for the validation of CFD software.
The sequential-computing solver is first used to obtain the computed results. It is found that the computed
results are in good agreement with the DNS data available. The parallel solver is also implemented here
and tested against the sequential-computing one. It is found that the parallel program can run faster than
the sequential one up to 7.37 times for the best case. Furthermore, the governing equations are solved on
the structured and body-fitted coordinates so that this computer program can be developed further for the
simulation of turbulent flow over or inside any object of complex geometry in the future.