In order to investigate the flow characteristics in circular settling tanks with or without temperature difference, a k-epsilon two layer turbulence model was developed for this problem. In the numerical process, staggered grid system was employed, pressure and velocity were decoupled using SIMPLE method, and the hybrid finite analytic method was used to solve the governing equations. The simulation results of circular settling tanks with inlet at the center and outlet at the skirt show that the simulation results were in line with the results available in literatures. Without temperature difference and low suspended solids SS