任意决口水激波的理论解模型及计算方法解剖.pdf

维普资讯 第3卷第2期 四川 联 合 大 学 学 报 (工 程 科 学 版 ) v0l【3No．z 1999年 3月 JOURNALOFSICHUANUNIONUNIVERSITY {ENGINE；F．,RLNGSCIENCE；EDITION) Mar1999 博-k生论文 。 ． ／13-- ‘ 任意决 口水激波的理论解模型及计算方法 黄黄国国富富伍伍超超郑郑永永红红 丁JV、／D67／(1-厶7／ (~Jti联合太学高速水力学国家重点实验室 成都 61006_5) 摘 要 本文提出了任意决 口断面下溃墁水誊泣的理论鼻筷型。通过引人分离变量 法、相似性■理论等，文中不仅解决了复杂央口断面承墩墟不舶■折求解的难题，而 且还全面地克服了激被求懈 数值方法 中存在的闻断抹平、■的非物理振动、域址处 “狗嚏 及洪水前锋无穷太瘴速等缺陷。新模型可直接根据决 口断面形状和域上游 、 下瓣的初始水位求■溃央后的承力特性，对此文中给出了求解方法和计算实倒。算 倒表明。本文所提出的模型及其计算方法是成功而有效的。 关t调 任意决口；理论鲁；廿算方法 ． 蜮～ 一 士嘧 、徼，埙 P涝r^．． ’ 八 r^ ． 大坝突然破坏时，在下游河道有一定水体存在的条件下，水流突然下泻，势必产生水 激波。水激波是一种间断渡，它以立波的形式向下游演进。如 图i所示，溃坝流动可分为 4 区 【：(A)区是未受扰动的上游静水区；(B)区为非恒定渐变流区；(C)区为激波演进区。 在该区内水流的速度和水位保持不变 ；(D)区为未受扰动的下游静水区。图中存在3间断 面：(A)区和 (B)区的交界面4～4为弱间断，(B)区和 (C)区的交界面b～b为弱间断； (C)区和(D)区的交界面 c～c为强间断，亦即所谓的激波间断面，它以恒定的行进速度 向下游推进。显然，要完整表达整个溃坝流动，在给定上游库水深 ho和下游河道初始水 深 ^ 的条件下，不仅需求出(B)区和 (C)区的水探，还需确定3个间断的位置。 在现有的渍坝激波求解方法中，绝大多数数值算法因间断问题没有得到妥善处理以 及受到数值方法本身的制约，易出现解的非物理振动、间断被抹平、坝址处 “狗腿 “等诸多 不合理现象，难以进行溃坝水激波的精确模拟；而在解析解方面，长期以来，仅限于矩形渍 口断面下的求解 ]。故寻求任意决 口激波理论解是十分必要和迫切的。 本文提出了适用于任意溃 口断面的水激波理论解模型，并给 出了求解方法和计算实 倒 。 ‘ 收稿 日期 1998—06—10 维普资讯 l14 四』Il联合大学学摄 (工程科学版) 第 3卷 1 经典的溃坝水激波数学模型 1．1 基本方程 (B)区的流动可用圣维南方程来描述 -： + “ + A差=0 (1a) 害+“+g差=g(S。一st) (16) 式中：A为过水断面面积 ；B为过水断面水面宽；￡为时间；为沿河道的坐标 ，朝向水流方 向为正 ；“为过水断面平均流速 ；g为重力加速度 ；S。；为河底坡 降；Si为摩阻坡降。方程 (1)适用于任意断面。 当前有的数值算法从 (C)区到(D)区