Se han realizado simulaciones resueltas por partículas (PR-DNS) utilizando un método de penalización viscosa implícito de segundo orden (VPM) para estudiar la transferencia de calor entre un conjunto de partículas y un fluido portador incompresible. Se utiliza una extrapolación de Lagrange acoplada a una interpolación de Taylor de alto orden para estimar con precisión los coeficientes de transferencia de calor en una esfera aislada, un arreglo cúbico centrado en la cara fijo de esferas y un empaque aleatorio de esferas. Los coeficientes de transferencia de calor simulados se comparan con éxito con varias leyes de Nusselt existentes en la literatura.