Particle-Resolved simulations (PR-DNS) have been conducted using a second order implicit Viscous Penalty Method (VPM) to study the heat transfer between a set of particles and an incompressible carrier fluid. A Lagrange extrapolation coupled to a Taylor interpolation of a high order is utilized to the accurate estimate of heat transfer coefficients on an isolated sphere, a fixed Faced-Centered Cubic array of spheres, and a random pack of spheres. The simulated heat transfer coefficients are compared with success to various existing Nusselt laws of the literature.
