考虑转捩影响的翼型动态失速数值模拟

用数值求解雷诺平均NS方程的方法来对考虑转捩影响的中低雷诺数下振荡翼型动态失速进行数值模拟,计算中采用了k-ω SST两方程湍流模式,并加入Chen-Thyson转捩模型来模拟流动中的转捩效应.采用该方法分别对雷诺数Re=1.35×105和Re=7.7×104情况下NACA0012翼型的动态失速进行了数值模拟.计算结果显示:计算出的翼型动态失速气动力系数迟滞曲线与实测结果符合较好;对于Re=7.7×104的工况,实测的升力系数迟滞曲线中出现了高频振荡,计算结果有效地捕捉到了这一流动现象,并通过分析瞬时流线的计算结果揭示出尾缘涡的涡脱落是引起该高频振荡的主要原因.此外,通过数值计算分析了中低雷诺数下减缩频率对升力系数迟滞曲线的影响,从结果中看到,随着减缩频率的增加,翼型的失速攻角值会增加,升力系数峰值会增加;当减缩频率减小时,升力系数曲线中的高频振荡频率会增加.通过进一步计算分析知,升力曲线中高频振荡的产生不仅取决于减缩频率,还取决于流动的雷诺数,只有在中低雷诺数、较小的减缩频率下翼型动态失速的升力迟滞曲线中才有可能出现高频振荡.     

基于eN转捩预测方法的增升装置失速特性数值模拟研究

基于三维可压缩Navier-Stokes方程耦合e~N转捩预测方法,对30P30N二维多段翼型和DLR-F11三维增升装置进行了数值模拟研究。计算结果表明,e~N方法计算得到的转捩位置与试验值符合较好,增升装置翼面特别是主翼和襟翼下翼面存在大范围层流流动,为提高增升装置绕流计算准度,计算中考虑转捩的影响是非常有必要的。相比全湍计算,e~N转捩计算升力系数与试验值符合更好,提高了增升装置大迎角失速特性计算准度。     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

翼型跨声速动态失速过程的数值研究

本文发展了一种模拟翼型跨声速动态失速过程的数值方法。     

NS方程中耦合转捩判断的风力机翼型气动力数值模拟

在Reynolds-Averaged Navier-Stokes(RANS)方程计算中耦合Michel经验判据和Chen-Thy-son转捩模型对风力机翼型S809进行计算。由RANS方程求得翼型表面压力分布作为层流边界层方程求解的输入参数,然后使用Michel经验判据分析判断层流边界层的解得到转捩点的位置,使用Chen-Thyson转捩模型得到从层流到湍流的转捩过渡区,这样随着流场的迭代求解,求解器自动更新判断出转捩点位置。对S809翼型进行全湍流和耦合转捩判断的RANS方程计算,可以看到考虑转捩判断后得到的升阻力系数与实验值吻合较好,验证了方法的可行性。     

压缩系统过失速过程的数值模拟

基于描述轴流压气机及其系统的动态气动响应的一维非定常数学模型,通过数值求解,即可对系统的过失速过程进行数值模拟。以一台三级压气机为算例,数值模拟了发生于该系统中的喘振和旋转失速现象     

翼型动态失速的变下垂前缘控制数值模拟

用CFD方法对VR-12翼型可压缩动态失速的变下垂前缘控制概念进行了数值模拟研究.模拟结果显示,变下垂前缘控制能在最大升力下降不大的情况下,非常显著地降低最大阻力,减小俯仰力矩负峰值,明显改善动态失速的负面效应.在流动机理上,变下垂前缘控制完全消除了动态失速涡.对参数影响的模拟研究表明,变下垂前缘控制的优势对马赫数、简缩频率或下垂控制方式的变化不敏感,是一种健壮的翼型可压缩动态失速控制手段.      

雷诺数变化对翼型边界层发展及失速特性的影响

翼型失速及其边界层发展是飞行器设计中的基础科学问题,而雷诺数变化对其影响很大。针对后缘失速翼型,采用Menter k-ω SST模型及耦合扰动放大因子输运方程的转捩模型,进行雷诺数变化对层流-湍流转捩边界层特性和失速特性的影响分析。结果表明:雷诺数增大时,对于转捩边界层,当地涡量雷诺数增大,转捩前移且分离泡减小,流动能量耗散减小,翼型整体表面剪切效应增强,动能更充沛,流动自持能力增强,压力分布可以维持较长距离的梯度抵抗分离能力增强;因此雷诺数增大使翼型失速迎角提高、升力系数增加。     

翼型动态失速的数值研究

用不可压缩流动的求解算法，结合WilcoxDC提出的k－ω模式和k－ωSST湍流模式，对翼型的动态失速进行了数值模拟。通过对典型的振荡翼型轻失速和深失速算例的计算结果分析可以看出：（1）绕动态失速翼型的流场结构十分复杂，轻失速和深失速在流动特性上有很大区别。计算结果显示：轻失速主要是由于后缘分离引起，分离涡的影响范围主要是在后缘附近。而深失速则首先形成很大的前缘分离涡，该分离涡在翼型表面上运动，并诱发出二次分离涡，引起翼型升、阻力系数的显著变化。（2）对于动态失速的翼型绕流，k－ωSST湍流模式是较为有效的，计算出的气动力系数迟带曲线变化趋势与实验结果符合得比较好。     

边界层转捩的数值模拟

在数值模拟中将Abu-Ghannam&Shaw(AGS)转捩模型引入到湍流模型中,以考虑流动中的层流和转捩现象,并在此基础上,分别对低雷诺数情况下的平板和涡轮叶栅流动进行了模拟,并与实验结果进行了比较,表明加入转捩模型可以较好地描述边界层的转捩过程。同时,结果还表明AGS转捩模型对压力梯度变化较大情况下的模拟稍有偏差。      

失速限制型风轮叶片流场的数值模拟

采用双时间步方法在旋转坐标系下求解非定常方程,构造了零偏航条件下风轮叶片流场的数值模拟方法,采用该方法对失速限制型风轮叶片NREL Phase-Ⅵ进行了模拟,其中网格采用重叠网格进行模拟,控制方程为可压缩N-S方程,同时加入了低速预处理技术,模拟的风速范围从13m/s到25m/s,这时候风轮叶片发生失速以限制扭矩输出.本文采用数值方法对风轮叶片在发生失速以后的主要流场特点进行了分析,从风轮扭矩、载荷分布等方面对数值计算结果进行了对比验证.     

基于粒子群算法的海上风机叶片优化算法研究

提出了一种针对不同风场条件提高海上风力发电机叶片气动效率的方法.风力发电机叶片是风机效率的核心部件,叶片在不同海上风场条件下功率差别很大,这导致了风机效率的降低.通过引入改进粒子群优化算法(PSO),基于修正动量叶素理论(BEM),对叶片沿展向的弦长分布和扭转角进行针对特定风场条件的气动优化.以美国再生能源实验室提供的5MW海上风机叶片作为算例,结合我国东海风场条件,叶片优化后平均功率提高了6.7％,取得了理想效果.结果表明,该优化方法具有较高的气动优化效率.     

大型风力机复合材料叶片动态特性及气弹稳定性分析

采用参数化建模技术快速建立大型风力机复合材料叶片三维有限元壳模型,并在此基础上对叶片的固有动力学特性进行停机及以额定转速旋转两种工况下的模态分析,其中旋转工况考虑了离心力导致的应力刚化效应和旋转软化效应。通过编制插值程序,将CFD计算所得的叶片表面分布压力,导算到叶片结构计算的有限元壳模型上,并以此为载荷对叶片进行静气弹稳定性分析。以某初步设计的1.5MW风机为例的计算结果表明:在参数化三维壳模型建模基础上进行的模态分析技术与结合CFD的静气弹稳定性分析技术,有利于快速、准确地计算大型复合材料叶片的动态特性、识别叶片结构的薄弱部位,并预测叶片发生局部屈曲的可能性及其发生的位置。     

尾缘改型对风力机翼型性能的影响研究

针对风力机工程常用翼型,采用有限元法的SIMPLE算法,对NACA4412翼型、加装2%弦长Gurney襟翼的NACA4412翼型及对应尾缘厚度为2%弦长的钝尾缘翼型进行了以来流风速为9.5m/s、攻角为0°～25°的气动性能数值计算。结果表明,改型后的钝尾缘翼型在流场中产生了强烈的下洗作用,明显改变了翼型表面的压力分布,使其升力系数、升阻比等较原型有显著改善,同时大大推迟了翼型的失速现象。通过翼型噪声机理的研究,对上述三种翼型分别采用在时间域上积分的方式进行了频谱特性的分析,气动声学研究表明,翼型噪声具有很强的指向性,改型后的翼型声级有明显降低,为低噪声风力机的优化设计和噪声预测提供了可靠的理论依据。     

大型风力机三维空气动力学数值模拟

采用数值模拟的方法,对三种风力机叶片的空气动力学性能进行研究.计算得到了旋转功率随着来流速度的变化曲线,并将所得结果与基于动量叶素理论的工程设计方法结果、风洞实验结果和风场实测结果进行对比.计算结果体现了CFD方法解决这种问题的有效性.所用的数值模拟方法可以广泛地应用到风力机设计和气动性能评价中,所得结果可作为叶片动态特性及气弹稳定性分析的载荷而被应用于风力机性能和可靠性的评估当中.     

旋翼翼型非定常动态失速响应的计算

基于旋翼非定常翼型气动模型,给出了计算分离流和深度失速状态下的翼型非定常升力、俯仰力矩的数值计算方法。该方法采用半经验指数响应公式,利用数值离散方法来求解翼型的非定常法向力和俯仰力矩。分别计算了NACA0012和SC-1095翼型上的非定常气动载荷,并与可得到的试验结果进行了对比,验证了方法的有效性。文中还讨论了缩减频率和马赫数对动态失速响应的影响;然后,这个模型被改进以适用于后掠流下的翼型动态失速响应计算,分析了后掠角对翼型动态失速响应的影响。最后,得出了一些结论。     

考虑气动弹性的风力机叶片性能分析

考虑气动弹性对风力机叶片的影响,采用叶素-动量理论计算气动力,采用盒形梁理论计算结构变形,耦合静气动弹性平衡方程,建立了风力机叶片静气动弹性分析程序。本文运用该程序进行了多种风速下叶片载荷及风轮性能的计算,分析了气动弹性对原设计的影响。结果表明,对于兆瓦级风力机,在大风速情况下,气动弹性对风轮性能有着明显影响,并会造成气动载荷的重新分布,影响结构设计的准确性。该方法可用于对叶片气动设计与载荷计算方法进行气动弹性修正。     

单级轴流压气机失速起始的数值模拟

基于某亚声单级轴流压气机试验台孤立转子环境下预失速行为的研究成果,进一步数值研究了该试验台单级环境下的失速起始过程.数值探针的监测结果表明,失速起始过程中单级轴流压气机的流动失稳始于转子叶尖部位,这与前期针对该试验台流动失稳的测量结果是一致的.对数值模拟的瞬态流场进行分析表明,先兆波的型式已从孤立转子环境下的“突尖”型转变为单级环境下的“多个小失速团”型,说明静叶的存在改变了近叶尖载荷的周向分布规律.然而,“叶尖二次涡”的形成与运动依然是“多个小失速团”型先兆波浮现的流体动力学机理.      

叶顶吹气对低速轴流压气机间隙流动 影响的数值研究

为了评价叶片顶部吹气对间隙流动的影响,采用NUMECA软件对带有叶顶吹气的叶栅流场进行了计算,并将其与不带吹气的基准叶栅流场结果进行了对比、分析.结果表明:叶顶吹气的位置和角度都对间隙流动有着较大影响.按周向和径向方向吹气时对泄漏流量密度的分布有一定的影响,但对稳定工况范围的影响不大;按一定的轴向方向吹出时可抑制间隙引起的流动分离,进而增大压气机的失速裕度.