发动机旋转整流帽罩的水滴撞击特性分析

基于计算流体力学软件FLUENT,采用拉格朗日法,提出了一种求取旋转部件表面水滴收集系数的方法,并采用用户自定义函（UDF）予以实现。将该方法运用于发动机旋转整流帽罩的水滴撞击特性研究,结果表明：随着帽罩锥体锥角及转速的增大,帽罩表面的局部水收集系数沿其表面外形线下降较快,且撞击区域变小。     

旋转帽罩电加热防冰瞬态过程研究

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟,提出了旋转部件表面对流换热计算方法,并基于改进Messinger结冰模型开发了旋转帽罩电热防冰计算程序,对旋转帽罩瞬态防冰过程进行了数值模拟研究.结果表明:旋转运动对旋转帽罩表面对流换热起到增强作用,且在供给相同加热热流密度时转速越大,防冰表面温度越低;防冰系统启动阶段,旋转帽罩表面会发生结冰和冰脱落现象,考虑结冰过程后系统的响应时间缩短;当电加热功率相同时,周期电加热防冰方式更为节能;当电加热能耗相同时,周期电加热方式系统响应更快.     

发动机旋转表面水滴撞击特性数值研究

介绍了发动机旋转机械内部过冷水滴轨迹的三维数值计算方法,阐述了旋转坐标系下气流及粒子的运动规律。采用ANSYS-CFX软件及其粒子输运模型,对某发动机风扇转子叶片外围空气及水滴流场进行了数值模拟。利用水滴的速度矢量、撞击区域等参数表征水滴的撞击特性,获得了转速、水滴直径等对风扇叶片表面水滴撞击特性的影响:水滴撞击区域集中在风扇叶片迎风面叶盆侧,且水滴撞击区域随着转速的增加而减小;水滴在叶片表面的撞击范围随着水滴直径的增大而减小。     

Roe格式在水滴欧拉方程数值求解中的应用

目前普遍采用欧拉法求解水滴运动控制方程以模拟结冰表面水滴撞击特性。水滴欧拉方程近似没有压力项的空气欧拉方程,但已应用于空气欧拉方程的空间离散格式多数不能直接应用于水滴欧拉方程,有研究者将经典Roe格式应用于水滴欧拉方程时出现了奇点问题。通过对水滴欧拉方程进行形式改造,基于中心格式验证改造方法的有效性,并应用Roe格式结合线性重构与限制器函数实现对水滴欧拉方程二阶精度的求解。将NACA23012翼型计算结果与FENSAP-ICE软件进行对比,两者水滴收集系数分布基本一致。     

调节角度对整流支板水滴撞击特性的影响

基于欧拉-拉格朗日法计算空气-水滴两相流场,获得了表面水滴撞击特性及撞击量,并对可调整流支板(调节角度0°和30°)的水滴撞击特性进行研究。计算结果表明:无论偏折角为0°还是30°,支板前缘处都有水滴收集,且沿着弦长方向局部收集率呈越来越小的趋势;而30°支板在偏折处也有水滴收集。同时分析了马赫数、液态水含量、水滴直径对水滴撞击特性的影响,随着来流马赫数、水滴直径的增大,水滴收集系数、撞击水量和撞击极限都是增加的;而随着液态水含量的增加,水滴收集系数和水滴撞击极限是不变的,而撞击水量却是增加的。     

考虑遮蔽区影响的旋翼三维水滴撞击特性计算新方法

针对直升机旋翼三维黏性流场特有的复杂环境,建立了一种基于欧拉法的旋翼三维水滴撞击特性计算的新方法。首先,在旋翼桨叶嵌套网格的基础上,发展了一套用于预测旋翼绕流流场的计算流体力学(CFD)模拟方法。然后,为克服传统直升机旋翼二维水滴撞击特性计算方法的不足,充分考虑旋翼流场的三维效应,在嵌套网格中基于欧拉法求解旋翼三维水滴撞击流场。其中,为解决尾流等区域的密度脉冲现象所引起的稳定性和收敛性问题,提出并建立了遮蔽区扩散模型。该模型通过判断遮蔽区变量,在计算域中动态生成遮蔽区域,并随迭代步数逐渐扩散。最后,通过与NACA0012翼型及国外UH-1H桨叶的试验和计算结果的对比,验证了旋翼三维水滴撞击特性计算新方法的可靠性,并进行了温度和水滴当量直径(MVD)对旋翼三维水滴撞击特性的影响分析。结果表明:遮蔽区扩散模型的加入,使二维情况的计算时间减少了22%,并增加了三维情况的计算稳定性,显著提高了旋翼三维水滴撞击特性的计算效率;沿着旋翼桨叶展向位置增大的方向,旋翼桨叶剖面水滴撞击范围有所增大,最大水滴局部收集系数呈先增加后减少再增加的变化趋势,其变化幅度接近50%;旋翼桨叶表面的水滴撞击区域和水滴局部收集系数随水滴当量直径的增加而增加。     

水滴撞击特性的数值计算方法研究

 对水滴撞击特性的数值计算方法进行了研究。在对防冰部件流场进行计算的基础上,通过对水滴所在单元的判断和单元内流场速度的插值求解以及水滴运动轨迹计算边界的判断,采用差分法对水滴运动方程进行了数值求解,得到了水滴运动轨迹,从而确定了水滴的撞击极限、总收集系数和局部水收集系数等水滴撞击特性参数,为飞机防冰系统的设计奠定了基础;此外,以发动机进气道的水滴撞击特性的计算为算例,研究了飞行高度、飞行速度及水滴半径对水滴撞击特性的影响。     

三维复杂表面水滴撞击特性计算

为预测多段翼和发动机进气道等三维复杂形状表面在结冰气象条件下的水滴撞击特性,提出了一种基于欧拉两相流模型求解水滴收集系数的方法.空气相和水滴相认为是单相耦合的,空气流场由Euler或Navier-Stokes(N-S)方程独立求得.给出了一种水滴相撞击壁面边界条件的处理方法.为避免因局部水滴容积分数异常而导致计算发散,提出了一种自适应的数值扩散模型以增加高阶计算的稳定性.验证计算了圆柱、MS(1)-0317翼型、球面等典型二三维物面的收集系数,并计算了水滴直径呈某种分布时的情形.与公开发表的试验数据对比表明,计算方法准确有效,能很好的应用于三维复杂表面的水滴撞击特性预测.      

过冷水滴撞击三维机翼的数值模拟

基于欧拉两相流理论,提出一种数值模拟过冷水滴撞击三维机翼的方法.根据水滴拟流体模型建立三维水滴控制方程;提出一种可穿透型壁面边界模拟水滴对机翼表面的撞击;由水滴流场的求解结果得到机翼表面的水滴收集特性.对ONERA M6机翼表面的水滴收集特性进行了研究,计算结果表明所提出的数值模拟方法是有效、可行的,可以为三维积冰的数值模拟研究奠定良好的基础.      

旋转帽罩结冰相似准则的研究

为了研究航空发动机旋转帽罩结冰相似理论,在静止部件结冰相似准则的基础上,忽略次要因素的影响,建立了一种适用于旋转帽罩的结冰试验相似准则。利用建立的旋转结冰相似准则,通过数值模拟方法,对旋转帽罩缩比前后的水滴收集特性和结冰冰形进行了误差分析,并对此准则的可靠性进行评估。计算结果表明:用建立的相似准则所确定的参数对于全尺寸旋转帽罩转速分别为1kr/min和5kr/min时进行数值模拟。在全尺寸旋转帽罩及其1/2缩比模型上的对应区域得到的水滴收集特性和结冰冰形,误差约为2.7%和5.3%,表明选取缩比参数的方法是可靠的,并且提出的相似准则是有效的,对结冰风洞试验的理论研究和参数选取有一定的参考价值。     

A review of recent studies on rotating internal cooling for gas turbine blades

Gas turbines have been used extensively for aircraft and marine propulsions as well as land-based power generation because of their high thermal efficiency and large power to weight ratios. To further increase the thermal efficiency, numerous prior researches on gas turbine blade internal cooling have been intensively carried out, majorly under stationary conditions. However, the stationary studies neglect the effects of Coriolis and buoyancy forces, which should change the velocity, turbulence and temperature distribution under rotating conditions. To elucidate the rotational effects on gas turbine internal cooling, the extensive results collected from recent investigations are discussed, which include the rotation and buoyancy effects on the rib turbulated cooling, pin fin cooling, jet impingement cooling, dimple/protrusion cooling, latticework cooling as well as swirl cooling. The rotational effects on the friction factors and the most employed experimental and numerical methods are also presented. Moreover, recommendations for future research are outlined. Therefore, this review article provides extensive literature information for the design of the next-generation high-efficiency internal cooling for rotating turbine blades.     

一种补偿惯性平台基座不水平的转位算法

针对静基座平台系统调平进行了研究, 详细分析了四轴平台系统台体坐标 系、框架坐标系及基座坐标系间的空间关系,推导出了台体不同位置调平时框架角间的 关系,设计了补偿平台基座不水平的转位方法。当基座不水平角度较大时,通过该方法 可将平台直接转位至调平角度附近,不需经过粗调平便可直接进入精调平,从而可以使 平台系统的工作过程得到简化,对于提高平台的任意位置转位及自标定、自瞄准效率有 重要意义。     

旋转坐标系N-S方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明：所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

一种用于飞机结冰的海绵冰模型及应用

基于两相流理论,采用欧拉法对飞机结冰进行数值模拟,在空气流场速度分布的基础上,求解水滴流场,得到水滴收集率和撞击极限.结果表明:欧拉法比拉格朗日法在计算水滴收集效率和撞击极限方面有优势,尤其是在较大的水滴平均等效直径情况下.对翼型的冰形和圆柱表面的温度预测数据和实验数据进行对比,圆柱冰形厚度误差在15.0%之内,计算结果与实验结果吻合较好,同时表明结冰模型是有效可行的.      

瞬态法测量高转速旋转盘表面传热系数

将实际发动机涡轮盘冷却系统简化为中心进气旋转盘,用瞬态实验的方法对该结构的换热特性进行了研究.以高转速旋转换热实验台为基础,建立了转静系盘腔中瞬态换热的实验流程和数据处理方法,得到了转盘表面的努塞尔数,并研究了流量系数、旋转雷诺数和出气间隙比对努塞尔数的影响.研究结果表明:对于中心进气冷却结构,在转盘低半径处,转盘表面的表面传热主要由冷气冲击控制,随着半径的的增大冲击对换热的影响减弱.转盘表面的努塞尔数随冷气流量系数的增大而增大,随旋转雷诺数的增大而增大,随出气间隙比的增大而减小.      

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

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

体积分数法测量几种转静系中的最小封严流量

为防止涡轮盘腔转静系统高温“燃气入侵”现象发生,需要确定最小封严流量.在不同的转速和不同出气间隙下,采用体积分数法分别对6种典型封严结构进行实验研究,以确定它们在不同工况下的最小封严流量并得出其变化规律:在固定的出气间隙值下,最小封严冷气流量随转盘转速的增加而增大,且几乎呈线性变化;对于固定的转速,随着出气间隙的增加,最小封严冷气流量也随之增加.拟合出不同封严结构的级间封严准则函数关系式,并探寻出封严效果最优的封严结构.      

旋转状态下矩形微小通道流动与换热试验研究

为了研究微小通道结构在航空发动机涡轮叶片中应用的前景和可行性,以空气为冷却介质,在Re＝1000～3000、转速为0～500 r/min、Ro＝0～3.5×10-3条件下,对水力直径为1 mm的旋转微小通道组的流动和换热特性进行试验研究。结果表明：微小通道流阻系数呈现粗糙壁通道特征,通道临界Re≈2350,流阻系数以及临界Re随转速增加未见明显改变。在静止状态下,通道组综合换热系数随Re增大而增大,换热系数分布沿流动方向逐渐减小;在旋转状态下,通道组平均综合换热系数略有增大,旋转对换热特性的影响随着流动的发展而增大     

动力对全机水滴收集率的影响计算

针对飞机在飞行中遭遇过冷水滴撞击并结冰现象,建立了适合于发动机带动力情况下结冰过程水滴收集率计算的三维数值方法和计算程序.其基本思路为:采用多块技术与SIMPLE方法计算空气流场,以流场分布的计算结果为基础,求解水滴相的控制方程,进而获得物体表面的水滴收集率.空气相控制方程和水滴相控制方程均写成典型输运方程的形式,采用一致的有限体积法离散求解,方便了计算程序的编制.对某型运输机巡航构型有/无动力条件的水滴收集率进行了比较计算,获得了不同直径水滴在飞机表面的撞击特征以及水滴收集率在飞机机翼、平尾、垂尾和发动机进气道唇口上的分布规律.研究表明:(1)发动机是否带动力对机翼、平尾、垂尾的水滴收集率基本无影响;(2)飞机带动力主要影响发动机进气道唇口处的水滴收集率,带动力后唇口的收集率比无动力情况高,水滴撞击范围增大,在进行防除冰研究和设计时需引起重视.     

轴向通流旋转盘腔换热特性

针对航空发动机压气机盘腔中流动与换热现象,采用该实验的方法对等温轴向通流旋转盘腔进行研究。通过对比不同工况下局部努塞尔数和平均努塞尔数的变化情况分析旋转腔中各力对流动与换热的影响,并总结努塞尔数与各无变量参数的经验公式。结果表明:盘腔内流动与换热主要由哥氏力、惯性力与离心力控制。腔内流动区域可划分为低半径位置惯性力占主导地位的惯性对流区与高半径位置哥氏力占主导地位的旋转对流区。轴向通流流量增大使惯性力增大,转速增高使哥氏力增大。惯性力与哥氏力的增加都会增强换热,两个力对彼此对换热的影响有削弱作用,两个力的综合作用使不同工况中不同半径位置的换热变化情况不同。平均努塞尔数与局部努塞尔数随变量的变化情况基本一致。