考虑稀薄效应的再入飞行器前缘热化学反应与烧蚀模型研究

An investigation of thermochemical reaction and aerodynamic ablation model on leading edges of reentry vehicles considering rarefaction effect

  • 摘要: 高温气流与飞行器表面材料的多物理、化学作用导致飞行器表面形貌显著变化,影响流动结构演化与飞行器气动力、热特性。准确预测再入过程的表面材料烧蚀过程对于飞行器热防护系统的设计至关重要。现有的气动烧蚀数值模拟研究主要针对固定壁面温度条件的流场变化,忽视了烧蚀过程中复杂化学反应与壁面材料性质差异导致的升温过程及烧蚀形貌的影响。本研究基于直接模拟蒙特卡罗(DSMC)方法,耦合壁面能量守恒方程,采用开源程序SPARTA对飞行器再入时的气动加热过程开展解耦分析。以柱体模型为例,综合激波后的气-气与气-固化学反应,对解耦后的壁面升温和模型烧蚀两个过程分别建立了相应的控制方程,并分析了在二维条件下的热化学反应与气动烧蚀机制。研究结果表明,本研究开发的可计算烧蚀模型不仅能够提升激波后气体分子内能采样的准确性,也能够复现文献中已有的烧蚀形貌。该模型不仅复现了圆柱烧蚀时前后缘的差异烧蚀形貌,还在引入表面粗糙度对气动加热的放大效应后,使球锥烧蚀退缩距离的预测结果与实验数据的相对误差控制在5%以内。该方法首次在DSMC框架内综合体现了壁面材料特性、变壁温效应以及表面粗糙度对烧蚀过程的物理影响,能够为深入理解变壁面温度下的复杂热化学非平衡现象提供理论依据与数据支撑。

     

    Abstract: High-temperature airflow interacts with the surface material of the aircraft through multiple physical and chemical processes, significantly altering the aircraft's surface morphology, which consequently affects the evolution of flow structures as well as the aerodynamic and thermal characteristics of the vehicle. Accurate prediction of the ablation process during re-entry is crucial for designing thermal protection systems. Existing numerical simulations of aerodynamic ablation primarily focused on flow fields under fixed wall temperature conditions, neglecting the influences of complex chemical reactions and material property differences on the heating process and ablation morphology during ablation. This study employed the Direct Simulation Monte Carlo (DSMC) method, coupled with the wall energy conservation equation, and utilized the open-source program SPARTA to conduct a decoupled analysis of the aerodynamic heating process during vehicle re-entry. Using a cylindrical model as an example, corresponding governing equations were established separately for the decoupled processes of wall heating and model ablation. By integrating gas-gas and gas-solid chemical reactions behind the shock wave, the thermochemical reaction and aerodynamic ablation mechanisms under two-dimensional conditions were analyzed. The results indicate that the developed computable ablation model not only improves the accuracy of internal energy sampling for gas molecules behind the shock wave but also successfully reproduces ablation morphologies documented in existing literature. This model not only reproduces the difference in erosion morphology between the front and rear edges of a cylinder but also, after incorporating the amplification effect of surface roughness on aerodynamic heating, controls the relative error between the predicted retreat distance of the ball cone erosion and experimental data within 5%. This method represents the first comprehensive integration of wall material properties, variable wall-temperature effects, and surface roughness into the physical modeling of the ablation process within a DSMC framework, providing a theoretical basis and data support for deepening the understanding complex thermochemical non-equilibrium phenomena under variable wall temperatures.

     

/

返回文章
返回