Simulating the complete pyrolysis and charring process of phenol–formaldehyde resins using reactive molecular dynamics

We examine the mechanism of pyrolysis and charring of large (> 10,000 atom) phenol–formaldehyde resin structures produced using pseudo-reaction curing techniques with formaldehyde/phenol ratios of 1.0, 1.5 and 2.0. We utilise Reactive Molecular Dynamics (RMD) with a hydrocarbon oxidation parameter set to simulate the high-temperature thermal decomposition of these resins at 1500, 2500 and 3500 K. Our results demonstrate that the periodic removal of volatile pyrolysis gasses from the simulation box allows us to achieve near complete carbonisation after only 2 ns of simulation time. The RMD simulations show that ring openings play a significantly larger role in thermal decomposition than has previously been reported. We also identify the major phases of phenolic pyrolysis and elucidate some of the possible mechanisms of fragment formation and graphitisation from the RMD trajectories and compute the thermal and mechanical properties of the final pyrolysed structures. Graphical abstract