基于Stereo-PIV 技术的三维发卡涡结构定量测量研究1)

田海平 伊兴睿 钟山 姜楠 张山鹰

∗(太原理工大学力学国家级实验教学示范中心,太原 030024)

†(曼彻斯特大学机械航空与土木工程学院,英国曼彻斯特M13 9PL)

∗∗(天津大学机械工程学院力学系,天津 300350)

引言

相干结构的发现被认为是20 世纪湍流研究的重大进展之一[1],深入研究相干结构的组织形态、发展演化及其与周围流场的动力学特性规律[2],是理解更复杂流动[3-4]及对湍流实现有效控制[5-7]的基础.在众多形式的相干结构中,发卡涡[8-9]引起了人们极大的关注:一是发卡涡的产生、发展和演化涉及壁湍流边界层的各个分区;二是发卡涡与高低速条带、喷射扫掠等猝发事件、发卡涡涡包等不同尺度、不同类型相干结构间均存在紧密关联[10-11],从而对雷诺应力分布、摩擦阻力的产生起着重要作用.

诚然,在实验研究领域,对发卡涡的组织形态和动力学机制的研究已经取得了丰硕的成果[12-15].但大部分实验结果是基于二维平面的瞬态结果经线性随机估计[15]、拓扑平均[4,16]等手段处理得到的,而直接通过实验手段在湍流边界层中定量地测量发卡涡的三维空间结构并分析其动力学特性仍然具有一定的挑战.目前通过3D-3C-PIV 实验技术直接对真实湍流边界层流场进行测量并提取到发卡涡三维结构还不够成熟,难度是普遍存在的,仅在湍斑[17]、人造发卡涡[18]的测量中观测到了发卡涡结构.究其原因是复杂的:一是,目前测量三维流场的主流技术是3D-3C-PIV,这类设备价格高昂,对相机和测试过程要求较高,且数据量大、分析处理耗时较长;二是,湍流边界层流场近壁面速度梯度大,现有PIV 技术在精细重构湍流边界层的缓冲层及以下的区域的流场时存在较大难度[19-21];三是,在真实壁湍流场中存在的某一个发卡涡结构的形状和尺度特征及流场特性均受很多因素的影响,完整的典型发卡涡实际很少,更难以捕捉、提取和分析.

尽管发卡涡来自壁湍流,但对其形态结构及流体动力学机制的研究,并不是只有在壁湍流场中才有答案.反而,在层流边界层中生成完整、可控的发卡涡结构,并通过控制发卡涡的产生位置、时间、频率、尺度、强度等,用以研究发卡涡的演化过程、相邻发卡涡的耦合机制、二次涡的诱导生成机理及控制发卡涡以实现减小壁面摩阻[22]等更具可操作性,层流场中的研究结果对于评估湍流边界层中相关的流体动力学行为,也具重要的借鉴价值.

用人工方式在边界层流动中制造发卡涡,实质上就是要在边界层近壁局部区域产生有形的、无形的、或时有时无的阻碍,来达到流体在近壁面局部区域积聚或受到障碍、绕流卷曲成发卡涡的动力过程.目前,为实现这一过程主要有4 种方式:其一,通过在边界层壁面局部区域不断注入低速流体,低速流体积聚震荡产生发卡涡结构[18,23],其特点是实验时需不断注入外界流体,且产生发卡涡的频率、尺度均不可控;其二,通过绕流固定障碍物(半球[24-25]、球[26]、圆柱[27]、其他形状[28]等)诱导生成发卡涡结构,其特点是:发卡涡的产生频率、形态和尺度与被绕流物体密切相关,除非替换固定障碍物,否则不可调节;其三,通过绕流可动结构(动粗糙元技术[29],压电陶瓷片[30])诱导生成发卡涡结构,其特点是发卡涡的产生频率、形态和尺度与可动结构密切相关,通过控制可动结构的频率、振幅等可实现对发卡涡结构的部分控制;其四,合成射流装置(synthetic jet)产生涡环后在边界层横流流场中形成发卡涡结构[31-33],其特点是发卡涡的产生频率、形态和尺度与合成射流参数密切相关,容易控制.综上所述,合成射流技术生成发卡涡的过程可控程度最高,即便在无3D-3C-PIV的情况下,只需与Stereo-PIV 系统密切配合即可定量测量发卡涡的三维结构,且该方法属于零质量射流技术,对原有的边界层流场无流量进出、无结构性影响,是一种较为理想的人造发卡涡方式,且有利于系统化实验的开展.Koo 和Kang[22]就是通过控制这种人造发卡涡来进行减小壁面摩擦阻力的研究的,三维发卡涡结构的生成过程通过流动显示做了定性分析.Wen 和Tang 等[34-36]也在层流场中通过合成射流产生人造发卡涡的方式研究了不同相位、不同空间位置的两个发卡涡相互耦合的流体动力学特性,以期发现壁湍流中发卡涡融合和诱导的机制,研究中发卡涡的三维结构也用流动显示的方法呈现,而定量分析只能采用2D-PIV 数据或数值模拟的结果.可见,对可控发卡涡三维结构进行定量测量是壁湍流领域一个亟需突破的课题.

本文通过Stereo-PIV 锁相实验技术定量测量了合成射流装置在层流边界层中所产生发卡涡的三维结构,介绍了一套实验研究发卡涡的新方案,对发卡涡结构特征和动力学行为进行了深入刻画,以期为研究发卡涡发展演化、涡系融合、二次诱导等问题提供思路.

1 实验系统及数据处理方法

1.1 实验系统介绍

实验是在英国曼彻斯特大学机械航空与土木工程学院的低速循环水槽实验室进行的.实验系统主要由低速循环水槽、合成射流装置、Stereo-PIV 系统等三大部分组成,实验平台布局如图1 所示.

图1 实验平台布局示意图Fig.1 Layout of experimental platform

1.1.1 低速循环水槽及边界层平板

实验用低速循环水槽长3.64 m,实验段长0.8 m,截面积为0.305 m×0.305 m,流速在0 ∼0.4 m/s 内连续可调.自由来流速度0.1 m/s 时,水槽湍流度为0.8%.实验段侧壁和底部开有玻璃窗口,便于开展流动显示和PIV 实验.实验时水的自由来流速度U∞=0.12 m/s,水流方向为自左向右(图1).

在距离水槽底面0.31 m 高度处,放置一块长1.38 m,宽0.3 m,厚5 mm 的铝合金板材,倒扣在水槽中,其前缘进行椭圆形修型,作层流边界层平板使用.边界层平板上距离前缘0.715 m 处,开有圆孔,能与合成射流装置的射流出口圆盘平面完美配合.平板和合成射流装置作为一个整体可沿流向在水槽上方实现平移.平板前缘不加扰动、合成射流装置不工作时,平板下表面自然发展成层流边界层.

图2 为分别用2D-PIV 测量流法向平面和Stereo-PIV 测量展法向平面获得的层流边界层外尺度无量纲化的平均速度剖面,二者与Blasisus 曲线符合较好,验证了流场和Stereo-PIV 实验的可靠性.

图2 层流边界层平均速度剖面Fig.2 Mean velocity profile of the laminar boundary layer

1.1.2 合成射流装置及其控制系统

合成射流是一种由激励器振动导致空腔内交替吹吸外界流体而产生的非连续射流.实验中用到的合成射流激励器可以在空腔中产生活塞运动,空腔通过孔口与外界流体连通.吹出的流体由于孔口剪切作用形成射流涡环,遇边界层后涡环沿流向拉伸并远离壁面向外抬升,进而形成具有展向涡涡头及两个反向旋转涡腿结构的发卡涡.实验用到的合成射流装置主要包括合成射流腔体、电磁激振器和位移传感器3 部分.

合成射流腔体是一个长方形腔体(如图3 所示),长H=100 mm,宽W=25 mm.腔体上有流体进出的一端连接在一个圆盘平板上;而另一端是隔膜片,可在电磁激振器的作用下,振动片使薄膜在腔体内上下做活塞运动,带动流体进出腔体.圆盘上有N=5个沿展向整齐排列的圆形射流出口,孔深h=5 mm(即圆板厚度),孔径为D0=5 mm,相邻两个射流出口的展向间距∆z=9 mm.圆盘恰好嵌入边界层平板的圆形开孔中,且两者配合的表面光滑无台阶.

图3 合成射流腔体结构Fig.3 Schematic of synthetic jet cavity structure

在合成射流研究中,αA和RV是两个非常重要的无量纲参量[31],其中αA是空腔截面积与射流出口总截面积的比值,定义为

式(1) 中,H为射流腔体的长,W为腔体的宽,N为圆形射流出口的数目,D0为孔径大小.本实验中αA≈25.5.RV是射流出口平均速度和来流速度U∞的比值,定义为

其中,∆是电磁激励器的振幅,f是激励器振动频率.实验中,电磁激励器振动的频率f=2 Hz,通过调节激励器的振幅∆,使RV=0.1,此时合成射流产生的涡环在层流边界层横流运动中产生典型发卡涡结构.

图4 合成射流激励装置Fig.4 Synthetic jet actuator

图4 是合成射流激励装置,展现了电磁激励器的结构模型及位移传感器的安装位置.电磁激励器的振幅∆、频率f等振动参数由Labview 软件控制,带动驱动连杆与射流空腔一端的隔膜片一起运动.受射流腔体内压力变化影响,隔膜片的运动无法与激励器之间实现完全同步、同幅振动.位移传感器能把真实的振动信号及动作延迟参数及时回传到Labview 软件界面上,用以实时监督实验过程及调整振动幅值、相位等信息,以确保控制效果.

1.1.3 Stereo-PIV 系统

实验中使用的是低频 Stereo-PIV 系统,购自TSI 公司,最高采集频率15 Hz,主要包括照明系统 (Nd:YAG LASER,200 mJ/pulse)、双相机系统(2352×1768 pixels)及同步控制系统(TSI Synchronizer 610036).如图1 所示:激光片光经激光反射镜反射后,从水槽底部窗口垂直入射,照亮法展向(y−z)平面.两相机置于水槽两侧同一高度,与被采集平面满足Scheimpflug 成像原则,为避免光在气液界面产生折射,实验过程中分在水槽两侧各加装了一个45◦的水棱镜,也间接确定了实验中两相机镜头光轴的夹角为90◦.双相机、激光装置是安装在三维立体坐标架上的,可对其移动产生精准控制.标定时,考虑标定靶放置和近壁面标定阻碍等实际情况,先标定“目标待测平面”正下方的“虚待测平面”,标定完成后再将双相机同步竖直提升到“目标待测平面”进行数据采集.实验时用的示踪粒子购自Dantec 公司,平均粒径5µm.

实验过程中,Labview 软件将空腔隔膜片与激励器振动的延迟情况通过数据采集卡传给PIV 系统的同步控制器信号触发端,从而实现隔膜片振动、激光光源系统出光和相机采集等工作的协调同步,保证了锁相采集工作的精度.

为遍历发卡涡发展的整个周期,低频 Stereo-PIV 系统相继采集了25 个相位的数据,相位间隔15◦.实验中采样频率f=2 Hz,相机曝光时间间隔为1000 ns.每个相位重构得到的2D-3C 流场共有235×108 (y−z) 个数据点,物理空间大小为95.9 cm×43.9 cm,空间分辨率为0.41 mm.每个相位采集的样本量为600 个粒子图像对,考虑实验中参数调节及恢复流场稳定等耗时,采集整个周期所需的时间跨度约为3.5 h.在控制参数RV=0.1 下,实验中选取了位于圆出口下游的多个不同的测量位置,本文只对射流出口下游x=D0=5 mm 处的实验结果进行分析.

1.2 三维流场重构方法

合成射流装置的运动周期为T,相邻相位间的时间分辨率为T/24.为得到一个完整周期内的三维流场,流向坐标应当将时间参数转换为空间流向坐标,遵循“泰勒冻结假设”的思想,相位间隔时间与相干结构流向迁移速度的乘积可以看作流向上的空间分辨率.众所周知,边界层流动近壁速度梯度大,合成射流所产生流动结构的迁移速度随空间位置变化而变化.本文主要探究Stereo-PIV 锁相实验技术测量三维发卡涡结构的可行性,并研究发卡涡的结构形态和动力学特性,而经时间坐标转化得到的空间流向长度对所关心的问题并不敏感.为简便,这里以来流速度U∞与相位时间间隔T/24 的乘积作为流向上的空间分辨率,即∆x=U∞·T/24.本文在重构发卡涡三维结构时,将25 个相位的平均流场结果按照相位顺序进行倒序排列,得到一个周期T内的三维流场(见图5).因为流向尺度问题,研究结果中未对结果进行无量纲化处理.同时,也应注意到:流向实质上是时间序列的反方向,x轴上的零点没有实际意义.

图5 三维流场重构方法示意图Fig.5 Schematic diagram of 3D flow field reconstruction method

2 实验结果与分析

Stereo-PIV 锁相实验数据,先进行相位平均,再按照图5 方法重构三维流场.最终获取的三维流场物理空间大小为62.5 mm×43.9 mm×95.9 mm,流向、法向和展向上的分辨率分别为2.5 mm、0.41 mm 和0.41 mm,有25 × 108 × 235(x×y×z)个速度矢量.

2.1 全局流场结果及实验效果评估

图6 Q 值等值面显示的发卡涡结构(等值面上是流向涡量ωx 云图)Fig.6 Iso-surfaces of Q-criterion color coded with streamwise vorticity ωx

2.2 单发卡涡结构所在流场的空间结构形态

为了深入研究发卡涡的结构特征和周围流场的特性,下面取最中间发卡涡所在的相位重构三维流场进行细致分析.选定的空间范围为:−60 mm

图7 Q 值等值面显示的单个发卡涡结构(等值面上是流向涡量ωx 云图)Fig.7 Single hairpin vortex structure(iso-surfaces of Q-criterion color coded with streamwise vorticity ωx)

发卡涡结构上红色区域代表的是正向旋转的涡腿结构,蓝色区域代表反向旋转的涡腿结构,与流向涡量空间分布(图8(a))相吻合.同时,涡腿结构位置与展向涡量分布高度一致(图9).从流向涡量切片分布图(图8(b))中易知,与一对反向旋转涡腿相对应的流向涡下面还诱导出一个与之符号相反的二次流向涡或壁面强剪切区域,其法向尺度非常小.两对流向涡量结构呈双层分布,其符号沿法向上下相反,在展向上左右相反.Chen[2,37]和Li 等[5]在数值模拟研究中也发现了相同结构.与3D-3C PIV 相比,也体现出Stereo-PIV 在近壁流场测量上的优势.

图8 一个合成射流周期内相位平均的流向涡量ωx 空间分布Fig.8 Spatial distribution of phase-averaged streamwise vorticity ωx during a synthetic jet actuation cycle

发卡涡的涡头结构位于俯视图(图7(c))中红蓝相间的位置,结合左视图(图7(b)),涡头中心位置坐标约为(x=−38,y=4,z=0).从展向涡量的分布图(图10) 看,射流出口的下游存在一个负的展向涡量集中区域(图10(a)),且沿着流动方向逐渐远离壁面(图10(b)和图10(c)),形成倾斜结构(图10(c)).

图9 一个周期内相位平均的法向涡量ωy 空间分布Fig.9 Spatial distribution of phase-averaged wall-normal vorticity ωy during a synthetic jet actuation cycle

图11 ∼图13 分别展示了单个发卡涡所在三维流场中3 个速度分量等值面的空间结构及法展向(y−z)切平面速度云图的分布情况.在射流出口下游−55 mm

图10 一个周期内相位平均的展向涡量ωz 空间分布Fig.10 Spatial distribution of phase-averaged spanwise vorticity ωz during a synthetic jet actuation cycle

可见,发卡涡的涡头及两个反向旋转的涡腿,使得低速流体在发卡涡涡头上游的中心区域产生积聚形成低速流体,同时低速流体还有远离壁面的趋势,即形成“喷射流体区域”,发生“喷射”事件.从切片云图上看,喷射流体区域两侧是高速流体(图11(b))、有向壁运动的趋势(图12(b)),即发生“扫掠” 运动.同时,发卡涡的涡腿结构位置靠近壁面一侧,流体的展向运动指向发卡涡结构中心轴线,较为强烈,是主要运动(图13(a));而远离壁面一侧,流体的展向运动趋势相反.

图11 一个周期内相位平均的流向速度u 空间分布Fig.11 Spatial distribution of phase-averaged streamwise velocity u during a synthetic jet actuation cycle

2.3 发卡涡及周围流场的运动规律

图12 一个周期内相位平均的法向速度v 空间分布Fig.12 Spatial distribution of phase-averaged wall-normal velocity v during a synthetic jet actuation cycle

图13 一个周期内相位平均的展向速度w 空间分布Fig.13 Spatial distribution of phase-averaged spanwise velocity w during a synthetic jet actuation cycle

本节以发卡涡结构涡头中心处的法向速度和流向速度作为参考基准对三维重构流场进行了重建,着重分析平面上的流动规律及与空间分布之间的关联.图14 分析了z=0 mm 处流法向x−y平面上带速度矢量的展向涡量云图.图中圆圈位置为展向涡头位置,虚线所示区域是负的展向涡量集中的区域,两者没有明显界限.从速度矢量上看,虚线区域上下两部分流体的运动方向发生明显变化,因此虚线所示区域是一个强剪切的区域,该区域与图10(c)中展向涡量结构相对应.图中黑色圆点是流动滞止点的大致位置.强剪切区域的右下方流体是低速流体,结合速度矢量看,该区域发生喷射事件,这与图11 和图12所呈现的关于发卡涡中心区域发生喷射事件的结果相同.而强剪切层外,尤其是滞止点左上方流体区域,u′> 0 且v′< 0,发生扫掠事件.因而,强剪切层是自由来流与涡结构产生的低速喷射流体共同作用的结果.

图14 流法向截面上(z=0 mm)带速度矢量的展向涡量ωz 云图Fig.14 Contour of spanwise vorticity ωz with velocity vectors on x −y plane at z=0 mm

图15 是4 个带有速度矢量的流向速度云图,分别位于x=−55 mm,x=−47 mm,x=−38 mm 和x=−30 mm 的法展向(y−z)截面.对图7 的分析已知,−55 mm 0),到达上部空间后向两侧缓缓流动.这种流动在平面流场上形成了两个反向旋转的涡旋运动,两个涡旋运动的外侧,流体向下运动(v′< 0),形成两个高速流体区域(u′> 0).这两个涡旋运动对应的就是发卡涡结构的两个反旋的涡腿结构.再一次证明,发卡涡存在的区域,展向中心内部发生喷射事件,而两侧发生扫掠事件.在发卡涡涡头位置x=−38 mm 的平面流场仍有相似趋势,而发卡涡存在区域之外位于x=−30 mm 的平面流场上则无此趋势.

图16 是3 个带速度矢量的法向速度云图,分别来自位于y=1.69 mm,y=3.19 mm 和y=4.19 mm 的流展向(x−z)平面.从图16(a)上看,发卡涡所在区域展向中心处的流体有正的法向速度分量(v′> 0),速度矢量显示流体从两侧汇聚,并向上游运动(u′< 0).而展向两侧位置具有两个明显的剪切层(图16(a)中白色虚线所示区域),位置上与流向涡量集中的区域相对应,因而它的形成与准流向涡或发卡涡涡腿结构密切相关.剪切层外侧的流体有负的法向速度(v′<0),速度矢量显示其向下游运动(u′> 0).随着法向位置的增高,剪切层逐渐向内收缩(图16(b)),并消失(图16(c)).同时,发卡涡所在中心区域流体的喷射运动也逐渐减弱.另外,也注意到同一流展向平面的低速流体区域沿着流向其展向的宽度是变化的,在发卡涡涡头附近最窄.

图15 法展向截面上带速度矢量的流向速度u′ 云图Fig.15 Contour of streamwise velocity u′ with velocity vectors on y −z plane

图16 流展向截面上带速度矢量的法向速度v′ 云图Fig.16 Contours of wall-normal velocity v′ with velocity vectors on x −z plane

经3 个平面上的速度矢量分析,可知以发卡涡涡头迁移运动速度为基准,流法向平面内展向涡量集中区域内速度矢量为零的坐标可指示发卡涡结构涡头的位置;法展向平面上两个反向旋转的旋涡中心(速度矢量为零) 的坐标可以指示发卡涡涡腿的法向中心位置;流展向平面上剪切层的位置,可指示发卡涡涡腿在该平面的投影.因此,利用二维流场的相关特征来重构发卡涡三维流场,在理论上具有可操作性.

2.4 关于主发卡涡下游的一对流向涡结构的讨论

从单个发卡涡三维结构的俯视图(图7(c))和侧视图(图7(a)) 上还可以看到发卡涡结构的下游,还有一对反向旋转的流向涡结构,位于发卡涡涡头的下方且未与主发卡涡结构完全分离,流向上大致位于−40 mm

Suponitsky 等[38]对射流涡环横流运动的数值模拟工作中也发现了相同现象,并认为是涡环的在横流运动中会同时生成主、次两个发卡涡.同时,鉴于锁相技术,该位置产生的这对流向涡也有可能是下游发卡涡的近壁部分,因为发卡涡结构在近壁的迁移速度较慢,在法向上也可能与后一个发卡涡主体部分重叠.因此,它的成因需要进一步探究.

3 结论

本文得到以下几点主要结论:

(1)基于Stereo-PIV 锁相实验技术可以在层流边界层中实现对合成射流装置所产生发卡涡三维结构的定量测量.同时,相较于3D-3C PIV,Stereo-PIV 在流场近壁结构精细测量上占有优势.

(2)所测发卡涡结构三维空间流场符合目前对发卡涡、高低速条带、喷射和扫掠事件的基本认知.

(3) 本文还对二次流向涡、展向涡量集中区域的展向涡头、与低速喷射流体相关的汇聚流动和发散流动等相关问题进行了深入的讨论.特别注意到:(a) 流向涡量和展向速度空间结构呈现双层分布特性; (b) 展向涡量集中区域,一部分是发卡涡展向涡头,另一部分是边界层流场来流与发卡涡涡腿间的喷射流体相互作用形成的强剪切区域,且二者无明显区分界限; (c)低速喷射流体的来源是两侧反向旋转涡腿带动近壁流体(涡腿下层流体)向中心产生的“汇聚” 流动,喷射流体遇到强剪切层后在展向上向两侧作“发散”流动,经涡腿结构上方,受涡腿旋转带动成为高速扫掠流体的主要来源.

(4)以发卡涡涡头的迁移运动为基准的三维空间流场,其各平面矢量图上呈现出与发卡涡结构有明显相关性的信息,可作为“用二维脉动流场的相关特征实现重构三维发卡涡”的实验依据,以期在有限的实验条件下得到质量更高、信息更全的发卡涡特征数据.

本文采用低频Stereo-PIV 锁相实验方案定量测得的发卡涡三维结构,其流向尺度不够严谨.因而严格意义上讲,它仍不是真实的发卡涡结构.尽管如此,该技术方案仍然为定量研究可控发卡涡、乃至进一步研究发卡涡发展演化、涡系融合、二次诱导等问题提供了一些思路.为实现严格意义上的发卡涡三维结构的定量测量,下一步考虑用更新颖的技术手段或更先进的测试设备完成对流向尺度的精确测量,以实现对结果进行无量纲化工作,增加研究工作的普适性.与此同时,通过改变无量纲化的控制参数,按需实现对生成发卡涡结构的有效调控.