单喷嘴燃烧室传热特性的测量方法

为开展燃烧室的传热特性等研究,采用了以紫铜作为材料的容热式推力室身部,通过测量热试车过程中身部壁温,来获得推力室内壁温度、热流沿轴向的分布。为测量数据设计了数据处理方法,以利用测得的温度值获得内壁温和热流值,且进行了主要误差估计。将该方法在气气喷注器推力室上进行了应用,研究了其传热特性。     

壁温对压缩拐角流动影响的数值模拟研究

压缩拐角是高超声速飞行器上的典型非连续区域,其分离/再附结构对局部热环境有较大影响。本文采用基于Navier-Stokes方程的自研程序,开展了典型压缩拐角外形的气动热环境数值模拟研究,分析了不同壁温条件下压缩拐角的热环境分布规律。计算结果表明:随着壁温升高,流场整体结构变化不大,但由于近壁面流体物性参数变化较大,导致拐角分离涡分离点向前移动、再附点向后移动,拐角处干扰区扩大;压缩拐角大部分区域的热流随着壁温的升高而减小,且干扰区内热流减小的幅度比无干扰区更大,但热流并不完全遵循随壁温升高而减小的规律。另外,通过对比变壁温计算热流和热壁修正公式修正热流,发现热壁修正公式在干扰区存在精度降低、适用性不足的问题。本文的研究进一步加深了壁温对压缩拐角流动影响的认识,缩小了热壁修正公式的适用范围。     

航空发动机试验峰值时间对示温漆测温的影响

采用示温漆和热电偶两种方式,在全温全压试验中对航空发动机阵列预混燃烧室火焰筒壁温进行测试,获得火焰筒的壁温分布,并对不同峰值时间下示温漆判读结果与热电偶测试结果进行对比分析。研究表明:阵列预混燃烧室火焰筒壁温小于所用材料的许用工作温度;火焰筒试验峰值时间是影响示温漆测温判读的主要原因,为保证测试结果的准确性,应严格要求试验峰值时间与示温漆标定的峰值时间一致。     

内流热环境的热流模拟方法

对内流热环境的热流测量方法进行了试验研究和模拟,对比分析了内流热环境与外流热环境的不同.采用超声速矩形湍流导管试验技术,利用等离子电孤加热器进行了内流与外流热环境的表面冷壁热流测量和材料考核试验.结果表明:在来流条件相同的情况下,与外流热环境下的测量值相比较,内流热环境下的冷壁热流值高出40%多,相同材料的背面温度高出400℃.通过直接测量的冷壁热流换算得到的热壁热流值却不到外流热环境下的热壁热流值的50%.因此,进行发动机防热设计时,必须进行材料考核试验,必须考虑材料的使用热环境,否则对材料的烧蚀性能的评估影响很大.     

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

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

传感器表面温度对热流测量的影响

高超声速飞行器热环境测量数据一直是防热设计和考核的基准数据,也是热环境计算方法的考核性数据,其准确性至关重要。针对热流测量中遇到的传感器表面和周边防热材料温度差异而导致的测量数据偏差问题,采用基于Navier-Stokes方程的自研程序开展了详细的气动热环境数值模拟,得到了不同温差条件下传感器表面热环境的分布规律,并根据场协同理论分析了局部热流变化的成因机理,研究了影响热流变化幅值的主要因素。结果表明：①当传感器和其周边材料的温度存在一定的差异时,导致该区域近壁面流场中的压力、密度等特征量梯度增大,改变了传感器当地的法向速度和温度分布,造成了局部热流的剧烈变化。②相同来流马赫数和高度下,来流攻角主要影响法向速度的分布,从而影响气动加热量,攻角越大,相同温差下加热量上升的幅度越小;来流总温主要影响法向温度梯度的分布,从而影响气动加热量,来流总温越大在相同温差条件下加热量上升幅度越小。所开展的研究工作可加深对传感器局部热环境分布规律的认识,避免对测量热流的误判,提升数据判断和分析的可靠性。     

激波风洞温敏热图技术初步试验研究

介绍了可用于激波风洞等高超声速脉冲设备复杂外形模型表面热流分布显示与测量的温敏热图技术。该技术利用温敏材料发光特性随温度变化的特点实现对模型表面的温度测量,进而分析获得其表面热流分布。由于有效运行时间一般只有10ms左右,在激波风洞中采用温敏热图技术在材料温敏响应、图像采集等方面均存在困难,导致此技术还很不成熟。介绍了近期在CARDC 0.6m激波风洞中开展温敏发光热图技术研究的情况,包括发光材料的研制、材料喷涂、图像处理及标定等方面的研究工作。验证试验采用了平板加钝舵模型,获得了与理论分析及传感器测量结果较一致的平板干扰区热图结果。热图与薄膜传感器测量数据差异大致在5%~25%之间。     

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

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

氢氧烃三组元发动机燃烧室传热、流阻初步分析

本文对几种发动机的冷却作了分析比较,指出使用铜锆合金是高压大热流发动机的必然要求.对大推力氢氧烃三组元试验发动机燃烧室的传热分析表明:热流仅相当于相同推力氢氧双组元发动机的1/3～1/2,热壁温低于800K,完全满足设计要求,但流阻较大,对氢泵不利.     

CARDC 激波风洞 TSP 技术研究进展

从基本原理、关键技术和验证应用三个方面总结了近两年在中国空气动力研究与发展中心激波风洞中开展的温敏涂层(TSP)技术相关研究工作。通过解决快速响应温敏发光材料研制、模型研制、数据处理等一系列关键技术,完成图像采集系统、光学系统及标定系统的配套和系统集成,建立了一套适于激波风洞试验的高速 TSP 测量及标定系统。该技术可在激波风洞试验中获取模型被测面温敏涂层的发光图像,基于该图像可以直接观察模型表面热流分布和捕捉峰值热流的准确位置。结合温敏发光材料的物性参数标定数据,能够实现对模型表面热流的定量测量。不同于传统的传感器点热流测量技术只能得到模型表面有限数量的离散点的热流值,TSP 技术能够以高空间分辨率得到较大面积区域的详细热流分布信息,可更加全面的测量模型外表面的热环境,并且可以据此进一步分析和辨别边界层流态以及确定边界层转捩位置。试验对比表明,TSP 技术的测量结果与点热流传感器的测量结果具有良好的一致性。目前该技术已趋于成熟,在Φ2 m 和Φ0.6 m 激波风洞上成功应用于边界层转捩研究、局部干扰区热环境研究和复杂外形飞行器热环境研究等领域,已成为激波风洞除点测热技术之外又一重要测热技术。     

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

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

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

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

燃烧室辐射换热一维区域法模型研究

目前用于工程计算燃烧室辐射热流的Lefebvre零维辐射模型,不考虑燃气区域对附近壁面的辐射影响,过于简化,不能反映实际燃烧室火焰筒内的辐射换热过程。本文发展了一种机理严密,精度高的用于燃烧室壁温计算的区域法辐射换热模型。文中采用Lefebvre模型和区域法模型计算了一个圆筒型不冷却火焰筒结构的换热;对比分析了两种方法计算的壁面辐射热流和壁面温度分布情况。结果表明,使用Lefebvre模型计算时,由于将燃气按一维划分为定向辐射,在燃烧室前后壁面辐射热流较小,存在一定不足;而使用区域法模型计算时,辐射热流分布更趋合理。     

带主动冷却的超声速燃烧室传热分析

介绍了流体、固体传热耦合的一维分析方法,对带主动冷却系统的马赫数2.5超声速燃烧室进行了传热分析。该分析以实验测量的燃烧室壁面静压以及超临界煤油换热特性数据为基础,考虑了燃气的高温离解效应,燃烧特性以及碳氢燃料的高温热物理特性,对不同燃烧状态、冷却条件下的主动冷却过程进行了分析。结果表明有燃烧时壁面热流可高达1MW/m2以上,是无燃烧时的2～3倍。当煤油流量较小时(当量比为0.45),冷却后的壁面温度仍偏高,而且冷却壁内温度分布不均匀。随着煤油流量的增加,冷却效果明显提高,冷却壁内温度分布趋于均匀;并且煤油的出口温度也显著减小。     

运用区域法模型计算燃烧室内一维辐射换热

采用一种机理严密,精度高的区域法辐射换热模型用于燃烧室壁温计算.采用Lefebvre模型和区域法模型计算了一个圆筒型不冷却火焰筒结构的换热;对比分析了两种方法计算的壁面辐射热流密度和壁面温度分布情况.结果表明,使用Lefebvre模型计算时,由于将燃气按一维划分为定向辐射,在燃烧室前后壁面辐射热流密度较小,存在一定不足;而使用区域法模型计算时,辐射热流密度分布更趋合理.      

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范围内分布.      

高温圆管内多孔介质耦合换热的辐射作用

通过数值模拟圆管内填充多孔介质的辐射对流耦合换热,研究高温固体骨架辐射效应对温度分布及换热的影响.基于局部非热平衡模型分别建立流体和固相能量方程,采用蒙特卡罗法求解固体骨架的辐射换热;对不考虑热辐射引起的温度场偏差和管壁发射率以及多孔结构参数的影响进行讨论.结果表明:多孔固体骨架的辐射效应对入口段温度场的影响明显,不考虑辐射将导致较大偏差,壁温为1500K时最大偏差为16%.管壁发射率对温度场影响较小,壁温为1500K时影响小于3%.孔隙率或孔径增大,壁面辐射热流密度比例增加,辐射效应体现明显.      

低速高温燃气流热模拟试验方法和设备

对比分析了两种气流状态参数和两种加热情况下典型前缘部件表面热流密度的相似性,论证了利用亚声速高温燃气流加热方式进行近地空间高超声速飞行工况气动热模拟试验的可行性.针对高超声速飞行器典型钝头锥结构提出“小喷口低速高温燃气流+石英灯”组合热试验方案.通过采用新型高效双腔蒸发管型燃气发生器、新型带保温夹层和耐高温陶瓷内衬的水冷不锈钢高温管道结构,同时引入电加热器预热及燃烧室两路供油方案,使所建低速高温燃气流热试验设备产生燃气流温度达到2100K,φ250mm喷口处平均径向温度分布梯度约3K/mm,具有线性温度控制功能且稳态控制温差约46K,满足24km、马赫数为6典型高超声速飞行器工况驻点区域高温/大热流密度气动热试验要求.      

典型尖楔结构燃气流热试验工况确定方法

在传统冷壁热流模拟方法的基础上,进一步提出以热壁温度及热流密度的时序变化曲线为控制目标的燃气流热试验工况确定方法,即利用壁温控制目标与实测值的偏差对热壁热流控制目标做一定修正,以尽可能消除和弥补前期试验误差,同时利用300K冷壁边界热流密度数据库插值迭代方法,快速确定一定气动热模拟所需燃气流温度,解决了沿飞行轨迹瞬态热试验技术难题之一。利用CFD数值模拟方法,建立了典型尖楔结构高/中温双路燃气流组合热试验300K冷壁边界热流密度数据库,并针对典型尖楔结构沿某飞行轨迹9个典型状态气动热模拟需求,确定相应双路燃气流热模拟参数。相关数值计算结果显示,驻点区域热流密度平均模拟偏差为4.5%,平板区热流密度平均模拟偏差为4.6%,两者最大模拟偏差均不大于8%,满足工程试验精度要求。同时,瞬态热分析结果显示第45s时,距驻点1mm处最大温度梯度达到21K/mm,距驻点10.1mm处最大温度梯度达到18K/mm,满足气动热大温度梯度效应需求。