高超声速流动壁面催化复合气动加热特性

针对高超声速流动壁面催化特性,计算了不同壁面催化复合系数条件下的球锥驻点热环境。引入了经验证的数值求解Navier Stokes方程的方法,在不同壁温500K～2500K的条件下分别分析了O 2和N 2气体在壁面处的催化复合气动加热特性,得到如下结论：（1） 原子复合放热将提高近壁面温度梯度,改变近壁面组分分布;原子复合放热一部分加热飞行器形成组分扩散热流,一部分加热近壁气体提高近壁温度梯度。（2） 在壁面催化复合系数较小时,原子复合放热主要转化为组分扩散加热,对于不同壁面温度,壁面催化复合系数α<0.1时,单一气体反应组分扩散热流小于总热流的20%。     

再入返回器极端热载荷预测方法

为快速预测返回过程再入器的极端气动热载荷情况,文章以充气式再入器为研究对象,基于动力学运动方程及Kemp-riddell气动热工程公式,采用龙格-库塔方法开展了136组工况的返回过程数值计算,获得了充气式再入系统返回过程的轨迹弹道与驻点热流密度变化情况,研究了驻点热流密度峰值和峰值出现高度与弹道系数、球头半径及再入角度的关系,发现驻点热流密度随弹道系数、再入角度的增加而增加、与球头半径的二次方成反比;但极端热载荷出现高度随弹道系数增加而降低,与球头半径和再入角度无关。文章提出了航天器以第一宇宙速度返回再入时极端热载荷的工程经验公式,采用公式对飞船返回舱、返回式卫星的极端热载荷进行预测,所得结果和试验数据基本一致,表明该预测公式具有较高的准确性和较好的通用性。文章的预测方法适用于再入返回器的设计初期阶段,可快速预测返回器再入过程的极端气动热载荷,满足气动热估算需求,为再入器气动热防护方案的选择提供支持与参考。     

高超声速飞行器主动气膜冷却热防护数值仿真研究

为研究高超声速飞行器头部全覆盖保护情况下的气动热分布特征,文章提出一种微孔射流的主动气膜热防护方案,并对射流微孔分布进行优化;通过数值求解N-S方程,得到高超声速飞行器头部的驻点压力及表面、附近温度分布。研究结果表明,在主动气膜冷却热防护下,高超声速飞行器壁面温度可以降到1000 K以下。该方案可为未来高超声速飞行器的外壳设计提供参考。     

可降低气动热效应的类凹腔外形优化设计

针对飞行器上由凸起物形成的类凹腔气动加热问题,采用数值方法求解三维N-S方程,研究了类凹腔外形结构的高超声速气动加热规律,获得了三维高超声速流场和局部热流分布,并详细分析了局部的流场结构和气动加热机理。针对凹腔前壁面热流密度过高的问题,提出并验证了一种降低前壁面边缘热流密度的优化外形,将前壁面的热流密度降低至优化前的20%左右。     

火星探测器表面材料催化特性数值模拟研究

针对火星探测器高超声速进入过程中的表面材料催化作用及其对气动热环境影响这一问题,建立了变壁面温度的火星大气表面材料催化作用模型,并基于火星大气物理化学模型和求解三维热化学非平衡N-S方程的数值方法,对典型火星探测器防热大底进行了数值模拟,获得了不同催化特性下的高超声速非平衡流场和气动热数据,分析了表面材料催化特性对气动热环境影响的规律性。研究结果表明：表面催化特性对壁面附近组分分布影响很大,催化反应进程主要受O原子浓度限制;有限催化热流随催化效率增大而增大,完全催化峰值热流比催化效率为1的有限催化峰值热流高25%~64%;表面温度随催化特性的变化规律与热流变化规律类似。有限催化模型能根据表面材料的催化特性精细化预测表面热流和温度,为防热设计提供更精确合理的参考标准。     

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

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

气动加热与热响应耦合分析及试验研究

文章介绍了Ma（马赫数）＜6的飞行器典型部位的气动加热与热响应耦合分析及试验验证技术。气动加热计算模块采用的是工程算法,同时将其集成到自主开发的结构有限元三维温度场计算ASTSA软件平台上,实现了气动加热与热响应耦合分析功能。在结构地面热试验方面,成功研发了全方程热流密度PLC控制热模拟试验系统,并利用该系统进行了高速飞行器结构气动加热与热响应耦合地面模拟试验。试验与数值分析的结果对比显示,二者吻合较好,验证了本数值分析方法的正确性。     

高超声速气流非结构化网格DSMC算法研究

高超声速气动热力学环境的研究是直接涉及飞行器轨道控制、热防护设计的关键问题之一。文章通过研究稀薄气体热化学非平衡态中的热力学环境,采用非结构化DSMC程序对“火星探路者号”（Mars Pathfinder）探测器的Ballute减速装置在地球大气层和火星大气层中的高超声速飞行进行了数值模拟,计算得到了流场的温度分布、探测器壁面的热流密度分布,分析表明稀薄气体热化学非平衡态对飞行器流场有影响。将仿真结果与NASA兰利研究中心的计算结果作了比较,二者吻合很好。研究结果可用于飞行器热防护设计。     

传感器催化特性差异对气动热影响的计算分析

针对高温气体效应及壁面催化效应对气动热环境影响,考虑高温空气各种化学反应、分子振动激发、流动中的非平衡效应以及壁面催化复合反应,通过数值求解热化学非平衡N-S方程和壁面能量、质量平衡方程,完善了有限催化条件下高温热化学非平衡流场气动热环境数值模拟方法和程序。在此基础上,针对廉金属热电偶传感器热流测量问题,开展了不同条件下高超声速热化学非平衡流场气动热环境数值模拟,分析了催化特性差异对局部气动热环境(传感器表面热流)的影响规律,为试验数据的误差带分析、修正处理和使用提供参考。研究表明:1)催化特性差异会给局部区域带来很强质量扩散热流,使总热流发生跳变,给传感器热流测量带来不可忽视的误差;2)材料催化特性差异越大,热流跳变量越大,某些条件下时,局部热流值将远远高于全表面FCW模拟的结果,其影响量可达100%以上;3)本文计算条件下,飞行马赫数越大、飞行高度越低,催化特性差异的影响越大;4)催化特性差异带来的影响还与飞行攻角、飞行器表面温度等因素存在一定关联,在催化复合系数相同情况下,表面温度越高,影响量越大。     

光子晶体热控结构防护性能研究

光子晶体在红外波段同时具有高反射率和高辐射率,可作为高超声速飞行器的防热材料。文章利用FLUENT软件和辐射模型中的离散坐标算法对有/无光子晶体防热层结构进行了气动加热数值模拟计算,得到了在不同热流密度条件下基体的温度分布;通过计算结果的对比发现,采用光子晶体防热层结构使基体温度下降约360～515 K,具有良好的热控性能。光子晶体是一种具有广阔应用前景的热防护材料,作为高超声速飞行器的防热层结构的可行性与有效性还有待深入研究。     

基于代理模型技术的返回舱外形优化

针对返回舱再入过程中需要满足的复杂气动力热务件,在典型高超声速条件下以升阻比和阻力系数为目标,以驻点热流和容积率为约束条件进行返回舱外形多目标优化设计。气动特性分析采用基于三维Navier-Stokes方程和结构网格的计算流体力学数值仿真（CFD）方法,采用Fay-Riddell经验公式计算驻点热流,采用NSGA－II优化算法进行全局寻优,为了增强多目标优化设计的计算效率,利用Kriging代理模型替代CFD计算,并引入改进的EI函数加点策略,大大减少了构建代理模型时所需的样本点数目。优化计算结果表明,代理模型计算结果与CFD计算结果误差可以控制在7％以内,基于代理模型技术的优化过程可以节省95％以上的计算开销。因此,该方法能够较好地满足初步设计的要求,为返回舱外形优化提供有益的借峪。     

喷管尺寸对电弧风洞流场和热环境的影响分析

为获得电弧风洞喷管尺寸对试验流场以及模型表面热流的影响规律,针对某特定模拟参数试验状态,采用高焓流动数值模拟方法对不同尺寸锥形喷管下的球柱校核模型试验流场进行了模拟和比较分析。研究发现,在模拟气流焓值和模型驻点热流的条件下,采用出口尺寸小的喷管所需电弧加热功率更低,同时单位流向截面上气流能量转化为模型驻点气动热的比例更低。不同喷管出口尺寸下,试验流场喷管出口区域热力学非平衡程度、波后氧原子质量分数、模型驻点区域压力以及表面传导热流和扩散热流占比都比较接近,但相较飞行状态存在明显差异;不同喷管出口尺寸下来流速度、激波脱体距离以及驻点线上平动温度之间的差异明显,喷管出口尺寸越大,其与飞行状态越接近。     

The optimization of earth orbit reentry vehicle configurations

The optimization design of aerodynamic configurations for earth orbit reentry vehicles is treated. Choosing the total heat transfer as the objective function and taking account of the design requirement of the longitudinal aerodynamical stability for reentry vehicles with the given reentry weight and the fixed fineness ratio, we utilize the interior point penalty function method to determine the minimum aerodynamic heat transfer configuration which meets the specified requirement of the longitudinal stability. The results indicate that the near double cone with a flat nose is an attractive aerodynamic shape for earth orbit reentry vehicles. The total reentry heat transfer and the weight of reentry thermal shielding can be reduced to a great extent. The longitudinal aerodynamic stability could be raised effectively.     

高超声速化学平衡流辐射流场数值解

本文用数值方法研究了高超声速大钝头体绕流辐射流场。采用粘性激波层控制方程和平衡化学反应气体模型。辐射部份包括连续辐射和线辐射。采用切平板辐射输运模型解决了空间推进的问题，计算结果表明：辐射传热与飞行器再入飞行速度，高度及头部半径密切相关，在计算的条件下，当再入速度小于７．６公里／秒时，辐射传热的影响可以忽略，当再入速度为１１公里／秒左右时，辐射传热可以与对流传热相比拟。当再入速度更大时，辐射传热远     

航天器再入全过程轴对称烧蚀热防护数值仿真研究

对航天器再入全过程轴对称烧蚀热防护进行了全过程数值仿真研究.采用修正Lees驻点热流密度方法和参考焓方法计算再入热流密度.采用JANAF模型计算烧蚀率.利用有限元法计算钝锥体再入航天器烧蚀层在移动边界条件下的轴对称温度场.采用碳化层-热解面-原始材料的轴对称碳化烧蚀模型;推导了热解气体流量计算方法.针对再入飞行大热流密度条件下,用有限元方法求解瞬态温度场时会产生的时间和空间上解的振荡问题.通过分析温度振荡现象产生的原因,采用集中热容矩阵向后差分方法解决振荡问题.计算结果表明,在时间步长选择合适的情况下,求解集中热容矩阵能够很好地解决数值振荡问题,同时烧蚀率和温度场计算比较准确.     

具有复杂围护结构光学窗口的传热特性分析

针对具有复杂围护结构光学窗口传热特性,建立了高热流作用下的光学窗口及其复杂围护结构的辐射-导热传热模型,分别讨论了窗口的辐射半透明、光谱选择性、外部加热热流以及窗口尺寸对其瞬态温度响应的影响。结果表明,忽略窗口半透明特征会产生最大31.4%误差,忽略光谱选择性会产生最大40.4%误差;窗口外部气动加热热流影响显著,温度响应随热流的增大而急剧增大;窗口厚度对其传热特性的影响较大,随着厚度的增大窗口温度响应减小;而窗口半径产生的影响可忽略不计。     

基于气动辅助变轨的变气动外形飞行器概念研究--最优控制律问题

从一个天地往返飞行器的上升轨道和再入返回轨道的优化，以及适用不同飞行任务的变轨要求的气动外形问题，提出一项基于气动力辅助变轨的变气动外形飞行器的新概念研究。对于一个固定气动外形飞行器要同时满足上升轨道有效载荷最大和再入轨道热流峰值、过载峰值及机动性能约束下的成本最低往往是困难的。若同时满足不同飞行任务：飞往太空站的运输任务，空间拦截和交会机动巡航任务及星际探测任务，则更为困难，实际上是不可能的。文章研究基于气动力辅助变轨，在热流约束下，气动外形参数变化对最优控制的影响。其结论为：热流约束下的最优控制解，包括考虑推力协同变轨，除了在非约束弧的滚转角不直接受气动外形影响外，其余的控制律，升力系数和滚转角都是气动外形参数和攻角的函数。因而变气动外形可作为一项新技术，即通过气动外形参数变化和相应的变轨策略而获得性能和成本都最佳的用途很广的一种新型飞行器。     

高超声速飞行器鼻锥迎风凹腔结构防热效能研究

对高超声速飞行器鼻锥使用迎风凹腔结构作为热防护系统时,凹腔结构的防热效能进行了数值研究。通过与相关实验对比,验证了本文数值方法的可靠性,获得了鼻锥的流场参数,外表面、凹腔内壁面的热流分布,分析了不同的凹腔尺寸参数选择对鼻锥冷却效果的影响。结果表明迎风凹腔结构能够有效的对高超声速飞行器的鼻锥尤其是驻点区域进行冷却,凹腔越深,其冷却效果越好。鼻锥气动加热的最大热流并不在尖锐唇缘的顶点,而是位于凹腔内的侧壁面上,凹腔的深度(L)变化对最大热流的出现位置影响很小。除非凹腔很浅(L/D<0.5),凹腔底面的热流值都非常小,基本可以忽略。     

进入火星大气的高温真实气体效应与气动加热研究

针对火星和地球大气分子热力学和化学行为的差异性,采用理论分析和数值模拟两种手段,研究探测器进入过程高超声速流动的分子振动激发、离解反应及热力学和化学非平衡等真实气体效应,获得不同气体模型条件下的高超声速气动加热规律,探究引起地火差异的根本原因。分析认为,探测器进入火星大气层的稀薄气体效应明显;激波层内发生CO 2气体为主的大规模离解,在极高温环境下O 2和CO也将离解;沿进入轨道的高超声速流动基本处于化学非平衡但热力学平衡状态;激波层内能量储存和分配模式因分子振动激发和化学反应而改变,分子振动激发会增强气动加热量,但均介于化学反应模型的完全非催化和完全催化壁结果之间;相同来流条件下CO 2介质高超声速气动加热强于空气介质,但真实的火星进入热载荷因大气稀薄而弱于地球再入环境。相关研究为我国未来火星探测器热防护系统设计提供技术支持。     

化学非平衡效应对返回舱再入气动力特性的影响

高空高马赫数条件下,化学非平衡效应将对飞行器气动特性产生影响,影响飞行器气动布局优化和飞行弹道设计。文章通过三维化学非平衡流动求解程序,针对再入返回器开展数值研究与机理分析,通过对比完全气体模型和化学非平衡气体模型获得的气动力参数,揭示化学非平衡效应对流场结构和气动力特性的影响和规律。结果表明,对Apollo的气动力计算结果验证了模型和计算方法;化学非平衡效应影响下,激波层内化学反应消耗大量能量,致使激波脱体距离减小,气体压缩性增强;典型状态高度为70 km,Ma=30条件下,化学非平衡效应导致返回器升力系数增大约6%、阻力系数增大约1.3%～3.3%、升阻比增大3%左右、俯仰力矩系数增大,从而使配平攻角减小约2.5＆#176;;通过机理分析,发现化学非平衡效应影响下表面压力系数发生变化的原因是飞行器周围激波形状及驻点压力改变,表现为气体沿流线经激波层、压缩区和膨胀区的历程变化;对于钝体形状的返回器,迎风面前体压力系数增加和后体压力系数降低,造成轴向力和法向力系数增大。