65°三角翼亚音速复杂流场计算和数据可视化

提出一种基于非结构混合网格和有限体积法的有效计算策略,对第二期国际涡流试验项目(second international vortex flow experiment,VFE-2)的尖前缘65°三角翼在马赫数0.4,迎角20.3°,雷诺数2×10~6条件下的亚音速复杂流场结构进行数值模拟,重点探讨了基于计算数据进行该类型复杂涡系干扰表面和空间流场关键特征提取和数据可视化问题.通过与相关试验类比,建立了与先进试验流动显示技术相比拟的定性和定量分析方法,为三角翼这类复杂流场结构的精细分析奠定了技术基础.采用上述方法,细致分析了亚音速三角翼的大迎角复杂旋涡流场结构,得到了与试验一致的结论.研究证实:在大迎角条件下,三角翼流动物理复杂,黏性效应耦合严重,只有通过N-S方程计算才能准确地捕捉主涡和二次涡的发展.     

矢量喷流对细长体大迎角非对称流动影响研究

采用测压方法研究了矢量喷流对细长旋成体大迎角非对称流动的影响特性.实验结果表明:矢量喷流对细长旋成体大迎角非对称侧向力有明显的抑制作用,该抑制作用是通过喷流诱导作用,改变其空间绕流涡系结构的分布来实现的,但是矢量喷流的存在并不能改变大迎角机身空间绕流涡系的本质结构;随着迎角的增大,矢量喷流对细长旋成体大迎角非对称流动的影响区域不断前移,甚至影响到头部;随着喷流落压比的增加,矢量喷流对细长旋成体大迎角非对称侧向力的抑制作用加强,但当喷流落压比达到临界落压比后(即喷管出口处达到设计马赫数时),喷流影响作用将不会随喷流落压比的增加而改变.     

流动显示技术及其现代发展概况

本文回顾了流动显示发展的历史,以及它在流体力学研究中所起的作用。根据获取流动信息的方式不同,流动显示技术可分为两大类型。作为近代技术发展的例子,三维非对称流干涉图的定量计算、片光技术、散斑照相术、激光诱导萤光法、多普勒干涉图象显示速度场、光电信号转换激光干涉仪、以及图象处理和计算机产生数字流谱技术均作了简要的介绍。     

飞机大迎角非对称涡组合扰动主动控制研究

采用测压、测力以及流动显示方法研究了头部微三角扰动块对飞机大迎角非对称背涡的主控作用和背风侧单孔位微吹气对背涡空间位置及相应侧向力的控制作用,以此为基础提出了基于微三角块扰动和单孔位微吹气扰动的组合扰动主动控制新技术,并在某飞机模型上进行了验证.实验是在北京航空航天大学D4风洞中进行的.研究结果表明:该组合扰动控制技术能够实现对飞机大迎角非对称侧向力的有效主动控制.     

应用小干扰柱体控制角区马蹄涡结构的实验研究

本文提出一种简单的抑制和控制角区马蹄涡的被动控制方法.即在角区平板上游放置一个远小于主柱体的小干扰柱体,用其产生的弱马蹄涡来抑制和控制角区的马蹄涡结构.目的是使角区原来的马蹄涡结构由强变为弱、由大变为小、由多变为少、由非定常变为定常,以获得减小冲刷、抑制湍流、降低噪声、避免振动的工程效果.作者在风洞中采用烟线法和激光片光流动显示的方法开展研究,实验表明,在平板上游适当位置放置小干扰柱体的确可以有效抑制和控制角区马蹄涡结构.实验发现,当小干扰柱体放置在原角区马蹄涡生成区时,其抑制和控制效果最佳;当小干扰柱体放置在上游区或下游区时,控制效果不好.本文讨论了小干扰柱体控制角区马蹄涡的机理.此外,实验还研究了小干扰柱体相对尺度和截面形状对角区马蹄涡结构抑制和控制的影响.     

带涡襟翼翼型流场的数值模拟

以数值计算为手段,分析了带涡襟翼的翼型的流场特性,分别对迎角及扰流板偏角对翼型气动性能的影响做了分析。结果表明,在小迎角来流情况下,保持迎角不变,涡襟翼偏转角度越大,升力越小,阻力越大,呈现较好的线性关系。在大迎角情况下,绕翼型的流动发生分离,通过适当控制涡襟翼的偏转角度,能够有效的改善翼型的失速特性,从而达到流动控制的目的,迎角越大,涡襟翼所需偏转的角度越大。     

前体非对称涡高风速烟雾显示技术

通过在细长旋成体模型机身的不同位置设置出烟孔,使用新型发烟装置对模型注入雾化油 滴,采用激光片光法在常规风洞中高风速(60m/s)条件下得到清晰、完整而稳定的截面涡 结构图像. 并通过同状态标定法和等比例网格法发展了简易涡位定量测量技术,为前体非对 称涡Re数效应的研究提供了重要的技术平台.     

湍流边界层中马蹄形涡的形成方式

本文使用一种新的流动显示方法——激光片光运动法和几种实验技巧对湍流边界层中的马蹄形涡进行了观测,发现并描述了其形成的四种方式:二次不稳定式、组合式、变形式和突发式。对这四种马蹄形涡的形成及发展进行了研究和比较。实验结果表明,这些马蹄形涡在尺度、运动速度和变形上是有差别的。     

纳米示踪平面激光散射技术在激波复杂流场测量中的应用

在激波以及激波边界层相互作用这类含激波的复杂流场中,流场结构具有明显的三维特征.研究这类流场,采用纹影、阴影和干涉等传统流动显示技术空间分辨率较低,难以分辨流场的三维特性.基于纳米示踪的平面激光散射技术(nano-tracer planar laser scattering,NPLS),是作者近年来开发的一种新的研究超声速流场的测试与显示技术,可对超声速复杂三维流场进行高时空分辨率流动显示与测量.NPLS技术的特点使其成为测量激波复杂流场的有力手段.近年来,作者以NPLS技术为主要手段,对航空航天领域典型的激波复杂流场进行了试验研究,包括超声速弹头绕流、超声速混合层、超声速边界层,以及激波边界层相互作用流场,显示出NPLS技术在激波复杂流场精细测试与流动显示中优势.本文简要介绍NPLS技术在激波复杂流场测量中应用的研究进展.      

湍流边界层拟序结构的实验研究

20世纪60年代后, 先后从流动显示发现了快慢斑、猝发、上升流、下扫流和多种涡结构等湍流边界层的拟序结构. 它们对湍流边界层的摩阻、传热传质和湍动能的产生等特性有重要影响. 涡结构是上述拟序结构的核心, 它影响其它拟序结构的发展和演变. 发卡涡通常被认为是基本涡结构. 发卡涡等涡结构的再生, 是湍流边界层拟序结构能够自持续的必要的因素.壁面低速流上升产生猝发, 是湍流边界层湍能的主要来源; 条件采样是测量猝发频率和其它拟序结构出现频率的重要手段. 流动显示对湍流边界层拟序结构作了大量定性观察, 有许多减阻和增加传热率等应用性研究在此基础上发展起来. 80年代后, 出现了测量湍流边界层的瞬时流速矢量场的多热线法和PIV技术, 三维PIV技术可望将来为湍流边界层的实验研究带来重大进展. 本文评述了流动显示法、多热线法和PIV技术的优点和不足之处, 以及它们在对湍流边界层拟序结构的研究中的贡献.      

Vortex structure of separated flows on model wings at low freestream velocities

Flow past model wings is experimentally investigated in a subsonic wind tunnel at large angles of attack at which the laminar boundary layer separates near the leading edge of the wing (flow stall). The object of the study was the flow structure within the separation zone. The carbon-oil visualization of surface streamlines used in the experiments showed that in the separation zone there exist one or more pairs of large-scale vortices rotating in the wing plane. Certain general properties of the vortex structures in the separation zone are found to exist, whereas the flow patterns may differ depending on the model aspect ratio, the yaw angle, and other factors.     

翼反角对高压捕获翼构型亚声速气动特性影响分析研究

高压捕获翼新型气动布局在高超声速设计状态下具有较好的气动性能, 新升力面的引入使其在亚声速条件下也具有较大的升力, 但在亚声速下的稳定特性还有待研究. 基于高压捕获翼气动布局基本原理, 在机身-三角翼组合体上添加单支撑和捕获翼, 设计了一种参数化高压捕获翼概念构型. 以捕获翼和机体三角翼上/下反角为设计变量, 采用均匀试验设计、计算流体力学数值计算方法及Kriging代理模型方法, 研究了0° ~ 10°攻角状态下不同翼反角对高压捕获翼构型亚声速气动特性的影响, 重点分析了升阻特性、纵向和横航向稳定性的变化规律以及流场涡结构等. 结果表明, 小攻角状态下翼反角对升阻比的影响比大攻角更加显著, 捕获翼上反时, 升阻比略微增大, 下反则升阻比减小; 三角翼上反时, 升阻比减小, 下反则升阻比先略微增大后缓慢减小; 翼反角对纵向稳定性的总体影响较小, 捕获翼上反会稍微提高纵向稳定性, 而三角翼上反则会降低纵向稳定性; 捕获翼或三角翼上反都会增强横向稳定性, 下反则减弱横向稳定性, 但大攻角状态时, 三角翼上反角过大对提升横向稳定性作用有限; 捕获翼上反航向稳定性增强, 下反航向稳定性则减弱, 而三角翼下反对提升航向稳定性的整体效果比上反更加显著.      

飞机翼身干扰流动分离机理研究

采用CFD(计算流体力学)方法,结合风洞试验结果,对飞机翼身干扰流动分离机理进行了研究。分析了机翼根部马蹄涡、边角涡的内部结构和形成过程;研究了机身和机翼的压力分布及附面层,进而总结给出了飞机翼根流动分离机理;同时也验证了CFD方法对翼身干扰计算的准确性。研究结果表明:机翼前缘对气流的阻挡作用形成较高的逆压梯度,使气流产生回流形成马蹄涡,马蹄涡影响翼根表面附面层,使机翼机身附面层掺混失稳,向下游流动发展形成边角涡。     

The effect of phase angle and wing spacing on tandem flapping wings

In a tandem wing configuration, the hindwing often operates in the wake of the forewing and, hence, its performance is affected by the vortices shed by the forewing. Changes in the phase angle between the flapping motions of the fore and the hind wings, as well as the spacing between them, can affect the resulting vortex/wing and vortex/vortex interactions. This study uses 2D numerical simulations to investigate how these changes affect the leading dege vortexes (LEV) generated by the hindwing and the resulting effect on the lift and thrust coefficients as well as the efficiencies. The tandem wing configuration was simulated using an incompressible Navier-Stokes solver at a chord-based Reynolds number of 5 000. A harmonic single frequency sinusoidal oscillation consisting of a combined pitch and plunge motion was used for the flapping wing kinematics at a Strouhal number of 0.3. Four different spacings ranging from 0.1 chords to 1 chord were tested at three different phase angles, 0°, 90° and 180°. It was found that changes in the spacing and phase angle affected the timing of the interaction between the vortex shed from the forewing and the hindwing. Such an interaction affects the LEV formation on the hindwing and results in changes in aerodynamic force production and efficiencies of the hindwing. It is also observed that changing the phase angle has a similar effect as changing the spacing. The results further show that at different spacings the peak force generation occurs at different phase angles, as do the peak efficiencies.     

偏转矢量喷流对双三角翼前缘涡破裂影响的实验研究

本文应用染色液流动显示技术对后缘偏转喷流情况下76°/40°双三角翼前缘涡破裂位置的变化进行了观测,实验结果表明偏转喷流主要推迟与喷流方向相同一侧前缘涡的破裂,而使另一侧前缘涡破裂略有提前.随着喷流偏转角度的增大,喷流使两前缘涡破裂位置差逐渐增大.另外,随着模型攻角的增大,前缘涡涡核与双三角翼翼面的夹角逐渐增大,导致偏转喷流的作用逐渐减弱.     

Calculation of the flow of an ideal fluid past a wing of finite thickness

The problem of irrotational flow past a wing of finite thickness and finite span can be reduced by Green's formula to the solution of a system of Fredholm equations of the second kind on the surface of the wing [1]. The wake vortex sheet is represented by a free vortex surface. Besides panel methods (see, for example, [2]) there are also methods of approximate solution of this problem based on a preliminary discretization of the solution along the span of the wing in which the two-dimensional integral equations are reduced to a system of one-dimensional integral equations [1], for which numerical methods of solution have already been developed [3–6]. At the same time, a discretization is also realized for the wake vortex sheet along the span of the wing. In the present paper, this idea of numerical solution of the problem of irrotational flow past a wing of finite span is realized on the basis of an approximation of the unknown functions which is piecewise linear along the span. The wake vortex sheet is represented by vortex filaments [7] in the nonlinear problem. In the linear problem, the sheet is represented both by vortex filaments and by a vortex surface. Examples are given of an aerodynamic calculation for sweptback wings of finite thickness with a constriction, and the results of the calculation are also compared with experimental results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 124–131, October–December, 1981.     

Effect of apex strake incidence-angle on the vortex development and interaction of a double-delta wing

The vortex flow characteristics of a sharp-edged delta wing with an apex strake was investigated through the visualization and particle image velocimetry (PIV) measurement of the wing-leeward flow region, and the wing-surface pressure measurement. The wing model was a flat-plate, and 65°-sweep cropped-delta wing with sharp leading edges. The apex strake was also a flat-plate wing with a cropped-delta shape of 65°/90° sweep, and it can change its incidence angle. The flow Reynolds number was 2.2 × 10⁵ for the flow visualization and 8.2 × 10⁵ for the PIV and wing-surface pressure measurements. The physics of the vortex flow in the wing-leeward flow region and the suction-pressure distribution on the wing upper-surface were interrelated and analyzed. The effect of a positive (negative) strake incidence-angle was the upward movement of the strake and wing vortices away from (downward movement of the strake and wing vortices toward) the wing-upper surface and the delayed (enhanced) coiling interaction between them. This change of vortex flow characteristics projected directly on the suction pressure distribution on the wing upper-surface.     

Reynolds number effects on the vortical-flow structure generated by a double-delta wing

An experimental investigation of the high-incidence vortical flowfield over a 76/40° double-delta wing model with sharp leading edges was conducted in the Naval Postgraduate School water tunnel facility at three nominal flow Reynolds numbers of 15000, 45000, and 75000 (based on centerline chord). Extensive flow visualization studies were performed with the dye-injection technique, followed by laser Doppler velocity measurements. The primary objective of this investigation was the determination of the influence of Reynolds number on vortex interactions/trajectories, and breakdown. It was found that there is a significant influence of Reynolds number. Specifically, with the increase of flow Reynolds number the strake and wing vortex trajectories tend to move outboards and closer to the model surface, and the vortex breakdown location moves forwards toward the apex of the model. The intertwining or coiling-up feature of the vortex interaction phenomenon becomes less dominant and disappears altogether at high Reynolds numbers. These trends in the vortex interaction and bursting data are found to be in good agreement with previous wind tunnel data. Received: 26 March 1998/Accepted: 2 February 1999     

A DPIV study of the trailing vortex elements from the blades of a horizontal axis wind turbine in yaw

Digital particle image velocimetry (DPIV) has been used in a wind tunnel study to measure the velocity field of the trailing vortices from the blades of a horizontal axis wind-turbine (HWAT) in yaw. The creation of the trailing vortex circulation is shown to vary as a function of the phase angle of the rotor and the angle of yaw between the wind and turbine rotor. The strength of the convecting vorticity was also shown to vary with time. The initial formation of the vortex is shown determined by the flow expansion angle while in yawed flows the developing vortex is then influenced by the vortex sheet shed from the inboard blade trailing edge. This interaction is shown to significantly affect the roll up of the tip vortex. Received: 17 December 1997/Accepted: 16 June 1999