本文介绍了一种利用高分辨率、高速粒子图像测速技术(PIV)研究边界附近瞬态流动的方法。PIV是一种非侵入式测量技术,通过优化图像与记录特性、激光片光源特性以及分析算法等多个参数约束,可应用于任何光学可访问的流动研究。
本文介绍了一种利用高分辨率、高速粒子图像测速技术(PIV)研究边界附近瞬态流动的方法。PIV是一种非侵入式测量技术,通过优化图像与记录特性、激光片光源特性以及分析算法等多个参数约束,可应用于任何光学可访问的流动研究。
多维瞬态流动在科学、工程和健康科学的许多领域中起着关键作用,但往往尚未被充分理解。这类复杂流动的特性可通过粒子图像测速技术(PIV)进行研究,PIV是一种基于激光的成像技术,适用于光学可访问的流动。尽管存在多种形式的PIV,能够将该技术扩展到超越最初二维平面双分量速度测量的能力,但基本的PIV系统由光源(激光器)、相机、示踪粒子和分析算法组成。通过调整成像与记录参数、光源以及算法,可优化目标流动的记录过程,从而获得有效的速度数据。
常规的粒子图像测速(PIV)研究通常在某一平面内测量二维速度分量,速率为每秒几帧。然而,近年来仪器技术的进步推动了高帧率(> 1 kHz)测量的发展,能够以高时间分辨率解析瞬态流动。因此,高帧率测量使得研究人员能够探究高度瞬态流动结构与动力学的演化过程。这些研究在理解复杂流动基本物理机制方面发挥着关键作用。
本文详细介绍了如何进行高分辨率、高速平面粒子图像测速(PIV),以研究平板表面附近瞬态流动的方法。内容包括调整图像与记录参数、激光片光源特性以及数据处理算法等参数约束的详细步骤,以便将PIV技术适用于任何感兴趣的流动研究。
速度的多维测量以及对流场进行时间追踪的能力为科学、工程和健康科学的许多领域提供了关键信息。在流动成像技术中,应用最广泛的技术之一是粒子图像测速(particle image velocimetry, PIV)。PIV最初是一种平面测量技术,用于获取两个面内速度分量的瞬时图像,现已发展出多种变体,能够实现三维速度分量和体积分量的测量。所有PIV系统均由示踪粒子、一个或多个光源以及一个或多个相机组成。常用的示踪粒子包括固体颗粒或液滴,但流动中固有的气泡也可用作示踪粒子。当示踪粒子被光源照射后,相机将捕获其散射或发射的光。在众多技术变体中1,2,最常见的是以每秒数帧的速率测量平面内的两个速度分量。近年来,新型仪器实现了高帧率测量(> 1 kHz),可在kHz量级的湍流时间尺度上追踪流动过程。
PIV 通过追踪由已知时间间隔分隔的一对图像中粒子群的平均运动来确定速度场。每幅图像被划分为由规则间隔的询问窗口组成的网格。最常见的询问窗口大小为 32 × 32 像素。一种算法计算所有询问窗口的互相关函数,从而得到每个询问窗口对应的一个位移矢量,因此生成一个规则的矢量网格。将位移矢量场除以时间间隔即可确定速度矢量场。
在规划粒子图像测速(PIV)测量时,需要认识到实验参数的选择通常是在相互冲突的要求之间做出的权衡。换句话说,必须仔细设计实验条件,以捕捉与当前研究相关的关键流动特征。Raffel et al. 1 以及 Adrian 和 Westerweel 2 的著作对这些限制条件提供了深入而详尽的讨论。本文重点介绍在当前背景下最为关键的几个方面。
视场(FOV)的选择将决定此处参数设置的起点。在确定相机芯片上的像素数量后,若在互相关处理过程中选择使用32×32像素的 interrogation window(通常具有50%的重叠率),则可进一步确定空间分辨率以及所获得的矢量数量。通常希望每个 interrogation window 内有8–10个示踪粒子,以利于互相关函数的计算。然而,也存在一些特殊算法,例如粒子追踪测速法(PTV)和时间平均相关方法,可用于处理示踪粒子密度较低的情况(每个 interrogation window 仅含1–3个粒子),这在靠近表面成像时尤为常见。需要注意的是,每个 interrogation window 内的速度梯度应尽可能小,以避免该窗口代表性矢量出现偏差。
一个既定的经验法则是,第一帧与第二帧之间示踪粒子的位移不应超过8个像素(即 interrogation window 尺寸的¼),以减少相关过程中的配对损失(即示踪粒子图像在 interrogation window 内从第一帧到第二帧丢失)。因此,两次连续激光脉冲之间的时间间隔(dt)必须相应调整。然而,若将 dt 减小至低于相当于8像素位移的值,将降低速度动态范围,因为速度分辨率的下限约为0.1像素的位移。
与成像平面内8像素的位移类似,速度最高的粒子移动不应超过 ¼ 光片厚度的调整同样是为了减少配对损失。由于两束激光脉冲之间的时间延迟用于确保光片平面内最佳的相关性,因此在此背景下光片厚度是一个可变参数。尽管光强均匀性对于基于强度的测量(如平面激光诱导荧光成像)至关重要,但在本应用中其重要性相对较低。 3近平顶光束轮廓有助于提高PIV质量,尤其是在高分辨率成像时。
通常,可以将关于所研究流动特性的某些假设作为选择实验参数的出发点。随后,可能需要进行探索性实验以优化设置。
本文介绍了如何建立一种粒子图像测速(PIV)实验,以实现对两个速度分量的高帧率成像测量,其空间分辨率足以解析边界层结构。该方法采用高重复频率的TEM 00模态二极管泵浦固体激光器、长距离显微镜和高帧率CMOS相机来实现。文中还包含了一些关于近表面成像的细节说明。
1. 实验室安全
2. 台式装置设置
3. 流路设置
4. 优化实验设置
5. 运行实验
6. 数据处理
其中,P1 和 P2 分别为相关峰中的第一和第二最高值,min 为相关平面中的最小值。 Q 是评估矢量质量的一个指标。Q 将产生最佳矢量的最高相关峰与由第二高相关峰代表的普通相关背景进行比较。Q 接近 1 的矢量表明最高相关峰为假峰。接下来,中值滤波器确定中值矢量u中位数,v中位数一组矢量的平均方向与相邻矢量偏离程度(u均方根, vrms)。如果中间向量(u,v)不符合以下标准,则中值滤波器将拒绝该向量: u中位数 - u均方根 ≤ u ≤ u中位数 + u均方根 和 v中位数 - v均方根 ≤ v ≤ v中位数 + v均方根此外,如果在计算速度矢量时设定了较大的重叠区域,则可能获得一些虚假矢量组。因此,可以移除矢量数量少于指定阈值的矢量组。去除虚假矢量后,可采用矢量填充方法,利用非零的相邻矢量进行插值,以填补空白区域。最后,永久应用掩膜将删除掩膜范围之外的所有矢量。实验装置的照片如图1所示。从连续采集的两幅图像中,靠近壁面的32 × 32像素 interrogation window 内的原始粒子图像如图2所示。图2a中的粒子在图2b中向右移动了2-3个像素,满足“四分之一规则”,即面内和面外粒子位移不应超过 interrogation window 尺寸的四分之一。此外,每个 interrogation window 内的粒子密度应约为8-10个粒子,因为PIV相关算法通过追踪粒子群进行计算。然而,在近壁区PIV研究中,示踪粒子密度通常仅为1-3个粒子。因此,对于低粒子密度的情况,应采用特殊算法,例如粒子追踪测速(particle tracking velocimetry, PTV)算法,该算法可追踪单个粒子1,2,4-6。也可采用时间平均相关方法7,8来应对低粒子密度问题,但这通常会导致时间分辨率的损失。此外,靠近壁面成像时会受到强烈激光反射的影响,可能对PIV相关计算产生不利影响并生成错误的速度矢量。这些强反射还会限制壁面法向方向上第一个有效速度矢量的位置。必须对原始粒子图像进行预处理,以降低来自激光反射等来源的背景噪声影响。在本演示中,第一个有效矢量位于距壁面23 μm处。
在使用PIV相关算法处理原始粒子图像后,应对所得速度矢量场的质量和有效性进行评估。原始矢量场中不可避免地会出现异常矢量,但这些矢量具有一些可辨识的特征。错误的矢量通常出现在表面附近、光片边缘以及流动边界的区域。此外,无效矢量的大小和方向与邻近矢量存在显著差异,且在物理意义上不合理。以本例中的边界层流动为例,有效的速度矢量应指向从左到右,因为粒子位移方向也应如此 图2 表明。此外,由于无滑移条件,壁面附近的流速应降低 9. 所示的瞬时速度场 图3 同时满足这两个物理标准。评估PIV结果有效性的另一个有用指标是确定速度矢量场中每个矢量的矢量选择情况。通常情况下,矢量场应由 >= 95% 首选载体 即 那些无需后期处理的数据,因此可采用稳健的后期处理算法来检测并替换虚假矢量,而不会产生明显的伪影 2. 所示的瞬时矢量场 图3 完全由首选载体组成。
通过检查一系列流动图像的时间序列,高速或电影式PIV测量的重要性变得显而易见。图3展示了记录序列开始、中间和结束时刻的瞬时速度(Vi)和速度脉动(V')矢量场。根据雷诺分解,Vi 是平均速度场(
)与 V' 之和10。在本实验中,
通过在时间上对序列中所有图像进行平均得到。在整个记录序列中,瞬时矢量场非常相似,显示流动从左向右运动。这些结果还表明,流动主要沿水平方向,因为水平速度分量(u)远大于垂直速度分量(v)。脉动矢量场也显示,水平速度脉动(u')大于垂直速度脉动(V')。然而,脉动结果还表明流动正在减速,因为在整个记录序列中,u' 的方向发生了反转。
在记录序列的多个不同时间点处,时间平均和瞬时的u速度剖面如图4所示,结果验证了流动随时间逐渐减慢。u速度剖面通过将靠近壁面的四个相邻矢量列进行平均得到,以提高近壁区域结果的统计显著性。该方法在先前的研究中已有应用6,8。误差棒表示四个相邻矢量列标准差的两倍。最大的误差棒出现在平板表面附近,再次说明在示踪粒子密度较低的区域使用PIV相关算法存在困难。目前已有多种分析算法用于应对低示踪粒子密度的情况,例如PTV方法5,6和时间平均相关算法7,8。

图 1. 台式装置组装。

图 2. 在壁面附近一个 32 × 32 像素的 interrogation 区域中,a) t = 0.2 msec 和 b) t = 0.4 msec 时刻的粒子图像。 interrogation 窗口的物理尺寸为 96 × 96 μm2。

图3. 左侧:瞬时(Vi)速度场,右侧:脉动(V')速度场,分别显示记录序列的起始、中间和结束时刻。 矢量场完全由首选矢量构成。为清晰起见,仅展示矢量场的一个较小子集。Vi 场显示流动方向从左向右,而 V' 场方向相反。请注意,为便于观察,水平方向上仅显示每第四个矢量列。此外,如每幅图像左上角所示,Vi 与 V' 场之间的速度标度不同。

图4. 水平速度(u流动过程中不同时间点的剖面图 时间平均的 u - 剖面以圆形表示。t = 0.1 毫秒剖面上的误差棒代表所有其他时刻的误差棒。时间历程的 u - 曲线显示随时间推移流速下降。
与任何光学流动测量技术一样,在规划高速粒子图像测速(PIV)实验装置时,必须评估各种限制条件,并针对具体的测量任务寻求最佳折中方案。图像放大倍数、帧率、激光片光源特性以及分析算法的选择,均取决于所研究流动的具体细节。如有必要,需先进行探索性测量,以确定能够实现高保真度测量的参数设置。
本文介绍了利用高速粒子图像测速技术(PIV)研究平板表面流动边界层的一般步骤及部分示例结果。实验中以5 kHz的频率连续采集了500帧图像。采用长距离显微镜系统,获得位于平板表面的2.4 × 1.8 mm2视场区域。示踪油滴的高质量照明由脉冲式二极管泵浦固体激光器产生的光束实现,该光束通过光束匀化器扩展为片状光。光束匀化器包含一个由多个小型柱面透镜组成的微透镜阵列以及一个额外集成的望远镜系统。微透镜阵列通过将入射光束分割为多个子光束,在垂直方向上扩展圆形光束。随后,望远镜系统将这些子光束叠加,从而在垂直于光束传播方向的平面内形成光强分布均匀的片状光。图像数据采用PIV互相关算法进行处理。需要注意的是,匀化光束在靠近表面区域工作时尤为有利,但对于本文所述应用而言并非关键条件。
本方法通过使用鲁棒的互相关算法,实现了对流动的非侵入式高分辨率、高速测量。这种高分辨率、高速测量技术的主要优势在于其高空间和时间分辨率,以及识别和追踪流场内结构演化的能力。利用这些技术,Alharbi6 和 Jainski et al.8 已成功实现了对内燃机边界层内涡旋结构的可视化与追踪。这些关键特性使得对高度瞬态流动的结构与动力学特性开展研究成为可能。此外,粒子图像测速技术(PIV)可进一步扩展,不仅限于本文所述的二维双分量(2D-2C)速度场测量,还可实现平面内的三分量(3C)测量(立体PIV)以及三维空间内的测量(层析PIV、扫描PIV、全息PIV)。同时,PIV还可与其他技术结合使用,例如平面激光诱导荧光(PLIF)、滤波瑞利散射(FRS)和热致磷光技术,以实现速度场与其他标量场(如温度、组分浓度、当量比)的同步二维测量11-14。这些基于激光的光学方法可直接应用于研究多种场景下的质量与能量交换过程,例如内燃机中的近壁面流动。
作者无任何利益冲突需要披露。
本材料基于美国国家科学基金会资助项目(资助号:CBET-1032930)支持的研究工作,以及在美国密歇根大学定量激光诊断实验室完成的研究工作。
| 姓名 | 公司 | 目录编号 | 评论 |
|---|---|---|---|
| 高速532 nm Nd:YAG激光器 | Quantronix | 型号:Hawk I | |
| 长距离显微镜(QM-100) | Questar | 型号:QM-100 | |
| 高速CMOS相机(Phantom v7.3) | Vision Research | 型号:Phantom v7.3 | |
| 气溶胶发生器(TSI 9306) | TSI | 型号:9306 | |
| 硅油 | Dow Corning CST 510 | CST 510 流体 | |
| 光束匀化器 | Fraunhofer | 定制部件 | |
| 45 ° 高反射率(HR)532 nm 反射镜 | Laser Optik | 多家供应商 | |
| 光阑 | 多家供应商 | ||
| 校准靶标 | 定制部件 | ||
| PIV记录与处理软件 | LaVision | 软件:Da Vis | |
| 高速控制器(HSC) | LaVision | ||
| 光学导轨及滑块 | 多家供应商 | ||
| 激光束阻挡块与吸收阱 | 多家供应商 | ||
| 光学元件支架 | 多家供应商 | ||
| 平移台 | Newport | ||
| 用于产生射流的金属管 | McMaster-Carr | 多家供应商 | |
| 组合方尺与对中方尺 | 多家供应商 |