Finite difference calculations of buoyant convection in an enclosure. Part 1: The basic algorithm

A novel mathematical model of buoyant convection in an enclosure, developed earlier, is solved by finite difference techniques in the two-dimensional case. This model has been developed as a principal analytical tool for the prediction of the movement of smoke and hot gases in fires. Effects of large density variations caused by substantial heating are retained while acoustic (high-frequency) waves, which are unimportant to buoyant convection, are analytically filtered out. No viscous or thermal conduction effects are included in the model. These two characteristics (filtering and no dissipative effects) distinguish the model from all others describing buoyant convection. The mathematical model consists of a mixed hyperbolic and elliptic set of non-linear partial differential equations: the problem is a mixed initial, boundary value one. An explicit time-marching algorithm, second-order accurate in both space and time, is used to solve the equations. The computational procedure uses a software package for solving a nonseparable elliptic equation developed especially for this problem.