In the present study, the turbulent flow around a square cylinder with uniform injection or suction through various surfaces has been numerically studied. The wall functions of the viscous sublayer and log-law layer for velocity and temperature on the cylinder have been modified to allow for the effect of the injection or suction through the cylinder. The finite-volume method has been used for solving the ensemble averaged Navier-Stokes equations for incompressible flow in conjunction with the k- turbulence model equations that include the Kato-Launder's modification. A nonuniform staggered grid arrangement has been used. The results show that the parameters taken into account such as drag coefficient, and Strouhal and Nusselt numbers have been influenced by the uniform injection and suction through certain surfaces of the square cylinder, as well as the flow field parameters.