圣维南方程组

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

2、领域中的应用。方程组的形式一维单宽水流情况下，圣维南方程组的典型形式为：式中t为时间；s为距水道某固定断面沿流程的距离；h、v、Z0分别为相应于s处 过水断面的水深、断面平均流速和水底高程;Hf为由于摩阻损失而引起的能量比降；g 为重力加速度；t和s为自变量；h和v为因变量；Z0、Hf可由s、h和v确定。（1）式 为连续方程，反映了水道中的水量平衡，即蓄量的变化率（第一项）应等于沿程流量的变 化率（第二项）。（2）式为运动方程。其中第一项反映某固定点的局地加速度，第二项反 映由于流速的空间分布不均匀所引起的对流加速度。以上两项称为惯性项。第三项反映 由于底坡引起的重力作用，称为重力项。第四项反

3、映了水深的影响，称为压力项。第 三、四项可合并为一项，即水面比降。第五项为水流内部及边界的摩阻损失。该式表达 了重力与压力的联合作用使水流克服惯性力和摩阻引起的能量损失而获得加速度。圣维南方程组还有许多其他形式。例如：以断面流量代替流速，以面积代替水深作 为因变量；也可考虑河道两侧的沿程入流、地转力和水面风力的影响；还可把垂线平均 流速作为因变量，写出二维水体渐变不恒定明流的运动方程。基本假定建立圣维南方程组的基本假定是：流速沿整个过水断面（一维情形）或垂线（二维情形）均匀分布，可用其平均值 代替。不考虑水流垂直方向的交换和垂直加速度，从而可假设水压力呈静水压力分布，即 与水深成正比；河床比降

4、小，其倾角的正切与正弦值近似相等；水流为渐变流动，水面曲线近似水平。此外，在计算不恒定的摩阻损失Hf时， 常假设可近似采用恒定流的有关公式，如曼宁公式（见河水运动）。圣维南方程组描述的不恒定水流运动是一种浅水中的长波传播现象，通常称为动力 波。因为水流运动的主要作用力是重力，属于重力波的范畴。如忽略运动方程中的惯性 项和压力项，只考虑摩阻和底坡的影响，简化后方程组所描述的运动称为运动波。如只 忽略惯性项的影响，所得到的波称为扩散波。运动波、扩散波及其他简化形式可以较好 地近似某些情况的流动，同时简化计算便于实际应用。求解方法圣维南方程组在数学上属于一阶拟线性双曲型偏微分方程组。联解方程组并使其

5、符 合给定的初始条件和边界条件，就可得出不恒定水流的流速和水深（或其他因变量）随流 程和时间的变化，即v = v（s, t）和h = h（s, t）o初始条件为某一起始时刻的水流状 态，如水道沿程各断面的水深和流速。边界条件为所计算的水体的边界水流状态，如某 一河段上、下游边界断面处的水位过程、流量过程或水位流量关系等。给定的初始条件 和边界条件的数目和形式必须恰当，符合水流的性质，才能保证方程组的解存在和唯 一，保证不致因数据的微小变化而使方程的解发生很大的变化。此时，问题称为是适定 的，求解才有意义。除特殊情况外，很难用解析方法求得圣维南方程组的解析解。一般只能通过数值计 算获得个别情况的

6、近似解。常用的数值计算方法主要有以下三类：有限差分法。将所 计算的水体按照一定的网格划分，每个网格点处的微分形式的圣维南方程组，用某种形 式的差分方程组来逼近。边界条件也写成差分形成。然后逐时段地求解差分方程组，得 出各网格点（如断面）处的水深及流速。根据所采用的差分计算方法的不同，对每一计 算时段来说，或可逐个算出各网格点处的水力要素，或是必须联立求解各网点处的水力 要素。前者称为显式差分法，后者称为隐式差分法。克莱茨提出的瞬态法就属于一种简 化的显式差分法。特征法。把圣维南方程组由偏微分方程组变换为在所谓“特征”上 成立的常微分方程组，通常称为特征方程组。在空间为一维的情况下，“特征”的几

7、何 表示称为特征线，而在二维则为特征面。不恒定水流中的波动和干扰是沿“特征”传播 的。用有限差分法联立求解表达“特征”儿何位置的方程和特征方程组，即可求得所需 的数值解。有限单元法。把水体划分成几何形状简单的单元（如一维的直线段，二维的 矩形、直边或曲边三角形等），在每一单元内，解用数学处理比较简单的内插函数来逼 近。把圣维南方程组应用于每个单元，变换为积分形式，并根据某种准则（如逼近的残 差最小）来确定内插函数中的待定系数便可定解。常用的是伽辽金半离散有限单元法。除了求解完全或简化形式的圣维南方程组的上述解法，在水文学中多年来还对一维 流动发展出许多简化计算方法。例如，把运动方程简化为计算时

8、段内计算河段的蓄水量 与出流量之间关系的方程，然后联立求解。同时，已对水文学中常用的方法与求解圣维 南方程组的关系进行了研究。如应用广泛的马斯金格姆（曾译“马斯京根”），可列为 扩散波中的特殊情形。水文学方法简单，而且能较好地适用于某些情况，今后仍将长期 广泛地被应用。对于非渐变的流动，水流通过激波把两部分渐变流连接起来。如通过水跃实现由急 流（超临界流）到缓流（次临界流）的过渡。在涨潮和溃坝波中也常出现近乎垂直的波 前。此时，两边的渐变流仍可用圣维南方程组来描述。只要补充激波处的跳跃条件和用 以判别物理上是否许可的某种准则（如熠条件等）即可求解。圣维南方程组所描述的具有自由表面的水体的渐变不恒定流动的计算具有重要的实 际意义。洪流演进计算是洪水预报、堤防