Real-time or faster-than-real-time flow simulation is crucial for studying airflow and heat transfer in buildings, such as building design, building emergency management and building energy performance evaluation. Computational Fluid Dynamics (CFD) with Pressure Implicit with Splitting of Operator (PISO) or Semi-Implicit Method for Pressure Linked Equations (SIMPLE) algorithm is accurate but requires great computational resources. Fast Fluid Dynamics (FFD) can reduce the computational effort but generally lack prediction accuracy due to simplification. This study developed a fast computational method based on FFD in combination with the PISO algorithm. Boussinesq approximation is adopted for simulating buoyancy effect. The proposed solver is tested in a two-dimensional case and a three-dimensional case with experimental data. The predicted results have good agreement with the experimental results. In the two test cases, the proposed solver generates lower Root Mean Square Error (RMSE) compared to the FFD and at the same time, the proposed method reduces computational cost by a factor of 10 and 13 in the two cases compared to CFD.