High-Temperature Oxidation of SiC-Based Composite: Rate Constant Calculation from ReaxFF MD Simulations, Part II
详细信息 查看全文
- 作者：David A. Newsome ; Debasis Sengupta ; Adri C. T. van Duin
- 刊名：The Journal of Physical Chemistry C
- 出版年：2013
- 出版时间：March 14, 2013
- 年：2013
- 卷：117
- 期：10
- 页码：5014-5027
- 全文大小：662K
- 年卷期：v.117,no.10(March 14, 2013)
- ISSN：1932-7455
文摘
Space vehicles often encounter very high temperature and harsh oxidative environments. To ensure proper thermal protection, layers composed of SiC and EPDM polymer are placed on the outer surface of the space vehicle. The O2 and H2O molecules are able to oxidize the SiC network, creating SiO2-type structures that may form a protective layer, while also pyrolyzing and burning the EPDM polymer, causing ablation. Reactive molecular dynamics simulations nicely complement experiment, as they provide direct observation and information to calculate physical parameters such as transport diffusivities and reaction constants. In this study, rate models were developed and molecular dynamics simulated trajectories were used to extract Arrhenius parameters that describe the initial stages of transport and kinetics of SiC oxidation by O2 and H2O and the combustion and pyrolysis of EPDM. The simulations showed that O2 was able to oxidize SiC more efficiently than H2O, with the transport activation barrier of O2 in the range of 40鈥?0 kJ/mol, and for H2O at 125鈥?50 kJ/mol. The oxidizer molecules created a continuous surface of SiO2 that grew on top of the SiC layer. The O atoms insert themselves in between the Si鈥揅 bonds, causing the C atoms to migrate into a carbonaceous phase in the center of the simulation box. A model of the EPDM polymer was used for the combustion and pyrolysis simulations. An Arrhenius analysis gave activation barriers of 183 kJ/mol for combustion and 213 kJ/mol for pyrolysis. By comparison, experimental observations of EPDM and similar polymers conclude a range of 100鈥?50 kJ/mol for both combustion and pyrolysis, indicating that ReaxFF can give reasonable quantitative predictions for these thermal studies, demonstrating the possibility of direct calculation of Arrhenius parameters from reactive molecular dynamics simulations.