The direct reduction of iron ore pellets with syngas or hydrogen is a promising technology to reduce the CO2 emissions of the iron and steel industry. The conversion rate of single iron ore pellets to iron is extensively investigated. In most of these studies, a shrinking core model is employed to reproduce the experimental observations. However, this model presents an inherent bias by assuming a sharp separation between a fully converted region and a fully unreacted one. Herein the present study, an improved porous solid model is proposed. This model solves the mass balances of the individual gas species and the solid ones assuming spherical symmetry. The governing equations, the main algorithm, and validation cases are presented. The present model also offers wide flexibility to incorporate complex phenomena such as porosity changes or carbon deposition. Furthermore, the proposed model is integrated into a computational fluid dynamics (CFD) environment. It is verified that identical input parameters yield almost identical results in both frameworks, opening the gate toward reliable CFD simulations of industrial‐scale reactors.