Numerical study of forced convection heat transfer from arrays of prolate particles is performed using the second-order Immersed Boundary-Lattice Boltzmann Method (IB-LBM). Prolate particle is studied with aspect ratio of 2.5 with solid volume fraction variation from 0.1 to 0.3. For each solid volume fraction, arrays of prolate particles are generated and simulations have been performed to calculate Nusselt number for four different Hermans orientation factors and various Reynolds numbers. From the simulation results, it has been observed that, for any specific value of Hermans orientation factor, Nusselt number increases with the increase of the Reynolds number and solid volume fraction. More importantly, it is found that the effect of orientations on Nusselt number is significant. Nusselt number correlation is developed for ellipsoidal particles as function of Reynolds number, Prandtl number, solid volume fraction, and orientation factors. This correlation is valid for and .