河网模拟JPWSPC方法和分级解法的对比.pdfVIP

第27卷 第10期 长 江 科 学 院 院 报 Vo1．27 No．10 2010年 10月 JournalofYangtzeRiverScientificResearchInstitute Oct． 2 0 l 0 文章编号：1001—5485(2010)10—0039—05 河网模拟 JPWSPC方法和分级解法的对比 朱德军，陈永灿，刘昭伟 ，王智勇 (清华大学 水利水电工程系，北京 100084) 摘要：比较了汊点水位预测一校正法 (junction—pointwaterstagepredictionandcorrection，JPWSPC)和经典的分级解 法，在处理缓流河网汊点处回流效应时的不同，并对 比了它们的计算效率。应用这2类方法时，Saint—Venant方程组 都采用Preissmann格式离散，生成的非线性离散方程用Newton—Raphson方法求解。比较表明：这2类方法都能处理 普适河网，JPWSPC法无需求解整体连接矩阵，同时不会增加每一时间步的迭代次数，因而节约了系统内存 ，提高了 计算效率；河网中河段数 目越多，JPWSPC法的效率优势越明显。 关 键 词：河网水动力数值模拟；汊点水位预测一校正法(JPWSPC)；分级解法；效率 中图分类号：TV131 文献标识码 ：A 尺寸大大降低，因而提高了计算效率。 Islam等人 ¨曾比较了2种分级解法，我国学者 1 概 述 提出的三级和四级等改进解法，都可以建立在这类 隐式差分法稳定性好，能允许较大的时间步长， 分级解法之上。但是，吕满英等人 发现，当初始 在一维水动力数值模拟中应用最为广泛，应用隐式差 误差 比较大时分级解法容易失稳。另外，当河网中 分法离散控制方程后，将形成离散方程组，在每一时 河段数 目很多时，第二步中的矩阵尺寸依然很大，由 间步都需要求解此方程组 。单一河道离散方程组的 于该矩阵是不规则的稀疏矩阵，求解相当耗时。文 系数矩阵有明显的带状特征，大大节省了计算量。当 献 [2]中，陈永灿和朱德军等人首次提出汊点水位 众多河道相互连接形成河网时，如果流动为缓流，由 预测校正 (JPWSPC)法，利用缓流的特征处理汊点 于回水效应，河网中各分支河道互相影响，必须对河 处的回水效应。该方法具有分级解法普适性强的优 网中的所有支流同时进行模拟¨J。与此同时，河网连 点，无需建立和求解总体矩阵。本文将JPWSPC法 接条件破坏了离散矩阵的带状特征，使得求解大型河 和Islam等人选取的2个分级解法进行对 比，重点 网系数矩阵所需的内存和计算量都非常巨大。 分析了他们的效率差别。 针对这一问题，诞生了很多处理方法 (详见文 献 [1]和[2]中的综述)，其 中分级解法 既能求 2 河网控制方程组 解树状河网，又能求解环状河网，被大量学者和工程 河网水动力过程由一维圣维南方程组描述，如 人员采用，尤其在我国，围绕此类方法开展了大量研 方程 (1)和 (2)所示 ： 究，并提出了一系列的改进形式，如三级、四级解法 OA／Ot+oQ／ox—q=o ， (1) 等 一 。分级解法主要分为3步 ：第一步，