计算水力学04.pptVIP

计算水力学 第四章 河道水流计算 天然河道水流:一维流动 基本方程:圣维南方程组 双曲线型微分方程有两类基本求解方法: 一是基于方程的特征线形式的特征线法 二是基于偏微分方程的有限差分法 介绍适用于河道水流计算的有限差分法 第一节 蛙跳格式 蛙跳方法是最先使用的方法之一 优点是接近二阶精度，质量守恒很好 B是常数时连续方程的近似解是二阶的 缺点是所得的结果为锯齿形线，原因是因变量是在奇、偶交替的点上计算的，两组网格点互不相关。 该格式的稳定性的必要条件是 其中U０为特征流速，C ０重力波速 第二节 Lax－Wendroff格式 Ｌａｘ和Wｅｎｄｒｏｆｆ提出了有限差分格式写成守恒形式的一阶微分方程系统为： 该方程首先应用于气体动力方程的系统。 应用台劳级数展开 保留表达式的头三项，并应用原方程 式中Ａ=A(Ｕ) 是Ｆ(Ｕ)关于Ｕ的雅可比行列式： 展开式变成 式中 该式是二阶精度的显式有限差分格式 方法可以简化，令Ａ为常数有: 此式可用两步法代替 1、用下式计算中间值： 2、用“校正”格式计算输出值： 格式的稳定性条件为： λ是矩阵的特征值 Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式的使用,需要初始条件尽可能满足流动方程,格式具有二阶精度,而且是不耗散的,这样,它不会使初始扰动光滑化。对于这一方程Ａ和Ｑ的边界值需用特征线法计算。 第三节 Ａｂｂｏｔｔ隐式格式 Ａｂｂｏｔｔ格式应用的微分方程形式与常用的有所不同，在连续方程中用“存储宽度 ”的概念，在动量方程中用“计算宽度Ｂ”的概念 离散格式 交错网格 差分方程采用追赶法求解。 Abbott格式是无条件稳定的，如果时步长过大时，它的收敛速度较慢。 不管怎样，这种格式可以在相当大的Courant数下计算仍然保持稳定。在潮汐波计算中， Courant数可取10~20，在明渠水流计算中， Courant数可取100。 第四节Preissmann隐式格式 由于采用显格式对时间步长ｄt要施加限制，从而促使有限差分的隐式方法得到发展。本节介绍１９６０年发展起来的Preissmann隐式方法，其实际是一个四点隐格式，对圣维南方程进行离散，得到以增量表达的非线性方程组，忽略二阶微量简化成为线性代数方程组，可以直接求解。 Preissmann格式 式中θ－－加权系数 四点线性隐格式 连续方程 连续方程 其中： 动量方程 动量方程 代入动量方程化简得 其中系数为 任一河段差分方程 为书写方便，忽略上标ｎ＋１，任一河段差分方程写成： 其中均由初值计算，所以方程组为常系数线性方程组。对一条具有Ｌ２-Ｌ１个河段的河道，有2·(L2-L1＋1）个未知变量，可以列出２· (L２-L１)个方程，加上河道两端的边界条件，形成封闭的代数方程组，可唯一求解未知量 。 格式的稳定条件和精度 追赶法求解 对于该方程组，根据不同的边界条件，可设不同的递推关系，用追赶法直接求解。对于河道的三种边界条件： 追赶法求解 １、水位边界条件的计算 设如下的追赶方程 因为 所以 追赶法求解 把Zj的表达式 代入差分方程得 追赶法求解 令 追赶法求解 追赶法求解 由 和 解得 追赶法求解 所以 追赶法求解 最后一个追赶方程 与边界条件 联立可求得 回代可求出各断面的水位和流量。 追赶法求解 ２、流量边界条件的计算 设如下的追赶方程 因为 所以 追赶法求解 把Qj的表达式 代入差分方程得 追赶法求解 以Zj+1为自由变量，解出Zj和Qj+1，得追赶系数表达式 追赶法求解 最后一个追赶方程 与边界条件 联立可求得 回代可求出各断面的水位和流量。 追赶法求解 ３、水位、流量关系边界的计算 对于水位流量关系，可线性化处理成同流量边界条件一样计算。 水位～流量关系边界 所以，可以把边界条件写成 与流量边界条件的表达式相同。 第五节 内边界的处理 在河道水流计算中，除了外部边界条件外，还可能遇到内部边界条件。 所谓内部边界条件是指河道的几何形状的不连续或水力特性的不连续点。例如，集中入流，过水断面突然放大，堰闸过流等等。 内边界处，圣维南方程组不再适用，必须根据其水力特性作特殊处理。内部边界条件通常包含两个相容条件，即流量的连续性条件和能量守恒条件(或动量守恒条件)。 一、集中旁侧入流 对于集中旁侧入流，可设一虚拟河段，这时基本的连接方程为 当上边界为水位边界条件 用上式计算虚拟河段的追赶系数，可同正常河道一样递推求解。 当上边界为流量条件 用上式计算虚拟河段的追赶系数，可同正常河道一样递推求解。 二、河道与贮水池汇合 相容方程 贮水池的连续