矢量喷管内燃气辐射与壁面温度的耦合计算

选取内、外调节片和隔热屏建立几何模型,基于封闭腔净辐射模型和壁面热平衡模型建立了燃气辐射与喷管壁面温度的耦合算法.波段为1~5μm的气体辐射采用窄波带模型计算,其他波段不考虑气体辐射,建立辐射净热流密度-有效辐射亮度-壁面温度的关联式求解燃气与壁面的辐射换热,采用牛顿-拉斐尔森迭代法求解壁面热平衡方程计算其温度.对某轴对称矢量喷管(偏转20°),计算了喷管壁面的红外光谱辐射和辐射净热流密度,以及各部分结构的温度.作为验证,还计算了文献中某液体火箭发动机轴对称矢量喷管壁面的辐射净热流密度,与文献的结果进行对比一致性较好.研究表明:轴对称矢量喷管偏转段沿周向的辐射热流密度和温度差异很大,沿偏转方向部位壁面的温度和辐射热流密度都较低,偏转方向壁面的温度比相反方向大约低10%,辐射热流密度大约低50%.      

考虑可压缩与热传导的壁面函数边界条件及其应用

考虑到可压缩和热传导效应的壁面函数边界条件,被耦合到了采用k-ω两方程湍流模型、用有限体积法求解N-S方程的程序中。壁面函数基于耦合的速度和温度型,并且在边界层内的粘性子区和对数区内一致有效。引入壁面函数边界条件后,通过算例验证在y <100的范围内,得到的物面压力、摩阻、热流与实验结果比较,结果可靠。而无壁面函数边界条件时,要得到相同精度的结果,要求y ≈1。壁面函数的引入,为工程上准确预测飞行器在湍流流动中表面受力与气动热提供了保障。     

滑移流区化学非平衡流中气动热环境的数值模拟

针对滑移流区化学非平衡流中的气动热环境预测问题，采用多组分化学非平衡纳维-斯托克斯（N-S）方程对轨道验证飞行器（OREX）进行数值模拟，对比了有限催化和非催化、滑移和无滑移壁面条件下的气动加热，研究了壁面催化和滑移效应对气动热的影响规律，并分析了影响气动热的主要机制。结果表明：壁面催化和滑移效应均会影响流场物理量分布，高度为92.82 km时对激波脱体距离的影响较为明显；随着高度减小，有限催化与非催化壁面间的热流偏差增大，而滑移与无滑移壁面间的热流偏差减小；高度低于92.82 km时有限催化壁面计算得到的驻点热流值与飞行数据吻合较好，偏差均在11%以内；机制分析发现，催化效应对壁面附近的组分分布影响较大，滑移边界条件中的温度跳跃条件对壁面附近的温度分布影响显著。     

高超声速再入飞船气体热辐射计算边界虚网格方法

针对再入飞行器激波层辐射加热,应用有限体积法计算辐射传输方程,引入了一种边界处理方式——虚网格方法,根据化学非平衡流场数据,计算出飞行器表面辐射加热热流。通过分析发现:采用虚网格方法处理辐射传输边界,能够确保真正的边界条件得到应用,满足真实的物理意义,且辐射传输计算的空间网格、立体角网格、初始温度值的收敛性很好,计算结果精度误差小于2.5%。利用虚网格方法计算飞船再入返回时的辐射加热,结果表明计算所得壁面辐射热流与其他研究者辐射热流计算结果符合良好。      

发动机燃烧室内壁材料热环境的气动热试验模拟技术研究

利用超声速矩形湍流导管和等离子电弧加热器模拟了发动机燃烧室内流和高超声速飞行器外壁面外流热环境,进行了平板表面冷壁热流测量和燃烧室内壁材料考核试验。结果表明:由于辐射换热的影响,在选取的两个典型来流条件下,发动机燃烧室内流热环境下的冷壁热流比外流热环境下的高出21%和40%,但是冷壁热流的增量基本相当,约为0.70～0.80MW/m2。随着冷壁热流的增加,辐射换热产生的热流增量的影响力会逐渐减小。材料考核时,相同配方的C/SiC复合材料在内流热环境下的表面温度高出约400℃,背面温度高出约90℃,这种差异对于发动机燃烧室内壁面材料考核至关重要,必须在材料考核试验中加以考虑。     

局部催化特性差异对气动热环境影响的计算分析

高温气体非平衡效应及其壁面催化效应对高超声速飞行器气动热环境造成显著影响,是当前高超声速飞行器气动热环境预测和热防护设计的关键问题之一。考虑高温空气离解与电离等化学反应、气体分子热力学激发、流动中的非平衡效应和壁面催化效应,通过数值求解三维热化学非平衡Navier-Stokes方程和壁面处质量、能量平衡关系,完善了高温气体热化学非平衡流场有限催化气动热环境数值计算方法和计算程序,采用典型算例进行了考核验证。在此基础上,开展了不同条件下高超声速飞行器热化学非平衡流场气动热环境数值模拟,分析局部催化特性差异对气动热环境的影响。研究表明：所建立的高超声速飞行器热化学非平衡流场有限催化气动热环境数值计算方法及程序,其数值模拟结果与飞行试验、文献符合;局部催化特性差异会导致热流跳变,其热流跳变量与催化特性差异量、材料分布方式等有关;催化特性差异较大时,局部区域热流可能远远高于飞行器全表面完全催化的热流结果,此时将飞行器在全表面完全催化（FCW）和完全非催化（NCW）条件下的数值模拟结果作为实际飞行过程中表面热流的上、下限这一简化处理方式,是不可取的。     

类乘波体防热壁板气动加热-温度场耦合特性分析

针对类乘波体防热壁板,采用松耦合算法推进气动热和传热迭代计算,研究气动加热-温度场的耦合特性.结果表明:在壁板达到辐射平衡前,冷壁热流和辐射平衡热流与真实气动加热的误差分别达+55.1%和-15.4%以上,必须将气动加热和温度场进行耦合分析;当地绝热壁面温度不随时间变化,表面传热系数是耦合效应的关键参数;采用平均表面传热系数进行瞬态气动加热-温度场耦合计算只进行2次气动加热计算,壁面温度预测误差在2.5%以内,可有效提高气动加热-温度场耦合计算的效率.      

基于PIC/MCC/DSMC方法霍尔推力器热分析

以提高霍尔推力器性能和使用寿命为目的对霍尔推力器进行了热分析。建立了等离子体在通道陶瓷壁面和阳极的能量沉积计算模型,并将模型耦合到PIC/MCC/DSMC流场计算程序中,计算加速通道壁面上的能量沉积分布。将计算得到的壁面能量沉积作为霍尔推力器温度计算的热流边界条件,考虑结构间的热传导与热辐射,计算温度分布。为了考察壁面的热流分布方式不同对温度场的影响,在总热流相同的情况下,将通道热流假定为线性分布与平均分布,计算霍尔推力器的温度分布。结果表明,壁面能量沉积占总功率的20.4%,陶瓷壁面能量沉积沿轴向位置先增大后减小,最大值在电离区。阳极的能量沉积,沿径向在中间位置达到最大,两端靠近壁面处较小。三种热流边界条件下的温度结果表明,三者最高温度都为700K左右,但高温位置不同,PIC热流边界条件下,最高温在电离区,而线性与平均热流边界条件下,高温区分别在出口区与近阳极区。通过与实验结果比较表明,PIC计算热流边界条件下高温区温度与测量者吻合更好,误差小于0.8%。     

发动机舱内部传热和蒙皮的降温规律

建立了简化的发动机舱物理模型,对通风冷却和辐射遮挡的蒙皮降温效果进行了数值模拟和实验研究,并对发动机舱内部传热规律进行了分析.结果表明:数值计算得到的发动机舱蒙皮表面沿程分布与实验结果吻合良好.在加热壁面分段的热边界条件下,无辐射遮挡且无通风时蒙皮内侧表面接受的热流占热壁面加热热流的比值约为78%,当发动机舱通风速度达到1m/s时该比值降至20%;当通风速度从0m/s增大至3.5m/s时,发动机舱蒙皮峰值温度的降幅约为10%;无通风条件下有辐射遮挡的蒙皮峰值温度相对于无辐射遮挡的情形降幅约为6%;在通风冷却和辐射遮挡的综合作用下,蒙皮峰值温度相对于无通风且无辐射遮挡的基准状况可降低约50K,相对降幅约为14%.      

TBCC燃烧室/喷管一体化壁面温度计算

基于Navier Stokes(N-S)方程组对包括隔热屏、隔热屏内外流、大气外流在内的涡轮基组合动力(TBCC)发动机燃烧室/喷管进行了一体化的气/热耦合数值模拟,考虑了燃气组分输运、辐射换热等影响,研究了其在某典型飞行状态下TBCC冲压发动机燃烧室/喷管筒体及隔热屏内外壁壁面温度、辐射换热热流及对流换热热流分布.结果表明:燃烧室/喷管筒体与对称面上下交线的壁面温度在轴向距离为0.5~2.6m内变化较小,在轴向距离为2.6~3.1m内急剧增加,在轴向距离为3.1~3.5m内急剧下降.之后,上交线筒体壁面温度沿流向减小,下交线筒体壁面温度先升高后降低.筒体壁面温度最高点在喷管下调节板收缩段,为1577K.隔热屏内壁面辐射热流在370~500kW/m2变化,上下交线处的辐射热流较外壁面的辐射热流约高300kW/m2,辐射热流沿流向先减小后增加.隔热屏外壁面辐射热流在50~200kW/m2范围内分布.