圣维南方程组求解.docxVIP

\o查看水动力学中的全部文章水动力学

圣维南方程组的求解

发布时间：2012年06月23日

分类：\o查看水动力学中的全部文章水动力学

自从Stoker（1953）首次尝试将完整的Saint-Venant方程组用于Ohio河流的洪水计算以来，出现了大量的针对完整的Saint-Venant方程组的数学模型（动力波模型）。求解圣维南方程组的数值方法很多，按离散的基本原理可分为特征线法、有限差分法、有限元法、有限体积法和有限分析法等。

有限差分法

显式方法的先驱是Stoker（1953），其后有Liggett和Woolhiser（1967），Martin和DeFazio(l969)及Strelkoff（1970）等人，Dronkers（1969），Balloffet（1969）及Johnson（1974）将显式方法用于分析河口的潮汐运动，Garrison等(1969)，Johnson(1974)将显式方法用于模拟河道及水库的洪水，Liggett和Cunge(l975)给出了数种显式差分格式的表达式及分析结果。对于每一计算时刻，关于计算断面的未知量，显式方法可直接从代数方程组中得出结果。

隐式方法的提出是出于显式方法由于计算的稳定性要求而存在时间步长限制的考虑。隐式方法首先是由Isaacson等(1953)建议的，其后在六、七十年代很多学者在隐式方法进行了大量的研究工作。隐式方法则要求解代数方程组。代数方程组又分为线性和非线性两种，前者既可用直接法又可用迭代法求解，而后者要用迭代法求解。在迭代法中，Newton-Raphson方法以其收敛速度快的特点而较为普遍地用于求解非线性代数方程组中，该方法首先由Amein和Fang(1970)应用于Saint-Venant方程组的数值解中，其后，国内外都有将这一方法用于河网水力数值模拟中。

直接法是用差商代替导数，将微分方程化为代数方程组，再求出区域网点上的解；而特征线法是首先将质量和动量方程进行等价变换，化为由四个常徽分方程构成的方程组，再用有限差分近似来求解。特征线法既可以是矩形网格，也可以是曲线网格，见Abbott(1966)；既有显式格式，如Liggett和Woolhiser(1967)，也有隐式格式，如Amein（1966）。

对于各种隐式差分格式的数值稳定性和精度问题，Cunge(1966)，Abbott和Ionescu(1967)、，Dronkers（l969），Gunaratnam和Perkins(1970)，Fread（1974），Liggett和Cunge(1975)以及Ponce和Simons（1977）都进行过研究。从线性化稳定性分析的角度，各种隐式差分格式都是无条件线性稳定的，即线性化的Saint-Venant方程组的隐式差分格式稳定性是与时间步长和空间步长无关的。然而Chaudhry和Contractor（1973），Fread(1974)以及Cunge(l975)发现，对于瞬间变化迅速的水流的模拟，若时间步长太大，则会出现数值不稳定；对于由截面在纵向及垂向变化迅速而引起的非线性，也会出现数值不稳定。同时时间步长还受精度、波形、Courant条件、空间步长及隐式差分格式类型的限制。

对Saint-Venant方程组的隐式差分格式，其线性方程组的求解技术是极为重要的。较为有效的是Fread(1971)的关于五对角元的压缩存贮消元法方法及Liggett和Cunge(1975)的双追赶方法。这两种方法在目前的河网水力数值模拟中使用较为普遍。然而，这两种方法均存在着绝对值较小的数作除数而引起计算中断或数值不稳定的隐患。

在众多离散格式中，Preissmann(1961)加权四点隐式格式是使用较为普遍的一种，该格式具有适用于非均匀空间步长、计算的稳定性和收敛性易于更改及处理边界条件较为简单的优点。另外还有六点格式风，但六点格式存在等空间变长网格的限制，且处理边界条件时要比四点格式复杂得多，因此使用并不普遍。

有限元法

对完整Saint-Venant方程组的求解方法还有有限元法，如Cooley和Moin(1976)，Gray等（l977）。

有限元法借鉴于固体力学，在六十年代用于流体力学中，其理论基础是极值原理和剖分插值，根据有限差分中离散处理思想，同时采用了变分计算中选择逼近函数及对任意形状（三角形或四边形）的许多微小单元进行积分处理的合理方法，因而具有很强的适应性，特别适合于几何、物理条件比较复杂的问题。常见的有限元计算方法有直接法、变分法、加权余量法及能量平衡法等。有限元法适应性强，计算精度较高，但存在计算格式复杂、计算量较大，大型系数矩阵求解困难等缺点，而且在误差估计、收敛性和稳定性等方面的理论研究与有限差分法相比还显得不够