Saturated fluid convective heat transfer in ordered three-dimensional porous media has been investigated numerically using a CFD software. Spherical particles with diameters dp equal to respectively 14, 9.4, 7 mm have orderly arrangements constituting the porous media skeleton. There is a copper plate of constant heat flux density above the porous skeleton. The distribution of temperature and local convective heat transfer coefficient in the channel, as well as factors affecting the convective heat transfer are investigated by the method of fluid-solid coupling. Results show that the thickness of the thermal penetration and the thermal boundary layer increase in the flow direction, and decrease with the flow rate. When the thermal conductivity of the skeleton is relatively high, the convective heat transfer increases slightly with decreasing particle diameter. The convective heat transfer coefficient decreases with the increase in mass concentration of polyacrylamide solution, because the viscous dissipation weakens convection heat transfer.