弯曲与分叉渠道溃坝波运动特性的数值预报.docxVIP

弯曲与分叉渠道溃坝波运动特性的 数值预报 (2)式中n为曼宁粗糙系数。3有限体 积TVD格式文献提出了卫星单元概念并定 义了单元的几何拓扑关系，在此基础上，构 造了任意四边形单元的有限体积高分辨率 格式。但采用的是一参数限量函数，本文采 用优选的限量函数，以避免在干河时出现非 物理扭曲。 针对单元i(内部区域 ，边 界)，对(1a)式进行积分，得到积分形式 的方程(3)式中A为区域 的面积，dl为边 界 的弧长，n为边界 的外法向单位向量。 边界 由四条线段组成，上式左端第2项可 以写成(4)其中 为边长， 为第k边外法 向通量。 设向量U在单元内部保持不变， 进一步离散(3)式，得到有限体积基本方程 (5)其中 满足(6)禾IJ用F(U)和G(U)的旋 转不变性，引入旋转变换矩阵 及其逆矩阵 可改写式为(7)其中(8) 记⑺式右端 为，则有(9)采用两步Runge-Kutta法离 散上式，时间精度也可以达到二阶，得到 在每一单元的每一边，其通量按组合 型TVD格式给出(例如对于单元i的第1边) 式中 为单元i和卫星单元1的第l个 右特征向量的算术平均。 采用组合型TVD格 式进行插值，(12)其中 代表单元i和卫 星单元1的第l个特征速度的算术平均， 为二单元特征变量差分的算术平均， 为限 量函数，采用文献优选得到的 MUSC型限量 函数，(13)其中正是计算溃坝波时，介于 缓流与急流之间的临界水深比。 当初始下游 与上游的水深比较大时，宜采用单参数限量 函数，否则，在干底时宜采用双参数限量 函数以避免在坝址附近的非物理扭曲现 象。(14) (15)其中(16) (17)时空步长 比为(18)其中 表示单元i和卫星单元1 形心之间的距离。 4边界条件为了使边界与内点保持同 样精度和采用统一的格式，本文增设域外虚 拟单元。(1)开边界，需视流态确定；(2)固 壁边界，与边界重合的边 .域外虚拟单 元；表示其法向通量和切线通量为零。 其中下标“一1”、“一 2”表示域外单元， “b”表示邻近边界的内部单元（如图1所 示）。不计域外单元的几何尺寸和形状，因 为计算中不需要对其量化。 5数值预报 180°弯道 设有一 180°弯曲渠道，入口与 出口直段长2400m,弯曲段内外半径分别是 300m和900m 渠道宽600m坝址位于2240 m处的直段。设上下游水深分别是 10m和， 曼宁粗糙系数为。一般认为，平均径宽比小 于3,则视为强弯曲，本例平均径宽比为， 显然属于强弯曲类型。 采用上述计算方法， 获得了该弯曲渠道全溃溃坝波的数值模拟 结果。计算历时100s、300s、470s后所得 到的自由水面和速度场分别示于图 2、图3 和图4。可以看出，刚开始时，溃坝波在平 直段均匀推进，进入弯道后，溃坝波的运动 受到两岸固壁的限制，流速方向和水位分布 发生变化。外侧水位略高于内侧水位，且最 前端的流体速度总是有偏离弯道方向推进 的趋势，这应是离心力作用以及渠道宽度较 大的缘故。但再进入平直段后，这种变化逐 渐减小，当经历更长一段时间后，这种差别 逐渐消失，两岸流体以几乎相同的波高和速 度推进，甚至外侧流体推进速度略快。这些 特征与已有的实验结果［14］是一致的。 图2弯道内溃坝流动的自由水面与速度场 图3弯道内溃坝流动的自由水面与速度场 图4弯道内溃坝流动的自由水面与速度场 90°双支分叉渠道有一 90°分叉渠道，坝址 位于920m处，考虑初始上游水位为 10m下 游水位为2m和，分别代表湿底和干底。曼 宁粗糙系数取为。图5（a）、（b）分别是计算 溃决100s时湿底和干底的溃坝波自由水面。 负波向上游传播，正波及下游流体在交汇处 自动进行质量和动量的再分配， 形成主流和 支流，且在角区水位变化剧烈，有旋涡。叉 点处的流动特征非常复杂，而且在不同水 深比情况下的特征也是有显着差别的。 干底 时波高显然比湿底时小而波形平坦。 图5双支分叉渠道溃坝波水面线 （左图为 湿底，右图为干底）45°三支分叉渠道一 45 三支分叉渠道，两条对称分布的分支渠道与 主河道成45°角，主河道长2300m,宽400m 分支河道宽300m坝址位于900m处。同样 取上游水深是10m下游水深是2m和，曼宁 粗糙系数。计算溃决 40s和80s时，水深比 为情况下水位等值线如图 6所示；水深比为 情况下水位等值线如图 7。可以看出，溃坝 波运动特征与上例类似，上游波形和波的推 进速度在不同水深比时具有相似的特征， 在 角区存在旋涡。但交汇处和下游波形及推进 速度却有较大差别，下游小水深情况下与大 水深相比，波形更平缓，而且在到达交汇处 之前推进速度略快。 图6三支分叉湿底渠道溃坝波的水位等值 线 图7三支分叉干底渠道溃坝波的水位等值 线6结语在作者近期研究的基础上，