舵几何特征对桨-舵系统尾流场演化的影响

发生在桨和舵之间的干扰会影响螺旋桨尾流的演化,导致尾流场中的湍流在下游增强,恶化船舶的振动和噪声性能,深入分析舵几何参数对桨-舵系统尾流场演化的影响能够为推进器尾流场的调节和减振降噪提供新思路.因此,从弦长、剖面和梯形舵入手分析不同的舵几何参数对螺旋桨尾流场演化特性的影响,使用大漩涡模拟方法模拟流场中的湍流结构,对不同舵弦长、剖面下的螺旋桨尾涡结构演化进行了分析,在舵弦长、剖面影响螺旋桨尾流场演化的研究的基础上分析了梯形舵对螺旋桨尾涡结构的影响,进一步分析了梯形舵影响下的螺旋桨尾流场中湍动能的分布.结果表明舵的弦长和剖面均会影响螺旋桨尾流场的演化,这种影响表现为更大的弦长和更厚的剖面会促进螺旋桨梢涡在舵压力面上的偏移,更薄的舵剖面会带来更强烈的螺旋桨毂涡偏移;涡管轮廓和舵表面脉动压力的对比均表明梯形舵会促进螺旋桨尾流场沿逆舵梯度方向偏移,从而导致螺旋桨的尾涡结构在舵两侧及下游呈现不对称分布,桨-舵系统下游的湍流结构与螺旋桨尾涡-舵碰撞过程、螺旋桨尾涡-舵随边涡干扰过程、螺旋桨梢涡-螺旋桨毂涡干扰有关,偏移更大的螺旋桨尾涡结构会在尾流场中更早地引起湍动能增强.     

Evolution of vortices in the wake of an ARJ21 airplane: Application of the lift-drag model

Wake separation is crucial to aircraft landing safety and is an important factor in airport operational efficiency. The near-ground evolution characteristics of wake vortices form the foundation of the wake separation system design. In this study, we analysed the near-ground evolution of vortices in the wake of a domestic aircraft ARJ21 initialised by the lift-drag model using large eddy simulations based on an adaptive mesh. Evolution of wake vortices formed by the main wing, flap and horizontal tail was discussed in detail. The horizontal tail vortices are the weakest and dissipate rapidly, whereas the flap vortices are the strongest and induce the tip vortex to merge with them. The horizontal tail and flap of an ARJ21 do not significantly influence the circulation evolution, height change and movement trajectory of the wake vortices. The far-field evolution of wake vortices can therefore be analysed using the conventional wake vortex model.     

螺旋桨梢涡不稳定性机理与演化模型研究

螺旋桨尾流场的涡流特性是一个基础但又十分复杂的流体力学问题, 它的复杂性源于其蕴含复杂的漩涡系统, 且该漩涡系统会在高速的剪切层流动中不断演化, 其流体动力学行为, 如由稳定态演变为不稳定态的机理以及复杂工况环境中的流动现象, 一直是流体力学领域的难点和备受关注的热点问题. 从工程应用的角度看, 桨后梢涡的演化特性与船舶结构物的宏观特性直接相关, 更好地理解多工况下螺旋桨尾流的动力学特性, 将有助于改善与振动、噪声以及结构问题等相关的推进器性能, 对综合性能优良的下一代螺旋桨的设计和优化有着重要的现实意义. 本文基于延迟分离涡模拟、大涡模拟和无湍流模型模拟方法以及粒子图像测速流场测试分别开展了螺旋桨尾流动力学特性的数值与试验研究, 对螺旋桨尾流不稳定性的触发机理进行了揭示. 基于均匀来流中螺旋桨梢涡的演化机理, 提出了螺旋桨梢涡演化模型. 该模型能够较为准确地模拟螺旋桨梢涡的演化过程, 预测螺旋桨梢涡融合的时间和位置, 对螺旋桨流噪声预报和控制以及性能优良的螺旋桨设计具有重要意义.      

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

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

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

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

单一水平轴风电机组尾迹的模拟方法与流动机理研究综述

为实现碳达峰、碳中和“3060”目标, 风能将在我国能源体系发挥重要作用. 风力机尾迹是影响风电性能和度电成本的关键因素, 需在风力机布置和控制设计中充分考虑. 本文首先介绍风力机尾迹的数值模拟方法, 包括解析模型、低阶模型、大涡模拟和来流湍流生成方法. 解析模型和低阶模型可快速计算风力机尾迹, 但依赖于模型参数, 且不能或不能准确预测尾迹湍流特性. 结合风力机参数化模型的大涡模拟可准确预测尾迹蜿蜒等湍流特征, 是流动机理研究的有力工具, 可为发展快速预测模型提供数据和理论支撑. 接着, 本文介绍了叶尖涡、中心涡和尾迹蜿蜒并讨论其产生机理. 对于湍流来流, 叶尖涡主要存在于近尾迹. 蜿蜒是远尾迹的主要特征, 影响下游风力机的来流特征. 尾迹蜿蜒的产生有两种机制: 来流大尺度涡和剪切层失稳. 数值和观测结果显示两种机制共同存在. 机舱和中心涡对尾迹蜿蜒有重要影响. 采用叶片和机舱的致动面模型可准确预测尾迹蜿蜒. 研究显示不同风力机尾迹间的湍流特征存在相似性, 为发展尾迹湍流的快速预测模型提供了理论依据. 当前研究多关注平坦地形上的风力机尾迹, 复杂地形和海洋环境下的大气湍流和风力机尾迹的机理复杂, 现有工程模型无法准确预测, 有待深入研究.      

Instantaneous and time-averaged flow fields of multiple vortices in the tip region of a ducted propulsor

The instantaneous and time-averaged flow fields in the tip region of a ducted marine propulsor are examined. In this flow, a primary tip-leakage vortex interacts with a secondary, co-rotating trailing edge vortex and other co- and counter-rotating vorticity found in the blade wake. Planar particle imaging velocimetry (PIV) is used to examine the flow in a plane approximately perpendicular to the mean axis of the primary vortex. An identification procedure is used to characterize multiple regions of compact vorticity in the flow fields as series of Gaussian vortices. Significant differences are found between the vortex properties from the time-averaged flow fields and the average vortex properties identified in the instantaneous flow fields. Variability in the vortical flow field results from spatial wandering of the vortices, correlated fluctuations of the vortex strength and core size, and both correlated and uncorrelated fluctuations in the relative positions of the vortices. This variability leads to pseudo-turbulent velocity fluctuations. Corrections for some of this variability are performed on the instantaneous flow fields. The resulting processed flow fields reveal a significant increase in flow variability in a region relatively far downstream of the blade trailing edge, a phenomenon that is masked through the process of simple averaging. This increased flow variability is also accompanied by the inception of discrete vortex cavitation bubbles, which is an unexpected result, since the mean flow pressures in the region of inception are much higher than the vapor pressure of the liquid. This suggests that unresolved fine-scale vortex interactions and stretching may be occurring in the region of increased flow variability.     

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.     

The effect of single-horn glaze ice on the vortex structures in the wake of a horizontal axis wind turbine

The present study experimentally investigated the effect of a simulated single-horn glaze ice accreted on rotor blades on the vortex structures in the wake of a horizontal axis wind turbine by using the stereoscopic particle image velocimetry(Stereo-PIV) technique. During the experiments, four horizontal axis wind turbine models were tested, and both "free-run" and "phase-locked" Stereo-PIV measurements were carried out. Based on the "free-run"measurements, it was found that because of the simulated single-horn glaze ice, the shape, vorticity, and trajectory of tip vortices were changed significantly, and less kinetic energy of the airflow could be harvested by the wind turbine. In addition, the "phase-locked" results indicated that the presence of simulated single-horn glaze ice resulted in a dramatic reduction of the vorticity peak of the tip vortices. Moreover,as the length of the glaze ice increased, both root and tip vortex gaps were found to increase accordingly.     

Coherent vorticity dynamics and dissipation in a utility-scale wind turbine wake with uniform inflow

The vorticity dynamics and its relationship to dissipation in the wake of a utility-scale wind turbine are investigated through large-eddy simulation. The vorticity dynamics is assessed through the enstrophy, which is related to the turbulent dissipation. The averaged enstrophy and turbulent dissipation are shown to be quantitatively similar in the wake. Using temporal phase averaging, the vorticity fluctuations are decomposed into coherent and random fluctuations with respect to the frequency of the tip vortices. The enstrophy in the tip vortices is dominated by coherent fluctuations, while the coherent fluctuations of root vortices are immediately saturated by the random vorticity fluctuations of the unstable hub vortex. The coherent strain rate has significant differences compared to the coherent enstrophy within one diameter downwind of blade tip, but the random enstrophy and strain rate are relatively similar. Differences in coherent enstrophy and strain rate decrease further from the rotor.     

Measurements of the aperiodic wake of a hovering rotor

 An experimental investigation was made to study the aperiodic flow characteristics of the tip vortices generated by one-bladed and two-bladed hovering rotors. Measurements of the tip vortex locations and accompanying aperiodicity statistics were established as a function of vortex age. Velocity field measurements were made using three-component laser Doppler velocimetry. The average amplitude of the aperiodicity was found to be a fraction of the measured viscous core radius, this being approximately 5% of blade chord or about 50% of the core radius for wake ages of less than two rotor revolutions. The aperiodicity appeared isotropic. A numerical analysis of the aperiodicity problem based on convolution with an assumed displacement probability function showed that for this experiment the measured tangential velocities in the tip vortices were underestimated by approximately 20% and the viscous core radii were overestimated by 20%. There was no evidence that the number of blades or the blade passage adversely influenced the aperiodicity of the rotor tip vortices. Received: 22 August 1997/Accepted: 4 February 1998     

Near-wake flow structure downwind of a wind turbine in a turbulent boundary layer

Wind turbines operate in the surface layer of the atmospheric boundary layer, where they are subjected to strong wind shear and relatively high turbulence levels. These incoming boundary layer flow characteristics are expected to affect the structure of wind turbine wakes. The near-wake region is characterized by a complex coupled vortex system (including helicoidal tip vortices), unsteadiness and strong turbulence heterogeneity. Limited information about the spatial distribution of turbulence in the near wake, the vortex behavior and their influence on the downwind development of the far wake hinders our capability to predict wind turbine power production and fatigue loads in wind farms. This calls for a better understanding of the spatial distribution of the 3D flow and coherent turbulence structures in the near wake. Systematic wind-tunnel experiments were designed and carried out to characterize the structure of the near-wake flow downwind of a model wind turbine placed in a neutral boundary layer flow. A horizontal-axis, three-blade wind turbine model, with a rotor diameter of 13 cm and the hub height at 10.5 cm, occupied the lowest one-third of the boundary layer. High-resolution particle image velocimetry (PIV) was used to measure velocities in multiple vertical stream-wise planes (x–z) and vertical span-wise planes (y–z). In particular, we identified localized regions of strong vorticity and swirling strength, which are the signature of helicoidal tip vortices. These vortices are most pronounced at the top-tip level and persist up to a distance of two to three rotor diameters downwind. The measurements also reveal strong flow rotation and a highly non-axisymmetric distribution of the mean flow and turbulence structure in the near wake. The results provide new insight into the physical mechanisms that govern the development of the near wake of a wind turbine immersed in a neutral boundary layer. They also serve as important data for the development and validation of numerical models.     

Mechanisms of wake deflection angle change behind a heaving airfoil

An immersed-boundary numerical method is applied to simulate the wake downstream of a two-dimensional heaving airfoil. A switch of vortex pattern is found to be the major reason that a deflected asymmetric wake reverses its deflection angle. Parameters of the heaving airfoil and flow that influence the onset and location of the vortex switching are discussed. While the previous literature deliberately discussed the wake deflection in the near wake region, this study shows that the deflection angle can change from the near wake to far wake regions. A cross-flow effective phase velocity is introduced to analyze the already-formed asymmetric wake behind the airfoil. A vortex dipole model and the related vortex dynamics analysis are developed to show that the change of the distance between the vortices is the primary factor that leads to the vortex pattern switching in the far wake.     

PIV study of near-field tip vortex behind perforated Gurney flaps

The impact of Gurney flaps, of different heights and perforations, on the growth and development of a tip vortex, both along the tip and in the near field of a finite NACA 0012 wing, at Re = 1.05 × 10⁵ was investigated by using particle image velocimetry (PIV). Wind-tunnel force balance measurements were also made to supplement the PIV results. This study is a continuation of the work of Lee and Ko (Exp Fluids 46(6):1005–1019, 2009) on the near-wake measurements behind perforated Gurney flaps. The present results show that along the tip, the overall behavior of the secondary vortices and their interaction with the primary, or tip, vortex remained basically unchanged, regardless of flap height and perforation. The peak vorticity of the tip vortex, however, increased with flap height and always exhibited a local maximum at x/c = 0.8 (from the leading edge). In the near field, the strength and structure of the near-field tip vortex were found to vary greatly with the flap height and perforation. The small flaps produced a more concentrated tip vortex with an increased circulation, while the large Gurney flaps caused a disruption of the tip vortex. The disrupted vortex can, however, be re-established by the addition of flap perforation. The larger the flap perforation the more organized the tip vortex. The Gurney flaps have the potential to serve as an alternative off-design wake vortex control device.     

Aperiodicity in the near field of full-scale rotor blade tip vortices

Blade tip vortices are the dominant vortical structures of the helicopter flow field. The inherent complexity of the vortex dynamics has led to an increasing interest in full-scale in situ experiments, where the near field, closely behind the blade, is of particular interest, since measures of vortex control mostly target this initial stage of development. To examine the near field, three-component particle image velocimetry (PIV) measurements of blade tip vortices of a full-scale helicopter in simulated hover flight in ground effect were conducted. A feasible and robust evaluation procedure was developed to minimise the shortcomings of full-scale PIV applications, such as a moderate spatial resolution and an elevated measurement noise level. At vortex ages ranging from y_v=1^°\psi_{\rm v}=1^{\circ} to 30°, a pronounced aperiodicity and asymmetry of the vortex were observed in -sections perpendicular to the vortex axes. At y_v=1^°\psi_{\rm v}=1^{\circ}, a preferential orientation of the vortex was observed. For increasing wake age, vortex wandering increased while the asymmetry of the vortex cores decreased. The high level of aperiodicity and core asymmetry must be taken into account when considering phase-averaged vortex characteristics in the near wake region.     

Symmetric vortex shedding in the near wake of a circular cylinder due to streamwise perturbations

Symmetric perturbations imposed on cylinder wakes may result in a modification of the vortex shedding mode from its natural antisymmetric, or alternating, to a symmetric one where twin vortices are simultaneously shed from both sides of the cylinder. In this paper, the symmetric mode in the wake of a circular cylinder is induced by periodic perturbations imposed on the in-flow velocity. The wake field is examined by PIV and LDV for Reynolds numbers about 1200 and for a range of perturbation frequencies between three and four times the natural shedding frequency of the unperturbed wake. In this range, a strong competition between symmetric and antisymmetric vortex shedding occurs for the perturbation amplitudes employed. The results show that symmetric formation of twin vortices occurs close to the cylinder synchronized with the oscillatory component of the flow. The symmetric mode rapidly breaks down and gives rise to an antisymmetric arrangement of vortex structures further downstream. The downstream wake may or may not be phase-locked to the imposed oscillation. The number of cycles for which the symmetric vortices persist in the near wake is a probabilistic function of the perturbation frequency and amplitude. Finally, it is shown that symmetric shedding is associated with positive energy transfer from the fluid to the cylinder due to the fluctuating drag.     

水下发射航行体尾涡不稳定性分析

在水下连续发射过程中前一发航行体尾流会对后一发航行体运动姿态稳定性产生流动干扰现象. 因此, 研究尾流中涡旋结构演变机理对解决多弹体水下连续发射流动干扰难题具有重要的意义. 本文采用改进型分离涡模型与能量方程, VOF多相流模型与重叠网格技术相结合方法, 对航行体水下发射尾流演变过程开展精细化模拟研究, 其中模拟结果和实验吻合度较好, 验证了本文数值方法的有效性. 以航行体尾流区域为重点研究对象, 分析了尾流区瞬态流场分布, 讨论了横流强度和雷诺数对尾涡结构演变以及脉动压力分布特性的影响. 结果表明: 由于尾流区高速流体核心区与低速自由流相互作用导致Kelvin-Helmholtz不稳定现象出现, 可以清晰地发现涡旋结构在剪切力的作用下发生脱落. 在横流条件下, 航行体尾端脱落的涡环与涡腿形成发卡涡, 而多个发卡涡沿轴向间隔排列组成发卡涡包存在于尾流中. 随着横流强度增大, 形成多级发卡涡包结构, 而导致脉动压力二次峰值均出现的主要原因是尾流涡旋流场演变引起的. 随着雷诺数的增大, 尾流中由圆柱形涡和U型涡组成的二次涡结构逐渐明显, 不稳定性加强.      

Flow above the free end of a surface-mounted finite-height circular cylinder: A review

The wake of a surface-mounted finite-height circular cylinder and the associated vortex patterns are strongly dependent on the cylinder aspect ratio and the thickness of the boundary layer on the ground plane relative to the dimensions of the cylinder. Above a critical aspect ratio, the mean wake is characterized by streamwise tip vortex structures and Kármán vortex shedding from the sides of the cylinder. Below a critical aspect ratio, a unique mean wake structure is observed. Recent experimental studies in the literature that used phase-averaged techniques, as well as recent numerical simulations, have led to an improved physical understanding of the near-wake vortex flow patterns. However, the flow above the free end of the finite circular cylinder, and its relationship to the near wake, has not been systematically studied. The effects of aspect ratio and boundary layer thickness on the free-end flow field are also not completely understood, nor has the influence of Reynolds number on the free-end flow field been fully explored. Common features associated with the free end include separation from the leading edge, a mean recirculation zone containing a prominent cross-stream arch (or mushroom) vortex, and reattachment onto the free-surface. Other flow features that remain to be clarified include a separation bubble near the leading edge, one or two cross-stream vortices within this separation bubble, the origins of the streamwise tip or trailing vortices, and various critical points in the near-surface flow topology. This paper reviews the current understanding of the flow above the free end of a surface-mounted finite-height circular cylinder, with a focus on models of the flow field, surface oil flow visualization studies, pressure and heat flux distributions on the free-end surface, measurements of the local velocity field, and numerical simulations, found in the literature.     

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.     

Effect of perforation on flow past a conic cylinder at <Emphasis Type="Italic">Re</Emphasis> = 100: vortex-shedding pattern and force history

The flow past a circular-section cylinder with a conic shroud perforated with four holes at the peak was simulated numerically at \(Re=100\), considering two factors, viz. the angle of attack and the diameter of the holes. The effects of the perforated conic shroud on the vortex shedding pattern in the near wake was mainly investigated, as well as the time history of the drag and lift forces. In the investigated parameter space, three flow regimes were generally identified, corresponding to weak, moderate, and strong disturbance effects. In regime I, the wake can mainly be described by alternately shedding Kármán or Kármán-like vortices. In regime II, the spanwise vortices are obviously disturbed along the span due to the appearance of additional vorticity components and their interactions with the spanwise vortices, but still shed in synchronization along the spanwise direction. In regime III, the typical Kármán vortices partially or totally disappear, and some new vortex shedding patterns appear, such as \(\Omega \)-type, obliquely shedding, and crossed spanwise vortices with opposite sign. Corresponding to these complex vortex shedding patterns in the near wake, the fluid forces no longer oscillate regularly at a single vortex shedding frequency, but rather with a lower modulation frequency and multiple amplitudes. An overview of these flow regimes is presented.