水库异重流一维水沙耦合模型.pdf

第２６卷 第 １期 水 科 学 进 展 Ｖｏｌ．２６，Ｎｏ．１ ２０１５年１月　 　 ＡＤＶＡＮＣＥＳ ＩＮ ＷＡＴＥＲ ＳＣＩＥＮＣＥ Ｊａｎ．，２０１５ ＤＯＩ： １０１４０４２／ ｊｃｎｋｉ３２１３０９２０１５０１０１０ 水库异重流一维水沙耦合模型 １ １ １，２ １，２ ２ 王增辉 ， 夏军强 ， 李　 涛 ， 郜国明 ， 张俊华 （１． 武汉大学水资源与水电工程科学国家重点实验室， 湖北 武汉　 ４３００７２； ２． 黄河水利科学研究院， 河南 郑州　 ４５０００３） 摘要： 建立了基于库区不规则断面的一维非恒定异重流数学模型， 并采用明流与异重流水沙输移模型交替运算的 两步模式， 即用潜入条件动态判别异重流计算的上游边界位置， 将潜入点上游的明流浑水段与下游异重流段计算 连接起来。 水流运动、 泥沙输移与河床变形过程完全耦合， 采用ＴＶＤ（Ｔｏｔａｌ Ｖａｒｉａｔｉｏｎ Ｄｉｍｉｎｉｓｈｉｎｇ）形式的ＭＵＳＣＬ⁃ Ｈａｎｃｏｃｋ格式进行数值求解。 将该模型应用于恒定流量与释放定量悬沙两种条件下的异重流水槽实验模拟， 比较了 有无水面梯度项对模拟精度的影响， 计算结果表明该模型能较为准确地预测异重流的厚度、 含沙量分布及传播 过程。 关键词： 水库异重流； 数值模拟； 耦合模型； 两步计算模式 中图分类号： ＴＶ１４５２　 　 　 文献标志码： Ａ　 　 　 文章编号： １００１⁃６７９１（２０１５）０１⁃００７４⁃０９ 水库内蓄水体与上游来流间密度差引起的分层流动被称为异重流。 在水库异重流中引起密度差异的原因 一般是温度或含沙量， 本文仅讨论后者引起的异重流。 维持异重流运动的动力是重力和惯性力， 所以入流流 ［１］ 量、 含沙量变化及库底形态等决定了异重流能否形成以及持续运动到坝前 。 异重流输沙过程与浑水明流 有很大不同， 所以异重流的模拟对于多沙河流上的水库泥沙管理具有重要意义。 利用异重流排沙已在一些水 ［２］ 库中得到初步运用， 如小浪底水库 。 但是在这些水库调度方案设计及排沙效果评估中使用较多的仍然是 经验关系或简化的异重流数学模型。 另外， 现有大部分一维异重流数学模型都是针对矩形断面， 仅对室内尺 ［３⁃４］ ［５］ 度的水槽实验资料进行验证 。 Ｂｏｎｎｅｃａｚｅ等 提出的异重流模型中假设雷诺数足够大而忽略了阻力， 对于 ［６］ 垂向泥沙交换也做了过多的简化处理。 Ｐａｒｋｅｒ 等 通过实验得出了水量掺混系数与理查德森数的经验关系 ［７］ 式， 然而在其模型 中清水掺入对下层异重流的影响只体现在连续性方程中， 忽略了动量的贡献项； 也没 ［８］ 有考虑河床冲淤对异重流质量和动量传输的影响。 李书霞等 在异重流运动方程中引入了反映含沙量与流 速沿垂向不均匀分布影响的修正系数， 建立了潜入点密度弗劳德数与体积比含沙量的关系式。 Ｈｕ 和Ｃａｏ［９］ 根据河床冲淤与异重流是否耦合以及床沙上扬通量是否用平均紊动能限制， 总结了６种异重流模型； 通过数 值模拟与时间尺度分析进行了系统的比较， 对不同模型的使用范围给出了建议。 异重流分层的本质特征是流速、 含沙量沿垂向的不均匀分布， 立面二维模型［１０⁃１２］可以较好地表现这种 垂向结构， 但是由于计算成本过高， 目前应用仍局限于对水槽实验过程的模拟。 另外实际水库地形具有明显 的平面二维特征， 对于