1. 上海交通大学 机械与动力工程学院,上海,200240
2. 空天飞行空气动力科学与技术全国重点实验室,绵阳,621000
李芹(1995—),女,山东淄博人,博士研究生,研究方向:高速化学非平衡流动及界面气固效应建模. E-mail:liqin01@sjtu.edu.cn
杨肖峰*,副研究员,研究方向:高温非平衡流动及气动热预测. E-mail: yxf.2006@tsinghua.org.cn
董威*,教授,研究方向:高速飞行器气动热预测. E-mail: wdong@sjtu.edu.cn
纸质出版:2024
移动端阅览
李芹, 杨肖峰, 董威, 等. 压力影响硅基防热材料界面多相催化的微观机理研究[J]. 空气动力学学报, 2024, 42(4): 84−95.
LI Qin, YANG Xiaofeng, DONG Wei, et al. Microscale mechanism study of pressure effects on interfacial heterogeneous catalysis of silica-based thermal protective materials[J]. 2024, 42(4): 84-95.
李芹, 杨肖峰, 董威, 等. 压力影响硅基防热材料界面多相催化的微观机理研究[J]. 空气动力学学报, 2024, 42(4): 84−95. DOI: 10.7638/kqdlxxb-2023.0195
LI Qin, YANG Xiaofeng, DONG Wei, et al. Microscale mechanism study of pressure effects on interfacial heterogeneous catalysis of silica-based thermal protective materials[J]. 2024, 42(4): 84-95. DOI: 10.7638/kqdlxxb-2023.0195.
反应分子动力学模拟是阐明高温壁面效应微观致热机理、深化高速飞行器非平衡气动热认知的重要途径之一,然而通过人为增加压力实现计算效率提升的微观模拟,往往会导致反应路径和速率系数差异,进而影响气动热,造成机理认知偏差。采用基于ReaxFF力场的分子动力学方法,以离解氧原子在硅基防热材料表面的催化复合反应体系为对象,计算分析了不同气相压力条件下的基元反应速率、表面覆盖率和复合系数,用于获得基元反应速率常数与压力的量化关系,明确通过增压提高计算效率的上限范围。结果表明,增压可导致主导反应路径从吸附相间作用至气相-吸附相作用的转变,且使基元反应速率常数-压力的关系偏离实验/飞行条件下的规律。在
1200
K条件下、单原子碰撞的压力范围内,各基元反应步骤的速率常数均随压力的降低而下降。其中,ER1~ER3复合反应的速率常数随压力呈近似线性变化,速率常数分别与压力的
1.10179
、
1.01686
和
0.91654
次方呈线性关系;LH1~LH3复合反应的速率常数与压力呈对数关系,且对数前因子显著小于非单原子碰撞区;热解附反应的速率常数与压力呈指数关系。根据气相压力影响催化反应机制的微观机理,以基元反应速率常数-压力关系可以稳定解析为判据,提出了人为增加压力的约束上限条件:以体系高度为特征长度的努森数应大于10
2
量级,以保证气固单原子碰撞。相关研究为气固界面反应的分子模拟方法和防热材料微观催化数据的累积提供了支撑。
Reactive molecular dynamics simulation provides an important approach to elucidate the microscale heating mechanisms involved in high-temperatu
re wall effects and to deepen our understanding of non-equilibrium aerothermodynamics of high-speed aircraft. However
microscale simulations that enhance the computational efficiency by artificially increasing the gaseous pressure often lead to discrepancies in reaction pathways and rate coefficients
thereby affecting the aerothermodynamics and causing the misunderstanding of reaction mechanisms. Using a molecular dynamics approach based on the ReaxFF force field
this investigation addresses the catalytic recombination reaction system of dissociated oxygen atoms on the surface of silica-based thermal protective materials. The primary objectives involve the computation and analysis of elementary reaction rates
surface coverages
and recombination coefficients under various gaseous pressure conditions. The purpose is to establish a quantitative correlation between the rate constants of elementary reactions and pressure
thereby elucidating the upper bounds of computational efficiency enhancement through pressure augmentation. The results indicate that pressure augmentation leads to a transition in the dominant reaction pathway from adsorbate-adsorbate interactions to gas-adsorbate interactions. Furthermore
it causes a deviation in the relationship between the rate constants of elementary reactions and pressure from the patterns observed under experimental/flight conditions. At
1200
K
within the pressure range associated with single-atom collisions
a consistent decrease in the rate constants of individual elementary reaction steps is observed as the pressure decreases. Among them
the rate constants of ER1~ER3 recombination exhibit a linear relationship with pressure. Specifically
the rate constants can be expressed as a power law function of pressure
with exponents of
1.10179
1.01686
and
0.91654
r
espectively. The rate constants of LH1~LH3 recombination reactions display a logarithmic dependence on pressure
with significantly smaller pre-logarithmic factors compared to those observed in the non-single collision region. The rate constant of thermal desorption is exponentially related to pressure. Based on the microscale mechanism of catalytic reaction influenced by gaseous pressure
and with the stable relationship between elementary reaction rate constants and pressure
a constraint upper limit for artificially increasing pressure is proposed. The Knudsen number with the system height as the characteristic length should be greater than the magnitude of
$ {10^2} $
<
img text_id='M2' class='formula-img' style='display:none;' src='2023-195_M2.jpg'>
to ensure the gas-solid single collision. This investigation provides support for molecular simulation methods of gas-solid interface reactions and the accumulation of microscale catalytic data for thermal protective materials.
MO F, GAO Z X, JIANG C W, et al. Progress in the numerical study on the aerodynamic and thermal characteristics of hypersonic vehicles: High-temperature chemical non-equilibrium effect[J]. Scientia Sinica (Physica, Mechanica Astronomica), 2021, 51(10): 22-37 (in Chinese).
GUI Y W, LIU L, WEI D . Combined thermal phenomena issues of long endurance hypersonic vehicles[J ] . Acta Aerodynamica Sinica, 2020 , 38 ( 4 ): 641 − 650 (in Chinese).
LIN L, WU B, WU C K . Studies on surface catalytic effect of materials in a high-temperature gas flow[J ] . Acta Aerodynamica Sinica, 2001 , 19 ( 4 ): 407 − 413 (in Chinese).
MENG S H, JIN H, WANG G L, et al . Research advances on surface catalytic properties of thermal protection materials[J ] . Acta Aeronautica et Astronautica Sinica, 2014 , 35 ( 2 ): 287 − 302 (in Chinese).
DING M S, DONG W Z, GAO T S, et al . Computational analysis of influence of differences in local catalytic properties on aero-thermal environment[J ] . Acta Aeronautica et Astronautica Sinica, 2018 , 39 ( 3 ): 44 − 54 (in Chinese).
ZHOU K, ZHU X J, OU D B, et al . Wall catalytic effects on aerodynamic heat of hypersonic flow[J ] . Acta Aerodynamica Sinica, 2020 , 38 ( 4 ): 660 − 667 (in Chinese).
DU Y X, YANG X F, XIAO G M, et al. Progress on multiscale coupling thermal effects of hypersonic vehicles[J]. Scientia Sinica (Technologica), 2023 (in Chinese).
SU S Y, SHI Y L, LIU S, et al . Finite-rate surface catalysis effects on aero-heating environment of a reentry capsule[J ] . Acta Aerodynamica Sinica, 2018 , 36 ( 5 ): 878 − 884 (in Chinese).
MOZETIČ M, ZALAR A . Recombination of neutral oxygen atoms on stainless steel surface[J ] . Applied Surface Science, 2000 , 158 ( 3-4 ): 263 − 267 .
BALAT-PICHELIN M, BÊCHE E . Atomic oxygen recombination on the ODS PM 1000 at high temperature under air plasma[J ] . Applied Surface Science, 2010 , 256 ( 16 ): 4906 − 4914 .
BALAT-PICHELIN M, BÊCHE E, SCITI D, et al . Emissivity, catalycity and microstructural characterization of ZrB 2 –SiC fiber based UHTC at high temperature in a non-equilibrium air plasma flow[J ] . Ceramics International, 2014 , 40 ( 7 ): 9731 − 9742 .
NASUTI F, BARBATO M, BRUNO C . Material-dependent catalytic recombination modeling for hypersonic flows[J ] . Journal of Thermophysics and Heat Transfer, 1996 , 10 ( 1 ): 131 − 136 .
BARBATO M, REGGIANI S, BRUNO C, et al . Model for heterogeneous catalysis on metal surfaces with applications to hypersonic flows[J ] . Journal of Thermophysics and Heat Transfer, 2000 , 14 ( 3 ): 412 − 420 .
NORMAN P, SCHWARTZENTRUBER T, COZMUTA I. A computational chemistry methodology for developing an oxygen-silica finite rate catalytic model for hypersonic flows[C]// 42nd AIAA Thermophysics Conference, Honolulu, Hawaii. Reston, Virginia: AIAA, 2011: AIAA2011-3644.
COZMUTA I. Molecular mechanisms of gas surface interactions in hypersonic flow[C]// 39th AIAA Thermophysics Conference, Miami, Florida. Reston, Virginia: AIAA, 2007: AIAA2007-4046.
ZHAO J, SUN X C, ZHANG J, et al . Research advances on heat and mass transfer coupling effect at gas-solid interface for thermal protection materials[J ] . Acta Aeronautica et Astronautica Sinica, 2022 , 43 ( 10 ): 527577 (in Chinese).
HE L C, CUI Z L, SUN X C, et al . Sensitivity analysis of the catalysis recombination mechanism on nanoscale silica surfaces[J ] . Nanomaterials, 2022 , 12 ( 14 ): 2370 .
VAN DUIN A C T, STRACHAN A, STEWMAN S, et al . ReaxFF SiO reactive force field for silicon and silicon oxide systems[J ] . The Journal of Physical Chemistry A, 2003 , 107 ( 19 ): 3803 − 3811 .
KULKARNI A D, TRUHLAR D G, GOVERAPET S S, et al. Oxygen interactions with silica surfaces: Coupled cluster and density functional investigation and the development of a new ReaxFF potential [J]. The Journal of Physical Chemistry C, 2013, 117(1): 258-69.
NORMAN P, SCHWARTZENTRUBER T. A computational chemistry methodology for developing an oxygen-silica finite rate catalytic model for hypersonic flows: Part Ⅱ[C]// 43rd AIAA Thermophysics Conference, New Orleans, Louisiana. Reston, Virginia: AIAA, 2012: AIAA2012-3097.
YANG Y. Multi-approach analysis of catalytic recombination in hypersonic flow[D]. Daejeon, Korea: Sejong University, 2020.
YANG Y, PEDDAKOTLA S A, KUMAR R, et al . Effect of argon gas in oxygen catalytic recombination on a silica surface: a reactive molecular dynamics study[J ] . Acta Astronautica, 2020 , 175 : 531 − 539 .
ZENG Q X. Spectral diagnosis and simulation of the interaction between silicon-based materials and oxygen atoms[D]. Harbin: Harbin Institute of Technology, 2020 (in Chinese).
CUI Z L, ZHAO J, YAO G C, et al . Competing effects of surface catalysis and ablation in hypersonic reentry aerothermodynamic environment[J ] . Chinese Journal of Aeronautics, 2022 , 35 ( 10 ): 56 − 66 .
HE L C, YE Z F, CAO Y F, et al . Influence of surface nanostructures on the catalytic recombination of hyperthermal non-equilibrium flow[J ] . Physics of Fluids, 2023 , 35 ( 4 ): 042108 .
YE Z F, ZHAO J, ZHU B J, et al . A hybrid CFD-RMD multiscale coupling framework for interfacial heat and mass simulation under hyperthermal ablative conditions[J ] . Internationa l Journal of Heat and Mass Transfer, 2023 , 213 : 124341 .
MASSUTI-BALLESTER B, PIDAN S, HERDRICH G, et al . Recent catalysis measurements at IRS[J ] . Advances in Space Research, 2015 , 56 ( 4 ): 742 − 765 .
CUI Z L, ZHAO J, YAO G C, et al . Molecular insight of the interface evolution of silicon carbide under hyperthermal atomic oxygen impact[J ] . Physics of Fluids, 2022 , 34 ( 5 ): 052101 .
CHEUNG T M, SCHRIJER F F J, PARK G . Nitrogen catalytic recombination on copper oxide in tertiary gas mixtures[J ] . Journal of Spacecraft and Rockets, 2016 , 53 ( 4 ): 644 − 653 .
MASSUTI-BALLESTER B, HERDRICH G . Experimental methodology to assess atomic recombination on high-temperature materials[J ] . Journal of Thermophysics and Heat Transfer, 2018 , 32 ( 2 ): 353 − 368 .
VAN DUIN A C T, DASGUPTA S, LORANT F, et al . ReaxFF: a reactive force field for hydrocarbons[J ] . Journal of Physical Chemistry A, 2001 , 105 ( 41 ): 9396 − 9409 .
PLIMPTON S . Fast parallel algorithms for short-range molecular dynamics[J ] . Journal of Computational Physics, 1995 , 117 ( 1 ): 1 − 19 .
NORMAN P, SCHWARTZENTRUBER T, COZMUTA I. Modeling air-SiO 2 surface catalysis under hypersonic conditions with ReaxFF molecular dynamics[C ] // 10th AIAA/ASME Joint Thermophysics and Heat Transfer Conference, Chicago, Illinois. Reston, Virginia: AIAA, 2010: AIAA2010-4320.
GARCIA A L, WAGNER W . Generation of the Maxwellian inflow distribution[J ] . Journal of Computational Physics, 2006 , 217 ( 2 ): 693 − 708 .
LI Q, YANG X F, DONG W, et al . Transition of catalytic recombination pathways on silica-based thermal protection materials at different temperatures using reactive molecular dynamics method[J ] . Physics of Fluids, 2023 , 35 ( 1 ): 012105 .
LIU L P, WANG G L, WANG Y G, et al . Test methods for determining surface catalytic properties of thermal protection materials in high enthalpy chemical non-equilibrium flows[J ] . Acta Aeronautica et Astronautica Sinica, 2017 , 38 ( 10 ): 96 − 104 (in Chinese).
CARLETON K L, MARINELLI W J . Spacecraft thermal energy accommodation from atomic recombination[J ] . Journal of Thermophysics and Heat Transfer, 1992 , 6 ( 4 ): 650 − 655 .
KOVALEV V L, KOLESNIKOV A F . Experimental and theoretical simulation of heterogeneous catalysis in aerothermochemistry (a review)[J ] . Fluid Dynamics, 2005 , 40 ( 5 ): 669 − 693 .
ZHANG S L, WANG Z H . Effects of chemical energy accommodation on nonequilibrium flow and heat transfer to a catalytic wall[J ] . Chinese Journal of Aeronautics, 2022 , 35 ( 10 ): 165 − 175 .
0
浏览量
0
下载量
0
CSCD
关联资源
相关文章
相关作者
相关机构
京公网安备11010802024621