航空学报 > 2017, Vol. 38 Issue (12): 121397-121397   doi: 10.7527/S1000-6893.2017.121397

再入弹头三维非对称烧蚀外形模拟

周述光1, 国义军2, 贺立新2, 刘骁2   

  1. 1. 中国空气动力研究与发展中心 空气动力学国家重点实验室, 绵阳 621000;
    2. 中国空气动力研究与发展中心 计算空气动力研究所, 绵阳 621000
  • 收稿日期:2017-05-08 修回日期:2017-06-29 出版日期:2017-12-15 发布日期:2017-06-29
  • 通讯作者: 周述光 E-mail:zhoushuguang111@163.com
  • 基金资助:
    国家重点研发计划(2016YFA0401200)

Simulation of 3D asymmetric ablation shape of reentry missile

ZHOU Shuguang1, GUO Yijun2, HE Lixin2, LIU Xiao2   

  1. 1. State Key Laboratory of Aerodynamics, China Aerodynamics Research and Development Center, Mianyang 621000, China;
    2. Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000, China
  • Received:2017-05-08 Revised:2017-06-29 Online:2017-12-15 Published:2017-06-29
  • Supported by:
    National Key Research and Development Program of China(2016YFA0401200)

摘要: 远程战略导弹再入大气层,热环境十分严酷,烧蚀量很大,外形呈现非对称性。烧蚀变形、质量引射和表面烧蚀花纹反作用于气动力、气动热及表面温度场。它们之间的相互作用十分复杂,影响落点精度。依据等效转换,将有迎角、侧滑角的绕流转换成总迎角绕流。对于不同子午面,将其变换到纵向坐标轴沿等效流动方向的气流坐标系,通过等效迎角化为零迎角绕流。采用Monte-Carlo统计方法,生成伪随机数,对弹头表面粗糙度随机抽样,表征不同子午线上转捩初始位置的随机分布。改进了零迎角烧蚀外形计算方法,拓展了其应用范围,并通过与能反映转捩点移动方向的松耦合计算策略相结合,采用数值求解烧蚀外形方程,成功模拟了小倾角再入过程中端头非对称烧蚀外形。

关键词: 再入弹头, 非对称体, 热防护, 气动热, 烧蚀, 碳/碳材料, 松耦合模拟

Abstract: Strategic missiles are heated by hypersonic gas while they are reenter into the earth's atmosphere at small attitude angles. The thermal environment is so harsh that the loss of surface material from spacecraft through evaporation or melting caused by friction with the atmosphere is a considerable amount of ablation, thereby forming a visible asymmetric shape. Shape change, matter inject to boundary layer and a number of different types of patterns by ablation will affect the aerodynamics, aerothermodynamics and temperature distribution in the matters. Thermal environment, ablation, and the shape are the cause and effect for each other, and therefore make the impact scatter to spread out. According to equivalent transform analysis, the flow field relative to the vehicle at some angles of attack and sideslip can be equal to that at a total angle of incidence. The vertical axis at a meridian plane is transformed to the flow direction, wherefore the flow flied is at the zero angle of attack. Random sampling of surface roughs is conducted with Monte-Carlo statistical measure to simulate the random distribution of the transition regions at different circumferential positions. The rough sizes are created by the randomizer. The calculation procedure has been improved with these techniques, and the loosely coupled simulation that can exhibit the movement of the transition point towards the vertex of the nosetip is also employed. The numerical calculation method is also used to sole the nonlinear partial differential equation for the ablative shape. The irregular shape of the nosetip of ballistic missiles at small incidence angle in reentry into the earth's atmosphere can be successfully simulated with the methods applied proposed.

Key words: reentry missile, asymmetric body, thermal protection, aerothermal, ablation, carbon/carbon material, loosely coupled simulation

中图分类号: