水电站过渡过程与仿真3.pptVIP

水电站过渡过程与仿真 武汉大学 程远楚 第三章 瞬变流的特征线计算方法 特征线方法是目前求解管道系统水力瞬变最常用的数值计算方法。特征线方法具有许多优点:①稳定性准则可以建立;②边界条件很容易编成程序，可以处理很复杂的系统;③可以适用于各种管道水力瞬变分析，包括汽液两相瞬变流;④在所有差分法中具有最好的精度。下面介绍特征线方法的计算原理(Wylie和streeter)。 求解基本方程式，必须有初始条件和边界条件。在水力过渡过程的分析中，我们总是以恒定流状态作为初始条件。至于边界条件，则因管道系统而异。本章还将介绍几种典型的边界条件。 3.1 特征线方程 连续方程和动量方程和组成了一对准线性双曲型微分方程。其中有两个因变量即速度和水力坡度线高度，两个自变量是沿管距离和时间。这两个方程可以用特征线方法变换成四个常微分方程。 水力瞬变的运动方程和连续方程可以写成下述形式: 将上述方程用一个未知因子几进行线性组合得 3.1 特征线方程 根据微分法则 且令 解得： 得： 将?的解代入上两式，有： 分别称为C+上成立的相容性方程和C+特征线方程 3.1 特征线方程 且令 解得： 得： 将?的解代入上两式，有： 分别称为C+上成立的相容性方程和C+特征线方程 分别称为C+上成立的相容性方程和C+特征线方程 3.1 特征线方程 在一般情况下，a》V，故可在特征线方程中略去V。另外，相容性方程中Vsin?/a项也可以忽略不计。 3.1 特征线方程 对于给定管道，a通常是常数，于是在图8所示的x-t平面上两个特征线方程画出来是两根直线。这些在x-t平面上的直线称为特征线。沿这些特征线，相容性方程成立。 3.2 有限差分方程 把一根管子等分为n段，每一段的长度为?x，如图3-2所示，算出时步为? t=?x/a。 在图示网格中，正向倾斜的对角线AP满足方程C+特征线方程。如果在A点的因变量V和H已知，则C+相容性方程(此方程沿C+线成立)可以在端点A和P间积分，从而可以用P点的未知量V和H把该方程表示出来； 图上朝负方向倾斜的对角线BP满足方程C-特征线方程，沿BP对C-相容性方程积分，用B点的已知条件和P点的未知条件，可以得出用在P点的同样两个未知变量表示的第二个方程。 联立求解这两个方程，可以得出xt平面中以P标出的点在特定时间和位置的参数。 3.2 有限差分方程 以adt=dx乘相容性方程，并引入管道流量Q＝AV，然后沿C+特征线积分得 即： 若假设上式积分号里被积函数在A点可导，则 当忽略上式中?x的二阶以上微量 3.2 有限差分方程 当忽略上式中?x的三阶以上微量，并取 整理后有： 3.2 有限差分方程 上两式可简化为： 求解有： 3.3 基本边界条件 截面i是x方向的任一网格交点，称为管内计算截面。 解液体瞬变流问题时，通常从时间为零时的定常状态开始。因此，每一个计算截面上的H及Q的起始值( t=0时)是已知的，解算首先就是沿着?t求每个网格点的H及Q，然后接着在t= 2?t上计算，依次类推，一直计算到所要的时间为止。 显然，对于管内计算截面，在任何一个时刻，它在前一时步的数值总是已知的，当管路两端的边界条件给定时，时刻t计算断面i的未知流量Qpi和水头Hpi可以联立C+和C-相容性方程得到。 单管的任何一端，只有一个两变量的相容性方程可用。对上游端，方程沿C-特征线成立，而对下游边界，方程沿C+特征线成立。这是两个关于Qp和Hp的线性方程，每一个方程把瞬变期间的管内流体的整个特性和响应传到相应的边界上去。这种情况，都要一个辅助方程来规定Qp和Hp，或者规定它们之间的某种关系。那就是说辅助方程应将边界的情况传给管子。求解时，每个边界条件和另一个边界条件无关，和内部点的计算也无关。 3.3 基本边界条件 1. 上游端为已知水位的水库: 对上游为大水库的情况，在很短的瞬变期间，该处的水力坡度线高度通常可以假定不变。这个边界条件可写成: 如果水库水面以一个已知的规律来变化，譬如说，正弦波，那末边界条件是: 2. 上游端流量是时间的已知函数：由排液泵输送的流量可以表示为时间的显函数，例 3.3 基本边界条件 3. 管道下游为盲端: 水轮机导叶或水泵出口阀门全关闭时，它们上游管道的末端就属于这种情况。在这样的条件下，盲