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