水面线误差论文.docVIP

分段计算的渠道恒定非均匀流能量方程稳定性探讨及验证 摘要：工程设计中求解渠道水面线常采用能量方程分段计算，但分段计算将导致计算误差向后续计算传播，因此需对能量方程的稳定性进行分析，以确保计算误差能在后续计算的传播过程中得到有效控制。本文对渠道恒定非均匀流能量方程的稳定性进行了初步探讨，并给出误差控制的建议，同时对天然河道水面线计算也有一定参考意义，可供工程设计人员参考。 关键词：恒定非均匀流；能量方程；稳定性； 1 导言 工程设计中常需求解渠道水面线，渠道水面线常采用能量方程分段计算。分段计算将会导致计算误差（譬如数字舍入）产生的小扰动向后续计算传播。如果不对计算误差扰动的传播进行有效控制，这些计算误差可能会在传播过程中恶性增长，甚至于能“淹没”能量方程的“真解”。因此对能量方程进行稳定性分析是十分必要的。 2 能量方程 渠道恒定非均匀渐变流的微分方程【1】为： （1） 式中： —渠道水位； —动能修正系数； —渠道流速； —沿程水头损失； —局部水头损失； g一重力加速度； 将（1）式改写为差分方程为： （2） 式中： —计算分段上游断面的水位； —计算分段下游断面的水位； —计算分段上游断面的平均流速； —计算分段下游断面的平均流速； 工程中常令，（2）式经变形可改写为： （3） 式中： ； —沿程水头损失，与计算分段上、下游的水位有关； —局部水头损失系数； 3 误差传播方程 设能量方程（3）式在(，)处达到平衡，即满足。此时如有一个扰动进入能量方程，使能量方程在z0+处产生新的平衡，将满足下式： （4） 式中： ε上—上游断面的水位误差扰动值； ε下—下游断面的水位误差扰动值； 将（4）式采用一阶泰勒公式在、处展开，并略去高阶微量，有： （5） 从（5）式中消去，即有： （6） 整理即可得误差传播方程： （7） 或 （8） 4 误差传播方程分析 我们先对（7）或（8）式中和的值域进行讨论。根据的定义，令，有: == （9） 式中： —通过计算段的流量； —计算断面过流面积； —计算断面的过流水面宽度； 显然，的值域为。 水力学计算中，常近似按下式计算【1】： （10） 式中： —计算分段的长度； —计算分段流量模数平方的平均值，工程计算中常采用以下三种公式之一计算[1]： （11） （12） （13） 式中： —上游断面的流量模数； —下游断面的流量模数； 对于（11）式，有。由于为单调增函数，且必有K0，于是必有；同理必有。采用同样的方法可分别验证，对于（12）式和（13）式，也必有和。 当从下游向上游推算水面线时，要使算法稳定需满足条件。在渠道水力条件变化不过分剧烈的情况下，根据（7）式中各分项的正负关系，可知当0时，容易得到满足，即水面线宜从下游向上游方向推算。当从上游向下游推算水面线时，要使算法稳定需满足条件。在渠道水力条件变化不过分剧烈的情况下，根据（8）式中各分项的正负关系，可知当0时，容易得到满足，即水面线宜从上游向下游方向推算。如忽略局部水头损失，即取时，此时0为缓流，0为急流。即缓流时渠道水面线推算宜从下游向上游方向推算、急流时宜从上游向下游方向推算。 5 能量方程稳定性验证算例 6、结语 能量方程分段求解渠道水面线是工程常用的方法。分段计算过程中产生的误差也将会逐段传播，工程计算中应注意对误差传播进行控制，防止误差传播“淹没”真解。水面线计算时，选择从下游向上游推算还是从上游向下游推算，对于误差的传播与控制，有时具有决定意义，因此在选择计算方向的时候，尤其需要审慎对待。本文虽然仅仅只针对渠道水面线进行了探讨，但对于天然河道水面线计算也具有参考价值，在进行天然河道水面线计算时，可用差分替代微分采用（7）式或（8）式对误差传播进行预估和控制。 参考文献： [1] 吴持恭.水力学[M].北京：高等教育出版社，2008。 [2] 李庆扬，王能超，易大义等.数值分析[M].武汉：华中科技大