基于加权本质无振荡格式二维溃坝水流数值模拟.pdfVIP

水 利 学 报 年 月 第 卷 第 期 文章编号 基于加权本质无振荡格式的二维溃坝水流数值模拟 魏文礼 郭永涛 西安理工大学水利水电学院 陕西西安 摘要 将加权本质无振荡 格式和 时间离散的思想应用于二维浅 水控制方程的求解中 建立了模拟大坝瞬间全溃或局部溃倒所致的洪水演进过程的数学模型 应用该模型对一维 矩形明渠中大坝瞬间全溃所致的水流运动进行了数值计算 并与理论解进行了比较 证实了方法的可靠性 最后 用该模型预测了矩形河道中大坝瞬间局部溃倒时的洪水演进过程 模拟结果与实际相符 算例表明采用 格 式所建立的高分辨率溃坝模型能够很好地模拟溃坝波的演进过程 关键词 格式 溃坝 数值模拟 中图分类号 文献标识码 研究背景 提出的 格式可以有效地处理间断问题 使之不产生新的极 值点 保持解的单调性 但由于其保持解的总变差不变 致使格式在局部极值点处降为一阶精度 为了 克服 格式的这一缺陷 出现了本质无振荡格式 这类格式 是守恒的 全场具有一致的高阶精度 并且总变差有界 其放松了对总变差施加的约束 可进一步提高 格式的精度 并且全场具有一致的高阶精度 即在极值点处精度不降低 该格式在计算中为避免跨越间 断插值引起数值振荡 在计算中自适应地选择模板 利用牛顿差商 绝对值的 大小判断单元内函数的光滑程度 以此来选择最光滑的模板 但 格式仍有许多不足 如 在进 行模板选择时 为得到 阶精度 需要 个单元 总体覆盖了 个单元 但只有一种形式的模板被使 用 模板选择程序包含的很多 如果 型逻辑判断语句对提高程序的计算速度非常不利 在模 板寻找中必须进行网格的扩展 欲得阶数越高 扩展的网格数越多 使实际计算中的界面 边界处理比较 复杂 格式 的出现则消除了以上 点不足 格式采用所有模板的凸组 合 使得包含间断解的模板的权重自动降为极小 而包含平滑解的模板权重自动加大 合理利用了所有 模板 大大提升了程序的运行速度 并简化了边界条件的处理 本文将 格式应用于二维浅水方程 组的求解中 同时引用 时间离散的思想 数值模拟二维溃坝洪水问题 溃坝模型与 格式求解 二维溃坝水流的控制方程为 式中 ? ? 其中 为水深 为 方向的流速 为 方向的流速 为 方向上的底坡 为 收稿日期 作者简介 魏文礼 男 陕西大荔人 教授 博士 主要从事计算水力学及环境水力学研究 方向上的底坡 ? 为 方向的摩阻坡降? ? 为 方向上的摩阻? 坡降 为曼宁系数 方程离散 将式 离散为 ? ? ? ? 整理后得 ? ? ? ? 式中 为时间步长 为空间步长 上标 表示时间层 下标 表示空间节点