河道崩岸、护岸、清淤、生态需水量等相关论文 河道非恒定流的数值模拟研究.pdfVIP

第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 河道非恒定流的数值模拟研究 1 2 2 2 龙吉平 张新华 李能川 严瑞平 (1. 四川大学 水电学院，成都 610065; 2. 四川大学 水力学与山区河流开发保护国家重点实验室，成都 610065) 摘要: 针对河道非恒定流控制方程简化后的运动波模型、扩散波模型和动力波模型，采用直接 差分法对三种模型进行数值离散，得出三种模型的数值模拟模型。通过数值模拟实验分析， 得出了三种波形的水深、流量的变化规律，结果表明：数值模型能够准确反映不同类型河道 中不同洪水波的波形特征，可以进行有针对性的实际应用。 关键词: 非恒定流;数值模拟;运动波;扩散波;动力波 1 引言 洪水的破坏性大，常给我们带来很大的经济损失和人员伤亡，所以对洪水的研究，了解 洪水的特性，了解洪水的运行规律，以便对洪水进行准确的预报，就显得很重要。洪水波在 河道中的运动属典型的非恒定流，河道非恒定流一般可用圣维南方程组进行描述，但是，由 于该方程组没有理论解，只能依赖数值方法求解。 河流有山区性河流和平原型河流之分，根据河流类型的不同，可将圣维南方程组进行简 化。为此，本研究针对简化后的运动波模型、扩散波模型和动力波模型，采用直接差分法对 三种模型进行数值离散，得出三种模型的数值模拟模型。 2 河道非恒定流的控制方程 对一维河道非恒定流计算的连续性方程和运动方程可以写成如下形式： ∂h ∂q + 0 （1） ∂t ∂x ∂q ∂vq ∂h gn 2 q + +gh −ghi0 + 7 3 0 （2） ∂t ∂x ∂x h 式中： x 为沿水流方向的坐标；t 为时间；h 为水深；q 为单位宽度的流量(等于vh)；v 为断 面平均流速；g 为重力加速度；n 为河道糙率系数。 - 948 - 第二十一届全国水动力学研讨会暨第八届全国水动力学学术会议 暨两岸船舶与海洋工程水动力学研讨会文集 2.1 非恒定流的三种模型 由于上述方程组没有可普遍适用的解析解，又由于这种方程的复杂性，在实际进行洪水 计算时，针对实际情况，对方程(1)和方程(2)进行简化，然后对简化方程进行数值计算。 2.1.1 运动波模型（Kinematic Wave Model） 如果将运动方程（2）中的惯性项和水深变化项的影响忽略，与连续方程（1）一起整理 则有 ∂h ∂h +V 0 （3） ∂t ∂x 式中： 5 5 1 5 3 1 2 V v ⋅ h i （4） 0 3 3 n 式（3）和式（4）就构成了常用运动波方程。 2.1.2