Fully resolved simulation in a Cartesian computational grid. Immersed boundary method is used to resolve particle and cylindrical column wall. Conjugate heat transfer between fluid and solid particles are computed numerically. Pressure drop and Nusselt correlations for wall-to-bed heat transfer are proposed.