NLR-7301高升力翼型绕流评估计算

探讨了在二维高升力多段翼型绕流计算中附面层网格采用各向异性三角形网格的技术。通过对经典的NLR-7301两段翼型绕流的计算,从数值计算角度评估了自研混合网格雷诺平均N-S方程求解软件对该类型附面层网格的解算能力。通过与确认试验数据对比,在附面层采用各向异性三角形网格进行二维高升力多段翼型绕流计算,其结果具有较高的可信度。在此基础上,采用同样策略,评估了多段翼型绕流计算中翼型后缘的几何处理、湍流模型选取等关键CFD要素对计算结果的影响。     

多段翼型大迎角下主翼、襟翼上的分离流及缝道流动

使用雷诺平均N-S方程、采用可用于较大分离区的Johnson-King紊流模型、嵌套网格和有限体积法研究大迎角下的多段翼型绕流,特别是主翼、襟翼上的分离流动及缝道流动。利用嵌合体技术对组合体每一部分生成高质量并适于高效求解的贴体网格。以具有17%相对厚度的GAW-1翼型带30%襟翼翼型为例进行了计算,计算结果与实验结果吻合很好,证实该方法可以较好地预计多段翼型上的分离流、缝道流动与最大升力。     

基于压力分布的开孔翼型气动特性数值模拟研究

以NACA0012为基准翼型,生成C型结构网格,采用SA湍流模型求解NS方程,获得基准翼型压力分布数值模拟结果。找出翼型上表面前缘低压区域和后缘高压区域,在两个区域进行开孔,让后缘高压区域气流通过导流管引流到前缘区域,加速翼型上表面气流流动,增大翼型速度环量。通过基准翼型与开孔翼型气动特性对比分析可知:小迎角下,基于压力分布的开孔翼型阻力系数略有增大,但能有效提高翼型的升力系数。     

求解平面翼型亚,跨声速绕流的一个新方法

本文提出了求解平面翼型亚、跨声速绕流的一个新方法。引入流函数和Ｖｏｎ Ｍｉｓｅｓ变换后，亚、跨声速绕机翼无旋流动的基本方程组被化为以流线纵坐标ｙ为未知量的单个二阶偏微分方程－流线控制方程。并通过变换将物理平面上的无限域变为计算平面上有限的矩形域，而后在计算平面采用有限差分线松弛迭代法求解。作为算例，计算了对称翼型ＮＡＣＡ００１２－３４和非对称翼型ＮＡＣＡ４４１２的亚、跨声速有攻角绕流，所得数值结果     

二维弹性扑翼沉浮运动流动特性的数值研究

采用非定常势流板块法和Euler-Bernoul1i梁振动微分方程,对后缘带弹性薄板的翼型在作沉浮运动时翼型绕流和薄板弹性振动之间的流-固耦合过程进行数值模拟,并计算流动特性、非定常气动载荷,尤其是沉浮运动的推力效应.用有限差分法求解薄板的弹性变形运动,气动弹性分析采用弱耦合迭代方式.通过对翼型非定常气动力(升力和推力)、推进效率、尾迹脱落涡计算结果的分析发现,通过调整弹性薄板的刚度,可以使得翼型产生最佳的非定常推力.     

旋转坐标系N-S 方程无网格/网格混合算法

为了模拟涡桨发动机等绕固定轴旋转的桨叶流场,发展了一种求解旋转体非定常黏性绕流的无网格/网格混合算法。算法基于求解旋转坐标系纳维尔-斯托克斯方程展开,避免了旋转角速度特征物理量的插值运算;计算域采用整体网格和物面附近局部无网格离散,通过引入无网格和网格对偶点,实现了混合算法绕流信息的跨区交换,并给出了一种对偶点调整选取的方法;基于无网格点云空间导数逼近方法离散控制方程,并结合隐式LU-SGS算法,给出了混合算法求解非定常问题的双时间步隐式推进格式。结果表明:所提算法通过2维振荡翼型非定常绕流、旋转圆柱黏性绕流和模拟发动机桨叶旋转运动的3维悬停旋翼绕流进行了考核,所得升力系数等重要气动数据的变化趋势与文献试验或计算值一致,典型截面处所示的激波强度和位置亦与文献值吻合,展现出算法在模拟旋转部件绕流问题方面具有广阔的工程应用前景。     

Gurney Flap绕流气动特性研究

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

超临界翼型尾缘噪声影响因素研究

使用大涡模拟(LES)和FW-H相结合的方法对超临界翼型尾缘噪声进行了研究。针对超临界翼型后缘的不同设计参数(后缘厚度、后缘角和几何形状)对尾缘噪声的影响,使用ANSYS FLUENT的LES湍流模型计算声源,采用FW-H积分方法求解远场噪声总声压级。首先完成二维非定常流圆柱绕流的流场及其噪声的验证计算,计算结果与实验值符合,证明了求解器设置和网格生成的合理性。然后基于此正确的求解器设置和网格生成,使用LES/FWH对比了典型超临界翼型的不同后缘设计参数对远场总声压级的影响,所得结论对低噪声超临界翼型优化设计有参考意义,同时为进一步的噪声控制优化设计提供了基础。     

超临界翼型射流环量控制的数值计算

为提高飞行器的升力特性,改善起飞着陆和机动性能,研究了射流环量控制对超临界翼型气动特性的影响。基于计算流体力学,对影响计算效率和精度的网格参数进行了优化设计,得到一种高效高精度网格的剖分方法;并以此为基础,建立超临界翼型环量控制的计算模型,求解雷诺平均N-S方程,分析了环量控制在不同动量系数和迎角下的作用效果,数值计算结果与实验结果基本一致。结果表明:环量控制大大提高了翼型的气动特性,且在小迎角或负迎角就可以产生较大的升力,但是动量系数并不是越大越好。计算结果可为进一步研究环量控制对飞行器的影响奠定理论基础。     

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

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

机翼和叶栅非定常流的Kutta条件

本文在现有的实验研究和数值计算的基础上，对经典Ｋｕｔｔａ条件在非定常流动中的适用性作了进一步探讨，并对文献〔１，２〕中提出的一种新的Ｋｕｔｔａ条件－Ｇｉｅｓｉｎｇ－Ｍａｓｋｅｌｌ（Ｇ－Ｍ）模型进行了分析与比较，本文给出了二维振荡机翼和叶栅绕流的变域变分有限元解法中Ｋｕｔｔａ条件的处理方法。本文旨在为机翼和叶栅非定常绕流的数值计算提供更合理可靠的尾缘边界条件。     

NPU翼型的气动力分析和改进设计

 在飞行器设计中用计算方法设计超临面翼型已完全取代了选用现成翼型的设计方法。为考察已设计出的NPU翼型是否满足飞行器设计要求我们对其进行了全面气动分析,发现这些翼型尚有不足之处,有必要进行改进设计。     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

不可压缩势流绕多段翼型流动的快速迭代解法

利用保角变换计算不可压缩势流绕翼型的流动,具有计算速度快和要求计算机内存小的优点,但主要用于单段翼型问题,而对于多段翼型,尤其是三段和三段以上的翼型则很难应用。本文在单段翼型的保角变换的基础上,提出了一种利用迭代算法计算任意多段翼型的不可压缩势流的绕流问题的方法。本文方法和目前常用的边界元法相比,具有计算速度快的优点。本方法还可进一步推广到考虑无粘流/边界层的迭代计算中。     

二维振荡机翼跨声速非定常绕流的变域变分有限元解

本文首先对文献（１－３）提出的二维振荡机翼含激波跨声速非定常绕流的变分原理进行了改造，使之能适用于时间推进法；构造了三维时空的可变节点有限元来捕获自由尾涡面，而跨声流中的激波采用人工密度法捕获。在远场边界上使用简化的无反射条件和新型非定常Ｋｕｔｔａ条件。用本文方法求解了作俯爷振荡的ＮＡＣＡ６４Ａ０１０翼型的跨声速绕流，计算结果令人满意。本文方法可以推广到三维机翼及二维和三维叶栅的同类气动力问题上去     

旋翼桨-涡干扰气动特性计算及参数影响研究

建立了一个适用于旋翼桨-涡干扰(BVI)气动特性分析的数值方法。在该方法中,控制方程采用惯性坐标系下的非定常Euler方程,以适合于分析BVI流场的特点;为便于前飞流场分区求解和信息传递,使用了嵌套网格方法;对于流场中涡线的模拟,建立了一种适用于有限体积格式的涡引入方法——广义网格速度法,以简化计算。应用上述方法对旋翼BVI流场进行了计算,并与可得到的试验数据进行对比,验证了方法的有效性。通过对比桨叶弦向位置压强的变化得出,在旋翼BVI过程中,气动载荷主要来自桨叶10%弦长内的前缘部分的压强突变。文中还进一步分析了涡强、干扰距离和干扰角度对BVI气动特性的影响。本文的BVI数值计算表明:当涡接近桨叶前缘时,升力达到最大;而涡位于后缘位置时,诱导速度改变桨叶环量分布,破坏了后缘的库塔条件,但随着涡的远离,桨叶表面环量开始重构。     

Self-sustained shock oscillations on airfoils at transonic speeds

Self-sustained shock wave oscillations on airfoils at transonic flow conditions are associated with the phenomenon of buffeting. The physical mechanisms of the periodic shock motion are not yet fully understood even though experiments performed over fifty years ago have demonstrated the presence of oscillatory shock waves on the airfoil surfaces at high subsonic speeds. The unsteady pressure fluctuations generated by the low-frequency large-amplitude shock motions are highly undesirable from the structural integrity and aircraft maneuverability point of view. For modern supercritical wing design with thick profiles, the shock-induced fluctuations are particularly severe and methods to reduce the shock wave amplitudes to lower values or even to delay the oscillations to higher Mach numbers or incidence angles will result in expanding the buffet boundary of the airfoil. This review begins with a recapitulation of the classical work on shock-induced bubble separation and trailing edge separation of a turbulent boundary layer. The characteristics of the unsteady pressure fluctuations are used to classify the types of shock-boundary layer interaction. The various modes of shock wave motion for different flow conditions and airfoil configurations are described. The buffet boundaries obtained using the standard trailing edge pressure divergence technique and an alternative approach of measuring the divergence of normal fluctuating forces are compared to show the equivalence. The mechanisms of self-sustained shock oscillations are discussed for symmetrical circular-arc airfoils at zero incidence and for supercritical airfoils at high incidence angles with fully separated flows. The properties of disturbances in the wake are examined from linear stability analysis of two-dimensional compressible flows. The advances in high-speed computing make predictions of buffeting flows possible. Navier–Stokes solvers and approximate boundary layer-inviscid flow interaction methods are shown to give good correlation of frequencies and other unsteady flow characteristics with experiments. Finally, passive and active methods of shock oscillation control show promising results in delaying buffet onset to higher Mach numbers or incidence angles, thus enhancing the transonic performance of airfoils.     

尾缘锯齿结构的降噪物理机制实验

对比分析了常规尾缘翼型与锯齿尾缘翼型尾缘湍流流场的基本特征,并通过线阵列的方法测量了两种尾缘结构的噪声.结果表明:锯齿尾缘翼型尾缘湍流流场的湍流强度以及3个方向上的湍流强度都相比于常规尾缘翼型有显著减少,声场结果显示锯齿尾缘翼型对尾缘噪声有显著减小,对前缘噪声影响很小.锯齿结构加宽了尾迹区域并加快了大涡的破碎,产生了额外的马蹄涡,湍流脉动衰减率沿着流动方向变大.      

机翼跨声速绕流的Galerkin有限元解法

本文给出了翼型和机翼跨声速绕流的Galerkin带权余数有限元解法。利用一系列适当形状的元素将计算区域离散,在元素内采用线性插值函数。对于超临界流动,文中采用了适合于有限元法的超声速上风技术,具有捕捉激波的能力。将这种技术应用于线松弛迭代解法中,在计算钝头翼型和机翼绕流时,都获得了成功。