NACA4412翼型低速绕流数值计算中湍流模型对比

使用Spalart-Allmaras(S-A)、SSTk-ω、Gao-Yong 3种湍流模型对NACA4412翼型低速绕流进行了数值计算,研究了尾迹流动松弛效应。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。计算中对翼型尾缘流松弛效应进行了分析,比较了翼型表面压力系数、速度剖面、雷诺应力等的分布,3种湍流模型总体上能够较好地模拟NACA4412翼型低速绕流。SSTk-ω模型对流动细节及升力系数计算最好,Gao-Yong模型对翼型平均速度剖面及雷诺剪切应力分布计算最准确。     

平均攻角和振幅对振荡翼型气动特性的影响

采用一种基于k-ωSST模型和γ-θ转捩模型的雷诺平均N-S方程数值方法,对雷诺数Re=1.35×105下的NACA0012振荡翼型和静态翼型非定常流场和升力特性进行模拟,在缩减频率K=0.1的条件下研究了翼型振荡运动中平均攻角和振幅对平均升力系数的影响,并与静态翼型的升力特性及实验结果进行了对比。结果表明:当平均攻角小于临界攻角时,翼型的振荡运动会降低平均升力系数,当平均攻角大于临界攻角同时最小攻角小于临界攻角时,翼型的振荡可以提高平均升力系数。在平均攻角为12°~17°时,翼型振幅为6°左右时可获得最大平均升力系数,与静态翼型相比,平均升力系数可提高30%~45.7%。当振荡过程中最小攻角对应静态翼型轻失速攻角时,翼型上仰阶段前缘涡的产生和集中涡的稳定附着是平均升力系数大幅度阶跃式提升的原因,静态翼型与振荡翼型的组合可提高升力并拓宽攻角范围。     

Gurney Flap绕流气动特性研究

采用C-H型多块结构网格和SST k-ω湍流模型求解二维定常雷诺平均NS方程,对NACA4412单段翼型后缘附近和NACA63(2)-215B多段翼型襟翼后缘附近的Gurney Flap绕流进行了数值模拟。讨论了局部网格对计算结果的影响,重点研究了Gurney Flap位置和高度对气动特性的影响。计算结果表明:对于NACA4412单段翼型,Gurney Flap位置对气动性能的影响较小,但随着高度的增加,零升攻角和失速攻角减小,最大升力系数增大,同时阻力系数和力矩系数增加;对于NACA63(2)-215B多段翼型,Gurney Flap对气动性能的影响类似于NACA4412单段翼型情形。     

S1223翼型低雷诺数下分离泡及气动特性

采用表面测压技术,测量了低雷诺数下(Re=6.0×104、1.0×105、2.0×105)S1223翼型的表面压力分布,通过时均化处理及瞬态处理方法,分别获得了翼型稳态和瞬态压力系数、升力系数,分析了流场结构随雷诺数及攻角的变化规律,研究了雷诺数及攻角对翼型升力的影响机理.结果表明,从时均升力系数随攻角的变化规律来看,S1223翼型在低雷诺数下存在"静态滞回"效应.攻角由负逐渐增大至0°时,下翼面由完全分离转变为出现层流分离泡,随后分离泡逐渐减小直至消失,导致升力系数斜率呈现随攻角逐渐增大的非线性现象.当攻角超过临界攻角后,不同雷诺数下翼型流场结构随攻角的变化规律存在本质不同,Re=6.0×104和1.0×105时,翼型周围流场迅速发生大范围流动分离,升力系数迅速减小;而Re=2.0×105时,上翼面周期性生成短泡,引发低频振荡现象,升力系数呈现准周期性变化,α=16°时上翼面时均流场呈现40％弦长的长泡结构.     

使用GAO-YONG湍流方程组计算翼型分离流

采用SIMPLE方法求解GAO-YONG不可压湍流方程组, 对不同来流迎角下的NACA0012翼型绕流结构进行了数值模拟, 给出了翼型绕流分离流结构随迎角的变化特征和翼型在分离绕流中的气动力参数.与实验数据以及大涡模拟结果的比较表明, GAO-YONG不可压湍流方程组能够对翼型绕流的分离点、分离涡形态、表面压力分布、升阻特性做出较好的预测, 能够模拟翼型大攻角分离流动, 计算结果优于FLU-ENT软件中的k-ε、k-ω、k-ωsst模型的计算结果.      

低雷诺数下二维翼型绕流的流场特性分析

采用高精度有限差分格式,对低雷诺数下二维翼型绕流进行了直接数值模拟,计算了雷诺数为1.0×104,NACA0012翼型0°,4°以及10°攻角下的流场。计算结果表明:在0°和4°攻角条件下,翼型绕流尾迹区的统计特性相似,0°攻角下的统计量值具有很好的对称性;在距翼型尾缘0.3弦长以后的尾迹区,旋涡排列成类似涡街的结构,涡量的极值、压力的极小值都位于旋涡中心,沿着流向,涡量极值的绝对值逐渐减小,压力的极小值逐渐增大。10°攻角下,翼型上表面从前缘开始分离,尾迹区统计分析结果所得图象与0°和4°完全不同,数值上较后者结果大;在翼型尾缘处,涡量的卷吸,高压力区域的形成,是旋涡脱落的条件,正向和反向旋涡的交替脱落,形成了类似涡街的结构。      

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

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

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

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

使用GAO-YONG湍流方程组对翼型绕流的计算

将GAO-YONG湍流方程组应用到翼型复杂绕流场的计算中,研究了机械能方程中Du/Dt项和各向异性湍流黏性系数分量中的实度系数Cs的影响,给出了控制Du/Dt项作用区域和选择实度系数Cs值的判别函数,对NACA4412翼型最大升力迎角下的低速绕流算例、AS240翼型8°迎角下的转捩算例、AS240翼型19°迎角下的分离流算例进行了模拟,与实验结果比较,得到了比较好的结果.结果表明GAO-YONG湍流方程组能够准确预测翼型前缘层流到湍流的转捩问题,翼型尾迹流动松弛问题和翼型大迎角分离流.     

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

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

Numerical simulation of airfoil flows with a turbulence model

Turbulent flows over AS240 and NACA4412 airfoil were simulated numerically using a two-equation turbulence model named k-ξ model. The predictions of velocity profiles and the pressure coefficient of airfoil AS240 at 8°/19° attack angle and NACA4412 at 13.87° attack angle were calculated. The results were compared with those using k-ε and k-ω models,as well as experimental data. It indicates that the new k-ξ model offers more realistic prediction than the other two models. The main finding shows that the new k-ξ model is good at predicting separated flows around airfoils, and it captures the flow feature of pressure-induced separation adequately. All calculations are implemented as per openFOAM 1.7.1(open source field operation and manipulation).      

大湍流度对超低雷诺数下翼型受力及绕流的影响

实验测量了超低雷诺数(Re=5 300)下 NACA 0012翼型在自由来流下的升力系数和阻力系数,重点研究来流的湍流度对升力系数和阻力系数的影响,并进一步通过对翼犁流场的研究揭示湍流度对翼型受力产生影响的机理.结果显示,低湍流度下升力系数和阻力系数均无失速特征;当湍流度提高,升力系数和阻力系数在12°～15°迎角下表...     

翼型失速的非线性动力学特性及其控制

以NACA0012翼型为研究对象,应用CBS( Characteristic-Based Split)有限元方法对翼型绕流流场进行数值模拟,着重对翼型失速现象的动力学机理和特征进行研究.首先给出了静止翼型在不同迎角下流场的瞬时流线图以及静态翼型的迎角-升力系数关系曲线,并且分析出此曲线具有滞后、突跳等典型的非线性动力学...     

翼型结冰过程数值模拟验算与分析

应用FENSAP-ICE结冰计算软件,对NACA0012翼型进行了流动特性、水滴撞击特性以及冰型生成过程的计算;同时,对结冰前后的翼型进行气动力特性计算对比分析,其中包括升力特性对比、阻力特性对比、流场细节分析以及压力系数分布对比。计算结果表明：翼型前缘结冰后,导致翼型前缘气流提前分离,最大升力系数、失速攻角大幅减小,...     

翼型湿空气非平衡凝结流动的数值研究

 湿空气非平衡凝结会对机翼的气动特性造成显著影响。通过对Fluent软件进行二次开发,发展了一种求解湿空气非平衡凝结流动的数值方法并进行了验证。对NACA0012翼型在迎角为0°和1.5°时的湿空气非平衡凝结流动进行了分析。结果表明：与干空气流动相比,在30％～70％相对湿度范围内,湿空气非平衡凝结流动中升力系数显著降低,而阻力系数则可能增加或减小;0°迎角条件下,压差阻力系数最大增加了68.3％;在迎角为1.5°时,升力系数最大减小39.8％,而压差阻力系数最大减小了34.8％。造成翼型气动特性显著变化的原因在于：湿空气中水蒸气凝结放热对跨声速气流加热,导致翼型表面附近的流速、压力与流场结构发生了显著变化。     

Experimental study of ice accretion effects on aerodynamic performance of an NACA 23012 airfoil

In this paper,the effects of icing on an NACA 23012 airfoil have been studied.Experiments were applied on the clean airfoil,runback ice,horn ice,and spanwise ridge ice at a Reynolds number of 0.6 106 over angles of attack from 8° to 20°,and then results are compared.Generally,it is found that ice accretion on the airfoil can contribute to formation of a flow separation bubble on the upper surface downstream from the leading edge.In addition,it is made clear that spanwise ridge ice provides the greatest negative effect on the aerodynamic performance of the airfoil.In this case,the stall angle drops about 10° and the maximum lift coefficient reduces about50% which is hazardous for an airplane.While horn ice leads to a stall angle drop of about 4° and a maximum lift coefficient reduction to 21%,runback ice has the least effect on the flow pattern around the airfoil and the aerodynamic coefficients so as the stall angle decreases 2° and the maximum lift reduces about 8%.     

低雷诺数下协同射流关键参数对翼型气动性能的影响

协同射流是一种高效的新型主动流动控制技术,至今缺乏关键参数对翼型气动性能影响规律的系统研究。通过引入螺旋桨激励盘模型,发展了一种新的协同射流翼型流动模拟方法,使得射流反作用力计算更符合实际。在低雷诺数条件下,以NACA6415为基准翼型开展了射流动量系数、开口尺寸和位置等关键参数对翼型气动性能的影响规律研究,并探讨了相应的物理机制。结果表明：大迎角分离流状态下,射流动量系数对翼型气动性能的影响规律比小迎角附着流状态更复杂;随着吹气口尺寸增加,气泵功率系数先减后增,有效升阻比先增后减;随着吸气口尺寸增加,气泵功率系数逐渐减小,有效升阻比先增加后趋于平稳;吹/吸气口位置对翼型气动性能和气泵功率系数的影响很小。     

高速风洞等离子体流动控制实验技术研究

针对开展等离子体高速流动控制研究的技术需求,通过专用模型及实验机构设计、绝缘密封走线、多层电磁屏蔽等技术手段,建立了一套适用于高速风洞的等离子体流动控制系统,提出了等离子体高速流动控制风洞实验的技术规范和运行策略,并初步探索了等离子体激励对二元翼型绕流的控制规律。采用该技术后,解决了高压电缆的绝缘、密封走线问题,模型与实验机构的感应电压减小90％以上。风洞实验结果表明：实验系统运行稳定,实验数据可靠,等离子体激励对犕犪＝0．2的流动可实现有效控制;施加等离子体激励后,NACA0012翼型的流动分离明显减弱,升力增大,阻力减小,临界失速迎角增大2＆#176;,最大升力系数增大4％,总体气动性能得到显著提升。     

低雷诺数下吹吸气射流对NACA0012翼型气动性能的影响

基于Favre过滤的大涡模拟方法,对雷诺数Re=10^4,迎角α=6°下的NACA0012翼型上表面吹吸气射流进行了数值模拟,从翼型周围流场流线图、速度场云图、上下表面压力系数曲线以及上表面边界层位移厚度等多角度地分析了射流位置以及速度变化对翼型气动性能的影响。结果表明:射流位置对翼型气动性能影响较大,且吸气射流要明显优于吹气射流。对于吸气射流,前缘吸气要明显优于中后缘吸气,可有效增升减阻,并减小翼型尾部流动分离,抑制翼型气动参数扰动,其最佳吸气位置随着速度的增大逐渐向下游移动;而吹气射流对翼型气动系数的作用效果较差,但中后缘的吹气射流可减小飞行过程中的气动扰动量,且吹气越大,效果越明显。