In this work a detailed mathematical model for free radical suspension polymerization of methyl methacrylate (MMA) in water is developed. This model is based on sound principles such as the free volume theory to account for the diffusion limited reactions in suspension polymerization. Additionally, the complex polymerization kinetics process of the aqueous suspension polymerization of MMA is studied as a one-dimensional numerical experiment. For this purpose, the polymerization process is modeled as a moving boundary mass transfer problem coupled with polymerization reactions. The Galerkin finite element method is used to simultaneously solve the nonlinear governing equations. The model predictions for conversion and average molecular weights vs. time were found to be in close agreement with laboratory data. It is believed that this work, as it provides fundamental understanding of the process, it might contribute to a more rational design of polymerization reactors. POLYM. ENG. SCI., 51:670-678, 2011. (C) 2010 Society of Plastics Engineers