一种涵道螺旋桨性能预测的CFD方法

为适应高性能涵道螺旋浆精细化设计需求,基于雷诺平均Navier—Stokes控制方程和非结构网格旋转／静止滑移面技术,开展了涵道与螺旋桨动／静部件干扰的非定常流及性能预测数值模拟方法研究。提出了．静止与旋转域的多块网格布局策略及转／静滑移面技术,研究了涵道壁而与螺旋桨桨梢间隙的网格布局形式,有效解决了螺旋桨旋转流及螺旋...     

静偏心对挤压油膜阻尼器减振特性影响的 数值分析

由于装配误差或重力等原因,挤压油膜阻尼器在实际使用中难以避免产生静偏心。为了研究静偏心对挤压油膜阻尼器减振特性的影响,基于广义雷诺方程推导出静偏心条件下挤压油膜阻尼器雷诺方程并建立了静偏心挤压油膜阻尼器-转子系统的动力学模型,随后运用数值积分获取静偏心下转子系统的非线性响应。研究表明:转子进动半径随着静偏心增大而减小;静偏心显著影响转子轴心轨迹,静偏心时转子轴心轨迹在大多数转速下近似为椭圆而不再是圆,静偏心下突加不平衡响应瞬态振幅远高于无静偏心条件下瞬态振幅。     

叶轮机中的焓变化

通过研究粘性非定常可压流的能量方程,讨论了单级和多级环境下在叶轮机转静叶片中和叶排间容腔中焓的变化和守恒条件。结论表明在旋转坐标系中涉及非定常流设计的压力扰动,或粘性应力对相对流动中的流体做功等因素时,转子叶片槽道中相对流动的转焓会发生变化,由于转静叶排叶排干涉造成的压力扰动和尾迹输送等因素使静叶槽道中的总焓也不守恒,而且,转静叶排之间的容腔中,由于流动的非定常性和剧烈的掺混,总焓也是变化的。     

适用于倾转旋翼机气动分析的新型嵌套网格方法

为克服围绕复杂倾转旋翼机生成高质量一体化结构网格的困难,将分块网格策略与嵌套网格方法融合,提出并建立了一套适用于倾转旋翼机气动干扰分析的新型嵌套网格生成方法。该方法中,将整体网格按照倾转旋翼机结构特性进行了分块生成,在保证插值精度的基础上降低了网格规模。为高效连接不同部件网格,首先,引入Inverse Map方法来辅助部件网格定位,采用新颖的多向投影方法保证了部件间洞边界的连续性;其次,基于独立挖洞嵌套策略完成机身与背景网格的组合,并使用了相应的辅助计算坐标进行背景边界网格标识。在贴体网格区域采用基于S A（Spalart Allmaras）湍流模型可压缩RANS(雷诺平均Navier Stokes）方程求解流场,背景网格区域采用欧拉方程,并运用了SPMD(single program, multiple data)模式的并行加速技术,建立了适合倾转旋翼机气动干扰特性分析的高效混合CFD方法。采用Robin直升机机身作为数值算例,验证了CFD方法的有效性。在此基础上着重对倾转旋翼机在悬停、过渡和前飞巡航模式下的干扰流场进行了数值分析研究,分别得出了悬停状态下倾转旋翼机典型的“喷泉效应”干扰现象、过渡状态下的非定常气动干扰特性以及巡航状态下的部件气动变化特性,计算结果表明建立的新型嵌套网格方法能够较好地表征倾转旋翼机外形特性,并能够有效地用于倾转旋翼机的气动特性分析,加速比能超过5.0。      

临近空间桨梢小翼螺旋桨布局气动增效研究

基于雷诺平均Navier—Stokes控制方程和非结构网格旋彬静止滑移面技术,开展了临近空间螺旋桨桨梢小翼构型气动布局与增效作用机理研究,采用静止与旋转域的多块网格布局策略及转／静滑移面技术,发展了高效的螺旋桨流场与性能分析方法;通过开展上翘小翼、上下翘小翼等多种桨梢小翼螺旋桨布局研究,初步建立了基础参数的性能影响规律,研究表明：桨梢小翼螺旋桨布局是临近空间工况环境提高气动效率的一种合理可行的技术途径。     

一种导管螺旋桨非定常干扰流CFD方法

基于雷诺平均Navier-Stokes控制方程和非结构网格旋转/静止滑移面技术,开展了导管与螺旋桨动/静部件干扰的非定常流及性能预测数值模拟方法研究。提出了静止与旋转域的多块网格布局策略及转/静滑移面技术,研究了导管壁面与螺旋桨桨梢间隙的网格布局形式,有效解决了螺旋桨旋转流及螺旋桨与导管干扰计算问题。典型敞开桨与导管螺旋桨算例的数值模拟结果与实验值符合一致,研究结果表明：所发展的导管螺旋桨CFD分析方法是实用、有效的。     

机翼静气动弹性数值模拟研究

通过求解基于三维非结构网格的欧拉方程计算机翼气动特性,采用结构柔度影响系数方法进行机翼静弹性变形的计算,利用三维非结构弹性网格技术进行几何外形变形下的网格修正,将以上三部分进行耦合迭代.实现了机翼的静弹性特性的数值模拟.采用开发的数值分析程序,进行了某支线飞机机翼静弹性特性的计算分析,计算结果合理,表明所发展的数值模拟方法可用于机翼静弹性气动特性的分析研究.     

双转子系统碰摩有限元接触分析模型及故障诊断

针对常用碰摩力模型的不足,应用AN SY S软件在求解转静子之间的碰摩力与法向相对位移之间的关系的基础上得到了碰摩力模型,并与整体传递系数法相结合建立转静件碰摩的动力学方程。采用整体传递系数法计算转子系统的动力特性,进而求解双转子系统的瞬态响应。利用时域波形、频谱分析和小波变换对响应进行了分析,得出并分析了相关的故障特征。      

倾转过渡状态旋翼-机翼气动干扰特性

针对倾转旋翼机倾转过渡状态旋翼-机翼之间的气动干扰开展了高分辨率数值模拟。计算采用有限体积法对积分形式的非定常Navier-Stokes方程进行离散,使用重叠网格处理旋翼桨叶与机翼之间的相对运动。背景网格采用基于八叉树结构的笛卡尔结构化网格,并使用基于无量纲Q准则的网格自适应技术来获得高分辨率尾迹流场。首先,计算了孤立旋翼和旋翼-机翼组合体悬停状态,根据计算结果证明高分辨率求解器对于旋翼气动力预测的准确性和旋翼尾迹结构演化的高分辨率特性,对比了旋翼-机翼组合体悬停状态改进延迟分离涡模拟(IDDES)和RANS计算结果,IDDES所得流场更为精细且气动力结果与RANS结果也有所差异。随后,对倾转过渡状态不同倾转角下旋翼-机翼组合体飞行状态进行了仿真计算,得到并分析了旋翼与机翼之间气动干扰流场。结果表明,在倾转过渡状态中期旋翼产生诱导滑流对机翼有着一定的升力增益效果,而倾转初期与末期则无明显升力增益,初期甚至为负增益;旋翼与机翼气动变化情况证明倾转过渡状态旋翼-机翼气动干扰对倾转旋翼气动性能有着重要影响。     

倾转旋翼典型飞行状态气动特性的CFD分析

针对倾转旋翼特有的悬停、巡航、过渡3种飞行模式,基于CFD方法分别对倾转旋翼的上述3种典型飞行状态进行了数值模拟.首先,基于运动嵌套网格思想建立了适合倾转旋翼运动模拟的多层运动嵌套网格系统.然后,基于RANS(Reynolds averaged Navier-Stokes)方程建立了适用于非定常流场分析的CFD模拟方法.为了提高整个流场的计算效率,时间推进采用高效的隐式LU-SGS(lower upper symmetric Gauss Seidel)方法,并运用SPMD(single program multiple data)并行计算模式.最后,采用所建立的网格生成方法和流场求解器,进行了不同飞行状态下的倾转旋翼气动特性数值分析,验证了建立CFD方法的有效性并得出了多种飞行状态下的气动特性.结果表明:悬停状态下倾转旋翼相对于直升机旋翼有更加明显的下洗流影响区域.巡航状态两旋翼中心间距为2.5R时,倾转旋翼之间存在明显的气动干扰.过渡状态下,倾转旋翼提供的拉力在倾转角为45°前下降较快,且在倾转角较小时,倾转旋翼桨叶表面拉力分布仍能表现出直升机旋翼的特性.      

双级风扇三维粘性非定常流动的数值解

针对一台双级跨音速轴流风扇动静叶片排相互干扰的三维非定常流动,进行了N-S方程的数值求解。首先计算风扇定常流场,以定常流收敛结果为起始场计算非定常流动。计算结果表明,非定常流动对各叶片排进口气流角有较大的影响,其中静叶的变化大于动叶,与定常流结果相比,非定常的瞬态攻角增大的量大于减少的量,最多时增大近10°左右。静子叶片受到的非定常气动力变化幅度也大于转子叶片,叶片排所处的轴向位置不同,非定常流动的影响会有较大的差别。      

轴流压气机转静二维非定常薄层N-S数值模拟

采用非定常薄层Navier-Stokes紊流运动方程来描述一跨声轴流压气机中径处的二维转-静非定常流场，速度脉动的相关分析表明转子尾迹对静子的热流流场几乎没有什么影响，而粘性的影响取决于叶片排之间的轴向距离。     

多级压气机非定常流动的数值模拟

进行了多级压气机动静叶片排干扰的三维非定常无粘流动数值模拟研究，采用高分辨率格式求解非定常Euler方程组。对某以级压气机的非定常无粘流场的计算结果表明，即使在不考虑粘性影响的情况下，压气机动静叶排干扰的非定常效应仍然是比较显著的，各叶片排所受到的非定常影响是不相同的，计算得到的叶片表面压力周期性波动的结果可为叶片强度计算提供较为准确的数据。     

风扇/压气机转子、静子干涉噪声的预测方法

 本文采用管道声学模型计算了由转子尾迹与静子相互作用产生的干涉声音的大小。分析了转子、静子相对距离的变化对声音大小的影响,并与试验结果进行了对比。     

非均匀栅距对压气机叶片非定常气动力的影响

研究了轴流压气机叶片采用非均匀栅距对叶片非定常力影响的问题,利用高分辨率格式求解准三维雷诺平均N-S方程。非定常流场计算采用了双时间步法与LU-SGS隐式解法耦合的方法。研究了带进口导叶的单级三排叶栅非定常流动,其中动叶栅距不变,而导叶和静叶采用非均匀栅距。对不同结构的计算分析表明,在导叶和静叶采用一定的非均匀栅距组合型式后,在保证压气机级的性能变化不大的同时,动叶表面所受非定常气动力的状况得到改善,这将有助于提高动叶片的疲劳寿命。     

跨音速压气机动静叶排相干的三维非定常流动数值分析

采用时空二阶精度ＮＮＤ显式差分格式求解三维Ｅｕｌｅｒ方程，对跨音速压气机级中动静叶排相干形成的非定常流动进行了数值分析。除动静叶排流动的控制方程体系、数值方法外，着重讨论了动静叶排计算域、网格以及边界条件的处理。对某压气机第１级动静叶排相干的非定常流场的计算，给出了叶片表面的马赫数分布、压力波动幅值，以及叶片槽道中瞬态密度等值线等，并分析了结果的合理性。     

Unsteady Flow Variability Driven by Rotor-stator Interaction at Rotor Exit

Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.     

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

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

静叶时序对压气机叶片附面层流动影响的数值研究

采用数值方法对两级低速压气机中径处的非定常流场进行模拟,针对压气机第2排静叶两个典型周向位置对动静叶干扰下的叶片附面层流动进行研究.建立尾迹与附面层干扰分析模型,结合叶片壁面摩擦力和近壁面附面层湍动能,详细分析了尾迹和势流干扰下静叶时序改变对叶片附面层流动产生的影响.对第2排静叶附面层的研究结果表明:静叶时序改变了尾迹在其叶排中的输运特征,能够降低壁面摩擦力和近壁面湍动能及其非定常最大波动幅值,影响吸力面附面层内动叶尾迹后沉寂区的宽度.在非定常条件下,尾迹能够诱导静叶层流附面层在尾迹干扰的局部范围内转捩发展为湍流状态,同时高湍流度尾迹的干扰具有抑制逆压梯度下附面层分离的作用,并能够延长层流区的范围.     

时间倾斜法在叶轮机械非定常模拟中的应用

应用时间倾斜法处理非定常模拟中前后两排叶片栅距不相等的问题,通过对流动控制方程的时空转换,发展了相邻两叶排任意转速下非等栅距的非定常时间精确算法,使该方法拓展为适用于超常规动/动结构叶轮机械内的非定常流动计算.在此基础上编写了完整的时间倾斜法计算程序与约化算法计算程序.对某进口导叶/转子之间静/动干涉的非定常流场进行的计算表明:与约化算法相比,时间倾斜法采用原始单通道模型,真实地反映了流场的非定常性;同时,计算量仅为全通道模拟的2/41,大大节约了计算成本.