In our earlier papers we have presented a general algorithm for the solution of both compressible and incompressible Navier‐Stokes equations. The objective of the present work is to show the performance of this algorithm when it is used to solve thermal flow problems. Both natural and forced convection and transient problems are considered in this study. The semi‐implicit form of the algorithm has been used to deal with a variety of these problems.