The steam reformer for hydrogen production from methane is studied by a numerical method. Langmuir- Hinshelwood model is incorporated for catalytic surface reactions, and the pseudo-homogeneous model is used to take into account local equilibrium phenomena between a catalyst and bulk gas. Dominant chemical reactions are Steam Reforming (SR) reaction, Water-Gas Shift (WGS) reaction, and Direct Steam Reforming (DSR) reaction. The numerical results are validated with experimental results at the same operating conditions. Using the validated code, parametric study has been numerically performed in view of the steam reformer performance. As increasing a wall temperature, the fuel conversion increases due to the high heat transfer rate. When Steam to Carbon Ratio (SCR) increases, the concentration of carbon monoxide decreases since WGS reaction becomes more active. When increasing Gas Hourly Space Velocity (GHSV), the fuel conversion decreases due to the heat transfer limitation and the low residence time. The reactor shape effects are also investigated. The length and radius of cylindrical reactors are changed at the same catalyst volume. The longer steam reformer is, the better steam reformer performs. However, system energy efficiency decreases due to the large pressure drop.