Neste trabalho estudou-se a modelagem matemática de processo de polimerização em suspensão de estireno, envolvendo a equação de balanço populacional para a previsão da distribuição de tamanhos de partículas. Neste problema os termos referentes aos processos de quebra e coalescência de partículas são importantes. Foram estudados métodos de solução da equação de balanço populacional e a solução desta equação, para o problema de polimerização em suspensão, foi obtida pelo método discretização do domínio em classes, usando a técnica dos pivôs fixos. Foram investigados os efeitos no número de elementos de discretização (número de classes) e o modo de modelar a distribuição inicial de tamanhos de gotas de monômero, quando esta não é medida experimentalmente. Neste caso, adotou-se uma distribuição normal com média elevada e variância pequena, arbitrárias, e simulou-se o processo de agitação com ocorrência de quebra e coalescência, até que se atingisse uma distribuição que não variasse com o tempo. Esta distribuição era usada então como condição inicial para a simulação do processo. Os resultados da simulação foram validados por comparação com dados experimentais referentes à evolução da distribuição de tamanhos das partículas, obtidos da literatura. O primeiro conjunto de dados foi referente ao processo de dispersão de gotas de estireno em água (sem ocorrência de polimerização). O segundo conjunto de dados foi referente ao processo de polimerização em suspensão de estireno em reator de batelada de escala de laboratório, nas quais os efeitos da velocidade de agitação, quantidade de agente estabilizante e fração de fase dispersa foram estudados. Em ambos os casos, o modelo representou corretamente o comportamento dos dados experimentais.

The mathematical modeling of suspension polymerization of styrene accounting for the population balance equation to predict the dynamic behavior of particle size distribution was studied. In this process, the population balance is determined by the rates of breakage and coalescence. Different methods for solving population balance (BP) equations were studied and, for the case of suspension polymerization, the numerical solution was obtained using the discretized BP by the method of classes and the fixed pivot technique. The effect of the number of discretization elements (number of classes) was investigated. In order to model the initial particle size distribution (PSD), which was not measured experimentally, an initial normal distribution with high average and low variance was assumed and the agitation process with breakage and coalescence was simulated until a constant PSD is reached. This PSD is then used as the initial condition for the process simulation. Simulation results were validated by comparison with experimental data taken from the literature. The first data set involves the dispersion process of styrene droplets in water, without polymerization. The second data set refers to the process of suspension polymerization of styrene carried out in a laboratory-scale batch reactor, in which the effects of agitation speed, amount of stabilizer (surface-active agent) and disperse phase holdup were studied. In both cases, the model was able to represent correctly the trends of the experimental data.