复式河道漫滩水流计算方法研究.docVIP

　　复式河道漫滩水流计算方法研究|第1 内容加载中... 　 摘要：对复式顺直河道和复式游荡型河道中漫滩洪水的计算方法进行了研究.其中，对复式顺直河道，从Shiono-Knight方法出发，将二次流项归并入雷诺切应力项，并将主槽和滩地的流速分布曲线直接相接，推导出一套不需联立求解系数方程组的漫滩水流流速和床面切应力分布的计算公式以及床面平均切应力的解析计算公式；对复式游荡型河道，建立了一种新的上层水流计算方法.文中的计算结果与前人的实验资料基本一致，初步表明文中的方法是合理的. 关键词：复式河道；漫滩；弯道；流速；切应力 洪水是一种常见的自然灾害.在洪水灾害的防治方面，存在许多水力学和河流动力学问题.由于受洪水灾害威胁较大的往往是河流的中下游地区，而这些区域的河道断面型式又常呈现出主槽与滩地相结合的复式形态，因此，对复式河道漫滩洪水的研究有着重要的实际意义.天然的复式河道可概括性地分为复式顺直河道和复式游荡型河道.对这二者，都有许多学者进行过大量的研究，包括模型实验、解析计算和数值模拟等.就计算方法而言，虽然目前的紊流精细模拟技术(如k-ε模型、Reynolds应力模型等各种紊流模型)已可以较详细、准确地模拟出紊流运动情况，但其计算相对复杂，因此，解析方法(或附以简单的数值方法)仍为工程设计人员所乐用.这使得许多研究人员仍不断致力于对简单的解析方法的研究，新的研究成果不断涌现. 在复式顺直河道漫滩水流计算方法方面，不同的学者提出了许多重要的方法，如SKM方法[1]、COHM方法[2，3]等.其中，SKM(Shiono-Knight方法)以水深平均的Navier-Stokes方程为基础，更准确、详细地反映了漫滩水流现象，并且，还给出了漫滩水流的解析计算公式，公式中的积分常数通过联立方程组解出.就床面平均切应力的计算而言，目前还是以经验半经验公式为主，解析公式尚比较缺乏. 对于复式游荡型河道，虽然研究成果很多[4～8]，但大多以实验为主，在计算方法方面的成果相对较少，其难度也相对较高.Shiono等[6]曾提出将复式游荡型河道以滩地高程为界划分为上层和下层.本文将采用这一划分方法，重点研究上层流动的计算问题. 1 计算复式顺直河道漫滩水流的SKM方法的近似解 SKM方法[1]给出的漫滩水流计算公式为(河道断面形态及相应的参数如图1所示)： 对于主槽或滩地 其中Ud为垂线平均流速，y为横向距离，S0为坡降，H为主槽水深，f为摩阻因子，A1,A2,A3和A4为积分常数，γ、β、α1、ω、ξ和η为中间系数，由下式确定： γ=(2/λ)1/2(f/8)1/41/H；β=Γ/ρgS0H；α1=-1/2+1/2[1+s(1+s2)1/2λ(8f)1/2]1/2；ω=gS0/(1+s2)1/2/sf/8-λ/s2(f/8)1/2；ξ=H-y-b/s；η=-Γ/(1+s2)1/2/sρ(f/8)·其中λ为无量纲涡粘性系数，Γ为二次流项，其余符号同前或见图1.计算中横断面被分成若干子区.积分常数通过在子区边缘上给出的边界条件联立解出. 图1 漫滩水流示意图 当SKM方法中的二次流项被归并入雷诺切应力项，Ud的解析公式成为： 式中下角标mc和fp分别表示主槽和滩地.取Ud(1)和Ud(3)分布曲线在y=b+hs处直接相接，即以Ud(1)分布曲线延长线近似图1中区域2分布曲线，可得： 2复式顺直河道流速分布与床面切力的计算和比较 图2～3是计算的流速分布与实验结果的比较，其中实验资料来自内容加载中... 　 3复式游荡型河道漫滩水流上层流动的基本方程 沿垂向，游荡型复式河道可分为两个子区，即上层和下层(如图5所示).从最基本的连续性方程和动量方程出发，对于上层流动，在上层水深近似均匀的情况下，压力项可省略，当水平方向的动量扩散与对流相比可以忽略时，对x(纵向)和y(横向)方向的动量方程沿水深积分，可得如下“层平均”的方程组. 式中下脚标“up”表示上层的平均值，下脚标“o”表示在上层的底边界(滩地固壁边界或滩地高程处上、下层子区的交界面边界)处的值. 连续性方程(10)可以被视作交界面处垂向流速Wo的控制方程，式(11)、(12)可用于计算水平方向的速度分量Uup和Vup.引入“表观”底部切应力τa,式(11)和(12)可写成： 沿曲线Vup/Uup=dy/dx，式(13)和(14)成为常微分方程.将其中的τa表达成与床面切应力类似的形式，则该常微分方程为： (16) 式中下脚标“lo”表示下层.式(15)和(16)便是计算上层流动的基本方程. 4 复式游荡型河道漫滩水流上层流动的计算方法 式(15)和(16)中的下层流速近似给成具有主槽的方向和断面平均流速的大小.沿求解线Vup/Uup=dy/dx，用中心差分格式离散式(15)和(16)，