飞行器表面气动载荷的柔性智能蒙皮多参量测量

郭栋梁, 侯超, 朱臣, 熊文楠, 陈爽, 许晓斌, 杨华, 黄永安

郭栋梁,侯超,朱臣,等. 飞行器表面气动载荷的柔性智能蒙皮多参量测量[J]. 实验流体力学,2022,36(2):146-154. DOI: 10.11729/syltlx20210115
引用本文: 郭栋梁,侯超,朱臣,等. 飞行器表面气动载荷的柔性智能蒙皮多参量测量[J]. 实验流体力学,2022,36(2):146-154. DOI: 10.11729/syltlx20210115
GUO D L,HOU C,ZHU C,et al. Multi-parameter measurement of aerodynamic load via flexible sensing skin[J]. Journal of Experiments in Fluid Mechanics, 2022,36(2):146-154.. DOI: 10.11729/syltlx20210115
Citation: GUO D L,HOU C,ZHU C,et al. Multi-parameter measurement of aerodynamic load via flexible sensing skin[J]. Journal of Experiments in Fluid Mechanics, 2022,36(2):146-154.. DOI: 10.11729/syltlx20210115

飞行器表面气动载荷的柔性智能蒙皮多参量测量

基金项目: 国家重点研发计划资助(2020YFA0405700);国家自然科学基金(51635007,51925503);中央高校基本科研业务费专项资金(2019kfyRCPY017)
详细信息
    作者简介:

    郭栋梁: (1992—),男,河北邢台人,博士研究生。研究方向:飞行器智能蒙皮技术,柔性传感技术,风洞测试技术。通信地址:湖北省武汉市洪山区珞喻路1037号华中科技大学(430074)。E-mail:78312889@qq.com

    侯超: (1996—),男,四川南充人,博士研究生。研究方向:飞行器智能蒙皮技术,柔性传感技术,形状重构技术。通信地址:湖北省武汉市洪山区珞喻路1037号华中科技大学(430074)。E-mail:809762746@qq.com

    通讯作者:

    黄永安: E-mail:yahuang@hust.edu.cn

  • 中图分类号: TP212.9;V211.753

Multi-parameter measurement of aerodynamic load via flexible sensing skin

  • 摘要: 飞行器表面气动参数特征是飞行器结构设计和安全评估的重要依据,而风洞试验作为最有效的测试手段,通常面临破坏结构、测量物理量单一等问题。提出曲面共形的柔性智能蒙皮测量技术,集成了多种超薄柔性传感器阵列,通过剪纸–拼接的完全共形方式集成到飞行器结构表面,在不改变结构表面形貌的情况下同步实时测量壁面静态压力、脉动压力、温度、壁面剪应力等多种气动参数。在直流式风洞、射流平台和FL–9风洞中对NACA0012机翼和飞行器尾翼进行了变风速和变迎角试验,分析风洞试验中采集获得的多种气动参数,验证了该系统的可用性,为风洞试验中柔性智能蒙皮多参量同步测量气动特性研究提供参考。
    Abstract: Aerodynamic characteristics of aircraft are critical for the design and security evaluation of aircraft structures. Wind tunnel test is the most effective experimental method, which is faced with the problems of structure damage and fewer parameters measurement. Herein, we propose a conformal measurement technique based on flexible sensing skin integrated with variety of ultra-thin flexible sensors, that are attached on surface of aircrafts consistently by kirigami assembly strategy. The sensing skin can simultaneously obtain multiple aerodynamic parameters like static pressure, pulsating pressure, temperature and wall shear force, without changing the surface morphology of the structure. Experiments on NACA0012 wing and aircraft tail under variable wind speed and variable angle of attack are carried out in the direct wind tunnel, jet platform and FL–9 wind tunnel. The characteristics of parameters collected in the wind tunnel are analyzed, which proves the availability of the system, paving a way for simultaneous measurement of various aerodynamic characteristics by flexible sensing skin in wind tunnel experiments.
  • 旋涡作为一种常见的流动现象广泛存在于自然界以及近现代工程中,被喻为流体运动的“肌腱”[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   Flexible sensing skin and flow field perception

    图  2   柔性智能蒙皮系统

    Fig.  2   System of flexible sensing skin

    图  3   柔性传感器及测量原理示意图

    Fig.  3   Flexible sensors and schematic diagrams

    图  4   柔性智能蒙皮系统风洞试验

    Fig.  4   Wind tunnel experiment of sensing skin system

    图  5   NACA0012机翼小型风洞试验脉动压力结果

    Fig.  5   Results of pulsating pressure in wind tunnel test with different angles of attack

    图  6   静态差压标定和升力系数随NACA0012机翼迎角变化

    Fig.  6   Results of static pressure in wind tunnel test with different angles of attack

    图  7   柔性温度传感器试验结果图

    Fig.  7   Calibration of flexible temperature sensor and wind tunnel test results

    图  8   热膜传感器在尾翼不同位置的响应

    Fig.  8   Voltage signal vibration of hot film sensors in aircraft tail

    图  9   飞行器颤振信号联合分析

    Fig.  9   Analysis of multi-signals for aircraft flutter test

    表  1   传感器性能参数

    Table  1   Performance parameters of sensors

    传感器类型敏感区尺寸厚度/μm响应时间/ms采样频率/Hz分辨率量程
    动压ϕ1.5 mm25550010
    静压ϕ3 mm6010010005 Pa0~6 kPa
    热膜1 mm×0.08 mm2021000
    温度1.5 mm×3.5 mm1010020.5 ℃–20~120 ℃
    下载: 导出CSV
  • [1]

    TABRIZIAN A,TATAR M,MASDARI M,et al. An experimental study on boundary layer transition detection over a pitching supercritical airfoil using hot-film sensors[J]. International Journal of Heat and Fluid Flow,2020,86:108743. doi: 10.1016/j.ijheatfluidflow.2020.108743

    [2]

    RUDMIN D,BENAISSA A,POIREL D. Detection of laminar flow separation and transition on a NACA-0012 airfoil using surface hot-films[J]. Journal of Fluids Engineering,2013,135(10):101104. doi: 10.1115/1.4024807

    [3] 柴聪聪,易贤,郭磊,等. 基于BP神经网络的冰形特征参数预测[J]. 实验流体力学,2021,35(3):16-21. DOI: 10.11729/syltlx20200016

    CHAI C C,YI X,GUO L,et al. Prediction of ice shape characteristic parameters based on BP nerual network[J]. Journal of Experiments in Fluid Mechanics,2021,35(3):16-21. doi: 10.11729/syltlx20200016

    [4] 陈舒越,郭向东,王梓旭,等. 结冰风洞过冷大水滴粒径测量初步研究[J]. 实验流体力学,2021,35(3):22-29. DOI: 10.11729/syltlx20200104

    CHEN S Y,GUO X D,WANG Z X,et al. Preliminary research on size measurement of supercooled large droplet in icing wind tunnel[J]. Journal of Experiments in Fluid Mechanics,2021,35(3):22-29. doi: 10.11729/syltlx20200104

    [5] 高文涛,张武林,寇宝智. T型尾翼飞机抖振试飞研究[J]. 实验流体力学,2020,34(6):45-51. DOI: 10.11729/syltlx20200040

    GAO W T,ZHANG W L,KOU B Z. Study on buffet flight test of aircraft with T-tail[J]. Journal of Experiments in Fluid Mechanics,2020,34(6):45-51. doi: 10.11729/syltlx20200040

    [6] 余涛,张威,张毅锋,等. 一种非介入式高超声速边界层不稳定波的测量方法[J]. 实验流体力学,2019,33(5):70-75. DOI: 10.11729/syltlx20190076

    YU T,ZHANG W,ZHANG Y F,et al. Focused laser differential interferometry measurement of instability wave in a hypersonic boundary-layer[J]. Journal of Experiments in Fluid Mechanics,2019,33(5):70-75. doi: 10.11729/syltlx20190076

    [7]

    MO Z W,FU H Z,HO Y S. Global development and trend of wind tunnel research from 1991 to 2014: a bibliometric analysis[J]. Environmental Science and Pollution Research International,2018,25(30):30257-30270. doi: 10.1007/s11356-018-3019-6

    [8]

    HUANG Y A,ZHU C,XIONG W N,et al. Flexible smart sensing skin for “Fly-by-Feel” morphing aircraft[J]. Science China Technological Sciences,2022,65(1):1-29. doi: 10.1007/s11431-020-1793-0

    [9]

    KUESTER M S, BORGOLTZ A, DEVENPORT W J. Pressure tap effects on the lift measurement of an airfoil section[C]//32nd AIAA Aerodynamic Measurement Technology and Ground Testing Conference.2016: 3654. doi: 10.2514/6.2016-3654

    [10] 汪玉,邱雷,黄永安. 面向飞行器结构健康监测智能蒙皮的柔性传感器网络综述[J]. 航空制造技术,2020,63(15):60-69,80. DOI: 10.16080/j.issn1671-833x.2020.15.060

    WANG Y,QIU L,HUANG Y A. Review of flexible sensor networks for structural health monitoring of aircraft smart skin[J]. Aeronautical Manufacturing Technology,2020,63(15):60-69,80. doi: 10.16080/j.issn1671-833x.2020.15.060

    [11]

    LI X T,LIU X,LIU Y Z,et al. Experimental study of near-wall underexpanded jet impingement on a flat plate using temperature-insensitive semi-transparent pressure-sensitive paint[J]. Experiments in Fluids,2020,61(11):1-17. doi: 10.1007/s00348-020-03068-5

    [12] 霍俊杰, 易仕和, 牛海波, 等. 基于温敏漆技术的圆锥高超声速大攻角绕流背风面流动结构实验研究[J/OL]. 气体物理, [2022-03-17]. https://kns.cnki.net/kcms/detail/detail.aspx?dbcode=CAPJ&#38;dbname=CAPJLAST&#38;filename=QTWL20210420000&#38;uniplatform=NZKPT&#38;v=5V9HsVG-07ysOtbqFkzbUOBCzSI1rjf0gx8QsqySKs7E_OEB73Fp1Na-RDwblyLh. doi: 10.19527/j.cnki.2096-1642.0906.
    [13]

    ZHU C,GUO D L,YE D,et al. Flexible PZT-integrated, bilateral sensors via transfer-free laser lift-off for multimodal measurements[J]. ACS Applied Materials & Interfaces,2020,12(33):37354-37362. doi: 10.1021/acsami.0c10083

    [14]

    CHEN F R,HOU C,JIANG S,et al. Mechanically-compensated bending-strain measurement of multilayered paper-like electronics via surface-mounted sensor[J]. Composite Structures,2021,277:114652. doi: 10.1016/j.compstruct.2021.114652

    [15]

    DI LUCA M,MINTCHEV S,SU Y X,et al. A bioinspired Separated Flow wing provides turbulence resilience and aerodynamic efficiency for miniature drones[J]. Science Robotics,2020,5(38):eaay8533. doi: 10.1126/scirobotics.aay8533

    [16] 杨华,黄永安. 大型风洞全场智能感知的研究进展[J]. 中国科学基金,2017,31(5):450-460.

    YANG H,HUANG Y A. Advance in the full-field intelligent sensing for large-scale wind tunnels[J]. Bulletin of National Natural Science Foundation of China,2017,31(5):450-460.

    [17] 赵建国,马炳和,邓进军. 用于壁面剪应力测量的柔性热敏传感器阵列[J]. 微纳电子技术,2009,46(7):414-418. DOI: 10.3969/j.issn.1671-4776.2009.07.006

    ZHAO J G,MA B H,DENG J J. Flexible thermal sensor array for wall shear stress measurements[J]. Micronanoelectronic Technology,2009,46(7):414-418. doi: 10.3969/j.issn.1671-4776.2009.07.006

    [18]

    SUN B Y,WANG P B,LUO J,et al. A flexible hot-film sensor array for underwater shear stress and transition measurement[J]. Sensors (Basel, Switzerland),2018,18(10):3469. doi: 10.3390/s18103469

    [19]

    XIONG W N,GUO D L,YANG Z X,et al. Conformable, programmable and step-linear sensor array for large-range wind pressure measurement on curved surface[J]. Science China Technological Sciences,2020,63(10):2073-2081. doi: 10.1007/s11431-020-1642-4

    [20]

    BIAN J,CHEN F,YANG B,et al. Laser-induced interfacial spallation for controllable and versatile delamination of flexible electronics[J]. ACS Applied Materials & Interfaces,2020,12(48):54230-54240. doi: 10.1021/acsami.0c18951

    [21] 梁撑刚,张琳,刘华伟. 翼型动态实验中脉动压力信号处理方法研究[J]. 弹箭与制导学报,2021,41(2):111-113,119. DOI: 10.15892/j.cnki.djzdxb.2021.02.023

    LIANG C G,ZHANG L,LIU H W. Research on fluctuating pressure signal processing in airfoil dynamic experiment[J]. Journal of Projectiles, Rockets, Missiles and Guidance,2021,41(2):111-113,119. doi: 10.15892/j.cnki.djzdxb.2021.02.023

    [22]

    POELS A, RUDMIN D, BENAISSA A, et al. Localization of flow separation and transition over a pitching NACA 0012 airfoil at transitional reynolds numbers using hot-films[J]. Journal of Fluids engineering-transactions of the asme, 2015, 137(12). doi: 10.1115/1.4031008

  • 期刊类型引用(4)

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

    其他类型引用(2)

图(9)  /  表(1)
计量
  • 文章访问数:  2930
  • HTML全文浏览量:  413
  • PDF下载量:  174
  • 被引次数: 6
出版历程
  • 收稿日期:  2021-08-29
  • 修回日期:  2022-03-03
  • 录用日期:  2022-03-15
  • 网络出版日期:  2022-05-25
  • 刊出日期:  2022-05-18

目录

/

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