城市给排水管道水力瞬变现象数值计算研究.docxVIP

城市给排水管道水力瞬变现象数值计算研究

引言：

城市发生暴雨会在短时间内增加排水管道水力瞬变现象的发生。从设计角度考虑，传统管道设计标准为恒定无压流，对于2.0km2汇水面积后的排水管道，随着暴雨前后排水管的净流量增加至设计的3～4倍，设计的恒定无压流在降雨过程中，由定无压流、到非恒定无压流、再由明满流逐渐过渡到有压流。

目前针对有压流的管道高频水力瞬变分析过程中，多采用经典的拟稳态摩阻模型进行分析，此理论计算针对管壁切应力具有拟稳态性的缓慢流体瞬变进行分析是可行的，但是却无法真实的模拟高频流体的瞬变流动，无法精确的模拟分析水击压力波的畸变现象和衰减过程，因此研究并采用具有更好精度，且能够精确模拟分析高水击压力波的畸变现象和衰减过程的摩阻计算模型时当前研究的方向。基于此，文章首先对于Brunone非恒定摩阻模型与传统拟稳态摩阻模型进行了分析，然后通过有限差分方程进行了特征线方程的计算，依托于某排水管道试验段项目，进行了水力瞬变现象模拟分析，结果显示Brunone非恒定摩阻模型相比传统拟稳态摩阻模型能够更好的使用高强度的水力瞬变，且堵塞模拟分析结果显示应用Brunone非恒定摩阻模型能够较好的进行堵塞定位分析，且具有较高的定位精度。

1.管道水力瞬变及危害分析

1.1水力瞬变定义

给排水管道中的水流流量会由于降雨、或者其他原因造成流量发生较大的变化，这会在较大程度上导致管道内的压力发生一定的增减，称为水击或者水锤。由于过流液体中存在一定的气体，在高压水流中产生截留气体的示意图见图1。

根据上图可知，在尾水的填充过程中，驱动水流的压头保持一定的恒定。驱动水头由上游恒定压力流提供，截留气体位于管道末端，由于重力的作用，截留气体内不可避免存在一定深度的尾水。初始状态时，通过阀门的管壁保持左侧管道处于稳定状态，同时在右侧冲水，保持截留气体的压强处于稳定状态。在阀门打开之后，左侧水柱在上游水头压力作用下向右冲击截留气体，截留气体受到压缩后压力增大保持与上游水头压力平衡，此时截留气体的压强占据主要地位；左侧水柱的压力逐渐降低直至为零，此时截留气体的压力达到最大值，气液交界面开始向左运动直至最初的位置，如此方完成一个压力的震荡过程。

1.2伴随现象分析

由于水击现象会造成管道内压力发生较大的变化，这会诱发管道产生剧烈振动甚至发生破裂，严重会导致爆炸事故的发生。根据水力瞬变的定义，发生的情况下会产生较大的破坏作用，且在水力瞬变发生过程中伴随着下述现象：

（1）给排水管道中的水压力发生急剧的上升和下降。

（2）管道水压力发生上升和下降时，伴随着水流流动产生较强的管道水力共振；

（3）管道内水力急剧下降，会导致部分水变成水蒸气，从而导致水蒸气和液体发生分离，在液柱二次聚集时发生较大的冲击压力。

1.3破坏效应分析

由于管道内水流量的不同，且部分旧管道故障容易导致掺气现象，对于给排水管道形成一定的截留气团效应。截留气团在管道中的存在具有两面性，具体分析如下。

（1）管道水力瞬变发生的过程中，由于气体的压缩程度较大，截留气团的存在可以起到一定的缓冲作用，减小了水体因突然降速产生的水击压力，对于管道系统起到一定的保护作用。

（2）在给排水管道中，当高压水流冲击截留气体时，过高的水压产生较大的水体流通障碍，具体如下：

①截留气体的存在占据了管道一定的过流面积，造成过流能力发生一定程度的下降。

②气体和液体交界处的断面要素因为截留气体的存在引入一定的垂直动量，导致局部流量发生较大的变化，从而导致过流压力的分布也发生较大的变化。

③城市暴雨等造成管道流量急剧增大，管道内气体大量释放导致一定的井喷现象，对于管道本身和公共安全造成较大的威胁。

2.数学模型

2.1拟稳态摩阻模型

传统水力瞬变认为流体为拟稳态的流体，采用拟稳态的摩阻模型，此模型认为发生水力瞬变时的管道壁面切应力之间偏差较小，可认为是一致的。因此对于流速较低的拟稳态流体，可采用拟稳态摩阻模型对水力瞬变时刻的管壁摩阻进行近似计算，公式如下。与定长速度是的管道壁面切应力

2.2Brunone模型

相对于传统拟稳态摩阻模型，Brunone等人认为高频流体下的水力瞬变中的摩阻不仅与流体的平均速度V有直接关系，同时还与管道瞬时加速度和对流瞬时存在直接的函数关系。因此，Brunone等人提出了管道内部高频流体在水力瞬变状态下的管道内壁切应力的计算公式，具体如下所示。

D—给排水管道的内径；

3.数值计算方法

3.1特征线方程

根据特征线理论，用管内流量Q表示管内液体的平均流速，即可得到Brunone非恒定摩阻模型的特征方程，其中包括正特征值和负特征值。

3.1.1正特征值

正特征值方程计算如下式：

3.1.2负特征值

负特征值方程计算如下式：

3.2有限差分方程

根据上一节