Molecular‐Weight‐Distribution Modelling of Radical Polymerization in Batch and Continuous Reactors with Transfer to Polymer Leading to Gel Formation

Full chain‐length distribution (CLD) modelling applying the Galerkin finite‐element method^1] (FEM) to polymerization reactors featuring a certain degree of gel formation is confronted with extremely long computation times. The paper describes a new method to predict CLDs for systems where gel formation may occur. The new concept is to model a part of the CLD up to a cut‐off length L, while satisfying the full set of population balances. With transfer to polymer as the mechanism responsible for gelation, this gives rise to a closure problem, which has been solved by assuming the dead CLD beyond L to be represented by a part of a Flory distribution. The method could be proved to work by performing simulations and comparing cut‐off CLDs to full CLDs for non‐gelling systems and comparing results for different L for systems with gelation. The model is demonstrated for polymerization reactors, the batch reactor and the continuous stirred‐tank reactor (CSTR), with either disproportionation or recombination termination. Reliable results are obtained for systems with moderate gel formation. Comparing these results to those from moment models including balance equations up to the fourth moment, a number of interesting differences have been found.