流向凹曲率壁面湍流边界层的TRPIV实验研究

王轩，范子椰，陈乐天，唐湛棋, 2，姜楠, 2, *

1.天津大学 机械工程学院，天津 300354 2.天津市现代工程力学重点实验室，天津 300354

湍流边界层是工程中常见的流动现象，同时也是湍流领域研究的热点问题之一。在实际流动中通常会遇到在曲面上发展的湍流边界层现象[1]，而曲面上的流动往往也会受到由曲率引起的压力梯度影响[2-3]。一般来说，带有凹曲率和凸曲率的壁面会分别增加和减少湍流活动[4]，同时由曲率引起的压力梯度也会影响湍流活动[5]。针对曲率和压力梯度影响的研究也一直受到船舶、航天等领域的重视。

迄今为止，研究者们主要研究了曲率和压力梯度分别对湍流边界层的影响，其中针对压力梯度的研究一般采用在光滑平板上施加压力梯度的方法。Spalart[6]通过数值模拟发现顺压梯度在达到一定的强度后会使湍流边界层发生层流化。Aubertine等[7]对逆压梯度的影响进行了激光多普勒测速（LDV）实验，发现逆压梯度会扩大尾迹区域，使得对数律区向下收缩，雷诺应力剖面与平板类似，但外层的应力比平板更强。Joshi等[8]对顺压梯度的流场进行了粒子图像测速（PIV）实验，发现近壁区的涡结构受到了限制，而加速流动也使得涡结构的倾角减小。最近，Volino[9]对不同压力梯度下的湍流边界层进行了激光多普勒测速实验，认为顺压梯度抑制湍流，且使相干结构的流向尺度更大，而在逆压梯度中则会产生相反的现象。

曲率对湍流边界层单独影响的相关研究大多使用带有曲率的管道来尽可能地抑制由于过流截面变化产生的压力梯度。目前针对曲壁面的实验研究方式一般为单点速度测量，如热线等。Meroney等[10]使用毕托管对流向曲壁面进行了测量，发现凹曲面上的流速比平板更高，而凸曲率对湍流边界层的影响则完全相反。So等[11]在针对凹壁面的热线实验中发现凹曲率会增强湍流强度。Arolla等[12]通过数值模拟研究了凹曲率和由凹曲率引起的压力梯度对湍流边界层的影响，发现在曲率起始处会出现逆压梯度，使壁面摩擦减小，而凹曲率对流场中的法向分量影响最大。Matsubara等[13]对弯曲槽道进行了模拟研究，发现曲壁面会显著增加大尺度结构的强度。

曲壁面上的湍流边界层会不可避免地受到由曲率引起的压力梯度作用[12]，因此曲壁面湍流边界层相比平板湍流边界层是一个更加复杂的流动。目前针对曲率和压力梯度综合影响下湍流边界层的研究仍局限于最基本的统计量，而对于内部流动及其结构的实验研究较少。本文利用高时间分辨率粒子图像测速技术（Time-Resolved Particle Image Velocimetry，TRPIV），获取凹壁面湍流边界层上不同流向位置的二维速度矢量场，利用大小视场获取法向高度上完整的平均速度剖面以及雷诺应力剖面，选取不同的法向高度并应用可以体现涡旋转方向的Λci准则进行条件相位平均，提取流场中的顺向涡（即沿顺时针旋转的展向涡），然后利用空间两点相关方法研究不同法向高度的相干结构形态，分析凹壁面上湍流边界层内部流动的变化。

实验在天津大学流体力学实验室重力溢流式低湍流度循环水洞中进行。水洞实验段尺寸为4.1m×0.6m×0.7m（长×宽×高），流速最高可达 0.5m/s。实验模型为带有流向凹曲率的曲面模型，总长度为2.75m，其前半部分为平板，尺寸为1.20m×0.59m×0.01m（长×宽×高），后半部分为曲板，尺寸为1.55m×0.59m×（0.01～0.18）m（长×宽×高）。在实际工程中，湍流边界层流动常存在不同程度的曲率作用情况，而发生小曲率作用的情况较为普遍，加之实验段尺寸也存在一定限制，最终选定曲面的型线为函数y=0.01ex。由于该模型的曲率半径较大，相机视野中拍摄区域的壁面接近水平，因此在以互相关算法进行计算时可以选择矩形计算域。为获得充分发展的湍流边界层，在平板前缘下游 5 0mm处增加了直径d=2mm的绊线。实验装置如图1所示，图中x、y和z分别代表流向、法向和展向。

图1 实验装置示意图Fig.1 Experimental facility diagram

为了获得较为完整的流场速度分布以及近壁面处的精细流动，采用大视场和小视场相结合的方式进行实验：共使用3台高速相机进行拍摄，其中2台相机分别布置在实验段的两侧，用于大小视场的拍摄，第3台高速相机用于拍摄实验段前端的自由来流速度。实验中，在两个不同位置使用不同激光进行照射。

由于实验模型带有沿流向的凹曲率，使得湍流边界层在不同流向位置受到的影响程度各不相同，因此在模型曲面部分的上游和下游各选择一个位置进行测量（图2），两个位置与模型前端的距离定义为L。为便于后期计算，调整高速相机视场与来流呈一定夹角，使得测量时观察到的壁面尽可能平行于相机视场。图中θ为相机视场与来流方向的夹角，在上游和下游测量位置分别为θ1=5.5°和θ2=9.0°。

图2 实验测量示意图Fig.2 Schematic diagram of experimental measurements

由于模型的拍摄表面为曲面，且在拍摄范围内其曲率变化不明显，在分析曲面流场时选取矩形计算域进行计算。为便于分析，定义x*为在不同位置选取的流场计算区域内的流向，y*为法向，U*和V*、u*和v*分别为流向和法向的平均速度和脉动速度。

实验使用的示踪粒子为直径20 µm的空心玻璃微珠，密度为 1.04kg/m3。3台高速相机为丹麦Dantec公司的SpeedSense 9072 CMOS相机，分辨率为1280像素×800像素。大视场相机对应的物理空间尺寸为 8 0.39mm×50.25mm（流向×法向），小视场相机对应的物理空间尺寸为28.26mm×17.66mm（流向×法向），小视场包含于大视场中。采样频率f=1000Hz，激光片光的厚度约为1 mm。两个位置的测量分两次进行，因此入口处的自由来流速度存在一定波动，在测量上游和下游时分别为0.316m/s和0.328m/s。

实 验 水 温T=20°， 水 的 密 度ρw=998.2kg/m3，水的运动黏度ν=1.01×10-6m2/s。实验共采集了3组流场数据，每组实验获取了两个不同位置大小视场各8 216张流场照片。大视场使用32像素×32像素的查询窗口，小视场由于粒子浓度较低，使用 4 8像素×48像素的查询窗口；
大小视场查询窗口的重叠率均为 7 5%。以互相关算法进行计算，小视场相邻矢量的间距为 0.26mm，大视场相邻矢量的间距为 0.5mm。

目前，尚无法在实际实验中将流动中压力梯度和流向曲率的影响进行单独分析。此前相关研究常采用零压力梯度的光滑平板流动作为基准进行对比研究。本文亦开展了光滑平板实验。实验中同样使用2台高速相机进行大小视场的拍摄。来流速度为0.267m/s；

采样频率f＇=800Hz；
大视场相机对应的物理空间尺寸为 1 26mm×77mm（流向×法向），小视场相机对应的物理空间尺寸为 2 7mm×17mm（流向×法向），小视场位于大视场中；
其余参数与曲面流场相同。计算矢量场时同样选取矩形计算域，因此使用与前文相同的坐标系x*和y*定义。

流向凹曲率壁面的湍流边界层流动不仅会受到凹曲率的作用，还会受到由凹曲率引起的流向顺压梯度的影响，使得湍流边界层存在沿流向逐渐增大的加速度。由于模型型线曲率半径较大，在初期测试中发现该模型壁面上的湍流边界层受到的凹曲率作用较弱，而受到的顺压梯度作用相对较大。对于带有流向顺压梯度的流动，需要了解流向顺压梯度引起的加速度强弱，研究者普遍使用加速度参数K[14-17]来表征这一强弱程度：

式中，U0表示当前测量位置的自由来流速度。当K＞3.0×10-6时，湍流边界层会发生层流化现象[6]。此外，曲面会使流体在不同流向位置流过的截面大小非线性变化，从而使得流向加速度在流向上非均匀；
但本实验测量区域尺寸不足以使K产生明显变化，因此选取平均加速度参数作为测量位置的加速度参数。通过对K的计算发现，在上游和下游测量位置均满足K＜2.0×10-6，流向顺压梯度的作用较弱，因此可以认为在该曲面上的湍流边界层未发生层流化现象。为比较其与平板的湍流边界层流动情况，引入光滑平板模型的实验数据进行对比，平板数据同样被分为大视场和小视场数据。本实验顺压梯度和流向曲率对湍流边界层的作用较弱，传统的对数律方程仍然具有较强的鲁棒性[6,18-19]，因此本文利用平均速度剖面的对数律区与对数律公式U+=（1/κ）lny++B进 行拟合来获取壁面摩擦速度uτ，其中y+=yuτ/ν、U+=U/uτ分别为内尺度无量纲化的法向高度和平均速度，κ为卡门常数，B为积分常数。由于研究对象为曲面湍流边界层，其对数律方程可能产生变化，因此将卡门常数κ与积分常数B均设为变量进行拟合，此外，还应用Spalding曲线[20]进行了拟合验证。结果发现：κ和B均未发生明显变化，Spalding曲线拟合结果也与对数律方程拟合结果相似，可以认为拟合得到的壁面摩擦速度未发生较明显变化。表1列出了实验基本湍流参数，其中，U∞为入口处来流速度，uτ，large、uτ，small分别为大小视场的壁面摩擦速度（本实验将小视场的uτ，small用于后续数据处理），边界层厚度δ99按照0.99倍当前测量位置的平均来流速度确定，剪切雷诺数Reτ=uτδ/ν表示内尺度雷诺数（uτ随着边界层向下游的发展逐渐增大，且下游增长幅度大于上游）。

表1 边界层参数Table 1 Boundary layer parameters

2.1 基本统计量

图3为在曲面模型上两个不同流向位置得到的内尺度无量纲平均速度剖面，横轴为两个测量位置与壁面垂直的内尺度无量纲法向高度，定义U*+=U*/uτ，small，y*+=y*uτ，small/ν。为获取高法向空间分辨率及近壁处平均速度剖面，对小视场利用单行互相关算法（Single-Row Cross Correlation，SRCC）思想进行流场高法向空间分辨率速度矢量计算[21-23]。SRCC算法的思想是通过牺牲流向分辨率来提升法向分辨率，以获得更高法向空间分辨率的速度矢量。因此，基于该算法的思想将查询窗口大小调整为256像素×8像素（流向×法向），重叠率为 5 0%。经实际计算，速度矢量的法向间距可达 0.09mm，同时可测量到的流场矢量最高可以达到黏性底层。由于实验模型为曲面模型，大视场不同流向位置的平均速度剖面存在较为明显的变化，为了统一无量纲化后的结果，大小视场均使用小视场的壁面摩擦速度uτ，small进行无量纲化，同时选取大视场与小视场拍摄位置重合的流向位置进行绘制。

图3 平均速度剖面Fig.3 Mean velocity profile

从图3可以看出，凹壁面上游速度剖面的对数律区变化不明显，但下游速度剖面的对数律区发生了明显偏离，速度剖面在对数律区向上隆起且对数律的斜率小于传统对数律方程。在缓冲层及以下区域，凹壁面与光滑平板情况下的速度分布无明显差异，这可能是由于凹曲率和顺压梯度对边界层内流向速度的影响相互抵消所导致[10，24]。在对数律区以上，两个测量位置的流动均受到抑制[9]，说明顺压梯度在对数律区以上造成的影响更强。

图4为不同位置的内尺度无量纲化雷诺应力随y*+变化的曲线，雷诺应力各分量均使用壁面摩擦速度实现内尺度无量纲化。

图4 雷诺应力剖面Fig.4 Reynolds stress profile

从图中可以看到：凹壁面湍流边界层的湍流强度显著减小，且随着边界层向下游发展，流向雷诺正应力在近壁处的峰值逐渐沿法向向上移动，达到y*+=20左右；
在对数律区以上时，与光滑平板情况之间差距增大；
法向湍流强度和雷诺应力在近壁附近无明显变化，而在对数律区及以上区域明显下降。这说明流向顺压梯度对湍流有显著的削弱作用，且随着法向高度增大，削弱作用先增大后减小。

2.2 涡结构识别

从2.1节对凹壁面流场基本统计量的分析中可以发现，湍流边界层在受到凹曲率和流向顺压梯度的组合影响后，湍流强度整体呈现下降趋势，这一影响势必也会对边界层内的流动结构产生影响。为了研究凹壁面流场内流动结构的变化，本文利用λci准则进行流场内涡的识别[25-27]。λci为速度梯度张量的共轭复特征值，该值可以表示流场中某一矢量点处的旋涡强度。应用该准则可以有效避免将剪切运动识别为涡旋运动。若将旋涡强度与展向涡涡量符号sgn（ωz）结合，即可识别出流场内涡结构的旋转方向：

式中，Λci表 示获得涡量符号后的λci值。

图5给出了凹壁面湍流边界层不同流向位置的Λci识 别结果，x*与y*均使用内尺度进行无量纲化。当Λci为负时，检测到顺向涡，为正时，则检测到逆向涡。为减少环境噪声对识别结果的影响，使用阈值条件|Λci|≥T×Λci，rms（y）（T为检测涡结构边界的阈值，本文取0.5）来确定涡结构的边界，同时考虑到流场的空间分辨率影响，将仅有单个网格节点大小的涡结构过滤，以减弱背景噪声影响[28]。从图中可以看出，流场中Λci较大的区域与速度矢量中直观显示出的涡结构吻合程度较高，流场内存在顺向涡和逆向涡的涡对，以及由发卡涡组成的发卡涡包结构。

图5 Λci识别结果Fig.5 Λci recognition results

为进一步了解凹壁面流场中涡结构的强度沿法向分布情况，对大视场中的λci进行了时间和沿流向的平均，并绘制了大视场λci沿内尺度无量纲化法向高度y*+的变化曲线，如图6所示。从图中可以发现，凹壁面和光滑平板上的λci均显示出随着法向高度的变化先增大后减小的趋势，而凹壁面上不同流向位置的旋涡强度在各法向位置上均大于平板情况，且在对数律区附近相差更大。这说明虽然凹壁面边界层内的湍流强度被削弱，但由于凹壁面对湍流边界层的影响，流场内不同位置的速度梯度增强，进而使得旋涡强度增大，同时也使涡的分布更加集中于对数律区附近。

图6 大视场 λci随 y＊+的变化曲线Fig.6 Variation of λci in boundary layer

2.3 条件相位平均

顺时针旋转的展向涡通常被认为与发卡涡的涡头密切相关，而发卡涡在向下游迁移的过程中会对周围流场施加诱导作用，产生猝发事件[29]，对湍流内部的能量产生和输运起到了重要作用，发卡涡的形成与分布也与湍流脉动密切相关[30]。

条件相位平均是研究流场中涡结构的典型方法。为研究凹壁面上展向发卡涡头的变化，2.3小节使用该方法提取了小视场内顺向涡（即顺时针旋转的展向涡）的统计特征。展向涡的识别条件使用2.2小节所述的Λci准则[31]，可表示为：

当检测位置达到检测条件，提取该时刻检测位置附近的流场。本文在小视场内选取了2个法向参考高度进行条件采样。

图7 小视场yr＊+e f=30附近顺向涡流向脉动速度的条件平均等值线图Fig.7 Conditional phase-averaged contour of the fluctuation velocity in the direction of the downward vortex near y＊re+ f=30 in a small field

图8 小视场y＊re+ f=80附近顺向涡流向脉动速度的条件平均等值线图Fig.8 Conditional phase-averaged contour of the fluctuation velocity in the direction of the downward vortex near yr＊e+ f=80 in a small field

2.4 两点相关分析

从之前分析可知，流场内的顺向涡在凹曲率和流向顺压梯度的组合影响下，产生了较大变化。众所周知，顺向涡与发卡涡头有着密切联系，而发卡涡是壁湍流中普遍存在的相干结构[32]，因此，在凹曲率和流向顺压梯度的影响下，湍流边界层内相干结构的运动和发展也会产生一定变化。

针对相干结构在不同法向高度的流向分布特征，本文使用两点流向脉动速度相关法来体现凹壁面上湍流边界层内部相干结构在受到凹曲率和流向顺压梯度组合影响时产生的变化[13,33-34]。两点流向脉动速度相关计算公式如下：

式中，yref表示计算选取的参考法向高度， Δx表示两点之间的流向空间延迟，σu表示流向脉动速度的均方根。

为研究近壁面附近相干结构的空间形态，图9绘制了小视场两点相关结果，所选取的参考高度与近壁处条件平均所选取的高度一致，x*与y*均使用内尺度进行无量纲化。从图中可以看到，相较于光滑平板情况，凹壁面近壁面处相干结构的空间尺度呈现先缩小后增大的趋势。在两点相关结果中，相干结构呈现出椭圆外形。筛选出某一等值线上的所有点，并对这些点进行椭圆拟合，即可获得椭圆方程，基于此方程可以获取相干结构与壁面之间的夹角。本文选取=0.4等值线（红色虚线）上的点作为参考点，计算出的相干结构与壁面夹角已在图9中指出。可以发现：位于缓冲层附近的相干结构与壁面的夹角明显小于光滑平板情况；
随着边界层向下游发展，相干结构逐渐向壁面倾斜，其沿流向的尺度也被逐渐拉长，这与仅有流向顺压梯度影响时的趋势类似[10]。但由于流向凹曲率的影响，相干结构的空间尺度在上下游增长不均匀，下游的空间尺度增长幅度更大。

图9 小视场流向两点相关系数等值线图Fig.9 The contour of two points correlation coefficient in the streamwise direction in the small field

小视场可以呈现出流场近壁处细致的空间形态，而大视场则可以较为完整地表征相干结构的空间尺度大小。图10为不同工况下大视场两点相关系数的结果，选取的参考法向高度在附近，x*与y*均使用内尺度进行无量纲化。在该高度附近，湍流边界层受到凹曲率和流向顺压梯度的影响，对数律区相干结构的尺度明显增大，且随着湍流边界层向下游的发展，相干结构的空间尺度逐渐增大，在下游位置超过了1倍边界层厚度（内尺度无量纲化后的边界层厚度δ+=δuτ/ν=753.26）。

图10 大视场流向两点相关系数等值线图Fig.10 The contour of two points correlation coefficient in the streamwise direction in the large field

通过双相机大小视场的TRPIV实验，对带有纵向凹曲率的壁面进行了测量并获取了瞬时速度场，经过数据处理和分析得出以下结论：

1）在流向凹曲率壁面上得到的平均速度剖面随着向下游发展逐渐偏离传统对数律，同时尾迹区被抑制。在相同法向高度的湍流强度均弱于平板情况。这说明流向凹曲率壁面上的湍流活动被明显抑制，且尾迹区的抑制作用相对更强。

2）相比于平板情况，凹壁面湍流边界层中涡的强度在不同法向高度均明显增强，且随着向下游的发展而逐渐增大。对数律区及其下方区域与尾迹区旋涡强度的差值增大，这表明凹壁面上湍流边界层内的涡强度显著增强，且更集中于近壁附近。

3）对凹壁面湍流边界层对数律区内的顺向涡进行条件相位平均，发现凹壁面会增强顺向涡下方负的流向脉动速度。缓冲层附近顺向涡上方正脉动峰值被削弱，但空间尺度有所增大，而对数律区上侧正脉动更加集中于顺向涡上方较小的区域内，且脉动强度明显增大。随着向下游的发展，缓冲层附近顺向涡上方正脉动的峰值和空间尺度均弱于平板情况，而在对数律区上侧则均强于平板情况。

4）两点相关结果表明，凹壁面上相干结构向外层迁移的趋势减弱，其空间尺度随着向下游的发展呈现出不断增大的趋势，其中上游的空间尺度增长被抑制，而在下游则相反。