This report describes a numerical method based on a finite difference technique for simulating indoor airflows in a building using a JC-e turbulence method. The model treats three dimensional non-isothermal turbulent flows using the Boussinesq approximation for buoyancy. It solves the resulting nonlinear system of momentum, energy and turbulence equations by an explicit time marching technique to obtain a solution for either a steady state or transient flow.