A mathematical model for the pyrolysis reaction of polystyrene (PS) in a semi-batch reactor has been presented. The thermal degradation of PS was flexibly modeled by a combination of random and specific chain-end scissions. Numerical simulation was used to investigate the effect of operating conditions on the PS products spectrum, the results of which were validated by the experimental data. It was found that as the reaction temperature increased (decreased), the monomer fraction in the products became lower (higher) while the trimer higher (lower). No significant variation in the product composition was, however, observed while constant temperature was maintained. These results indicate the reaction temperature is an effective manipulated variable for the control of products composition of PS pyrolysis. The calculation of the optimum temperature trajectories through the optimization study can thus be of interest for achieving productivity enhancement in plastics pyrolysis processes.