Bulk styrene homopolymerization in cylindrical molds was studied in order to simulate the reactive injection process when the reaction kinetics was much slower than the filling process. Mathematical models were built and implemented, allowing the analysis of certain kinetic phenomena, such as the gel-effect, and the relevance of natural convection. Based on a detailed kinetic mechanism, the polymer properties (average molecular weight and polydispersity) and monomer conversion were calculated as functions of spatial coordinates and time. The results obtained show that polymer properties and monomer conversion may be very heterogeneous inside the cylindrical mold, although insensitive to variations of the diffusion coefficient, so that defects are expected to appear as a consequence of the heterogeneity of polymer properties inside the mold. Experiments confirmed the simulation results. Simulations indicate that heterogeneities may be expected to increase when the reaction rates are one order of magnitude higher, as in acrylic acid polymerizations.