圣维南方程组

《圣维南方程组》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、简介　　描述水道和其他具有自由表面的浅水体中渐变不恒定水流运动规律的偏微分方程组。由反映质量守恒律的连续方程和反映动量守恒律的运动方程组成。1871年由法国科学家A.J.C.B.de圣维南提出，故名。　　一百多年来，虽然为了考虑更多的因素和实际应用方便对它的基本假定作了某些简化或改进，产生出多种不同的表达形式，但其实质没有变化。主要进展表现在求解方法的改进和创新。1877年法国工程师克莱茨提出了瞬态法。1938年苏联С.А.赫里斯季安诺维奇提出另一类解法──特征线法。但均因计算量较大，不得不进行各种简化处理,使实际应用受到限制。自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)。初始条件为某一起始时刻的水流状态，如水道沿程各断面的水深和流速。边界条件为所计算的水体的边界水流状态，如某一河段上、下

6、游边界断面处的水位过程、流量过程或水位流量关系等。给定的初始条件和边界条件的数目和形式必须恰当，符合水流的性质，才能保证方程组的解存在和唯一，保证不致因数据的微小变化而使方程的解发生很大的变化。此时，问题称为是适定的，求解才有意义。　　除特殊情况外，很难用解析方法求得圣维南方程组的解析解。一般只能通过数值计算获得个别情况的近似解。常用的数值计算方法主要有以下三类：①有限差分法。将所计算的水体按照一定的网格划分，每个网格点处的微分形式的圣维南方程组，用某种形式的差分方程组来逼近。边界条件也写成差分形成。然后逐时段地求解差分方程组，得出各网格点（

7、如断面）处的水深及流速。根据所采用的差分计算方法的不同，对每一计算时段来说，或可逐个算出各网格点处的水力要素，或是必须联立求解各网点处的水力要素。前者称为显式差分法，后者称为隐式差分法。克莱茨提出的瞬态法就属于一种简化的显式差分法。②特征法。把圣维南方程组由偏微分方程组变换为在所谓“特征”上成立的常微分方程组，通常称为特征方程组。在空间为一维的情况下，“特征”的几何表示称为特征线，而在二维则为特征面。不恒定水流中的波动和干扰是沿“特征”传播的。用有限差分法联立求解表达“特征”几何位置的方程和特征方程组，即可求得所需的数值解。③有限单元法。把水

8、体划分成几何形状简单的单元(如一维的直线段,二维的矩形、直边或曲边三角形等),在每一单元内，解用数学处理比较简单的内插函数来逼近。把圣维南方程组应用于每个单元，变换