Gatheri, F. K.2015-06-112015-06-112005East African Journal of Physical Sciences 6(1): 33-48 (20051729-5106http://ir-library.ku.ac.ke/handle/123456789/12896Research ArticleA mesh generation procedure for the solution of coupled, non-linear, partial differential equation is described. The essence of the procedure is the use of a grid generating function, which results in a strong refinement along the walls. In viscous fluid flow, the problem arises where the solution varies rapidly over a small part of the domain but over the rest of the domain changes very slowly. At very large Reynolds number, viscous fluid flow pattern changes rapidly in narrow boundary layers close to walls where the fluid is brought to rest. In such convection dominated flows it is desirable to concentrate more mesh points in certain regions of the cavity. This is in order to place more mesh points within areas of steep velocity gradients. The effectiveness of the procedure is demonstrated by applying it to the problem of turbulent natural convection in an enclosure with colliding boundary layers.enGrid generatorNatural convectionTurbulent flowsArticle