Abbott六点中心格式应用研究与探讨.pdfVIP

Abbott六点中心格式应用研究与探讨 郑国栋1’2 荣小红3 黄 东2 (1．武汉大学水电学院武汉430072； 2，广东省水利水电科学研究院 广州5106lo 3．广州市水电建设工程有限公司 广州510600) 摘要 本文对目前国内外一些知名的明渠非恒定流计算商业软件所采用的Abbott六点 中心格式的求解方法、离散格式及边界条件的确定进行了详细的论述和推导。通过实例与其他 成熟算法比较，讨论其适用性和优缺点，以便为国内工程设计、研究等单位采购或引进商业软 件时提供一些有益的借鉴。 关键词 Abbot{六点中心格式，Preissmann格式，柯胡格式 1 前 言 明渠非恒定流计算是水利、交通及航运等部门在设计、规划、防洪减灾等工作中经常面临 的工作。过去，受计算工具的制约，对于菲恒定流计算，人们不得已被迫采取一些简化措施(对 模型方程简化，计算手段、方法上进行简化)来较粗略地近似求解所关心的问题。近20年来随 着计算机技术的快速发展和计算机的普及，计算水力学发展突飞猛进，数学模型(数值试验)也 发展成与原型观测、物理模型实验、理论分析等研究方法并驾齐驱重要手段之一。在国内研究 单位中，目前明渠非恒定流计算仍多以自主开发研制的程序为主，这些程序大都经过实践检 验，实用性较好，较为成熟；缺点是未经包装，没有形成商业化软件，短时间内设计或工程单位 人员接受、掌握困难较大。国外开发的一些商业水力学计算软件，由于其友善的界面，操作简 单，数据处理方便，易于被工程人员掌握，而受到设计和工程单位的青睐。缺点是这些软件一方 面运算属于“黑箱操作”，计算失稳或发散不容易查出原因，另～方面软件受其采用的计算格式 限制．本身适用范围和条件也受一定限制．使用不当。容易出错，而这个缺点是比较隐蔽的。本 文就目前国外明渠非恒定漉商用计算软件多采用的Abbott六点中心格式．详述其离散、求解 果进行比较，来研究探讨其适用性。 2 Abbott六点中心格式 2．】 控制方程 描述水体运动的一维明渠非恒定流Sain卜Venant方程由连续方程与动量方程构成： dr 罢抛鲁=g厦 ·78l· u (2) 警+掣州篓+渊；。i十—F十g以磊十℃蕊5 式中：z，f表示空间坐标和时间坐标；Q、^为断面流量和水位；A、“、尺分别为断面过水面积、 断面平均流速和水力半径；g为重力加速度；q为单位河长的旁侧人流量；e为谢才系数；口垂 向速度分布系数。 2．2 Abbott六点中心格式 Abbott六点中心格式如下图1所示。 图l Abbott六点中心格式 式的有限差分法交替计算Q一和^一points。 图2河道上断面(节点)布置示意图 在连续方程中，Q仅对工求偏导，所以方程容易写成以^一point为中心的形式(见图3)； 动量方程则以Q—points为中心(见图4)。 卜—竺—叫 n+1 n+l n+1 n n j-l j j+l j-l j j+l 图3连续方程Abbott六点中心格式 图4动量方程Abbott六点中心格式 在月+1／2时刻，式(1)中的导数项可以写成 !垒整i±垒垒!!一!垡竺i±垒：=!! (3)