Methanol as a hydrogen carrier can be reformed with steam over Cu/ZnO/Al2O3 catalysts. In this paper a comprehensive pseudo-homogenous model of a multi-tubular packed-bed reformer has been developed to investigate the impact of operating conditions and geometric parameters on its performance. A kinetic Langmuir-Hinshelwood model of the methanol steam reforming process was proposed based on the work of Peppley et al. . In addition to the kinetic model, the pressure drop and the mass and heat transfer phenomena along the reactor were taken into account. This model was verified by a dynamic model in the platform of ASPEN. The diffusion effect inside catalyst particles was also estimated and accounted for by the effectiveness factor. The simulation results showed axial temperature profiles in both tube and shell side with different operating conditions. Moreover, the lower flow rate of liquid fuel and higher inlet temperature of thermal air led to a lower concentration of residual methanol, but also a higher concentration of generated CO from the reformer exit. The choices of operating conditions were limited to ensure a tolerable concentration of methanol and CO in H2-rich gas for feeding into a high temperature polymer electrolyte membrane fuel cell (HT-PEMFC) stack. With fixed catalyst load, the increase of tube number and decrease of tube diameter improved the methanol conversion, but also increased the CO concentration in reformed gas. In addition, increasing the number of baffle plates in the shell side increased the methanol conversion and the CO concentration.