排沙漏斗悬移质泥沙运动数值模拟.pdf

水 利 学 报 2017年8月 SHUILI XUEBAO 第48卷 第8期 文章编号：0559-9350（2017）08-0986-07 排沙漏斗悬移质泥沙运动数值模拟 肖柏青，戎贵文 （安徽理工大学 地球与环境学院，安徽 淮南 232001） 摘要：排沙漏斗是一种利用立轴旋涡实现水沙分离的二级排沙设施。为了探明排沙漏斗悬移质泥沙运动特征和排 沙机理，采用欧拉-拉格朗日方法模拟了排沙漏斗水气沙三相流动，其中水相和气相采用雷诺应力模型（RSM）和 VOF方法描述，悬移质泥沙用颗粒轨道模型描述。通过对计算的悬移质含沙量垂向分布与实验结果的对比，验证 了数学模型的可靠性。计算结果表明：在漏斗室外侧区域含沙量分布整体上较均匀，在漏斗室内侧区域的近底部 存在一个明显的高浓度带。作者认为，上述排沙漏斗悬移质泥沙含沙量分布特征主要是由排沙漏斗内的二次流引 起的。 关键词：排沙漏斗；悬移质；数值模拟；颗粒轨道模型 中图分类号：TV673 文献标识码：A doi ：10.13243/ki.slxb 1 研究背景 在多沙河流上引水所带来的泥沙问题是水利工程建设需要面对的重大问题。排沙漏斗是一种利 用立轴旋涡实现水沙分离的二级排沙设施，在排沙漏斗中，含沙水流被切向引入底部中心设有排沙 底孔的漏斗室做旋流运动，在各种涡流和重力的作用下泥沙被输送至排沙底孔而排出，经过排沙后 的水流则从漏斗室外围溢出再进入引水渠道以资利用。排沙漏斗具有截沙率高、耗水率低和占地面 积小等优点，受到国内外学者的广泛关注和研究。 目前排沙漏斗的研究仍然以实验为主，主要涉及排沙漏斗的排沙效果影响因素、水流流场结构 ［1-6］ 和结构优化等方面 。随着计算机技术和计算流体力学的迅速发展，数值模拟已经成为水沙多相流 研究的一个重要手段，然而目前排沙漏斗数值模拟研究还很少，且基本上是排沙漏斗清水流场模 ［7］ ［8］ 拟。唐毅等 采用三维贴体坐标变换方法对排沙漏斗水流流场进行了数值计算，Ziaei等 使用涡量定 ［9］ ［10］ 义开边界条件的方法计算了简化的排沙漏斗内水流流动情况，肖柏青等 、Chapokpour等 采用 VOF方法模拟了排沙漏斗水气两相流流动。排沙漏斗的中心位置存在漏斗状的自由水面，这使得在 泥沙数学模型中常用的不考虑动水压力影响的自由水面计算方法难以应用，另外，在自由水面变化 复杂的情况下悬移质泥沙运动的计算过程更容易发散，因此，排沙漏斗的浑水模拟比一般的悬移质 ［11］ 泥沙三维数值模拟更加复杂。Athar等 通过求解溶质扩散方程对排沙漏斗内悬移质的浓度分布进行 了数值计算，水流流场和泥沙扩散系数均由试验数据得出的经验公式给出，这种处理方法避开了对 自由水面的求解，但同时也降低了计算精度，并导致所建泥沙数学模型存在较大的局限性。 本文采用欧拉-拉格朗日方法模拟排沙漏斗水气沙三相流动，水相和气相采用雷诺应力模型 （RSM）和VOF方法描述，悬移质泥沙用颗粒轨道模型（也常被称为离散相模型）描述，综合考虑了泥 沙颗粒受力分析、颗粒湍流扩散、颗粒粒径分布和水流流场结构等影响因素，根据计算结果分析了 收稿日期：2016-05-26；网络出版日期：2017-08-29 网络出版地址：http：///kcms/detail/11.1882.TV1100.004.html 基金项目：国家自然科学基金项目 51579002） 作者简介：肖柏青（1980-），男，湖南衡阳人，博士，讲师，主要从事水力学研究。E-mail：baiqingx@126.com — 986 — 排沙漏斗悬移质含沙量分布特征及其排沙机理。 2 物理模型及实验概况