基于层析PIV的椭圆水翼近尾迹梢涡实验研究

赵航, 佘文轩, 高琪, 邵雪明

赵航,佘文轩,高琪,等. 基于层析PIV的椭圆水翼近尾迹梢涡实验研究[J]. 实验流体力学,2022,36(2):82-91. DOI: 10.11729/syltlx20210108
引用本文: 赵航,佘文轩,高琪,等. 基于层析PIV的椭圆水翼近尾迹梢涡实验研究[J]. 实验流体力学,2022,36(2):82-91. DOI: 10.11729/syltlx20210108
ZHAO H,SHE W X,GAO Q,et al. TPIV study for near-field tip vortex from an elliptical hydrofoil[J]. Journal of Experiments in Fluid Mechanics, 2022,36(2):82-91.. DOI: 10.11729/syltlx20210108
Citation: ZHAO H,SHE W X,GAO Q,et al. TPIV study for near-field tip vortex from an elliptical hydrofoil[J]. Journal of Experiments in Fluid Mechanics, 2022,36(2):82-91.. DOI: 10.11729/syltlx20210108

基于层析PIV的椭圆水翼近尾迹梢涡实验研究

基金项目: 国家自然科学基金(91852204);国家重点研发计划资助(2020YFA0405700)
详细信息
    作者简介:

    赵航: (1994—),男,湖南安乡人,博士研究生。研究方向:实验流体力学流动测量。通信地址:浙江省杭州市西湖区浙江大学玉泉校区航空航天学院(310027)。E-mail:hangzhao@zju.edu.cn

    通讯作者:

    高琪: E-mail:qigao@zju.edu.cn

  • 中图分类号: O352

TPIV study for near-field tip vortex from an elliptical hydrofoil

  • 摘要: 梢涡空化作为一种常见的空化现象,广泛存在于水力机械及船舶推进领域。梢涡空化初生与桨叶梢部的旋涡流动密切相关,因此有必要深入研究梢涡流场,揭示其流动特征与空化的内在联系。基于高时间解析度的层析PIV技术,在高速空泡水洞中对椭圆水翼的近尾迹梢涡流场开展了实验研究。结果表明:梢涡在近尾迹区域内存在明显的摆动现象,未考虑旋涡摆动的时间平均会在时均流场中引入额外的误差,因此在梢涡特性的定量研究中有必要滤除旋涡摆动的影响;在水翼脱落剪切层的作用下,涡核中心两侧的切向速度分布明显不对称,且在剪切层与涡核之间存在高速轴向流动区域;梢涡流场中的湍流脉动能量主要集中在涡核内部,且由法向、展向速度脉动主导。结合前人研究,发现法向、展向速度脉动是涡核内部湍流压力脉动的主要来源。
    Abstract: Tip vortex cavitation (TVC) is a common type of cavitation in hydraulic machinery and marine propulsion. Since TVC inception is highly relevant to the vortical flow around the blade tip of turbines and propellers, it is essential to give more insights into the flow field of the tip vortex to reveal the inherent relationship between flow properties and TVC. Measurement for the tip vortex from an elliptical hydrofoil has been conducted in a high-speed cavitation tunnel utilizing tomographic particle image velocimetry (TPIV) with high time-resolution. The results show that the wandering motion of the tip vortex is noticeable in the near field. The time-averaging process without taking into account wandering motion can bring extra errors into the time-averaged flow field. Therefore, it is necessary to filter out the wandering motion for the quantitative analysis on vortex characteristics. The tip vortex is under roll-up process and can be greatly affected by the shear layer from the hydrofoil, which contributes to the asymmetric circumferential velocity distribution and a high-axial-velocity area between the shear layer sheet and the vortex core. The tip vortex contains the most of the turbulence energy within its core and the turbulence energy is dominated by the vertical and spanwise velocity fluctuations which are considered as the main source of the fluctuating pressure in the core center combining with previous researches.
  • 旋涡作为一种常见的流动现象广泛存在于自然界以及近现代工程中,被喻为流体运动的“肌腱”[1]。在水动力学领域,旋涡身影随处可见,如轴流泵等水力机械中的梢隙涡[2]和船舶螺旋桨的毂涡、梢涡[3]等。尤其是桨叶梢涡常诱导产生空化现象,即梢涡空化。梢涡空化现象会引起较大的脉动噪声及结构振动,阻碍了水力机械和船舶推进技术的发展[4-5]。因此,近几十年来,国内外学者对梢涡空化问题开展了大量的研究。

    在梢涡空化研究中,常使用空化数这一基本物理参数。该参数与流场中的压力高度相关,因而早期对梢涡空化的研究主要以基于压力的空化预测为主。早在20世纪60年代,McCormick[6]针对梢涡空化问题开展了系统性的研究,基于理论分析,提出以旋涡中心最小压力系数来预测梢涡的临界空化数,推导得出的最小压力系数能够较好地预测梢涡空化的趋势,但其具体数值与实验中实际观察到的临界空化数有较大差异。此后,有相当一部分研究基于实验获取的梢涡特性(如旋涡强度、涡核半径等),结合理论涡模型推导旋涡中心的压力系数,并以此建立空化预测模型[7-9]。基于简化的径向N-S方程则是另一种获取梢涡内部压力分布的方法[10-11],但该方法存在一定误导性,其前提假设为轴对称涡,而梢涡空化常发生于近尾迹区域,该区域内的旋涡处于发生发展阶段并且高度不对称。基于前人相关工作,Arndt[12]总结完善了梢涡初生空化数预测模型,该预测模型主要考虑了3类影响因素,分别为旋涡中心平均最小压力、当地湍流压力脉动以及水质影响。Zhang等[13]则在此基础上进一步强调了桨叶梢部的旋涡流动会显著影响预测模型中的相关参数,且局部的湍流压力脉动可能由某些流动诱导产生,因此随着流场测量技术和数值模拟方法的发展,有相当一部分研究工作聚焦于桨叶梢涡的流动结构。

    近年来,研究人员采用粒子图像测速技术(Particle Image Velocimetry,PIV)对梢涡流场开展了大量的实验研究,重点关注有无空化情况下的涡流场差异[14-16]。这部分研究工作极具开创性,重点揭示了处于空化状态下的梢涡流场特性:空泡附近的液相涡流场会受到阻滞作用,具体表现为空泡附近的旋涡切向速度减小,旋涡强度减弱。Bosschers[17]在基于空化Lamb-Oseen涡模型的数值模拟中也观察到了类似现象,其结果与实验吻合较好。在全湿流以及空化流状态下,基于大涡模拟(Large Eddy Simu-lation,LES)的梢涡流动预测也取得了不错的进展,相应模型能够有效地预测梢涡流动状态[18]。此外,还有部分工作针对抑制梢涡空化展开,采用流动控制的手段,如主动射流[19]、梢部柔性丝线[20-21]、翼梢小翼[22]、粗糙翼面[23-24]等,以减弱梢涡强度,从而达到推迟或抑制梢涡空化的目的,这方面研究工作同样也与梢涡流动密切相关。

    上述研究主要集中于梢涡时均二维流场结构,只有极少部分数值模拟涉及梢涡三维流场。迄今为止,关于梢涡三维流场的实验研究仍然十分欠缺,开展这方面工作不仅有利于从三维流动的角度进一步阐明梢涡空化的内在机理,准确描述梢涡的时空发展特性,还能在工程应用上指导船舶螺旋桨叶设计及梢涡空化控制策略,同时也可作为数值模拟及梢涡空化预测模型的重要参考及验证。因此,本文将从这一角度出发研究梢涡流动特性。

    本文基于高时间解析度的层析PIV技术,在高速空泡水洞中对椭圆水翼的梢涡流场开展实验研究,重点关注椭圆水翼梢涡的发生发展规律、三维流动结构及涡核内部脉动特性等,尝试从流动特征的角度说明其与空化的联系。

    实验在浙江大学流体工程研究所的小型高速空泡水洞中进行。实验段尺寸为1000 mm×200 mm×200 mm,来流速度0.3~10 m/s连续可调,实验段湍流度小于1%。实验段侧壁配备有角编码器和步进电机,可用于调节模型迎角,角度精度为0.1°。所用实验模型为椭圆水翼,如图1 所示,其截面翼型为NACA662-415,根弦长c = 83.66 mm,半展长b = 100 mm,安装在实验段侧壁时,模型梢部正好处于实验段的中线上。图中xyz方向分别代表流向(轴向)、法向及展向,坐标原点位于水翼梢部。

    图  1  椭圆水翼
    Fig.  1  Elliptical hydrofoil

    采用高频层析PIV技术对水翼单相梢涡流场开展实验研究,如图2 所示。4台分辨率为1024 像素×1024 像素的Photron FASTCAM高速相机呈“×”形布置在水洞一侧(2台型号为SA4 200K M3,满画幅最大帧率3600帧/s;另外2台型号为Mini AX100,满画幅最大帧率4000 帧/s)。每台相机配备有尼康200 mm定焦微距镜头,为满足景深需求,镜头光圈f#=24。激光器系统为Nd-YLF高频脉冲激光器(镭宝光电Vlite-Hi-527-50),激光波长527 nm,单脉冲最大能量50 mJ,重频范围0.2~1.0 kHz。通过光学镜组将激光扩展为厚度18 mm的体光源,用以照亮测量区域。相机与激光器系统通过同步器和上位机控制软件同步控制,层析PIV测量系统采样频率为1 kHz,两帧粒子图像的曝光时间间隔Δt分别为100 μs和75 μs,以适用于不同的来流速度。实验采用平均粒径为50 μm的空心玻璃珠作为示踪粒子,密度约为1.05 g/cm3,跟随性良好。

    图  2  实验布置示意图
    Fig.  2  Schematic of the experimental setup

    采用北京立方天地科技有限公司MicroVec4 V1.2商用软件对粒子图像进行三维粒子重构及互相关运算。基于灰度增强倍增代数重构算法(IntE-MART)[25]重构三维粒子灰度场,三维重构体大小为51 mm(x)× 39 mm(y)× 18 mm(z),放大率为0.06 mm/像素,总共包含850×650×300个重构体素。采用多层迭代体变形算法[26]对重构粒子场进行互相关运算,总共迭代2次,第一次查询窗口大小为48 voxel×48 voxel×48 voxel,第二次查询窗口大小为32 voxel×32 voxel×32 voxel,每次互相关运算的查询窗口重叠率为75%,速度矢量间隔为0.48 mm(8 voxel),空间分辨率为1.92 mm(32 voxel),重构体内包含106×81×37个速度矢量。

    本文采用λci判据[27]辨识流场中的三维旋涡结构,该判据是一种基于速度梯度的旋涡识别方法。由于水翼梢涡流场中仅存在单个强度较大的流向旋涡,剪切层对涡识别的干扰较小,因而在本文中采用λci-xλcix方向上的分量)辨识流场中的梢涡结构,认定λci-x的局部极值点为涡核中心。从水翼下游向上游观察,梢涡在y-z平面内的旋转方向为顺时针,因此λci以及涡量在x方向上的分量均为负值。

    在梢涡这种典型尾涡系统中,存在旋涡摆动(vortex wandering)现象,该现象被认为是旋涡的涡轴空间位置随时间存在某种随机波动[28]。未考虑旋涡摆动的梢涡时均流场测量会引入额外的不确定性因素,从而影响测量结果的准确性,尤其是在单点测量(如探针、热线等)中。一般而言,未考虑旋涡摆动往往会得到一个“平滑”的时均涡流场,所获取的旋涡尺寸比真实尺寸偏大,而其最大切向速度偏小[29-30]。在本文中,采用旋涡中心对齐平均的方法滤除旋涡摆动对时均流场的影响[31-32]:将相同空间位置的时间序列数据按所确定的旋涡中心(λci-x局部极值点)对齐,再进行时间平均,从而得到摆动滤除后的时均数据。该方法简单明了,适用于不同类型的涡流场研究[2, 33]。旋涡摆动滤除前后的时均涡流场如图3 所示,图中ωx为轴向涡量,U为来流速度。可以看到,与图3 (a)直接时间平均的结果相比,图3 (b)所涵盖的流场区域变小,包含的流场信息也明显减少,即摆动滤除方法会导致一定量的数据损失,但摆动滤除后的时均流场仍包含了完整的梢涡结构,不会影响后续对梢涡流场的分析。

    图  3  摆动滤除前后的时均涡流场
    Fig.  3  Time-averaged flow field before and after vortex wandering correction

    本文主要研究不同来流速度(U=3.01、5.70 m/s)、不同迎角(α=5°、10°)下的水翼梢涡流场。基于水翼根弦长的雷诺数Rec分别为2.48×105、4.70×105,此时水洞实验段内空化数σ分别为22.60、5.80,流动处于单相状态,无空化现象,流向测量范围x/c =0.40~0.89。

    为说明旋涡摆动滤除的必要性,本节分析了梢涡的摆动特性。图4 给出了不同工况下梢涡向下游发展的运动轨迹,其中误差棒代表梢涡在当地位置的摆动幅度。可以看到,梢涡轨迹及摆动幅度在yz方向上呈各向异性:测量区域内,梢涡在展向上沿z轴负向偏移明显,而在法向上沿y轴正向略有偏移。对比不同实验状态,发现改变雷诺数不会对梢涡的运动轨迹产生明显影响,但不同迎角下梢涡在法向上的相对位置存在差异,α=10°状态下梢涡在法向上的空间位置更低,这与不同迎角下梢涡脱离水翼的空间位置相关。对于摆动幅度而言,梢涡在法向上的摆动幅度大于其在展向上的摆动幅度。

    图  4  近尾迹区域内涡核中心运动轨迹
    Fig.  4  Trajectory of the tip vortex core center in the near field

    图5 给出了不同工况下梢涡在近尾迹区域内的摆动幅度,图中AyAz分别为yz方向的摆动幅度。可以看到,梢涡在法向上的摆动幅度明显大于其在展向上的摆动,说明梢涡在法向上更不稳定,且向下游发展引起梢涡摆动的不稳定性会逐渐增强,摆动幅度沿流向放大,这与前人关于翼尖涡摆动的研究结论基本一致[34]。随着迎角增大,梢涡强度增大,摆动幅度有所减小,二者呈负相关,但雷诺数变化对梢涡摆动幅度的影响较小。

    图  5  近尾迹区域内梢涡摆动幅度
    Fig.  5  Amplitude of the vortex wandering in the near field

    为定量对比摆动滤除前后的时均梢涡速度场,定义相对偏差Δ如下:

    $$ \varDelta = \frac{{{E_{{\rm{uncor}}}} - {E_{{\rm{cor}}}}}}{{{E_{{\rm{cor}}}}}} \times 100\% $$ (1)

    式中,EuncorEcor分别代表摆动滤除前后的时均速度场,包括轴向速度Ux(速度方向沿流向,即x轴正方向)和切向速度Ut(极坐标下周向速度,即为yz方向速度的合速度,以顺时针方向为正)。

    图6 给出了摆动滤除前后的时均速度场。此处选取迎角α=10°状态进行分析,原因在于一定迎角范围内,梢涡强度随着水翼迎角增大而增强,而旋涡的摆动幅度与其强度呈弱负相关,讨论旋涡小幅度摆动状态对时均流场结果的影响更具说服力。图6 (e)、(f)中的洋红色曲线代表涡核边界(以λci-x局部极值的12.5%所绘制的等高线,2.2节将具体阐述涡核边界和涡核尺寸)。结果表明:旋涡的摆动运动对轴向速度分布影响不明显,摆动滤除前后相对偏差不大于1%(图6(f));相比而言,滤除旋涡摆动的影响后,涡核内部时均切向速度分布发生了较大变化,相对偏差在20%左右(图6(e)),涡核边界上相对偏差约为4%,其原因可能在于旋涡内部切向速度的空间梯度较大,直接平均会引入较大偏差。此外,摆动滤除后时均切向速度场分布更加规整,涡核形态接近圆形,而摆动滤除前的涡核呈“拉长”的椭圆形。

    图  6  旋涡摆动滤除前后时均速度场,α = 10°
    Fig.  6  Time-averaged velocity field before and after vortex wandering correction, α = 10°

    综上可知,梢涡在近尾迹区域内所表现出的摆动特性十分明显,未考虑旋涡摆动的直接时间平均会引入较大偏差,尤其是对时均切向速度场而言;而采用旋涡中心对齐平均方法能够有效滤除旋涡摆动所带来的影响,获取更为准确的梢涡流场。因此,后文对梢涡流场的研究均采用此方法获取时均涡流场。

    本节主要对椭圆水翼梢涡的三维流场结构展开分析。图7 给出了不同实验状态下梢涡的瞬时三维流场,其中,蓝色等值面代表测量流场中λci-x最小值的35%,用以展示梢涡的三维结构,黑、红色圆点分别代表所在平面与x-y平面、x-z平面平行的流线起始点,相应颜色虚线代表流线,用以说明梢涡附近的流动状态。结果表明:近尾迹区域内的梢涡尚未完全卷起,处于发生发展阶段,呈不规则的柱状形态,且梢涡在向下游发展过程中会向翼根方向(即z轴负方向)偏移。从流线看,周围的流体不断向梢涡卷绕,并向下游传播。随着迎角和来流速度增大,所产生的梢涡强度提升,其结构也更加稳定,影响区域也有所增大。

    图  7  三维瞬时梢涡流场结构
    Fig.  7  Three-dimensional instantaneous tip vortex structure

    图8 给出了不同工况下梢涡轴向速度三维分布,可观察到,梢涡所表现出的轴向速度特性存在明显差异。小迎角、低雷诺数状态下所产生的梢涡强度较弱,涡流场内的轴向速度普遍小于来流速度,呈尾迹型(wake-like profile)轴向流动(图8(a));当进一步增大迎角或来流雷诺数时,梢涡轴向流动由尾迹型转变为射流型(jet-like profile),即梢涡内部轴向流动速度大于来流速度(图8(b)、(c))。这种转变的原因在于旋涡强度增大,将梢部的分离剪切层卷绕在旋涡周围,从而使得旋涡免受尾迹流动的破坏性影响,使其内部轴向流动保持较高速度[35]图8 中青色和红色等值面分别代表2种不同轴向流动类型涡流场中的高速区域,结果表明:无论梢涡的轴向流动表现为尾迹型还是射流型,流场中的高速轴向流动区域并不完全出现于涡核内部,其空间位置与涡核中心存在一定的偏差,主要分布于梢涡的一侧,后文将对该现象进行具体说明。

    图  8  时均轴向速度三维分布
    Fig.  8  Three-dimensional distribution of time-averaged axial velocity

    α = 10°、Rec = 2.48×105为例,图9给出了梢涡在不同流向位置的轴向涡量ωx及轴向速度分布云图。图中洋红色圆点代表涡核中心。在x/c=0.4处,可明显观察到梢涡与水翼脱落剪切层的相互作用,剪切层逐渐卷绕在梢涡周围,裹挟梢涡向下游发展,此时涡量主要集中分布于涡核内部,且在向下游发展的过程中存在一定衰减,同时涡核尺寸膨胀,梢涡作用区域逐渐增大;高速轴向流动区域主要分布于涡核与水翼脱落剪切层之间,而不是涡核中心。出现该现象的原因在于:近尾迹区域内,水翼脱落剪切层对梢涡作用明显,涡核周围旋转流体受涡核和剪切层挤压存在轴向加速[36];在向下游发展过程中,涡核与剪切层之间的相互作用减弱,该区域内的轴向速度也逐渐减小。

    图  9  时均轴向涡量及轴向速度云图,α = 10°
    Fig.  9  Contour of time-averaged axial vorticity and velocity, α = 10°

    为说明水翼脱落剪切层对梢涡的影响情况,图10 给出了过涡核中心沿法向及展向的切向速度分布(α = 10°)。法向上,梢涡受剪切层作用明显,受作用侧的切向速度峰值明显较低,且速度分布带有低速拐点等不规则特征,另一侧则为单峰值曲线,旋涡中心切向速度处于较低水平。梢涡向下游发展,低速拐点逐渐向涡核外侧移动,说明剪切层对梢涡的影响逐渐减弱,同时旋涡中心两侧的速度峰值差减小,减小量约为来流速度的9.5%,涡核中心切向速度持续降低,以上特征均说明梢涡逐渐发展至完全卷起状态。在展向上,剪切层的作用不明显,涡核中心两侧速度分布基本对称。此外,涡核尺寸同样也是描述梢涡特性的重要参数之一,常采用涡核中心两侧切向速度峰值间距作为涡核的等效直径d。但近尾迹区域内梢涡处于高度不对称状态,仅通过等效直径无法准确给出涡核区域,因而本文采用以λci-x局部极值的12.5%所绘制的等高线确定涡核边界,对比发现法向和展向上的切向速度峰值均处于此边界上。而对于涡核尺寸,仍通过等效直径d定量描述(本文采用法向和展向上切向速度峰值间距的均值)。x/c=0.40时,其等效直径d约为4 mm,发展至x/c=0.89,d约为4.7 mm,同时法向、展向上的切向速度均随之降低,结合图9 中涡核内部涡量变化趋势可知,上述现象为旋涡在黏性涡核旋转作用下的内部涡量再分配过程。

    图  10  梢涡时均切向速度分布,α = 10°
    Fig.  10  Distribution of time-averaged circumferential velocity, α = 10°

    影响梢涡空化初生的因素主要有3类,分别为旋涡中心的平均最小压力、当地湍流压力脉动以及水质(含气量、气核分布)[5, 12-13]。其中,当地湍流压力脉动往往与局部的不稳定流动高度相关,且梢涡空化普遍初生于旋涡中心,因而本节主要探讨涡核内部的脉动特性,以说明涡核内部流动特征与梢涡空化的联系。

    图11 给出了α = 10°时不同流向位置的湍动能(Turbulence Kinetic Energy,TKE)云图,湍动能计算公式如下:

    图  11  湍动能云图,α = 10°
    Fig.  11  Contour of turbulence kinetic energy, α = 10°
    $$ k = \frac{1}{{2U_\infty ^2}}({u'^2} + {v'^2} + {w'^2}) $$ (2)

    式中,u'v'w'分别代表xyz方向上的均方根速度脉动。图中洋红色圆点代表涡核中心,青色曲线代表涡核边界(λci-x局部极值的12.5%)。梢涡流场的湍动能主要集中分布于涡核内部,梢涡向下游发展,涡核内部的湍流脉动增强,其涡核内部的湍动能分布不再呈单一峰值,且湍动能峰值不再出现于涡核中心。

    图12给出了湍动能沿法向的分布曲线(α = 10°),以说明涡核内部的湍动能分布情况。图12 (a)展示了未考虑旋涡摆动的湍动能分布情况,湍动能呈单一峰值的正态分布,不同流向位置的湍动能峰值均出现于涡核中心,这与Peng等[16]通过LDV获取的结果一致。滤除旋涡摆动影响后的结果如图12 (b)所示,在靠近水翼尾缘的流向位置,湍动能仍呈正态分布,但较图12 (a)中的分布曲线更加“扁平化”,涡核内部湍流能量分布更为分散,说明滤除了旋涡摆动对湍流能量分布的影响。当发展至下游位置(x/c = 0.89)时,涡核内部湍动能继续增强,但曲线在涡核中心位置出现“凹陷”,此时涡核中心区域的湍动能较其周围区域偏低,该位置的湍流能量衰减增强,梢涡内核存在再层流化(relaminarization)趋势,即在黏性涡核旋转作用下,旋涡内核中的湍流脉动受到抑制进而衰减,呈向未受扰自由来流转化的迹象[36]

    图  12  旋涡摆动滤除前后的法向湍动能分布,α = 10°
    Fig.  12  Distribution of turbulence kinetic energy along y axis before and after vortex wandering correction, α = 10°

    图13 给出了不同工况下涡核中心各速度分量脉动强度沿流向的变化情况,图中不同颜色虚线代表相应速度脉动沿流向平均值。在所有实验工况下,涡核中心速度脉动呈各向异性:轴向速度保持低水平稳定脉动,而展向速度脉动略强于法向速度脉动,约为轴向速度脉动的2倍,在测量区域内各速度分量脉动强度沿流向略有增强。此外,迎角改变会对脉动强度产生影响,迎角越大,各速度分量的脉动越强,其原因在于较大迎角下,从水翼梢部卷起的剪切层以及次生旋涡等流动结构与主涡相互作用更为剧烈,导致主涡涡核内部的湍流脉动增强;但迎角改变对各速度分量脉动强度分布规律无明显影响。同时,增大来流雷诺数也会使得各速度分量脉动强度有所增大。结合2.2节中的结论发现,无论水翼梢涡的轴向流动形式是高速射流型还是低速尾迹型,其轴向速度脉动较其他2个方向的速度脉动均处于较低水平,说明轴向流动形式同样对速度脉动规律无明显影响。

    图  13  不同流向位置涡核中心各速度分量脉动强度
    Fig.  13  Velocity fluctuation in the core center at different streamwise position

    在不考虑水质的前提下,水翼梢涡空化初生主要与涡核中心平均最小压力以及当地湍流压力脉动相关。而空化初生的第一要素始终是涡核中心平均最小压力,只有当该压力明显低于空化初生的压力阈值时才会出现稳定的空化现象,如图14 (a)所示;而当涡核中心平均最小压力与初生压力阈值接近时,湍流压力脉动则可能是不稳定空化现象的成因所在,如图14 (b)所示的间歇空化。因此,为探究梢涡空化与流动特征的联系,下面从涡核中心平均最小压力和当地湍流压力脉动两方面展开讨论。

    图  14  梢涡空化的不同类型
    Fig.  14  Different type of tip vortex cavitation

    涡核中心平均最小压力主要取决于旋涡强度,旋涡强度越大,其切向速度在径向上的梯度越大,涡核中心压力越低。涡核中心的平均切向速度接近于零,远小于当地的轴向速度,而当地速度一定程度上能够反映局部静压,即通过涡核中心的平均轴向速度沿流向的分布情况能够判断当地静压的大致分布。Asnaghi等[18]通过数值模拟方法也证明了这一观点,即涡核中心空化数与当地轴向速度高度相关,最大轴向速度与最小空化数(最低压力)所在位置相对应,但Asnaghi没有进一步考虑脉动因素的影响。对于湍流压力脉动,Stinebring等[11]总结了梢涡内部无量纲脉动压力公式:

    $$ \frac{{\sqrt {{{p}}{{{'}}^2}} }}{{\rho U_\infty ^2}} = K\frac{{u{'^2} + v{'^2} + w{'^2}}}{{2U_\infty ^2}} = K \cdot {{k}} $$ (3)

    式中,p'为涡核中心的脉动压力,K为正常数,与流场湍流特性相关,k为涡核中心湍动能。可以发现,涡核内部的无量纲压力脉动与湍动能呈正相关。结合图13 中的结果,表1进一步给出了不同工况下,各速度分量脉动对湍动能的贡献情况。从表中可知,法向和展向速度脉动贡献了90%以上的涡核中心湍动能,说明水翼梢涡中心的湍流脉动由这2个方向上的速度脉动主导,而流向速度脉动的影响十分有限,进一步结合公式(3)可知,在本文所有工况下,根据涡核中心法向和展向的速度脉动即可基本确定当地湍流压力脉动大小。

    表  1  涡核中心各速度分量脉动对湍动能的贡献
    Table  1  The contribution of each velocity component fluctuation to TKE in the core center
    α/(° )Rec0.5u' 2U–2k–10.5v' 2U–2k–10.5w' 2U–2k–1
    52.48×1057.6%40.4%52.0%
    4.70×1055.7%41.6%52.7%
    102.48×10510.3%40.6%49.1%
    下载: 导出CSV 
    | 显示表格

    基于上述结论,通过涡核中心平均轴向速度以及当地法向、展向速度脉动等流场信息,能够有效估计局部压力变化情况,这对预测水翼梢涡空化初生、完善空化初生预测模型具有重要的应用价值。

    本文通过水洞实验及高频层析PIV技术对椭圆水翼的梢涡流动开展了实验测量,清晰、完整地描述了梢涡的三维流动特性,不仅有助于典型尾涡流场的流动特性研究,同时也为梢涡空化数值模拟研究提供了三维数据参考。研究首先讨论了梢涡的摆动特性及其影响情况,突出了摆动滤除对涡流场定量分析的重要性;其次对梢涡三维流场结构展开分析,重点强调梢涡与剪切层间的相互作用及其对梢涡流场的影响情况;最后讨论了涡核内部的脉动特性,结合前人工作明晰了涡核中心湍流脉动压力的主要来源,阐明了流动特征与空化的联系。结果表明:

    1)近尾迹区域内,水翼梢涡存在各向异性的偏移及摆动(偏移:法向 < 展向;摆动:法向 > 展向)。未考虑旋涡摆动的直接时间平均会引入较大偏差,影响梢涡时均流场、脉动流场结果的准确性。采用旋涡中心对齐平均方法能有效滤除摆动带来的影响。

    2)在近尾迹区域内,水翼梢涡具有典型的三维特性,其瞬时三维形态不规则,且与水翼脱落剪切层相互作用明显,具体表现为梢涡流场中的高速轴向流动区域出现于涡核与剪切层之间,受剪切层作用侧的切向速度明显较低。向下游发展,剪切层的作用明显减弱,高速轴向流动也随之衰减,涡核中心两侧切向速度峰值差降低,梢涡尺寸增大。

    3)梢涡流场中湍流脉动集中分布于涡核内部,且该湍流脉动主要来自法向、展向上的速度脉动。基于无量纲脉动压力公式进一步发现,这2个方向上的速度脉动是涡核内部湍流脉动压力的主要来源。

  • 图  1   椭圆水翼

    Fig.  1   Elliptical hydrofoil

    图  2   实验布置示意图

    Fig.  2   Schematic of the experimental setup

    图  3   摆动滤除前后的时均涡流场

    Fig.  3   Time-averaged flow field before and after vortex wandering correction

    图  4   近尾迹区域内涡核中心运动轨迹

    Fig.  4   Trajectory of the tip vortex core center in the near field

    图  5   近尾迹区域内梢涡摆动幅度

    Fig.  5   Amplitude of the vortex wandering in the near field

    图  6   旋涡摆动滤除前后时均速度场,α = 10°

    Fig.  6   Time-averaged velocity field before and after vortex wandering correction, α = 10°

    图  7   三维瞬时梢涡流场结构

    Fig.  7   Three-dimensional instantaneous tip vortex structure

    图  8   时均轴向速度三维分布

    Fig.  8   Three-dimensional distribution of time-averaged axial velocity

    图  9   时均轴向涡量及轴向速度云图,α = 10°

    Fig.  9   Contour of time-averaged axial vorticity and velocity, α = 10°

    图  10   梢涡时均切向速度分布,α = 10°

    Fig.  10   Distribution of time-averaged circumferential velocity, α = 10°

    图  11   湍动能云图,α = 10°

    Fig.  11   Contour of turbulence kinetic energy, α = 10°

    图  12   旋涡摆动滤除前后的法向湍动能分布,α = 10°

    Fig.  12   Distribution of turbulence kinetic energy along y axis before and after vortex wandering correction, α = 10°

    图  13   不同流向位置涡核中心各速度分量脉动强度

    Fig.  13   Velocity fluctuation in the core center at different streamwise position

    图  14   梢涡空化的不同类型

    Fig.  14   Different type of tip vortex cavitation

    表  1   涡核中心各速度分量脉动对湍动能的贡献

    Table  1   The contribution of each velocity component fluctuation to TKE in the core center

    α/(° )Rec0.5u' 2U–2k–10.5v' 2U–2k–10.5w' 2U–2k–1
    52.48×1057.6%40.4%52.0%
    4.70×1055.7%41.6%52.7%
    102.48×10510.3%40.6%49.1%
    下载: 导出CSV
  • [1]

    KÜCHEMANN D. Report on the I. U. T. A. M. symposium on concentrated vortex motions in fluids[J]. Journal of Fluid Mechanics,1965,21(1):1-20. doi: 10.1017/s0022112065000010

    [2]

    DREYER M,DECAIX J,MÜNCH-ALLIGNÉ C,et al. Mind the gap: a new insight into the tip leakage vortex using stereo-PIV[J]. Experiments in Fluids,2014,55(11):1-13. doi: 10.1007/s00348-014-1849-7

    [3]

    FELLI M,FALCHI M. Propeller tip and hub vortex dynamics in the interaction with a rudder[J]. Experiments in Fluids,2011,51(5):1385-1402. doi: 10.1007/s00348-011-1162-7

    [4] 潘森森, 彭晓星. 空化机理[M]. 北京: 国防工业出版社, 2013.

    PAN S S, PENG X X. Physical mechanism of cavitation[M]. Beijing: National Defense Industry Press, 2013.

    [5] 刘玉文,徐良浩,宋明太,等. 水翼叶梢涡空化实验研究进展[J]. 实验流体力学,2020,34(5):1-11. DOI: 10.11729/syltlx20190083

    LIU Y W,XU L H,SONG M T,et al. Experimental research progress of hydrofoil tip vortex cavitation[J]. Journal of Experiments in Fluid Mechanics,2020,34(5):1-11. doi: 10.11729/syltlx20190083

    [6]

    McCORMICK B W Jr. On cavitation produced by a vortex trailing from a lifting surface[J]. Journal of Basic Enginee-ring,1962,84(3):369-379. doi: 10.1115/1.3657328

    [7]

    BILLET M L,HOLL J W. Scale effects on various types of limited cavitation[J]. Journal of Fluids Engineering,1981,103(3):405-414. doi: 10.1115/1.3240800

    [8]

    PAUCHET A,BRIANGON-MARJOLLET L. Recent results on tip vortex cavitation scale effects at high Reynolds numbers[J]. WIT Transactions on The Built Environment,1993(1):1-8. doi: 10.2495/NEVA930131

    [9]

    FRUMAN D H,CERRUTTI P,PICHON T,et al. Effect of hydrofoil planform on tip vortex roll-up and cavitation[J]. Journal of Fluids Engineering,1995,117(1):162-169. doi: 10.1115/1.2816806

    [10]

    POGOZELSKI E, SHEKARRIZ A, KATZ J, et al. Three dimensional near field behavior of a tip vortex developing on an elliptic foil[C]//Proc of the 31st Aerospace Sciences Meeting. 1993: 865. doi: 10.2514/6.1993-865

    [11]

    STINEBRING D R,FARRELL K J,BILLET M L. The structure of a three-dimensional tip vortex at high Reynolds numbers[J]. Journal of Fluids Engineering,1991,113(3):496-503. doi: 10.1115/1.2909524

    [12]

    ARNDT R E A. Cavitation in vortical flows[J]. Annual Review of Fluid Mechanics,2002,34(1):143-175. doi: 10.1146/annurev.fluid.34.082301.114957

    [13]

    ZHANG L X,ZHANG N,PENG X X,et al. A review of studies of mechanism and prediction of tip vortex cavitation inception[J]. Journal of Hydrodynamics,2015,27(4):488-495. doi: 10.1016/s1001-6058(15)60508-x

    [14]

    PENNINGS P C,WESTERWEEL J,TERWISGA T J C. Flow field measurement around vortex cavitation[J]. Expe-riments in Fluids,2015,56(11):1-13. doi: 10.1007/s00348-015-2073-9

    [15]

    DREYER M. Mind the gap: tip leakage vortex dynamics and cavitation in axial turbines[D]. Lausanne: École Poly-technique Fédérale de Lausanne, 2015. doi: 10.5075/epfl-thesis-6611

    [16]

    PENG X X,XU L H,LIU Y W,et al. Experimental measurement of tip vortex flow field with/without cavitation in an elliptic hydrofoil[J]. Journal of Hydrodynamics,2017,29(6):939-953. doi: 10.1016/s1001-6058(16)60808-9

    [17]

    BOSSCHERS J. An analytical and semi-empirical model for the viscous flow around a vortex cavity[J]. International Journal of Multiphase Flow,2018,105:122-133. doi: 10.1016/j.ijmultiphaseflow.2018.03.021

    [18]

    ASNAGHI A,SVENNBERG U,BENSOW R E. Large eddy simulations of cavitating tip vortex flows[J]. Ocean Enginee-ring,2020,195:106703. doi: 10.1016/j.oceaneng.2019.106703

    [19]

    CHANG N,GANESH H,YAKUSHIJI R,et al. Tip vortex cavitation suppression by active mass injection[J]. Journal of Fluids Engineering,2011,133(11):111301. doi: 10.1115/1.4005138

    [20]

    LEE S J,SHIN J W,ARNDT R E A,et al. Attenuation of the tip vortex flow using a flexible thread[J]. Experiments in Fluids,2017,59(1):1-12. doi: 10.1007/s00348-017-2476-x

    [21]

    AMINI A,SEO J,RHEE S H,et al. Mitigating tip vortex cavitation by a flexible trailing thread[J]. Physics of Fluids,2019,31(12):127103. doi: 10.1063/1.5126376

    [22]

    AMINI A,RECLARI M,SANO T,et al. Suppressing tip vortex cavitation by winglets[J]. Experiments in Fluids,2019,60(11):1-15. doi: 10.1007/s00348-019-2809-z

    [23]

    ASNAGHI A,SVENNBERG U,GUSTAFSSON R,et al. Investigations of tip vortex mitigation by using roughness[J]. Physics of Fluids,2020,32(6):065111. doi: 10.1063/5.0009622

    [24]

    SVENNBERG U,ASNAGHI A,GUSTAFSSON R,et al. Experimental analysis of tip vortex cavitation mitigation by controlled surface roughness[J]. Journal of Hydrodynamics,2020,32(6):1059-1070. doi: 10.1007/s42241-020-0073-6

    [25]

    WANG H P,GAO Q,WEI R J,et al. Intensity-enhanced MART for tomographic PIV[J]. Experiments in Fluids,2016,57(5):1-19. doi: 10.1007/s00348-016-2176-y

    [26]

    SCARANO F. Iterative image deformation methods in PIV[J]. Measurement Science and Technology,2002,13(1):1-19. doi: 10.1088/0957-0233/13/1/201

    [27]

    ZHOU J,ADRIAN R J,BALACHANDAR S,et al. Mecha-nisms for generating coherent packets of hairpin vortices in channel flow[J]. Journal of Fluid Mechanics,1999,387:353-396. doi: 10.1017/s002211209900467x

    [28]

    BHAGWAT M J,RAMASAMY M. Effect of tip vortex aperiodicity on measurement uncertainty[J]. Experiments in Fluids,2012,53(5):1191-1202. doi: 10.1007/s00348-012-1348-7

    [29]

    DEVENPORT W J,RIFE M C,LIAPIS S I,et al. The structure and development of a wing-tip vortex[J]. Journal of Fluid Mechanics,1996,312:67-106. doi: 10.1017/s0022112096001929

    [30]

    IUNGO G V,SKINNER P,BURESTI G. Correction of wandering smoothing effects on static measurements of a wing-tip vortex[J]. Experiments in Fluids,2009,46(3):435-452. doi: 10.1007/s00348-008-0569-2

    [31] 薛栋,潘翀,袁先士,等. 低雷诺数下翼尖涡统计特性实验研究[J]. 实验流体力学,2019,33(5):36-41. DOI: 10.11729/syltlx20180129

    XUE D,PAN C,YUAN X S,et al. Experimental investiga-tion on the characteristics of wingtip vortex at low Reynolds number[J]. Journal of Experiments in Fluid Mechanics,2019,33(5):36-41. doi: 10.11729/syltlx20180129

    [32]

    HEYES A L, JONES R F, SMITH D A R. Wandering of wing-tip vortices[C]//Proceedings of 12th international sym-posium on the applications of laser techniques to fluid me-chanics. 2004.

    [33]

    BERESH S J,HENFLING J F,SPILLERS R W. Meander of a fin trailing vortex and the origin of its turbulence[J]. Experiments in Fluids,2010,49(3):599-611. doi: 10.1007/s00348-010-0825-0

    [34] 邱思逸,程泽鹏,向阳,等. 基于线性稳定性分析的翼尖涡摇摆机制[J]. 航空学报,2019,40(8):122712.

    QIU S Y,CHENG Z P,XIANG Y,et al. Mechanism of wingtip vortex wandering based on linear stability analysis[J]. Acta Aeronautica et Astronautica Sinica,2019,40(8):122712.

    [35]

    LEE T,PEREIRA J. Nature of wakelike and jetlike axial tip vortex flows[J]. Journal of Aircraft,2010,47(6):1946-1954. doi: 10.2514/1.c000225

    [36]

    SKINNER S N,GREEN R B,ZARE-BEHTASH H. Wingtip vortex structure in the near-field of swept-tapered wings[J]. Physics of Fluids,2020,32(9):095102. doi: 10.1063/5.0016353

  • 期刊类型引用(4)

    1. 陈翰,何聪,许祖锋,张远媛,张为,朱城磊. 柔性压力传感阵列弯曲应变与温变联合补偿方法. 传感技术学报. 2024(03): 423-429 . 百度学术
    2. 季承,刘俊,程进,赵天祥. 基于位置补偿的蒙皮数字波束形成方法. 信息记录材料. 2024(06): 211-213 . 百度学术
    3. 全栋梁,赵雄涛,张宇鹏. 空天飞行器结构与热防护系统健康监测技术的需求和挑战. 空天技术. 2023(01): 84-97 . 百度学术
    4. 黄俊,郭雨欣,冀晶晶,黄永安. 基于数据同化的气动压力稀疏重构方法. 实验流体力学. 2023(05): 9-17 . 本站查看

    其他类型引用(2)

图(14)  /  表(1)
计量
  • 文章访问数:  2463
  • HTML全文浏览量:  314
  • PDF下载量:  90
  • 被引次数: 6
出版历程
  • 收稿日期:  2021-08-25
  • 修回日期:  2021-11-29
  • 录用日期:  2021-12-01
  • 网络出版日期:  2022-05-25
  • 刊出日期:  2022-05-18

目录

/

返回文章
返回
x 关闭 永久关闭