旋转爆震发动机燃烧室壁面烧蚀热防护研究

为了研究碳化烧蚀材料对旋转爆震燃烧室的被动热防护作用效果,建立了旋转爆震燃烧室热环境计算模型和多层复合壁一维烧蚀热响应分析模型。基于旋转爆震燃烧过程的数值模拟,确定了沿燃烧室轴向的壁面平均温度与热流的分布特征;针对典型的壁面热负荷输入条件,采用动边界隐式差分格式,对碳化型烧蚀材料的烧蚀过程和传热过程进行了耦合计算,分析了碳化层剥蚀、热解热、热解气体质量流率、材料厚度等因素对燃烧室壁面热防护效果的影响。研究结果表明,旋转爆震燃烧室温度和热流密度分布不均匀,斜激波扫掠区域热负荷更为剧烈;碳化层剥蚀后会使烧蚀材料的烧蚀率增加,壁面温度迅速上升;增大热解热及热解气体质量流率都有利于延长工作时间、降低壁面温度; 旋转爆震燃烧室不同区域的烧蚀有一定差异,增加烧蚀材料厚度无疑有利于壁面热防护,但是在实际应用中应综合考虑壁面重量和空间尺寸。     

液体火箭发动机燃烧室壁面热流测量方法研究

为了研究液体火箭发动机燃烧室的壁面热载荷,提出了一种简单、准确的燃烧室壁面热流测量方法。通过发动机热试,测得了距燃烧室壁面不同径向深度的瞬态温度数据,采用一种基于分离变量的数据处理方法,最终得到燃烧室壁面热流和壁面温度。研究表明,提出的热流测量方法能够准确获得燃烧室壁面热流和壁面温度且不会放大误差。发动机热试的瞬态过程中,温度的响应存在明显的迟滞性,5s时间不足以使温度达到稳态。而燃烧室壁面热流响应迅速,达到稳态所需时间在2 s左右,且热流的变化规律与燃烧室压力变化规律一致。     

旋转爆震燃烧室梯度复合热防护结构热分析模型及验证

针对旋转爆震燃烧室高热流密度的热防护需求,提出了一种碳化硅耐烧蚀层-高硅氧烧蚀层-气凝胶隔热层-不锈钢金属基体层的梯度复合热防护结构,建立了考虑烧蚀的多层平壁一维瞬态热分析模型,结合旋转爆震燃烧室的典型热环境,采用动边界隐式差分计算格式求解获得了壁面输入热流密度和高硅氧烧蚀层主要参数对热防护结构内部温度分布的影响;同时开展了旋转爆震燃烧室梯度复合热防护结构热考核试验,对热分析模型进行了试验验证。研究结果表明：旋转爆震燃烧室壁面沿轴向温度分布存在着内在的不均匀性,前端由于预混气的及时补充而得到有效冷却,温度峰值出现在位于中截面和燃烧室出口的尾端区域。基于径向一维传热简化,从实测的旋转爆震燃烧室壁面轴向温度分布反演出时均热流密度沿程分布,并以此对特定轴向截面的热防护结构温度瞬态变化进行了分析,与试验结果的对比验证了所建立地考虑烧蚀过程的一维瞬态热分析模型可以较好地预测梯度复合热防护结构的温度变化特性。     

气氢/气氧同轴剪切单喷嘴燃烧室壁面热流计算与分析

以气氢/气氧为推进剂,采用数值模拟方法,研究了同轴剪切喷嘴设计参数——氢氧速度比和氧压降比对单喷嘴燃烧室内燃烧过程和壁面热载的影响,并将绝热壁面条件、等温壁面条件的计算结果与试验结果进行了对比分析,结果表明:氢氧速度比增大,燃烧性能提高,壁面热载增加;氧压降比增大,燃烧性能下降,壁面热载减小;相比采用壁面绝热燃气温度,采用热流预示燃烧室壁面热载与真实情况更为接近.      

气气喷注器推力室传热特性研究

设计了单喷嘴气气喷注器容热式推力室,进行了0.92～6.1 MPa范围内7个燃烧室压力工况,共17次热试车;采用壁面测温方法获得推力室沿轴向方向的内壁热流分布,得到各室压工况下燃烧室内壁热流分布曲线,不同工况的热流曲线显示出相似的分布;并从中分离出对流传热热流,得到对流传热热流与室压的关系。为拓展范围,采用多组分湍流N-S方程描述推力室内燃烧流动,采用6组分9步反应模型来描述氢氧反应,反应速率由Arrhenius公式计算,进行了5～20 MPa更高室压范围内的燃烧内流场的数值模拟,并耦合计算了各工况燃气与室壁之间的传热,获得了与试验相同规律的结果。     

热容燃烧室基于单点温度的热流测量方法研究 *

为了研究推力室内壁面与燃气间的传热特性,获得多喷嘴燃烧室内壁热载分布,从而为推力室设计及热结构分析提供参考,开展了热容燃烧室壁面热流测量方法研究。以位于距内壁一定距离位置的单点测量温度作为输入条件,文章提出了三种热流计算方法,并进行火箭发动机试验研究,得到了对应位置燃气对内壁的热流密度。可以看出,运用单点方法得到的热流密度与传统两点法结果符合较好,两者计算结果误差在10%以内。结果显示:不同点火时长工况下在同一位置相同时间点上重复性很好,在燃烧室圆柱段热流密度随着点火时长的增加而减小,而在喷管段变化规律则相反,热流密度随时间增加而增加。在提出的计算平均热流的三种方法中,热积累方法的计算结果最高,热平衡方法次之,瞬态法最低。在火箭发动机热容试验中,本文提出的方法可以用于测量燃烧室内壁面热流密度。     

三级压缩楔前缘半径对壁面静压及热流分布影响研究

为了研究三级压缩楔前缘半径对壁面静压和热流分布的影响,在0.6m激波风洞上开展了试验测量,模型长约0.6m前缘半径为0mm3mm试验名义马赫数为5.98。研究结果表明,试验得到的壁面静压和热流数据重复性很好。采用Fluent软件进行了二维和三维流场参数辅助模拟分析,三种不同湍流模型获得的壁面静压分布差别不大,均与试验结果吻合较好;不同湍流模型获得的壁面热流分布差异较大,采用标准k-ε模型得到的结果与试验吻合较好。试验和数值模拟结果均表明,在第二和第三压缩面上,经过激波后,壁面静压逐渐上升到一个压强平台;壁面热流逐渐上升到一个局部极大值,然后在同一压缩面持续下降。随前缘半径增加,壁面静压和壁面热流整体减小,压强平台值和热流局部极大值也减小,而达到压强平台和热流局部极大值需要的长度增加,显示激波边界层干扰影响区域增大。     

基于传热反问题方法的N_2O/C_2H_4预混推进剂燃烧室热流测量研究

测量液体火箭发动机的热载荷是获取燃烧室内部信息的重要方法。为了获取N_2O/C_2H_4预混推进剂燃烧室内壁的热载荷,建立了液体火箭发动机的热流计算的反问题方法,该方法基于对燃烧室壁面温度场的直接求解,通过对轴向多个位置测量温度的反演计算得到燃烧室内壁热流和温度。研究表明:应用文中建立的传热反问题方法能够较为准确地获得热流随时间及空间的分布;热电偶的位置对计算准确性有明显的影响,与理论深度偏差在0.2mm以内的随机深度偏差可导致超过4%热流反演误差;N_2O/C_2H_4预混推进剂燃烧室热流及温度沿轴向逐渐降低,表明燃烧室内的反应释热过程主要在燃烧室头部附近发生。     

热沉式燃烧室热载分析

采用了以紫铜为材料的热沉式燃烧室身部设计,测量了不同氢/氧速度比下同轴剪切气-气喷注器燃烧室壁面的温度分布和燃烧室压力,分析了不同速度比下推进剂的特征速度效率和燃烧室壁面温度分布的不同,计算了沿燃烧室轴向的热流强度和燃烧室吸热量占推进剂完全反应所应释放能量的比例分数.研究方法对工程应用具有参考价值.     

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

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