4. 对标前的实验方案设计

PIV实验与CFD对标的成败,往往在实验开始之前就已经决定了大半。如果测量区域选在回流区边缘、相机视角斜视严重、片光过厚导致离面粒子污染相关峰、示踪粒子跟随性不足或采样频率低于流场特征频率,那么后续无论用多么精细的后处理和对标算法,都很难得到可信的对比结论。本节围绕六个关键设计参数展开:测量区域、相机布置、激光片光厚度、示踪粒子浓度、采样频率以及测量不确定度预估,并给出可执行的估算方法与示例参数。

一、测量区域选择:先对齐CFD的“可解释窗口”

测量区域不是越大越好,而是要与CFD中你真正关心的流动结构对齐。选择时先明确三个问题:要捕捉的主导结构是什么(剪切层、回流泡、尾迹涡街、分离再附点);该结构的特征尺度L是多少;CFD网格在该区域的分辨率是否足以支撑逐点对比。经验上,测量区域应至少覆盖2~3倍特征尺度,并在主流方向留出上下游各约0.5L的缓冲,避免边界效应污染统计量。

以典型钝体尾迹为例,若钝体特征尺寸D=50 mm,回流区长度约1.2D,则测量窗口建议取3D×2D(约150 mm×100 mm),使回流区、剪切层和部分尾迹涡都落在视场内。同时要保证视场内没有强反光壁面、没有遮挡光路的支架,否则局部矢量会出现系统性缺失。测量区域一旦确定,应记录其相对于模型坐标系的原点与范围,便于与CFD提取相同空间窗口的数据。

二、相机布置:视场、分辨率与离面位移的三角平衡

相机布置的核心是让“视场—分辨率—景深”三者匹配。设相机传感器为2048×2048像素,镜头焦距f,工作距离Z,则视场宽度W≈(传感器宽度×Z)/f。若要求矢量空间分辨率达到0.5 mm,则每个查问区至少需要16×16像素,对应放大率约为32 pixel/mm,即视场宽度约64 mm。这意味着你不能同时拥有大视场和高分辨率,必须根据对标目标取舍。

对于二维PIV,相机光轴应垂直于激光片光平面,偏差控制在±3°以内,否则面内位移会被投影压缩,产生系统偏差。两台相机做立体PIV时,夹角通常取30°~60°:夹角过小则离面位移灵敏度不足,夹角过大则视场重叠区变小、透视畸变加剧。相机还应尽量靠近片光法线方向布置,并避免镜头正对强反光面。若使用Scheimpflug调节,需保证片光平面、镜头平面和传感器平面满足成像清晰条件。

三、激光片光厚度:在光强与离面污染之间取平衡

片光厚度直接决定离面粒子对相关峰的贡献。厚度越薄,离面粒子越少,但片光能量密度下降,粒子散射光变弱,信噪比降低。一般原则是:片光厚度应小于最小涡结构尺度的1/4,同时保证片光最厚处不超过查问区深度的1/3。对于常规二维PIV,片光厚度常取0.5~2 mm;若测量边界层或小尺度剪切层,可取0.3~0.8 mm。

片光在传播方向上会因高斯光束发散而变厚,因此要记录片光腰位置和厚度分布。若片光腰在视场中心,边缘厚度可能增加30%~50%。建议在实验前用感光纸或光束分析仪沿光路测量3~5个位置的片光厚度,并在数据处理时评估离面粒子比例。若离面粒子比例超过10%,相关峰会出现明显展宽或双峰,需要重新调整片光或缩小查问区。

四、示踪粒子浓度:既要相关峰,也要跟随性

粒子浓度需同时满足两个条件:查问区内粒子数足够多(通常每个查问区至少8~15对粒子),且粒子跟随性足够好。粒子跟随性用斯托克斯数St=τp/τf衡量,其中τp为粒子松弛时间,τf为流场特征时间。一般要求St<0.1,最好小于0.05。对于空气流场,常用DEHS或橄榄油粒子,粒径1~3 μm;对于水流场,常用空心玻璃微珠或聚苯乙烯粒子,粒径10~50 μm。

浓度估算可用粒子图像密度法:若查问区为32×32像素,对应物理尺寸0.5 mm×0.5 mm,片光厚度1 mm,则查问区体积约0.25 mm³。若希望每个查问区有10个粒子,则粒子数浓度约为40个/mm³。实际布撒时,应通过试拍调整浓度,使粒子图像既不过稀(相关峰弱)也不过密(粒子重叠导致匹配错误)。经验上,粒子图像覆盖率在5%~15%之间较为合适。

五、采样频率:不能低于流场特征频率的两倍

采样频率决定了时间分辨能力。若只做时均统计,采样频率可较低,但样本数要足够,通常要求独立样本数不少于1000~2000。若要解析涡脱落等非定常现象,采样频率应至少为特征频率的2倍,工程上常取5~10倍。例如,若涡脱落频率约为200 Hz,则采样频率建议不低于1000 Hz,最好达到2000 Hz。

对于时间分辨PIV,采样频率还需与激光脉冲频率、相机帧率和数据吞吐量匹配。若使用双脉冲激光,脉冲间隔Δt应根据最大预期位移选择:一般要求粒子在查问区内的位移不超过查问区的1/4。若最大流速为10 m/s,放大率为32 pixel/mm,查问区32像素,则最大位移应小于8像素,对应Δt<8/(10×32)=0.025 ms,即脉冲间隔约25 μs。这个计算必须在实验前完成,否则会出现速度矢量截断或相关峰丢失。

六、测量不确定度预估:把误差预算写进方案

PIV测量不确定度主要来自粒子位移估计、标定误差、片光厚度引起的离面误差、粒子跟随性误差和采样统计误差。位移估计误差通常为0.1~0.2像素,对应速度误差约为0.1~0.2像素/Δt。若放大率为32 pixel/mm,Δt=25 μs,则0.1像素对应速度误差约0.125 m/s。标定误差若为0.5%,在10 m/s量级上贡献约0.05 m/s。离面误差与片光厚度和视角有关,立体PIV可将其控制在1%~3%。

建议在实验前填写一份不确定度预算表,把各项误差按平方和根合成,得到总不确定度。若总不确定度超过CFD与实验对比允许的阈值(例如5%),则需要调整相机布置、片光厚度或采样参数。下表给出一个示例预算,数值为合理假设,用于说明方法而非实测结果。

误差来源典型量级对速度的影响控制措施
位移估计0.1~0.2 pixel0.05~0.10 m/s提高放大率、优化查问区
标定误差0.3%~0.8%0.03~0.08 m/s高精度标定板、多位置标定
离面运动1%~3%0.10~0.30 m/s减小片光厚度、立体PIV
粒子跟随性St<0.05<0.05 m/s减小粒径、匹配密度
采样统计1000样本0.02~0.05 m/s增加样本数、延长采集
实验方案设计完成后,建议先做一次“空场试拍”:不开启主流,只布撒粒子并采集100对图像,检查粒子浓度、片光均匀性、相机对焦和背景反光。这一步能在正式实验前暴露大部分布置问题。
PIV实验方案设计:六参数与对标窗口的耦合关系 ① 测量区域选择 覆盖 2~3 倍特征尺度 L 上下游各留 0.5L 缓冲 示例:D=50 mm → 150×100 mm 避开反光壁面与支架遮挡 记录相对模型坐标系原点 ② 相机布置 光轴垂直片光,偏差 ≤ ±3° 立体PIV夹角 30°~60° 放大率 32 pixel/mm 示例 视场 64 mm @ 2048 px Scheimpflug 保证全视场清晰 ③ 激光片光厚度 厚度 < 最小涡尺度 / 4 常规 0.5~2 mm 边界层 0.3~0.8 mm 离面粒子比例 < 10% 沿光路测 3~5 点厚度 ④ 示踪粒子浓度 每查问区 8~15 对粒子 St = τp/τf < 0.05~0.1 空气:DEHS 1~3 μm 水:空心玻璃微珠 10~50 μm 图像覆盖率 5%~15% ⑤ 采样频率 时均:独立样本 ≥ 1000~2000 非定常:≥ 5~10 倍特征频率 示例:200 Hz → 1000~2000 Hz Δt 使位移 < 查问区 1/4 10 m/s → Δt ≈ 25 μs ⑥ 测量不确定度预估 位移估计 0.1~0.2 pixel 标定误差 0.3%~0.8% 离面误差 1%~3% 平方和根合成总不确定度 阈值示例:总不确定度 < 5% 六参数耦合约束 视场 ↔ 分辨率 ↔ 片光厚度 ↔ 采样率 ↔ 不确定度 输出:与CFD同窗口、同尺度、同统计量的对标数据集 含坐标原点、采样参数、不确定度预算与质量标记 设计顺序建议: 先定测量区域 → 再定相机与片光 → 再定粒子与采样 → 最后做不确定度预算 若预算超标,回退调整前序参数

把上述六项参数写成实验方案后,还应做一次“对标可行性检查”:CFD提取的窗口是否与PIV视场一致;CFD的统计时间是否与PIV采样时长匹配;CFD的网格分辨率是否足以解析PIV矢量尺度;双方的速度参考系和坐标原点是否统一。只有这些条件同时满足,后续的对标曲线和误差分析才有意义。

不要等到实验结束才计算不确定度。若在方案阶段发现总不确定度接近或超过对标阈值,应优先调整相机放大率、片光厚度或采样频率,而不是寄希望于后处理修正。PIV的系统误差一旦进入图像,后处理只能有限缓解,无法完全消除。