高温高速稀薄流的DSMC算法与流场传热分析

在非结构网格下,研究了热力学碰撞传能的6种方式和化学反应的8种类型,详细给出了上述碰撞传能和化学反应碰撞类型的子程序计算框图.在可变硬球(VHS)分子模型、Borgnakke-Larsen唯象模型、Bird的化学反应几率模型以及壁面CLL(Cercignani-Lampis-Lord)反射模型的基础上,用Fortran语言编制了能够模拟内能松弛、热力学非平衡和化学非平衡的稀薄气体DSMC(direct si mulation Monte-Carlo)源程序,在地球大气层和火星大气层中完成了Ballute减速装置的8个工况计算(其中Knudsen数从0.05变到30.0,飞行Mach数从26.3变到11.2),并与NASA Langley研究中心2007年发表的计算结果作了比较,本文的结果令人满意.另外,这里用|T-Tv|/T去分析热力学非平衡,用Damk hler数去分析化学反应非平衡,用Stan-ton数去分析飞行器壁面的传热效果,所有这些分析对高温部件的热防护设计十分有益.

激波后高温高速流场中的传热特性研究

激波后高温、高速流场中的热力学与传热特性分析是直接涉及到飞行器热防护设计与传热分析的关键问题之一,借助于多组分、考虑非平衡态气体的振动以及激波与热化学非平衡态效应的守恒积分型Navier-Stokes方程组,并用高分辨率总变差减小(TVD)格式进行求解,计算与研究了Apollo工程AS-202返回舱再入地球大气层的6个飞行工况(飞行马赫数15.52～22.63)以及Huygens飞行器再入土卫六大气层的6个工况(飞行马赫数17.29～24.47),分析了不同工况下弓形脱体激波后高温高速流场的热力学与传热特性,计算得到了沿壁面的热流密度分布、温度分布以及Stanton数分布,并与国外相关飞行数据进行了比较,两者吻合较好.相关计算可以指导有关飞行器的热防护设计.

四种飞行器绕流的三维DSMC计算与传热分析

在地球大气层与火星大气层中,使用自己编制的DSMC(direct si mulation Monte Carlo)源程序完成了四种飞行器(即Apollo,Orion,Mars Pathfinder以及Mars Microprobe)高超声速穿越稀薄气体时的三维绕流计算,给出了上述飞行器42个典型飞行工况(其中包括在地球大气层中,飞行高度从250 km变到90km,飞行攻角从45°变到-45°,Knudsen数从111.0变到0.0057,飞行速度从7.6km/s变到9.6km/s;在火星大气层中,飞行高度从141.8km变到80.28km,飞行攻角从45°变到0°,Knudsen数从100.0变到0.0546,飞行速度从7.47km/s变到6.908km/s)时详细气动力与气动热的数值结果,并与国外的飞行数据以及美国NASA(National Aeronautics and Space Administration)Langley研究中心发表的计算结果进行了比较,所得结果令人满意.文中采用了三种无量纲参数分别刻画这四种典型飞行器绕流流动的热力学非平衡、化学反应非平衡以及壁面热流分布的特征,这些结果对于指导空间飞行器的热防护气动设计十分有益.