基于MATLAB图像处理算法的风沙两相流PIV图像的分割方法.pdfVIP

第31卷第2期 中 国 沙漠 V01．31No．2 2011年3月 RESEARCH Mar．2011 JOURNALOFDESERT 文章编号：1000-694X(2011)02—367—05 基于MATLAB图像处理算法的风沙两相流 PIV图像的分割方法 蒋缠文，梅凡民’，王晓艳 (西安工程大学环境与化学工程学院，陕西西安710048) 摘要：为了深入理解沙粒起动的力学机制，需要从风洞实验中同步地获得风沙边界层中沙粒群和气流的瞬时速 度信息。基于MATLAB平台的风沙两相流PIV图像分割的方法，利用沙粒群和气流示踪粒子特征粒径的显著差 异，对风沙两相流PIV图像中的气固两相流进行图像分割，可分别获得沙粒群和气流示踪粒子的图像，并进一步对 分离后的两相粒子按粒径进行分级分离。得到各级粒径粒子PIV图像，为进一步恢复气固两相流的流场奠定了基 础。 关键词：图像处理；PIV；风沙两相流；沙粒群；示踪粒子；MATLAB 中国分类号：X169；TPl91-41 文献标识码：A 作为风沙物理学的基本关键的科学问题，学术 粒群和湍流的运动速度、瞬时运动应力等参数，而目 界对于沙粒起动的力学机制还存在争议[1q]，前日益成熟的粒子图像测速技术(PIV)为实现风沙 Bagnold[31提出流体起动(后人常称为aerodynamic两相流同步测速各地提供了可能[10-”]。王大伟 entrainment，空气动力学起动或称为空气动力学卷等[14]以颗粒粒径为主要特征，采用基于K一均值的 吸)和冲击起动机制，并认为跃移系统达到稳态时， DigitalMask自动生成算法和互相关PIV算法来分 沙粒起跳以碰撞冲击作用为主，此观点被目前大多 离风沙两相流。代钦等[15]主要通过气固两相流的 数学者所接受[4—5]。根据沙粒冲击速度很小的事灰度差异实现气固两相流的自动分离。赵玉峰 实，董治宝等[11对冲击起动的作用提出了质疑。 等[163利用小波变换和图像分析技术，成功地计算出 Doorschot等[6]利用数值模拟方法模拟了风雪输运分析图像视场内平均粒径大小，该方法为PIV图像 过程中空气动力学起动和冲击起动的作用。结果表 中分割风沙两相流流场内粒子粒径信息提供了借 明，低风速下雪粒起跳以流体起动为主(壁面的湍流 鉴。 切应力大于或接近流体起动的切应力)，高风速下以 为了进一步认识沙粒运动过程中气固两相流的 冲击起动为主，中等风速条件下二者作用相当，且随 动量输送过程和沙粒起跳的动力学机制，需要从 着雪粒粒度增加，空气动力学起动的作用愈大。此 PIV图像气固两相流的图像中同时能够获得不同粒 结论是否适合风沙流，尚需要实验工作的验证。赵 级的颗粒群的粒度和速度信息，而这需要发展一种 建华等[73通过理论推导提出了冲击起动和气动力起 新的PIV图像气固两相流分离方法。该方法不但 动交替循环的机理。至于升力起动说、压差起动说、 能够分离风沙两相流的运动图像，还需要同时识别 振动起动说、斜面飞升说、猝发起动说、湍流起动说、 不同粒级沙粒群和运动速度。笔者提出一种基于颗 负压起动说和涡旋起动说等各种假说，仅是从某一 粒群粒径分级分离的气固两相流分离方法，该方法 侧面反映了沙粒起动的原因，还不能够系统地揭示 不但可以精细地分离PIV图像中风沙两相流，而且 影响沙粒起跳的主要因素和主要的动力学过 可以获得运动颗粒群粒径和速度场，为深入地理解 程E8-9]。笔者认为，风沙颗粒起动是普通的空气动风沙颗粒起动的力学机制提供了技术基础。 力学起动、典型的冲击起动·、湍流猝发起动等多种起 MATLAB是一种高级的面向对象的编程矩阵 动方式的循环交替的动态