由PIV的速度场获取压强分布的数值方法.pdfVIP

由 PIV 的速度场获取压强分布的数值方法 王大伟，王元 西安交通大学能动学院流体工程系，陕西西安（710049） E-mail: wangdawei@ 摘 要：论述了 PIV 所测得速度场与同位网格形式的相似性，按照有限容积法对 N－S 方程 进行离散，将由 PIV所获得速度场的速度分量代入离散后的方程，求得测量区域的压强分布， 解决了由于非接触测量所导致的 PIV无法同时获得压强分布的问题。同时，在求解过程中发 现，加密网格和采用高阶精度格式方法可以有效减小由于局部梯度过大引起的离散误差。 关键词：PIV；压强分布；数值方法 中图分类号：O351 1．引言 粒子图像测速（PIV ）技术是上个世纪 80 年代末出现的一种激光测速方法[1,2] 。由于其 具有瞬态、非接触式、整场测量的优点，近年来得到飞速发展，其功能也由最初的单相二维 速度场的测量向多相多维的方向发展[3-8] 。然而正是由于 PIV 采用了非接触测量的方法，无 法获得测量区域的压强值，为这种新的测量技术带来了遗憾。本文利用 PIV 测得流场与同 位网格的相似性，对 N －S 方程进行离散，将 PIV 测得流场的速度值代入方程，求得测量区 域的压强分布。同时，在求解过程中发现，加密网格和采用高阶精度格式方法可以有效减小 由于局部梯度过大引起的离散误差。 2. 原理 2.1 PIV 测得流场与同位网格形式的相似性 图 1 PIV 测得流场与同位网格相似性示意图 Fig.1 Comparability of flow field from PIV and collocated grid PIV 测速可分为两部分，拍摄流场图像和图像处理得到流场中各点的速度分布。在后一 部分的图像处理中，首先就是将图像分为若干个大小相等的正方形查问域。最终得到的速度 是各个查问域的平均流动速度。这恰好与采用内节点法的长方形区域均分同位网格有异曲同 工之妙。两者都是将一个小区域内各种物理量的值赋予其中心节点之上。PIV 图像中的每个 基金项目：教育部博士点基金（20030698004）； 国家、新疆自治区“十五”攻关重大项目（2004BA901A21,200433121）. - 1 - 查问域相当于同位网格中所对应节点的控制容积。在图 1 中，两者的区别在于，与网格计算 结果相比，PIV 得到的结果缺少压强分布，但我们恰好可以利用 PIV 的速度场和网格的计算 格式，得出 PIV 测量区域的压强分布。 2.2 根据 PIV 速度场求解测量区域压力分布 对于稳态不可压缩二维流动，N －S 方程可以简化为： ∂ρ u2 ∂ρ uv ∂p ∂ ⎛ ∂u ⎞ ∂ ⎛ ∂u ⎞ + =− + ⎜ μ ⎟ + ⎜ μ ⎟ (1) ∂x ∂y ∂x ∂x ⎝ ∂x ⎠ ∂y ⎝ ∂y ⎠ ∂ρ uv ∂ρ v2 ∂p ∂ ⎛ ∂v ⎞ ∂ ⎛ ∂v ⎞ + =− + ⎜ μ ⎟ + ⎜ μ ⎟ (2) ∂x ∂y ∂y ∂x ⎝ ∂x ⎠ ∂y ⎝