卫星回收舱再入过渡流区气动热数值计算

简介

返回式卫星进入地球大气层95km高度左右时,速度达到Ma>20,回收舱被高超声速稀薄来流形成的弓形激波环绕,气动加热问题非常明显,准确预测过渡流区的气动热成为一个十分突出的问题。过渡区由于气体空气分子仍然比较密集,基于分子动力学的直接蒙特卡洛模拟方法(DSMC)极其耗费计算资源,而求解Navier-Stokes(N-S)方程的方法误差较大。采用添加二阶滑移条件的N-S方程求解过渡流区气动热,并与开源DSMC2V程序计算结果对比,研究了舱体母线热流、压强系数变化;通过对比两种方法下舱体前缘弓形激波及流场参数变化,分析了滑移条件影响壁面和流场参数的机理。结果表明:回收舱再入过渡区时,钝头驻点区为高压、高热流区,锥身区气动热和壁面压强保持在较低水平。带滑移条件的计算机流体力学方法(ComputationalFluidDynamics,CFD)计算得到的热流、压强系数与DSMC结果吻合良好,具有计算效率高、精度较高的优势。对流场压力、温度、速度等参数分析显示,滑移条件中壁面速度滑移和温度跳跃的加入,改变了壁面流动参数,进而改善了壁面热流和压强的准确模拟能力。滑移条件对外流场参数的影响极小,没有改善稀薄流区流场的模拟能力和激波捕获能力,模拟得到的温度、速度等流场参数及激波位置、激波层厚度等与DSMC结果仍有差别。可以认为,文章采用方法能够满足工程上快速高效预测卫星回收舱再入过渡流区气动热需要。