北京航空航天大学五系流体力学实验PIV实验报告.docxVIP

研究生《流体力学实验》—— 粒子成像测速（PIV）技术实验报告 班 级 姓 名 实验日期3月23日 指导教师北京航空航天大学流体力学研究所实验目的1. 利用粒子成像测速技术测量二维流场速度分布。2. 利用matlab中FFT实行互相关运算记录流场速度分布。基本原理粒子成像测速技术（Particle Image Velocimetry）是一种全流场测速技术，可测得流场中某一截面上的瞬时二维速度矢量分布，体视PIV可获取三维速度分量。三、实验步骤1 .实验流程如下图2 .系统构成：粒子及投放装置；双脉冲激光器；图像记录设备；信息处理系统①光源系统：激光器+片光系统=激光片光要求：要求短时间内(脉冲宽度 5ns)保证大量的光能 (20 mJ--500 mJ)；脉冲时间间隔能够视流速大小而变化，范围要求：1μs --若干 ms；脉冲激光器典型重复频率为 10-30Hz (激光器蓄能需要时间)；双脉冲Nd:YAG激光器，脉冲时间间隔Δt = 1-150 μs ,重复频率15 Hz，适合高速气流速度的测量。其中：?t: 脉冲时间间隔 T: 单个激光器脉冲重复时间 τ: 脉冲宽度②图像记录设备（CCD 相机）CCD — Charge Coupled Device：电荷耦合元件或CCD图像传感器，CCD是一种半导体器件，能够把光学影像转化为数字信号。CCD上植入的微小光敏物质称作像素（Pixel）。一块CCD上包含的像素数越多，其提供的画面分辨率也就越高。CCD的作用就像胶片一样，但它是把图像像素转换成数字信号。CCD 阵列的空间分辨率至少比摄影胶片低两个量级。CCD按加工工艺分为两种：TTL和CMOS，TTL工艺成像质量要优于CMOS工艺；CMOS成像和信息存储、传输速度快，用于高速高频响PIV。而CCD的信号存储和传输：像素→ 存储区: 500 ns，存储区 → PC: 33 ms，跨帧技术如下图所示因此对CCD相机的性能要求：记录图像序列，一帧帧连续排列；高空间分辨率；高速捕获多帧图像； 高感光度。③信息处理系统：计算机实验注意事项①示踪粒子选择两个依据：散射光强度和粒子跟随性粒子光散射特性如下图：部分光在四周方向散射光强:由于光的漫射大约正比于1/r2，故希望粒径大。粒子跟随性粒子跟随性粒子的动态响应，包括阶梯响应时间常数，假设粒子为球形，低雷诺数假设下，Stocks方程求解得到球形粒子的受力公式为设粒子从静止到达来流速度U? 的响应时间为，则及与流体密度不同引起的粒子沉浮速度：故从此方面希望粒径小点。故要求液体中粒径 10-20 μm综合性能好，常用粒子有聚苯乙烯 (10-100 μm); 铝粉(2-7 μm); 空心玻璃球 (10-100 μm).②相机调节和标定在PIV 实验中，会同时有多台相机分别以一定视角对实验区域成像。当相机以一定视角对着被摄区域时，被摄平面的两端到镜头的距离是不同的。由于镜头具有景深，对焦时可能被摄平面的两端无法同时落入镜头的景深内，就会导致成像总有一部分区域是模糊的。这时就需要调节镜头的角度，使被摄平面、镜头平面和像平面满足Scheimpflug 条件，就可以解决景深对有视角成像的不利影响，看到整幅清晰的图像。③matlab实现平面互相关算法互相关函数的计算有两种途径：直接互相关计算，计算量 O(N2*N2)和利用FFT进行互相关计算，计算量 O(N2log2N)，本实验利用FFT计算， 要求两个询问区必须是(N×N)的正方形，且N为2的整数幂。当N选取的比所要求的大，会引起空间分辨率的降低，此处用两帧/单重曝光，如图所示要求粒子在脉冲时间间隔内运动的位移必须小于查询区边长的1/4；为提高空间分辨率，可取相邻查询区50%的重叠；每个询问区的粒子数量: 至少3-4个，10-15较佳。互相关的优化包括利用高斯拟合进行亚像素内插。下面是matlab程序：t=1/180;i=1;fangdabeishu=1;type=.tif;for number=1:20if number10filename=C001H001S000100000;elsefilename=C001H001endP(:,:,i)=imread([filename,num2str(number),type]);if i==1i=i+1;elseP1=P(:,:,i-1);P2=P(:,:,i);[height,width]=size(P1);W=32; bushu=32;m=32;n=32;M=(height-W-2*m)/bushu+1;N=(width-W-2*n)/bushu+1;[u,v]=huxiangguanyunsuan(P1,P2,t,fang