圣维南方程组

描述水道和其他具有自由表面的浅水体中渐变不恒定水流运动规律的偏微分方程组。由反映质量守恒律的连续方程和反映动量守恒律的运动方程组成。1871年由法国科学家A.J.C.B.de圣维南提出,故名。

一百多年来,虽然为了考虑更多的因素和实际应用方便对它的基本假定作了某些简化或改进,产生出多种不同的表达形式,但其实质没有变化。主要进展表现在求解方法的改进和创新。1877年法国工程师克莱茨提出了瞬态法。1938年苏联С.А.赫里斯季安诺维奇提出另一类解法──特征线法。但均因计算量较大,不得不进行各种简化处理,使实际应用受到限制。自50年代以来,随着电子计算机的普及,研究和提出了一整套解法,并研究出若干个通用性较强的应用软件(即程序系统),促进了圣维南方程组在水文和其他工程领域中的应用。

方程组的形式

一维单宽水流情况下,圣维南方程组的典型形式为:

式中t为时间;s 为距水道某固定断面沿流程的距离;hvZ0分别为相应于s处过水断面的水深、断面平均流速和水底高程;Hf为由于摩阻损失而引起的能量比降;g重力加速度ts为自变量;hv为因变量;Z0Hf可由shv确定。(1)式为连续方程,反映了水道中的水量平衡,即蓄量的变化率(第一项)应等于沿程流量的变化率(第二项)。 (2)式为运动方程。其中第一项反映某固定点的局地加速度,第二项反映由于流速的空间分布不均匀所引起的对流加速度。以上两项称为惯性项。第三项反映由于底坡引起的重力作用,称为重力项。第四项反映了水深的影响,称为压力项。第三、四项可合并为一项,即水面比降。第五项为水流内部及边界的摩阻损失。该式表达了重力与压力的联合作用使水流克服惯性力和摩阻引起的能量损失而获得加速度。

圣维南方程组还有许多其他形式。例如:以断面流量代替流速,以面积代替水深作为因变量;也可考虑河道两侧的沿程入流、地转力和水面风力的影响;还可把垂线平均流速作为因变量,写出二维水体渐变不恒定明流的运动方程。

基本假定

建立圣维南方程组的基本假定是:

(1)流速沿整个过水断面(一维情形)或垂线(二维情形)均匀分布,可用其平均值代替。不考虑水流垂直方向的交换和垂直加速度,从而可假设水压力呈静水压力分布,即与水深成正比;

(2)河床比降小,其倾角的正切与正弦值近似相等;

(3)水流为渐变流动,水面曲线近似水平。此外,在计算不恒定的摩阻损失Hf时,常假设可近似采用恒定流的有关公式,如曼宁公式(见河水运动)。

圣维南方程组描述的不恒定水流运动是一种浅水中的长波传播现象,通常称为动力波。因为水流运动的主要作用力是重力,属于重力波的范畴。如忽略运动方程中的惯性项和压力项,只考虑摩阻和底坡的影响,简化后方程组所描述的运动称为运动波。如只忽略惯性项的影响,所得到的波称为扩散波。运动波、扩散波及其他简化形式可以较好地近似某些情况的流动,同时简化计算便于实际应用。

求解方法

圣维南方程组在数学上属于一阶拟线性双曲型偏微分方程组。联解方程组并使其符合给定的初始条件和边界条件,就可得出不恒定水流的流速和水深(或其他因变量)随流程和时间的变化,即vv(st)和hh(st)。初始条件为某一起始时刻的水流状态,如水道沿程各断面的水深和流速。边界条件为所计算的水体的边界水流状态,如某一河段上、下游边界断面处的水位过程、流量过程或水位流量关系等。给定的初始条件和边界条件的数目和形式必须恰当,符合水流的性质,才能保证方程组的解存在和唯一,保证不致因数据的微小变化而使方程的解发生很大的变化。此时,问题称为是适定的,求解才有意义。

除特殊情况外,很难用解析方法求得圣维南方程组的解析解。一般只能通过数值计算获得个别情况的近似解。常用的数值计算方法主要有以下三类:

(1)有限差分法。将所计算的水体按照一定的网格划分,每个网格点处的微分形式的圣维南方程组,用某种形式的差分方程组来逼近。边界条件也写成差分形成。然后逐时段地求解差分方程组,得出各网格点(如断面)处的水深及流速。根据所采用的差分计算方法的不同,对每一计算时段来说,或可逐个算出各网格点处的水力要素,或是必须联立求解各网点处的水力要素。前者称为显式差分法,后者称为隐式差分法。克莱茨提出的瞬态法就属于一种简化的显式差分法。

(2)特征法。把圣维南方程组由偏微分方程组变换为在所谓“特征”上成立的常微分方程组,通常称为特征方程组。在空间为一维的情况下,“特征”的几何表示称为特征线,而在二维则为特征面。不恒定水流中的波动和干扰是沿“特征”传播的。用有限差分法联立求解表达“特征”几何位置的方程和特征方程组,即可求得所需的数值解。

(3)有限单元法。把水体划分成几何形状简单的单元(如一维的直线段,二维的矩形、直边或曲边三角形等),在每一单元内,解用数学处理比较简单的内插函数来逼近。把圣维南方程组应用于每个单元,变换为积分形式,并根据某种准则(如逼近的残差最小)来确定内插函数中的待定系数便可定解。常用的是伽辽金半离散有限单元法。

除了求解完全或简化形式的圣维南方程组的上述解法,在水文学中多年来还对一维流动发展出许多简化计算方法。例如,把运动方程简化为计算时段内计算河段的蓄水量与出流量之间关系的方程,然后联立求解。同时,已对水文学中常用的方法与求解圣维南方程组的关系进行了研究。如应用广泛的马斯金格姆(曾译“马斯京根”)流量演算法,可列为扩散波中的特殊情形。水文学方法简单,而且能较好地适用于某些情况,今后仍将长期广泛地被应用。

对于非渐变的流动,水流通过激波把两部分渐变流连接起来。如通过水跃实现由急流(超临界流)到缓流(次临界流)的过渡。在涨潮和溃坝波中也常出现近乎垂直的波前。此时,两边的渐变流仍可用圣维南方程组来描述。只要补充激波处的跳跃条件和用以判别物理上是否许可的某种准则(如熵条件等)即可求解。

圣维南方程组所描述的具有自由表面的水体的渐变不恒定流动的计算具有重要的实际意义。洪流演进计算是洪水预报、堤防设计和防洪系统运用的重要依据。水电厂引水渠、下游尾水、灌溉和通航水道中的不恒定流计算,是确定堤岸、尾水管出口高程,论证发电、航运、给水等工程设施的安全和效益的根据。此外,潮汐河口的潮流计算、溃坝决堤造成的洪水灾害的估计等也都具有十分重要的经济意义。

参考书目
  1. V.Yevjevich,K.Mahmood,ed.,Unsteady Flow in Open Channels,Vol.1,Water Resources Publ.,Fort Collins,Colorado,1975.