The present research is focused on the two-step ZnO/Zn thermochemical water splitting cycle for hydrogen production. In the present paper, the numerical modeling of the first step, which involves endothermic reduction of zinc oxide (ZnO), is carried out in a cylindrical reactor using Computational Fluid Dynamics (CFD). The parametric study shows that the fractional conversion of ZnO increases with an increase in the flow rate of ZnO, while it decreases with an increase in the ZnO particle diameter and carrier gas mass flow rate. Six different reactor configurations are also assessed comprehensively. It is observed that a cylindrical reactor with a tangential inlet at the top plane and a tangential outlet at the bottom plane has higher robustness to the variation of various operating parameters with consistently high ZnO fractional conversion. Copyright (c) 2013 John Wiley & Sons, Ltd.