Abstract
A finite-difference method formulated in boundary-fitted nonorthogonal general curvilinear coordinates is used for the prediction of transient free convection in enclosures of arbitrary geometry. The governing system of equations in primitive variable form, simplified by the Boussinesq approximation, is transformed into a general curvilinear system in which all physical boundaries are described by constant coordinate lines. The finite-difference form of the equations is obtained by employing the control volume formulation and solved by a segregated numerical algorithm using a skewed time-marching procedure. Several test problems are considered in order to evaluate the performance of the method in terms of reliability, accuracy, stability, and efficiency. A problem of transient cooling of a square cavity is then considered in order to determine the effect of time step and under-relaxation coefficient on the accuracy and efficiency of the solution technique. The results of all test problems are presented and analyzed.