1. 课程导入与对标意义
在流体力学研究与工程设计中,数值计算与物理实验长期构成两条并行的技术路线。CFD(Computational Fluid Dynamics,计算流体力学)以离散求解Navier–Stokes方程为核心,能够在计算机中重建流场,给出全场、瞬态、可重复的定量结果;PIV(Particle Image Velocimetry,粒子图像测速)则以示踪粒子与激光片光源为手段,通过高速相机连续拍摄粒子图像,经互相关算法反演瞬时速度场。两者一个来自“方程的解”,一个来自“光学的观测”,各有优势也各有盲区。本课程要解决的核心问题,就是让这两条路线在同一物理对象、同一工况、同一坐标系下真正“对得上”,从而互相校验、互相补充。
为什么必须做对标
CFD的可信度取决于模型假设、网格分辨率、湍流模型、边界条件与数值格式。若不做实验校验,仿真结果可能“自洽但错误”——例如圆柱绕流中湍流模型对分离点与涡脱频率的预测偏差,或叶栅通道中端壁二次流被过度耗散。反过来,PIV虽然直接测量真实流动,但受限于示踪粒子跟随性、片光厚度、相机分辨率、采样频率与互相关窗口,存在峰值锁定、速度欠估、空间平均等系统误差,且通常只能获得有限平面或体积内的数据。对标的意义在于:用CFD补全PIV难以覆盖的全场与壁面附近信息,用PIV约束CFD的模型选择与参数设置,最终形成“实验定标—仿真外推—再验证”的闭环。
典型应用场景
- 钝体绕流与涡脱:以圆柱、方柱为对象,对比斯特劳哈尔数St、阻力系数Cd与尾迹涡结构,检验湍流模型与网格。
- 叶轮机械内流:压气机/涡轮叶栅通道中,PIV测量叶片间速度场,CFD预测总压损失与二次流涡系。
- 燃烧与热流:PIV结合PLIF获取速度与组分场,CFD校验燃烧模型与湍流—化学反应相互作用。
- 微流控与生物流动:低雷诺数下μPIV与CFD对比,验证滑移边界与电渗流模型。
- 风工程与建筑绕流:风洞PIV与LES大涡模拟对比,评估钝体气动力与尾流干扰。
整体学习路径
本课程按“先建立共同语言,再打通数据链路,最后落到案例与排障”的顺序展开。第一阶段理解对标的基本概念、误差来源与评价指标;第二阶段掌握几何、坐标系、工况与数据格式的统一方法;第三阶段学习PIV图像处理与CFD后处理的关键操作,包括互相关参数、矢量场修正、网格插值与时间平均;第四阶段进入典型工程案例,完成从实验布置到仿真设置再到定量对比的全过程;第五阶段聚焦不确定度量化、误差溯源与常见问题排查。每一阶段都强调“可复现”,即读者能依据给出的参数、步骤与判断依据,独立完成一次对标。
核心术语与约定
为避免后续章节出现歧义,这里先统一几个高频术语。瞬时速度场指某一时刻的空间速度分布,PIV单次拍摄与CFD单个时间步均可给出;时均速度场指对足够长时间序列做平均,需保证采样时长覆盖主要低频脉动;脉动均方根反映湍流强度,PIV需足够样本量,CFD需足够统计时间。互相关窗口是PIV速度估计的基本单元,窗口越小空间分辨率越高但信噪比越低;网格分辨率是CFD的空间离散尺度,需与PIV矢量间距匹配才能公平比较。坐标系约定:本课程统一采用右手系,主流方向为+x,展向为+z,法向为+y,所有速度分量与涡量符号均以此为准,若实验或仿真采用不同约定,必须在对比前完成转换并明确记录。
一个可操作的入门示例
以二维圆柱绕流为例,圆柱直径D=20 mm,来流速度U∞=1.0 m/s,水为介质,运动黏度ν≈1.0×10⁻⁶ m²/s,则Re=U∞D/ν≈2.0×10⁴,处于亚临界区,尾迹存在明显涡街。PIV布置:激光片光沿圆柱展向中心平面,跨帧时间Δt取50 μs,查询窗口32×32像素、重叠率50%,采样频率200 Hz,采集2000对图像。CFD设置:二维计算域上游10D、下游30D、侧向各10D,壁面加密,采用SST k-ω模型,时间步长0.001 s,统计时长20 s。对标时先比较时均流向速度沿y方向的剖面,再比较涡脱频率对应的St数。若PIV给出St≈0.20,CFD给出St≈0.21,且速度剖面在尾迹区偏差小于5%,则可认为该工况下模型可信;若偏差显著,应优先检查边界条件、网格分辨率与PIV示踪粒子跟随性,而非直接否定某一方。
需要强调的是,对标的目标不是让两条曲线完全重合,而是让差异可解释、可量化、可追溯。当差异落在合成不确定度范围内,或能归因于已知的物理/数值机制时,对标就是成功的。本课程后续将围绕这一判断标准,逐步展开几何统一、数据处理、案例实操与排障方法。