基于双流体模型的三维溃坝洪水模拟.pdfVIP

第 37卷第3期 人 民 黄 河 Vo1．37．No．3 2015年 3月 YELL0W RIVER Mal．，2015 【防洪治河 】 基于双流体模型的三维溃坝洪水模拟 史志鹏 r，张根广 ，王国平 ，何婷婷 ，刘关艳。 (1．江苏建筑职业技术学院，江苏 徐州221116；2．西北农林科技大学，陕西 杨凌 712100； 3．江苏省水利勘察设计研究院有限公司徐州分公司，江苏 徐州221009) 摘 要 ：在溃坝洪水模拟中，泥沙因素的影响不可忽略。针对坝体溃决中的水沙流动，溃 口采用DAMBRK模型中的概化 梯形模型，通过欧拉模型即双流体模型成功模拟了溃坝洪水中的悬移质泥沙颗粒。对比来流分别为清水和含悬移质泥 沙的水沙流、在下游有村庄阻水和无阻水情况下坝面剪切力的大小，证明悬移质泥沙加快 了坝体溃决，下游阻水建筑物 对水沙流的破坏有抑制作用，最后给出了河道中泥沙的淤积分布情况。所有模拟结果均符合实际渍坝规律和环境。 关 键 词：溃坝洪水模拟；溃 口；欧拉模型；水沙流模拟；紊流模型 中图分类号：TV131 文献标志码：A doi：10．3969／j．issn．1000—1379．2015．03．011 近年来，随着极端天气、暴雨的频发，溃坝事故时 用，给出泥沙在坝体附近的分布情况，为溃坝洪水精确 有发生。2014年，江西龙溪镇栖源村华标山塘水库溃 模拟和溃坝洪水破坏能力的分析提供科学依据。 坝，部分村庄与农 田被淹；2013年，山西曲亭水库因管 厂————————————] 涌而溃坝，库水几乎全部溃泄；2013年，新疆联丰水库 溃坝，村庄淹没，多人死伤。据解家毕等…对全国水 I } 1—．——．—．。．．—．—．．．18．．．—．．．．．— 库溃坝统计及分析，在所统计的溃坝坝型中，土坝比例 最高，达到93．00％，其 中均质土坝占85．85％；溃坝坝 高主要集中在 10—20m；溃坝原因中洪水漫顶是最主 要的原因，占47．85％。溃坝事故一旦发生，必然造成 人员、环境和财物的巨大损失，专家学者对此进行了大 量相关的研究工作。曹志先等 对溃坝水流进行了 (e】俯视图 (d)效果图 定床物理模型试验，并通过V0F模型进行了三维数值 图1 溃坝模型三视图及效果图(单位：m) 模拟，捕捉到了溃坝水流三维流场、湍动能等成果；曹 1 数学模型及数值求解方法 洪建等 通过 自主研发的非定长流动求解器，对溃坝 水流和直立方柱的相关作用进行了数值模拟，给出了 1．1 控制方程 三维溃坝波的演化过程。王晓玲等 通过 VOF模型 假设溃坝洪水属于不可压缩黏性流体，以连续性 和三维 紊流模型，对复杂区域的溃坝洪水进行了三维 方程和雷诺时均N—S方程为控制方程 ¨，连续性方 数值模拟，绘制 了溃坝洪水风险图。此外，DAVID 程为 一 等 L5j、WANG等 J、夏军强等 、张红武等 均对三 +V( )： (1) 维溃坝洪水演进进行了数值模拟研究，取得 了一定的 雷诺时均动量方程为 成果，也验证了数值模拟方法的可行性和正确性 ，但是 这些学者的研究多利用的是表达水汽两相流的VOF +V( M)=一 Vp+ g+V