天津大学《实验流体力学》数据处理Matlab示例程序（PIV、热线）

粒子图像测速技术（PIV，Particle Image Velocimetry）、热线测速（HWA，Hot Wire Anemometry）是两种广泛应用的流场测速手段，在天津大学流体力学实验室中经常利用这两种测速技术测得各种流场数据，本文主要介绍利用Matlab对这些数据进行基本处理的方法（数据导入、数据处理、数据导出）。

目录

一、PIV数据处理

二、HWA数据处理

三、文件下载

四、注意（尤其是速度剖面拟合）

---

一、PIV数据处理

PIV数据处理程序有四个部分：

初始化：给定基本信息。

流场信息：计算基本的湍流边界层流场信息，如平均速度场、雷诺数、边界层厚度等，并输出到文件。

统计分析：包括平均速度场、平均湍动能、雷诺应力，并输出到文件。 

函数：Spalding速度剖面函数和它的导数，用于拟合壁面磨擦速度。

---

二、HWA数据处理

HWA数据处理程序有六个部分：

读取数据：读入存放在文件夹中的数据

统计量：利用原始数据求统计量，包括平均速度、二阶矩（平均湍动能）、三阶矩、四阶矩、平坦因子。

降低采样率：后续的自相关、频谱的计算不需要过高的采样率，降低采样率以降低计算量。

计算自相关、频谱：自相关和频谱均可以反应数据的周期性。

绘图：对于以上的数据处理进行可视化。

函数：包括三个函数：对自相关结果求其周期、n阶矩可视化（线性坐标）、n阶矩可视化（对数坐标）。

---

三、文件下载

以上程序（附带小量的示例数据和程序说明），整理在“示例程序.rar”中：

链接：https://pan.baidu.com/s/1i4YNXf9NE3Xz6ndpiHxd6w?pwd=1895

提取码：1895

---

四、注意（尤其是速度剖面拟合）

本例给出的所有数据均为示范数据，尤其是PIV边界层数据，测量时没有把控精度。且由于数据量较少，边界层的数据并未收敛。

拟合边界层速度剖面时，例子给出的是Spalding剖面（同时含有粘性底层和对数律区的剖面），但实际上PIV精度并不能测到粘性底层，甚至不能够测量到缓冲层，因此使用Clauser对数律剖面进行拟合更为适宜。拟合时选取的法向高度范围也不能是全部完整的速度剖面，应当去除内区不准的1~2个点、去除外区尾迹律区的点，只选取对数律的点进行Clauser剖面。

更多关于边界层速度拟合的细节请参阅论文，或与老师进行沟通。