In the present article, radiation effect on mixed convection boundary layer flow of a viscoelastic fluid over a horizontal circular cylinder with constant heat flux has been numerically analyzed. The governing boundary layer equations are transformed to dimensionless nonlinear partial differential equations. The equations are solved numerically by using Keller-box method. The computed results are in excellent agreement with the previous studies. Skin friction coefficient and Nusselt number are emphasized specifically. These quantities are displayed against the curvature parameter. The effects of pertinent parameters involved in the problem namely effective Prandtl number and mixed convection parameter on skin friction coefficient and Nusselt number are shown through graphs and table. Boundary layer separation points are also calculated with and without radiation and a comparison is shown. The presence of radiation helps to decrease or increase the skin friction coefficient for the negative or positive values of the mixed convection parameter accordingly. The decrease in value of effective Prandtl number helps to increase the value of skin friction coefficient and Nusselt number for viscoelastic fluids
