复杂河网水动力数值模型 Hydrodynamic model for complex flow in open-channel networks.pdfVIP

第22卷第2期 水科学进展 V01．22。No．2 201 l 1年3月 ADVANCESINWATERSCIENCE Mar．，201 DOI：CNKI：32—1309／P．20l10313．2236．017 复杂河网水动力数值模型 朱德军，陈永灿，王智勇，刘昭伟 (清华大学水沙科学与水利水电工程国家重点实验室，北京100084) water and 摘要：利用汉点水位预测一校正(junction-pointstageprediction 建了非恒定河网水动力模型。模型既无需特殊的河道编码。又避免了建立和求解总体矩阵。与常用的分级解法模型 相比，保留了其既适用于树状又适用于环状河网的优点，同时克服了分级解法需要建立和求解河网总体矩阵的缺 点，显著提高了稳定性和计算效率。 关键词：树状河网；环状河网；水动力数值模拟；汉点水位预测-校正法 中图分类号：TVl31．2 文献标志码：A 文章编号：lOOl-6791(2011)02-0203-05 隐式差分格式具有良好的稳定性，在河网数值模拟中得到了最广泛的应用。当水流为缓流时，扰动在河 道中双向传播。上下游流态互相影响，这就是所谓回水效应。由于回水效应，河网中汊点处变量值一方面受 到所有汇聚于该汉点的河道流态影响。另一方面又反过来影响着这些河道流态。也就是说，通过汉点的传递 作用，整个河网中各河道互相影响，所以必须对河网中的所有河道同时进行模拟…。 与单一河道易于求解的带状系数矩阵不同，河网离散方程组的系数矩阵是稀疏的非带状。最初的求解方 法，是直接求解整个河网的离散方程组，需要很大的计算机存储量，同时，由于计算耗时长，往往不能满足 适时预测的需要。针对这一问题，目前最普遍采用的方法主要有两类。第一类方法¨刮是通过特殊的河道编 码，将河网的系数矩阵化为类似单一河道的带状矩阵，能有效减少存储量，提高计算效率。这类方法只适用 于树状河网，而不适用于环状河网。第二类方法¨一刮首先提取出河网各河段的首末断面，建立首末断面变 量之间的联系，并与汊点处的连接条件和边界条件联合组成整体矩阵，求解此整体矩阵，再将所得首末段面 变量值作为边界条件，求解各单一的组成河段，最终得到河网中所有断面的变量值。这类方法在我国常称为 分级解法，国内学者围绕此类方法开展了大量研究，并提出了一系列的改进形式¨引。 分级解法既能求解树状河网，又能求解环状河网，但是，杨开林等一1和吕满英等¨训都发现在建立整体 矩阵的过程中，当河道内断面数很多，而且初始误差较大时，由于误差累积作用，首、尾断面变量之间的关 系过于不准确，将导致分级解法计算崩溃。与上述分级解法不同，杨开林和白正裕采用线性变换的方法建立 整体矩阵。对这一问题进行了处理。除此以外，采用分级解法时，每一迭代步都需要建立和求解整体矩阵， 这一过程很复杂耗时，影响了求解效率。Zhu等…1利用缓流的特征处理汊点处的回水效应，首次提出了汊 water and 点水位预测校正(junction-pointstageprediction 网，而且无需建立和求解整体矩阵。本文在JPWSPC法的基础上，构建了一个普适的河网水动力模型，模型 l-03-13 收稿日期：2010-05-07：网络出版时问：201 10313．2236．017．html 网络出版地址：http：／／www．cnki．net／kcms／detail／32．1309．P．201 ；水利部公益性行业专项经费资助项目(200801008) 作者简介：朱德军(1980一)．男．江苏盐城人．博士后．博士．主要从事水环境和水生态问题研究。 E-mail：zhudj02@mails．tsinghua．edu．c