Asymmetric vortex shedding flow past an inclined flat plate at high incidence

This paper reports an experimental investigation of the vortex shedding wake behind a long flat plate inclined at a small angle of attack to a main flow stream. Detailed velocity fields are obtained with particle-image velocimetry (PIV) at successive phases in a vortex shedding cycle at three angles of attack, α=20°, 25° and 30°, at a Reynolds number Re≈5,300. Coherent patterns and dynamics of the vortices in the wake are revealed by the phase-averaged PIV vectors and derived turbulent properties. A vortex street pattern comprising a train of leading edge vortices alternating with a train of trailing edge vortices is found in the wake. The trailing edge vortex is shed directly from the sharp trailing edge while there are evidences that the formation and shedding of the leading edge vortex involve a more complicated mechanism. The leading edge vortex seems to be shed into the wake from an axial location near the trailing edge. After shedding, the vortices are convected downstream in the wake with a convection speed roughly equal to 0.8 the free-stream velocity. On reaching the same axial location, the trailing edge vortex, as compared to the leading edge vortex, is found to possess a higher peak vorticity level at its centre and induce more intense fluid circulation and Reynolds stresses production around it. It is found that the results at the three angles of attack can be collapsed into similar trends by using the projected plate width as the characteristic length of the flow.     

基于涡旋识别方法的高速列车尾涡结构的讨论

利用改进型延迟分离涡模拟方法对缩尺比例1:30的高速列车简化模型的绕流流场进行数值计算,主要针对近尾流区的涡旋结构展开具体讨论. 通过不同的涡旋识别方法,发现在尾涡结构中,高涡量的强涡旋主要聚集于尾车附近,而涡量较低但处于相对稳定状态的涡旋分布在大部分尾流空间中. 对此,主要基于最新提出的涡旋定义及其物理意义认为,由于边界层在尾部发生的流动分离,剪切变形以及高涡量的扩散对强涡旋的形成发挥着重要的作用,而涡旋会被较强的剪切旋转拉伸,使得局部复杂的流动表现出突出的湍流特性;另一方面,尽管涡强度明显下降,但是在强剪切应变迅速衰减的情况下,流向涡核中的涡旋涡量是主要的,此时,在较接近地面的情况下,流体微团以涡核为中心的旋转运动使得涡旋与地面之间的相互作用成为主导的流动机制. 虽然涡强度会相对缓慢地衰减,但是从湍流能量产生的角度,该机制对涡旋的自维持发挥重要的作用,从而使尾涡结构能够相对稳定地存在于尾流流动中.      

Unsteady RANS and detached-eddy simulations of flow around a circular cylinder in ground effect

Unsteady Reynolds-averaged Navier–Stokes (URANS) simulations and detached-eddy simulations (DES) were performed of flow around a circular cylinder placed near and parallel to a moving ground, on which substantially no boundary layer developed to interfere with the cylinder. The results were compared with experiments previously reported by the authors to examine how accurately the URANS and DES can predict the cessation of von Kármán-type vortex shedding and the attendant critical drag reduction of the cylinder in ground effect. The DES, which were performed in a three-dimensional domain with spanwise periodicity imposed, correctly captured the cessation of the vortex shedding, whereas both two- and three-dimensional URANS also predicted it but at a much smaller gap-to-diameter ratio compared with the experiments. The wake structures of the cylinder predicted by the DES were in good agreement with the experiments in both large- and small-gap regimes, and also in the intermediate-gap regime, where the DES captured the intermittence of the vortex shedding in the near-wake region. Based on the results obtained, further discussions are also given to the reason why the von Kármán-type vortices in the URANS solutions incorrectly ‘survived’ until the cylinder came much closer to the ground.     

长编组高速列车的列车风动模型实验研究

列车风是高速列车运行时诱导产生的气流流动, 是列车空气动力学重要的研究内容和保证列车运行安全的重要方面. 本文利用缩尺比例为1\bh8的八编组高速列车模型进行了列车风的动模型实验,测试了明线运行状态下列车周围的流动参数,突破了短编组列车风动模型测试所带来的局限性.列车风系综平均曲线和标准差曲线说明: 列车头部会引起稳定的列车风,在车身和车尾处的列车风具有非常明显的非定常特征.列车风所反映的车身周围的气流扰动在第二节车厢开始显现, 列车风速振荡上升,在第七节车厢达到局部最大值.转向架舱和车厢间风挡间隙的气流干扰并没有在列车风曲线上表现出来.利用本征正交分解法分析列车风尾迹区的实验结果,发现列车风的扰动能量集中于近尾迹区, 次之是车身发展区.以各次实验结果中列车风的峰值位置距离车尾远近为条件,对列车风实验结果进行条件平均分析,表明列车尾涡生成时与列车风探针间的相对位置关系会影响列车风尾迹区的峰值形态.      

Influences of marshalling length on the flow structure of a maglev train

In this paper, transient numerical simulations of maglev trains of different marshalling lengths (2, 4, and 8-car group trains) were conducted in the open air and without wind. This was done by solving the three-dimensional incompressible Navier-Stokes equations using an SST K-ω double-equation IDDES turbulence model. The results were compared with the results of wind tunnel experiments to verify the feasibility of numerical simulation. The results show an increase in the marshalling lengths of the train affects the flow above and below the train. With the increase of the marshalling length, the position of the flow separation in the tail car is advanced. The turbulence generated by the average shear on the x, y-plane and the x, z-plane as a component of the turbulence of the wake region increases. The region that produces non-vorticial vorticity in the main vortex becomes narrow and moves towards tail car. The structural analysis of the wake indicates that the wake structure of the 8-car group train is quite different from the other two groups. Both the time-averaged slipstream and the gust analysis show that the maximum expected slipstream velocity at the track-side increases as the train marshalling length increases. At the platform height, the change in vertical position of the wake vortex structure of the 2, 4, and 8-car group trains caused the difference of the shear flow regions. This is why the maximum expected slipstream velocity generated by the 4-car group train at this position is largest. As the marshalling length of the train increases, the time-average drag and lift force coefficient of the tail car have a significant negative correlation.     

一种风力机气动计算的全自由涡尾迹模型

采用全自由方式建立风力机尾流场的涡尾迹模型,引入“虚拟周期”的概念,并发展一种自适应松弛因子方法,从而改善了自由尾迹迭代的稳定性,提高了迭代收敛速度。利用建立的自由涡尾迹模型,计算了风力机叶片的尾流场结构、气动性能及叶片载荷,并与实验结果进行了对比分析。结果表明,尖速比越大,自适应松弛因子方法对缩小模型计算时间越有效;全自由涡尾迹模型能准确给出风力机尾流场的结构,包括尾迹的扩张以及叶尖涡和叶根涡的产生、发展和耗散的过程,风轮扭矩与实验数据吻合;叶片载荷分布的计算结果在低风速下与实验值基本一致,但是在大风速下差别较大,说明需要一个准确的失速模型。     

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.     

高速列车空气动力学研究技术综述

随着我国高速铁路的快速发展,高速列车运行速度越来越快,包括气动阻力、横风效应、会车效应、隧道效应和气动噪声等一系列空气动力学问题日益突出. 利用模型试验、实车测量和数值计算等不同的研究手段,开展全面的高速列车空气动力学研究显得十分重要. 本文比较全面和系统地介绍了国内外高速列车空气动力学研究在模型试验、实车测量、数值计算等方面的技术现状和进展情况,对今后的发展方向和内容进行了展望.     

Detached eddy simulation of flow around two wall-mounted cubes in tandem

We investigate the performance of unsteady Reynolds-averaged Navier–Stokes (URANS) computation and various versions of detached eddy simulation (DES) in resolving coherent structures in turbulent flow around two cubes mounted in tandem on a flat plate at Reynolds number (Re) of 22,000 and for a thin incoming boundary layer. Calculations are carried out using four different coherent structure resolving turbulence models: (1) URANS with the Spalart–Allmaras model; (2) the standard DES [Spalart, P.R., Jou, W.H., Strelets, M., Allmaras, S.R., 1997. Comments on the feasibility of LES for wings, and on a hybrid RANS/LES approach. In: Liu, C., Liu, Z., (Eds.), Advances in DNS/LES. Greyden Press, Columbus, OH]; (3) the Delayed DES (DDES); and (4) the DES with a low-Re modification (DES-LR) [Spalart, P., Deck, S., Shur, M., Squires, K., Strelets, M., Travin, A., 2006. A new version of detached eddy simulation, resistant to ambiguous grid densities. Theor. Comput. Fluid Dyn. 20 (3), 181–195]. The grid sensitivity of the computed solutions is examined by carrying out simulations on two successively refined grids. The computed results for all cases are compared with the experimental measurements of Martinuzzi and Havel [Martinuzzi, R., Havel, B., 2000. Turbulent flow around two interfering surface-mounted cubic obstacles in tandem arrangement. ASME J. Fluids Eng. 122, 24–31] for two different cube spacings. All turbulence models reproduce essentially identical separation of the approach thin boundary layer and yield an unsteady horseshoe vortex system consisting of multiple vortices in the leading edge region of the upstream cube. Significant discrepancies between the URANS and all DES solutions are observed, however, in other regions of interest such as the shear layers emanating from the cubes, the inter-cube gap and the downstream wake. Regardless of the grid refinement, URANS fails to capture key features of the mean flow, including the second horseshoe vortex in the upstream junction and recirculating flow on the top surface of the downstream cube for the large cube spacing, and underestimates significantly turbulence statistics in most regions of the flow for both cases. On the coarse mesh, all three DES approaches appear to yield very similar results and fail to reproduce the second horseshoe vortex. The standard DES and DDES solutions obtained on the fine meshes are essentially identical and both suffer from premature switching to unresolved DNS, due to the mis-interpretation of grid refinement as wall proximity, which leads to spurious vortices in the inter-cube region. Numerical solutions show that the low-Re modification (DES-LR) is critical prerequisite in DES on the ambiguously fine – not fine enough for full LES – mesh to prevent excessive nonlinear drop of the subgrid eddy viscosity in low cell-Re regions like in the inter-obstacle gap. Mean flow quantities and turbulence statistics obtained with DES-LR on the fine mesh are in good overall agreement with the measurements in most regions of interest for both cases.     

Coherent structure in the turbulent wake behind a circular cylinder

A wake behind a circular cylinder at Reynolds number 850–1700 was visualized by the smoke-wire method. The observations of the How together with the results of quantitative measurements, such as various velocity correlation coefficients, illustrated the formation process of spoon-shaped large eddies in the region 90 ⩽x/d⩽ 230 attained through the deformation and rearrangement of the regular Karman vortices. A spoon vortex was likely to pair with the counterpart on the opposite side of the wake. The large-scale bulges of the turbulent and non-turbulent interface of the wake were shown to correspond to these spoon vortices.These results indicate that some coherent structures are organized by rearrangement and deformation of initially regular vortices in turbulent flow.     

Flow topology and unsteady features of the wake of a generic high-speed train

The unsteady wake of a high-speed train is investigated experimentally. From a practical point of view, the wake region is of considerable importance as it is where slipstream velocities—velocities induced by the vehicles movement through air—are largest. In turn, this can create a considerable risk for passengers and track-side workers as the train passes. The flow is quantified in a 1:10 scale wind-tunnel experiment using high-frequency 4-hole dynamic pressure cobra probes, surface-pressure measurements and flow visualisation. The dominant feature of the time-average wake topology consists of a clearly identifiable counter-rotating streamwise vortex pair. Although the wake structure and evolution should perhaps be considered as a whole, the near wake exhibits periodic unsteadiness, at a Strouhal number of 0.2, that could be attributed to periodic shedding from the sides and to a lesser extent the top surface. This periodicity feeds into the trailing vortices, consistent with lateral and vertical displacement of the cores as they advect downstream and thus affecting maximum slipstream velocities.     

Secondary wake instabilities of a blunt trailing edge profiled body as a basis for flow control

Flow in the wake of a blunt trailing edge profiled body, comprised of an elliptical leading edge and a rectangular trailing edge, has been investigated experimentally, to identify and characterize the secondary instabilities accompanying the von Kármán vortices. The experiments, which involve laser-induced fluorescence for visualization and particle image velocimetry for quantitative measurement of the wake instabilities, cover Reynolds numbers ranging from 250 to 2,150 based on thickness of the body, to include the wake transition regime. The dominant secondary instability appears as spanwise undulations in von Kármán vortices, which evolve into pairs of counter-rotating vortices, with features resembling the instability mechanism predicted by Ryan et al. (J Fluid Mech 538:1–29, 2005). Feasibility of a flow control approach based on interaction with the secondary instability using a series of discrete trailing edge injectors has also been investigated. The control approach mitigates the adverse effects of vortex shedding in certain conditions, where it is able to amplify the secondary instability effectively.     

Vortex shedding flow behind a slowly rotating circular cylinder

This paper investigates flow past a rotating circular cylinder at 3600?Re?5000 and α?2.5. The flow parameter α is the circumferential speed at the cylinder surface normalized by the free-stream velocity of the uniform cross-flow. With particle image velocimetry (PIV), vortex shedding from the cylinder is clearly observed at α<1.9. The vortex pattern is very similar to the vortex street behind a stationary circular cylinder; but with increasing cylinder rotation speed, the wake is observed to become increasing narrower and deflected sideways. Properties of large-scale vortices developed from the shear layers and shed into the wake are investigated with the vorticity field derived from the PIV data. The vortex formation length is found to decrease with increasing α. This leads to a slow increase in vortex shedding frequency with α. At α=0.65, vortex shedding is found to synchronize with cylinder rotation, with one vortex being shed every rotation cycle of the cylinder. Vortex dynamics are studied at this value of α with the phase-locked eduction technique. It is found that although the shear layers at two different sides of the cylinder possess unequal vorticity levels, alternating vortices subsequently shed from the cylinder to join the two trains of vortices in the vortex street pattern exhibit very little difference in vortex strength.     

Nonlinear characteristics of a thin wing in a shock-free flow around the nose

An effective means of controlling wing leading-edge stall at high angles of attack is deflection of the nose in order to assure shock-free entrance of the stream. A numerical method of computing the angles of nose deflection and the aerodynamic characteristics of a thin wing of arbitrary planform for a shock-free entrance of the steady ideal incompressible fluid stream is elucidated in this paper on the basis of nonlinear wing theory [1]. The problem is solved by the method of discrete vortices. In the computations, the wing and its wake, replaced by a vortex sheet, are modeled by a system of discrete vortices which are nonlinear segments with constant circulation along the length. The angles of deflection of the nose and the aerodynamic characteristics of the wing, including shunting of the free vortices shed from the side and trailing edges, are determined during the computation. Examples of an electronic digital computer are presented.     

The Vortical Structures in the Rear Separation and Wake Produced by a Supersonic Micro-Ramp

The vortical structures in the rear separation and wake region produced by a micro-ramp that immersed in a supersonic turbulent boundary layer are investigated. The small scale separation close to the trailing edge was revealed and this confirms the previous experimental observation. Between the reverse region and surrounding fast moving flow, a three-dimensional shear layer was formed, and vortices are generated. By using vortex line method, the spiral points were illustrated as the cross-sections of the Ω-shaped vortices that follow the shape of the separation. The vortical structure was analogous to that in the wake region, where similar Ω-shaped vortex which follows the deficit region caused by the micro-ramp can be observed. Finally, the revealed flow topology was conceived new and beneficial to the studying of wall bounded turbulence which involves similar vortical structures but in a smaller scale, compared with the vortical pattern in the current micro-ramp wake.     

Three-component velocity field measurements of propeller wake using a stereoscopic PIV technique

A stereoscopic PIV (Particle Image Velocimetry) technique was used to measure the three-dimensional flow structure of the turbulent wake behind a marine propeller with five blades. The out-of-plane velocity component was determined using two CCD cameras with an angular displacement configuration. Four hundred instantaneous velocity fields were measured for each of four different blade phases, and ensemble averaged in order to find the spatial evolution of the propeller wake in the region from the trailing edge up to one propeller diameter (D) downstream. The influence of propeller loading conditions on the wake structure was also investigated by measuring the velocity fields at three advance ratios (J=0.59, 0.72 and 0.88). The phase-averaged velocity fields revealed that a viscous wake formed by the boundary layers developed along the blade surfaces. Tip vortices were generated periodically and the slipstream contracted in the near-wake region. The out-of-plane velocity component and strain rate had large values at the locations of the tip and trailing vortices. As the flow moved downstream, the turbulence intensity, the strength of the tip vortices, and the magnitude of the out-of-plane velocity component at trailing vortices all decreased due to effects such as viscous dissipation, turbulence diffusion, and blade-to-blade interaction.     

The wake flow downstream of a propeller-rudder system

We report wall-resolved, large-eddy simulations for the case of a propeller operating upstream of a hydrofoil, mimicking a rudder. Our primary objective is the identification of wake features that are unique to this coupled system, when compared to open-water cases, which are usually the focus of experiments and computations in the literature. We were able to achieve unprecedented levels of numerical resolution, which capture the dynamics of all energetic eddies in the flow by using a scalable, conservative, structured solver in cylindrical coordinates. The boundary conditions on the rotating propeller and hydrofoil were enforced via an immersed boundary formulation. The largest values of turbulent stresses in the wake of the hydrofoil are achieved outwards from the radial coordinate of the tip of the propeller blades. This is due to spanwise gradients across the hydrofoil (in the direction parallel to the span of the hydrofoil), producing a displacement of the pressure side legs of the tip vortices towards outer coordinates, where they experience shear with the wake of the hydrofoil. The evolution of turbulence is non-monotonic across the streamwise direction. This is a consequence of the growing shear resulting from the complex interactions involving the shear layers from the trailing edge, the tip vortices and the two branches of the hub vortex coming from the two sides of the hydrofoil. Such a shear is reinforced by the spanwise velocities developed by the two branches of the propeller wake across the hydrofoil. Compared to an isolated propeller, these phenomena enhance turbulence production. The present results highlight that a downstream hydrofoil, typical of surface ships, is able to significantly intensify the wake signature of a propeller.     

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.     

A nonlinear model for aerodynamic configuration of wake behind horizontal-axis wind turbine

Determination of the aerodynamic configuration of wake is the key to analysis and evaluation of the rotor aerodynamic characteristics of a horizontal-axis wind turbine.According to the aerodynamic configuration, the real magnitude and direction of the onflow velocity at the rotor blade can be determined, and subsequently, the aerodynamic force on the rotor can be determined. The commonly employed wake aerodynamic models are of the cylindrical form instead of the actual expanding one. This is because the influence of the radial component of the induced velocity on the wake configuration is neglected. Therefore, this model should be called a "linear model". Using this model means that the induced velocities at the rotor blades and aerodynamic loads on them would be inexact. An approximately accurate approach is proposed in this paper to determine the so-called "nonlinear" wake aerodynamic configuration by means of the potential theory,where the influence of all three coordinate components of the induced velocity on wake aerodynamic configuration is taken into account to obtain a kind of expanding wake that approximately looks like an actual one. First, the rotor aerodynamic model composed of axial(central), bound, and trailing vortexes is established with the help of the finite aspect wing theory. Then, the Biot-Savart formula for the potential flow theory is used to derive a set of integral equations to evaluate the three components of the induced velocity at any point within the wake. The numerical solution to the integral equations is found,and the loci of all elementary trailing vortex filaments behind the rotor are determined thereafter. Finally, to formulate an actual wind turbine rotor, using the nonlinear wake model, the induced velocity everywhere in the wake, especially that at the rotor blade,is obtained in the case of various tip speed ratios and compared with the wake boundary in a neutral atmospheric boundary layer. Hereby, some useful and referential conclusions are offered for the aerodynamic computation and design of the rotor of the horizontal-axis wind turbine.     

Numerical simulation of flapping‐wing insect hovering flight at unsteady flow

A computational fluid dynamics (CFD) analysis was conducted to study the unsteady aerodynamics of a virtual flying bumblebee during hovering flight. The integrated geometry of bumblebee was established to define the shape of a three‐dimensional virtual bumblebee model with beating its wings, accurately mimicking the three‐dimensional movements of wings during hovering flight. The kinematics data of wings documented from the measurement to the bumblebee in normal hovering flight aided by the high‐speed video. The Navier–Stokes equations are solved numerically. The solution provides the flow and pressure fields, from which the aerodynamic forces and vorticity wake structure are obtained. Insights into the unsteady aerodynamic force generation process are gained from the force and flow‐structure information. The CFD analysis has established an overall understanding of the viscous and unsteady flow around the virtual flying bumblebee and of the time course of instantaneous force production, which reveals that hovering flight is dominated by the unsteady aerodynamics of both the instantaneous dynamics and also the past history of the wing. A coherent leading‐edge vortex with axial flow and the attached wingtip vortex and trailing edge vortex were detected. The leading edge vortex, wing tip vortex and trailing edge vortex, which caused by the pressure difference between the upper and the lower surface of wings. The axial flow, which include the spanwise flow and chordwise flow, is derived from the spanwise pressure gradient and chordwise pressure gradient, will stabilize the vortex and gives it a characteristic spiral conical shape. Copyright © 2006 John Wiley & Sons, Ltd.