Numerical investigations into the asymmetric effects on the aerodynamic response of a pitching airfoil

The effects of asymmetric sinusoidal motion on pitching airfoil aerodynamics were studied by numerical simulations for 2-D flow around a NACA0012 airfoil at Re=1.35×10⁵. Various unsteady parameters (amplitude of oscillation, d; reduced frequency, k) were applied to investigate the effect of asymmetry parameter S on the instantaneous force coefficients and flow patterns. The results reveal that S has a noticeable effect on the aerodynamic performance, as it affects the instantaneous force coefficient, maximum lift and drag coefficient, hysteresis loops and the flow structures.     

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.     

Investigation of flow field of clap and fling motion using immersed boundary coupled lattice Boltzmann method

This paper deals with the investigation of flow field due to clap and fling mechanism using immersed boundary coupled with lattice Boltzmann method. The lattice Boltzmann method (LBM), an alternative to Navier–Stokes solver, is used because of its simplicity and computational efficiency in solving complex moving boundary problems. Benchmark problems are simulated to validate the code, which is then used for simulating flow over two elliptic wing of aspect ratio 5 performing clap and fling flapping motion for different flow parameters such as Reynolds number (Re=75, 100, 150), advance ratio (J=10E−3.10E−2,0.2) and frequency (f=0.05 Hz, 0.25 Hz). Numerical simulation is able to capture typical low Reynolds number unsteady phenomena such as, ׳wake vortex wing interaction׳, ׳Kramer effect׳ and ׳delayed stall׳. The results are both qualitatively and quantitatively consistent with experimental observation. The parametric study involving different combinations of Re, f and J depict distinctly different aerodynamic performances providing physical insights into the flow physics. It is observed that a combination of low f, low J and high Re flow results in better aerodynamic performance. Pronounced lift enhancement via leading edge vortices are obtained in unsteady regime (J<1) compared to quasi-steady regime (J>1). The role of leading edge vortices in enhancing lift are investigated by studying the size and strength of these vortices for different flow conditions. For a given Re, the magnitude of maximum lift coefficient decreases with increasing f irrespective of the value of J; while the same is enhanced with the increasing Re.     

Influence of heat transfer on the aerodynamic performance of a plunging and pitching NACA0012 airfoil at low Reynolds numbers

The unsteady low Reynolds number aerodynamics phenomena around flapping wings are addressed in several investigations. Elsewhere, airfoils at higher Mach numbers and Reynolds numbers have been treated quite comprehensively in the literature. It is duly noted that the influence of heat transfer phenomena on the aerodynamic performance of flapping wings configurations is not well studied. The objective of the present study is to investigate the effect of heat transfer upon the aerodynamic performance of a pitching and plunging NACA0012 airfoil in the low Reynolds number flow regime with particular emphasis upon the airfoil's lift and drag coefficients. The compressible Navier–Stokes equations are solved using a finite volume method. To consider the variation of fluid properties with temperature, the values of dynamic viscosity and thermal diffusivity are evaluated with Sutherland's formula and the Eucken model, respectively. Instantaneous and mean lift and drag coefficients are calculated for several temperature differences between the airfoil surface and freestream within the range 0–100 K. Simulations are performed for a prescribed airfoil motion schedule and flow parameters. It is learnt that the aerodynamic performance in terms of the lift C_L and drag C_D behavior is strongly dependent upon the heat transfer rate from the airfoil to the flow field. In the plunging case, the mean value of C_D tends to increase, whereas the amplitude of C_L tends to decrease with increasing temperature difference. In the pitching case, on the other hand, the mean value and the amplitude of both C_D and C_L decrease. A spectral analysis of C_D and C_L in the pitching case shows that the amplitudes of both C_D and C_L decrease with increasing surface temperature, whereas the harmonic frequencies are not affected.     

Numerical study of large amplitude,nonsinusoidal motion and camber effects on pitching airfoil propulsion

The effects of large amplitude and nonsinusoidal motion on pitching airfoil aerodynamics for thrust generation were numerically studied with a 2-D NACA0012 airfoil used, and various 2-D NACA asymmetric airfoils were applied for camber effect study. The large amplitude effect study has been undertaken over a wide range of reduced frequency k (from 6 to 18) and pitching amplitude θ (from 5° to 30°) at Re=1.35×10⁴ with sinusoidal pitching profile used. It is shown that the large pitching amplitude results in much more thrust generated than that at low pitching amplitude and the increase of thrust with amplitude becomes slow when the amplitude reaches some degree. However, the propulsive efficiency noticeably decreases with the increase of θ at a fixed k.An adjustable parameter K was employed to realize various nonsinusoidal motions and the effect of nonsinusoidal motion was investigated with various unsteady parameters (θ, k) applied. The results reveal that nonsinusoidal motion has a noticeable effect on the aerodynamic performance, as it affects the instantaneous force coefficients, maximum thrust coefficients and flow structures. An increase in K results in a better thrust generation performance at fixed θ and k, especially for K>0. It is also shown that the larger K noticeably influences the wake pattern and induces a stronger reverse von Karman vortex street in the wake, which in turn leads to the increased thrust. The camber study was performed on various 2-D NACA airfoils with different cambers and camber locations undergoing sinusoidal pitching motion at θ=5° and Re=1.35×10⁴. It is found that varying camber offers little improvement in thrust generation performance.     

Numerical simulation of Reynolds number effects on velocity shear flow around a circular cylinder

Three-dimensional Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) are performed to investigate the shear effects on flow around a circular cylinder at Reynolds numbers of Re=60–1000. The shear parameter, β, which is based on the velocity gradient, cylinder diameter and upstream mean velocity at the center plane of the cylinder, varies from 0 to 0.30. Variations of Strouhal number, drag and lift coefficients, and unsteady wake structures with shear parameter are studied, along with their dependence on Reynolds number. The presented simulation provides detailed information for the flow field around a circular cylinder in shear flow. This study shows that the Strouhal number exhibits no significant variation with shear parameter. The stagnation point moves to the high-velocity side almost linearly with shear parameter, and this result mainly influences the aerodynamic forces acting on a circular cylinder in shear flow. Both the Reynolds number and shear parameter influence the movement of the stagnation point and separation point. Mode A wake instability is suppressed into parallel vortex shedding mode at a certain shear parameter. The lift force increases with increasing shear parameter and acts from the high-velocity side to the low-velocity side. In addition, a simple method to estimate the lift force coefficient in shear flow is provided.     

Dynamic stall behavior from unsteady force measurements

A direct force measurement technique employing piezoelectric load cells is used to experimentally investigate a two-dimensional airfoil (NACA 0012) undergoing dynamic stall. The load cells are installed at each end of the airfoil and give the force response in two directions in the plane normal to the airfoil axis during oscillations. Experiments are carried out at a Reynolds number based on the airfoil chord equal to 7.7×10⁴, and at four reduced frequencies, k=0.005, 0.01, 0.02, and 0.04. Phase-averaged lift of the airfoil undergoing dynamic stall is presented. It is observed that hysteresis loops of the lift occur both when the airfoil is pitched to exceed its static stall limit and when it is still within its static stall limit, and they grow in size with increasing k at the same pitching mean angle of attack and pitching amplitude. Both the lift and the drag induced by the pitching motion are further analyzed using the methods of higher order correlation analysis and continuous wavelet transforms to undercover their nonlinear and nonstationary features, in addition to classical FFT-based spectral analysis. The results are quantitatively illustrated by an energy partition analysis. It is found that the unsteady lift and drag show opposite trends when the airfoil undergoes transition from the pre-stall regime to the full-stall regime. The degree of nonlinearity of the lift increases, and the lift show a nonstationary feature in the light-stall regime, while the nonlinearity of the drag decreases, and the drag shows nonstationary feature in both the light-stall and the full-stall regimes. Furthermore, the lift and the drag have significant nonlinear interactions as shown by the correlation analysis in the light-stall regime.     

Experimental benchmark of a free plunging wing with imposed flap oscillations

The validation of fluid–structure interaction solvers is difficult since there is a lack of experimental data. Therefore, in this work an aeroelastic experiment is presented. The focus is on the temporal coupling between fluid and structure dynamics. Issues in the spatial coupling are eliminated by using a rigid wing. The wing, with a harmonically actuated 0.2c trailing edge flap, has a degree of freedom in the plunge (vertical) direction. The wing has a chord of 0.5 m and is suspended with springs. The wing motion is constrained by a vertical rail system.For simplicity attached flow is desired and therefore the set angle of attack is α=0°. The Reynolds number is approximately Re=700 000 and the flap deflects over a range of about ±2°. The damped natural frequency of the structure expressed as a reduced frequency is about k=0.194 and measurements are performed for reduced flap frequencies ranging from k=0.1 to k=0.3. Displacements and time dependent aerodynamic forces are measured and for k=0.198 2-D PIV measurements are performed. The planar PIV measurements are used to intrinsically determine the unsteady loads using Noca׳s method.As expected the aeroelastic problem shows similarities with a viscously damped mass–damper–spring, meaning the maximum excursion of the wing is found near the system eigenfrequency. The lift is dominated by the flap motion and the effective angle of attack due to the motion introduces phase shifts of the lift signal with respect to the flap phase angle.The experiment has been set up and executed with the necessary precision, but small ambiguities are found in the lift and drag disqualifying the data for validation. Nevertheless the data set provides a clear insight into typical loads and motions and can be used for comparative studies. It can also be used to (re)design future experiments to improve the quality of the data to the desired level of accuracy for validation.     

CONTROL OF FLOW PAST BLUFF BODIES USING ROTATING CONTROL CYLINDERS

Computational results for control of flow past a circular cylinder using small rotating cylinders are presented. A well-proven stabilized finite-element method, that has been applied to various flow problems earlier, is utilized to solve the incompressible Navier–Stokes equations in the primitive variables formulation. The formulation is first applied to study flow past an isolated rotating cylinder. Excellent match with experimental results, reported earlier, is observed. It is found that in purely two-dimensional flows, very high lift coefficients can be realized. However, it is observed, via three-dimensional Navier–Stokes simulations, that the end-effects and centrifugal instabilities along the cylinder span lead to a loss of lift and increase in drag. The aspect ratio of the cylinder plays an important role. The flow past a bluff body with two rotating control cylinders is studied using 2-D numerical simulations. The effect of the Reynolds number is studied by carrying out simulations for Re=10²and 10⁴. Finite element meshes with an adequate number of grid points are employed to resolve the flow in the gap between the main and control cylinders. Two values of the gap are considered: 0·01D and 0·075 D, where D is the diameter of the main cylinder. It is observed that when the control cylinders rotate at high speed, such that the tip speed is 5 times the free-stream speed, the flow at Re=100 achieves a steady state. For Re=10⁴, even though the flow remains unsteady, the wake is highly organized and narrower compared to the one without control. The results are in good agreement with the flow-visualization studies conducted by other researchers for bluff bodies using similar control concepts. In all the cases, a significant reduction in the overall drag coefficient and the unsteady aerodynamic forces acting on the main cylinder is observed. Results are also presented for the power requirements of the system for translation and rotation. It is found that the coefficient of power required for the rotation of control cylinders is significant for Re=100 but negligible for Re=10⁴flow. The size of the gap is found to be more critical for the Re=10⁴flows. This study brings out the relevance of the gap as a design parameter for such flow control devices.     

Numerical modelling of confined flow past a cylinder of square cross-section at various orientations

Results are presented for the unsteady, two-dimensional flow and heat transfer due to a square obstruction of diameter d located asymmetrically between the parallel sliding walls of a channel with length-to-height ratio W/H = 6·44. Analysis is based on the numerical solution of spatially and temporally second-order accurate finite difference approximations of the transport equations expressed in curvilinear co-ordinates. Laminar, constant property flow is assumed for obstruction configurations in which the blockage ratio is d/H = 0·192, the nearest-wall distances are g/d = 0·2, 0·5 and 1, the orientation angles are α=0°, 10° and 20° and the Reynolds numbers are Re=100, 500, and 1000. Preparatory testing of the numerical procedure was performed for a variety of documented flows to verify its physiconumerical accuracy and obtain estimates of the residual grid-dependent uncertainties in the variables calculated. Heat transfer, drag and lift coefficients and Strouhal numbers for the present flow were finally calculated to within 4%–7% of their grid-dependent values using non-uniformly spaced grids consisting of (x=99, y=55) nodes. Above a critical value of the Reynolds number, which depends on the geometrical parameters, the flow is characterized by alternate vortex shedding from the obstruction top and bottom surfaces. Streamline, vorticity and particle streakline plots provide qualitative impressions of the unsteady vortical flow. Especially noteworthy are the extremes in the lift coefficient which ranges from large positive values for an obstruction with g/d=0·2 and α=10° to negative values for one with g/d=0·5 and α=0°. Both the drag and lift coefficients as well as the Strouhal number exhibit non-monotonic variations with respect to the parameters explored. Asymmetries in the obstruction location and orientation account for relatively large vortex-induced periodic variations in heat transfer, especially along the wall nearest the obstruction. Notable differences are also predicted for the heat transfer coefficients of the individual obstruction surfaces as a function of the orientation angle.     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

Reynolds number, thickness and camber effects on flapping airfoil propulsion

The effect of varying airfoil thickness and camber on plunging and combined pitching and plunging airfoil propulsion at Reynolds number Re=200, 2000, 20 000 and 2×10⁶ was studied by numerical simulations for fully laminar and fully turbulent flow regimes. The thickness study was performed on 2-D NACA symmetric airfoils with 6-50% thick sections undergoing pure plunging motion at reduced frequency k=2 and amplitudes h=0.25 and 0.5, and for combined pitching and plunging motion at k=2, h=0.5, phase ?=90°, pitch angle θ_o=15° and 30° and the pitch axis was located at 1/3 of chord from leading edge. At Re=200 for motions where positive thrust is generated, thin airfoils outperform thick airfoils. At higher Re significant gains could be achieved both in thrust generation and propulsive efficiency by using a thicker airfoil section for plunging and combined motion with low pitch amplitude. The camber study was performed on 2-D NACA airfoils with varying camber locations undergoing pure plunging motion at k=2, h=0.5 and Re=20 000. Little variation in thrust performance was found with camber. The underlying physics behind the alteration in propulsive performance between low and high Reynolds numbers has been explored by comparing viscous Navier-Stokes and inviscid panel method results. The role of leading edge vortices was found to be key to the observed performance variation.     

飞行器纵向俯仰气动特性计算研究

应用有限体积方法求解三维可压缩雷诺平均N-S方程,计算了巡航导弹外形飞行器作小振幅俯仰运动时的动态绕流流场和空气动力特性,开展了导弹绕不同转轴、以不同频率和在不同迎角范围内进行俯仰运动的非定常气动力迟滞特性研究。计算结果表明,当导弹作快速俯仰运动时,在上仰和下俯过程中的同一迎角瞬间,绕导弹流场流动明显不同,表现出明显的非定常迟滞特性。导弹的非定常气动力迟滞特性随俯仰运动频率的增大明显增强,且气动力迟滞曲线随着俯仰轴位置的变化而变化。在同一减缩频率下,导弹在不同迎角范围内作周期俯仰运动时,相同的运动相位角所对应的升力系数对迎角的导数是一致的,而不同减缩频率下升力系数对迎角的导数随运动相位角变化曲线明显不同。     

Effect of flow regime change from subsonic to transonic on the air loads of an oscillating airfoil

In this research, the effect of flow regime change from subsonic to transonic on the air loads of a pitching NACA0012 airfoil is investigated. To do this, the effect of change in flow regime on the lift and pitching moment coefficients hysteresis cycles is studied. The harmonic balance approach is utilized for numerical calculation due to its low computational time. Verifications are also made with previous works and good agreements are observed. The assessment of flow regime change on the aforementioned hysteresis cycles is accomplished in the Mach number range of M=0.65–0.755. The reduced frequency and pitch amplitude also vary from k=0.03 to 0.1 and α₀=1–2.51°, respectively. Results show that the effect of increase in Mach number is to increase and decrease the lift coefficient during downstroke and upstroke, respectively, whereas at low reduced frequencies, the effect of increase in Mach number may lead to a reverse manner when airfoil moves toward its extremum angle of attack. Results also reveal that as the pitch amplitude varies, the shape of lift coefficient hysteresis cycle depends more on the pitch amplitude than on the appearance of shock. It is shown that as the Mach number increases, the incidence angles correspond to the extremum pitching moment, and depending on the reduced frequency, lie between zero and extremum angle of attack. These incidence angles shift toward the extremum angle of attack as the reduced frequency decreases. Results also show that the increase in pitch amplitude at low Mach number, in such a way that leads to the formation of shock around the extremum angle of attack, causes the extremum pitching moment to appear around these angles and at high Mach number, depending on the reduced frequency, the extremum pitching moment incidence angles would be between zero and extremum incidence angle.     

Effect of flexure on aerodynamic propulsive efficiency of flapping flexible airfoil

The aim of present study is to investigate the effect of chord-wise flexure amplitude on unsteady aerodynamic characteristics for a flapping airfoil with various combinations of Reynolds number and reduced frequency. Unsteady, viscous flows over a single flexible airfoil in plunge motion are computed using conformal hybrid meshes. The dynamic mesh technique is applied to illustrate the deformation modes of the flexible flapping airfoil. In order to investigate the influence of the flexure amplitude on the aerodynamic performance of the flapping airfoil, the present study considers eight different flexure amplitudes (a₀) ranging from 0 to 0.7 in intervals of 0.1 under conditions of Re=10⁴, reduced frequency k=2, and dimensionless plunge amplitude h₀=0.4. The computed unsteady flow fields clearly reveal the formation and evolution of a pair of leading edge vortices along the body of the flexible airfoil as it undergoes plunge motion. Thrust-indicative wake structures are generated when the flexure amplitude of the airfoil is less than 0.5 of the chord length. An enhancement in the propulsive efficiency is observed for a flapping airfoil with flexure amplitude of 0.3 of the chord length. This study also calculates the propulsive efficiency and thrust under various Reynolds numbers and reduced frequency conditions. The results indicate that the propulsive efficiency has a strong correlation with the reduced frequency. It is found that the flow conditions which yield the highest propulsive efficiency correspond to Strouhal number St of 0.255.     

Control of the aerodynamic characteristics of a profile oscillating with respect to the incidence angle in subsonic viscous flow

The results of a numerical simulation of the unsteady subsonic viscous gas flow around a two-dimensional profile oscillating with respect to the incidence angle are presented and the possibility of controlling the nonstationary aerodynamic characteristics is considered. The hysteresis phenomena typical of oscillatory profile motions are investigated, the dependence of the lift force and drag is found for various laws of periodic variation of the incidence angle with time, and the effect of the frequency and amplitude of the angular profile oscillations on the shape of the hysteresis curves is studied. The calculations were based on the numerical solution of the nonstationary Navier-Stokes equations averaged in the Reynolds sense (Reynolds equations) which were closed using the k-ω turbulence model with modeling of the laminar/turbulent transition.     

Investigation of steady plasma actuation effect on aerodynamic coefficients of oscillating airfoil at low Reynolds number

In this work, numerical study of two dimensional laminar incompressible flow around an oscillating NACA0012 airfoil is proceeded using the open source code Open FOAM. Oscillatory motion types including pitching and flapping are considered. Reynolds number for these motions is assumed to be 12000 and effects of these motions and also different unsteady parameters such as amplitude and reduced frequency on aerodynamic coefficients are studied. For flow control on airfoil, dielectric barrier discharge plasma actuator is used in two different positions on airfoil and its effect is compared for the two types of considered oscillating motions. It is observed that in pitching motion, imposing plasma leads to an improvement in aerodynamic coefficients, but it does not have any positive effect on flapping motion.Also, for the amplitudes and frequencies investigated in this paper, the trailing edge plasma had a more desirable effect than other positions.     

APPLICATION OF A MODIFIED k–ε MODEL TO THE PREDICTION OF AERODYNAMIC CHARACTERISTICS OF RECTANGULAR CROSS-SECTION CYLINDERS

Aerodynamic characteristics of rectangular cylinders of infinite length and various breadth-to-depth ratios ranging from B/D=0·6 to 8·0 were investigated numerically by a two-layerk –ε model with a modified k -production term. Although the numerical method is two-dimensional (2-D), a physically reasonable smooth, periodic vortex shedding was obtained, even in the range of high Reynolds number. This kind of behaviour cannot normally be simulated by ordinary 2-D analyses which do not incorporate a turbulence model. Various typical aerodynamic features were successfully obtained, particularly including the discontinuity in Strouhal number at the critical section of B/D=2·8 and 6·0. Also, drag coefficients and distributions of mean surface pressure were in good agreement with results of experiments and 3-D analyses for the range of B/D ratios. However, as regards the prediction of pressure and force fluctuations, since the currently available Reynolds-Averaged Navier–Stokes models are able to analyse only periodic components and not stochastic components, the total fluctuations in surface pressure and aerodynamic lift force were considerably underestimated in some cases, compared with those measured in experiments and calculated from 3-D analyses.     

Flow structure of knuckling effect in footballs

The flight trajectory of a non-spinning or slow-spinning soccer ball might fluctuate in unpredictable ways, as for example, in the many free kicks of C. Ronaldo. Such anomalous horizontal shaking or rapid falling is termed the ‘knuckling effect’. However, the aerodynamic properties and boundary-layer dynamics affecting a ball during the knuckling effect are not well understood. In this study, we analyse the characteristics of the vortex structure of a soccer ball subject to the knuckling effect (knuckleball), using high-speed video images and smoke-generating agents. Two high-speed video cameras were set at one side and in front of the ball trajectory between the ball position and the goal; further, photographs were taken at 1000 fps and a resolution of 1024×512 pixels. Although in a previous study (Taneda, 1978), shedding of horseshoe vortices was observed for smooth spheres in the Reynolds number (Re) range of 3.8×10⁵<Re<10⁶, in the case of the soccer ball, the vortex structure, which consisted of distorted loop vortices, appeared in the wake behind the ball in the supercritical Re number region. Moreover, after the knuckleballs were airborne, large-scale undulations were observed in the vortex trail visualised with a smoke technique. On the other hand, aerodynamic forces acting on the ball were estimated from the data of the ball’s flight trajectory, and a statistically high correlation (r=0.94, p<0.01) between the fluctuation frequency of the lift and side forces and the undulation frequency of the vortex trail was shown to exist. This fact suggests that the phenomenon of large-scale undulations of the vortex trail is closely related to the cause of the unsteady aerodynamic forces acting on the knuckle ball.     

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

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