Numerical method for computing nonlinear, time dependent, buoyant circulation of air in rooms