柔性扑翼的气动特性研究

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

Flow field characteristics study of a flapping airfoil using computational fluid dynamics

The flow field of a flapping airfoil in Low Reynolds Number (LRN) flow regime is associated with complex nonlinear vortex shedding and viscous phenomena. The respective fluid dynamics of such a flow is investigated here through Computational Fluid Dynamics (CFD) based on the Finite Volume Method (FVM). The governing equations are the unsteady, incompressible two-dimensional Navier-Stokes (N-S) equations. The airfoil is a thin ellipsoidal geometry performing a modified figure-of-eight-like flapping pattern. The flow field and vortical patterns around the airfoil are examined in detail, and the effects of several unsteady flow and system parameters on the flow characteristics are explored. The investigated parameters are the amplitude of pitching oscillations, phase angle between pitching and plunging motions, mean angle of attack, Reynolds number (Re), Strouhal number (St) based on the translational amplitudes of oscillations, and the pitching axis location (x/c). It is shown that these parameters change the instantaneous force coefficients quantitatively and qualitatively. It is also observed that the strength, interaction, and convection of the vortical structures surrounding the airfoil are significantly affected by the variations of these parameters.     

N-S方程数值研究翼型对微型扑翼气动特性的影响

首先基于嵌套网格发展了一套适用于三维扑翼研究的非定常雷诺平均Navier-Stokes(RANS)方程数值模拟方法.为了解决微型扑翼在低马赫数下的收敛问题,使用了预处理方法,湍流模型为BL模型.在该方法的基础上,保持状态参数和扑翼表面形状一定的情况下,分别研究了一系列不同厚度、不同弯度的翼型对于微型扑翼气动特性的影响....     

Power extraction efficiency improvement of a fully-activated flapping foil: With the help of an auxiliary rotating foil

The power extraction efficiency improvement of a fully-activated flapping foil with the help of an auxiliary rotating foil is numerically examined in this work. A NACA0015 airfoil is placed in a two-dimensional laminar flow and synchronously performs the imposed pitching and plunging motions. An auxiliary smaller foil, which rotates about its center, is arranged below the flapping foil. As a consequence, the vortex interaction between the flapping foil and the rotating foil occurs. At a Reynolds number of 1100 and the position of the pitching axis at one-third chord, the effects of the distance between the flapping foil and the auxiliary foil, the phase difference between the rotating motion and the flapping motion as well as the frequency of flapping motion on the power extraction performance are systematically investigated. It is found that compared to the single flapping foil, the efficiency of power extraction for the flapping foil with an auxiliary device can be improved. Based on the numerical analysis, it is indicated that the enhanced plunging contribution, which is caused by the increased lift force owing to the vortex interaction, directly helps the efficiency improvement.     

Flow Characteristics of Flapping Motion of a Plane Water Jet Impinging onto Free Surface

A submerged turbulent plane jet in shallow water impinging vertically onto the free surface will produce a large-scale flapping motion when the jet exit velocity is larger than a critical one. The flapping phenomenon is verified in this paper through a large eddy simulation where the free surface is modeled by volume of fluid approach. The quantitative results for flapping jet are found to be in good agreement with available experimental data in terms of mean velocity, flapping-induced velocity and turbulence intensity. Results show that the flapping motion is a new flow pattern with characteristic flapping frequency for submerged turbulent plane jets, the mean centerline velocity decay is considerably faster than that of the stable impinging jet without flapping motion, and the flapping-induced velocities are as important as the turbulent fluctuations.     

Aerodynamic force and flow structures of two airfoils in flapping motions

Aerodynamic force and flow structures of two airfoils in a tandem configuration in flapping motions are studied, by solving the Navier-Stokes equations in moving overset grids. Three typical phase differences between the fore- and aftairfoil flapping cycles are considered. It is shown that: (1) in the case of no interaction (single airfoil), the time average of the vertical force coefficient over the downstroke is 2.74, which is about 3 times as large as the maximum steady-state lift coefficient of a dragonfly wing; the time average of the horizontal force coefficient is 1.97, which is also large. The reasons for the large force coefficients are the acceleration at the beginning of a stroke, the delayed stall and the “pitching-up” motion near the end of the stroke. (2) In the cases of two-airfoils, the time-variations of the force and moment coefficients on each airfoil are broadly similar to that of the single airfoil in that the vertical force is mainly produced in downstroke and the horizontal force in upstroke, but very large differences exist due to the interaction. (3) For in-phase stroking, the major differences caused by the interaction are that the vertical force on FA in downstroke is increased and the horizontal force on FA in upstroke decreased. As a result, the magnitude of the resultant force is almost unchanged but it inclines less forward. (4) For counter stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased. As a result, the magnitude of the resultant force is decreased by about 20 percent but its direction is almost unchanged. (5) For 90°-phase-difference stroking, the major differences are that the vertical force on AA in downstroke and the horizontal force on FA in upstroke are decreased greatly and the horizontal force on AA in upstroke increased. As a result, the magnitude of the resultant force is decreased by about 28% and it inclines more forward. (6) Among the three cases of phase angles, inphase flapping produces the largest vertical force (also the largest resultant force); the 90°-phase-difference flapping results in the largest horizontal force, but the smallest resultant force.     

昆虫拍翼方式的非定常流动物理再探讨

基于提出的理论模化方法来探讨昆虫拍翼方式的非定常流动物理. 以悬停飞行为 例,通过对拍翼运动的分析,不仅解释了昆虫利用高频拍翼的方式为何能够克服低雷诺数带 来的气动局限性(St \gg 1/Re),而且还指出高升力产生和调节的3个流动 控制因素：(1) 由于拍翼的变速运动即时引起了流体动力响应,这种附加惯性效应 可产生瞬时的高升力; (2) 保持前缘涡不脱离翼面有助于减少升力的下降; (3) 增大后缘涡的强度并加速其脱离后缘能够有效地提高升力.     

HUNTER: A conformational search program for acyclic to polycyclic molecules with special emphasis on stereochemistry

A new conformational search program, HUNTER, connected with the force fields MMP2 and MM3(92) is presented. The program accepts all types of molecules with most different substructures, considers stereochemical facts, and covers conformational space efficiently and completely. The most important facilities are an automated analysis of the stereochemistry including topographical facts, a separate perturbation of the acyclic and cyclic parts of the molecule using modified corner flapping, and an incremental rotation around single bonds with fixed flap and rotation angles, respectively; an exclusion of high energy structures by simulated annealing; the choice of the conformer lowest in energy, which is new as an initial structure for the next sampling run; and the use of a reduced set of dihedral angles to define a conformation. A specifically devised graphic interface, SERVANT, is used to feed in and control all informations necessary for a program run and to visualize the results. Most of the parameters are user-defined and thereby allow a flexible search, including a search for the most stable diastereomer. The efficiency of the different parameter sets was tested in calculation with cycloundecane ( 12 ), (Z)-oct-3-ene ( 13 ), and sipholenol-A monoacetate ( 14 ). The best performance regarding the number of different low-energy conformers was achieved with 60° ( 14 ) and 90° flaps ( 12 ), respectively, including substituent correction for the cyclic parts, and with 105° ( 14 ) and 120° rotations ( 13 ), respectively, for the acyclic parts. In comparison to the stochastic search routine implemented in MM3(92), HUNTER performed two ( 12 ) to six ( 14 ) times better. © 1997 John Wiley & Sons, Inc. J Comput Chem 18: 1264–1281, 1997     

Stream function–potential function coordinates for aeroacoustics and unsteady aerodynamics

‘Stream function as a coordinate approach’ (SFC) combined with compact high-order finite difference schemes has been developed and applied to aeroacoustics and unsteady aerodynamics problems. Straightforward implementation of SFC creates coarse grids at the vicinity of stagnation points that smears high-order numerical computations. Grid clustering is employed to resolve coarse grid near stagnations points. The agreement between numerical results and particle image velocimetry (PIV) measurements for flapping airfoil shows the robustness of the current approach for performing high-order computations.     

Spatio-temporal analysis of the turbulent flow in a ribbed channel

The turbulent flow in a channel with transverse ribs over one wall is studied experimentally. The height of the obstacles is about one tenth of the channel height, and the spacing is 10 times their height. The Reynolds number based on the channel hydraulic diameter and bulk flow velocity is 15,000. Velocity fields are obtained with high spatial and temporal resolution along the streamwise/wall-normal plane by means of time-resolved particle image velocimetry. Beside mean velocity and Reynolds stresses, the flow is investigated through two-point correlations, distributions of spanwise-swirling events, space–time velocity diagrams and power spectral density. Although the separated flow reattaches before the following obstacle is approached, a strong rib-to-rib interaction occurs. Spanwise vortices, 0.2 rib heights in size, are generated in the free shear layer, travel across the whole pitch, and may impact on the next rib. The large scale motions triggered by the separation grow in size until they reach the following obstacle. Flapping of the separated shear layer is observed at frequencies consistent with previous studies, causing the instantaneous reattachment point to fluctuate. The flapping initiates at the downstream edge of the obstacle tip, rather than at the upstream edge where the instantaneous separation occurs.