Abstract
As societal interest in recycling of plastics increases, modeling thermochemical recycling of vinyl polymers, e.g., via pyrolysis or reactive extrusion, becomes increasingly important. A key aspect remains the reliability of the simulation results with fewer evaluation studies regarding convergence as in the polymerization or polymer reaction engineering field. Using the coupled matrix-based Monte Carlo (CMMC) framework, tracking the unzipping of individual chains according to a general intrinsic reaction scheme consisting of fission, β-scission, and termination, it is however illustrated that similar convergence demands as in polymerization benchmark studies can be employed, i.e., threshold values for the average relative error predictions on conversion and chain length averages can be maintained. For this illustration, three theoretical feedstocks are considered as generated from CMMC polymer synthesis simulations, allowing to study the effect of the initial chain length range and the number of defects on the convergence demands. It is shown that feedstocks with a broader chain length distribution and a long tail require a larger Monte Carlo simulation volume, and that the head−head effects play a key role in the type of degradation mechanism and overall degradation rate. A minimal number of chains around 5 × 105 is needed to properly reflect the degradation kinetics. A certain degree of noise can be allowed at the higher carbon-based conversions due to the inevitable decrease in number of chains.