应用Godunov格式模拟复杂河网明渠水流运动.pdfVIP

中国水利水电科学研究院 第13届青年学术交流会论文 文章编号： 应用Godunov格式模拟复杂河网明渠水流运动 张大伟1,2 (1 中国水利水电科学研究院，北京 100038；2 水利部防洪抗旱减灾工程技术研究中心，北京 100038) 摘要：为克服传统河网数学模型无法准确模拟急流或存在水面间断的复杂明渠水流等问题，采用具有良好激波捕捉能 力的 Godunov 格式建立了新的一维河网水动力数学模型。模型采用基于 HLL 格式的有限体积法离散 Saint-Venant 方程 组，采用数值重构的方法使模型整体获得时空二阶精度。在汊点处由质量、能量守恒方程和特征线理论构造汊点方程 组，求解该方程组为各河段提供内边界条件。模型计算无需特殊的河道编码，每个汊点和河段分开单独求解。通过四 个经典的数值算例对模型进行了验证，结果表明模型能够自动处理流态变化，具有良好的精度和稳定性，同时具备很 好的模拟复杂河网水流运动的能力，说明新研发的河网数学模型是成功的。 关键词：河网；数值模拟；Godunov格式 中图分类号：TV 131.2 文献标识码：A 河网是指由一系列河段相互连接形成的整体，河段和河段之间的连接点称为汊点，根据河网的几 何形状，一般可将河网分为树状河网和环状河网。在工程上，河网问题比较常见，如自然流域的水系 和灌渠系统等。河网水流数值模拟问题一直是水利工程界和学术界关注的热点和难点问题，从上世纪 五、六十年代至今，很多学者对该问题进行过研究[1-3] 。在这些数值方法中，隐格式有限差分法由于其 较好的稳定性应用较为广泛。 最初的隐格式是对整个河网的所有断面水力要素值同时求解，耗时较长，计算效率较低。后来逐 渐发展出分级解法、汊点分组解法和汊点的水位预测校正方法等[4-7] ，将隐格式求解效率大大提升。但 同时需要注意的是，这类求解方法在缓流状态下一般应用效果较好，但是当河段中出现急流或激波现 [8] 象时，往往会遭到失败。如 Preissmann 四点格式，在不做修正的情况下不能模拟跨临界流 ,为了保持 模型的稳定，需要进行专门的近似处理。在 FLDWAV 模型中使用了 LPI （Local Partial Inertia ）方法来 [9] 处理跨临界流动，虽然保证了计算稳定，但在水流非恒定性较强时会引起较大的误差 。 Godunov 格式不仅适用于光滑的古典解，同时具备模拟大梯度流动和自动捕捉激波的能力，由于 这种独特的优势，近年来，该格式在一、二维浅水数值模拟领域得到了广泛的应用[10,11] 。目前，Godunov 格式较多用于单河道复杂明渠水流数值模拟中，尚未见到用于河网数值模拟的报道。为了突破传统河 网数学模型在处理复杂水流流态时遇到的困难，本文采用 HLL 格式的有限体积法来离散 Saint-Venant 方程组，利用特征线法构造汊点方程组提供内边界条件，构造了全新的基于 Godunov 格式的河网水流 数学模型。文中采用的计算格式在处理单河道的复杂水流流态时取得了较好的效果，详见文献[12]， 本文重点论述基于该格式开发的新的数学模型处理复杂河网时的适应能力。 1 控制方程及数值离散 1.1 控制方程 描述天然河网一维浅水流动的 Saint-Venant 方程组的向量形式如下[13] ： U F D  S （1） t x 200  Q   0  B 0 Z f