气膜孔倾角对层板隔热屏冷却性能影响

为了获得气膜孔倾角对层板隔热屏(冲击/发散复合冷却隔热屏)冷却性能的影响规律,基于加力燃烧室真实工况,对0°到90°范围内的十种不同气膜孔倾角的层板隔热屏进行了三维流热耦合数值模拟研究,得到了层板隔热屏冲击壁面Nu数、层板隔热屏气膜冷却表面的冷却效果、层板隔热屏冷流体热负荷及气膜孔流量系数的变化规律。结果表明,气膜孔倾角的变化对冲击壁面Nu数的影响较小;气膜冷却表面的综合冷却效果随气膜孔倾角的增大而减小,15°倾角模型比10°倾角模型的平均综合冷却效果降低2.8%;单位面积冷流体热负荷随气膜孔倾角的增大而增大,最小值比最大值低30.7%;气膜孔倾角对层板隔热屏平均流量系数的影响不大,但上游气膜孔的出流会对下游气膜孔的流量系数产生影响。     

冲击/发散冷却层板隔热屏冷却性能及对比

为了分析对比新型冲击/发散冷却层板隔热屏冷却性能,论证其应用于加力燃烧室的可行性,在3种不同主次流总压比条件下对其进行了三维流固耦合传热数值模拟研究,并与某型波纹板隔热屏和单层平板隔热屏进行了相同工况的对比分析,得到了冷却效果、冷气用量、冷气热负荷和次流总压损失系数等的对比结果和变化规律.结果表明:冲击/发散冷却层板隔热屏具有较好的冷却效果,但其受主次流总压比变化的影响较大.相比某型波纹板,冲击/发散冷却层板隔热屏的冷气消耗量平均减少41.6%,单位面积冷气热负荷平均降低65.9%.     

流向/周向曲率对加力尾部双层壁隔热屏冷却性能影响

为了得到流向/周向曲率对双层壁隔热屏冷却性能的影响规律,基于加力燃烧室尾部圆转方双层壁隔热屏弯曲型面的工况条件,建立了两组共6种不同流向曲率半径和5种不同周向曲率半径的双层壁隔热屏冷却模型,采用三维流热耦合数值模拟方法对流向和周向曲率半径对双层壁隔热屏冷却性能的影响进行了研究,得到了双层壁隔热屏气膜冷却表面的综合冷却效果、冲击靶面的Nusselt数Nu、气膜冷却表面的单位面积冷却气量等的变化规律。结果表明,随着正流向曲率半径的增大,综合冷却效果先增大后减小,变化幅度达3.18%,冲击靶面Nu先减小后增大,单位面积冷却气量增大,增大幅度达9.78%;随着负流向曲率半径的增大,综合冷却效果先减小后增大,变化幅度为6.56%,冲击靶面Nu先减小后增大,单位面积冷却气量减小,减小幅度达15.82%;随周向曲率半径的增大,气膜冷却表面的综合冷却效果先减小后增大,变化幅度为1.04%,冲击靶面Nu减小,单位面积冷却气量减小,减小幅度为3.48%;相比于周向曲率半径对双层壁隔热屏冷却性能的影响程度,正流向曲率半径的影响程度为3倍左右,负流向曲率半径的影响程度为6倍左右。相比于周向曲率半径,流向曲率半径的变化对双层壁隔热屏冷却性能的影响更强烈。     

发散孔纵向波纹隔热屏气膜冷却特性

对燃烧室内开有发散孔的纵向波纹隔热屏进行了数值模拟.研究隔热屏的四种结构参数开孔率、波纹板高度、孔径和冷却通道高度的改变对隔热屏冷却效果的影响.研究表明:在气膜孔出流总量相同的情况下,3%开孔率比6%开孔率的隔热屏平均冷却效率较高;波纹板高度对隔热屏冷却效果影响较大,波纹板无量纲高度为1%比波纹板无量纲高度为2%和3.33%的隔热屏平均冷却效率高;孔径和冷却通道高度的改变对隔热屏冷却效果几乎没有影响.      

发散孔结构参数对横向波纹表面气膜绝热冷却效率的影响

针对加力燃烧室特定横向波纹结构的隔热屏发散冷却进行数值模拟,主要研究发散孔孔节距、孔径以及开孔率等3个结构参数对发散冷却效率的影响。结果表明:由于处于波纹波谷的相邻气膜射流更易于形成相干,从而在波谷形成更强的集聚效应,造成波谷附近的绝热壁面温度低于波峰区域;在相同的壁面单位面积冷气用量下,减小孔径、增大开孔率均显著改善气膜冷却效率,尤其是在发散气膜的前排起始段;发散孔流向节距大于展向节距的长菱形排布相对较优,但在小吹风比下发散孔排布节距比对展向平均绝热冷却效率的影响非常微弱。      

扫掠冲击-气膜冷却结构的气动传热特性数值研究

为了揭示“扫掠冲击-气膜”冷却结构的换热机理,采用气热耦合方法和SST k-ω湍流模型,对比分析了吹风比为1, 2, 3, 4和气膜孔角度为30°, 45°, 55°, 65°等条件下“直接冲击-气膜”组合方式和“扫掠冲击-气膜”组合方式在平板模型上的气动传热特性。结果表明,流体激振器的扫掠频率、冲击靶面上的Nu数随吹风比增大而增大,并且几乎不受气膜孔角度影响。两种组合方式的总压损失系数和综合冷却效率随吹风比增大而增大,并且随气膜孔角度的增大而略微减小。尽管在使用相同冷气流量时“扫掠冲击-气膜”组合方式的冷气进口静压较高,但是其具有冲击靶面上Nu数分布均匀、综合冷却效率更高且分布面积更大的优势。     

气膜孔与冲击孔面积比和动量比对加力燃烧室双层壁隔热屏综合冷却效率的影响

基于综合冷却效率模化理论和匹配原则,采用红外测温技术测量了加力燃烧室双层壁隔热屏的综合冷却效率分布,分析了气膜孔与冲击孔面积比( 1、2、3、4)和动量比(I=0.02-0.88)等参数对加力燃烧室双层壁隔热屏冷却特性的影响规律。研究结果表明：综合冷却效率分布由冲击冷却、气膜冷却及冲击孔和气膜孔的相对位置决定。高冷效区域集中在冲击驻点附近和气膜覆盖的区域。沿着主流的流动方向,气膜板下游的综合冷却效率高于上游的综合冷却效率。随着动量比的增大,综合冷却效率增大。随着气膜孔与冲击孔面积比的增大,气膜板上游的综合冷却效率主要是由于冲击换热的增强而增大,气膜板中下游的综合冷却效率是由于冲击换热和气膜冷却效果的增强而增大。气膜孔与冲击孔面积比由1增大到3,面平均综合冷却效率提高72%,气膜孔与冲击孔面积比由3增大到4时,面平均综合冷却效率提高了12%。     

开孔率对加力燃烧室隔热屏冷却特性的影响

为了分析开孔率对航空发动机加力燃烧室隔热屏流动和冷却特性的影响,建立了双层壁隔热屏结构。在给定冲击距、 发散板与冲击板开孔面积比的条件下,对开孔率为0.6%和0.8%的冲击发散冷却结构进行数值模拟,获得了发动机真实工作状态 下的气动参数、几何参数对流量系数和综合冷却效率的影响规律。结果表明：冲击孔、发散孔流量系数沿主流方向略微增大,随着 开孔率由0.6%增大到0.8%,不同吹风比下冲击孔的流量系数均减小,而气膜孔的流量系数变化不明显,在吹风比由0.2增加到0.8 时,0.6%开孔率结构的综合流量系数增幅明显高于0.8%的;在吹风比为0.2时,0.6%开孔率结构上游区综合冷效优于0.8%开孔率 结构的,而在下游区则正好相反,在吹风比为0.3～0.8时,0.8%开孔率结构各区域的综合冷效均高于0.6%开孔率结构的。     

纵向波纹隔热屏气膜冷却特性数值研究

以加力燃烧室纵向波纹隔热屏为研究对象,采用数值方法,对不同吹风比、不同展向间距以及不同流向间距的纵向波纹隔热屏进行计算,获得了纵向波纹隔热屏流场流动特征和气膜冷却效率的变化。结果表明:纵向波纹隔热屏气膜冷却与常规的平板气膜冷却有着本质的区别,表现出气膜射流脱离壁面、较小吹风比下的冷却效率受结构影响大等特征;较低吹风比(M1.3)下气膜冷却效率沿着流向呈现出起伏变化,较高吹风比(M≥1.3)下的气膜冷却效率沿程逐渐增加,最终趋于平缓;对于几何结构,吹风比约为2.5时效率达到最大值;随着展向间距的减小,气膜冷却效率逐渐增加,但增加较为缓慢;气膜冷却效率并不是一直随流向间距的减小而增加,当流向间距从4.2倍孔径减小到3.25倍孔径时,在第二周期波峰区域,冷却效率反而降低,随着流向间距进一步减小,局部冷却效率降低的区域甚至扩大,并占据了波纹隔热屏的第二个波峰。     

弯曲段壁面冲击发散冷却效率的实验

为了研究回流燃烧室弯曲段采用冲击发散冷却结构形式对冷却效率的影响规律,设计了多种不同几何尺寸的计算模型,采用试验方式分别对其冷却效果进行了试验研究,得出了如下结论：（1）吹风比对冷却效率影响显著,随着吹风比的增大冷却效率升高。沿整个弯曲段冷却效率呈现先增加后减小的趋势,但是变化幅度很小。（2）弯曲段冲击发散冷却结构发散孔倾角对冷却效率影响很大,冷却效率随着发散孔倾角的减小而增大。（3）发散孔纵向间距小的发散孔板冷却效率高于纵向间距大的发散孔板的冷却效率。     

带冲击的两相气膜冷却特性数值模拟

用数值模拟方法研究了带冲击的气膜冷却系统的流场结构，示出了气流和冲击腔内的螺旋流动、气膜孔下游的肺状低压区、横截而对称双漩涡以及由双漩涡形成的交叉双锥形水蒸气体积分数等值面。比较了体积分数为0．5的水蒸气膜与纯空气膜的冷却效率，得到水蒸气膜平均冷却效率高于空气膜。在x=3d处高出31．6％。在x=26d处高山19％。水蒸气膜与空气膜相比换热系数变化不大。     

冲击距与气膜孔方位角对旋流气膜冷却性能影响

为了获得冲击距Hi与气膜孔方位角α对旋流气膜冷却性能的影响规律,以六边形供气腔圆形气膜孔平板气膜冷却结构为研究对象,对五种冲击距参数(0.74D,1.14D,1.54D,1.94D,2.34D)(D为气膜孔直径)和五种气膜孔方位角参数(0°,10°,15°,20°,25°)进行了三维数值计算研究,得到了绝热壁面气膜冷却效率、展向平均气膜冷却效率、流场空间无量纲浓度分布等随冲击距与气膜孔方位角的变化规律,分析了肾形涡对旋流气膜冷却性能的影响机理。结果表明:冲击距对绝热壁面气膜冷却效率展向分布规律影响不大,而方位角增加能够明显提高绝热壁面气膜冷却效率及展向气膜覆盖面积,方位角0°模型展向气膜冷却效率最大值为0.42,方位角25°模型的最大值为0.48,相比前者增加14.3%;绝热壁面同一流向位置的展向平均气膜冷却效率随冲击距的增加而增大,随方位角的增加而增大,Hi=2.34D时的展向平均气膜冷却效率最佳,α=20°时的展向平均气膜冷却效率曲线最佳。方位角的增加能够明显破坏流场中存在的肾形涡结构。     

叶片前缘两相流冲击冷却的耦合数值模拟

为了更好地强化叶片前缘的换热,对叶片前缘进行两相流强化冲击冷却的耦合数值模拟研究。利用水滴/蒸汽冲击曲面冷却系统进行数值模拟程序的校核。在此基础上,利用欧拉-拉格朗日颗粒追踪方法对叶片前缘的两相流冲击冷却进行耦合数值研究。研究了水滴加湿量、水滴直径和水滴与壁面的作用边界条件对换热性能的影响。结果表明当加湿量从3%增加到8%时,最大的平均传热增强系数从1.671增加到4.913;水滴直径从5μm增加到20μm时,最大的平均传热增强系数从2.128降低到1.164;反弹和反弹+破碎边界条件预测的传热效果好于粘附和粘附+破碎边界条件。     

水滴撞击特性的重力影响分析

以NACA0012翼型为研究对象,采用拉格朗日法分析了重力对水滴撞击特性的影响,获得了过冷大水滴（SLD）的撞击特性。计算结果表明：当水滴直径超过50μm时,重力对水滴轨迹以及极限轨迹的初始位置影响很大;SLD的水收集系数和收集区域都较直径为10～40μm的水滴（美国联邦航空条例FAR25部附录C规定的防冰系统设计用水滴直径范围）有较大幅度的提高,该结论可为飞机防冰系统的设计以及冰风洞中喷嘴安装位置的确定提供参考。     

具有微小W型肋的结构化表面冲击冷却实验

将冲击冷却技术与肋化表面相结合，研究了一种具有微小W型肋的表面射流冲击冷却结构。通过稳态实验和瞬态热色液晶测试技术来探究光滑靶板和微小W型肋靶板的传热特性。测试时的冲击间距比为1.5，3和5，基于水利直径的雷诺数为15000~40000。结果表明:两种靶板的平均努塞尔数和压力损失均随雷诺数的增加而增加，随冲击间距比的增加而减小。当冲击间距比为1.5时，与光滑靶板相比，微小W型肋靶板的平均努塞尔数提高了5.1%~7.3%，压力损失却几乎不变。但当冲击间距比大于3时，由微小W型肋带来的强化传热效果并不显著。      

基于SPH方法的射流撞击仿真

为了研究液体推进剂的撞击雾化问题,采取光滑粒子流体动力学方法 (SPH)进行了射流撞击形成液膜过程的数值模拟。修正了基于连续表面力模型的表面张力算法,推导了考虑粘性与表面张力作用的SPH形式的控制方程,对撞击雾化模型进行了合理简化,采用水作为推进剂模拟液。得到了不同撞击速度、角度下的液膜图像及液膜的形成过程图像。证明了SPH方法适于解决存在自由表面及大变形的射流撞击问题。     

非均匀横流作用下冲击射流冷却的数值研究

运用数值计算的方法对非均匀横流作用下冲击射流冷却的流动和传热过程进行了三维数值研究,揭示出射流冲击孔和横流孔相对排列方式、射流雷诺数和横流雷诺数相对大小、射流孔中心至横流孔出口截面距离等因素对冲击射流冷却流动传热特性的影响。结果表明:若保持总的冷却流量不变,当横流雷诺数与射流雷诺数之比等于或大于0.5后,有横流时的冲击射流冷却局部努塞尔数比无横流时有较为显著的下降;在非均匀的横流作用下,对于叉排和顺排方式,由于射流和横流的相互作用,壁面温度分布呈现出复杂的特征;射流孔和横流孔呈顺排布置时,峰值努塞尔数所对应的射流驻点区向下游偏移。      

过冷水滴撞击结冰表面的数值模拟

为了计算过冷水滴撞击结冰表面,本文应用雷诺平均Navier-Stokes方程和k-ε两方程紊流模型计算结冰表面外的空气流场。用欧拉法建立了空气中过冷水滴的运动方程,并采用有限控制容积法对方程进行了数值计算。为了提高计算精度,方程中的对流—扩散项采用MUSCL格式进行了离散。计算结果表明过冷水滴的大小和机翼攻角的变化对结冰区的大小和水收集系数的影响很大。      

飞机水平尾翼水滴撞击特性及防冰热载荷计算

建立了利用欧拉法求解水滴撞击特性的方法,并基于Messinger质量和能量守恒方程建立了热载荷计算模型,以某型飞机三维水平尾翼为研究对象,展开了撞击特性和热载荷的计算.结果显示水平尾翼的局部水收集系数极值从翼根到翼尖逐渐增加,不同截面水滴收集系数的分布为设计防冰系统范围提供了依据.驻点附近及机翼上、下表面各存在一个热载荷较大的区域,进而确定该尾翼电加热防冰分为3个区域:中间连续加热区和上、下表面驻点附近的各一个间断加热区.