Measurement of parallel blade–vortex interaction at low Reynolds numbers

In this study parallel blade–vortex interaction for a Schmidt-propeller configuration has been examined using particle image velocimetry (PIV). This tandem configuration consists of a leading airfoil (forefoil), used to generate a vortical wake of leading-edge vortices (LEVs) and trailing-edge vortices (TEVs) through a pitching or plunging motion, and a trailing airfoil (hindfoil), held fixed with a specified angle of attack and vertical spacing in its wake. The hindfoil incidence (loading) and not the vertical spacing to the incoming vortical wake has been found to dictate the nature of the interaction (inviscid vs. viscous). For cases where the vortex–blade offset is small and the hindfoil is loaded, vortex distortion and vortex-induced separations are observed. By tracking the circulation of the LEV and TEV, it has been found that the vortices are strengthened for the tandem arrangement and in certain cases dissipate quicker in the wake when interacting with the hindfoil. Time-averaged forces obtained using a standard control-volume analysis are then obtained and used to evaluate these vortex-interaction cases. A subsequent analysis of the varying pressure distribution over the suction side of the hindfoil is performed by integrating the Navier–Stokes equations through the velocity field. This allows for a direct comparison of the vortex-induced loading for the various configurations.     

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.     

On the unsteady motion and stability of a heaving airfoil in ground effect

This study explores the fluid mechanics and force generation capabilities of an inverted heaving airfoil placed close to a moving ground using a URANS solver with the Spalart-Allmaras turbulence model. By varying the mean ground clearance and motion frequency of the airfoil, it was possible to construct a frequency-height diagram of the various forces acting on the airfoil. The ground was found to enhance the downforce and reduce the drag with respect to freestream. The unsteady motion induces hysteresis in the forces’ behaviour. At moderate ground clearance, the hysteresis increases with frequency and the airfoil loses energy to the flow, resulting in a stabilizing motion. By analogy with a pitching motion, the airfoil stalls in close proximity to the ground. At low frequencies, the motion is unstable and could lead to stall flutter. A stall flutter analysis was undertaken. At higher frequencies, inviscid effects overcome the large separation and the motion becomes stable. Forced trailing edge vortex shedding appears at high frequencies. The shedding mechanism seems to be independent of ground proximity. However, the wake is altered at low heights as a result of an interaction between the vortices and the ground.     

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.     

Time-resolved scanning tomography PIV measurements around a flapping wing

The three-dimensional flow that develops around a finite flapping wing is investigated using a tomographic scanning PIV technique. The acquisition and correlation processes employed to achieve such measurements have been carefully validated. Issues regarding the relevant timescales of the flow and the spanwise space-resolution are addressed. Results obtained on a hovering flapping wing whose plunging phase is described by a rectilinear motion highlight the influence of the free end condition and the formation of the tip vortex on the leading edge vortices behavior, wing/wake interactions, and wake stabilization.     

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.     

High fidelity numerical simulation of airfoil thickness and kinematics effects on flapping airfoil propulsion

High-fidelity numerical simulations with the spectral difference (SD) method are carried out to investigate the unsteady flow over a series of oscillating NACA 4-digit airfoils. Airfoil thickness and kinematics effects on the flapping airfoil propulsion are highlighted. It is confirmed that the aerodynamic performance of airfoils with different thickness can be very different under the same kinematics. Distinct evolutionary patterns of vortical structures are analyzed to unveil the underlying flow physics behind the diverse flow phenomena associated with different airfoil thickness and kinematics and reveal the synthetic effects of airfoil thickness and kinematics on the propulsive performance. Thickness effects at various reduced frequencies and Strouhal numbers for the same chord length based Reynolds number (=1200) are then discussed in detail. It is found that at relatively small Strouhal number (=0.3), for all types of airfoils with the combined pitching and plunging motion (pitch angle 20°, the pitch axis located at one third of chord length from the leading edge, pitch leading plunge by 75°), low reduced frequency (=1) is conducive for both the thrust production and propulsive efficiency. Moreover, relatively thin airfoils (e.g. NACA0006) can generate larger thrust and maintain higher propulsive efficiency than thick airfoils (e.g. NACA0030). However, with the same kinematics but at relatively large Strouhal number (=0.45), it is found that airfoils with different thickness exhibit diverse trend on thrust production and propulsive efficiency, especially at large reduced frequency (=3.5). Results on effects of airfoil thickness based Reynolds numbers indicate that relative thin airfoils show superior propulsion performance in the tested Reynolds number range. The evolution of leading edge vortices and the interaction between the leading and trailing edge vortices play key roles in flapping airfoil propulsive performance.     

Dynamic wind loads and wake characteristics of a wind turbine model in an atmospheric boundary layer wind

An experimental study was conducted to characterize the dynamic wind loads and evolution of the unsteady vortex and turbulent flow structures in the near wake of a horizontal axis wind turbine model placed in an atmospheric boundary layer wind tunnel. In addition to measuring dynamic wind loads (i.e., aerodynamic forces and bending moments) acting on the wind turbine model by using a high-sensitive force-moment sensor unit, a high-resolution digital particle image velocimetry (PIV) system was used to achieve flow field measurements to quantify the characteristics of the turbulent vortex flow in the near wake of the wind turbine model. Besides conducting “free-run” PIV measurements to determine the ensemble-averaged statistics of the flow quantities such as mean velocity, Reynolds stress, and turbulence kinetic energy (TKE) distributions in the wake flow, “phase-locked” PIV measurements were also performed to elucidate further details about evolution of the unsteady vortex structures in the wake flow in relation to the position of the rotating turbine blades. The effects of the tip-speed-ratio of the wind turbine model on the dynamic wind loads and wake flow characteristics were quantified in the terms of the variations of the aerodynamic thrust and bending moment coefficients of the wind turbine model, the evolution of the helical tip vortices and the unsteady vortices shedding from the blade roots and turbine nacelle, the deceleration of the incoming airflows after passing the rotation disk of the turbine blades, the TKE and Reynolds stress distributions in the near wake of the wind turbine model. The detailed flow field measurements were correlated with the dynamic wind load measurements to elucidate underlying physics in order to gain further insight into the characteristics of the dynamic wind loads and turbulent vortex flows in the wakes of wind turbines for the optimal design of the wind turbines operating in atmospheric boundary layer winds.     

An investigation into the effects of unsteady parameters on the aerodynamics of a low Reynolds number pitching airfoil

The growing applications of low Reynolds number (LRN) operating vehicles impose the need for accurate LRN flow solutions. These applications usually involve complex unsteady phenomena, which depend on the kinematics of the vehicle such as pitching, plunging, and flapping of a wing. The objective of the present study is to address the issues related to LRN aerodynamics of a harmonically pitching NACA0012 airfoil. To this end, the influence of unsteady parameters, namely, amplitude of oscillation, d, reduced frequency, k, and Reynolds number, Re, on the aerodynamic performance of the model is investigated. Computational fluid dynamics (CFD) is utilized to solve Navier–Stokes (N–S) equations discretized based on the Finite Volume Method (FVM). The resulting instantaneous lift coefficients are compared with analytical data from Theodorsen’s method. The simulation results reveal that d, k, and Re are of great importance in the aerodynamic performance of the system, as they affect the maximum lift coefficients, hysteresis loops, strength, and number of the generated vortices within the harmonic motion, and the extent of the so-called figure-of-eight phenomenon region. Thus, achieving the optimum lift coefficients demands a careful selection of these parameters.     

Investigation of the unsteady flow velocity field above an airfoil pitching under deep dynamic stall conditions

The unsteady flow field above a NACA 0012 airfoil pitching under deep dynamic stall conditions has been investigated in a low-speed wind tunnel by means of particle image velocimetry. The measurements of the instantaneous flow velocity field show the characteristic features of the dynamic stall process: formation and development of an organized vortex structure for increasing incidences and the subsequent separation. Vorticity and divergence estimated from the measured data give a good insight into the complex flow behaviour during the downstroke motion. Furthermore, small-scale structures could be observed in the separated flow field and even within the dynamic stall vortex.The authors would like to thank Dr. Schäfer (ISL) for his support in organizing the cooperative measurements, Mr. Seyb (DLR) for his help during the recording of PIV images, Dr. Bretthauer (DLR) and Mr. Vollmers (DLR) for his assistance during the phase of evaluation and post processing of the PIV recordings and Dr. Geißler (DLR) for helpful discussions on the dynamic stall problem.     

PIV study of flow around unsteady airfoil with dynamic trailing-edge flap deflection

The flow around an oscillating NACA 0015 airfoil with prescheduled trailing-edge flap motion control was investigated by using particle image velocimetry (PIV). Aerodynamic load coefficients, obtained via surface pressure measurements, were also acquired to supplement the PIV results. The results demonstrate that upward flap deflections led to an improved negative peak pitching moment coefficient C _m,peak, mainly as a consequence of the increased suction pressure on the lower surface of the flap. The behavior of the leading-edge vortex (LEV) was largely unaffected. Its strength was, however, reduced slightly compared to that of the uncontrolled airfoil. No trailing-edge vortex was observed. For downward flap deflection, the strength of the LEV was found to be slightly increased. A favorable increase in C _l,max, as a consequence of downward flap-induced positive camber effects, accompanied by a detrimental increase in the nose-down C _m,peak, due to the large pressure increase on the lower surface of the flap, was also observed.     

High-fidelity simulations of moving and flexible airfoils at low Reynolds numbers

The present paper highlights results derived from the application of a high-fidelity simulation technique to the analysis of low-Reynolds-number transitional flows over moving and flexible canonical configurations motivated by small natural and man-made flyers. This effort addresses three separate fluid dynamic phenomena relevant to small fliers, including: laminar separation and transition over a stationary airfoil, transition effects on the dynamic stall vortex generated by a plunging airfoil, and the effect of flexibility on the flow structure above a membrane airfoil. The specific cases were also selected to permit comparison with available experimental measurements. First, the process of transition on a stationary SD7003 airfoil section over a range of Reynolds numbers and angles of attack is considered. Prior to stall, the flow exhibits a separated shear layer which rolls up into spanwise vortices. These vortices subsequently undergo spanwise instabilities, and ultimately breakdown into fine-scale turbulent structures as the boundary layer reattaches to the airfoil surface. In a time-averaged sense, the flow displays a closed laminar separation bubble which moves upstream and contracts in size with increasing angle of attack for a fixed Reynolds number. For a fixed angle of attack, as the Reynolds number decreases, the laminar separation bubble grows in vertical extent producing a significant increase in drag. For the lowest Reynolds number considered (Re _c  = 10⁴), transition does not occur over the airfoil at moderate angles of attack prior to stall. Next, the impact of a prescribed high-frequency small-amplitude plunging motion on the transitional flow over the SD7003 airfoil is investigated. The motion-induced high angle of attack results in unsteady separation in the leading edge and in the formation of dynamic-stall-like vortices which convect downstream close to the airfoil. At the lowest value of Reynolds number (Re _c  = 10⁴), transition effects are observed to be minor and the dynamic stall vortex system remains fairly coherent. For Re _c  = 4 × 10⁴, the dynamic-stall vortex system is laminar at is inception, however shortly afterwards, it experiences an abrupt breakdown associated with the onset of spanwise instability effects. The computed phased-averaged structures for both values of Reynolds number are found to be in good agreement with the experimental data. Finally, the effect of structural compliance on the unsteady flow past a membrane airfoil is investigated. The membrane deformation results in mean camber and large fluctuations which improve aerodynamic performance. Larger values of lift and a delay in stall are achieved relative to a rigid airfoil configuration. For Re _c = 4.85 × 10⁴, it is shown that correct prediction of the transitional process is critical to capturing the proper membrane structural response.     

Experimental study on the flow regimes past a confined prism undergoing self-sustained oscillations

An experimental study based on Particle Image Velocimetry (PIV) is presented with the objective of studying the flow regimes that appear in the flow past a confined prism undergoing self-sustained oscillations at low Reynolds numbers (Re). The square-section prism, placed inside a 3D square cross-section vertical channel with a confinement ratio of 1/2.5, was tethered to the channel walls and, therefore, it was allowed to move freely transverse to the incoming flow. Re (based on the prism cross-section height) was varied in the range from 100 to 700. Three different prism to fluid density ratios (m¹) were considered: 0.56, 0.70, and 0.91. These two parameters, Re and m¹, were used to map the results obtained. In particular, it was found that five different regimes appear: (1) steady prism with steady recirculation bubble, (2) steady prism with unsteady vortex shedding wake, (3) large amplitude low frequency oscillating prism with unsteady vortex shedding wake, (4) small amplitude high frequency oscillating prism with unsteady vortex shedding wake, and (5) irregular/chaotic motion of both the prism and the wake. The PIV results and associated numerical simulations were used to analyze the different prism and wake states.     

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.     

Parameter dependence of vortex interactions on a two-dimensional plunging plate

The structure and dynamics of the flow field created by a plunging flat-plate airfoil are investigated at a chord Reynolds number of 10,000 while varying plunge amplitude and Strouhal number. Digital particle image velocimetry measurements are used to characterize the shedding patterns and the interactions between the leading- and trailing-edge vortex structures (LEV and TEV), resulting in the development of a wake classification system based on the nature and timing of interactions between the leading- and trailing-edge vortices. The streamwise advancement of the LEV during a plunge cycle and its resulting interaction with the TEV is primarily dependent on reduced frequency; however, for Strouhal numbers above approximately 0.4, significant changes are observed in the formation of vortices shed from the leading and trailing edges, as well as the circulation of the leading-edge vortex. The functional form of the relationship between leading-edge vortex circulation and Strouhal number suggests that the Strouhal number dependence is more specifically a manifestation of the effective angle of attack. Comparison with low-Reynolds-number studies of plunging airfoil aerodynamics reveals a high degree of consistency and suggests applicability of the classification system beyond the range examined in the present work.     

A study on the mechanism of high-lift generation by an airfoil in unsteady motion at low reynolds number

The aerodynamic force and flow structure of NACA 0012 airfoil performing an unsteady motion at low Reynolds number (Re=100) are calculated by solving Navier-Stokes equations. The motion consists of three parts: the first translation, rotation and the second translation in the direction opposite to the first. The rotation and the second translation in this motion are expected to represent the rotation and translation of the wing-section of a hovering insect. The flow structure is used in combination with the theory of vorticity dynamics to explain the generation of unsteady aerodynamic force in the motion. During the rotation, due to the creation of strong vortices in short time, large aerodynamic force is produced and the force is almost normal to the airfoil chord. During the second translation, large lift coefficient can be maintained for certain time period and , the lift coefficient averaged over four chord lengths of travel, is larger than 2 (the corresponding steady-state lift coefficient is only 0.9). The large lift coefficient is due to two effects. The first is the delayed shedding of the stall vortex. The second is that the vortices created during the airfoil rotation and in the near wake left by previous translation form a short “vortex street” in front of the airfoil and the “vortex street” induces a “wind”; against this “wind” the airfoil translates, increasing its relative speed. The above results provide insights to the understanding of the mechanism of high-lift generation by a hovering insect. The project supported by the National Natural Science Foundation of China (19725210)     

Measurements of the flow over a low-aspect-ratio cylinder mounted on a ground plane

The flow over a finite-height cylinder of aspect ratio 1, with one end mounted on a ground plane and the other end free, has been studied by means of surface flow visualisation, particle image velocimetry (PIV) and surface pressure measurements. The diameter-based Reynolds number was 200,000. The mean flow topology has been identified in three areas: the horseshoe vortex system, the separated flow over the free-end and the wake region. Evidence is shown for the existence of a three-horseshoe vortex system, while the mean flow over the free-end consists of an arch vortex with its bases on the forward half of the free-end. There are two tip vortices coming off the free-end. The wake region is found to be highly unsteady, with considerable variation from the mean flow.     

Scaling of flow separation on a pitching low aspect ratio plate

We use two different dye injection approaches, in two different water tunnels, to visualize the formation and subsequent evolution of leading-edge vortices and related separated structures, for a pitching low aspect ratio plate. The motion is a smoothed linear pitch ramp from 0° to 40° incidence, brief hold, and return to 0°, executed at reduced pitch rates ranging from 0.1 to 0.35 and about various pivot locations. All cases evince a leading edge vortex with pronounced axial flow, which leads to formation of large-scale, three-dimensional flow structures, culminating in a large vortical structure centered at the wing symmetry plane. Pitch is also compared to plunge, where the plunge-induced angle of attack is taken as the geometric pitch incidence angle, ignoring pitch-rate effects. At successively increasing values of convective time C/U, the three-dimensional patterns of the flow structure are remarkably similar for the pitching and plunging motions. The similarity of these patterns persists, though they are shifted in time, for variation of either the location of the pitching axis or the dimensionless pitch rate.     

Assessment of 2D/3D numerical modeling for deep dynamic stall experiments

The results of computational fluid dynamics (CFD) simulations in two and three spatial dimensions are compared to pressure measurements and particle image velocimetry (PIV) flow surveys to assess the suitability of numerical models for the simulation of deep dynamic stall experiments carried out on a pitching NACA 23012 airfoil. A sinusoidal pitching motion with a 10° amplitude and a reduced frequency of 0.1 is imposed around two different mean angles of attack of 10° and 15°. The comparison of the airloads curves and of the pressure distribution over the airfoil surface shows that a three-dimensional numerical model can better reproduce the flow structures and the airfoil performance for the deep dynamic stall regime. Also, the vortical structures observed by PIV in the flow field are better captured by the three-dimensional model. This feature highlighted the relevance of three-dimensional effects on the flow field in deep dynamic stall.     

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.