The two-way coupling of temperature and flow fields in the indoor environment makes the conventional CFD simulation computationally intensive, especially when the transient evolution of temperature and velocity fields is of interest in areas such as unsteady ventilation. This study, by applying different splitting CFD methods (fast fluid dynamics (FFD) and fractional step (FS) method) to solve the momentum equation to decouple velocity and pressure fields, further explored different strategies to tackle the coupled temperature equation, and validated the different sequences of operation based on different combinations of temperature/velocity solution schemes. Each operational sequence of different combined numerical schemes was validated on a mixed convection airflow benchmark case with a heated floor. The results showed that different operational sequences using the numerical schemes which considered the strong coupling effect of temperature and flow fields could predict the final steady-state result well, while for the transient simulation that need to track the physical flow/thermal fields variation process precisely, the operational sequence which introduce the effect of buoyancy into volumetric flux to solve pressure Poisson equation showed better accuracy and owned the capability of predicting the same transient evolution process as conventional CFD method. Based on the comprehensive comparison, several promising numerical schemes under proper operational sequence were recommended and the highest computing efficiency was provided by sequence D of RIPC + Iterative which saved 68% computing time compared to the conventional CFD method.