考虑河岸变形的三维水沙数值模拟研究.pdfVIP

第20卷 第3期 水 科 学 进 展 V01．2O．No．3 2OO9年 5月 ADVANCESIN WATER SCIENCE May，2OO9 考虑河岸变形的三维水沙数值模拟研究 假冬冬，邵学军，王 虹，周 刚 (清华大学水沙科学与水利水电工程国家重点实验室，北京 100084) 摘要：将粘性河岸崩塌模拟力学方法与水沙模型相结合，构建了考虑河岸变形的三维数值模型。基于非正交网格， 采用局部网格可动技术处理由河岸崩塌引起的河道摆动过程。模型采用有限体积法对方程进行离散，采用与动量 插值技术相结合的SIMPLEC算法进行模型求解。通过弯道发展过程的计算表明，模型能够有效地模拟出河道的复 杂演变过程。 关 键 词：河岸崩塌；三维模型 ；水沙数值模拟；河床变形 中图分类号：TV147．3 文献标识码：A 文章编号：1001．6791I2OO9)03．0311．07 天然冲积河流总是处于不断地冲淤变化过程中，其冲淤变形既包括沿垂向上的冲淤变化，也表现在平面上 的摆动。据统计，长江荆江河段崩岸速度可达88．4H ，可见河道的平面变形是河床演变重要的一方面，其 研究具有重要意义。沿垂向的冲淤变化及边岸的淤长过程可由传统的水沙模型确定，但对于河道平面上的冲刷 后退过程则需另加考虑才能实现。对于河道平面变形的数值模拟虽已取得了一定的进展l22J，但由于问题的复杂 性，这类研究仍存在一些不足。采用沿水深平均的二维数值模型研究 』，不能准确地反映出水沙运动的三维 性，在弯道处该问题尤为突出，虽然可引进经验公式来反映二次流以及泥沙的横向输移对弯道演变所起的作 用_66_，但其局限性较大；M0ssel—，瑚_7j在其研究中亦指出，采用三维模型来模拟河道的横向变形，会得到比较 满意的结果。近年，一些学者开始对该问题的三维数值模拟进行了探索性研究，但由于河岸崩塌简化程度较 大，采用泥沙水下休止角和内摩擦角来确定边岸崩塌 8【,9J，仅适用于非粘性河岸，而未考虑河岸崩塌的力学模 型¨，不能真实反映河岸的崩塌情况。为此本文借鉴 0s咖 -l崩岸模型，将粘性河岸崩塌力学模型与水沙模 型相结合，构建了能够考虑河岸变形的三维数值模型。结合非正交网格采用局部网格可动技术来处理由河岸崩 塌引起的河道摆动过程。模型采用有限体积法对方程进行离散，采用与动量插值技术相结合的SIMPLEC算法 进行模型求解。通过弯道发展过程的计算表明，模型具备了模拟河道河势演变的能力。 1 三维水沙模型 水流模型求解及边界条件详见文献 [12]，泥沙模型介绍如下。 1．1 模型方程 1．1．1 悬沙输移方程 三维悬移质输沙方程[9]： as Z a + 。 (uS)+(s)+(．s)：( )+(e雾)+(e。)+(叫 ()· 式中 、 、W为沿 、Y、 方向的流速；S为含沙量；cU为相应泥沙沉速，采用根据阻力叠加原则得到的 统一公式 1【3j进行计算；￡为泥沙扩散系数。 收稿 日期：2OO8．11-11 基金项 目：国家自然科学基金资助项 目50221903)；国家科技支撑计划资助项 目(20(~BAB05B03；2006BAB05B02) 作者简介：假冬冬(1982一)，男，广西桂林人 ，博士研究生。主要从事水力学与河流动力学研究。 E—mail：jdd05@mails．tsinghua．edu．erl 312 水 科 学 进 展 第20卷 1．1．2 推移质不平衡输沙方程 推移质不平衡输沙方程_lj： 去(qb)+ + =0 (2) 式中 qb及g 分别为实际推移质输沙率和平衡推移质输沙率；‰及q6分别表示沿 和Y方向的推移质输沙 率；Ls为推移质不平衡输沙距离。 平衡推移质输沙率公式采用vanRUN(1984)公式[1