Fluid flow and heat transfer of mixed convection from a constant wall temperature circular cylinder in zero-mean velocity oscillating cooling flows have been simulated based on the projection method with two dimensional exponential stretched staggered cylindrical meshes. Cycle mean temperature and secondary streaming are obtained by the method of partial sums of the Fourier series. Present numerical results are validated by comparing the heat transfer results of free convection and the secondary streaming of pure oscillating flow over a circular cylinder to published experimental and numerical results. The complete structures of the cycle mean temperature and secondary streaming patterns are provided by numerical simulations over wide ranges of the Reynolds number, the Keulegan-Carpenter number and the Richardson number. Based on turning points of the curves of the overall Nusselt numbers versus Reynolds numbers and the characteristics of the cycle averaged temperature and flow patterns, the heat transfer can be divided into three linear regimes (conduction, laminar convection, and turbulent convection dominated regimes) and two non-linear transition regimes. The effects of wave directions, amplitudes, frequencies, and buoyancy forces on the enhancement of heat transfer are also investigated. The effective ranges of the governing parameters for heat transfer enhancement are identified. © 2013 Elsevier Inc.