复杂边界及实际地形上溃坝洪水流动过程模拟 Two-dimensional modelling of dam-break floods over actual terrain with complex geometries using a finite volume method.pdfVIP

第2l卷第3期 水科学进展 V01．21．No．3 2010年5月 ADVANCESINWATERSCIENCE May，2010 复杂边界及实际地形上溃坝洪水流动过程模拟 夏军强1’2，王光谦1，LIN Bin—lian93，谈广鸣2 (1．清华大学水沙科学与水电工程国家重点实验室，北京100084；2．武汉大学水资源与水电工程科学国家重点实验室， Research Centre，Schoolof 3AA。UK) 潮北武汉430072；3．Hydro·EnvironmentalEngineering，CardiffUniversity，CF24 摘要：建立了基于无结构三角网格下采用有限体积法求解的二维水动力学模型，用于模拟溃坝洪水在复杂边界及 格式及时问方向的预测一校正格式，可使模型在时空方向具有二阶计算精度。模型中引入最小水深概念，提出了有 效的于湿界面处理方法。模拟了理想条件下溃坝水流过程。研究不同最小水深取值对于河床上洪水演进的影响，并 用两组简单溃坝水流的水槽试验资料对模型进行验证。采用该模型模拟了实际溃坝洪水的流动过程，所得计算结 果与实测资料及已有模型计算结果较为符合。 关键词：复杂边界；实际地形；二维浅水方程；有限体积法；无结构网格；溃坝洪水；干湿界面 中图分类号：TVl22．4 文献标志码：A 文章编号：1001-6791l2010)03-0289．10 溃坝洪水在于河床上演进时，一般以顺行正波向下游演进，以逆行负波向上游传播。溃坝洪水的演进过 程，一般可以用守恒形式的二维浅水动力学模型来描述。多年来研究者采用数值模拟及室内试验研究溃坝洪 水的流动过程¨引。Bradford等旧’提出了模拟溃坝洪水的二维水动力学模型，该模型能较为准确地模拟干湿 界面的变化。Zhao等¨虬详细地评述了现有的溃坝水流模型，具有TVD特性的有限体积算法被认为是模拟溃 坝洪水流动较为成功的计算方法之一。因此目前该方法已广泛应用于溃坝洪水流动过程的数值模拟。 在有限体积算法中，水深平均的二维浅水方程在每一个控制单元内离散，能较好地保持质量与动量守 恒，即使在水流条件不连续处。通过界面的法向水流通量一般采用近似Riemann解计算，这种格式的计算效 率较高，又能准确地捕获激波界面。此外该格式通过引入合适的通量限制器，还能消除水跃等间断流附近的 分裂格式及Osher格式，采用不规则四边形网格，模拟了理想条件下溃坝洪水的流动过程。模拟结果表明， 这3种格式在计算精度、计算时间及数值稳定性方面差别不大。Brufau等…。建立了基于单元中心方式的一阶 中，采用一阶精度的有限体积格式难以准确地模拟溃坝洪水在干河床上的演进过程，因此通常采用在时空方 程。因具有守恒性及容易处理实际问题中的复杂边界及地形，目前多采用基于无结构三角网格下具有TVD 特性的有限体积算法模拟溃坝洪水的流动过程。 本文采用无结构三角网格，建立模拟复杂边界及实际地形上溃坝洪水流动的二维有限体积算法模型。该 模型采用Roe．MUSCL格式及时间方向上的预测校正格式，使模型在时空方向具有二阶计算精度，同时采用 收稿日期：2009-05—11 (NECT一10-0619) 作者简介：夏军强(1974一)，男，浙江绍兴人，副研究员，博士，主要从事河流动力学方面的研究。 E-mail：xiajq@tsinghua．edu．Cli 万方数据 290 水科学进展 第2I卷 非常有效的千湿界面处理方法。然后采用该模型模拟理想条件及试验水槽内溃坝洪水的流动过程，以及复杂 边界及实际地形上溃坝洪水的流动过程，并研究了处理干湿界面时不同最小水深取值对计算结果的影响。 1 二维浅水控制方程 溃坝洪水在河流、渠道及洪泛区内的演进，一般可以用二维浅水控制方程来描述，但通常必须满足以下 假定，如静水压力分布、自由水面、垂向加速度