Combined forced and free convection around horizontal cylinders as well as spheres has been numerically investigated, with the emphasis placed on the pressure solutions. The stream function (ψ)-vorticity (ω) formulation is adopted in conjunction with hybrid (central and upwind) difference and nonuniform grids. For recovery of the pressure distribution, which seems unavailable in the literature, the pressure Poisson's equation is solved after convergence of ψ and at solutions has been achieved. Computed results of Nusselt numbers are compared with several reported literature values. Good agreement has been found between these values.