圣维南方程组 de Saint-Venant system of equations
中文名称 | 圣维南方程组 | 外文名称 | de Saint-Venant system of equations |
---|---|---|---|
提出时间 | 1871年 | 提出者 | 法国科学家A.J.C.B.de圣维南 |
含义 | 描述水道渐变不恒定水流运动规律 |
建立圣维南方程组的基本假定是:
①流速沿整个过水断面(一维情形)或垂线(二维情形)均匀分布,可用其平均值代替。不考虑水流垂直方向的交换和垂直加速度,从而可假设水压力呈静水压力分布,即与水深成正比;
②河床比降小,其倾角的正切与正弦值近似相等;
③水流为渐变流动,水面曲线近似水平。此外,在计算不恒定的摩阻损失Hf时,常假设可近似采用恒定流的有关公式,如曼宁公式(见河水运动)。
圣维南方程组描述的不恒定水流运动是一种浅水中的长波传播现象,通常称为动力波。因为水流运动的主要作用力是重力,属于重力波的范畴。如忽略运动方程中的惯性项和压力项,只考虑摩阻和底坡的影响,简化后方程组所描述的运动称为运动波。如只忽略惯性项的影响,所得到的波称为扩散波。运动波、扩散波及其他简化形式可以较好地近似某些情况的流动,同时简化计算便于实际应用。
V.Yevjevich,K.Mahmood,ed.,Unsteady Flow in Open Channels,Vol.1,Water Resources Publ.,Fort Collins,Colorado,1975.
圣维南方程组方程组的形式
一维单宽水流情况下,圣维南方程组的典型形式为: 式中t为时间;s 为距水道某固定断面沿流程的距离;h、v、Z0分别为相应于s处过水断面的水深、断面平均流速和水底高程;Hf为由于摩阻损失而引起的能量比降;g为重力加速度;t和s为自变量;h和v为因变量;Z0、Hf可由s、h和v确定。(1)式为连续方程,反映了水道中的水量平衡,即蓄量的变化率(第一项)应等于沿程流量的变化率(第二项)。 (2)式为运动方程。其中第一项反映某固定点的局地加速度,第二项反映由于流速的空间分布不均匀所引起的对流加速度。以上两项称为惯性项。第三项反映由于底坡引起的重力作用,称为重力项。第四项反映了水深的影响,称为压力项。第三、四项可合并为一项,即水面比降。第五项为水流内部及边界的摩阻损失。该式表达了重力与压力的联合作用使水流克服惯性力和摩阻引起的能量损失而获得加速度。
圣维南方程组还有许多其他形式。例如:以断面流量代替流速,以面积代替水深作为因变量;也可考虑河道两侧的沿程入流、地转力和水面风力的影响;还可把垂线平均流速作为因变量,写出二维水体渐变不恒定明流的运动方程。
1. 一次篮、排球比赛,共有48个队,520名运动员参加,其中篮球队每队10名,排球队每队12名,求篮、排球各有多少队参赛? 2. 某厂买进甲、乙两种材料共56吨,用去9860元。若甲种材料每吨190...
比甲金属的40分之一重7克,求两种金属各重多少克? (和差倍问题)某厂第二车间的人数比第一车间的人数的五分之四少30人.如果从第一车间调10人到第二车间,那么第二车间的人数就是第一车间的四分之三.问这...
曲线y=2x²+1,在点(1.3)处的切线方程是?求解解:y′=4x+1,故y′(1)=5,∴在点(1,3)处的切线方程为y=5(x-1)+3=6x-2.
圣维南方程组在数学上属于一阶拟线性双曲型偏微分方程组。联解方程组并使其符合给定的初始条件和边界条件,就可得出不恒定水流的流速和水深(或其他因变量)随流程和时间的变化,即v=v(s,t)和h=h(s,t)。初始条件为某一起始时刻的水流状态,如水道沿程各断面的水深和流速。边界条件为所计算的水体的边界水流状态,如某一河段上、下游边界断面处的水位过程、流量过程或水位流量关系等。给定的初始条件和边界条件的数目和形式必须恰当,符合水流的性质,才能保证方程组的解存在和唯一,保证不致因数据的微小变化而使方程的解发生很大的变化。此时,问题称为是适定的,求解才有意义。
除特殊情况外,很难用解析方法求得圣维南方程组的解析解。一般只能通过数值计算获得个别情况的近似解。常用的数值计算方法主要有以下五类 :①有限差分法。将所计算的水体按照一定的网格划分,每个网格点处的微分形式的圣维南方程组,用某种形式的差分方程组来逼近。边界条件也写成差分形成。然后逐时段地求解差分方程组,得出各网格点(如断面)处的水深及流速。根据所采用的差分计算方法的不同,对每一计算时段来说,或可逐个算出各网格点处的水力要素,或是必须联立求解各网点处的水力要素。前者称为显式差分法,后者称为隐式差分法。克莱茨提出的瞬态法就属于一种简化的显式差分法。②特征法。把圣维南方程组由偏微分方程组变换为在所谓"特征"上成立的常微分方程组,通常称为特征方程组。在空间为一维的情况下,"特征"的几何表示称为特征线,而在二维则为特征面。不恒定水流中的波动和干扰是沿"特征"传播的。用有限差分法联立求解表达"特征"几何位置的方程和特征方程组,即可求得所需的数值解。③有限单元法。把水体划分成几何形状简单的单元(如一维的直线段,二维的矩形、直边或曲边三角形等),在每一单元内,解用数学处理比较简单的内插函数来逼近。把圣维南方程组应用于每个单元,变换为积分形式,并根据某种准则(如逼近的残差最小)来确定内插函数中的待定系数便可定解。常用的是伽辽金半离散有限单元法。 ④有限元。常见的有限元计算方法有直接法、变分法、加权余量法及能量平衡法等。⑤有限分析法。在局部单元上线性化微分方程和插值近似边界的条件下,在局部单元上求微分方程的解析解,而构成整体的线性代数方程组。
除了求解完全或简化形式的圣维南方程组的上述解法,在水文学中多年来还对一维流动发展出许多简化计算方法。例如,把运动方程简化为计算时段内计算河段的蓄水量与出流量之间关系的方程,然后联立求解。同时,已对水文学中常用的方法与求解圣维南方程组的关系进行了研究。如应用广泛的马斯金格姆(曾译"马斯京根")流量演算法,可列为扩散波中的特殊情形。水文学方法简单,而且能较好地适用于某些情况,今后仍将长期广泛地被应用。
对于非渐变的流动,水流通过激波把两部分渐变流连接起来。如通过水跃实现由急流(超临界流)到缓流(次临界流)的过渡。在涨潮和溃坝波中也常出现近乎垂直的波前。此时,两边的渐变流仍可用圣维南方程组来描述。只要补充激波处的跳跃条件和用以判别物理上是否许可的某种准则(如熵条件等)即可求解。
圣维南方程组所描述的具有自由表面的水体的渐变不恒定流动的计算具有重要的实际意义。洪流演进计算是洪水预报、堤防设计和防洪系统运用的重要依据。水电厂引水渠、下游尾水、灌溉和通航水道中的不恒定流计算,是确定堤岸、尾水管出口高程,论证发电、航运、给水等工程设施的安全和效益的根据。此外,潮汐河口的潮流计算、溃坝决堤造成的洪水灾害的估计等也都具有十分重要的经济意义。
Biot固结有限元方程组的病态规律分析——针对Biot固结有限元方程组的病态问题,采用正交试验和量纲分析法探讨病态性的变化规律及其影响因素。主要考虑单元平均尺寸、时间步长、压缩性和渗透性等4个因素的影响。分别提出2组相互独立的相似准则,在此基础上得到2...
洪水入库后,其运动是属于不稳定流。水库沿程的水位、流速和过水断面等均随时变化,可用圣维南方程组表示。由于圣维南方程组一般很难求得精确的解析解,实用中多采用瞬态法,直接差分法及特征线法等近似解法。水库调洪计算的实用方法可分为静库容法及动库容法。静库容法假定库内流速趋近于零,库水面为水平的,即库容与坝前水位成单值函数关系,忽略动力方程对调洪的影响,将连续方程写成有限差的水量平衡方程并建立水库下泄流量与水库蓄水量之间的关系。动库容法指水库尾部地形开阔或回水尾端库水面曲线上翘显著,动库容(实际库水面与水平库面间的库容)不容忽视时,按不稳定流计算方法逐段进行洪水演算。这种方法的工作量很大。实际应用采用回水曲线法和近似法 。
库区水流形态受入库洪水和坝址下泄量变化的影响,属于非恒定流范畴。为此,和水库调洪计算一样,可通过圣维南方程组求解,严格推求不同时间库区沿程各断面的水位变化。为进行某一洪水标准下的水库回水计算,通常可采用入库洪水过程线为其上边界条件,采用由调度方式规定的坝前水位和泄量过程,或水位与泄量关系为其下边界条件,并取调洪开始时的入库流量与坝址泄量相等,即库区沿程处于恒定流状态下的流量及水位为初始条件,求出整个洪水过程中水库库区的流态,然后连接各断面的最高水位,即是所求洪水标准的水库回水线。
自下而上逐河段计算,即可求得整个库区的回水线。实用上,也常采用图解法或半图解法代替试算过程。中国20世纪70年代以前,常采用的这类方法有艾斯考福法和H.M.别尔纳德斯基的控制曲线法。这些方法在规划设计中应用很广。用于推求库区淹没水面线时,一般先由规定频率的调洪成果给出起始条件,并视情况拟定必要的计算条件,如最高洪水位高于正常蓄水位,通常可取下列数值作为计算条件分别推求两条水面线,然后取其外包线作为所求频率的回水线:
①坝前最高水位与相应时刻的入库流量;
②入库最大洪峰流量与相应时刻的坝前水位。
当洪水位低于正常蓄水位时,则还要以正常蓄水位与汛末相应频率的洪峰流量作为计算条件,推求一条水面线,而后再取包线。对于为分析枯水季节航运和引水条件所需的低水位回水线,一般可采用由死水位及相应于供水设计保证率的枯水期流量进行推算求得 。
水库淤积将使库区沿程过水断面积减小,引起回水上延。通常可采用水库淤积计算方法,先求出不同淤积水平(年限)的库区淤积量与分布位置,据以求得淤积后的河道断面,然后再按上述方法推求淤积后的水库回水线 。
明渠不恒定流研究始于18世纪。1788年J.L.拉格朗日(J.L.Lagrange)创立浅水波理论。1871年法国D.圣维南(De Saint Venant)根据质量守恒与动量守恒原理提出一维渐变不恒定流方程组。
明渠不恒定流的研究方法有水力模型法和数值计算法等。目前解决一维不恒定流问题多用数值计算法。圣维南方程组为一阶拟线性双曲型偏微分方程,难于直接求解。工程规划设计中一般用特征差分或直接差分等数值计算法求解。直接差分法又分为显式差分法与隐式差分法两类。显式差分法假定时段末瞬流态的沿程变率用时段初瞬值代替,流态时间变率用流段中断面数值表示,方程中非偏导数因子以时段初瞬流段中断面数值代入,因此可根据前一瞬时各点的已知量逐点分别求解下一瞬时的未知量,计算较简便。当时间步长Δt与距离步长ΔS足够小时,这种假定引起的误差是很小的。隐式差分法不作上述假定,而是联解n个差分方程,同时求出下一瞬时各点的未知量,计算工作量较大,但隐式差分法的稳定性比显式差分法好,近年来应用范围较广。当水力因素随时间和地点变化比较缓慢时,不恒定流问题(如洪水传播)可采用忽略圣维南方程组的惯性项与流速水头项的简化法,如瞬态法或简化差分法。
上列方程组的求解除需要纵横断面及糙率资料外,还需给定起始条件和边界条件。起始条件为计算初瞬河道沿程的水深(或水位)和流速(或流量)状态的数据。边界条件为河段两端断面的流量过程、水位过程或水位流量关系数值。