A high-resolution laser Doppler anemometer: design, qualification, and uncertainty

A new high-resolution laser Doppler anemometer (LDA) has been developed with a working distance of 350 mm, allowing operation in lab-scale wind tunnels. The measurement volume size is 35 μm in diameter by 60 μm in length, allowing resolution of the smallest turbulence scales even at fairly high Reynolds numbers. The controversial question of velocity and validation bias in LDA data is resolved with an experimental method for measuring and removing those effects. Uncertainty estimates are also derived for all the mean and Reynolds stress measurements. Received: 27 June 1999/Accepted: 30 August 2000     

Experimental study of flow past a square cylinder at high Reynolds numbers

 Measurements of two-components of velocity in the wake of a square cylinder using a hot-wire anemometer are reported. Two Reynolds numbers, namely 8700 and 17,625, have been considered. The measurements were carried out in a low-speed, low-turbulence wind tunnel. Benchmark experiments at much lower Reynolds numbers show good agreement between the present experiments and published results. At higher Reynolds numbers, the experimental data reveal anticipated trends in terms of wake recovery and turbulence decay. Both velocity and velocity fluctuations show symmetry about the wake axis. The experimental data have been compared with the large eddy simulation (LES) calculation reported by Wang et al. [University of Illinois at Urbana – Champaign (1996) Report CFD 96-03] and LDV measurements of Lyn et al. [J Fluid Mech (1995) 304: 285–319]. The agreement among the three sets is generally acceptable in terms of the time-averaged velocity components, but not the velocity fluctuations. The turbulence fluctuations in the present experiments are seen to be lower than in the referred work. The differences have been traced to factors such as the aspect ratio, blockage ratio and upstream turbulence. Experiments with increased upstream turbulence did show a reduction in the discrepancy between the present experiments and the published data. An assessment of the experimental data in terms of physical mechanisms revealed that (a) streamwise normal stresses were correlated with the vortex centers, and (b) the turbulence kinetic energy profiles are similar to the turbulence shear stress. Spectral analysis of the velocity signals was carried out in the present work. Energy transfer from the mean flow to the streamwise velocity fluctuation was confirmed in the near wake. A redistribution of the kinetic energy between the streamwise and transverse components of velocity over a longer distance downstream was subsequently observed. Received: 17 May 1999/Accepted: 29 December 1999     

Unsteady PTV velocity field past an airfoil solved with DNS: Part 2. Validation and application at Reynolds numbers up to <Emphasis Type="Italic">Re</Emphasis> ≤ 10<Superscript>4</Superscript>

Hybrid unsteady-flow simulation combining particle tracking velocimetry (PTV) and direct numerical simulation (DNS) is introduced in the series of two papers. Particle velocities on a laser-light sheet acquired with time-resolved PTV in a water tunnel are supplied to two-dimensional DNS with time intervals corresponding to the frame rate of the PTV. Hybrid velocity fields then approach those representing the PTV data in the course of time, and the reconstructed velocity fields satisfy the governing equations with the resolution comparable to numerical simulation. In part 2, by extending the capabilities of the hybrid simulation to higher Reynolds numbers, we simulate flows past the NACA0012 airfoil over ranges of Reynolds numbers (Re ≤ 10⁴) and angles of attack (−5° ≤ α ≤ 20°) and validate the proposed technique by comparing with experimental results in terms of the lift and drag coefficients. We also compare the results with unsteady Reynolds-averaged Navier–Stokes (URANS) simulation in two-dimensions and show the advantages of the hybrid simulation against two-dimensional URANS.     

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.     

基于涡粘性模型的翼型湍流流场熵产率计算

利用有限体积法实现了基于非正交同位网格的SIMPLE算法。基于熵分析方法,采用涡粘性模型求解湍流熵产方程,系统研究了湍流模型对二维翼型绕流流场熵产率的影响。通过计算NACA0012翼型在来流雷诺数为2.88×10⁶时,0°攻角~16.5°攻角范围内的翼型表面压力系数分布和升阻力特性,验证了算法及程序的正确性。结果表明,选择不同湍流模型时,翼型流场熵产的计算结果存在差异,湍流耗散是引起流场熵产的主要原因;翼型流场的熵产主要发生在翼型前缘区、壁面边界层和翼型尾流区域,流场熵产率与翼型阻力系数线性相关;当产生分离涡时,粘性耗散引起的熵产下降。     

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.     

The influence of free-stream disturbances on low Reynolds number airfoil experiments

The results of an investigation of the influence of free stream disturbances on the lift and drag performance of the Lissaman 7769 airfoil are presented. The wind tunnel disturbance environment is described using hot-wire anemometer and sound pressure level measurements. The disturbance level is increased by the addition of a ‘turbulence screen’ upstream of the test section and/or the addition of a flow restrictor downstream of the test section. For the Lissaman airfoil it was found that the problems associated with obtaining accurate wind tunnel data at low chord Reynolds numbers (i.e., below 200,000) are compounded by the extreme sensitivity of the boundary layers to the free stream disturbance environment. The effect of free stream disturbances varies with magnitude, frequency content, and source of the disturbance.     

Effects of Body Cross-sectional Shape on Flying Snake Aerodynamics

Despite their lack of appendages, flying snakes (genus Chrysopelea) exhibit aerodynamic performance that compares favorably to other animal gliders. We wished to determine which aspects of Chrysopelea’s unique shape contributed to its aerodynamic performance by testing physical models of Chrysopelea in a wind tunnel. We varied the relative body volume, edge sharpness, and backbone protrusion of the models. Chrysopelea’s gliding performance was surprisingly robust to most shape changes; the presence of a trailing-edge lip was the most significant factor in producing high lift forces. Lift to drag ratios of 2.7–2.9 were seen at angles of attack (α) from 10–30°. Stall did not occur until α > 30° and was gradual, with lift falling off slowly as drag increased. Chrysopelea actively undulates in an S-shape when gliding, such that posterior portions of the snake’s body lie in the wake of the more anterior portions. When two Chrysopelea body segment models were tested in tandem to produce a two dimensional approximation to this situation, the downstream model exhibited an increased lift-to-drag ratio (as much as 50% increase over a solitary model) at all horizontal gaps tested (3–7 chords) when located slightly below the upstream model and at all vertical staggers tested (±2 chords) at a gap of 7 chords.     

Experimental and numerical investigation on sound generation from airfoil-flow interaction

Aerodynamic noise due to interaction between incoming turbulence and rotating blades is an important component in the wind turbine noise. The rod-airfoil configuration is used to investigate the interactive phenomenon experimentally and numerically. Distribution of unsteady pressure on the airfoil surface is measured for different rod positions and airfoil attack angles. Two National Advisory Committee for Aeronatics (NACA) airfoils, NACA0012 and NACA0018, and two wind turbine airfoils, S809 and S825 are investigated. In addition, for low angles of attack, the flow field around the airfoil's leading edge is investigated with the particle image velocimetry (PIV). The experimental results indicate that unsteady pressure disturbances on the airfoil surface are related to the rod vortex shedding. Meanwhile, the interaction flow field of the rod and NACA0012 airfoil is simulated with the unsteady Reynolds averaged Navier-Stokes method (URANS), and the obtained pressure spectra are compared with the experimental results.     

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.     

PIV measurements of the asymmetric wake of a two dimensional heaving hydrofoil

Particle image velocimetry is used to examine the flow behind a two-dimensional heaving hydrofoil of NACA 0012 cross section, operating with heave amplitude to chord ratio of 0.215 at Strouhal numbers between 0.174 and 0.781 and a Reynolds number of 2,700. The measurements show that for Strouhal numbers larger than 0.434, the wake becomes deflected such that the average velocity profile is asymmetric about the mean heave position of the hydrofoil. The deflection angle of the wake, which is related to the average lift and drag on the hydrofoil, is found to lie between 13° and 18°. An examination of the swirl strength of the vortices generated by the hydrofoil motion reveal that the strongest vortices, which are created at the higher Strouhal numbers, dissipate most rapidly. This research article was submitted for the special issue on Animal locomotion: The hydrodynamics of swimming (Vol. 43, No. 5).     

Orthogonal grid generation for Navier–Stokes computations

Body conforming orthogonal grids were generated using a fast hyperbolic method for aerofoils, and were used to solve the Navier–Stokes equation in the generalized orthogonal system for the first time for time accurate simulation of incompressible flow. For grid generation, the Beltrami equation and the definition equation for the orthogonality are solved using a finite difference method. The grids generated around aerofoils by this method have better orthogonality than the results published by earlier investigators. The Navier–Stokes equation at Reynolds numbers of 3000 and 35 000 for NACA 0012 and NACA 0015 respectively, have been solved as an application. The obtained results match quite well with the corresponding experimental results. © 1998 John Wiley & Sons, Ltd.     

Anisotropic Organised Eddy Simulation for the prediction of non-equilibrium turbulent flows around bodies

The unsteady turbulent flow around bodies at high Reynolds number is predicted by an anisotropic eddy-viscosity model in the context of the Organised Eddy Simulation (OES). A tensorial eddy-viscosity concept is developed to reinforce turbulent stress anisotropy, that is a crucial characteristic of non-equilibrium turbulence in the near-region. The theoretical aspects of the modelling are investigated by means of a phase-averaged PIV in the flow around a circular cylinder at Reynolds number 1.4×10⁵. A pronounced stress–strain misalignment is quantified in the near-wake region of the detached flow, that is well captured by a tensorial eddy-viscosity concept. This is achieved by modelling the turbulence stress anisotropy tensor by its projection onto the principal matrices of the strain-rate tensor. Additional transport equations for the projection coefficients are derived from a second-order moment closure scheme. The modification of the turbulence length scale yielded by OES is used in the Detached Eddy Simulation hybrid approach. The detached turbulent flows around a NACA0012 airfoil (2-D) and a circular cylinder (3-D) are studied at Reynolds numbers 10⁵ and 1.4×10⁵, respectively. The results compared to experimental ones emphasise the predictive capabilities of the OES approach concerning the flow physics capture for turbulent unsteady flows around bodies at high Reynolds numbers.     

Supersonic leeside flow topology on delta wings revisited

The leeside vortex structures on delta wings with sharp leading edges were studied for supersonic flow at the Institute of Theoretical and Applied Mechanics of the Russian Academy of Sciences in Novosibirsk. The experiments were carried out with three wings with sweep angles of χ=68°, 73°, and 78° and parabolic profiles in the 0.6 × 0.6 m² test section of the blow-down wind tunnel T-313 of the institute. The test conditions were varied from Mach numbers M=2 to 4, unit Reynolds numbers from Re _l=26 × 10⁶ to 56 × 10⁶ m⁻¹, and angles of attack from α=0° to 22°. The results of the investigations revealed that for certain flow conditions shocks are formed above, below, and between the primary vortices. The experimental data were accurate enough to detect the onset of secondary and tertiary separation as well as other boundaries. The various flow regimes discussed in the literature were extended in several cases. The major findings are reported. Received: 6 September 1999/Accepted: 24 January 2000     

On the generation of large-scale homogeneous turbulence

An active turbulence generating grid, based on the rotating-vane design of Makita (1991), was developed for a large wind tunnel. At 2.14 m square, the grid is the largest of this type ever developed. To improve the isotropy of the turbulence generated, the grid was placed in the wind tunnel contraction. Measurements show that the grid produces a closely uniform mean flow and homogeneous isotropic turbulence to within two integral scales from the wall. By systematically varying the flow speed and parameters controlling the random motion of the vanes, grid turbulence with a wide variety of characteristics was produced and the dependence of those characteristics on the operating parameters of the grid revealed. Taylor Reynolds numbers of the grid turbulence varied from 100 to 1,360 and integral scales from 5 to almost 70 cm. The extreme cases represent some of the highest Reynolds number and largest scale homogeneous turbulent flows ever generated in a wind tunnel.     

A comparison of PIV measurements of canopy turbulence performed in the field and in a wind tunnel model

Particle image velocimetry (PIV) has been used to compare between turbulence characteristics just within and above a mature corn canopy and those of a model canopy setup in a wind tunnel (WT). The laboratory normalized mean velocity profile is adjusted using variable mesh screens to match the normalized mean shear of the corn field (CF) data. The smallest resolved scale in the field is about 15 times the Kolmogorov length scale (η_CF ≈ 0.4 mm), whereas in the WT it is 5 times η_WT (η_WT ≈ 0.15 mm). In both cases, the mean velocity and turbulence statistics are consistent with those measured using single point sensors. However, the profiles of normalized Reynolds shear stress in the field and the laboratory differ. Turbulent spectral densities calculated from PIV spatial and time series in the field display an inertial range spanning three decades. In the laboratory due to lower Reynolds numbers, the inertial range shrinks to two decades. Quadrant-Hole analysis is applied to Reynolds shear stress, vorticity magnitude and dissipation rates. In quadrants 1–3, the WT and field conditionally sampled stresses show similar trends. However, a conflicting trend is found in the sweep quadrant. The analysis confirms that sweep and ejections dominate the momentum flux and dissipation rate.The content of this paper, entitled “Applying PIV for Measuring Turbulence just within and above a Corn Canopy,” was presented at the 6th International Symposium on Particle Image Velocimetry at Pasadena, CA, USA, September 21–23, 2005.     

低雷诺数俯仰振荡翼型等离子体流动控制

针对低雷诺数翼型气动性能差的特点, 通过介质阻挡放电(dielectric barrier discharge, DBD)等离子体激励控制的方法, 提高翼型低雷诺数下的气动特性,改善其流场结构. 采用二维准直接数值模拟方法求解非定常不可压Navier-Stokes方程,对具有俯仰运动的NACA0012翼型的低雷诺数流动展开数值模拟.同时将介质阻挡放电激励对流动的作用以彻体力源项的形式加入Navier-Stokes方程,通过数值模拟探究稳态DBD等离子体激励对俯仰振荡NACA0012翼型气动特性和流场特性的影响.为了进行流动控制, 分别在上下表面的前缘和后缘处安装DBD等离子体激励器,并提出四种激励器的开环控制策略,通过对比研究了这些控制策略在不同雷诺数、不同减缩频率以及激励位置下的控制效果.通过流场结构和动态压强分析了等离子体进行流场控制的机理. 结果表明,前缘DBD控制中控制策略B(负攻角时开启上表面激励器,正攻角时开启下表面激励器)效果最好,后缘DBD控制中控制策略C(逆时针旋转时开启上表面激励器,顺时针旋转时开启下表面激励器)效果最好,前缘DBD控制效果会随着减缩频率的增大而下降, 同时会导致阻力增大.而后缘DBD控制可以减小压差阻力, 优于前缘DBD控制,对于计算的所有减缩频率(5.01~11.82)都有较好的增升减阻效果.在不同雷诺数下, DBD控制的增升效果较为稳定, 而减阻效果随着雷诺数的降低而变差,这是由流体黏性效应增强导致的.      

Tailoring turbulence with an active grid

Using an active grid in a wind tunnel, we generate homogeneous shear turbulence and initiate turbulent boundary layers with adjustable properties. Homogeneous shear turbulence is characterized by a constant gradient of the mean velocity and a constant turbulence intensity. It is the simplest anisotropic turbulent flow thinkable, and it is generated traditionally by equipping a wind tunnel with screens which have a varying transparency and flow straighteners. This is not done easily, and the reachable turbulence levels are modest. We describe a new technique for generating homogeneous shear turbulence using an active grid only. Our active grid consists of a grid of rods with attached vanes which can be rotated by servo motors. We control the grid by prescribing the time-dependent angle of each axis. We tune the vertical transparency profile of the grid by setting appropriate angles of each rod such as to generate a uniform velocity gradient, and set the rods in flapping motion around these angles to tailor the turbulence intensity. The Taylor Reynolds number reached was R _λ = 870, the shear rate S = ∂U/∂y = 9.2 s⁻¹, the nondimensional shear parameter S ^*≡ Sq ²/ε = 12 and u = 1.4 ms⁻¹. As a further application of this idea we demonstrate the generation of a simulated atmospheric boundary layer in a wind tunnel which has tunable properties. This method offers a great advantage over the traditional one, in which vortex-generating structures need to be placed in the wind tunnel to initiate a fat boundary layer.     

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.     

An experimental comparison of two adaptation strategies in an adaptive-walls wind-tunnel

In the present work an experimental comparison is made between two adaptation strategies: the Judd's method and the Everhart's method. A NACA 0012 airfoil has been tested at Mach numbers up to 0.4: models with chords up to 200 mm have been tested in a 200 mm × 200 mm test section. The two strategies, though based on different theoretical approaches, show a fairly good agreement as far as c _p distribution on the model, lift and drag curves and residual interference are concerned and agree, in terms of lift curve slope and drag coefficient at zero lift, with the McCroskey correlation.A version of this paper was presented at the 2nd International Conference on Experimental Fluid Mechanics, Torino, Italy, July 4–8, 1994.This work has been supported by the Italian Ministry for University and Scientific and Technological Research (MURST 60%).