水力学及河流动力学研究展望

水力学及河流动力学研究展望河流动力学的发展具有悠久的历史，但采用现代科学体系进行系统的研究则是20世纪才开始的。

河流动力学是以流体力学、地学、海洋和环境科学等为基础的交叉学科，其趋势仍是采用各学科之长，在理论探索、科学实验和数学模拟等方面深入发展。

1研究发展趋势展望河流动力学的研究，它应包含两个方面的内容，一是在传统理论张现代化量测技术的基础上，对已有的研究成果进行系统的总结、归纳和提高，对一些假定和近似处理给出更严密的论证，对一些经典的试验成果重新进行检验。

二是开拓新的研究领域和研究方向，特别要注重与其它学科和最新的科学技术融会贯通。

在上世纪的30年代至50年代，以Shields 曲线、Rouse悬沙公式、Meyer-Peter及Einstein推移质公式为代表，基本奠定了泥沙运动力学的理论体系，半个世纪以来，主要是进行补充和完善的工作，除在工程应用方面取得巨大的进展外，在理论体系上没有重大的突破。

通过数十年来的理论积蓄和量测技术的时代跨越，有望在近些年内在理论体系上取得突破性进展，在试验科学上获得重大的成果。

1.1.1基础理论研究河流动力学基础理论研究包括泥沙运动力学基本理论和河流过程原理及调整规律的研究。

早在30年代，Rouse应用扩散理论导出了悬移质泥沙浓度分布公式，即扩散方程，它是进行输沙计算的基本方程。

在现代两相流理论中，扩散模型只是宏观连续介质理论的一种简单模型。

更一般的模型是双流体模型，两相流中关于固液两相流的基本方程、作用力分析及其应力本构关系的理论，极大地促进了泥沙运动力学理论的发展。

但泥沙运动理论与固液两相流理论又有所区别，其内容更丰富，更独具创新性。

悬移质、推移质、水流挟沙力、动床阻力等等都是一般两相流理论中没有的概念。

这些概念是泥沙运动力学理论体系的基础，使得泥沙运动力学理论纰固液两相流理论更生动、更便于在生产实际中应用。

悬移质和推移质输沙理论、非平衡输沙理论、水流挟沙力、床面形态和动床阻力等都是泥沙运动力学基础理论研究的重要内容，而且在80年代以前已经发展得比较成熟，之后除了引入固液两相流的双流体模型外，并没有重大的进展，许多理论研究是低水平重复。

因此，该领域的理论研究应集中在两个方面：1）对现有的理论成果或成果或公式进行认真总结，去伪存真，归纳提高。

如钱宁（1980）关于推移质公式比较的研究堪称范例，几家著名的推移质输沙率公式尽管基于不同的理论，但都能转化为统一的结构形式，便于比较各家公式的适用范围及优缺点。

倪晋仁（1987）导出了悬移质泥沙浓度分布的统一公式，其它著名的公式都是其特例，并论证了不论从哪一种理论出发，最后的结果都与扩散理论具有相同的形式。

各公式在推导过程中都不可避免的要引入一些假设，因而理论上并不完善，适用范围也不尽相同。

关于动床阻力、挟沙力等，都已经取得不少的成果，也应该进行类似的归纳总结工作。

2）对不成熟的理论进行深入研究，争取取得理论上的突破。

这些方面包括：非均匀沙不平衡输沙理论、高含沙水流运动理论、床面形态的空间结构及动术阻力、管道输送固体物料的减阻机理、水流相干结构对泥沙输移的影响等等。

河流过程原理主要是指河流的自动调整原理。

“水往低处流”即形象地揭示了河流发源于高山峻岭，奔流到海不复返的自然规律。

来水来沙作用于不同的边界条件，形成了丰富多彩的河道演变现象。

河床演变学就是研究河道演变过程，它不仅仅停留在对现象的描述，而且更重要的是探讨控制河道演变的规律。

如不同河型的形成、演变及转化条件，河流的自动调整原理等。

在固壁边界条件下，水流泥沙运动参量可以通过动力学方程求解得出。

但对于不断蜿蜒展宽（或缩窄）的松散边界的冲积河流来说，还缺少一个能反映河流（横向）调整规律的动力学方程。

80年代以来，以杨志达（C.T.Yang）为代表的一批学者提出“能量耗散率极值”的条件（见Yang et al,1996），建立补充方程来封闭动力学方程组，取得了明显的进展，成为河流过程原理研究的重要方向。

此外，自80年代以来，黄河频繁断流，河道断流引起河道萎缩，加重了黄河下游洪水灾害的危险。

断流条件下黄河下游河道演变规律亦是一个全新课题，是值得高度重视的研究方向。

1.1.2不平衡输沙和非恒定流输沙1.1.2.1非恒定流输沙对于一条天然的冲积河流，在恒定水流的作用下，其河床的冲淤变化总是趋向于平衡，但在非恒定流的作用下，冲刷或淤积的变化可能向单一的方向发展而造成灾害。

河道的冲淤变化不仅取决于水流能量的大小，而且与其能量的变化率有直接的关系，河床的剧烈变化一般都是在洪水陡涨陡落的过程中发生的，这也是边岸坍塌甚至溃决的最危险的时期。

王兆印（1998）认为“非恒定流中的挟沙力、沙波运动和河床演变都有其特有的规律，需要专门研究”宋天成和Graf（1996）的文章“明渠非恒定流的流速和紊动分布”因其“在水流研究中具有卓越的价值”而在第27届国际水力学大会上被美国土木工程师协会（ASCE）评为1997年的Hilgard水力学奖（每两年从全世界水利类的学术论文中评选一篇）。

颁奖公告认为“在洪水（非恒定流）条件下的泥沙输移可能带来灾害性的后果和对水利工程（如大坝和水库）以及环境的实际损坏。

迄今为止，人们主要进行均匀流的研究。

现在，量测仪器和数据采集系统的进展使得非恒定挟沙水流的研究成为可能。

该文的研究奖有望开创一个新的研究领域。

”迄今在恒定均匀流的研究方面已取得了丰硕的成果，在清水非恒定流的研究方面也有较大的进展（Nezu,1997），而非恒定挟沙水流的研究则处于刚起步阶段，代表性成果见于-德非恒定流输沙研究成果论文集（IJSR,1997,1997,2001）.在非恒定流的条件下，泥沙输移一定是不平衡的，即不平衡输沙是该课题研究的核心。

如Cellino 和Graf (1990)的水槽试验结果表明，在饱和与非饱和条件下，泥沙的输移规律是不相同的。

1.1.2.2不平衡输沙窦国仁（1963）最早提出了在矩形均匀断面条件下的不平衡输沙公式：()*vv L V U a S S h ςςωααχ∂∂+=-- (9-1)式中：α为泥沙恢复饱和系数；*S 为垂线平均的水流挟沙能力。

韩其为等（Han 1980）将方程（9-1）进一步扩展应用于天然河道，在恒定流的条件下将上式沿垂线积分，并采用在床面的泥沙扩散率和沉降率为零的条件得出()*v V k ds aS a S dx qω=-- (9-2) 式中：α为底部含沙浓度与断面平均含沙浓度的比值，K α为底部饱和含沙浓度与断面平均含沙浓度的比值，若近似认为K αα=，即（9-1）式中的恢复饱和系数，将上式改写成：()()***V V d S S ds a S S dx q dxω-=--- (9-3) 对上式积分可得：()*0*01a La L q q V q S S S S e e a L ωωω⎛⎫=+-+- ⎪ ⎪⎝⎭(9-4) 上式即为恒定流动中平均含沙浓度沿程的变化，出口断面的含沙浓度取决于进口断面的含沙浓渡0S 、进口断面的饱和含沙浓*0S 、出口断面的挟沙能力*S 、河段长度Ｌ及恢复饱和系数α等参数。

韩其为(1997)对非均匀沙的二维不平衡输沙方程及边界条件进行了深入的研究，较严密地推导了恢复饱和系数的表达式，能较好地概括已有的研究成果。

周建军（Zhou1990，1997）在假定的不平衡垂线浓度分布剖面条件下，得到了不平衡输沙方程和恢复饱和系数的近似计算公式，采取侧向积分的方法推导了适用于天然河道总流的不平衡输沙方程。

研究结果表明，在二维数学模型和一维数学模型计算中，应采用不同的恢复饱各系数。

1.1.3 颗粒流研究研究河流动力学的理想方法应是分别写出两相各自的控制方程和和建立两相之间的本构关系，从数学上求解方程组，以获得对两相流运动的完整描述。

动理学的方法为这方面的研究提供了新的思路（傅旭东，王光廉2002）。

1.1.3.1动理学的理论基础固液两相流动问题在自然界和工程应用中广泛存在，相应的研究方法多种多样。

这些方法基本可分为宏观描述的连续介质方法和微观描述的动理方法。

其中，连续介质方法由于在流体力学中的成功而应用较早。

在近二三十年内，随着在气固两相流和快速颗粒的研究中取得长足进展，动理学方法在固液两相流的研究中也有一定的应用，如Wang 和Ni(1991)、Aragon(1995)等的研究。

基于微观单颗粒分析的动理学方法，不仅能够提供单颗粒尺度上的微观信息，还可以通过对颗粒运动信息的统计平均，导出颗粒相连续介质形式的守恒型方程，并提供相应的宏观输运系数。

从描述颗粒相运动的Boltzmann 方程出发，在一定的流动条件下求解出均、弹性、无粘性的球形颗粒的速度分布函数f ，那么该条件下的颗粒相运动的宏观特征参量也就随之确定下来。

若单颗粒的质量为m ，则单位几何体积内的颗粒数目n 为：i n fdV =⎰ (9-5) 颗粒相的平均速度i U 为：1i i i U U fdV n =⎰(9-6) 颗粒相的何种比浓度C 、颗粒相的分密度s P 和颗粒相脉动速度i U 分别为：s C nV = 、s P nm =、i i i u U U =- (9-7)其中：s V 为单颗粒休积。

在经典的气体分子动理学理论中，通常将颗粒速度相空间内的随机速度坐标i U 变换为脉动速度坐标i u ，并有下面形式的Boltzman 方程：()i i i i i ci i i i i F f dU U df f f f f u u dt x dt u u U x t ∂∂∂∂∂∂⎛⎫+-+-= ⎪∂∂∂∂∂∂⎝⎭ 9-8 定义颗粒属性ϕ的平均量1in fdV ϕϕ=⎰，并将ϕ乘以方程（9-8）的两边，在整个脉动速度空间内对方程（9-8）积分，就有颗粒属性ϕ的输运方程：()()i i i i i i i i U dU d d n n n u n u dt x x dt x dt u ϕϕϕϕϕϕ⎧⎫∂∂∂∂⎪⎪++-+-⎨⎬∂∂∂∂⎪⎪⎩⎭ i i i i c i j i U f n F u du u u x t ϕϕϕ⎧⎫∂∂∂∂⎪⎪⎛⎫--=⎨⎬ ⎪∂∂∂∂⎝⎭⎪⎪⎩⎭⎰ 9-9 在颗粒碰撞弹性、无摩擦的假定下，粒间碰撞并不改变颗粒相的数量、动量和能量。

分别令m ϕ=，i mu ϕ=,/2i i mu u ϕ=,由输运方程就得到颗粒相的守恒方程如下：连续方程 ： 0s i s idp V dt x ρ∂+=∂ （9-10） 动量方程： ji i ss i j p dV F dt x ρρ∂=-∂ （9-11） 脉动能方程： 32i i s s i i ij j iV q dT Fv dt x x ρρρ∂∂=--∂∂ （9-12） 式中，ij s i j P v v ρ=为颗粒相的脉动应力张量、212i s s i q v v ρ=为颗粒相的脉动能传导通量、13i iT v v =为颗粒相脉动能。