In the present work, the free radical polymerization of styrene is modeled by considering the phenomenology of the process (a simplified model, which does not include the diffusional effects, gel, and glass effects) in combination with an empirical model represented by an artificial neural network. Differential evolution (DE) algorithm, belonging to the class of evolutionary algorithms, is applied for developing the neural models in optimal forms. For improving the results—predicted conversion and molecular weights as function of time, temperature, and initiator concentration—different combinations between phenomenological model and neural network are tested; also, individual and stacked neural networks have been developed for the polymerization process. This methodology based on hybrid models, including neural networks aggregated in stacks, provides accurate results.