基于非结构网格的分蓄洪区水沙演进数学模型研究.pdf

第28卷 第4期 长 江 科 学 院 院 报 VolJ28 No．4 2011年 4月 JournalofYangtzeRiverScientfiicResearchInstitute Apr．2 0 1 1 文章编号：1001—5485(2011)04—0075—05 基于非结构网格的分蓄洪区水沙演进数学模型研究 张细兵 ，欧治华。，崔 占峰 。孙 先 (1．长江科学院 河流研究所 ，武汉 430010；2．武汉大学 水利水电学院，武汉 430072； 3．荆州市长江河道管理局 洪湖分局，湖北 洪湖 433200) 摘要：基于非结构三角形网格，采用有限体积算法 ，建立了分蓄洪区水沙演进数学模型。模型网格能很好适应分蓄 洪区复杂的边界条件，应用了DEM数字模型插值技术 ，能快速获取计算域的地形数据，实现了洪水演进过程中泥 沙输移和河床冲淤变形的同步模拟。模型中还对网格划分、数值离散、动边界及干河床处理等关键技术进行 了研 究，并提供 了相应的处理方法。该模型在荆江分洪区得到成功运用 ，较好地复演了1954年分洪过程。洪水计算成 果可为分蓄洪区防洪决策与抢险提供参考，分洪闸后最大冲刷深度计算成果可为工程设计提供参考。 关 键 词 ：非结构网格 ；分蓄洪区；洪水演进；泥沙冲淤；数值模拟 中图分类号 ：TV143 文献标识码 ：A 泥沙输移和河床冲淤变形，从而为分蓄洪区的洪水 调度及防洪抢险提供更为丰富和准确的数据。 1 概 述 河流在 自然或人工溃堤情况下，常常会有水沙 2 非结构网格生成技术 的急变运动发生 ，洪水在洪泛区内的演进将给洪泛 区内人 民生命财产安全构成巨大威胁 。因此，正确 复杂流动问题数值计算 中采用的网格，根据拓 模拟溃坝洪水在洪泛区内的演进过程 ，将给防洪抢 扑结构大致分为结构化网格和非结构化网格 2大 险和防洪决策提供重要的依据。由于 问题的复杂 类。非结构化网格是一种无规则随机 的网格结构， 性，国内外许多学者开展了分蓄洪区洪水演进 的研 其优点为：复杂区域的网格生成 自动化程度高、贴体 究。如武汉大学 的余 明辉 、张小峰…采用矩形 网格 性好 ；流动方程直接在物理空间上的非结构化控制 和有限体积算法，建立 了分蓄洪区洪水演进数学模 体上离散和求解 ，不需要坐标变换和方程变换，物理 型 ；王嘉松 、王志刚 基于矩形网格 ，分别采用有 空间即为计算空间；能够较方便地求解复杂边界的 限差分算法和有限体积算法，对局部溃坝水流演进 流动且更容易实现 网格 自适应处理。 问题进行 了研究；Macehione和 Morelli 提 出了用 非结构化网格生成方法主要有 2种 ：Delaunay 于求解溃坝水流的MaeComack格式。国外许多商 三角形化方法和阵面推进法 。这 2种方法均能实现 业化的软件也具备了分蓄洪区洪水演进功能，如丹 网格生成的自动化；能够实现局部加密；能够加 以推 麦 DHI水环境研究所研制的MIKE—Flood模型能进 广生成三维非结构化网格 。笔者基于阵面推进法， 行分蓄洪区洪水演进的一 、二维嵌套模拟 ；美 国密西 引入 了背景 网格 (BackgroundMesh)概念 ，提出了河 西 比大学研制 的CCHE—Flood能进行溃坝洪水演进 道有 限元网格 自动剖分方法 。 的模拟等。但以上模型大都采用结构化 网格 ，一般 未考虑河床冲淤变形 问题。 3 非结构网格上的数值求解方法 为此