Experimental study of TRPIV for turbulent boundary layer of longitudinal concave curvature wall
-
摘要: 带有曲率和压力梯度作用的湍流边界层广泛地存在于各类工程应用中。以带有流向凹曲率的壁面(简称凹壁面)模型为研究对象,利用双相机大小视场的TRPIV系统分别测量了上游和下游两个不同流向位置的瞬时速度场。以光滑平板模型的实验数据为基准,研究了凹壁面上湍流边界层的变化规律。通过对比光滑平板的平均速度剖面和雷诺应力剖面,发现凹壁面情况下的平均速度剖面逐渐偏离了传统对数律,且湍流强度比平板情况下更弱。以
$ {\varLambda _{ci}} $ 准则识别出的顺向涡为条件进行条件相位平均,发现凹壁面中涡上方正脉动的峰值与空间尺度的变化相反,而下方的负脉动比平板情况更强。进一步使用空间两点相关方法提取相干结构,同时利用椭圆拟合的方法计算近壁区相干结构与壁面的倾角,发现凹壁面上相干结构沿流向的空间尺度逐渐增大。结果表明:当湍流边界层受到凹曲率和顺压梯度的共同作用时,近壁附近与尾迹区内湍流强度的差异增大,同时边界层内相干结构的旋涡强度增加;随着向下游的发展,缓冲层附近顺向涡上方正脉动的空间尺度增大,但峰值降低,而在对数律区的上侧则产生了相反的现象;凹壁面湍流边界层内相干结构向更高法向高度迁移的趋势减弱,且相干结构的空间尺度在下游一侧的增长幅度更大。Abstract: Turbulent boundary layers with curvature and pressure gradient are widely used in various engineering applications. A wall model with concave curvature in streamwise was used as the research object, and the instantaneous velocity fields at two different streamwise positions upstream and downstream were measured by TRPIV experiment with dual camera size field of view. Using experimental data from a smooth flat plate as a benchmark, the variation pattern of the turbulent boundary layer on a concave wall is investigated. By comparing the mean velocity profile and the Reynolds stress profile with those of the smooth flat plate, it was found that the mean velocity profile in the concave wall case gradually deviated from the conventional logarithmic law and that the turbulence intensity was weaker than that in the flat plate case. Conditional averaging is carried out under the condition of the prograde vortex identified by$ {\varLambda _{ci}} $ criterion. It is found that the peak of the positive fluctuations above the vortices in the concave wall varies oppositely to the scale range, while the negative fluctuations below the vortices is stronger than that of flat plate. The coherent structures were further extracted using a spatial two-point correlation method and the inclination of the coherent structures in the near-wall region to the wall was calculated using an ellipse fitting method. It was found that the scale range of the coherent structure on the concave wall gradually increases along the streamwise. The results show that when the turbulent boundary layer is subjected to a combination of concave curvature and favorable pressure gradient, the difference in turbulence intensity near the near-wall region and within the wake region increases, along with an increase in the intensity of vortices in the coherent structures within the boundary layer. With development downstream, the scale range of positive fluctuations above the prograde vortex on the buffer layer increases, but the peak value decreases, while the opposite phenomenon occurs on the upper side of logarithmic region. The trend of coherent structure migration to higher normal height in the curved wall turbulent boundary layer is weakened, and scale range of the coherent structures grows mainly towards the downstream side. -
0 引言
激波干扰现象,因其产生的激波、剪切层、膨胀波、超声速射流以及伴随的旋涡和湍流等复杂流动结构,不仅是高超声速飞行器研究与发展过程中需要重点关注的空气动力学问题之一,也是困扰工程应用的一大难题[1, 2]。以斜激波与圆柱/圆球弓形激波之间的相互作用为典型模型,根据激波干扰位置和激波强度不同,Edney[3]总结出6种具有普适性的激波干扰类型(Ⅰ~Ⅵ类)。其中,第Ⅳ类激波干扰因其产生的超声速射流冲击壁面会带来严峻的热/力载荷而受到设计和研究人员的广泛关注[4-7]。在某些情况下,激波干扰产生的超声速射流高度弯曲,不会直接冲击壁面,这种激波干扰被称为第Ⅳa类[8]。
高超声速进气道流动存在多种形式的激波干扰现象,尤其以设计工况下前体激波入射进气道唇口[9]形成的激波干扰最为突出。针对前体激波与二元进气道唇口弓形激波的干扰问题,Wieting等[10]进行了系统性研究,获得了丰富的压力和热流数据。近年来,前缘高度后掠的内转式进气道[11-13]逐渐成为研究热点之一。然而,其普遍采用的V形溢流口(唇口)与二元进气道唇口具有明显的区别。最近,Xiao等[14]在不考虑前体激波的情况下,将内转式进气道溢流口简化为V形钝前缘,研究发现:在一定的几何参数下,V形钝前缘自身就会形成类似于Edney第Ⅳ类激波干扰的流动结构。显然,当前体激波入射V形钝前缘溢流口时,将产生更加复杂的波系结构,很可能加剧溢流口的气动热。
另一方面,根据来流条件的不同,进气道会出现亚额定、额定和超额定等3种工作状态[15]。针对二元进气道的研究表明[16-17]:在超额定状态下,前体激波入射在唇口内部,当前体激波与唇口激波发生Edney第Ⅱ类激波干扰时,形成的亚声速区可能导致进气道局部不起动。相对而言,内转式进气道唇口更加复杂,前体激波与其相互作用的研究尚未见系统报道,亟需针对性的机理揭示,以补充认识上的不足。
采用完整的内转式进气道对V形溢流口激波干扰特性进行研究,不仅会耗费大量的计算和实验资源,更会带来多种复杂因素的影响。抓住内转式进气道V形溢流口的主要流动特征,有侧重地进行简化,继而揭示其流动机理和关键影响因素,当属高效简洁的研究方式。因此,本文将内转式进气道V形溢流口简化为V形钝前缘平板。鉴于内转式进气道的前体激波为复杂的三维曲面激波,为便于准确描述和控制,本文暂不考虑入射激波的展向曲率,而是采用简化的斜激波,以突出V形钝前缘溢流口自身带来的主要特征。继而采用激波风洞实验和数值模拟相结合的方法,揭示入射激波与V形钝前缘溢流口相对位置变化引起的激波干扰形态演变,以期为内转式进气道溢流口设计提供有价值的参考。
1 实验和计算方法
1.1 实验设备和模型
V形钝前缘平板由较长的平直段和交叉位置倒圆部分组成,如图 1所示。交叉位置倒圆半径R=6.5mm,前缘钝化半径r=2mm,前缘扩张角β=48°。前缘平直段的长度L=31r,以使其产生的脱体激波充分发展。以V形前缘驻点位置为坐标原点,定义流向为x方向,平板展向为z方向,垂直于平板方向为y方向。
为了更加清楚地观察波系结构,对V形钝前缘平板实验模型进行了2.5倍放大。实验在中国科学技术大学的KDJB330反射型激波风洞[18-20]中进行。实验来流马赫数为6,静压为891Pa,静温为101K。纹影系统采用Z字形光路布置,如图 2所示。通过高速摄影记录流场纹影图像, 拍摄帧率为15000Hz,曝光时间为2μs。入射斜激波由气流偏转角为12°的斜劈产生,斜劈流向长度为200mm,展向宽度为260mm。为了获得丰富的流场信息,考虑到模型对光路的遮挡,分别从平行于平板(z方向,见图 2)和垂直于平板(y方向,见图 3)2个方向对流场进行拍摄。通过前后移动斜劈的方式,改变入射斜激波与V形钝前缘溢流口的相对位置。为了便于分析,斜激波入射位置采用不考虑激波干扰时无粘斜激波在平板上的入射点坐标。
1.2 计算方法
数值计算采用Fluent软件求解Navier-Stokes方程,并使用具有二阶精度的格式离散。对流通量采用AUSM格式,湍流模型采用k-ω SST模型,气体分子粘性采用Sutherland公式,气体定压比热使用CHEMKIN温度多项式进行计算。对于V形钝前缘平板流场的计算,由于来流和几何的对称性,为了节省计算资源,采用1/4计算域进行模拟。采用结构化六面体网格,计算域及局部网格如图 4和5所示,近壁面采用等比加密,壁面y+ < 1,总网格量约3000万。有斜激波入射时的计算域如图 6所示,此时六面体网格总数在5000万左右。为了与风洞实验进行对比,来流条件设置与实验一致。在数值模拟中监视各方程残差,监测出口位置的流量及质量平均马赫数、驻点位置的壁面压力和热流量。待残差收敛并且各监测参数稳定,认为流场收敛。
由于三维流动的复杂性,数值模拟获得的流场结构难以直接与具有光路积分效应的实验纹影进行对比。为此,本文按照纹影法的原理对数值模拟流场进行处理。如为了获得z方向穿过全流场的数值纹影,将流场沿z方向划分为k个等间距的区域,假设每个小区域内的流场参数沿z方向均匀分布,则经过流场的光线偏转角εx由式(1)[21]得到:
(1) 式中,n为折射率,ξ1和ξ2分别为光线进入流场和离开流场位置的z坐标。由此获得的光强变化可以与z方向穿过流场的实验纹影进行直观对比。
2 结果与分析
V形钝前缘自身的三维绕流会导致复杂的激波干扰,因此在研究斜激波入射V形钝前缘平板激波干扰前,有必要对V形钝前缘基本波系结构进行研究。
2.1 V形钝前缘平板基本流场
从y方向拍摄的流场纹影与流场x-z对称面数值纹影的对比如图 7(a)所示,可以看到数值模拟较好地捕捉到了基本的波系结构。图 7(b)为上半部分波系结构的示意图,其中DS代表V形钝前缘斜直段产生的脱体激波,BS代表V形前缘交叉位置形成的弓形激波,脱体激波DS入射弓形激波BS,相交于三波点T,在三波点位置产生透射激波TS和滑移线Σ。脱体激波DS波后的超声速气流在V形前缘交叉倒圆位置被压缩,压缩波汇聚形成弯曲激波CS。透射激波TS与弯曲激波CS发生异侧激波干扰,形成马赫杆MS。激波BS以及马赫杆MS波后气流均为亚声速,而透射激波TS波后气流为超声速,超声速气流在亚声速区域之间不断地压缩和膨胀,交替出现激波和膨胀波,形成超声速射流结构。超声速射流沿壁面运动,在驻点前发生碰撞,并卷起形成反向旋转的涡对结构。涡对结构与弓形激波BS相互作用,使得弓形激波波面向上游凸起,并出现Kink。将脱体激波DS视为入射激波,则V形钝前缘激波干扰结构与Edney第Ⅴ类激波干扰类似。流场中存在射流失稳和剪切层脱涡等流动过程,数值模拟虽然不能捕捉这些非定常过程,但能够刻画基本波系结构,流动的非定常性没有显著影响本文对波系干扰类型的研究。
为了探讨流场波系的三维结构,从z方向也进行了纹影拍摄。使用本文1.2节所描述的三维流场后处理方法得到的数值纹影与实验纹影图像进行了对比,如图 8(a)所示。由于流场的三维性,在z方向纹影图像上可以看到3道激波,从左到右依次为弓形激波BS上凸起的激波面(Bulge)、弓形激波BS和弯曲激波CS。图 8(b)进一步定量比较了数值模拟与实验获得的激波BS及其凸起的激波面的位置。结果表明:激波BS的位置符合较好,而凸起的激波面(Bulge)略有差别。这是由于激波BS相对稳定,而凸起的激波面受涡对结构的影响(见图 7),具有一定的非定常性。整体而言,数值模拟与实验纹影吻合较好,可用于本文后续的激波干扰研究。
图 9比较了相同钝化半径下二维圆柱产生的弓形激波和V形钝前缘平板弓形激波。其中,Δ1和Δ2分别代表V形钝前缘和圆柱弓形激波的脱体距离,V形钝前缘弓形激波由x-y对称面的数值模拟给出,圆柱产生的弓形激波形状由Billig[22]公式给出,而弓形激波面上声速点位置的激波角β*由以下公式[23]给出:
(2) 式中,Ma∞为来流马赫数, γ为气体比热比。从图 9可以看出,由于自身的激波干扰,V形钝前缘交叉位置弓形激波的脱体距离较大,并且驻点前存在大范围的亚声速区。进一步采用马赫数为5.95的等值面给出V形钝前缘平板的三维激波结构(图 10),可以看出:前缘脱体激波DS和弓形激波BS均为三维曲面激波,并且BS波面有较大范围的凸起。这表明V形钝前缘弓形激波的近正激波区域较大,当斜激波与弓形激波干扰时,出现第Ⅲ和第Ⅳ类干扰的范围更大[24]。
2.2 斜激波入射V形钝前缘溢流口激波干扰流场
当入射斜激波距离V形溢流口上游较远时,对应于亚额定工况下内转式进气道前体激波与唇口激波相交于唇口外的情况。图 11(a)为斜激波无粘入射平板位置坐标为x/r=-10.5时的流场纹影图片,图 11(b)为波系示意图。斜激波IS从平板下侧入射弓形激波BS超声速区域,发生同侧激波相交,出现Edney第Ⅵ类激波干扰。激波相互作用导致弓形激波的激波角减小,入射激波干扰后形成透射激波TS和剪切层Σ。
向下游移动斜劈,当斜激波无粘入射位置为x/r=-3时,斜激波入射在V形前缘驻点前的亚声速区,此时对应于内转式进气道额定工况下前体激波入射唇口的情况。仅从z方向拍摄纹影时,由于模型遮挡,难以获得干扰点附近激波干扰形态。因此,需要结合y方向的纹影来判断干扰类型。从图 12(a)流场z方向的纹影可以看出, 入射斜激波IS分为了两道斜激波,分别为IS和IS′。其中,IS′是入射斜激波IS与V形钝前缘的脱体激波DS干扰产生的透射激波。IS′入射V形前缘平板壁面形成反射激波RS,而入射斜激波IS与弓形激波BS的干扰点被模型遮挡。从图 12(b)流场y方向的纹影可以看到,在驻点前存在两道弓形激波BS和BS′。结合图 12(a)和(b)两个方向的纹影,推测激波干扰类型可能为Edney第Ⅳ类激波干扰,其流场结构如图 13所示。此时,由于斜激波IS入射,原有的弓形激波BS转变为通过透射激波TS连接的两段弓形激波BS和BS′,并伴随着超声速射流的产生,该流场结构很复杂,将在2.3节进一步分析。因此,在内转式进气道额定工况下,前体激波入射溢流口驻点前亚声速区域,可能带来复杂的激波干扰流场,对进气道流动带来负面影响。
进一步向下游移动斜劈,使得斜激波无粘入射位置为x/r=0.2,此时对应于内转式进气道超额定状态下前体激波入射溢流口内部的情况。流场纹影如图 14(a)所示,图 14(b)为其流场示意图。入射激波IS与弓形激波BS干扰形成的剪切层Σ并没有撞击壁面,而是形成一种有别于经典的Edney第Ⅲ类激波干扰的结构。其原因在于:V形钝前缘脱体激波DS与弓形激波BS干扰形成的脱体距离和亚声速区域较大,使得发生第Ⅲ类激波干扰的区域增大[23]。
当斜激波无粘入射位置为x/r=3时,斜激波与弓形激波BS的超声速区域干扰,流场纹影如图 15(a)所示,流场如图 15(b)所示。入射斜激波IS与弓形激波BS相交,形成长度较短的马赫杆MS、透射激波TS以及剪切层Σ1和Σ2。这种由Edney第Ⅱ类激波干扰形成的亚声速区,可能会带来进气道局部不起动[16-17],应避免出现。
2.3 斜激波/弓形激波第Ⅳ类激波干扰流场分析
当斜激波入射V形钝前缘时,三维波系干扰的情况如图 16所示(这里采用斜激波在上)。根据激波位置和特点,可能会出现3种干扰情况:第一种是入射斜激波IS与脱体激波DS相交,第二种是入射斜激波IS与弓形激波BS相交,第三种是入射斜激波与激波BS和DS的交线相交。鉴于三维激波干扰的复杂性,为了深入理解斜激波入射V形前缘时的激波干扰流场,本节以斜激波无粘入射壁面位置坐标为x/r=-5时作为典型工况,结合数值模拟进行分析。与入射位置为x/r=-3的工况(图 12)相比,此时激波干扰点位置更接近BS的正激波部分。
图 17(a)采用流向截面和展向对称面(z/r=0)展示了斜激波IS与V形钝前缘波系干扰三维结构的演化情况。其中,流向截面的马赫数云图仅显示了一半流场,而展向对称面通过马赫数等值线显示。在x/r=-5截面上(局部放大见图 17(b)),入射斜激波IS与脱体激波DS的较弱部分发生异侧相交,形成Edney第Ⅰ类干扰,产生的透射激波TS相对于激波IS向DS激波面偏折。在x/r=-2.5截面上(局部放大见图 17(c)),入射斜激波IS与脱体激波DS的较强部分发生异侧相交,形成Edney第Ⅱ类激波干扰,并产生马赫杆MS,此时的透射激波TS1入射钝前缘壁面。图 17沿流向发展的波系结构表明:由于入射斜激波IS与前缘脱体激波DS相交处DS激波强度的变化,激波干扰类型由Edney第Ⅰ类转变为第Ⅱ类。而在V形钝前缘溢流口处,如图 17(a)中的蓝色线框所示,入射激波IS与弓形激波BS发生干扰并产生超声速射流,这一区域的流场结构将在图 18中进行分析说明。
在斜激波入射条件下,多股气流共同汇集在V形钝前缘交叉位置,形成更加复杂的三维激波干扰流场。图 18给出了V形钝前缘交叉位置z/r=0和y/r=-0.4两个正交截面上的流场结构。从z/r=0截面可以看出,入射斜激波与弓形激波BS干扰后,形成通过透射激波TS4连接的两段弓形激波BS和BS′以及超声速射流。该超声速射流并未直接入射壁面,而是发生高度弯曲沿壁面方向溢流,此干扰类型是第Ⅳ a类干扰。而y/r=-0.4截面展示了斜激波IS同时与弓形激波BS和脱体激波DS发生多波系干扰的情况,由于IS与DS干扰已经形成了马赫杆MS(见图 17(c)),MS再与弓形激波BS干扰形成透射激波TS3以及超声速射流区域。可见,入射激波IS、弓形激波BS与脱体激波DS三者之间的干扰互相耦合。
将数值模拟流场沿y方向进行三维纹影处理,如图 19(a)所示,与图 19(b)给出的实验纹影图像符合较好。从三维纹影结果中可以看到,V形钝前缘溢流口弓形激波BS被分为两道激波BS和BS′。实验纹影图像与数值计算结果相互印证了斜激波与V形钝前缘弓形激波发生了第Ⅳa类激波干扰。
3 结论
通过简化构型对内转式进气道V形钝前缘溢流口激波干扰特性进行了研究,得到以下结论:
(1) V形钝前缘自身形成的激波干扰,导致其驻点前弓形激波的脱体距离较大,波后存在大范围的亚声速区,因此在估算内转式进气道溢流口激波脱体距离时,需慎用圆柱脱体激波距离计算公式。
(2) 斜激波入射在V形钝前缘驻点前弓形激波接近正激波的位置时,发生Edney第Ⅳa类激波干扰。此时,入射激波与弓形激波BS的干扰以及弓形激波BS与脱体激波DS的干扰互相耦合,形成多处超声速射流区域。因此,在内转式进气道额定工况下,前体激波入射V形溢流口驻点前亚声速区域可能带来复杂的流场结构,在实际应用中应当予以关注。
(3) 当斜激波入射在V形钝前缘弓形激波亚声速区的声速点附近时,出现一种有别于Edney第Ⅲ类激波干扰的波系结构,此时干扰产生的剪切层未入射壁面。
(4) 入射斜激波与V形钝前缘脱体激波DS、弓形激波BS的超声速区域干扰时的形态与Edney第Ⅰ、Ⅱ、Ⅵ类激波干扰类似,其中第Ⅱ类干扰形成的马赫杆后局部亚声速区域,可能造成内转式进气道局部不起动。
鉴于V形溢流口三维流动的复杂性,在激波干扰引起的气动热以及三维前体激波入射方面,还有待进一步的细致研究。
-
表 1 边界层参数
Table 1 Boundary layer parameters
上游(L=2.0 m) 下游(L=2.3 m) 平板实验 U∞/(m·s−1) 0.316 0.328 0.267 uτ,large/(m·s−1) 0.0175~0.0180 0.0196~0.0203 0.0116 uτ,small/(m·s−1) 0.0179 0.0202 0.0116 δ99/mm 36.85 37.29 60.16 Reτ 659.6 753.3 697.9 K 9.03×10−7 1.68×10−6 — -
[1] YOU J,BUCHTA D A,ZAKI T A. Concave-wall turbulent boundary layers without and with free-stream turbulence[J]. Journal of Fluid Mechanics,2021,915:A4. doi: 10.1017/jfm.2021.12
[2] HOLLOWAY A G L,ROACH D C,AKBARY H. Combined effects of favourable pressure gradient and streamline curvature on uniformly sheared turbulence[J]. Journal of Fluid Mechanics,2005,526:303-336. doi: 10.1017/s0022112004003088
[3] LOPES A S,PIOMELLI U,PALMA J. Large-eddy simula-tion of the flow in an S-duct[J]. Journal of Turbulence,2006,7(11):1-24. doi: 10.1080/14685240500331900
[4] PATEL V C,SOTIROPOULOS F. Longitudinal curvature effects in turbulent boundary layers[J]. Progress in Aero-space Sciences,1997,33(1-2):1-70. doi: 10.1016/S0376-0421(96)00001-2
[5] HARUN Z,MONTY J P,MATHIS R,et al. Pressure gradient effects on the large-scale structure of turbulent boundary layers[J]. Journal of Fluid Mechanics,2013,715:477-498. doi: 10.1017/jfm.2012.531
[6] SPALART P R. Numerical study of sink-flow boundary layers[J]. Journal of Fluid Mechanics,1986,172(1):307-328. doi: 10.1017/s0022112086001751
[7] AUBERTINE C D,EATON J K. Turbulence development in a non-equilibrium turbulent boundary layer with mild adverse pressure gradient[J]. Journal of Fluid Mechanics,2005,532:345-364. doi: 10.1017/s0022112005004143
[8] JOSHI P, LIU X F, KATZ J. Turbulence in accelerating boundary layers[C]//Proceedings of ASME/JSME/KSME Joint Fluids Engineering Conference. 2012. doi: 10.1115/AJK2011-25010
[9] VOLINO R J. Non-equilibrium development in turbulent boundary layers with changing pressure gradients[J]. Journal of Fluid Mechanics,2020,897:A2. doi: 10.1017/jfm.2020.319
[10] MERONEY R N,BRADSHAW P. Turbulent boundary-layer growth over a longitudinally curved surface[J]. AIAA Journal,1975,13(11):1448-1453. doi: 10.2514/3.7014
[11] SO R M C,MELLOR G L. Experiment on turbulent boundary layers on a concave wall[J]. Aeronautical Quarterly,1975,26(1):25-40. doi: 10.1017/s0001925900007174
[12] AROLLA S K,DURBIN P A. LES of spatially developing turbulent boundary layer over a concave surface[J]. Journal of Turbulence,2015,16(1):81-99. doi: 10.1080/14685248.2014.959126
[13] MATSUBARA K,MUROMOTO T. Two-point correlation and integral scale of spatially advancing curved channel flow at friction-velocity-based Reynolds number 550[J]. Interna-tional Journal of Heat and Fluid Flow,2019,77:31-39. doi: 10.1016/j.ijheatfluidflow.2019.03.003
[14] SALTAR G,ARAYA G. Reynolds shear stress modeling in turbulent boundary layers subject to very strong favorable pressure gradient[J]. Computers & Fluids,2020,202:104494. doi: 10.1016/j.compfluid.2020.104494
[15] UZUN A,MALIK M R. Simulation of a turbulent flow subjected to favorable and adverse pressure gradients[J]. Theoretical and Computational Fluid Dynamics,2021,35(3):293-329. doi: 10.1007/s00162-020-00558-4
[16] COHEN E,GLOERFELT X. Influence of pressure gradients on wall pressure beneath a turbulent boundary layer[J]. Journal of Fluid Mechanics,2018,838:715-758. doi: 10.1017/jfm.2017.898
[17] VOLCHKOV E P,MAKAROV M S,SAKHNOV A Y. Boundary layer with asymptotic favourable pressure gradient[J]. International Journal of Heat and Mass Transfer,2010,53(13-14):2837-2843. doi: 10.1016/j.ijheatmasstransfer.2010.02.014
[18] VOLINO R J,SCHULTZ M P. Determination of wall shear stress from mean velocity and Reynolds shear stress profiles[J]. Physical Review Fluids,2018,3(3):034606. doi: 10.1103/physrevfluids.3.034606
[19] JONES W P,LAUNDER B E. Some properties of sink-flow turbulent boundary layers[J]. Journal of Fluid Mechanics,1972,56(2):337-351. doi: 10.1017/s0022112072002903
[20] SPALDING D B. A new analytical expression for the drag of a flat plate valid for both the turbulent and laminar regimes[J]. International Journal of Heat and Mass Transfer,1962,5(12):1133-1138. doi: 10.1016/0017-9310(62)90189-8
[21] WILLERT C E. High-speed particle image velocimetry for the efficient measurement of turbulence statistics[J]. Ex-periments in Fluids,2015,56(1):1-17. doi: 10.1007/s00348-014-1892-4
[22] 王建杰, 潘翀, 王晋军. 湍流边界层壁面剪切应力的光学测量及统计特性分析[C]//中国力学大会(CCTAM 2019)论文集. 2019. [23] NGUYEN T D,WELLS J C,NGUYEN C V. Wall shear stress measurement of near-wall flow over inclined and curved boundaries by stereo interfacial particle image velocimetry[J]. International Journal of Heat and Fluid Flow,2010,31(3):442-449. doi: 10.1016/j.ijheatfluidflow.2009.12.002
[24] DIXIT S A,RAMESH O N. Large-scale structures in turbulent and reverse-transitional sink flow boundary layers[J]. Journal of Fluid Mechanics,2010,649:233-273. doi: 10.1017/s0022112009993430
[25] ZHOU J,ADRIAN R J,BALACHANDAR S,et al. Mechanisms for generating coherent packets of hairpin vortices in channel flow[J]. Journal of Fluid Mechanics,1999,387:353-396. doi: 10.1017/s002211209900467x
[26] LEHEW J A,GUALA M,McKEON B J. Time-resolved measurements of coherent structures in the turbulent boundary layer[J]. Experiments in Fluids,2013,54(4):1-16. doi: 10.1007/s00348-013-1508-4
[27] 刘铁峰,王鑫蔚,唐湛棋,等. 超疏水表面对湍流边界层相干结构影响的TRPIV实验研究[J]. 实验流体力学,2019,33(3):90-96. DOI: 10.11729/syltlx20180101 LIU T F,WANG X W,TANG Z Q,et al. TRPIV experimental study of the effect of superhydrophobic surface on the coherent structure of turbulent boundary layer[J]. Journal of Experiments in Fluid Mechanics,2019,33(3):90-96. doi: 10.11729/syltlx20180101
[28] WU Y,CHRISTENSEN K T. Population trends of spanwise vortices in wall turbulence[J]. Journal of Fluid Mechanics,2006,568:555-76. doi: 10.1017/s002211200600259x
[29] ADRIAN R J,MEINHART C D,TOMKINS C D. Vortex organization in the outer region of the turbulent boundary layer[J]. Journal of Fluid Mechanics,2000,422:1-54. doi: 10.1017/s0022112000001580
[30] 苏健,田海平,姜楠. 逆向涡对超疏水壁面减阻影响的TRPIV实验研究[J]. 力学学报,2016,48(5):1033-1039. DOI: 10.6052/0459-1879-16-140 SU J,TIAN H P,JIANG N. Trpiv experimental investi-gation of the effect of retrograde vortex on drag-reduction mechanism over superhydrophobic surfaces[J]. Chinese Journal of Theoretical and Applied Mechanics,2016,48(5):1033-1039. doi: 10.6052/0459-1879-16-140
[31] TANG Z Q,WU Y H,JIA Y X,et al. PIV measurements of a turbulent boundary layer perturbed by a wall-mounted transverse circular cylinder element[J]. Flow, Turbulence and Combustion,2018,100(2):365-389. doi: 10.1007/s10494-017-9852-8
[32] ADRIAN R J. Hairpin vortex organization in wall turbu-lence[J]. Physics of Fluids,2007,19(4):041301. doi: 10.1063/1.2717527
[33] MARUSIC I. On the role of large-scale structures in wall turbulence[J]. Physics of Fluids,2001,13(3):735-743. doi: 10.1063/1.1343480
[34] JAISWAL P,MOREAU S,AVALLONE F,et al. On the use of two-point velocity correlation in wall-pressure models for turbulent flow past a trailing edge under adverse pressure gradient[J]. Physics of Fluids,2020,32(10):105105. doi: 10.1063/5.0021121