明渠紊流中泥沙颗粒输移的大涡模拟研究.pdfVIP

泥 沙 研 究 2014年 8月 JournalofSedimentResearch 第4期 明渠紊流中泥沙颗粒输移的大涡模拟研究 及春宁，陈威霖，宋晓阳，许栋 (天津大学 水利工程仿真与安全 国家重点实验室 ，天津 300072) 摘要 ：应用大涡模拟、浸入边界法和有限元一离散元联合法 ，数值模拟 了明渠紊流 中2—3层紧密堆积 的离散 泥沙颗粒 的起动和输移 ，计算了不 同输沙强度下 ，泥沙颗粒运动和受力的统计特征、泥沙浓度 、泥沙体积通量 密度和输沙率 ，以及流场的速度剖面和紊动强度等。研究发现 ：泥沙颗粒 的起动和输移对紊流的统计特征有 明显影响；主流区的泥沙浓度和体积通量密度随高度的增加呈指数下降；颗粒所受平均升力与输沙强度关系 不明显，其最大值约为 0．6倍的颗粒浮容重 ；升力的紊动强度随输沙强度等增大而增大，且数值远大于其平均 值 ，紊流结构 引起 的脉动升力是泥沙颗粒起动的主要原因。 关键词 ：离散泥沙颗粒 ；大涡模拟 ；浸入边界法 ；明渠紊动 ；泥沙输运 中图分类号：TV142．1 文献标识码 ：A 文章编号：0468—155X(2O14)04 o01_o9 1 研究背景 泥沙输运研究对于预测河流和河 口泥沙演变具有十分重要 的价值 。虽然学术界对近壁区域紊 流的 “猝发”现象是引起泥沙颗粒起动和输运的主要因素达成共识 ，但对其机理仍未清楚 。当泥沙 粒径较小时 ，泥沙颗粒 的运动对周围流场不产生明显的影 响，近壁处的泥沙颗粒汇集在低流速区域 ，呈 条带状分布 ，并被喷发 (Ejection)结构 的两个反向旋转的漩涡带到流场的上部 ，形成跃移 ，进而被主 流区水流挟带输运。然而 ，当泥沙粒径较大时 ，泥沙颗粒的运动对周 围紊流的影响不可忽略。具体为 ： 泥沙颗粒对于尺度与其粒径相似或较小的漩涡具有明显的耗散效应，近壁处的水体紊动强度明显降低 ， 紊流的高、低流速条带结构上移 。泥沙颗粒被扫掠 (Sweep)结构形成 的向下俯冲的高速水体 冲 起 ¨ ，之后或重新 回到床面 (有效重力大于脉动升力)，形成跃移质 ；或悬浮在水 中(有效重力小于脉动 升力)，形成悬移质。 借助于高性能计算技术的发展 ，研究者们逐渐摆脱 了 “采用浓度表示泥沙分布 ，并通过求解对流扩 散方程研究泥沙输运过程”的数值研究方法，转而将泥沙视为单个颗粒，通过求解数 以千万的泥沙颗粒 的运动 (考虑泥沙颗粒与紊流之间的相互作用 ，泥沙颗粒之间、泥沙与床面之间的碰撞 ，泥沙颗粒之间 的遮蔽和咬合作用等)，来获取泥沙运动 的整体特性 ，如输沙率、泥沙浓度分布、泥沙颗粒跳跃长度和高 度 、泥沙颗粒运动速度 ，甚至沙纹、沙垄和沙浪的形成和演变过程 、泥沙颗粒分选等。此方面的数值模拟 研究分为两类 ，一类是泥沙粒径远小于网格长度 ，泥沙颗粒被 当作 “点球 ”处理，泥沙颗粒所受的流体力 通过经验模型来模拟 ，泥沙颗粒运动对周 围流场的影响通过在流体控制方程 中增加附加体积力 (浸入 边界法)实现 ，紊流多通过大涡模拟求解 ，泥沙颗粒之间的接触力通过颗粒离散元计算。此类方法可模 收稿 日期 ：2013一l1 5 基金项 目：天津市青年科学基金资助项 目(12JcQNJcO2600，12JcQNJc0560O)；国家 自然科学基金创新研究群体科学基金资助项 目 ；国家 自然科学基金项 目51009105) 作者简介 ：及春宁 (1978一)，男 ，河北人 ，博士，副教授 ，主要从事泥沙动力学数值模拟研究。E-mail：cnji@tju．edu．Cn 通讯作者 ：许栋。E-mail：xudong@tiu．edu．cn 1 拟大量 (千万量级)泥沙颗粒的运动 ，适合研究较大空间尺度的泥沙演变过程。Nabi等 ¨ ”利用大涡 模拟和颗粒离散元方法数值模拟了明渠紊流中沙纹 、沙垄和沙浪的发展和演变，以及桩柱周 围的局部冲 刷和淤积问题。另一类研究中，网