一维河网水流数学模型数值计算方面的探讨.pdfVIP

第9期(总第277期) 吉 林 水 利 2005年9月 [文章编号]1009—2846(2005)09’0038—04 一维河网水流数学模型数值 计算方面的探讨 黄惠明。施春香 (河海大学，江苏 南京210098) [摘要] 本文在考虑了实际计算河段的特性的基础上．采用直接联解一维圣雏南方程组的方法建立了一堆河网水流数学 模型，并以此撂讨了数值计算中的关于求解失型稀疏矩阵方法的选择及糙率的选择等方面的问题。 [关键词] 一雏河网}稀疏矩阵；糙率 [中图分类号]P333．9 [文献标识码] B 稳定的离散方法。本文在满足以上所提出的一些 1 概述 一维水流模型的基本要求的前提下，采用普林斯 曼的差分离散方法进行离散，并考虑实际计算河 20世纪50年代以来，前苏联、北美、西欧段的特性，采用直接联解的方法求解一维圣维南 和我国就已经展开了水流泥沙数值模型的研究与 方程组，并以此探讨了数值计算中的一些问题。 应用。主要是用于计算大型水库的淤积和坝下游 冲刷及河口挡潮闸下游潮波变形引起的泥沙淤积 2研究方法 等方面。然限于当时的条件，计算精度较差，工 作量大，较难推广应用。随着计算机机技术的提 2．1 一维河网水流模型的基本控制方程 高及计算手段的多样化，数值计算的运用才受到 连续方程： 了重视。 ㈣ 一维水流数学模型经过数十年的发展和应用 冀+BⅣ警=g 已较为成熟，国内外已经普遍使用于研究大型水 运动方程： 利枢纽上下游长距离长时段的冲淤变化及河口感 甏+2“豢+(gA—Bu2)鬈一乱2磬+ 潮河段，如珠江三角洲航道网的水沙变化研究。 一维模型便捷、灵活，其可为二维数学模型提供 g掣一0(2) 必要的边界条件，亦可以作为独立的研究手段， 式中：￡为时间坐标，z为空间坐标，Q为 如对长时段、长河段的水流及河床进行来水、来 流量，Z为水位，乱为断面平均流速，，z为糙率， 沙条件和侵蚀基准点条件变化的研究。 A为过流断面面积，B为主流断面宽度，为水面 针对一维数学模型的主要研究对象的地形特 宽度(包括主流断面宽度B及仅起调蓄作用的 性，如长距离长时段多节点河道的水流变化及冲 附加宽度)，R为水力半径，q为旁侧流量。 淤变化、河口感潮河段等，其地形均属于滩潮交 由于此方程组为二元一阶双曲拟线性方程 错、河床沿断面方向起伏不平、水沙交换频繁的 组，因此采用有限差分法求其数值解。为了加大 地段，因此一维数学模型要较为准确的模拟河道 计算时间步长，节省计算时间，在求解上述方程 断面由于水位的涨落变化而导致的地形变化等情 组中，采用了常用的四点全隐式差分格式，即普 况，且要具有足够的计算精度，就需要采用较为 林斯曼格式。 [收稿日期]2005—05—1s [作者简介]