河网模拟JPWSPC方法和分级解法的对比-长江科学院.PDFVIP

　 第２７卷第１０期 长　江　科　学　院　院　报 Ｖｏｌ．２７　Ｎｏ．１０　　 　　２０１０年１０月 ＪｏｕｒｎａｌｏｆＹａｎｇｔｚｅＲｉｖｅｒＳｃｉｅｎｔｉｆｉｃＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ Ｏｃｔ．２０ １０　 　 　　文章编号：１００１－５４８５（２０１０）１０－００３９－０５ 河网模拟ＪＰＷＳＰＣ方法和分级解法的对比 朱德军，陈永灿，刘昭伟，王智勇 （清华大学 水利水电工程系，北京　１０００８４） 摘要：比较了汊点水位预测校正法（ｊｕｎｃｔｉｏｎ－ｐｏｉｎｔｗａｔｅｒｓｔａｇｅｐｒｅｄｉｃｔｉｏｎａｎｄｃｏｒｒｅｃｔｉｏｎ，ＪＰＷＳＰＣ）和经典的分级解 法，在处理缓流河网汊点处回流效应时的不同，并对比了它们的计算效率。应用这２类方法时，ＳａｉｎｔＶｅｎａｎｔ方程组 都采用Ｐｒｅｉｓｓｍａｎｎ格式离散，生成的非线性离散方程用ＮｅｗｔｏｎＲａｐｈｓｏｎ方法求解。比较表明：这２类方法都能处理 普适河网，ＪＰＷＳＰＣ法无需求解整体连接矩阵，同时不会增加每一时间步的迭代次数，因而节约了系统内存，提高了 计算效率；河网中河段数目越多，ＪＰＷＳＰＣ法的效率优势越明显。 关　键　词：河网水动力数值模拟；汊点水位预测校正法（ＪＰＷＳＰＣ）；分级解法；效率 中图分类号：ＴＶ１３１　　　文献标识码：Ａ 尺寸大大降低，因而提高了计算效率。 ［１］ １　概　述 Ｉｓｌａｍ等人 曾比较了２种分级解法，我国学者 提出的三级和四级等改进解法，都可以建立在这类 隐式差分法稳定性好，能允许较大的时间步长， ［８］ 分级解法之上。但是，吕满英等人 发现，当初始 在一维水动力数值模拟中应用最为广泛，应用隐式差 误差比较大时分级解法容易失稳。另外，当河网中 分法离散控制方程后，将形成离散方程组，在每一时 河段数目很多时，第二步中的矩阵尺寸依然很大，由 间步都需要求解此方程组。单一河道离散方程组的 于该矩阵是不规则的稀疏矩阵，求解相当耗时。文 系数矩阵有明显的带状特征，大大节省了计算量。当 献［２］中，陈永灿和朱德军等人首次提出汊点水位 众多河道相互连接形成河网时，如果流动为缓流，由 预测校正（ＪＰＷＳＰＣ）法，利用缓流的特征处理汊点 于回水效应，河网中各分支河道互相影响，必须对河 处的回水效应。该方法具有分级解法普适性强的优 网中的所有支流同时进行模拟［１］。与此同时，河网连 点，无需建立和求解总体矩阵。本文将ＪＰＷＳＰＣ法 接条件破坏了离散矩阵的带状特征，使得求解大型河 和Ｉｓｌａｍ等人选取的２个分级解法进行对比，重点 网系数矩阵所需的内存和计算量都非常巨大。 分析了他们的效率差别。 针对这一问题，诞生了很多处理方法（详见文 ［１，３］ 献［１］和［２］中的综述），其中分级解法 既能求 ２　河网控制方程组 解树状河网，又能求解环状河网，被大量学者和工程 人员采用，尤其在我国，围绕此类方法开展了大量研 河网水动力过程由一维圣维南方程组描述，如 方程（１）和（２）所示： 究，并提出了一系列的改进形式，如三级、四级解法 Ａ／ｔ＋ Ｑ／ｘ－ｑ＝０　， （１） 等［４－７］。分级解法主要分为３步：第一步，利用各河