基于混合CFD方法的旋翼桨-涡干扰噪声计算

建立了一个基于Euler方程和自由尾迹分析的新的、高效的混合计算方法,用于旋翼桨-涡干扰(BVI)气动和噪声特性的研究.在该方法中,仅针对单片桨叶进行计算,通过求解Euler方程来捕捉流场中旋翼近尾迹的影响以及桨叶表面的非定常载荷,流场空间中旋翼远尾迹的位置由自由尾迹模型给出;提出了广义网格速度方法以将旋翼远尾迹的影响...     

共轴式双旋翼悬停流场和气动力的CFD计算

考虑到共轴式双旋翼流场特征高度复杂、桨叶承受非定常气动载荷的特点,为更好地预测共轴式双旋翼的气动特性,把非结构嵌套网格方法和网格的自适应技术相结合,发展了一套适合于共轴式双旋翼流场数值计算的求解器.在该求解器中,采用非结构嵌套网格方案来模拟桨叶之间存在的相对运动,自适应网格技术用来捕捉尾迹对流场和气动特性的影响,求解惯性坐标系下的非定常N-S主控方程来模拟流场的非定常特性.应用该求解器,首先计算了有试验结果可供对比的一副试验旋翼的诱导速度场分布,在此基础上,计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性,并与单旋翼的计算结果进行了对比和分析,得到了一些有意义的结论.     

共轴刚性旋翼悬停状态地面效应气动特性

为了研究地面效应下共轴刚性旋翼的气动特性,建立了一套基于非定常雷诺平均Navier-Stokes方程的气动干扰数值方法,采用运动嵌套网格模拟双旋翼的反转运动。地面采用无滑移边界条件,并对旋翼和地面附近的网格进行加密,以更好地捕捉旋翼的流场细节和尾迹特征。计算结果与Lynx尾桨试验结果进行对比,验证了所建立方法的有效性。对地面效应下共轴刚性旋翼的气动性能和流场进行分析,结果发现：相对于单独的上下旋翼而言,共轴旋翼地面效应下的拉力增益更大,这是由于上下旋翼桨叶表面的压强干扰受地面高压的影响而减弱;地面的干扰主要影响双旋翼尾迹的径向位置,对其轴向位置影响不大,上下旋翼尾迹在地面附近相互融合、分裂,形成复杂的桨尖涡尾迹;双旋翼在地效下的尾迹径向扩张半径比单旋翼大,这是由于双旋翼的径向射流速度更大;随着旋翼距地面高度的增加,双旋翼间的气动干扰强度逐渐恢复,因此下旋翼拉力增益的下降速度比上旋翼更大;共轴旋翼桨尖涡相对卷起高度和扩张半径均随离地高度增加而减小。     

基于Navier-Stokes方程的旋翼前飞状态气动特性数值模拟研究

建立了一个适用于旋翼前飞状态非定常流场和气动载荷计算的数值方法。在该方法中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上将三阶逆风格式(MUSCL)与通量差分分裂方法相结合。为较好地模拟流动分离,采用了一方程的Spalart-Allmaras湍流模型,时间方向则采用双时间法推进求解。为计入桨叶的旋转、周期性挥舞和变距运动,采用了运动嵌套网格方法。此外,为了更真实地反映旋翼桨叶的实际运动,发展了一个旋翼配平分析模型及求解方法。应用所建立的方法,首先计算了具有先进气动外形的四片桨叶的UH-60A直升机旋翼悬停状态的流场及气动性能,验证了该方法对悬停状态旋翼气动性能数值模拟的有效性。在此基础上,着重对Caradonna模型旋翼无升力前飞状态和AH-1G直升机旋翼桨-涡干扰前飞状态的气动特性(表面压强分布、扭矩、气动载荷)进行了计算,数值结果和试验值吻合良好,表明本文建立的CFD方法对旋翼前飞状态流场的模拟和气动载荷的计算是有效的。     

基于IBC方法的旋翼BVI噪声主动控制机理

为揭示单片桨叶控制(IBC)主动控制技术抑制旋翼桨-涡干扰(BVI)噪声的降噪机理,建立了一套基于CFD/CSD/FW-H_pds方程的综合噪声分析方法。旋翼桨-涡干扰噪声与旋翼桨叶载荷特性、气动变形以及旋翼桨尖涡结构等密切相关,为有效模拟旋翼桨叶的载荷特性及桨尖涡结构,将Navier-Stokes方程作为前飞流场的主控方程,空间离散上采用三阶MUSCL插值格式与通量差分裂Roe格式相结合;时间方向上采用双时间法,使用隐式LU-SGS格式在伪时间方向上进行推进;湍流模型采用对分离流动具有较好捕捉能力的Spalart-Allmaras模型。为提高旋翼桨叶弹性变形运动的模拟精度,建立了基于Hamilton变分原理的CSD模型,并与高精度的CFD求解器结合,发展了适合旋翼桨叶变形及载荷特性模拟的流固耦合分析方法。在CFD/CSD耦合方法分析流场基础上,使用可穿透空间积分面的FW-H_pds方法对旋翼气动噪声特性进行计算。首先,对流场及噪声数值方法进行验证;然后,着重针对UH-60A旋翼的斜下降飞行状态,分别对有/无IBC噪声主动控制条件下的旋翼BVI气动噪声特性进行了模拟,相位角、幅值和频率等不同控制参数的影响对比分析结果表明:IBC主动控制减小了前行侧桨叶表面尤其是桨叶尖部的负压峰值,降低了桨-涡干扰发生位置附近的桨叶气动载荷;同时主动控制后的桨尖涡集中程度变弱,并且增加了桨叶与桨尖涡之间的相遇距离,从而显著降低了桨-涡干扰噪声;选取合理的相位角、幅值和频率等主动控制参数组合,BVI噪声降低可达5~7dB。     

悬停状态共轴刚性双旋翼非定常流动干扰机理

基于运动嵌套网格方法,建立了一套适合于悬停状态下共轴刚性双旋翼非定常干扰流场分析的计算流体力学(CFD)方法。首先,基于高效的运动嵌套网格技术,采用积分形式的可压雷诺平均Navier-Stokes(RANS)方程作为双旋翼非定常流场求解控制方程,湍流模型选用Baldwin-Lomax模型,时间推进采用双时间方法。在CFD方法的验证基础之上,对干扰过程中的桨尖涡涡核位置及强度演变规律进行了细致分析,揭示了共轴双旋翼非定常干扰流场中上、下旋翼桨尖涡与双旋翼桨叶之间的贴近干扰、碰撞现象,以及上、下旋翼桨尖涡之间的相互干扰机理。然后,进一步研究了不同总距角下的共轴旋翼系统中上、下旋翼的非定常气动特性以及影响规律。计算结果表明:上旋翼桨叶的桨尖涡会直接与下旋翼桨叶发生碰撞,导致下旋翼桨叶拉力损失;上旋翼桨叶的桨尖涡和下旋翼桨叶的桨尖涡相互干扰,改变了桨尖涡的强度和轨迹;上、下旋翼桨叶相互靠近时,上、下旋翼桨叶的拉力均会上升,之后相互远离时上、下旋翼桨叶拉力均会先下降再上升。     

基于滑移网格的倾转旋翼机全机干扰流场研究

为了对倾转旋翼机全机干扰规律及机理进行研究,对某倾转旋翼机的孤立机身流场特性、悬停及巡航状态下全机干扰流场特性进行模拟分析.首先,分别生成围绕机身和旋翼桨叶的贴体网格,然后通过网格搭接得到全机干扰流场计算网格;基于滑移网格方法,建立一套适用于倾转旋翼机全机气动干扰特性分析的数值方法,并以有试验结果对比的ROBIN机身进...     

基于气弹耦合的旋翼桨涡干扰噪声预测方法

建立了基于气动/弹性耦合的旋翼桨涡干扰（BVI）气动和噪声分析方法。气动模型包括修正Beddoes尾迹模型和CFD模型,噪声计算采用基于声学类比法推导出的FW H（Ffowcs Williams Hawkings）方程,弹性桨叶动力学建模采用有限元方法。应用所建立的方法,对刚性的OLS（operational load survey）旋翼桨涡干扰状态的气动和噪声特性进行了计算,对比了两种气动模型在研究桨涡干扰问题的有效性;以弹性的HART Ⅱ旋翼为研究对象,分析了桨叶弹性、时间步长对桨涡干扰气动载荷和噪声的影响。结果表明:进行桨涡干扰计算时所采用的时间步长不宜超过2°。CFD方法由于固有的数值耗散,计算出的OLS旋翼噪声声压峰值仅为试验值的60%,而修正Beddoes尾迹模型能够避免数值耗散,且具有高效率的优势。考虑桨叶气动弹性能够提高旋翼桨涡干扰噪声的预测精度。      

直升机旋翼桨-涡干扰状态非定常气弹载荷高精度预估

为准确计算直升机旋翼在复杂的桨-涡干扰(BVI)状态下的气弹载荷,在刚性旋翼计算流体力学方法中引入桨叶弹性变形的影响,建立了一套适合于弹性旋翼BVI状态气动特性分析的计算流体力学/计算结构力学(CFD/CSD)耦合方法.CFD模块对Reynolds averaged Navier-Stokes(RANS)/Euler方程进行求解,并采用双时间法推进和Baldwin-Lomax(B-L)湍流模型.CSD模块采用中等变形梁假设的有限元模型,通过Newmark-Beta方法求解桨叶运动方程.通过代数变换方法进行桨叶网格变形,并建立一个适于流场/结构信息交换的CFD/CSD耦合方法.在分别验证CFD和CSD模块的有效性的基础上,开展UH-60A直升机旋翼的BVI 状态载荷分析,并与飞行测试数据进行了对比.计算结果表明:相比于旋翼综合分析中的升力线理论和刚性旋翼CFD方法,耦合的CFD/CSD方法可以更准确地预测BVI状态气弹载荷,并有效地模拟桨叶前行侧方位角和后行侧方位角附近的BVI现象,对BVI导致的升力波动幅值和相位的计算结果均与试验值吻合良好.      

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。     

空调器室外机轴流风机内部流动的PIV研究

采用粒子图像测速仪PIV,对具有半管道式结构特点的空调器室外机轴流风机内部流场进行了实验研究,并结合实验结果分析了叶片顶部的叶尖涡和叶片出口尾缘涡的流动特性。实验结果显示在轴流风机流道内部叶顶区域存在与叶轮旋转方向相反的叶尖涡结构。叶尖涡产生于叶片前缘叶顶近吸力面侧,在流道内部与主流发生干涉后朝向周向和出口传播并逐渐耗散。叶尖涡涡心轨迹与叶顶弦长方向的夹角为10°,在叶高方向上叶尖涡的径向位置并不固定。与普通管道内部流动不同,叶片顶部与导风罩间的间隙中未捕捉到明显的叶顶泄漏涡现象。叶片出口近尾缘处30％以上叶高明显捕捉到尾缘涡结构,叶片压力面和吸力面侧的径向速度存在明显的方向变化,切向速度在尾迹区增加。     

跨声速压气机转子旋转失速现象的非定常模拟

为研究跨声速压气机转子失速机理,全周非定常数值模拟了某跨声速压气机单转子的失稳过程。结果表明:该转子由叶尖Spike扰动诱发旋转失速。在小流量稳定工作状态,压气机转子叶尖区域存在"旋转不稳定"(Rotating Instability,RI)流动现象。压气机节流过程中,转子进出口的流量降低,叶尖区流场非定常波动幅值增大。近失速状态时,RI扰动团的典型流场结构"径向涡"在叶尖区域形成堵塞,导致相邻叶片前缘间歇性地出现溢流现象。随着压气机进一步节流,转子叶尖的负荷达到极值,叶片通道尾缘逆压力梯度过大,出现倒流。尾缘倒流的出现又进一步增加通道内的堵塞,最终形成Spike扰动。失速先兆对应的流场结构是沿叶片前缘额线向相邻叶片压力面周向运动的"径向涡"结构。     

斜下降旋翼桨涡干扰声源定位及表面压力

旋翼桨-涡干扰（BVI）是直升机在进场和离场等近地飞行时后行桨叶切割前行桨叶脱落桨尖涡产生的气动扰动，该扰动不仅对桨叶表面压力载荷产生激励作用，同时也会引起旋翼噪声出现激增，旋翼噪声激增的主要成分为桨-涡干扰噪声。本文首先对斜下降桨-涡干扰状态桨叶表面载荷进行数值计算；然后分别阐述了基于改进整周期同步平均旋翼噪声去噪方法、基于多层小波包分解的桨-涡干扰声源识别和分离方法以及基于bartlett时延计算和球面插值的声源定位方法，设计并开展了风洞斜下降状态桨-涡干扰桨叶表面压力和声源定位试验，给出了开发的声源定位软件界面、声源定位图像畸变校准方法及声阵列现场校准方法；最后对比分析了不同试验状态的桨-涡干扰噪声声源特性以及和桨叶表面压力之间的关联性，并给出了声源定位及表面压力试验数据分析结果。结果表明典型斜下降状态后行侧桨-涡干扰主要出现在方位角310°~330°、径向位置1.6~1.8 m桨盘平面区域。     

动叶尾迹对静子非定常气动性能影响的研究

本文利用南京航空航天大学单级低速轴流压气机试验器,对动叶尾迹对静子非定常气动性能的影响开展了系统的实验研究。通过静叶表面动态压力分布的测量,获得了在不同的动—静叶排轴向间距、不同的工作转速和不同的工作流量下,动叶尾迹对静子非定常气动性能影响的大量基础实验数据     

低速轴流压气机叶尖非定常流动试验

利用动态压力传感器对一低速轴流压气机转子的叶顶间隙流场进行详细的试验测量,通过对信号特征的分析,对压气机节流过程中叶顶间隙的非定常流动发展演变规律进行了研究。结果表明：压气机完全失速时,叶尖存在一以46.5%转子转速周向传播的失速团;节流过程中,叶尖前缘处的动态压力信号中存在非定常波动的特征频率带,其变化规律与叶顶流场压力非定常波动的能量迁移有关;随着压气机流量减小,叶顶泄漏流影响区域向前缘移动,失速团在叶顶前缘附近产生,并向尾缘方向扩展,最终覆盖叶片全部弦长;近失速工况时,叶顶间隙相邻通道内泄漏流相互作用,造成通道中的低压区“一前一后”交替分布从而形成一个空间上周期约2个叶片通道的扰动波。     

耦合进排气壳的增压器涡轮流动干扰及激振力弱化研究

为探究船用增压器涡轮在耦合进排气壳条件下的流动特性,采用SST湍流模型和FFT方法对某增压器轴流涡轮与进排气壳耦合的流场进行非定常数值模拟,重点分析排气壳非对称流场对涡轮动叶片压力的干扰特性和扰动来源,在此基础上提出基于导叶非对称布局的激振力弱化改型方案并进行了验证。结果表明：进气壳影响静叶前缘静压沿周向的分布,排气壳流场导致的非对称背压会给动叶表面造成低频的压力波动,约为3.6%的转子通过频率,在尾缘吸力面处该低频波动幅值大于动静叶排干涉导致的高频波动幅值;排气壳内部复杂的分离流动和涡系结构是产生上述流动干扰的主要原因。4种非对称布局方案都能在几乎不影响涡轮性能的前提下分散并弱化高频波动幅值;其中Case2尾缘高频幅值减小了98.6%,低频幅值减小了52.6%,而且涡轮效率还有略微提高,为最优布局。     

不同转速下跨声速轴流压气机内部流动失稳的机理

以跨声速轴流压气机转子NASA Rotor 67为研究对象,采用数值模拟方法,开展100%、80%及60%转速下跨声速轴流压气机内部流动失稳触发机制的机理研究。数值结果与实验数据的对比分析表明:在3个转速下,数值总性能曲线的变化趋势与实验数据符合一致。通过压气机内部流场的详细分析,得出其基本流动机理。在3个转速下,随着压气机节流,叶顶泄漏涡(TLV)的起始位置逐渐向叶片前缘移动,叶顶泄漏涡也逐渐向相邻叶片压力面偏转,相比近峰值效率点,近失速点时在100%、80%以及60%转速下叶顶泄漏涡的偏转角度分别为3°、6°和9°。在100%和80%转速下,叶顶泄漏涡与激波相互作用所导致的堵塞是触发压气机内部流动失稳的机制,并且在80%转速下,叶顶泄漏涡发生破碎;而在60%转速下,泄漏涡在相邻叶片出现的叶顶前缘溢流(LESF)是触发压气机内部流动失稳的主要机制,叶片吸力面尾缘出现的小尺度附面层气流分离(BLFS)不是主要机制。      

上游叶片尾迹对转子叶片非定常表面压力频谱特性影响的研究

为了研究上游叶片尾迹对轴流压气机转子叶片非定常表面压力的影响,采用在叶片表面埋设微型压力传感器的方法,测量了一个单级低速轴流压气机转子叶片的非定常压力分布。对转子叶片非定常压力分布的分析显示:转子叶片压力面非定常表面压力主要受上游叶片尾迹影响,其主导频率为上游叶片尾迹频率及其倍频,压力波动幅度随上游叶片尾迹的衰减而沿流向减弱。转子叶片吸力面非定常压力受叶片附面层和上游叶片尾迹共同影响,其主导频率为上游叶片尾迹频率及其倍频,但是在吸力面前部气流加速区压力波动幅度沿流向增大,而在吸力面后部气流减速区压力波动幅度沿流向减小。在存在尾缘分离的情况下,在分离区附近产生较大压力波动。尾缘气流分离产生的压力波动频谱中含有上游叶片尾迹频率及其倍频的分量,但其频谱可能比较复杂,并非单一地受上游叶片尾迹影响。      

压气机动静叶干扰的非定常流动特性

采用大涡模拟(large eddy simulation,LES)方法计算动叶尾迹对静叶干扰的流场信息。利用涡量分布揭示动叶尾迹在静叶通道内的演化过程,利用压力梯度识别激波结构及波振源,运用动力学模态分解(dynamic mode decomposition,DMD)方法对静叶通道内流场的时空结构进行模态分解。结果表明:流场中存在3处波振源,分别位于动叶尾缘、静叶前缘和静叶尾缘处;发现静叶通道内流场的频谱具有多峰值现象,模态分解的第1阶流动代表动叶尾迹在通道内随时间迁移,对应频率为动叶通过频率(blade passing frequency,BPF)是通道内旋涡非定常波动的主导频率;第2阶流动是动叶通过频率的2倍频流动,旋涡的空间尺度为1阶模态的1/2,为更小尺度的扰动。