考虑泥沙冲淤对水流控制方程的改进与模拟-水科学进展.PDFVIP

第 17 卷 第 2 期 水 科 学 进 展 Vol17 ,No2 2006 年 3 月 　　 ADVANCES IN WATER SCIENCE 　 Mar . , 2006 　 考虑泥沙冲淤对水流控制方程的改进与模拟 崔占峰 , 张小峰 , 冯小香 , 谢作涛 (武汉大学水资源与水电工程科学国家重点实验室 , 湖北 武汉 　430072) 摘要 : 以一维数学模型为例 , 通过数值实验及渭河实际资料的验证 , 发现在含沙量较大的情况下 , 水流连续方程里 的浑水与床面的交换项对流量影响是明显的 , 计算中应予考虑 。解释了在像黄河这样的多沙河流中存在的 “水变 沙”、“沙变水”的现象 。由此可见 , 在模拟含沙量较大河流的水沙数学模型中 , 应完整考虑水流连续方程的微元河 段的水量变化率 、进出通量变化及浑水与床面的交换通量这三项 。 关 　键 　词 : 水流控制方程 ; 高含沙 ; 浑水连续方程 ; 数值试验 ; 数学模型 ( ) 中图分类号 : TV1423 　　　文献标识码 : A 　　　文章编号 : 2006 目前 , 人们开发的大多数泥沙数学模型是在含沙量较小的条件下建立的 , 此时浑水与床面的交换通量相对 于过水断面面积的变化率是一个低阶量 , 可以忽略其影响 , 这种处理在大多数情况下都能满足精度要求[ 1～7 ] 。 但在高含沙量情况下 , 这些模型不再适用 。为解决高含沙问题 , 人们提出了许多新的模型 : Zeng 和 Beck 开发 了 STAND 模型 , 并成功模拟渭河上的高含沙过程[8 ] ; 张红武等分析了高含沙洪水的特点 , 提出通过增加经验 公式耦合水沙方程来模拟黄河下游的高含沙过程 , 取得了较好的成果[9 ] 。Cao 等[ 10～12 ] 采用完整浑水连续方程 , 在矩形河槽上进行了数值模拟 , 研究发现在河床变形剧烈的情况下 , 浑水与床面的交换是水流连续控制方程中 不可忽视的一项 。本文从三维浑水连续方程出发 , 积分得出考虑浑水与床面交换的一维浑水连续方程 , 得到水 流连续方程中A 0/ t 项即为浑水与床面交换项 , 并通过数值实验 , 分析了浑水与床面交换通量的影响因素及 1 其对计算结果的影响 。发现高含沙水流发生大量淤积时会导致流量明显减少 , 即出现黄河上所说的 “水变沙” 现象 。计算结果表明浑水与床面的交换通量考虑与否 , 对计算结果的影响是明显的。 1 　一维水流基本方程及其离散 ( ) 11 　一维水流 浑水 连续微分方程 [ 1] 一维浑水连续方程可以由不同的方法得到 , 如谢鉴衡 通过控制体法推导得到 。本文通过对三维连续方程 积分得到一维浑水连续微分方程 。 ρ