水电站非恒定流计算-课程及设计指导书.docx

水电站非恒定流计算课程设计指导书简新平河北工程大学水电学院前言水电站非恒定流，包括动力渠道涌波、调压室水位波动、压力管道水击等过程。这些过程都属于水力非恒定流，是水电站设计中水力计算的难点，也是相应的水电站引水建筑物设计的控制工况。课程设计的目的，是通过锻炼使同学们把握这些水电站非恒定流计算的方法，结合计算机的应用，锻炼解决比较复杂的计算问题的能力。本指导书主要讲述动力渠道涌波、调压室水位波动、压力管道水击等过程的计算方法。希望同学们在课程设计工程中，充分理解这些方法，并在此基础上灵活应用，利用计算机编程或者Excel表格计算。最后，要求同学们不能在课程设计说明书中照抄本指导书中的内容，课程设计说明书应该根据自己的计算过程，独立编写。扩散差分求解明渠一维非恒定流圣维南（Saint-Venant）方程组规范中的计算公式运动方程：（1）连续方程：（2）以上为《水电站引水渠道及前池设计规范》（SL/T205—97）附录公式。式中，A——横断面积，Q——流量，v——平均流速，h——水深，m——渠底纵坡——摩擦坡度t——时间x——沿渠底度量的距离向下游为正g——重力加速度q——横向进流量，入流为正，出流为负，因次为——横向进流流速沿下游方向的分量对于求解的水电站引水渠道中的涌波，属于弱解，其差分格式应满足相容性、收敛性、稳定性及幅度耗散性。计算的初始条件为渠道恒定流时的流速和水深。上游边界条件，一般假定上游水位为常数，对于自动调节渠道是适宜的；对非自动调节渠道（通常设有侧堰）或有调节池布置的情况时，宜按实际情况建立其上游边界条件。下游边界条件一般为出流量变化条件，此时忽略压力管道中的水弹性现象，假定机组过流量的变化就是前池出流量的变化。公式的变化对梯形断面，，。此时，B为水面宽度。对矩形断面，，。B为水面宽度，也是矩形断面宽度。事实上，对于非棱柱体，都有，B为水面宽度。当时，考虑到（摩阻坡度），用s代替x表示沿渠底向下游的长度，上述方程可简化为：（3）（4）对于（3）式，结合上图理解如下：由，得到。当地加速度(，即时变加速度)是液体在重力的分力（）、摩擦阻力（）、水压力差（）这三个力的合力下产生的总加速减去位变加速度(，即迁移加速度)后的加速度值。非恒定流的本质是，非恒定流在断面上的加速度随时间变化。恒定流不是没有加速度，而是流场内各点的加速度分布不随时间发生变化。当地加速度为0、迁移加速度也为0时，为（恒定）均匀流。从欧拉视角，如果停留在某一断面观察，则是某一断面上速度随时间的变化率。可以这样理解，当时间增加时，速度的增量为，即为当时间变化时液体沿长度方向移动，而不同s位置的速度不同（非恒定流或非均匀流时水深沿程变化，流量和动能沿程变化。均匀流时水深、流量和动能沿程不变化。）流速沿程变化。部分加速度导致液体之间交换能量，导致液体自身的势能和动能相互转化。当水流为恒定流时，任一断面上，速度不随时间变化，。但是加速度的部分可以不为0，考虑到，此时得到恒定渐变流的方程：上式说明，恒定非均匀流时，势能沿程的变化，一方面转化为动能，一方面用来克服摩擦阻力。加上流量沿程不变的连续性方程，就可以计算恒定非均匀流水面线。对于（4）式，可以利用质量守恒来理解，进入水体和流出水体的流量不相等，会导致水水深变化。扩散差分法几种差商导数定义，则差商是导数的近似，而导数是差商的极限。一阶差商可有向前、向后和中心差商三种情况。一阶向前偏差商在B点处对s的一阶偏差商：一阶向后偏差商在B点处对s的一阶偏差商：一阶中心偏差商在B点处对s的一阶偏差商：一阶中心偏差商，也可以看成是前后两偏差商的平均值。显然，用中心偏差商来代替导数，较之用前后偏差商的精度为高。差分格式收敛的必要条件是库朗（Courant）条件[]：常用的显示差分格式有扩散格式、蛙跳格式等。扩散格式在时间上取带权的差商代替微商，在空间上取中心差商。蛙跳格式在时间和空间上均取中心差商。扩散差分格式用G代表函数，在节点上取微商：域内节点采用（5）（6）式中称为权因子，它在[0,1]区间内取值。显然，取值不同式（7）右端值不同。实践表明，时，结果比较光滑；当时，格式不稳定；当时为拉克斯格式（Lax格式，纯扩散格式，稳定格式）。经验表明，取较小值时，如，常可得较好的成果。当时，当时，如果在每个瞬时计算流段中，函数成线性变化，无论取什么样的值，式（7）右端则不变。显而易见，若函数呈凹或凸曲线变化，则或。因此，说明函数呈非线性变化时（河渠中非恒定流绝大多数如此），的偏差商随的取值而改变，而且随着时间的推移不断地发展下去。由此可见，的取值是否符合实际的水流情况这是数值计算中影响稳定性的原因。域上边界节点为（7）（8）域下边界节点为（9）（10）方程（5）、（6）中的变系数及非导数项、、等，用、、、计算，或者用、、… 等计