Three-dimensional characteristics of fluid flow and heat transfer around a wavy circular cylinder having sinusoidal variation in cross sectional area along the spanwise direction are numerically investigated using the immersed boundary method. The three different wavelengths of ${pi}4$, ${pi}3$ and ${pi}2$ at the fixed wavy amplitude of 0.1 have been considered to investigate the effects of waviness especially on the forced convection heat transfer around a wavy cylinder when the Reynolds and Prandtl numbers are 300 and 0.71, respectively. The present computational results for a wavy cylinder are compared with those for a smooth cylinder. The time- and total surface-averaged Nusselt number for a wavy cylinder with ${lambda}={pi}/2$ is larger than that for a smooth cylinder, whereas that with ${lambda}={pi}/4$ and ${pi}/3$ is smaller than that for a smooth cylinder. However, because the surface area exposed to heat transfer for a wavy cylinder is larger than that for a smooth cylinder, the total heat transfer rate for a wavy cylinder with different wavelengths of ${lambda}={pi}/4$, ${pi}/3$ and ${pi}/2$ is larger than that for a smooth cylinder.