环状河网水流的数值计算方法

《环状河网水流的数值计算方法》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、环状河网水流的数值计算方法钟娜王玲玲河海大学水电学院2006.5.27HohaiUniversity模型建立的必要性河网水流的水动力状况，对于平原地区特别是我国东部沿海经济发达地区的经济发展和社会生活具有重要的意义。长江及珠江三角洲城市水域的流通以及水质问题已越来越被人们所重视，对于河网水流水动力条件数值模拟的精度要求也越来越高。由于城市河网交错复杂，大部分以环状形式出现，故本文对环状河网进行建模和水力计算方法研究。模型建立的理论依据对于环状河网的求解，可以利用显式差分求解，但工程应用中为了满足收敛及计算稳定性要求，一般倾向于利用隐式求解。早期曾出现过对断面所有未知量建立方程组直

2、接求解的一级解法，以及以河道的首、末断面水力要素为基本未知量的二级解法，但这二类方法的方程组系数矩阵都过于庞大，很难应用于大型河网，因此，本文采用以节点水位为基本未知量的三级解法。为了使模型具有可扩充性、可移植性和通用性，对于系数矩阵的求解采用矩阵标识法。河网水流模型单一河道水流的模拟计算环状河网模型的建立河网的概化节点编码河道的编码河道流动方向河道断面编码环状河网的水力计算流程算例结论1.河网水流模型1.1单一河道水流的模拟计算河网模型控制方程组为圣维南方程组连续方程（1）动量方程（2）式中，为流量，；为水位，；为过水断面面积，；为动量校正系数；为水面宽，；为水力半径，；为谢才

3、系数，，；为河床的造率；为时间，；为流程，；为单位流程上侧向出流量，，负号表示流入，为重力加速度，。Pressimann格式离散图pressimann四点线性隐格式（3）式中，为权重系数，离散后的圣维南方程组（4）（5）由方程(4)和(5)组成的方程组为常系数线性方程组。对于任意河道，设首断面号为，末断面号为，则每条河道有个河段，有个未知变量，可以列出个方程，结合河道两端的边界条件，形成封闭的代数方程组，可以得到唯一的解。1.2环状河网模型的建立对于环状河网，其节点分为内节点和外节点。为了减少方程组的阶数，通常以内节点构成基本河网。对于环状河网，内河道和外河道因为边界条件的不同而

4、计算的形式不同。外河道的计算由于外边界条件的不同而分为两类：水位型边界和流量型边界。这两种形式采用追赶法求解的递推公式如下：1.2.1水位流量递推关系的建立水位型边界条件：(6a)流量型边界条件：(6b)1.2.1水位流量递推关系的建立由上面两式得出，末断面的流量都可表达成(7)这里，表示与该外河道相连接的基本节点水位。内河道的计算采用双追赶方程求解。(8a)(8b)设为首节点水位，为末节点水位，式(8)表示首、末断面流量可以用首、末断面水位的线性组合求得。当首、末断面的水位已知时，利用式(8a)和(8b)，可以推出下式：(9)求得后代入式（8a）中即可求出。1.2.2交叉节点水

5、位方程的建立节点水位方程建立的依据是水量平衡原理和能量平衡原理。水量平衡原理表示流入某节点的水量之和等于该节点的蓄水量的变化，即：(10)能量平衡原理表示节点所连接的各个河道断面上水位和流量与节点平均水位之间，必须符合实际的能量衔接要求。1.2.2交汊节点水位方程的建立这种衔接一般分为三种情况：如果节点可以概化为一个几何点，出入该节点的水流平缓，不存在水位突变的情况，则各节点断面水位相等，均等于该点的平均水位，即(11a)如果各断面的过水面积相差悬殊，流速有明显的改变，但仍属于缓流的情况，按照Bernouli方程，当略去节点的局部损失时，各个断面之间的能量水头应该相等，即(11b

6、)一般情况下，当节点建有堰、闸等水工建筑物时，节点的断面动力衔接条件可写为：(11c)将流量与节点水位关系式（11）代入水量平衡方程（10），得到与节点相邻的节点水位为未知变量的线性代数方程：(12)对于河网每一个节点，都可建立上述的节点水位方程，形成以河网的节点水位为基本未知变量的线性代数方程组：(13)其中，为系数矩阵，为节点水位列阵，为右端项列阵。为使模型具有可扩充性、可移植性和通用性，节点水位方程的求解采用矩阵标识法。根据节点水位方程系数矩阵的高稀疏性，对矩阵中非零元素进行代码标识。按照代码指示，将非零元素用一维数组存储，排除零元素，节约内存。求解时，由代码指示，只对非零

7、元素进行运算，从而大大提高方程组求解计算的效率。2.河网的概化节点编码对河道中节点进行统一编码，编码顺序可以任意。但为了编程的方便，需对节点的类型进行分类：一类为普通节点，即几条河流的交汇处；另一类节点为建有水工建筑物，如闸，堰等的特殊节点。河道的编码对于外河道和内河道进行统一编码，通常编码先从内河道开始。河道流动方向对于内河道，流向是任意的，计算水流方向为正表示沿箭头所指方向，为负表示流向逆箭头方向。河道断面编码每条河道独立编码，首节点对应于首断面，末节点对应于末断面，首断面向