Turbulence modelling of the flow past a pitching NACA0012 airfoil at and Reynolds numbers

This paper provides a study of the NACA0012 dynamic stall at Reynolds numbers 10⁵ and 10⁶ by means of two- and three-dimensional numerical simulations. The turbulence effect on the dynamic stall is studied by statistical modelling. The results are compared with experiments concerning each test case. Standard URANS turbulence modelling have shown a quite dissipative character that attenuates the instabilities and the vortex structures related to the dynamic stall. The URANS approach Organised Eddy Simulation (OES) has shown an improved behaviour at the high Reynolds number range. Emphasis is given to the physical analysis of the three-dimensional dynamic stall structure, for which there exist few numerical results in the literature, as far as the Reynolds number range is concerned. This study has shown that the downstroke phases of the pitching motion are subjected to strong three-dimensional turbulence effects along the span, whereas the flow is practically two-dimensional during the upstroke motion.     

Self-sustained aeroelastic oscillations of a NACA0012 airfoil at low-to-moderate Reynolds numbers

A wind tunnel experimental investigation of self-sustained oscillations of an aeroelastic NACA0012 airfoil occurring in the transitional Re regime is presented. To the authors’ knowledge this is the first time that aeroelastic limit cycle oscillations (LCOs) associated with low Re effects have been systematically studied and reported in the public literature. While the aeroelastic apparatus is capable of two-degree-of-freedom pitch-plunge motion, the present work concerns only the motion of the airfoil when it is constrained to rotate in pure pitch. The structural stiffness is varied as well as the position of the elastic axis; other parameters such as surface roughness, turbulence intensity and initial conditions are also briefly discussed. In conjunction with the pitch measurements, the flow is also recorded using hot-wire anemometry located in the wake at a distance of one chord aft of the trailing edge. It is observed that for a limited range of chord-based Reynolds numbers, 4.5×10⁴Re_c1.3×10⁵, steady state self-sustained oscillations are observed. Below and above that range, these oscillations do not appear. They are characterized by a well-behaved harmonic motion, whose frequency can be related to the aeroelastic natural frequency, low amplitude (θ_max<5.5°) and some sensitivity to flow perturbations and initial conditions. Furthermore, hot-wire measurements for the wing held fixed show that no periodicity in the undisturbed free-stream nor in the wake account for the oscillations. Overall, these observations suggest that laminar separation plays a role in the oscillations, either in the form of trailing edge separation or due to the presence of a laminar separation bubble.     

Numerical and experimental research of flow control on an NACA 0012 airfoil by local vibration

A flow control technique by local vibration is proposed to improve the aerodynamic performance of a typical airfoil NACA 0012. Both wind-tunnel experiments and a large eddy simulation(LES) are carried out to study the effects of local vibration on drag reduction over a wide range of angles of attack. The application parameters of local vibration on the upper surface of the airfoil are first evaluated by numerical simulations.The mounted position is chosen at 0.065–0.09 of chord length from the leading edge.The influence of oscillation frequency is investigated both by numerical simulations and experiments. The optimal frequencies are near the dominant frequencies of shear layer vortices and wake vortices. The patterns of shear vortices caused by local vibration are also studied to determine the drag reduction mechanism of this flow control method. The results indicate that local vibration can improve the aerodynamic performance of the airfoil. In particular, it can reduce the drag by changing the vortex generation patterns.     

Stochastic estimation of flow near the trailing edge of a NACA0012 airfoil

A stochastic estimation technique has been applied to simultaneously acquired data of velocity and surface pressure as a tool to identify the sources of wall-pressure fluctuations. The measurements have been done on a NACA0012 airfoil at a Reynolds number of Re _c  = 2 × 10⁵, based on the chord of the airfoil, where a separated laminar boundary layer was present. By performing simultaneous measurements of the surface pressure fluctuations and of the velocity field in the boundary layer and wake of the airfoil, the wall-pressure sources near the trailing edge (TE) have been studied. The mechanisms and flow structures associated with the generation of the surface pressure have been investigated. The “quasi-instantaneous” velocity field resulting from the application of the technique has led to a picture of the evolution in time of the convecting surface pressure generating flow structures and revealed information about the sources of the wall-pressure fluctuations, their nature and variability. These sources are closely related to those of the radiated noise from the TE of an airfoil and to the vibration issues encountered in ship hulls for example. The NACA0012 airfoil had a 30 cm chord and aspect ratio of 1.     

Computation of unsteady viscous flow around a locally flexible airfoil at low Reynolds number

A numerical method for fluid–structure interaction is presented for the analysis of unsteady viscous flow over a locally flexible airfoil. The Navier–Stokes equations are solved by ALE–CBS algorithm, coupling with a structural solver with large deformation. Following the validation of the method, a numerical example for the flight of micro-air vehicles at low Reynolds number is chosen for the computation. The coupling effect of flexible structure with different elastic stiffness on aerodynamic performance is demonstrated. A noticeable camber effect is induced by the deflection of the structure as the elastic stiffness of the structure goes smaller. Moreover, when the vibrating frequencies of the structure with smaller elastic stiffness have a close correlation with the shedding frequencies, the positive impact of the vibration of local flexible surface on the lift of the airfoil is highlighted, which results from the formation of the coherent vortices.     

Flow control over a NACA 0012 airfoil using dielectric-barrier-discharge plasma actuator with a Gurney flap

Flow control study of a NACA 0012 airfoil with a Gurney flap was carried out in a wind tunnel, where it was demonstrated that a dielectric-barrier-discharge (DBD) plasma actuator attached to the flap could increase the lift further, but with a small drag penalty. Time-resolved PIV measurements of the near-wake region indicated that the plasma forcing shifted the wake downwards, reducing its recirculation length. Analysis of wake vortex dynamics suggested that the plasma actuator initially amplified the lower wake shear layer by adding momentum along the downstream surface of the Gurney flap. This enhanced mutual entrainment between the upper and lower wake vortices, leading to an increase in lift on the airfoil.     

Measurement of tonal-noise characteristics and periodic flow structure around NACA0018 airfoil

The characteristics of tonal noise and the variations of flow structure around NACA0018 airfoil in a uniform flow are studied by means of simultaneous measurement of noise and velocity field by particle-image velocimetry to understand the generation mechanism of tonal noise. Measurements are made on the noise characteristics, the phase-averaged velocity field with respect to the noise signal, and the cross-correlation contour of velocity fluctuations and noise signal. These experimental results indicate that the tonal noise is generated from the periodic vortex structure on the pressure surface of the airfoil near the trailing edge of the airfoil. It is found that the vortex structure is highly correlated with the noise signal, which indicates the presence of noise-source distribution on the pressure surface. The vorticity distribution on the pressure surface breaks down near the trailing edge of the airfoil and forms a staggered vortex street in the wake of the airfoil.     

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.     

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.     

Simulation of dynamic stall for a NACA 0012 airfoil using a vortex method

The unsteady, incompressible, viscous laminar flow over a NACA 0012 airfoil is simulated, and the effects of several parameters investigated. A vortex method is used to solve the two-dimensional Navier–Stokes equations in the vorticity/stream-function form. By applying an operator-splitting method, the “convection” and “diffusion” equations are solved sequentially at each time step. The convection equation is solved using the vortex-in-cell method, and the diffusion equation using a second-order ADI finite difference scheme. The airfoil profile is obtained by mapping a circle in the computational domain into the physical domain through a Joukowski transformation. The effects of several parameters are investigated, such as the reduced frequency, mean angle of attack, location of pitch axis, and the Reynolds number. It is observed that the reduced frequency has the most influence on the flow field.     

Physics of forced unsteady flow for a NACA 0015 airfoil undergoing constant-rate pitch-up motion*

The unsteady Navier-Stokes (NS) analysis of Osswald, Ghia and Ghia in velocity-vorticity variables is modified to study the dynamic stall phenomenon for a NACA 0015 airfoil undergoing constant Ω₀ pitch-up maneuvers at Reynolds number Re =10 000 and 45000. The use of third-order accurate biased upwind differencing for the nonlinear convective terms in the vorticity transport equation removes the spurious oscillations observed in the earlier studies by the authors for these values of Re. The fully implicit and vectorized ADI-BGE method of the authors is used to solve the unsteady NS equations. Instantaneous inertial surface vorticity, which is an invariant of the choice of reference frame selected, is employed to determine the location of separation of the boundary-layer flow on the suction surface; also a separation bubble embedded within the boundary layer is observed for both cases somewhere between the leading edge and the quarter-chord point. Primary, secondary, tertiary and quarternary vortices have been observed before the dynamic-stall vortex evolves and gathers its maximum strength.     

Direct numerical simulation of a NACA0012 in full stall

This work aims at investigating the mechanisms of separation and the transition to turbulence in the separated shear-layer of aerodynamic profiles, while at the same time to gain insight into coherent structures formed in the separated zone at low-to-moderate Reynolds numbers. To do this, direct numerical simulations of the flow past a NACA0012 airfoil at Reynolds numbers Re = 50,000 (based on the free-stream velocity and the airfoil chord) and angles of attack AOA = 9.25° and AOA = 12° have been carried out. At low-to-moderate Reynolds numbers, NACA0012 exhibits a combination of leading-edge/trailing-edge stall which causes the massive separation of the flow on the suction side of the airfoil. The initially laminar shear layer undergoes transition to turbulence and vortices formed are shed forming a von Kármán like vortex street in the airfoil wake. The main characteristics of this flow together with its main features, including power spectra of a set of selected monitoring probes at different positions on the suction side and in the wake of the airfoil are provided and discussed in detail.     

Numerical evaluation of passive control of shock wave/boundary layer interaction on NACA0012 airfoil using jagged wall

Shock formation due to flow compressibility and its interaction with boundary layers has adverse effects on aerodynamic characteristics, such as drag increase and flow separation. The objective of this paper is to appraise the prac-ticability of weakening shock waves and, hence, reducing the wave drag in transonic flight regime using a two-dimensional jagged wall and thereby to gain an appropriate jagged wall shape for future empirical study. Different shapes of the jagged wall, including rectangular, circular, and triangular shapes, were employed. The numerical method was validated by experimental and numerical studies involving transonic flow over the NACA0012 airfoil, and the results presented here closely match previous experimental and numerical results. The impact of parameters, including shape and the length-to-spacing ratio of a jagged wall, was studied on aerodynamic forces and flow field. The results revealed that applying a jagged wall method on the upper surface of an airfoil changes the shock structure significantly and disinte-grates it, which in turn leads to a decrease in wave drag. It was also found that the maximum drag coefficient decrease of around 17%occurs with a triangular shape, while the max-imum increase in aerodynamic efficiency (lift-to-drag ratio) of around 10%happens with a rectangular shape at an angle of attack of 2.26?.     

Duty cycle and directional jet effects of a plasma actuator on the flow control around a NACA0015 airfoil

This paper reports on the effects of a series of fluid-dynamic dielectric barrier discharge plasma actuators on a NACA0015 airfoil at high angle of attack. A set of jet actuators able to produce plasma jets with different directions (vectoring effect) and operated at different on/off duty cycle frequencies are used. The experiments are performed in a wind tunnel facility. The vectorized jet and the transient of the flow induced by unsteady duty cycle operation of each actuator are examined and the effectiveness of the actuator to recover stall condition in the range of Reynolds numbers between 1.0 × 10⁵ and 5.0 × 10⁵ (based on airfoil chord), is investigated. The actuator placed on the leading edge of the airfoil presents the most effective stall recovery. No significant effects can be observed for different orientations of the jet. An increase of the stall recovery is detected when the actuator is operated in unsteady operation mode. Moreover, the frequency of the on/off duty cycle that maximizes the stall recovery is found to be a function of the free stream velocity. This frequency seems to scale with the boundary layer thickness at the position of the actuator. A lift coefficient increase at low free stream velocities appears to linearly depend on the supply voltage.     

Numerical simulation of the flow around a simplified vehicle model with active flow control

Large-eddy simulation (LES) was used to study the influence and the resulting flow mechanisms of active flow control applied to a two-dimensional vehicle geometry. The LES results were validated against existing Particle Image Velocimetry (PIV) and force measurement data. This was followed by an exploration of the influence of flow actuation on the near-wake flow and resulting aerodynamic forces. Not only was good agreement found with the previous experimental study, but new knowledge was gained in the form of a complex interaction of the actuation with the coherent flow structures. The resulting time-averaged flow shows a strong influence of the extension of the actuation slots and the lateral solid walls on the near-wake flow structures and thereby on the resulting drag.     

低Reynolds数NACA0012翼型绕流的流动特性分析

在水槽中应用PIV测速技术研究了NACA0012翼型在Reynolds数为8200时的流动特性,重点关注了翼型绕流结构中主频和扰动增长速率随迎角的变化。结果表明,分离剪切层的扰动增长符合指数规律;且随着迎角的增大,转捩过程加速,表现为扰动增长率逐渐增大,转捩的起始位置逐渐向上游移动。在所有实验迎角情况下,流场均由脱落旋涡主导,但其主导作用随着迎角的增大而削弱。     

Numerical simulation on the NACA0018 airfoil self-noise generation

Airfoil self-noise is a common phenomenon for many engineering applications. Aiming to study the underlying mechanism of airfoil self-noise at low Mach number and moderate Reynolds number flow, a numerical investigation is presented on noise generation by flow past NACA0018 airfoil. Based on a high-order accurate numerical method, both the near-field hydrodynamics and the far-field acoustics are computed simultaneously by performing direct numerical simulation. The mean flow properties agree well with the experimental measurements. The characteristics of aerodynamic noise are investigated at various angles of attack. The obtained results show that inclining the airfoil could enlarge turbulent intensity and produce larger scale of vortices. The sound radiation is mainly towards the upper and lower directions of the airfoil surface. At higher angle of attack, the tonal noise tends to disappear and the noise spectrum displays broad-band features.     

Computational aeroelastic simulations of self-sustained pitch oscillations of a NACA0012 at transitional Reynolds numbers

The phenomenon of low amplitude self-sustained pitch oscillations in the transitional Reynolds number regime is studied numerically through unsteady, two-dimensional aeroelastic simulations. Based on the experimental data, simulations have been limited in the Reynolds number range 5.0×10⁴<Re_c<1.5×10⁵. Both laminar and URANS calculations (using the SST k–ω model with a low-Reynolds-number correction) have been performed and found to produce reasonably accurate limit cycle pitching oscillations (LCO). This investigation confirms that the laminar separation of the boundary layer near the trailing edge plays a critical role in initiating and sustaining the pitching oscillations. For this reason, the phenomenon is being labelled as laminar separation flutter. As a corollary, it is also shown that turbulence tends to inhibit their existence. Furthermore, two regimes of LCO are observed, one where the flow is laminar and separated without re-attachment, and the second for which transition has occurred followed by turbulent re-attachment. Finally, it is established that the high-frequency, shear instabilities present in the flow which lead to von Kármán vortex shedding are not crucial, nor necessary, to the maintaining mechanism of the self-sustained oscillations.