Simulating the complete pyrolysis and charring process of phenol–formaldehyde resins using reactive molecular dynamics
Source abstract
Abstract
Abstract 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
Source: crossref · Quality: baseline
Research context
OpenAlex topics: Thermal and Kinetic Analysis · Catalysis and Oxidation Reactions · Zeolite Catalysis and Synthesis
SciVX discovery tags are terminology-based. OpenAlex topics are source-provided classifications.
References (64)
Show source-provided references
1. Allen HJ (1958) Hypersonic flight and the re-entry problem: the twenty-first wright brothers lecture. J Aerosp Sci 25:217–227. https://doi.org/10.2514/8.760010.2514/8.7600 ↗
2. Adams MC (1959) Recent advances in ablation. ARS J 29:625–632. https://doi.org/10.2514/8.485210.2514/8.4852 ↗
3. Fyfe CA, McKinnon MS, Rudin A, Tchir WJ (1983) Investigation of the mechanism of the thermal decomposition of cured phenolic resins by high-resolution carbon-13 CP/MAS solid-state NMR spectroscopy. Macromolecules 16:1216–121910.1021/ma00241a033 ↗
4. Mathieu RD (1964) Mechanical spallation of charring ablators in hyperthermal environments. AIAA J 2:1621–1627. https://doi.org/10.2514/3.262910.2514/3.2629 ↗
5. Sykes Jr GF (1967) Decomposition characteristics of a char-forming phenolic polymer used for ablative composites. National aeronautics and space administration Hampton va Langley research
6. Uyanna O, Najafi H (2020) Thermal protection systems for space vehicles: a review on technology development, current challenges and future prospects. Acta Astronaut 176:341–356. https://doi.org/10.1016/j.actaastro.2020.06.04710.1016/j.actaastro.2020.06.047 ↗
7. ASTM E (2008) 285–80. Standard test method for oxyacetylene ablation testing of thermal insulation materials. Annu B ASTM Stand 15:
8. Vojvodich NS (1969) Hypervelocity heat protection- a review of laboratory experiments. J Macromol Sci Part A Chem 3:367–394. https://doi.org/10.1080/1060132690805381910.1080/10601326908053819 ↗
9. Park C, Lundell JH, Green MJ et al (1984) Ablation of carbonaceous materials in a hydrogen-helium arcjet flow. AIAA J 22:1491–1498. https://doi.org/10.2514/3.4858910.2514/3.48589 ↗
10. Milos FS, Chen Y-K, Mahzari M (2018) Arcjet tests and thermal response analysis for dual-layer woven carbon phenolic. J Spacecr Rocket 55:712–722. https://doi.org/10.2514/1.A3414210.2514/1.A34142 ↗
11. Tran HK (1994) Development of lightweight ceramic ablators and arc-jet test results
12. Tsuge S, Ohtani H (1997) Structural characterization of polymeric materials by pyrolysis—GC/MS. Polym Degrad Stab 58:109–130. https://doi.org/10.1016/S0141-3910(97)00031-110.1016/S0141-3910(97)00031-1 ↗
13. Car R, Parrinello M (1985) Unified approach for molecular dynamics and density-functional theory. Phys Rev Lett 55:2471–2474. https://doi.org/10.1103/PhysRevLett.55.247110.1103/PhysRevLett.55.2471 ↗
14. Gooneie A, Schuschnigg S, Holzer C (2017) A review of multiscale computational methods in polymeric materials. Polymers 9:16. https://doi.org/10.3390/polym901001610.3390/polym9010016 ↗
15. van Duin ACT, Dasgupta S, Lorant F, Goddard WA (2001) ReaxFF: a reactive force field for hydrocarbons. J Phys Chem A 105:9396–9409. https://doi.org/10.1021/jp004368u10.1021/jp004368u ↗
16. Castro-Marcano F, Kamat AM, Russo MF Jr et al (2012) Combustion of an Illinois No. 6 coal char simulated using an atomistic char representation and the ReaxFF reactive force field. Combust Flame 159:1272–128510.1016/j.combustflame.2011.10.022 ↗
17. Gao M, Li X, Guo L (2018) Pyrolysis simulations of Fugu coal by large-scale ReaxFF molecular dynamics. Fuel Process Technol 178:197–20510.1016/j.fuproc.2018.05.011 ↗
18. Zheng M, Li X, Guo L (2018) Investigation of N behavior during coal pyrolysis and oxidation using ReaxFF molecular dynamics. Fuel 233:867–87610.1016/j.fuel.2018.06.133 ↗
19. Bhoi S, Banerjee T, Mohanty K (2014) Molecular dynamic simulation of spontaneous combustion and pyrolysis of brown coal using ReaxFF. Fuel 136:326–33310.1016/j.fuel.2014.07.058 ↗
20. Liu X, Li X, Liu J et al (2014) Study of high density polyethylene (HDPE) pyrolysis with reactive molecular dynamics. Polym Degrad Stab 104:62–7010.1016/j.polymdegradstab.2014.03.022 ↗
21. Zhao T, Li T, Xin Z et al (2018) A ReaxFF-based molecular dynamics simulation of the pyrolysis mechanism for polycarbonate. Energy Fuels 32:2156–216210.1021/acs.energyfuels.7b03332 ↗
22. Liu Q, Liu S, Lv Y et al (2021) Atomic-scale insight into the pyrolysis of polycarbonate by ReaxFF-based reactive molecular dynamics simulation. Fuel 287:11948410.1016/j.fuel.2020.119484 ↗
23. Zheng M, Wang Z, Li X et al (2016) Initial reaction mechanisms of cellulose pyrolysis revealed by ReaxFF molecular dynamics. Fuel 177:130–14110.1016/j.fuel.2016.03.008 ↗
24. Diao Z, Zhao Y, Chen B et al (2013) ReaxFF reactive force field for molecular dynamics simulations of epoxy resin thermal decomposition with model compound. J Anal Appl Pyrol 104:618–62410.1016/j.jaap.2013.05.005 ↗
25. Zhang X, Wu Y, Chen X et al (2017) Theoretical study on decomposition mechanism of insulating epoxy resin cured by anhydride. Polymers 9:34110.3390/polym9080341 ↗
26. Li G, Hu P, Luo W et al (2021) Simulation of pyrolysis of crosslinked epoxy resin using ReaxFF molecular dynamics. Comput Theor Chem 1200:11324010.1016/j.comptc.2021.113240 ↗
27. Odegard GM, Jensen BD, Gowtham S et al (2014) Predicting mechanical response of crosslinked epoxy using ReaxFF. Chem Phys Lett 591:175–17810.1016/j.cplett.2013.11.036 ↗
28. Natali M, Kenny JM, Torre L (2016) Science and technology of polymeric ablative materials for thermal protection systems and propulsion devices: a review. Prog Materials Sci 84:192–275. https://doi.org/10.1016/j.pmatsci.2016.08.00310.1016/j.pmatsci.2016.08.003 ↗
29. Jiang D, van Duin ACT, Goddard WA, Dai S (2009) Simulating the initial stage of phenolic resin carbonization via the ReaxFF reactive force field. J Phys Chem A 113:6891–6894. https://doi.org/10.1021/jp902986u10.1021/jp902986u ↗
30. Qi T, Bauschlicher CW Jr, Lawson JW et al (2013) Comparison of ReaxFF, DFTB, and DFT for phenolic pyrolysis. 1. Molecular dynamics simulations. J Phys Chem A 117:11115–1112510.1021/jp4081096 ↗
31. Bauschlicher CW Jr, Qi T, Reed EJ et al (2013) Comparison of ReaxFF, DFTB, and DFT for phenolic pyrolysis. 2. Elementary reaction paths. J Phys Chem A 117:11126–1113510.1021/jp408113w ↗
32. Desai TG, Lawson JW, Keblinski P (2011) Modeling initial stage of phenolic pyrolysis: graphitic precursor formation and interfacial effects. Polymer 52:577–585. https://doi.org/10.1016/j.polymer.2010.11.01810.1016/j.polymer.2010.11.018 ↗
33. Ng TY, Yeo JJ, Liu ZS (2012) A molecular dynamics study of the thermal conductivity of graphene nanoribbons containing dispersed Stone–Thrower–Wales defects. Carbon 50:4887–4893. https://doi.org/10.1016/j.carbon.2012.06.01710.1016/j.carbon.2012.06.017 ↗
34. Zou J-H, Ye Z-Q, Cao B-Y (2016) Phonon thermal properties of graphene from molecular dynamics using different potentials. J Chem Phys 145:134705. https://doi.org/10.1063/1.496391810.1063/1.4963918 ↗
35. Zhang X, Chen Z, Chen H, Xu L (2021) Comparative studies of thermal conductivity for bilayer graphene with different potential functions in molecular dynamic simulations. Results Phys 22:103894. https://doi.org/10.1016/j.rinp.2021.10389410.1016/j.rinp.2021.103894 ↗
36. Diao C, Dong Y, Lin J (2017) Reactive force field simulation on thermal conductivities of carbon nanotubes and graphene. Int J Heat Mass Transf 112:903–912. https://doi.org/10.1016/j.ijheatmasstransfer.2017.05.03610.1016/j.ijheatmasstransfer.2017.05.036 ↗
37. Müller-Plathe F (1997) A simple nonequilibrium molecular dynamics method for calculating the thermal conductivity. J Chem Phys 106:6082–6085. https://doi.org/10.1063/1.47327110.1063/1.473271 ↗
38. Fourier J-B-J (1955) The analytical theory of heat. Dover Publishers, New York
39. Sprik M, Impey RW, Klein ML (1984) Second-order elastic constants for the Lennard-Jones solid. Phys Rev B 29:4368–4374. https://doi.org/10.1103/PhysRevB.29.436810.1103/PhysRevB.29.4368 ↗
40. Theodorou DN, Suter UW (1986) Atomistic modeling of mechanical properties of polymeric glasses. Macromolecules 19:139–15410.1021/ma00155a022 ↗
41. (2021) BIOVIA, Dassault Systèmes, BIOVIA Materials Studio
42. Hall SA, Howlin BJ, Hamerton I et al (2012) Solving the problem of building models of crosslinked polymers: an example focussing on validation of the properties of crosslinked epoxy resins. PLoS ONE 7:e4292810.1371/journal.pone.0042928 ↗
43. Purse M, Edmund G, Hall S et al (2019) Reactive molecular dynamics study of the thermal decomposition of phenolic resins. J Compos Sci 3:3210.3390/jcs3020032 ↗
44. Mungiguerra S, Di Martino GD, Cecere A et al (2020) Ultra-high-temperature testing of sintered ZrB2-based ceramic composites in atmospheric re-entry environment. Int J Heat Mass Transf 156:119910. https://doi.org/10.1016/j.ijheatmasstransfer.2020.11991010.1016/j.ijheatmasstransfer.2020.119910 ↗
45. Thompson AP, Aktulga HM, Berger R et al (2022) LAMMPS-a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comput Phys Commun 271:108171. https://doi.org/10.1016/j.cpc.2021.10817110.1016/j.cpc.2021.108171 ↗
46. Aktulga HM, Fogarty JC, Pandit SA, Grama AY (2012) Parallel reactive molecular dynamics: numerical methods and algorithmic techniques. Parallel Comput 38:245–259. https://doi.org/10.1016/j.parco.2011.08.00510.1016/j.parco.2011.08.005 ↗
47. Chenoweth K, van Duin ACT, Goddard WA (2008) ReaxFF reactive force field for molecular dynamics simulations of hydrocarbon oxidation. J Phys Chem A 112:1040–1053. https://doi.org/10.1021/jp709896w10.1021/jp709896w ↗
48. Wang J, Jiang H, Jiang N (2009) Study on the pyrolysis of phenol-formaldehyde (PF) resin and modified PF resin. Thermochim Acta 496:136–142. https://doi.org/10.1016/j.tca.2009.07.01210.1016/j.tca.2009.07.012 ↗
49. Zheng F, Ren Z, Xu B et al (2021) Elucidating multiple-scale reaction behaviors of phenolic resin pyrolysis via TG-FTIR and ReaxFF molecular dynamics simulations. J Anal Appl Pyrolysis 157:10522210.1016/j.jaap.2021.105222 ↗
50. Jiang H, Wang J, Wu S et al (2012) The pyrolysis mechanism of phenol formaldehyde resin. Polym Degrad Stab 97:1527–1533. https://doi.org/10.1016/j.polymdegradstab.2012.04.01610.1016/j.polymdegradstab.2012.04.016 ↗
51. Paton K (1969) An algorithm for finding a fundamental set of cycles of a graph. Commun ACM 13:514–518. https://doi.org/10.1145/363219.36323210.1145/363219.363232 ↗
52. Wirnsberger P, Frenkel D, Dellago C (2015) An enhanced version of the heat exchange algorithm with excellent energy conservation properties. J Chem Phys 143:124104. https://doi.org/10.1063/1.493159710.1063/1.4931597 ↗
53. Ikeshoji T, Hafskjold B (1994) Non-equilibrium molecular dynamics calculation of heat conduction in liquid and through liquid-gas interface. Mol Phys 81:251–261. https://doi.org/10.1080/0026897940010017110.1080/00268979400100171 ↗
54. Chung D-H (1963) Elastic moduli of single crystal and polycrystalline MgO. Philos Mag J Theor Exp Appl Phys 8:833–841. https://doi.org/10.1080/1478643630821384010.1080/14786436308213840 ↗
55. Chen Z, Chen Y, Liu H (2013) Pyrolysis of phenolic resin by TG-MS and FTIR analysis. Adv Mater Res 631–632:104–109. https://doi.org/10.4028/www.scientific.net/AMR.631-632.10410.4028/www.scientific.net/AMR.631-632.104 ↗
56. Moldoveanu SC (2019) Chapter 2-pyrolysis of hydrocarbons. In: Moldoveanu SC (ed) Moldoveanu SCBT-P of OM, 2nd edn. Elsevier, pp 35–16110.1016/B978-0-444-64000-0.00002-0 ↗
57. Costa L, di Montelera LR, Camino G et al (1997) Structure-charring relationship in phenol-formaldehyde type resins. Polym Degrad Stab 56:23–35. https://doi.org/10.1016/S0141-3910(96)00171-110.1016/S0141-3910(96)00171-1 ↗
58. Cui J, Yan Y, Liu J, Wu Q (2008) Phenolic resin-MWNT nanocomposites prepared through an in situ polymerization method. Polym J 40:1067–1073. https://doi.org/10.1295/polymj.PJ200811810.1295/polymj.PJ2008118 ↗
59. Cai J, Liu R (2007) Research on water evaporation in the process of biomass pyrolysis. Energy Fuels 21:3695–3697. https://doi.org/10.1021/ef700442n10.1021/ef700442n ↗
60. Lousada CM, Sophonrat N, Yang W (2018) Mechanisms of formation of H, HO, and water and of water desorption in the early stages of cellulose pyrolysis. J Phys Chem C 122:12168–12176. https://doi.org/10.1021/acs.jpcc.8b0217310.1021/acs.jpcc.8b02173 ↗
61. Niu B, Ren S, Liu Y et al (2011) Low-temperature oxidation of oil components in an air injection process for improved oil recovery. Energy Fuels. https://doi.org/10.1021/ef200891u10.1021/ef200891u ↗
62. Sivaramakrishnan R, Brezinsky K, Vasudevan H, Tranter RS (2006) A shock-tube study of the high-pressure thermal decomposition of benzene. Combust Sci Technol 178:285–305. https://doi.org/10.1080/0010220050029234010.1080/00102200500292340 ↗
63. Clayton WA, Kennedy PB, Evans RJ, et al (1968) Thermal properties of ablative chars. Seattle
64. Makeev MA, Srivastava D (2009) Thermal properties of char obtained by pyrolysis: a molecular dynamics simulation study. Appl Phys Lett 95:181908. https://doi.org/10.1063/1.324963210.1063/1.3249632 ↗