仿座头鲸鳍肢前缘单/双凸起翼段失速特性的实验研究

时德伟, 陈乐天, 范子椰, 姜楠, 左志钢, 唐湛棋

时德伟, 陈乐天, 范子椰, 等. 仿座头鲸鳍肢前缘单/双凸起翼段失速特性的实验研究[J]. 实验流体力学, doi: 10.11729/syltlx20230165.
引用本文: 时德伟, 陈乐天, 范子椰, 等. 仿座头鲸鳍肢前缘单/双凸起翼段失速特性的实验研究[J]. 实验流体力学, doi: 10.11729/syltlx20230165.
SHI D W, CHEN L T, FAN Z Y, et al. Experimental study on stall characteristics of wing sections with leading-edge protuberances inspired by humpback whale flipper[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20230165.
Citation: SHI D W, CHEN L T, FAN Z Y, et al. Experimental study on stall characteristics of wing sections with leading-edge protuberances inspired by humpback whale flipper[J]. Journal of Experiments in Fluid Mechanics, doi: 10.11729/syltlx20230165.

仿座头鲸鳍肢前缘单/双凸起翼段失速特性的实验研究

基金项目: 国家工程研究海上风电研究中心项目(HSFD22001);国家自然科学基金资助项目(12272265, 12332017)
详细信息
    作者简介:

    时德伟: (1996—),男,辽宁鞍山人,硕士研究生。研究方向:流动分离与控制。通信地址:天津市津南区海河教育园区雅观路135号天津大学36号力学大楼412室(300354)。E-mail:13512993936@163.com

    通讯作者:

    左志钢: E-mail: Zhigang200@tsinghua.edu.cn

    唐湛棋: zhanqitang@tju.edu.cn

  • 中图分类号: O353.5

Experimental study on stall characteristics of wing sections with leading-edge protuberances inspired by humpback whale flipper

  • 摘要:

    近年来,相关研究表明仿座头鲸鳍肢前缘凸起可以有效延缓翼段失速或者改善失速后翼段气动性能。为研究仿生前缘凸起翼型的流动分离现象,针对二维平直翼段 NACA 634-021以及带有单个和两个前缘凸起的改型翼段开展了粒子图像测速(PIV)实验研究。对比了不同迎角下各翼段吸力面的平均速度、湍动能等统计量的分布,并对实验数据进行了本征正交分解(POD)模态分析,从能量的角度对失速过程中的流动结构进行了研究。研究发现,单、双凸起翼段都具有两步失速特性,当发生单侧失速时,凸峰处存在稳定的附着流动,抑制了失速区向凸起另一侧的扩展,起到了类似于翼刀的作用;相比于单凸起翼段,双凸起翼段在更小的迎角发生失速模式的转换;POD分解的各阶模态的能量大小与翼段所处的失速模式有关,失速模式的切换同时导致流场中流动结构的特征发生变化;单凸起翼段处于单侧失速与双侧失速转换的临界迎角时,流动分离模式具有不稳定性。本研究为仿生前缘凸起翼型的理论研究和工程应用提供了实验参考。

    Abstract:

    In recent years, relevant studies have shown that wing sections with humpback whale-like leading-edge protuberances can effectively delay stall and improve aerodynamic performance after stall. To investigate the flow separation phenomenon of the biomimetic leading-edge protuberances wing section, a particle image velocimetry (PIV) experimental study on the NACA 634-021 wing section and modified wing sections with single and double leading-edge protuberances is carried out. The distributions of statistics such as the average velocity of the wing section suction surface and the turbulent kinetic energy at different angles of attack are compared. The proper orthogonal decomposition (POD) method is adopted to analyze the experimental data. The flow structures during the stall process are analyzed from an energy perspective. The study finds that: both single and double leading-edge protuberances wing sections have two-step stall characteristics. When a unilateral stall occurs, there is a stable attached flow at the convex peak, which suppresses the extension of the stall zone to the other side of the convex peak and plays a role similar to a wing blade. Compared to the single leading-edge protuberance wing section, the wing section with double leading-edge protuberances undergoes stall mode transition at smaller angles of attack. The energy of each mode of POD decomposition is related to the stall mode at which the wing section is located. The characteristics of the flow structure in the flow field change with the switching of stall modes. The flow separation mode of the single leading-edge protuberance wing section exhibits instability during the transition from unilateral stall to bilateral stall. This study provides an experimental reference for the theoretical research and potential engineering application of biomimetic leading-edge protuberances wing sections.

  • 翼面绕流问题广泛存在于飞机和流体机械的运行中,翼型在大迎角下会出现失速现象,吸力面发生流动分离,导致升力陡然下降,阻力大幅升高,引起气动性能和机组效率的降低。因此对流动分离现象的研究与控制具有重要的工程应用价值。

    座头鲸是海洋中最大的生物之一,其鳍肢前缘具有独特的波状凸起。在捕食时,它先快速游过猎物身边,形成气泡云聚集猎物[1],再完成一个快速的转弯,回过头来穿过气泡云捕食猎物[2]。在转弯过程中,鳍肢依靠升力产生向心力,当迎角很大时,有可能出现失速现象。因此,座头鲸捕猎时的高速转弯要求其鳍肢在大迎角下具有良好的性能。通过仿生学研究,Fish等[3]发现座头鲸鳍肢的剖面形状与翼型相似,前缘凸起轮廓类似于正弦曲线,提出这种凸起结构具有在大迎角下延缓失速的作用[4]。Miklosovic等[5]在风洞中分别对带有前缘凸起的和平滑的座头鲸鳍肢模型进行了测力实验,对比结果发现,具有前缘凸起的模型最大升力提高了6%,失速迎角推迟了40%。 Pedro[6]、Webber [7]、Johari[8]和Zhang[9]等人分别采用数值和实验方法研究了仿鲸鱼鳍前缘凸起翼段的气动性能,验证了前缘凸起对翼段失速特性的改善作用。Cai等[10]对平滑翼段、具有单个前缘凸起的翼段进行了实验研究,发现平滑翼段在达到失速迎角后发生升力骤降,而单凸起翼段具有特殊的两步失速特性。在平滑翼段失速迎角之前,单凸起翼段先发生第一次失速,翼段仍有相当一部分区域处于未失速状态,升力系数介于平滑翼段失速前后之间;在平滑翼段失速迎角之后,单凸起翼段发生第二次失速,升力系数降至与平滑翼段失速后相当,并在凸峰位置观测到一股附着流动。同时Cai等[11]对具有多个前缘凸起的翼段进行了实验和数值分析,发现其绕流流场在小迎角下呈周期性分布,在大迎角下则具有非周期性且流场出现分区效应,并发现多凸起翼段的失速过程较平缓,在大迎角时升力系数明显高于平滑翼段。李德友等[12]针对Cai等学者的研究成果做了更进一步的数值分析,发现双凸起翼段也呈现两步失速特性,同时提出前缘凸起引发的不对称流向涡是产生不同失速现象的主要原因,不同前缘凸起数,对流动的作用机理也不同。Zhao等[13]采用本征正交分解方法研究了平滑翼段和多个前缘凸起翼段的绕流流场,提出卡门涡在平滑翼段尾迹中占主导地位,而前缘凸起结构抑制了这种旋涡脱落过程。

    Yang等[14]使用PIV实验方法研究了前缘凸起对低展弦比机翼模型流动的影响,发现前缘凸起可以有效缓解模型的流动分离,并降低尾流涡对流动脉动能量的贡献。Ghazi等[15]使用RANS方法对具有前缘凸起的涡轮机叶片进行了数值模拟,发现对比平滑前缘叶片,改型叶片的功率和扭矩系数显著增加。Dong等[16]研究了前缘凸起结构的振幅和波长对压气机叶片流动控制效果的影响,并分析了前缘凸起结构轮廓的设计方法。除此之外,李润泽等[17]分析了凹凸前缘对压气机叶栅的气动性能提升,陈柳等[18]和杨伟琦[19]发现凹凸前缘具有良好的空化抑制作用,侯宇飞等[20]发现凹凸前缘对翼段动态失速有较好的控制效果。这表明基于仿生学的前缘凸起结构在流动控制领域引起了越来越广泛的关注与应用。

    尽管上述研究工作已对前缘凸起翼段的流动控制机理获得了一定的认识,但主要是基于数值模拟进行分析,并且没有深入地对流场分离现象进行系统的研究。本文采用PIV 实验手段,以NACA 634-021翼型作为研究对象,对平直翼段、单凸起和双凸起翼段的绕流流场进行高精度测量,系统地分析不同迎角下各翼段绕流流场的分离现象。

    实验在天津大学低湍流度回流式风洞中进行,风洞实验段长2300 mm,有效横截面积1000 mm × 1000 mm。风洞实验段来流速度可在2~60 m/s范围内连续调节,背景湍流度小于0.2%。实验对象为以NACA 634-021作为基准翼型的平直翼段,以及具有单个、双个前缘凸起的改型翼段(图1),其基准弦长c为200 mm,翼展S为1000 mm,凸起部分在基准翼型的前缘基础上外扩,每个凸起的轮廓线是一个余弦曲线的周期,波长λ为50 mm,幅值a为20 mm,单个、双个前缘凸起都位于翼展中心位置。

    图  1  翼段模型和测量截面
    Fig.  1  Airfoil model and measurement cross-section

    平直翼段测量翼展中心截面,单凸起翼段在凸起部分沿展向平均划分5个测量截面,双凸起翼段沿展向平均划分9个测量截面,相邻截面距离12.5 mm,将凸起中心处的翼段剖面(弦长最大的位置)称为凸峰截面,基准翼型和凸起交界处的翼段剖面称为根部截面。翼段展向沿竖直方向安装在风洞实验段底部的转盘上,可连续调节翼段迎角,翼段与转盘刚性连接,实验过程中翼段无振荡。激光沿水平方向照射,两台高速相机同步由上至下拍摄,实验装置如图2所示。

    图  2  实验装置图
    Fig.  2  Experimental setup diagram

    实验来流速度U为13.3 m/s,实验雷诺数为1.8 × 105,翼段迎角α在14°~24°范围内先增后减,每隔2°为一个测量工况。每个实验测量截面分为两个视场,视场1对应翼型前半部分,位置固定不动;视场2对应翼型后半部分,由于各工况下流动分离程度不同,所以视场2的位置随迎角变化发生调整。相机采用双帧模式拍摄,脉冲时间间隔为40 μs,每个工况拍摄两组数据,采样频率分别为200 Hz和800 Hz,采样时长分别为20.535 s和5.133 s,每组数据包含4107对图像。采用自适应互相关算法计算速度场,查询窗口大小为16像素 × 16像素,每个瞬时速度场包含159 × 99个速度矢量,相邻矢量间的实际距离约为1mm。

    通过同一瞬时时刻两个流场中出现的同一个涡结构,可以确定两个视场的相对坐标,将两个小视场进行拼接,能够获得更完整的大视场。由于两台相机到激光平面的拍摄距离不可能完全相同,所以两个视场中矢量点的间距存在细微的差别,因此在拼接前基于视场1的矢量点间距在视场2中重生成矢量点坐标,对视场2的数据进行插值,获得一个矢量点间距与视场1相同的新视场,再将这个新视场与视场1进行拼接。在拼接时两个小视场存在重叠的区域,拼接所得大视场中该区域内的数据取两个小视场数据的平均值。此外拼接得到的大视场中还存在两个无数据的矩形区域,这两个区域的数据用零填充。

    基于PIV获得的定量的速度场,图3中给出了平直翼段中心截面的流向平均速度Um云图,Inc表示增迎角工况,Dec表示减迎角工况,各工况拍摄时流场都已达到稳定状态。迎角递增工况中,14°迎角流场处于后缘分离状态,16°迎角时分离程度略有增大。当迎角增长到18°时,流场突然发生前缘分离,回流区域扩大至整个翼面,该流场状态一直持续至24°迎角。在迎角递减工况中,流场直至16°迎角时仍处于前缘分离状态,与增迎角工况对比具有明显的迟滞效应,直至14°迎角流场才恢复后缘分离状态。此现象表明该翼型属于尾缘-前缘混合型失速的翼型[18],其特点是在小迎角下先发生尾缘失速,当达到临界迎角时突然发生前缘分离,且在迎角减小时存在迟滞效应。

    图  3  平直翼段时均流向速度云图
    Fig.  3  Time-averaged flow velocity cloud diagram of the flat wing section

    单、双凸起翼段具有多个测量截面,为了便于观察整体的流动分离规律,给测得的二维速度场添加展向坐标,组合不同展向位置的测量截面,通过Tecplot软件生成三维速度场,并抽取流向平均速度等值面。 图4为单凸起翼段流向平均速度的等值面分布,参考平直翼段流线和流向速度云图,等值面取自由来流速度的0.3倍,该等值面与流线所显示的流动分离区域边界较为符合。云图显示的是使用来流速度无量纲化的湍动能$tke'$,其中湍动能表示为:

    图  4  单凸起翼段时均流向速度等值面
    Fig.  4  Mean flow velocity contour surface for the single leading-edge protuberance wing section
    $$ tke' = {{\frac{1}{2}\left( {{u_x}{'^2} + {u_{\textit{y}}}{'^2}} \right)} /{U_\infty }^2} $$ (1)

    式中:${u_x}'$表示流向脉动速度;${u_{\textit{y}}}'$表示法向脉动速度;${U_\infty }$表示自由来流速度。

    图中Um为流向平均速度,等值面底部为对应迎角的单凸起翼段模型。

    在迎角递增时,14°和16°迎角下流场呈现局部失速现象,各展向截面都处于后缘分离状态。16°迎角时,右侧根部截面分离程度最大,等值线显示分离点约在1/4弦长处,分离点后方高湍动能区逐渐增厚,说明流体剪切运动范围逐渐增大。18°迎角时流场出现单侧失速现象,右侧根部截面发生前缘分离,回流区域显著增大,高湍动能区集中在后缘部分,表明随着流动分离模式的切换,流场结构特征发生转变。此时凸起左侧截面仍为后缘分离状态,这一流场现象一直持续到22°迎角。

    当迎角增长到24°时,流场出现双侧失速现象,两侧截面均发生前缘分离,凸峰截面具有一股附着流动(如图6所示),且在分离点附近存在较强的脉动。此外,附着流动在向下游发展的过程中存在沿左侧的展向分量,在凸峰左侧截面中形成了流动再附着区域。此处要注意的是,由于凸峰左侧的两个截面的弦长小于凸峰截面,这两个截面在实验测量时无法拍摄到真实的前缘轮廓,只能拍摄到凸峰截面的前缘,双凸起翼段同样如此。减迎角工况下的流动分离规律与增迎角工况类似,流场失速状态切换的临界迎角与增迎角工况相同。

    图  6  凸峰及左侧截面的时均流向速度分布, α = 24°(Inc)
    Fig.  6  Mean flow velocity contour distribution of the left three sections

    这种特殊的分离规律与前缘凸起结构密切相关,目前的主流观点认为[21],前缘凸起具有涡流发生器的作用,其诱导产生一对反向旋转的流向涡对,在凸起根部产生上洗作用,增大了当地有效迎角,并将边界层内部的动量带入主流,降低了凸起根部区域的动量;在凸峰位置产生下洗作用,降低了当地有效迎角,同时将主流的动量带入边界层,补充了凸峰位置的动量。这一现象导致了迎角递增时根部截面最先发生前缘分离,并在前缘位置产生局部失速涡团,引起凸峰处的附着流动产生背离失速区域的展向分量,补充了对侧翼面的边界层动量,推迟了前缘失速的发生,直至迎角进一步增大,凸起两侧都发生前缘分离后,凸峰附着流动仍然存在。

    图  5  右侧根部截面的时均流向速度分布
    Fig.  5  Mean flow velocity contour distribution of the right root section

    双凸起翼段流向平均速度的等值面分布如图7所示,云图显示为无量纲湍动能$tke'$。双凸起翼段流场的演化规律与单凸起翼段相似,都具有从局部失速到单侧失速,再到双侧失速的转变过程。

    图  7  双凸起翼段时均流向速度等值面
    Fig.  7  Mean flow velocity contour surface for the double leading-edge protuberances wing section

    14°迎角时流场发生局部失速,右半部分凸起范围内的流场与单凸起翼段流场相同。但凸起数量的增多导致两侧凸峰诱导产生的流向涡在凹谷位置发生了叠加,大大减小了凹谷位置的边界层动量,形成一个局部失速涡团。迎角增大到16°后,流场发生单侧失速,右侧根部区域发生前缘分离,形成一个高压区,凸峰附着流动起到翼刀的作用,抑制了翼段右侧失速涡团的扩展,同时形成了较强的指向左侧的展向分量,补充了凹谷位置的能量,消除了局部失速涡团。迎角由20°先增后减至18°的过程中,流场始终处于双侧失速状态,凸峰附着流动始终偏向翼段左侧,表明流场结构受到流动历史的影响。此外,迎角减小时流场的双侧失速状态一直保持到18°,在16°迎角切换回单侧失速,表现出迟滞效应。

    对比发现,双凸起翼段三种失速模式转换的临界迎角小于单凸起翼段,同时凸起数量的增加使得凸峰附着流动的影响区域扩大,因此迎角较小时气动性能降低,迎角较大时气动性能提高。

    本征正交分解是由Lumley 等[22,23]提出的一种模态分解方法,被广泛地应用于复杂流场的降阶和流动结构提取。这种方法将流场分解成若干空间正交模态,按照各模态的能量大小进行排序,从而提取出流动主要模态。

    使用$N$个离散时刻的脉动速度场构建矩阵$U$:

    $$ U = \left[ {{u_1},} \right.{u_2}, \cdots \left. {,{u_N}} \right] $$ (2)

    计算脉动速度场协方差矩阵的特征值矩阵$\mu $和特征向量矩阵$A$为:

    $$ \left( {U{U^T}} \right)A = \mu A $$ (3)

    通过特征值与对应的特征向量构建POD模态,其中第$i$阶模态${\phi _i}$为:

    $$ {\phi _i} = \frac{{\sum\limits_{n = 1}^N {A_n^i{u_n}} }}{{\left\| {\sum\limits_{n = 1}^N {A_n^i{u_n}} } \right\|}},\left( {i = 1,2, \cdots ,N} \right) $$ (4)

    式中:$A_n^i$为第$i$阶特征值的特征向量的第$n$阶分量。各阶模态的时间系数${a_i}$表示为:

    $$ {a_i} = {\psi ^T}{u_n} $$ (5)

    式中:$\psi $为模态矩阵。

    图8给出了在不同工况下,平直翼段POD分解的各阶模态能量占比,不同的流动分离状态,对应的模态能量占比曲线规律也不同。平直翼段在迎角递增时,α = 14°和16°工况的第1阶模态能量占比分别为28.42%和26.71%,第2阶模态能量占比分别为4.49%和4.85%,而α ≥ 18°的各工况第1阶模态能量占比约为15%,第2阶模态能量占比约为8%。首阶模态中,后缘分离流场的模态能量远高于前缘分离流场;从第2阶模态开始直至第10阶模态,后缘分离流场的模态能量占比则明显小于前缘分离流场。这表明流动分离程度较小时,流场中存在一个能量较高的主要模态结构,其对流场的影响程度远大于其余高阶模态;而发生前缘分离时,各阶模态携带的能量差距相对较小,流场受到多个流动结构的共同影响。迎角递减过程中,处于迟滞迎角α = 16°的工况第1阶模态能量占比为18.06%,介于14°和18°迎角能量占比中间,高阶模态能量占比与前缘分离流场的数值相似,表明迟滞迎角时,流动结构更接近于前缘分离流场。

    图  8  平直翼段POD模态能量占比
    Fig.  8  The POD modes' energy proportion of the flat wing section

    图9为平直翼段前2阶模态结构(mode1~2)的流向脉动速度u′云图,流场处于后缘分离和前缘分离状态的模态结构明显不同。小迎角流场的第1阶模态结构呈现扁长形,起始于流动分离点,贴近翼段表面,分布于整个回流区内,表示了明显的负速度波动;第2阶模态结构具有法向高度差,表明回流区内存在大尺度的旋涡,但其能量占比远低于第1阶模态,说明流场发生后缘分离时的主要特征是回流区内的减速。

    图  9  平直翼段POD分解前2阶模态
    Fig.  9  First two POD modes of the flat wing section

    大迎角流场的第1阶模态结构法向高度较高,分布比较集中,表明在流动分离区和主流的交界处流向速度的变化较大;第2阶模态都包含一对沿流向分布且脉动速度方向相反的模态结构,表明流场存在沿流向的振荡,流动加速与减速交替出现。这说明流场发生前缘分离时的特点是分离区和主流区交界处明显的速度波动。

    单凸起翼段在每个迎角下都测量了5个展向截面,取各个截面POD分解的模态能量占比的平均值,作为对应迎角下单凸起翼段整体的模态能量占比,如图10所示。迎角增减过程中,处于局部失速状态(α = 14°和16°)的流场前3阶模态能量相对较高,第1阶能量占比都超过了20%,处于单侧和双侧失速迎角范围的流场前3阶模态能量占比相对较小,这与平直翼段的规律相似。但α = 22°工况的模态能量较为突出,此时流场处于两种失速状态转变的临界迎角,流场存在不稳定性[18]

    图  10  单凸起翼段POD模态能量占比
    Fig.  10  The POD modes' energy proportion of the single leading-edge protuberance wing section

    图11为单凸起翼段不同工况下前2阶模态结构的等值面分布,流场的三种失速状态对应了不同特征的模态结构。第1阶模态中,α = 14°和16°工况流场发生局部失速,各截面都处于后缘分离状态,模态结构的起始位置与翼段分离线的趋势相符。迎角增长到18°后流场单侧失速,凸峰两侧的分离程度不同,但各截面模态结构非常接近,与平直翼段处于前缘分离状态下的模态类似。这表明单侧失速时,虽然凸峰附着流动补充了左侧翼面的动量,抑制了前缘分离的发生,但流场仍然具有前缘分离的结构特征。当迎角α = 24°时,附着流动已经不足以补充足够的能量,左侧翼面也发生了前缘分离,流场中的模态结构分为两部分,其中下游结构和单侧失速时相同,印证了单侧失速时对流场特征的判断;上游结构位于凸起左侧根部区域,与下游结构方向相反,这表明当下游流动减速时,凸起左侧区域会出现流动加速,可能引起流动分离模式的变化。第2阶模态中,α = 14°和16°流场包含的两部分模态结构流向间距较短,18°发生单侧失速后间距变大,表明随着流场转变为前缘分离,波动尺度范围增大。α =22°时流场处于分离模式转变的临界迎角,右侧根部截面模态结构发生突变,对应了该工况POD分解模态能量占比的异常。因此,有必要对临界迎角附近的流场进行深入分析。

    图  11  单凸起翼段POD分解前2阶模态
    Fig.  11  First two POD modes of the single leading-edge protuberance wing section

    POD分解得到的时间系数值表征了该阶模态结构在不同时刻对瞬时流场的影响程度大小。图12(a)为24°迎角时单凸起翼段5个截面的第1阶模态时间系数分布,其特点是系数围绕零值对称波动。若某时刻的时间系数绝对值较大,表明该瞬时流场受模态结构影响程度较高,时间系数处于不同水平的流场会存在差异,在时间系数围绕其均方根值发生波动时,对应流场的流动分离程度也围绕平均流场发生波动。通过对比发现,当选定时间系数小于均方根值的负1.8倍时,流场呈现单侧失速现象,与之相反,当时间系数大于均方根值的正1.8倍时,流场仍为双侧失速状态,如图12(b)所示。图中是对所有时间系数符合条件的瞬时流场做平均后的统计结果,各个截面被选出的瞬时流场数量占总体的比例标注在图12(a)中。其中单侧失速流场平均占比为3.7%,同时伴随着尾流区域的流动减速,分离面上抬;双侧失速流场平均占比为3.3%,相比于整体的时均流场其左侧前缘分离程度更大,凸峰附着流动更偏向下游的位置,并在尾流区发生流动加速,分离面更接近翼面。其余93%的流场形态介于二者之间,表明 24°迎角下的流场在单侧失速和双侧失速状态间发生波动,且绝大多数时刻发生双侧失速。同时要注意的是,各展向截面的实验测量是独立进行的,图12(b)所示的流场是统计的结果,仅代表某些时刻下可能出现的流态,时间系数的持续波动表明流场状态在不断地发生变化。

    图  12  单凸起翼段流场对比,α =24°
    Fig.  12  Flow field comparison of the single leading-edge protuberance wing section, α =24°

    22°迎角处于分离模式转换的临界状态,流场存在三种不同的流态,如图13(a)所示。增迎角工况时,截面1和2的时间系数都出现了较稳定的分段,两部分各自围绕不同的系数值波动,表明流场在实验过程中出现了两种不同的分布形态,且都能持续一定时间,如图13(b)、(c)所示。根据上述分析,当某个截面在前缘分离后发生流动再附着时,表明凸峰附着流动具有指向该侧截面的展向分量,则凸起的另一侧必然处于前缘分离状态。因此,通过单个截面的流动形态,可以推测整体翼段的分离模式。当截面1出现流动再附着时,对应时间系数占比为56.3%,表明在超过一半时间内,流场处于双侧失速且附着流动指向左侧;当截面2出现流动再附着时,系数占比为24.9%,表明流场双侧失速且附着流动指向右侧。减迎角工况时,截面2的时间系数发生了多次分段,说明流场状态在短时间内发生连续变化,其中双侧失速状态发生两次,合计占比为29.8%;单侧失速状态发生三次,合计占比为70.2%。

    图  13  单凸起翼段流场对比,α =22°
    Fig.  13  Flow field comparison of the single leading-edge protuberance wing section, α =22°

    对双凸起翼段每个迎角下的9个展向截面POD分解的各阶模态能量占比取平均值,得到对应迎角下双凸起翼段整体的模态能量占比,如图14所示。

    图  14  双凸起翼段POD模态能量占比
    Fig.  14  The POD modes' energy proportion of the double leading-edge protuberances wing section

    在增迎角工况中,流场处于单侧失速状态(α = 16°和18°)的第1阶模态能量占比较低,高阶模态能量变化均匀;局部和双侧失速流场第1阶模态能量更高,其中双侧失速流场的能量占比下降迅速,局部失速流场则相反。这表明局部失速和双侧失速的流场具有能量更集中的主要模态结构。α = 18°(Dec)工况的流场处于迟滞区,在迎角递减时能量占比曲线符合双侧分离流场的规律。

    图15为双凸起翼段不同工况下前2阶模态结构的等值面分布。当流场处于局部失速状态时,右半部分凸起区域的模态结构和单凸起翼段相似,左半部分凸起区域的模态结构对应了凹谷处的局部失速涡团。流场发生单侧失速后,随着局部失速涡团的消失,前两阶模态结构都与单凸起翼段单侧失速时相一致,流场具有前缘失速的特点,同时模态结构的不对称形状,反映了附着流动向下游发展时指向左侧的展向分量。双侧失速时前两阶模态结构都呈现不规则形状,各展向截面之间差距较大。第1阶模态结构在左半部分凸起区域流向分布较长,表明附着流动始终对左侧区域施加影响。第2阶模态结构更加破碎,表明流场具有很多小尺度的失速涡团,流动带有很强的三维性。

    图  15  双凸起翼段POD分解前2阶模态
    Fig.  15  First two POD modes of double leading-edge protuberances wing section

    为了验证双凸起翼段是否存在类似的分离不稳定现象,分别取各截面POD分解的时间系数均方根值的正负1.8倍作为阈值,对筛选出的瞬时流场做平均,统计结果如图16所示。当模态结构的影响相反时,流场形态也存在差异。在时间系数小于均方根值的负1.8倍时,小迎角下凹谷处的局部失速涡团相对较大,具有向下游发展的趋势,而大迎角下的凸峰附着流动相对更加集中。反之当时间系数大于均方根值的正1.8倍时,小迎角下局部失速涡团减小,大迎角下凸峰附着流动扩散。在临界迎角下,双凸起翼段保持了良好的分离稳定性,失速模式没有发生波动,表明凸起数量的增加提高了翼段的流动稳定性。

    图  16  双凸起翼段流场对比
    Fig.  16  Flow field comparison of double leading-edge protuberances wing section

    采用PIV实验测量了平直翼段、单凸起和双凸起翼段的瞬时流场,对比了各工况下的流动分离规律,并基于POD方法分析了流场的流动特性,结论如下:

    1)平直翼段达到失速迎角后直接发生前缘分离,单凸起和双凸起翼段都具有独特的两步失速特性,随着迎角增大,流场先由局部失速状态转变为单侧失速状态,再进一步转变为双侧失速状态,其中双凸起翼段失速模式转换的临界迎角小于单凸起翼段,但发生双侧失速后的流动附着区域大于单凸起翼段。

    2)POD模态能量占比表明,平直翼段和单凸起翼段分别发生后缘分离和局部失速时,流场中存在能量较高的主要模态结构,其对流场的影响程度大于其余高阶模态;两种翼段分别发生前缘分离和单侧失速后,各阶模态携带的能量差距相对较小,流场受到多个流动结构的共同影响。双凸起翼段流场具有对称性,在发生局部失速和双侧失速时首阶模态能量更高,单侧失速时能量更低,说明局部失速时凹谷处的失速涡团和双侧失速时两股凸峰附着流的汇合都导致了流动能量的集中。

    3)平直翼段的POD模态结构表明,流场发生后缘分离时的主要特征是回流区内的减速,发生前缘分离时分离区和主流区交界处存在明显的速度波动;前缘凸起翼段在发生单侧失速后,各截面的分离模式不同,但都具有前缘分离流场的特征;发生双侧失速后,双凸起翼段流场具有很多小尺度的失速涡团,流动带有很强的三维性。

    4)POD时间系数分析显示,单凸起翼段处于单侧失速与双侧失速转换的临界迎角时,流动分离模式具有不稳定性;凸起数量的增加提高了翼段的流动稳定性。

  • 图  1   翼段模型和测量截面

    Fig.  1   Airfoil model and measurement cross-section

    图  2   实验装置图

    Fig.  2   Experimental setup diagram

    图  3   平直翼段时均流向速度云图

    Fig.  3   Time-averaged flow velocity cloud diagram of the flat wing section

    图  4   单凸起翼段时均流向速度等值面

    Fig.  4   Mean flow velocity contour surface for the single leading-edge protuberance wing section

    图  6   凸峰及左侧截面的时均流向速度分布, α = 24°(Inc)

    Fig.  6   Mean flow velocity contour distribution of the left three sections

    图  5   右侧根部截面的时均流向速度分布

    Fig.  5   Mean flow velocity contour distribution of the right root section

    图  7   双凸起翼段时均流向速度等值面

    Fig.  7   Mean flow velocity contour surface for the double leading-edge protuberances wing section

    图  8   平直翼段POD模态能量占比

    Fig.  8   The POD modes' energy proportion of the flat wing section

    图  9   平直翼段POD分解前2阶模态

    Fig.  9   First two POD modes of the flat wing section

    图  10   单凸起翼段POD模态能量占比

    Fig.  10   The POD modes' energy proportion of the single leading-edge protuberance wing section

    图  11   单凸起翼段POD分解前2阶模态

    Fig.  11   First two POD modes of the single leading-edge protuberance wing section

    图  12   单凸起翼段流场对比,α =24°

    Fig.  12   Flow field comparison of the single leading-edge protuberance wing section, α =24°

    图  13   单凸起翼段流场对比,α =22°

    Fig.  13   Flow field comparison of the single leading-edge protuberance wing section, α =22°

    图  14   双凸起翼段POD模态能量占比

    Fig.  14   The POD modes' energy proportion of the double leading-edge protuberances wing section

    图  15   双凸起翼段POD分解前2阶模态

    Fig.  15   First two POD modes of double leading-edge protuberances wing section

    图  16   双凸起翼段流场对比

    Fig.  16   Flow field comparison of double leading-edge protuberances wing section

  • [1]

    REIDENBERG J S, LAITMAN J T. Blowing bubbles: an aquatic adaptation that risks protection of the respiratory tract in humpback whales (Megaptera novaeangliae)[J]. The Anatomical Record, 2007, 290(6): 569–580. doi: 10.1002/ar.20537

    [2]

    HAIN J, CARTER G, KRAUS S, et al. Feeding behavior of the humpback whale, Megaptera novaeangliae, in the western North Atlantic[J]. Fishery Bulletin, 1981, 80: 259–268.

    [3]

    FISH F E, BATTLE J M. Hydrodynamic design of the humpback whale flipper[J]. Journal of Morphology, 1995, 225(1): 51–60. doi: 10.1002/jmor.1052250105

    [4]

    FISH F E, WEBER P W, MURRAY M M, et al. The tubercles on humpback whales' flippers: application of bio-inspired technology[J]. Integrative and Comparative Biology, 2011, 51(1): 203–213. doi: 10.1093/icb/icr016

    [5]

    MIKLOSOVIC D S, MURRAY M M, HOWLE L E, et al. Leading-edge tubercles delay stall on humpback whale (Megaptera novaeangliae) flippers[J]. Physics of Fluids, 2004, 16(5): L39–L42. doi: 10.1063/1.1688341

    [6]

    PEDRO H C, KOBAYASHI M. Numerical study of stall delay on humpback whale flippers[C]//Proc of the 46th AIAA Aerospace Sciences Meeting and Exhibit. 2008. doi: 10.2514/6.2008-584

    [7]

    WEBER P W, HOWLE L E, MURRAY M M, et al. Computational evaluation of the performance of lifting surfaces with leading-edge protuberances[J]. Journal of Aircraft, 2011, 48(2): 591–600. doi: 10.2514/1.C031163

    [8]

    JOHARI H, HENOCH C, CUSTODIO D, et al. Effects of leading-edge protuberances on airfoil performance[J]. AIAA Journal, 2007, 45(11): 2634–2642. doi: 10.2514/1.28497

    [9]

    ZHANG M M, WANG G F, XU J Z. Aerodynamic control of Low-Reynolds-Number airfoil with leading-edge protuberances[J]. AIAA Journal, 2013, 51(8): 1960–1971. doi: 10.2514/1.J052319

    [10]

    CAI C, LIU S H, ZUO Z G, et al. Experimental and theoretical investigations on the effect of a single leading-edge protuberance on airfoil performance[J]. Physics of Fluids, 2019, 31(2): 027103. doi: 10.1063/1.5082840

    [11]

    CAI C, ZUO Z G, MAEDA T, et al. Periodic and aperiodic flow patterns around an airfoil with leading-edge protuberances[J]. Physics of Fluids, 2017, 29(11): 115110. doi: 10.1063/1.4991596

    [12] 李德友, 常洪, 左志钢, 等. 仿座头鲸前缘凸起数对NACA 634-021翼型失速控制机理影响研究[J]. 大电机技术, 2020(4): 1–9. DOI: 10.3969/j.issn.1000-3983.2020.04.001

    LI D Y, CHANG H, ZUO Z G, et al. The effect of humpback whalelike leading-edge protuberances number on stall control mechanism of NACA634-021 airfoil[J]. Large Electric Machine and Hydraulic Turbine, 2020(4): 1–9. doi: 10.3969/j.issn.1000-3983.2020.04.001

    [13]

    ZHAO M, ZHAO Y J, LIU Z X, et al. Proper orthogonal decomposition analysis of flow characteristics of an airfoil with leading edge protuberances[J]. AIAA Journal, 2019, 57(7): 2710–2721. doi: 10.2514/1.J058010

    [14]

    YANG P X, ZHU Y C, WANG J J. Effect of leading-edge tubercles on the flow over low-aspect-ratio wings at low Reynolds number[J]. Theoretical and Applied Mechanics Letters, 2023, 13(1): 100386. doi: 10.1016/j.taml.2022.100386

    [15]

    GHAZI M A, GHASSEMI H, GHAFARI H R. Hydrodynamic performance and the power production of the tidal turbine by new profiles at the leading-edge tubercles[J]. Ships and Offshore Structures, 2023, 18(6): 875–885. doi: 10.1080/17445302.2022.2078532

    [16]

    DONG J Z, CHU W L, ZHANG H G, et al. Flow control mechanism of compressor cascade: a new leading-edge tubercles profiling method based on sine and attenuation function[J]. Physics of Fluids, 2023, 35(6): 066118. doi: 10.1063/5.0151476

    [17] 李润泽, 王琦, 徐宁, 等. 凹凸前缘压气机叶栅旋涡结构探索[J]. 热能动力工程, 2021, 36(9): 69–78. DOI: 10.16146/j.cnki.rndlgc.2021.09.009

    LI R Z, WANG Q, XU N, et al. Exploration on the vortex structure of concave-convex leading edge compressor cascade[J]. Journal of Engineering for Thermal Energy and Power, 2021, 36(9): 69–78. doi: 10.16146/j.cnki.rndlgc.2021.09.009

    [18] 陈柳, 曹琳琳, 赵国寿, 等. 仿生前缘流动与空化控制机理的数值研究[J]. 工程热物理学报, 2019, 40(10): 2291–2298. DOI: CNKI:SUN:GCRB.0.2019-10-014

    CHEN L, CAO L L, ZHAO G S, et al. Numerical study on the mechanism of flow and cavitation control by leading-edge tubercles on hydrofoil[J]. Journal of Engineering Thermophysics, 2019, 40(10): 2291–2298. doi: CNKI:SUN:GCRB.0.2019-10-014

    [19] 杨伟琦. 仿座头鲸胸鳍翼段空化流动控制机理研究[D]. 哈尔滨: 哈尔滨工业大学, 2021.

    YANG W Q. Study on the control mechanism of cavitation flow in the pectoral fin segment of humpback whale[D]. Harbin: Harbin Institute of Technology, 2021. doi: 10.27061/d.cnki.ghgdu.2021.002838

    [20] 侯宇飞, 李志平. 仿生正弦前缘对翼面动态失速的影响[J]. 航空学报, 2020, 41(1): 123276. DOI: 10.7527/S1000-6893.2019.23276

    HOU Y F, LI Z P. Effect of bionic sinusoidal leading-edge on dynamic stall of airfoil[J]. Acta Aeronautica et Astronautica Sinica, 2020, 41(1): 123276. doi: 10.7527/S1000-6893.2019.23276

    [21] 蔡畅. 仿座头鲸鳍肢前缘凸起对翼型失速特性控制机理研究[D]. 北京: 清华大学, 2018.

    CAI C. Effect of leading-edge protuberances inspired by humpback whale flipper on airfoil stall control[D]. Beijing: Tsinghua University, 2018. doi: 10.27266/d.cnki.gqhau.2018.000320

    [22]

    LUMLEY J L. The structure of inhomogeneous turbulence[J]. Atmospheric Turbulence and Radio Wave Propagation, 1967: 166–178. doi: 10.1007/BF00271656

    [23]

    BERKOOZ G, HOLMES P, LUMLEY J L. The proper orthogonal decomposition in the analysis of turbulent flows[J]. Annual Review of Fluid Mechanics, 1993, 25: 539–575. doi: 10.1146/annurev.fl.25.010193.002543

图(16)
计量
  • 文章访问数:  94
  • HTML全文浏览量:  33
  • PDF下载量:  23
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-12-25
  • 修回日期:  2024-03-18
  • 录用日期:  2024-03-21
  • 网络出版日期:  2024-05-30

目录

/

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