防洪自适应网格生成技术(有限体积法有点).pdf

水 利 学 报 年 月 %’()( *’+,-. 第 期 !# $ $ !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 文章编号： （ ） 55202/5 !# $03106 防洪模拟中的地形自适应网格生成技术 马建明，陆吉康 （中国水利水电科学研究院 防洪减灾研究所，北京 1## ） 摘要：在洪水波演进过程中，地形是主要影响因素之一，根据计算区域内地形变化自动调整网格分布，可以使计算 模型和计算结果能更加细致和精确地反映客观地形地貌。本文分析了江河洪水演进模拟中网格自动生成方面的 需求，探讨了洪水模拟中两类地形自适应问题：突变地形和缓变地形。对于如堤防等阻水建筑物，概括为一种地形 突变，网格自动生成时此类地形的自适应可以转化为边界自适应问题，因为此类地形突变一般均能以连续折线来 表示。而对于缓变地形，本文提出了通过曲面样条函数拟合，形成计算区域内连续光滑的地形曲面，从而通过各点 高程梯度来控制网格尺寸，实现网格生成过程中的地形自适应，即在地势平坦区域加大网格尺寸，而在地势变化剧 烈的区域，加密网格。 关键词：防洪模拟；网格；地形；自适应；曲面样条函数 中图分类号： 文献标识码： 78$31 - 二维水动力学数值模拟技术目前已经广泛应用于防洪、治河、河口海岸工程等领域。网格生成技术 是计算流体力学发展的一个重要分支。 世纪 年代以来，随着有限元方法的发展，有限元网格自动 ! 5 生成成为一个研究热点，产生了许多网格自动生成算法，如：狄洛尼（ ）方法，推进阵面法（ 9:;=? -@0 ABCD EFGH ）等。有限元方法最早应用于固体力学，由于有限元网格属于非结构网格，能够处理不规则 边界，特别适合于几何、物理条件比较复杂的问题，因而逐渐在流体力学中也开始应用。但有限元方法 ［，］ 用于水流计算问题时存在一些局限，迄今在大尺度水体流动计算中尚未得到广泛应用 1 ! 。有限体积 法（ ）汲取了有限差分和有限元的优点，用有限体积法结合非结构网格模拟二维浅水流在 世纪 E8I ! 2 ［］ / 年代已呈明显优势。如何自动生成适用的网格，是有限体积法应用的关键之一 。 在江河防洪中，常常涉及蓄滞洪区运用，堤防溃决后洪水在堤防保护区内的演进分析等。在我国， 大型堤防保护区一般均位于流域中、下游平原区，人口、资产密度高，并且有城镇等重要保护对象，一旦 发生堤防溃决需要在很短的时间内分析洪水的走势，并作出决策，同时还需要实时评估抢险方案，如修 建隔堤等的效果。此时，为了应急指挥决策，需要快速动态生成计算网格，体现应急工程措施位置并分 析其效果。网格剖分属于模型计算前处理范畴，根据情况变化自动调整网格，这种能力叫适应性。一是 几何适应性，如在生成网格和存在动边界时根据域的几何形状作出调整。二是解适应性，在求解过程中 不断根据数值解重新分布节点和局部增减节点。在洪水波演进过程中，地形