基于两相流模型的挤压油膜阻尼器空化流场特性数值模拟

基于气液两相流理论中的Mixture模型和Zwart-Gerber-Belamri空化模型,采用动网格技术建立了两端开口中心槽供油型挤压油膜阻尼器三维非定常空化流场求解模型。数值模拟表明:随着阻尼器内环进动,考虑两相流动的油膜低压区在中心槽的两侧产生两道对称的条状负压带,条状带内具有较高的气相体积分数;阻尼器低压区的压力和气相体积分数对进油孔位置十分敏感,阻尼器油膜力和流场气穴比的变化频率与进油孔数密切相关。内环同心进动半径、进动频率对阻尼器空化流场影响的数值计算表明,进动半径和频率的增大均会使得流场内的空化现象加剧,同时气穴比的相位滞后现象愈加显著。      

对转压气机外伸激波对叶顶泄漏流影响

摘要:为了揭示对转压气机下游转子外伸激波对上游转子泄漏流的影响规律,针对上游转子叶顶间隙分别为0.2、0.5、0.8 mm的对转压气机开展了非定常数值模拟研究。研究发现:受下游转子外伸激波掠扫影响,上游转子尾缘附近压力面会形成弱压缩波,且随上游转子泄漏流增强而逐渐减弱;而该外伸激波在上游转子尾缘附近吸力面,会形成与型线切向相垂直的较强压缩波,且其位置基本不受叶顶间隙大小影响;外伸激波使上游转子尾缘附近吸、压力面压差增大,叶顶泄漏流增强,进而导致其损失增大;随着叶顶间隙增大,上游转子叶尖区弦长前半段压力波动的频率,由通道激波转为叶顶泄漏流主导,且呈现减小的趋势,而弦长后半段压力波动的频率主要由外伸激波主导,且基本不变。      

离心压缩机密封及空腔流动的一维/三维耦合建模研究

为了在保证数值预测精度的同时降低离心压缩机的建模复杂度,以一维密封泄漏计算替代完整的密封及空腔建模,开展了离心压缩机盖盘密封及空腔流动的源项建模方法研究。结果表明:基于源项方法的密封及空腔简化建模仅需在主流道建模基础上施加泄漏质量流量、总温和摩擦损失源项边界条件,即可较高精度地还原完整的盖盘密封及空腔建模的性能和关键流动细节。与完整建模相比,在全工况下源项建模预测的等熵效率的绝对偏差小于0.65%,总压比和等熵效率的相对性能偏差小于1%。该方法可在满足下游部件匹配设计精度需求的同时,减少50%以上的网格数量,并显著提高计算的经济性和收敛性。泄漏流量的计算精确度是影响源项建模法性能预测准确性的最关键因素,耦合建模可以采用一维计算中推荐的密封泄漏系数值。      

回流燃烧室流动特性试验

为了揭示有/无燃烧状态下燃烧室热态和冷态流场的特征和流动特性,针对某型回流燃烧室单头部试验件,使用粒子图像测速仪（PIV）,测量燃烧室燃烧状态下不同截面处的热态流场,以及没有燃烧状态下不同截面处的冷态流场,探讨不同总压损失系数对回流燃烧室热态/冷态流场特征及流动特性的影响。研究表明:随着总压损失系数的增大,冷态条件下各截面流场结构基本保持不变,如射流孔穿透深度、射流角度、回流区位置及大小、流线等基本保持一致,但是各位置点速度大小逐渐增大。热态条件下各截面流场随着总压损失系数增大,流场结构也基本保持不变;相同总压损失系数时,热态流场与冷态流场存在差异,燃油喷射与气流的相对运动将会对燃烧室头部的流场结构造成影响,速度较冷态流动时略微增大。      

偏心撞击对撞击式喷嘴雾化特性的影响

为研究偏心撞击对撞击式喷嘴雾化特性的影响,建立了求解自燃推进剂冷态射流撞击雾化过程的数值模拟方案,计算了不同偏心度条件下的射流撞击雾化过程。采用树形自适应加密算法直接求解不可压Navier-Stokes方程组,由分段线性的流体体积(VOF)方法对流体界面进行捕捉。结果表明偏心撞击会导致雾场发生偏转,当无量纲偏心度E为1/8时,雾场偏转角度约为9.2°,应控制加工偏差小于该值。随着偏心度的增大,液膜的偏转角度增大,理论推导得到的液膜偏转角度要小于数值计算得到的液膜偏转角度。正心撞击时燃料与氧化剂流强峰值接近,雾场的流强分布呈单峰分布。当发生偏心撞击时,由于燃料与氧化剂部分射流未参与撞击导致流强峰值出现交错,雾场的流强分布呈双峰分布,混合比的空间分布发生较大改变。正心撞击时撞击点下游液滴的速度分布近似呈轴对称分布,而偏心撞击之后的速度分布则呈中心对称分布。偏心撞击导致的射流动量损失使得雾化性能变差,当无量纲偏心度E为1/8时,一甲基肼(MMH)的Sauter平均直径增大约4.8%,四氧化二氮(NTO)的Sauter平均直径增大约5.8%。      

大型水陆两栖飞机侧风起降的滑流影响分析

大型水陆两栖飞机AG600的动力装置为安装在机翼上的4台同向旋转涡轮螺旋桨发动机,针对1：15缩比模型带动力风洞试验显示的螺旋桨滑流对侧风起降状态的偏航力矩不稳定影响,对全机带动力风洞试验模型进行了大规模并行非定常数值计算,再现了风洞试验现象,通过流动机理分析明确其产生原因主要是左侧滑时右外翼分离和垂尾背鳍涡破裂,这些原因和数值模拟的准确性也为后期的风洞试验所证实。考虑到模型风洞试验中尺度限制造成的低雷诺数和高螺旋桨转速,为保证飞行安全,继续采用该非定常方法对全尺寸飞机真实侧风起降状态进行了详细数值分析和偏航稳定性评估。研究结果显示,在飞行雷诺数和螺旋桨转速下,相同侧风范围内风洞试验显示的流动不稳定因素基本消失,偏航稳定性允许的侧风范围明显增加。本研究实现了四发螺旋桨飞机起降状态横向气动特性的滑流影响非定常数值分析,建立了基于计算流体力学的风洞与飞行雷诺数效应的相互关系,进行了偏航稳定性的虚拟试飞评估,研究成果也为AG600飞机的首飞和飞行试验所验证。     

涡桨飞机螺旋桨滑流气动干扰效应及流动机理

螺旋桨滑流对飞机各气动部件的干扰是涡桨飞机气动设计中面临的难点之一。以某双发涡桨支线客机为对象,采用数值模拟手段分析了滑流对机翼、平尾、垂尾的影响及其流动机理。计算采用基于动态面搭接网格技术的非定常方法,气动力与表面压力分布的计算结果均与实验值吻合良好。研究结果表明：在滑流影响下,全机升力和阻力有所增加,升阻比和纵向静稳定度有所降低,并在无侧滑状态下产生了滚转力矩和偏航力矩;在不同展向位置,滑流对机翼表面流动分离的影响存在显著的差异。在螺旋桨下行运动一侧,滑流的旋转减小了当地迎角,同时桨后气流速度较高,翼面流动分离被有效抑制。在螺旋桨上行运动一侧,滑流的旋转增大了当地迎角,而且桨后气流速度较低,因而翼面流动分离并未得到明显改善;在中小迎角下,滑流对背景飞机平尾当地的动压没有产生影响,但增加了下洗角变化率,进而导致平尾效率降低;垂尾的侧力与偏航力矩是由滑流的侧洗引起的,而滑流的侧洗又与左右两侧机翼升力分布的不对称性有关。     

翼型低雷诺数层流分离现象随雷诺数的演化特征

层流分离现象是翼型低雷诺数条件下出现的典型流场特征。层流分离流动中包含流动分离、转捩、再附等非定常流动结构,层流分离流动的形成与演化会对翼型气动特性产生恶化作用。采用大涡模拟（LES）方法对低雷诺数范围内不同雷诺数下的翼型层流分离流动开展精细数值模拟,研究了雷诺数对翼型气动特性的影响规律及作用机理。LES方法采用隐式亚格子模型,基于结构化拼接网格,对流项离散和时间推进方法分别采用AUSM⁺格式以及双时间步方法。验证算例计算结果表明数值模拟方法的正确性及可靠性,雷诺数对翼型气动特性具有显著影响。随雷诺数降低,时均分离泡外形增大、位置后移,平均阻力系数增大,特别是在较低雷诺数下,翼型升阻力系数随时间出现振荡现象。进一步研究表明,造成不同时均分离泡形态和气动特性的原因在于翼型上表面分离剪切层的失稳与转捩特征。随雷诺数降低,流动黏性增大,导致分离剪切层速度梯度减小,流动发生转捩及再附位置后移,直至翼型表面不再发生转捩和再附。     

旋转整流罩积冰过程中水膜流动和脱离的数值研究

为了更准确地模拟旋转整流罩积冰过程,在现有三维积冰与冰层表面薄水膜流动耦合模型基础上,基于功平衡分析的方法引入了旋转部件表面水膜脱离模型,并发展了相应的计算方法,给出了水膜脱离的判定依据:当气流曳力做的功和由于离心力使水膜增加的潜能之和大于黏附功时整流罩表面的水膜会发生脱离。对旋转整流罩积冰进行数值模拟,计算结果与实验结果吻合得较好,验证了该模型的合理性和计算方法的可行性。之后分析了转速和来流速度对整流罩表面水膜脱离和积冰的影响,结果表明:转速和来流速度越大,水膜发生脱离的比例越大。在研究范围内,转速为3000r/min和6000r/min时,因水膜脱离导致积冰总量分别减少13.4%和15.8%;来流速度为40、50m/s和60m/s时,因水膜脱离导致积冰总量分别减少为12.2%、13.4%、14.2%。      

一种可变流量系数的通气短舱匹配方法

要获得与动力短舱流量系数（MFR）一致的通气短舱,可在飞机气动设计和风洞试验过程中使用通气短舱（TFN）代替动力短舱（PN）模拟发动机短舱效应和短舱进气,以简化设计流程和降低试验成本。在研究流量系数变化对短舱和进气道气动特性影响时,需要匹配多种不同流量系数的通气短舱方案,本文在保持短舱外罩和进气道外形不变的条件下,在短舱内增加锥形堵块,通过调整堵锥的位置可快速获得不同流量系数的通气短舱,相较传统设计通气内涵匹配流量系数的方式效率更高。应用所提出的思路设计的独立通气短舱风洞试验方案得到了验证,实测流量系数与设计值吻合良好,证明了方法的可行性。