风力机专用翼型的反设计

本文根据翼型表面压力分布对边界层的影响,分析翼型表面压力分布的特点,以XFOIL计算软件为设计平台,采用其中的混合反设计模块,通过合理改变翼型表面的速度分布来得到气动特性满意的翼型。所设计的翼型主要要求具有高的设计升力系数、高的最大升阻比、良好的前缘粗糙不敏感性及和缓的失速特性。在提高设计升力系数的同时限制最大升力系数,以减小两者的差值,减小叶片的极限载荷。经计算分析,采用反设计优化得到的翼型与DU翼型相比具有很好气动特性。     

凹槽对风力机叶片尾缘襟翼性能的影响

本文采用单方程S—A湍流模型计算了尾缘襟翼与主体翼型连接处的凹槽对二维翼型气动特性及流场的影响。选用带有30％弦长固定偏斜20°角度尾缘襟翼的NACA0015翼型作为研究对象,分析了流场的相关特性。数值计算结果表明：凹槽对于带有尾缘襟翼的二维翼型气动特性有一定影响,在负攻角及较小正攻角时,减小翼型升力系数;随着攻角增大...     

减缩频率和平均攻角对俯仰振荡翼型影响分析

本文以NACA0012翼型为研究对象,采用混合网格划分方法和SST κ-ω湍流模型,数值模拟了雷诺数Re=2.7×10~5条件下减缩频率和平均攻角对翼型俯仰振荡气动特性的影响。结果表明:翼型俯仰运动时的平均升力系数均低于静态条件下的升力系数;减缩频率对翼型下行段气动特性影响最为显著,当减缩频率较小时,翼型刚开始下行运动,出现流动分离越显著,这导致平均升力系数与静态条件下升力系数差值变大;平均攻角越大,俯仰运动时的最大升力系数越小;翼型俯仰运动上行段升力系数大,主要是因为前缘流动加速剧烈,增大了上下表面压差。     

基于AC-DBD等离子体激励的NACA0012翼型动态失速控制

以NACA0012翼型为研究对象,采用动态测压及PIV测量技术,研究了AC-DBD等离子体激励器对翼型俯仰及耦合运动动态失速的控制作用和机理.研究表明,等离子体激励能够显著推迟失速迎角,抑制失速后的升力系数陡降,提前流动再附和升力系数回升,减小升力及俯仰力矩系数曲线迟滞环面积,改善翼型气动特性.研究了不同运动参数及激励器设置参数对控制效果的影响,结果表明翼型俯仰运动频率及激励器激励频率分别对激励器控制效果影响最大,为后续相关研究提供了数据基础.      

风力机翼型尾缘襟翼动态特性分析

采用计算流体力学方法研究了带有运动尾缘襟翼的风力机翼型,考察了襟翼偏转角频率对翼型气动参数及非定常特性的影响。结果表明:多数情况下,翼型升力系数滞后于偏转角变化,且相位差随着角频率的增加先增大后减小;尾缘襟翼改变升力系数的能力随着角频率的增加而减小;以尾缘襟翼长度为特征尺度定义的襟翼折合频率可作为尾缘襟翼问题非定常特性的判断准则,当该折合频率大于或接近0.01时,流场具有明显的非定常特性。     

YLSG 107高升力翼型主动流动控制的数值模拟

数值模拟零质量射流与YLSG 107翼型绕流的干扰流场,探讨零质量射流在高升力翼型失速控制中的控制效果、控制特性及控制机理.数值模拟以积分形式雷诺平均Navier-Stokes(N-S)方程为控制方程,采用格心有限体积法进行求解.通过在喷口上施加非定常边界吹/吸边界条件模拟射流对翼型绕流的干扰.采用与风洞实验相同的来流状态和控制参数进行数值模拟,得到与实验相吻合的结果.为进一步研究控制特性和控制规律、提出改进的实验方案,研究不同动量系数、不同射流偏角对控制效果的影响,并对法向射流和近切向射流进行较深入的比较.研究表明,先前的风洞实验对应的射流动量系数(0.000 014)偏小是控制效果不显著的重要原因之一,必须达到0.001以上才有明显控制效果(射流动量系数为0.005时可使该翼型失速迎角增大2°,最大升力提高8.7%);近切向射流在失速控制方面明显优于法向射流.     

波状前缘对风力机翼型失速性能影响数值分析

为研究仿生波状前缘对翼型失速性能的影响,本文采用S-A湍流模型,对风力机翼型NACA634-021(光滑前缘)以及对应的正弦波状前缘仿生翼型的绕流流场进行了数值模拟。结果表明,光滑翼型在20°攻角附近发生深度失速,升力系数骤然下降;而波状前缘仿生翼型有效改善了失速特性,升力系数变化较平稳,在大攻角下高于光滑翼型。通过流场分析发现光滑翼型失速前后升力系数骤然下降的主要原因在于前缘压力面和吸力面的压差大幅度下降,而仿生翼型改变了前缘的压力分布特性,进而改变了大攻角下的分离特性,促进流向涡对的产生和发展,使得凸峰附近保持附着流动,进而提高升力。     

水平轴风力机专用翼型族试验分析及优化设计

对水平轴风力机专用翼型族—CAS-W1-XXX薄翼型族试验结果进行了分析,并将其与国外同等厚度翼型进行对比。试验结果表明,与国外同等厚度翼型相比CAS-W1-XXX薄翼型具有良好的前缘粗糙不敏感性、高的最大升力系数、设计升力系数和良好的失速特性。为进一步提高翼型的气动特性,在试验结果的基础上对CAS-W1-XXX薄翼型族进行再次优化。根据XFOIL计算结果,优化后翼型的最大升阻比得到提高,并且与DU翼型相比具有良好的气动特性。同时对CAS-W1-XXX厚翼型中出现的小攻角失速现象进行了优化改进。     

大厚度钝尾缘翼型的设计研究

本文研究了MW级水平轴风力机叶片内侧翼型的设计准则,采用混合设计方法得到了四种适用于该部位的大厚度钝尾缘翼型:依据一个5 MW风力机的运行特征和性能需求,在不同的设计雷诺数下,以翼型在运行攻角范围内的升力水平为主要气动目标,以升力随雷诺数变化的稳定性为主要约束。RFOIL预测表明,在15°~30°攻角之间,四种翼型的升力系数大致从1.0以近似线性方式增加到1.7,升力曲线随着雷诺数变化稳定,具有良好的变工况特性。     

基于简易风洞的薄翼模型升力测量实验

为加深学生对薄翼型升力理论相关知识的理解,使抽象的理论直观化,我们设计了低速简易风洞测量薄翼模型升力系数的实验:利用简单的弹簧测力装置测量不同迎角下翼型的升力系数。实验数据表明,给定薄翼模型的升力系数与来流速度无明显关系、和迎角近似线性。多次实验结果表明,该简易风洞实验可用于半定量验证薄翼型升力理论,也可用于研究其他自制翼型的升力系数。     

Numerical study on the correlation of transonic single-degree-of-freedom flutter and buffet

Transonic single-degree-of-freedom(SDOF) flutter and transonic buffet are the typical and complex aeroelastic phenomena in the transonic flow. In this study, transonic aeroelastic issues of an elastic airfoil are investigated using Unsteady Reynolds-Averaged Navier-Stokes(URANS) equations. The airfoil is free to vibrate in SDOF of pitching. It is found that, the coupling system may be unstable and SDOF self-excited pitching oscillations occur in pre-buffet flow condition, where the free-stream angle of attack(AOA) is lower than the buffet onset of a stationary airfoil. In the theory of classical aeroelasticity, this unstable phenomenon is defined as flutter. However, this transonic SDOF flutter is closely related to transonic buffet(unstable aerodynamic models) due to the following reasons. Firstly, the SDOF flutter occurs only when the free-stream AOA of the spring suspended airfoil is slightly lower than that of buffet onset, and the ratio of the structural characteristic frequency to the buffet frequency is within a limited range. Secondly, the response characteristics show a high correlation between the SDOF flutter and buffet. A similar "lock-in" phenomenon exists, when the coupling frequency follows the structural characteristic frequency. Finally, there is no sudden change of the response characteristics in the vicinity of buffet onset, that is, the curve of response amplitude with the free-stream AOA is nearly smooth. Therefore, transonic SDOF flutter is often interwoven with transonic buffet and shows some complex characteristics of response, which is different from the traditional flutter.     

基于进化算法和复合参数化方法的低速翼型优化

针对设计工况为Re=1×106的翼型在2°和8°两种迎角下的气动性能进行综合优化, 提升其在巡航与起飞着陆时的性能。利用气动估算软件Xfoil结合NACA 4位数参数化方法进行翼型初步优化, 缩小优化空间, 减少计算量。通过求解Reynolds平均Navier-Stokes方程的方式进行精确优化, 最后利用Hicks-Henne形函数法对翼型表面施加扰动进行细节优化, 弥补NACA 4位数参数化方法细节表现不足的缺点。实现了2°和8°两种迎角下气动性能普遍提升20%以上的效果, 并发展了一种复合参数化方法下基于改进进化算法的多目标、多步骤气动优化方法。      

任意马赫数非定常流动数值模拟的统一算法

发展适用于从低速到高速任意马赫数非定常流动数值模拟的统一算法.通过引入一个伪时间导数项和一个新的预处理矩阵,得到双时间非定常预处理可压缩Navier-Stokes方程.方程的对流项采用三阶Roe通量近似差分格式离散,粘性项采用二阶中心差分格式离散.基于数值通量的线性化技术,实现伪时间步的隐式ADI-LU格式迭代,进而获得物理时间步的二阶推进精度.重点以低马赫数流动为例,求解了圆柱绕流和NACA0015翼型等速上仰动态失速问题.计算结果表明该统一算法能够较好地模拟低马赫数乃至任意马赫数非定常流动.     

Preconditioned characteristic boundary conditions for solution of the preconditioned Euler equations at low Mach number flows

Preconditioned characteristic boundary conditions (BCs) are implemented at artificial boundaries for the solution of the two- and three-dimensional preconditioned Euler equations at low Mach number flows. The preconditioned compatibility equations and the corresponding characteristic variables (or the Riemann invariants) based on the characteristic forms of preconditioned Euler equations are mathematically derived for three preconditioners proposed by Eriksson, Choi and Merkle, and Turkel. A cell-centered finite volume Roe’s method is used for the discretization of the preconditioned system of equations on unstructured meshes. The accuracy and performance of the preconditioned characteristic BCs applied at artificial boundaries are evaluated in comparison with the non-preconditioned characteristic BCs and the simplified BCs in computing steady low Mach number flows. The two-dimensional flow over the NACA0012 airfoil and three-dimensional flow over the hemispherical headform are computed and the results are obtained for different conditions and compared with the available numerical and experimental data. The sensitivity of the solution to the size of computational domain and the variation of the angle of attack for each type of BCs is also examined. Indications are that the preconditioned characteristic BCs implemented in the preconditioned system of Euler equations greatly enhance the convergence rate of the solution of low Mach number flows compared to the other two types of BCs.     

进入轨道偏差对火星科学实验室气动力特性的影响

针对火星科学实验室(MSL)高超声速进入过程,利用三维并行程序求解流体力学Navier-Stokes方程,耦合真实气体模型,分析火星大气中真实气体效应对进入器气动力特性的影响量在进入轨道发生偏差时的变化规律.结果表明:对海盗号的计算结果与飞行数据符合很好,验证了火星大气真实气体模型和计算方法;真实气体效应影响下,激波层厚度大为减小,温度下降明显,进入器阻力系数明显增加,升力系数变化不大,俯仰力矩系数增加,基准状态下配平攻角较完全气体减小约2.2°;高度不变,Ma数增加导致阻力系数和俯仰力矩系数增大,配平攻角和完全气体的差值由1.6°增加到2.6°,表明Ma数变大时真实气体效应引起的气动力变化增强;Ma数不变,高度增加略微减弱波后化学反应,对进入器气动力特性基本没有影响.     

Simulation of Viscous Flows Around a Moving Airfoil by Field Velocity Method with Viscous Flux Correction

Based on the field velocity method, a novel approach for simulating unsteady pitching and plunging motion of an airfoil is presented in this paper. Responses to pitching and plunging motions of the airfoil are simulated under different conditions. The obtained results are compared with those of moving grid method and good agreement is achieved. In the conventional field velocity method, the Euler solver is usually used to simulate the movement of the airfoil. However, when viscous effect is considered, unsteady Navier-Stokes equations have to be solved and the viscous flux correction must be taken into account. In this work, the viscous flux correction is introduced into the conventional field velocity method when non-uniform grid speed distribution is occurred. Numerical experiments for the flow around NACA0012 airfoil showed that the proposed approach can well simulate viscous moving boundary flow problems.     

A discontinuous Galerkin method for inviscid low Mach number flows

In this work we extend the high-order discontinuous Galerkin (DG) finite element method to inviscid low Mach number flows. The method here presented is designed to improve the accuracy and efficiency of the solution at low Mach numbers using both explicit and implicit schemes for the temporal discretization of the compressible Euler equations. The algorithm is based on a classical preconditioning technique that in general entails modifying both the instationary term of the governing equations and the dissipative term of the numerical flux function (full preconditioning approach). In the paper we show that full preconditioning is beneficial for explicit time integration while the implicit scheme turns out to be efficient and accurate using just the modified numerical flux function. Thus the implicit scheme could also be used for time accurate computations. The performance of the method is demonstrated by solving an inviscid flow past a NACA0012 airfoil at different low Mach numbers using various degrees of polynomial approximations. Computations with and without preconditioning are performed on different grid topologies to analyze the influence of the spatial discretization on the accuracy of the DG solutions at low Mach numbers.     

Optical-velocimetry, wake measurements of lift and induced drag on a wing

Particle Tracking Velocimetry was used, in a low-speed wind-tunnel study, to obtain simultaneous cross-flow velocity measurements in three planar regions downstream of an airfoil having a NACA 0015 profile. In order to measure both the total lift and induced drag, and their distributions on the airfoil, a wake integration technique based on the control volume approach was used. The airfoil model was tested in clean wing configuration and with a 3.33% Gurney flap attached to the trailing edge. The Gurney flap was found to increase the complexity of the wake and cause a systematic increase in lift and induced drag values. The changes in total lift and induced drag resulting from the fitting of a Gurney flap compared very favourably to investigations using traditional aerodynamic techniques demonstrating the validity of the method.     

二维地效翼及地效流动特性数值研究

用数值模拟的方法,对二维NACA0012翼型在地面效应下的空气动力特性和地效流动特性进行研究,得到地效翼的升力、阻力和翼型表面压力分布随攻角及相对飞行高度的变化情况.通过对计算结果的分析,可以看出,在一定的攻角,靠近地面飞行,机翼的升力得到提高;随着飞行高度的降低,地面效应增强,机翼的失速攻角减小;地面附近的粘性流动对机翼的空气动力特性影响很小;当相对飞行高度小于0.1时,应该考虑空气的可压缩性.