溃坝流量过程数值探究.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 溃坝流量过程数值研究* ** 王鑫 ，岳志远，曹志先 武汉大学水资源与水电工程科学国家重点实验室，武汉 430072 摘要：本研究运用一维水沙耦合数学模型研究坝体瞬间全溃条件下的溃坝洪水流量过程，并 对以往经验性溃坝流量公式的精度和可靠性进行了初步论证。模型的基本控制方程采用一维 浅水方程；运用能够准确捕捉激波和泥沙运动不连续性的 WAF TVD 数值格式离散方程。本研 究模型无需要求对计算区域进行划分，坝址处流量过程作为模型的数值计算结果。研究表明： 坝体溃决过程中，坝址处上下游水位、糙率等因素均会显著地影响坝址处溃坝流量过程，现 有经验性溃坝流量公式忽略了其中一个或诸多影响因素，难以获得可靠、准确的结果。 关键词：溃坝洪水数学模型；溃坝流量过程；激波；瞬间全溃 1 前言 作为一类典型的洪水灾害，溃坝洪水灾害能够给人们的生命和财产造成巨大损失。历史 上发生的溃坝洪水灾害不胜枚举，仅 20 世纪 50 年代以来，印度、意大利、法国、美国、中 国等共有60 余座大、中型水库大坝失事。特别是我国 1975 年 8 月板桥、石漫滩等水库的溃 决，其损失非常惨重。因此，研究溃坝洪水的演进特性进行具有重要意义。 数学模型作为洪水预报的重要工具，几十年以来得到了很大发展。然而，基于限于渐变 流数值方法的传统溃坝数学模型难以捕捉急剧的水流变化过程，在模拟过程中需要假定坝址 处流量过程作为模型的内部边界条件。采用这种处理方法，溃坝流量的可靠性和精度能够严 重影响模型的数值结果，因此有必要对其进行深入研究。物理上，坝址处溃坝流量过程与诸 多因素有关 (如上游水位、下游水位、初始上下游水流流速、河床比降、糙率、坝体溃决过程 等) 。但是，现有的经验性溃坝流量计算公式大多仅考虑水头差的影响[1] ，忽略了一个或诸多 影响因素，其精度和可靠性均难以保证。 近年来，基于能够捕捉激波的高性能数值格式发展了许多溃坝洪水数学模型[2-5] ，其无需 对计算区域进行划分，坝址处流量过程可作为模型的数值计算结果。但是，这些模型大多局 限于定床条件下的清水流动过程，且都没有对坝址处溃坝流量过程进行专门研究。本研究运 用一维水沙耦合数学模型研究坝体瞬间全溃情况下的溃坝洪水流量过程，集中考虑了下游水 位、糙率及床面变形对坝址处溃坝流量的影响，并对以往经验性流量计算公式的精度和可靠 性进行了初步论证。该数学模型曾用于研究溃坝水–沙–河床演变动力学过程[6] ，其性能亦得到 *国家重点基础研究发展(973)计划(2007CB714106)资助. ** wangxincb@ - 980 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 验证。 2 数学模型 2.1 基本控制方程 考虑矩形断面明渠流，根据流体力学基本守恒律[7](质量守恒定律和动量守恒定律)可以获 得一维浅水动力方程，将其整理成守恒形式如下 ∂ ∂ U F + S (1) ∂ ∂ t x 其中U 为守恒向量；F 为通量向量；S 为源项向量，分别定义如下：  ( E −D ) /( 1− p )  h  hu    