This paper presents a mathematical model for pulping that considers the effect of alkali diffusion and chemical reaction. An approximation for 3-dimensional diffusion in wood chips is tested against pulping data. A method for predicting pulp viscosity is also presented. The model simulates the effect of key pulping variables on pulp kappa number, viscosity, yield, residual alkalinity and pulping uniformity. Experimental data from a number of sources are successfully predicted by the model.