基于雨燕翅膀的仿生三角翼气动特性计算研究

针对低雷诺数微型飞行器的气动布局, 设计出类似雨燕翅膀的一组具有不同前缘钝度的中等后掠($\varLambda =50^{\circ}$)仿生三角翼. 为了定量对比研究三角翼后缘收缩产生的气动效应, 设计了一组具有同等后掠的普通三角翼. 为了深入研究仿生三角翼布局的前缘涡演化特性以及总体气动特性, 采用数值模拟方法详细地探索了低雷诺数($Re=1.58\times 10^{4})$流动条件下前缘涡涡流结构和气动力随迎角的变化规律. 分析结果表明, 前缘钝度和后缘收缩对仿生三角翼前缘涡的涡流强度和涡破裂位置有显著影响. 相对于钝前缘来说, 尖前缘使仿生三角翼上下表面的压力差增大, 涡流强度也更大, 增升作用也更显著. 相对于普通三角翼构型, 仿生三角翼的前缘斜切使其阻力更大, 但后缘的收缩使涡破裂位置固定在此位置, 因此整个上翼面保持低压, 总的升力更大. 由于小迎角时升力增大更明显, 因此仿生三角翼的气动效率在小迎角时明显大于普通三角翼. 这些结论对于揭示鸟类的飞行机理以及未来微型仿生飞行器的气动布局设计具有重要的研究价值.      

HEAT-PIPE-COOLED THERMAL PROTECTION FOR THE LEADING EDGE OF THE WING

针对高超声速飞行器飞行时翼前缘存在着严重的气动加热问题，为了保证翼前缘的尖锐外形，提出疏导式热防护结构，利用内置高温热管结构为翼前缘提供热防护。采用数值模拟和电弧风洞试验的方法对翼前缘疏导式结构进行了分析，得到翼前缘内置高温热管具有的防热效果。数值模拟结果表明在一定热环境条件下，翼前缘驻点温度下降了304K，尾部最低温度升高了130K，实现了热流从高温区到低温区的疏导，减弱了翼前缘的热载荷，强化了翼前缘的热防护能力。通过电弧风洞试验可以获得相同的热防护结果，并且在一定飞行条件下高温热管可以自适应启动，验证了数值模拟方法的准确性以及翼前缘内置高温热管疏导式热防护结构的可行性。     

Wings of minimum drag for given lift

The variational problem of determining the optimal shape (camber and twist) of the midsurface of a wing having minimum wavedrag is examined in the linear formulation. It is shown that for wings with supersonic leading edge and straight trailing edge, whose shape is given in the form of a double polynomial, the over-all aerodynamic characteristics can be simply expressed in terms of the equation for the leading edge of the wing. This makes it possible not only to solve the variational problem by the Ritz method and obtain the minimum wave drag [1] but also to find the optimal shape of the wing. As examples we consider delta and double-delta wings.     

扑翼飞行器带翼刀机翼气动特性实验研究

通过进行微型扑翼飞行器低速风洞试验,研究了带弯度机翼下翼面翼刀对扑翼飞行器升阻特性的影响。文中进行了带翼刀机翼和不带翼刀机翼在不同迎角下的风洞吹风试验。试验结果表明,带翼刀机翼升力系数大于不带翼刀机翼升力系数,从而证明了翼刀可以阻止机翼下表面气流展向流动,起到增加机翼升力的作用。当扑翼在小迎角飞行时,带翼刀机翼可以有效地提高扑翼的气动效率,改善扑翼的飞行性能。研究结果可为带翼刀机翼在扑翼飞行器上的应用提供技术支持。     

三角翼的双襟翼控涡作用的数值模拟研究

对装有“前端襟翼”和“前缘襟翼”的７４°后掠三角翼的不可压缩流场作了数值模拟,以研究襟翼的旋涡控制作用．数值模拟是用拟压缩性方法求解一般曲线坐标系下的三维不可压缩Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,时间离散用向后Ｅｕｌｅｒ差分,空间无粘项的离散用二阶迎风ＴＶＤ格式,所得的离散方程用对角化形式的近似隐式因子分解格式求解．湍流模型用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数模式．计算了三种平面形状的机翼在迎角范围为１０°～５０°的绕流和气动特性．计算和实验的比较表明,襟翼向下偏转可以推迟旋涡破裂,且对提高机翼的减阻能力、升阻比和改善失速前后的气动特性有明显效果,双襟翼具有更佳的控涡效果．     

Nondimensional frequency scaling of aerodynamically-tensioned membranes

Membrane wings have applications that involve low Reynolds number flyers such as micro air vehicles. The time-averaged and time-dependent deformations of the membrane affect the aerodynamic characteristics of the wing, primarily in the region beyond the maximum aerodynamic efficiency of the wing. This paper investigates an appropriate nondimensional vibration frequency scaling of a spanwise tensioned membrane with free (unattached) leading and trailing edges at low Reynolds numbers relative to nondimensional aeroelastic parameters. Silicone rubber membranes with varying spanwise pre-tension, aerodynamic tension (due to wing angle-of-attack and flow dynamic pressure), modulus of elasticity, span, and thickness are studied. Experimental results are compared to a proposed scaling that simplifies the aerodynamic loading as a uniform pressure distribution acting on the membrane. Data is further compared and discussed relative to previous published results of membrane wings with finite wing spans (three-dimensional flow) and fixed (rigid) leading edges.     

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

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

柔性扑翼的气动特性研究

以往扑翼的气动力计算研究都很少考虑扑翼的柔性,而在鸟的扑翼动作中,在外加气动力和鸟自身的扑动力作用下,扑翼的柔性变形相当大。本文在原有匀速刚性模型的基础上,提出考虑了扑翼扑动速率变化和形状变化的扑翼分析模型,使之更接近鸟翼柔性扑动真实情况。通过计算分析气动特性发现,控制适当的话,柔性变形能大大改善扑翼的气动性能。本文通过模拟鸟扑翼的柔性运动,计算了时柔性扑翼气动力以及平均升力系数和平均推力系数随着扑动角、倾斜角等参数变化的情况,从而从气动的角度解释了为什么鸟在不同的飞行阶段扑翼规律各不相同,并为柔性扑翼飞行器的设计提供了理论依据。     

Flow past a delta wing with a sinusoidal leading edge: near-surface topology and flow structure

The near-surface flow structure and topology on a delta wing of low sweep angle, having sinusoidal leading edges of varying amplitude and wavelength, are investigated using a stereoscopic technique of high-image-density particle image velocimetry at a Reynolds number of 15,000. Identification of critical points, in conjunction with surface-normal vorticity and velocity, provides a basis for determining the effectiveness of a given leading edge. At high angle of attack, where large-scale three-dimensional separation occurs from the wing with a straight leading edge, an amplitude of the leading-edge protuberance as small as one-half of one percent of the chord of the wing can substantially alter the near-surface topology. When the amplitude reaches a value of four percent of the chord, it is possible to completely eradicate the negative focus of large-scale, three-dimensional separation, in favor of a positive focus of attachment. Moreover, alteration of the near-surface topology is most effective when the ratio of the wavelength to amplitude of the sinusoidal leading edge is maintained at a small value.     

ACTIVE WING FLUTTER SUPPRESSION USING A TRAILING EDGE FLAP

The aeroservoelastic behaviour of a thin rectangular wing with a controllable trailing edge flap is investigated. A rather high aspect ratio motivates a numerical model based on linear beam theory for the structural dynamics and strip theory for the unsteady aerodynamic loads. Experimental flutter testing shows good agreement with the numerical stability analysis, and the impact of the trailing edge flap on the dynamics is verified by open-loop testing. The problem of stabilizing the wing utilizing the trailing edge flap is posed, and the design of a fixed-structure feedback controller is performed using numerical optimization. The problem of maximizing closed-loop modal damping with constraints on actuator performance is solved for a sequence of flow speeds and the obtained controller is synthesized using gain scheduling. The fairly large predicted increase in critical speed is experimentally verified with satisfactory accuracy.     

On the aerodynamic characteristics of hovering rigid and flexible hawkmoth-like wings

Insect wings are subjected to fluid, inertia and gravitational forces during flapping flight. Owing to their limited rigidity, they bent under the influence of these forces. Numerical study by Hamamoto et al. (Adv Robot 21(1–2):1–21, 2007) showed that a flexible wing is able to generate almost as much lift as a rigid wing during flapping. In this paper, we take a closer look at the relationship between wing flexibility (or stiffness) and aerodynamic force generation in flapping hovering flight. The experimental study was conducted in two stages. The first stage consisted of detailed force measurement and flow visualization of a rigid hawkmoth-like wing undergoing hovering hawkmoth flapping motion and simple harmonic flapping motion, with the aim of establishing a benchmark database for the second stage, which involved hawkmoth-like wing of different flexibility performing the same flapping motions. Hawkmoth motion was conducted at Re = 7,254 and reduced frequency of 0.26, while simple harmonic flapping motion at Re = 7,800 and 11,700, and reduced frequency of 0.25. Results show that aerodynamic force generation on the rigid wing is governed primarily by the combined effect of wing acceleration and leading edge vortex generated on the upper surface of the wing, while the remnants of the wake vortices generated from the previous stroke play only a minor role. Our results from the flexible wing study, while generally supportive of the finding by Hamamoto et al. (Adv Robot 21(1–2):1–21, 2007), also reveal the existence of a critical stiffness constant, below which lift coefficient deteriorates significantly. This finding suggests that although using flexible wing in micro air vehicle application may be beneficial in term of lightweight, too much flexibility can lead to deterioration in flapping performance in terms of aerodynamic force generation. The results further show that wings with stiffness constant above the critical value can deliver mean lift coefficient almost the same as a rigid wing when executing hawkmoth motion, but lower than the rigid wing when performing a simple harmonic motion. In all cases studied (7,800 ≤ Re ≤ 11,700), the Reynolds number does not alter the force generation significantly.     

Aerodynamic mechanism of forces generated by twisting model-wing in bat flapping flight

The aerodynamic mechanism of the bat wing membrane Mong the lateral border of its body is studied. The twist-morphing that alters the angle of attack （AOA） along the span-wise direction is observed widely during bat flapping flight. An assumption is made that the linearly distributed AOA is along the span-wise direction. The plate with the aspect ratio of 3 is used to model a bat wing. A three-dimensional （3D） unsteady panel method is used to predict the aerodynamic forces generated by the flapping plate with leading edge separation. It is found that, relative to the rigid wing flapping, twisting motion can increase the averaged lift by as much as 25% and produce thrust instead of drag. Furthermore, the aerodynamic forces （lift/drag） generated by a twisting plate-wing are similar to those of a pitching rigid-wing, meaning that the twisting in bat flight has the same function as the supination/pronation motion in insect flight.     

Discontinuous solutions of an interacting three-dimensional hypersonic boundary layer

The paper is a mathematical study of the three-dimensional flow of viscous gas in a hypersonic boundary layer that develops along a flat wing whose leading edge has a step shape. The flow interacts with a flap on the wing set at a small angle. A linear solution to the problem is constructed under the assumption that the deflection angle of the flap is small and the difference between the length of the plates is of order unity. It is shown that an important part in the formation of the flow near and behind the flap may be played by the change in the pressure along the span of the wing due to the step shape of the leading edge. It is significant that although the pressure and displacement thickness are continuous functions of the transverse coordinate, the longitudinal and transverse components of the friction force have discontinuities.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 19–26, March–April, 1991.I thank V. V. Sychev and A. I. Ruban for suggesting the problem, for valuable advice, and assistance.     

A novel MAV with treadmill motion of wing

The Magnus effect is well known phenomena for producing high lift values from spinning symmetrical geometries such as cylinders, spheres, or disks. But, the Magnus force may also be produced by treadmill motion of aerodynamic bodies. To accomplish this, the skin of aerodynamic bodies may circulate with a constant circumferential speed. Here, a novel wing with treadmill motion of skin is introduced which may generate lift at zero air speeds. The new wing may lead to micro aerial vehicle configurations for vertical take-off or landing. To prove the concept, the NACA0015 aerofoil section with circulating skin is computationally investigated. Two cases of stationary air and moving air are studied. It is observed that lift can be generated in stationary air although drag force is also high. For moving air, the lift and drag forces may be adopted between the incidence angles 20° to 25° where lift can posses high values and drag can remain moderate.     

Combined application of pressure sensitive paint and particle image velocimetry to the flow above a delta wing

The pressure sensitive paint (PSP) and particle image velocimetry (PIV) techniques are applied to the flow around a delta wing. These investigations are a contribution to the International Vortex Flow Experiment 2 (VFE-2). The delta wing is equipped with rounded as well as with sharp leading edges. To demonstrate the advantages of a combined application of PSP and PIV for flow investigation the results of a specific case are discussed, where an inner and outer primary vortex above the delta wing with rounded leading edges develop.     

高速飞行器中空翼结构高温热振动特性试验研究

远程高速飞行器飞行速度快, 滞空时间长, 飞行过程中翼、舵等结构会出现长时间的剧烈振动, 由气动加热产生的高温还会使飞行器材料和结构的弹性性能发生变化, 从而引起翼、舵等结构振动特性的改变.因此获得高温与振动复合环境下的远程高速飞行器翼、舵等结构的振动特性参数对于高速飞行器的安全设计具有非常重要的意义.将高温热环境试验系统与振动试验系统相结合, 在对中空翼面结构进行振动激励的同时使用红外辐射加热方式对翼面结构生成可控的热环境, 并通过自行设计的耐高温引伸装置将中空翼结构的振动信号传递到非高温区进行数据采集与分析的方式, 实现了高达800℃~900℃的力热复合环境下的翼结构固有频率、模态等振动特性参数的试验测试, 其试验结果为远程高速飞行器中空翼结构在高温振动环境下的动特性分析和安全可靠性设计提供了重要依据.     

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

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

The PELskin project-part V: towards the control of the flow around aerofoils at high angle of attack using a self-activated deployable flap

During the flight of birds, it is often possible to notice that some of the primaries and covert feathers on the upper side of the wing pop-up under critical flight conditions, such as the landing approach or when stalking their prey (see Fig. 1) . It is often conjectured that the feathers pop up plays an aerodynamic role by limiting the spread of flow separation . A combined experimental and numerical study was conducted to shed some light on the physical mechanism determining the feathers self actuation and their effective role in controlling the flow field in nominally stalled conditions. In particular, we have considered a NACA0020 aerofoil, equipped with a flexible flap at low chord Reynolds numbers. A parametric study has been conducted on the effects of the length, natural frequency, and position of the flap. A configuration with a single flap hinged on the suction side at 70 % of the chord size c (from the leading edge), with a length of \(L=0.2c\) matching the shedding frequency of vortices at stall condition has been found to be optimum in delivering maximum aerodynamic efficiency and lift gains. Flow evolution both during a ramp-up motion (incidence angle from \(\alpha _0=0\) to \(\alpha _s=20^\circ\) with a reduced frequency of \(k= 0.12\, U_{\infty }/c\), \(U_{\infty }\) being the free stream velocity magnitude), and at a static stalled condition (\(\alpha =20^\circ\)) were analysed with and without the flap. A significant increase of the mean lift after a ramp-up manoeuvre is observed in presence of the flap. Stall dynamics (i.e., lift overshoot and oscillations) are altered and the simulations reveal a periodic re-generation cycle composed of a leading edge vortex that lift the flap during his passage, and an ejection generated by the relaxing of the flap in its equilibrium position. The flap movement in turns avoid the interaction between leading and trailing edge vortices when lift up and push the trailing edge vortex downstream when relaxing back. This cyclic behaviour is clearly shown by the periodic variation of the lift about the average value, and also from the periodic motion of the flap. A comparison with the experiments shows a similar but somewhat higher non-dimensional frequency of the flap oscillation. By assuming that the cycle frequency scales inversely with the boundary layer thickness, one can explain the higher frequencies observed in the experiments which were run at a Reynolds number about one order of magnitude higher than in the simulations. In addition, in experiments the periodic re-generation cycle decays after 3–4 periods ultimately leading to the full stall of the aerofoil. In contrast, the 2D simulations show that the cycle can become self-sustained without any decay when the flap parameters are accurately tuned.     

Numerical simulation of flapping‐wing insect hovering flight at unsteady flow

A computational fluid dynamics (CFD) analysis was conducted to study the unsteady aerodynamics of a virtual flying bumblebee during hovering flight. The integrated geometry of bumblebee was established to define the shape of a three‐dimensional virtual bumblebee model with beating its wings, accurately mimicking the three‐dimensional movements of wings during hovering flight. The kinematics data of wings documented from the measurement to the bumblebee in normal hovering flight aided by the high‐speed video. The Navier–Stokes equations are solved numerically. The solution provides the flow and pressure fields, from which the aerodynamic forces and vorticity wake structure are obtained. Insights into the unsteady aerodynamic force generation process are gained from the force and flow‐structure information. The CFD analysis has established an overall understanding of the viscous and unsteady flow around the virtual flying bumblebee and of the time course of instantaneous force production, which reveals that hovering flight is dominated by the unsteady aerodynamics of both the instantaneous dynamics and also the past history of the wing. A coherent leading‐edge vortex with axial flow and the attached wingtip vortex and trailing edge vortex were detected. The leading edge vortex, wing tip vortex and trailing edge vortex, which caused by the pressure difference between the upper and the lower surface of wings. The axial flow, which include the spanwise flow and chordwise flow, is derived from the spanwise pressure gradient and chordwise pressure gradient, will stabilize the vortex and gives it a characteristic spiral conical shape. Copyright © 2006 John Wiley & Sons, Ltd.