溃坝水流数值计算的非结构有限体积模型 Unstructured finite volume model for numerical simulation of dam-break flow.pdfVIP

第22卷第3期 水科学进展 V01．22．No．3 2011年5月 ADVANCESINWATERSCIENCE May，2011 DOI：CNKI：32-1309／P1449．012 溃坝水流数值计算的非结构有限体积模型 宋利祥1，周建中1，王光谦2，王玉春1，廖力1，谢 田1 (1．华中科技大学水电与数字化工程学院，湖北武汉430074；2．清华大学水利水电工程系，北京100084) 摘要：针对溃坝洪水数值计算面临不规则边界和复杂地形等问题，建立了三角形网格下求解二维浅水方程的高精 时间上，采用Hancock预测．校正法。将底高程定义于单元顶点，结合单元水位一体积关系，提高了干湿界面处理 能力。采用单元中心型底坡项近似，并通过构造通量修正项，保证了格式的和谐性。分析并采用改进的半隐式格式 解决了摩阻项可能引起的刚性问题。算例研究结果表明，模型稳定，和谐性好，具有较好的推广应用价值。 关键词：溃坝；二维浅水方程；Godunov格式；源项；干湿界面；非结构网格 中图分类号：TVl31．2 文献标志码：A 文章编号：1001-6791(2011)03-0373-09 中国病险水库多且分布广，其溃决已成为国家防灾减灾工作中的一个突出问题。为分析溃坝洪水风险， 实现致灾过程预报与灾害评估，溃坝水流数值计算¨刮问题已成为研究的热点和重点。溃坝水流的流态较为 复杂，且存在大梯度解，属于间断流动问题。此外，溃坝洪水淹没模拟要求格式具备和谐性，即保持流速为 零、水位为常数的恒定状态。Godunov格式通过构造和求解单元界面处的黎曼问题计算数值通量，能有效处 理复杂混合流态问题和间断流动问题¨。1，结合线性重构和限制器技术，可在光滑解处达到空间二阶精度， 在间断解附近避免数值振荡。为保证格式的和谐性，Liang等¨1在浅水方程中用水位变量替换水深变量，推 预测一校正法只需要计算一次黎曼问题，因此，Hancock预测．校正法的计算效率相对较高；从格式稳定性看， 针对溃坝洪水数值计算面临不规则边界和复杂地形等问题，本文采用三角形网格划分计算域，并在河道 虑了初始条件和不同波型的影响，并在出现干底时对波速计算进行了改进；将底高程定义于单元顶点，解决 了底高程定义于单元中心时需要人为定义单元界面底高程估算方式的不足，也避免了干湿界面计算时干单元 底高程重构，同时，使单元中心型底坡项近似具有二阶精度。进一步，引入斜底单元的水位一体积关系，并 改进相应的水深蕈构方法，提高了干湿界面处理能力，有利于水量守恒。为保证格式的和谐性，构造了通量 修正项。分析了强不规则地形条件下摩阻项可能引起的刚性问题，进而采用改进的半隐式格式处理摩阻项， 该格式既能保证不改变流速分量的方向，也能避免小水深引起的非物理大流速问题。 收稿日期：2010-05-07；网络出版时间：201l-05—14 (200701008) 通信作者：周建中，E·mail：jz．zhou@hust．edu．cn 万方数据 374 水科学进展 第22卷 1 控制方程 溃坝水流控制方程为浅水方程， 其守恒形式为 ——+——+——=J型+塑+堕：s at a髫 ay hu hv 『 o ] 其中 ．2．1I 2 ，G= huv ，s=I u旧小hu+虿g^ gh(S眠一瓢)l huv 舻+扣^2 Lghis。，一％)J 式中 程；rg为糙率；t为