The calculation of heat transfer between fluid and a set of particles is a complex task as soon as the simulation mesh is not adapted to the shape of the particles. In the present work, a, Implicit Tensorial Viscous penalty method is used [1] to resolve the flow. A new Lagrange extrapolation coupled to a Taylor interpolation of high order [3] is extended for estimating heat transfer coefficients on isolated sphere, a fixed Faced-Centered Cubic arrays of spheres and a random packed spheres. The simulated heat transfers are compared to various existing Nusselt laws of the literature. The numerical convergence and accuracy of the method is also discussed. An illustration of obtained results is given in figure 1 for FCC arrangement of spheres. Figure 1: left: temperature field around a FCC arrays of spheres for Re=50, í µí»¼ í µí± = 0.5-right: dimensionless heat transfer coefficient for different solid volume fraction at Re=100: () Massol[4], () Gun[5], present work. () [1] S. Vincent, J.C. Brändle de Motta, A. Sarthou, J.-L. Estivalezes, O. Simonin, E. Climent, A Lagrangian VOF tensorial penalty method for the DNS of resolved particle laden flows, J. Comput. Phys., 256, pp. 582-614, 2014 [2] M.-A. Chadil, S. Vincent, J.-L. Estivalezes, Accurate estimate of forces on particles using particle-resolved direct numerical simulations, submitted to Acta Mechanica, 2018 [3] A. Massol. Simulations numériques d'écoulements à travers des réseaux fixes de spheres monodisperses et bidisperses, pour des nombres de Reynolds modérés. PhD thesis, 2004 [4] Gunn, D.J., 1972. Transfer of heat or mass to particles in fixed and fluidized beds. Int. J. Heat Mass Transfer 21, 467–476.< Réduire