A two-dimensional (2D) lattice gas automata model was developed to describe the behavior of flow and heat transfer around a heated stationary circular cylinder. To validate the feasibility of this model, heat transfer of a heated stationary circular cylinder, which is concentrically placed in a square enclosure, was simulated and compared with analytical results. Further investigations on flow and heat transfer past a heated stationary cylinder in a channel at Reynolds number of 100–200 were then carried out. The flow and heat characteristics were presented by streamlines and thermal nephograms. The dependences of drag coefficient and Nusselt number on the Reynolds number were further studied and compared to literature data. Simulation results of the heat transfer in the enclosure agree with the solution of the diffusive heat transfer equation. For the investigated Reynolds number range, time averaged drag coefficient decreases as Reynolds number increases, while Nusselt number increases with increasing Reynolds number. Reasonable agreements with previous investigations were achieved, demonstrating that the presented model provides an alternate method to simulate flow and heat transfer problems.