张奕，潘翀,2，窦建宇，张淼

1.北京航空航天大学 航空科学与工程学院，北京 100191

2.北京航空航天大学 宁波创新研究院先进飞行器与空天动力技术创新研究中心，宁波 315800

3.中国商用飞机有限责任公司 上海飞机设计研究院，上海 201210

0 引 言

附着在交通工具表面的边界层流态大部分为湍流，其所产生的摩阻较层流边界层高出约一个量级。通过控制湍流边界层内多尺度流动结构可有效减小边界层的摩阻并实现总阻力的降低，对节能减排具有重要意义。

随着雷诺数的升高，湍流边界层内小尺度流动结构对湍动能和摩阻的贡献逐渐降低，以近壁面小尺度结构为主要控制对象的传统主/被动控制技术的减阻效果随之降低。近期对湍流边界层内跨流层跨尺度湍流结构相互作用的研究表明，外区大尺度流动结构影响着湍流边界层内跨流层能量传递和近壁雷诺切应力的生成[1-2]。因此，通过控制较高流层中的大尺度流动结构，既可直接降低其在壁面留下的大尺度摩阻成分，亦可抑制其对近壁小尺度结构的调制作用，降低小尺度结构对壁面摩阻的贡献，进而有望在高雷诺数条件下降低壁面摩阻[3-5]。

目前针对湍流边界层的减阻流动控制技术包括沟槽壁面[6-9]、超疏水壁面[10]、局部吹吸气[11]、近壁面行波[12]和涡流发生器（Vortex Generator,VG）[13-15]等。其中，VG 具有结构简单、附加重量低、易于维护等优点，工程应用潜力较大。VG 的特征高度h 与边界层厚度δ相当（h/δ～1），存在附加压差阻力大的缺点。与之相比，微型涡流发生器（Micro Vortex Generator,MVG）几何尺寸更小（h/δ< 0.5）、附加阻力更低[16-17]，具有比常规VG 更好的减阻控制能力[13,17-18]。因其特征高度h 的差异，MVG 也被称为低形阻VG（Low-Profile VGs,LPVG）[18-20]、亚边界层VG（Sub-Boundary-Layer VGs,SBVG）[21]、浸没式VG（Submerged VGs）[16]等。MVG 有叶片形、楔形、车辙形及叉骨形等多种形状。当来流绕过MVG时，在其下游会形成较为稳定的时均流向涡对，从而对近壁流动施加控制。

目前，对MVG 流动控制技术的研究主要集中于3 个方向：1）分析不同形式MVG 在零压梯度或逆压梯度边界层中诱导产生的流动结构及其沿程演化规律[19,21-22]；2）应用MVG 进行边界层分离流动控制[18,20-21]；3）将MVG 应用于二维翼型、三维机翼、机身等处进行流动减阻试验[20-21,23-26]。为解决MVG 在湍流边界层减阻控制中的工程应用问题，需明确MVG 对流场产生有效控制的流向范围。文献[20]中提到，当基于边界层动量厚度θ的特征雷诺数Reθ≈ 3.5 × 105时，MVG 产生的时均流向涡对沿程抬升并在其下游约50h 处耗散殆尽。Zaman 等[22]的研究表明，当Reθ≈ 1.1 × 104时，MVG 所产生的流向高、低速流动区域沿主流方向逐渐融合，并在其下游约24h 处形成较厚的均匀低速流层。

本文针对上述第一个研究方向开展实验测量，旨在研究MVG 对经典零压梯度湍流边界层不同下游位置流动特性和摩阻的影响，为掌握MVG 控制效果的沿程演化规律提供新的认识，为MVG 在中高雷诺数湍流边界层中的减阻流动控制应用提供机理支撑。为此，本文使用体视粒子图像测速（Stereoscopic Particle Image Velocimetry,SPIV）技术和免标定双层热膜摩阻传感器，在低速风洞中测量中等雷诺数条件下平板湍流边界层内单排楔形MVG 阵列下游不同流向站位的法向-展向平面内的三分量速度场和当地壁面摩阻，使用统计分析、能谱分析和本征正交分解（Proper Orthogonal Decomposition,POD）等方法研究MVG 对湍流边界层内流场统计特性和流动结构的影响，并讨论MVG 下游摩阻的沿程变化及MVG 几何参数对减摩阻能力的影响。

1 实验设备与研究方法

1.1 实验模型布置

风洞实验在北京航空航天大学D6 低速风洞中进行。该风洞为直流吸气式，风速可调范围为0～50 m/s。实验段总长度为5 m，收缩比约为9∶1，采用变截面设计以减小洞壁体发展边界层对主流区流场均匀性的影响，其入口、出口横截面尺寸分别为0.4 m × 0.5 m、0.48 m × 0.5 m。

实验模型布置如图1 所示。在风洞实验段内放置长4.5 m、厚0.015 m、宽0.40～0.48 m 的有机玻璃平板，用于产生零压力梯度平板湍流边界层。平板距风洞底壁0.085 m，平板前缘距风洞实验段入口约0.1 m，前缘具有9.5°切角以防止流动在前缘产生局部分离。在距平板前缘0.300 m 的下游位置，沿展向布置一根直径1 mm 的圆柱形碳棒作为绊线，以触发边界层转捩。本文所使用的笛卡尔坐标系原点位于平板展向中心处MVG 后缘与壁面的接触点，x、y、z 分别代表流向、法向和展向，各方向的速度、时均速度和脉动速度分别记为（U、V、W）、（、、）和（u、v、w）。U∞为自由来流速度。如无特殊说明，在下文光滑壁面湍流和MVG 壁面湍流各工况结果中，具有上标“+”的速度和长度单位均已使用光滑壁面湍流边界层的内尺度（摩擦速度uτ及涡黏尺度l*）作无量纲化。

图1 低速风洞中实验布置示意图Fig.1 Schematic diagram of the experimental setup in the test section of a low-speed wind tunnel

在距绊线3.1 m 的下游位置，沿展向布置一排3D 打印加工的MVG 阵列。边界层在此流向位置已经完全转捩为湍流边界层。表1 为自由来流速度U∞=14 m/s 时该处湍流边界层的主要特征参数，其中Reτ为摩擦雷诺数（Reτ=uτδ/ν，ν为运动黏性系数），Reθ为动量厚度雷诺数，H 为形状因子，Δy+、Δz+为空间分辨率，uτT/δ为以边界层外尺度无量纲的采样总时长（T 为各工况实验的拍摄时间，用于衡量计算结果是否收敛）。上述特征参数系利用Chauhan 等[27]提出的复合速度型公式对SPIV 所测流向时均速度型函数(y+)进行拟合得到，综合考虑了内区及尾迹区内的流向时均速度的分布特点。

表1 光滑壁面湍流边界层主要特征参数Table 1 Main characteristic parameters of the studied smooth-wall TBL

1.2 MVG 尺寸参数

MVG 阵列特征尺寸如图2 所示，其中l 为流向长度，a 为展向宽度，s 为展向间距。本文以Ashill等[21]提出的前向三角楔形MVG（Forwards Wedge MVG）作为基准模型。该构型的MVG 在已有研究中表现出最小的附加阻力及较优的控制性能[21,28]。本文中，基准楔形MVG 记作MVG0，其特征高度h0=5 mm、展向间距s0=0 mm（即相邻MVG 无间隙）、展向宽度a0=10 mm、流向长度l0=20 mm。下文中的部分空间坐标将使用MVG0的特征高度h0进行无量纲化。此外，本文还对比研究了具有更大展向间距的MVG 阵列（记作MVGs）及更高特征高度的MVG 阵列（MVGh）。3 种MVG 阵列的主要几何参数如表2 所示。

表2 3 种MVG 阵列的主要几何参数Table 2 Main size of the three MVG array

图2 MVG 阵列特征尺寸示意图Fig.2 Schematic diagram of MVG array characteristic size

1.3 SPIV 测量

当U∞=14 m/s 时，对3 种MVG 阵列下游x/h0=10（近尾迹区）和x/h0=30（远尾迹区）这2 个流向站位开展SPIV 测量，获得该x 站位y-z 平面内的三分量速度场，并与无MVG 条件下的光滑壁面流场测量结果进行对比。SPIV 实验所使用的示踪粒子为平均粒径约2 µm 的癸二酸二异辛酯（DEHS）液滴，通过网格状的分布式粒子布撒装置注入直流式风洞的进气口，被吸入风洞后流经收缩段、稳定段与主流充分掺混。双曝光脉冲激光器（Beamtech Vlite-380）激发能量为380 mJ/pulse的激光束，经透镜组整形后形成片光，从风洞侧壁垂直入射照亮测量平面。片光厚度为3 mm，以降低流向速度的测量不确定度。2 台超高分辨率跨帧CCD 相机（Imperx CLB-B6620M，分辨率为6 600像素 × 4 400 像素）搭配180 mm 微距镜头（TAMRON SP AF 180 mm F/3.5 MACRO1∶1）记录测量平面内的示踪粒子信息，其物理分辨率是0.016 mm/像素。2 台相机分别布置于测量平面的上、下游两侧，以异侧布局的形式拍摄，相机透视角约为65°。使用经GPU 加速的光流法速度场求解器（Multiiterative Lucas-Kanade,MILK）[29]分别计算2 台相机像平面内的粒子图像，获得二维二分量速度场，计算时对粒子图像进行了4 次降采样，最终查询窗口尺寸为16 像素 × 16 像素，重叠率为75%，最大跨帧位移不超过12 像素。之后，将2 台相机像平面上的二维二分量速度场通过空间标定映射至测量平面，获得二维三分量速度场[30-32]。SPIV 测量的有效视野范围为80 mm（法向）× 60 mm（展向）。示踪粒子以14 m/s 的速度在跨帧时间20 µs 内沿流向最多前进2.8 mm（小于片光厚度3 mm），可保证在跨帧粒子图像对上计算出互相关峰。各工况均采样1 000对瞬时速度场，相机的采样频率为0.5 Hz。经检验，SPIV 测量得到的湍流二阶统计量在此样本量下达到收敛。

应用MILK 算法得到粒子位移的不确定度约为0.1 像素[29]，根据陈启刚和钟强[30]给出的估算方法，面内速度分量（V、W）的测量精度约为0.4%，面外速度分量（U）的测量精度约为1.0%。图3 为SPIV 测量得到的光滑壁面湍流边界层的流向时均速度型函数(y+)、雷诺正应力曲线（(y+)、(y+)、(y+)）和雷诺切应力曲线（(y+)），同时给出Watanabe 等[33]在Reτ=3 058 条件下开展时间发展湍流边界层直接数值模拟（Direct Numerical Simulation,DNS）得到的结果作为对比，图中κ为卡门常数，(i=j=1,2,3)中u1、u2、u3分别表示脉动速度u、v、w。可以看到，SPIV 实验测量得到的在近壁区y+=10 处仍然能够和DNS 的结果相吻合；在雷诺应力方面，实验测量的、及沿法向的变化趋势和DNS 的结果吻合，两者的数值在y+=30 以上的流层基本一致；但在y+=10～20 以内（内区）总体低于DNS 的结果，究其原因，一方面是本文研究的体发展湍流边界层和DNS 模拟的时间发展边界层[33]在流动特性上存在差异，另一方面则是本文采用的SPIV 实验布置方式导致测量面外的U 分量不确定度更大，且较高的片光厚度使得沿片光厚度方向的体平均效应较强，进一步降低了对U 分量小尺度脉动的解析能力，从而使内区测量结果相较于DNS 结果偏低。

图3 光滑壁面工况下SPIV 测量结果与 DNS 结果对比Fig.3 The SPIV measurement results are compared with the DNS results under the smooth wall condition

1.4 壁面摩阻测量

湍流边界层的壁面摩阻由高南和刘玄鹤[34]发展的一种新型免标定双层热膜摩阻传感器测量得到。该传感器装置如图4（a）所示，其工作原理为：在外接恒温开尔文电桥的控制下，双层热膜（镍箔）同温协同工作，此时下膜的存在将上膜向壁面的传热损失限制在总产热量的5%以内，使上膜热量几乎只传给流体，解决了传统单层热膜摩阻测量技术因壁面热损失而导致的测量偏差问题。Liu 等在此基础上发展了根据上层热膜发热量计算壁面摩阻的免标定测量技术，其计算原理详见文献[35-37]。本文使用的矩形镍箔流向长度为1.75 mm、展向宽度为30 mm、厚度约为40 µm，可测量矩形镍箔表面的平均摩阻。传感器电桥所测得的瞬时电压数据由数据采集卡（NI USB-6259）实时采集，采样频率为1 024 Hz。

图4 摩阻传感器及其测量结果Fig.4 The friction sensor and its measurement results

壁面时均摩阻τw理论值[38]的半经验关系式为：

式中：Rex表示以距边界层起点长度（绊线位置）为特征尺寸的雷诺数。

对各工况均开展7 次摩阻重复测量实验，多次测量所得不确定度用误差条表示。摩阻单次测量持续时间为40 s（满足Ashill 等[21]建议的收敛性准则）。如图4（b）所示，在不同自由来流速度下利用摩阻传感器测得的光滑壁面湍流边界层的时均摩阻τw（红色曲线）与式（1）估计的理论值（黑色曲线）相比较，二者最大误差约为4.7%。图4（b）中也给出了U∞=14 m/s 时经SPIV 所测(y+)估算所得的τw（蓝色点），比摩阻传感器的测量结果高14%。此外，本文还利用摩阻传感器测量了U∞=14 m/s 时光滑壁面及3 组MVG 阵列下游x/h0=10～80 范围内的壁面摩阻，研究发现：即便是大间距的MVGs阵列，摩阻传感器的展向测量尺寸仍然包含了2 个MVG 的展向周期尺寸，这意味着摩阻传感器获得的是具有展向平均意义的当地摩阻。

2 MVG 下游速度场统计量

图5 给出了3 种MVG 阵列下游近、远尾迹区2 个x 站位（x/h0=10 和30）的y-z 平面内的时均速度场。彩色云图表示垂直平面的流向速度分量，黑色箭头表示法向速度分量和展向速度分量合成的面内速度矢量。蓝色虚线为MVG 阵列下游时均速度场=17 等值线，红色点划线为光滑壁面工况下=17 等值线。可以看到，在3 种MVG 阵列的下游近尾迹区，边界层近壁区均出现时均流向涡对和时均流向速度亏损区，时均流场的展向不均匀性显著，与前人研究相符[21-22]；时均流向涡对结构强度及时均流向速度亏损区沿程逐渐衰减，在下游远尾迹区仅表现为近壁处的低动量层，展向不均匀性明显减弱。相较于基准MVG0工况（图5（b）），大间距MVGs工况（图5（a））在远尾迹区形成的近壁低动量层更薄、更接近光滑壁面水平；而更高的MVGh（图5（c））则显著增厚了远尾迹区的近壁低动量层。

图5 MVG 阵列下游的时均速度场Fig.5 Time averaged results of three velocity component fields behind MVG arrays

考虑到MVG 引入的展向不均匀性，对流场统计量沿展向进行了平均。图6 给出了3 种MVG 阵列的近、远尾迹区（x/h0=10 和 30）时均流向速度型和雷诺应力型的展向平均结果。图中，“〈〉”表示沿时间和展向的系综平均，蓝色和红色的点划线分别为MVG 阵列h=5、10 mm 所对应的法向位置。从图6（a）可以看到，MVG 阵列对近尾迹区（x/h0=10）流场统计量的影响主要体现在3 个方面：1）在MVG 高度的流层（y=h）附近，产生低平均动量及高湍流脉动区，主要由MVG 诱导产生的时均流向涡对和时均流向速度亏损区所致；2）在MVG 高度以下的近壁区（y < h），雷诺正应力（即湍流脉动强度）及雷诺切应力均有所降低，可归因于时均流向涡对近壁湍流结构的规则化作用，在一定程度上反映了MVG 的减摩阻潜力；3）MVG 的h 越高、s 越小，对边界层的扰动作用越强，诱导产生的时均流向涡对越强。从图6（b）可以看到，在远尾迹区（x/h0=30），MVG 对平均流场的诱导作用随时均流向涡对的耗散而减弱，但仍然在湍流脉动强度场上留下峰值“印记”，且峰值出现的流层相比近尾迹区更高，反映了时均流向涡对从壁面逐渐抬升的bottom-up 过程[21,39]。值得注意的是，在x/h0=30 处，3 种MVG阵列工况下近壁区（y+< 15）的流向雷诺正应力和雷诺切应力均低于光滑壁面湍流边界层的水平，降低的幅度仍和MVG 高度及间距相关，预示MVG 在远尾迹区具有减摩阻潜力。

图6 MVG 阵列下游的时均流向速度型与雷诺应力的展向平均结果Fig.6 Spanwise average results of wall-normal profiles of time-averaged streamwise velocity and Reynolds stress behind MVG arrays

3 展向能谱特性

为从尺度空间上厘清MVG 对湍流脉动特性的影响，对SPIV 测量得到的速度场进行展向能谱分析。图7 给出了MVG0阵列近、远尾迹区（x/h0=10 和30）u、v、w 分量的展向预乘能谱及u 和v 的互相关展向预乘能谱，即（i=j=1,2,3 或i=1,j=2），同时给出光滑壁面的情况作为对比。各预乘能谱均使用所对应的光滑壁面工况的能谱峰值进行无量纲化，即=kzϕuiuj/。其余工况（MVGh和MVGs）的能谱特性与MVG0工况类似，受限于篇幅不再给出。如图7（a）所示，SPIV测量能够解析到光滑壁面能谱的内区峰P1和外区峰P2，峰值的法向高度和峰值波长均与前人的DNS 结果[40]一致。此外，、的内区峰P1和的外区峰P2也和DNS 结果[40]吻合，但现有的SPIV 测量难以解析到在y+≈ 20 处的内区峰。图7（b）表明，在MVG 作用下，近尾迹区、、和的能谱中均出现一个独立的能谱峰PMVG（即第二外区峰），其所诱导的流动结构展向波长λ+z与MVG 的展向周期（相邻MVG 的中心距离，≈350，如图7（b）中竖直虚线所示）基本一致，其法向高度对应于图6显示的〈uu〉+、〈vv〉+、〈ww〉+和〈uv〉+在对数区的峰值位置。此外，近尾迹区和分量的外区峰P2有所增强，并能够持续到远尾迹区。从图7（c）可以看到，由MVG 所导致的能谱峰PMVG在远尾迹区各分量能谱中基本消失，且、和的内区峰相较于光滑壁面工况均有降低、能谱内区峰值区域有所收缩，表明MVG 对近壁小尺度结构的抑制及对内外区不同尺度结构间的传能与调制作用降低，导致图6 中近壁处雷诺应力低于光滑壁面水平。

图7 光滑壁面及MVG0 下游近、远尾迹区展向预乘能谱Fig.7 Spanwise pre-multiplied energy spectra of smooth-wall and MVG0 arrays

4 POD 模态分解

对SPIV 测得的光滑壁面及MVG 阵列下游速度场用POD 进行模态分解，以厘清MVG所引入的流动结构在湍动能层面的表征。使用Snapshot POD 算法对各工况的1 000 帧三分量平面速度场进行POD 分解[41-42]。图8 给出了光滑壁面及MVG0阵列工况下各阶模态的能量占比及能量积累曲线，其中i 为模态阶数。图9 给出了光滑壁面及MVG0工况下近尾迹区（x/h0=10）和远尾迹区（x/h0=30）的部分典型POD 模态的空间基Ψi（y,z），云图为流向分量Ψiu，黑色箭头表示法向分量Ψiv、与展向分量Ψiw的合成矢量。

图8 光滑壁面及MVG0 工况近、远尾迹区POD 分解所得各阶模态的能量占比及能量积累曲线Fig.8 Energy ratio of each rank and the cumulative energy of the POD result of smooth-wall case and near- or far-wake regions of MVG0 case

图9 光滑壁面及基准MVG0 工况近、远尾迹区流场第1、5、10、20 阶POD 模态的空间基（y,z）Fig.9 Rank 1,5,10 and 20 mode Ψi(y,z) of POD decomposition results of smooth-wall case and near- or far-wake regions of MVG0 case

从图9 可以看到，光滑壁面工况下，第1 阶POD模态Ψi=1以贯通内外区的超大尺度结构为典型特性，高阶POD 模态Ψi=5,10,20则表现出内外区分离的双层结构特征。MVG0阵列工况下，近尾迹区第1 阶POD 模态中出现了类似图5 近壁处半圆形时均流向速度亏损区的“印记”结构，从壁面延伸到对数区中部。值得注意的是，该“印记”一直持续至第20 阶模态，说明MVG 所诱导的流动结构（时均流向涡对）贡献了整个SPIV 测量视野范围内的大部分湍动能（即该流动结构与湍流边界层内的大、超大尺度结构对湍动能的贡献相当），并由此造成图8 中MVG0阵列工况下低阶POD 能量贡献率明显高于光滑壁面工况。这与Bai 等[43]在湍流边界层内利用斜置沟槽产生大尺度流向环流时所观察到的各阶POD 模态能量分布情况类似。在MVG0阵列下游的远尾迹区，其所诱导的流动结构不再明显呈现于低阶POD 模态中，结合〈uu〉+、〈vv〉+及〈ww〉+仍然在对数区存在峰值的观察，可以推测，MVG 所诱导的时均流向涡对结构在远尾迹区难以维持规则的完整外形，而是衰减为沿展向蜿蜒的涡量层并继续影响速度统计量的分布。对比图9 中的第20 阶模态可以发现，MVG0诱导的流动结构影响了近壁处含能结构分布，间接表征了MVG 对下游结构的控制作用。MVGs和MVGh阵列工况下的低阶POD 模态在近、远尾迹区的空间分布特性与MVG0阵列工况下类似，不再另行给出。

5 MVG 阵列的减摩阻特性

本节讨论MVG 阵列下游的摩阻沿程变化，主要评估各型MVG 阵列的减摩阻效果，并分析间距、高度对MVG 减摩阻特性的影响。减摩阻率RD（Drag Reduction）根据下式计算：

式中：τ光滑与τMVG分别为免标定双层热膜摩阻传感器所测的光滑壁面和MVG 工况下同一流向位置的当地时均摩阻。

当U∞=14 m/s 时，3 种MVG 阵列下游x/h0=10～80 范围内减摩阻率RD的沿程变化如图10 所示，其中误差条表示多次重复实验的不确定度。可以看到，在x/h0=10 处，3 种MVG 工况下的当地时均摩阻相较于光滑壁面工况均略有增大（RD<0），增阻量在5%以内。增阻的主要原因是MVG 产生的时均流向涡对在该处距壁面较近，其诱导产生的上洗、下洗运动增强了近壁区各流层的动量掺混。当x/h0≥ 20 后，3 种MVG 均呈现出明显的减摩阻效果，增大MVG 阵列高度能显著提升减摩阻率。已有研究表明[21]，高度越高的MVG 阵列能够诱导时均流向涡对产生更大的环量，因而对近壁处湍流结构的抑制和规则化作用也更强，最终表现为减摩阻率的提升，这与Cheng 等[4]以等离子体涡流发生器产生时均流向涡对减摩阻的研究结果相符（Cheng 等通过增大电极电压提升时均流向涡对的强度，实现了减摩阻率的提升）。大间距MVG 阵列所诱导的时均流向涡对展向排列稀疏（图5（a）），对边界层内流动结构的不规则运动影响作用有限，在摩阻测量中则体现为大间距MVGs的减摩阻率低于MVG0阵列。但更高更密的MVG 具有更大的型阻，在实际应用时需考虑MVG 所引起的附加阻力。高度更高、展向间距更小的MVG 阵列的最优减阻位置更靠近上游，在峰值减阻率后，各型MVG 阵列的减摩阻率均沿程逐渐衰减，但在x/h0=80 处MVGh及MVG0仍然具有超过5%的减摩阻率，说明MVG 阵列的减摩阻效果能够至少持续至下游80 倍自身特征高度处。

图10 各型MVG 阵列下游减摩阻率的沿程变化Fig.10 Drag reduction in different streamwise stations behind each MVG arrays

6 结 论

在U∞=14 m/s、Reτ=3 453 的湍流边界层中，利用SPIV 及免标定双层热膜摩阻传感器对3 种MVG 阵列下游速度场统计特性及减摩阻的效果开展实验测量，结论如下：

1）MVG 在近尾迹区诱导产生时均流向涡对和时均流向速度亏损区，并在MVG 特征高度附近流层产生低平均动量及高湍流脉动区。随着时均流向涡对的沿程抬升，MVG 对湍流统计量的峰值扰动也在下游抬升至更高的流层。此外，MVG 引起近壁区雷诺切应力的降低，说明了MVG 的减摩阻能力。

2）展向预乘能谱表明，MVG 导致近尾迹区出现了第二外区峰值，其所诱导流动结构的展向波长与相邻MVG 的中心距离一致，并显著影响了近尾迹区湍流边界层外区峰的强度与形貌。POD 模态分析表明，MVG 所诱导的流动结构与湍流边界层内的大、超大尺度结构对湍动能的贡献相当，并影响了近尾迹区近壁含能结构的分布。

3）实验证实MVG 阵列具有减摩阻效果，减摩阻有效范围至少可持续至下游80 倍自身特征高度处。

后续应在MVG 减摩阻机理与参数影响规律等方面开展进一步研究，并在减阻应用中综合考虑MVG 的减摩阻效果及其额外型阻。

致谢：北京航空航天大学龙彦光博士、朱熠辰博士生在2D-PIV 的基础上开发并完善了本实验所用的SPIV 基础程序。加拿大新不伦瑞克大学高南教授、大连航华科技有限公司刘玄鹤博士研发了本实验所用的免标定双层热膜摩阻测量设备，天津大学王昊博士生指导了热膜传感器的制备。北京航空航天大学张清福、于悦、李拓、徐颢文等多名研究生协助开展了实验测量工作。在此，对上述个人与机构提供的帮助表示感谢。