A drag force model for spheroids, referred as the spheroid model, was implemented in OpenFOAM, in order to better predict the thermochemical conversion of pulverized biomass. Our previous work has found that the spheroid model predicts more dispersed results in terms of particle velocities and local concentrations comparing to other conventional particle models under non-reactive conditions. This work takes the spheroid model one step further, by validating against experiments performed under reactive conditions with a newly implemented heat transfer model for spheroids as well as updated devolatilization kinetic parameters. In addition, simulations were conducted in a configuration similar to a pilot-scale entrained flow gasifier for more realistic scenarios. Particle mass and axial velocity development were compared accordingly using four different modelling approaches with increasing complexity. When compared with models of spheroidal shape assumptions, the sphere and simplified non-sphere model predict 61% and 43% longer residence times, respectively. The combination of the spheroid shape assumption with the heat transfer model for spheroids tends to promote drying and devolatilization. On the other hand, the traditional spherical approach leads to longer particle residence times. These opposing effects are believed to be a major contributing factor to the fact that no significant differences among modelling approaches were found in terms of syngas production at the outlet. Furthermore, particle orientation information was reported in both experiments and simulations under reactive conditions. Its dependency on gas velocity gradient under reactive conditions is similar to what was reported under non-reactive conditions. © 2020 The Authors