在水文学中,圣维南方程具有重要的意义,主要是两个方程:连续方程与运动方程。连续方程比较简单,动力方程相对比较复杂。不过,就是连续方程,在不同的水文学教材中,讲解的方式也有差异,但在直观理解上似乎还有改进之处。
我经过学习与比较,觉得可以用经过改写的式子,虽不很严谨、但很直观的方法帮助理解。我对连续性方程用“水涨面高、水落面绌”(洪水波的波前期,即涨水段,水位上升,过水断面水面抬升,面积增加;洪水波的波后期,即落水段,水位下降,过水断面水面下载,面积减少)来表达,公式仅用一行式子就可以直观解释(见下图左上角)。现在的连续方程,形状上就有点怪异。我是改用了等价、直观的表达。
附1:我与一位水文学同行的对话:
陈:水文学中,圣维南方程是理论上的基石。我对一些教材中的阐解总是感到不容易理解。因为我教学多以土地课程为主,涉及不到给学生讲解。有哪位老师也感到难以向学生直观解释清楚圣维南方程的连续性方程的?或者感到很容易理解的,不妨介绍一下经验。
彭:连续性方程还比较好解释。可以先假定河道宽度是一个定值,即等河宽河道,那么所研究河段的河道水位增加,是河道下断面泄水能力比上断面来水能力低造成的,即上下断面流量差造成。在某个时间段内,这样造成的来水总量与泄水总量差值,就积累在这个河道里了,造成水位增长,也就是造成这个河段总蓄水量增加。这样一个差分方程,进行微分化,就是连续性微分方程了。连续性方程的本质,就是水量不可能在这个河段内自己增值,必须是从河段的外部进入,然后又出到外部。即,研究河段内没有源汇项。
陈:谢谢彭老师的专业点评!这两本书都对连续性方程做了详细的图文解释,我经过学习与思考,终于找到了简单的理解方法。我将把其放入科学网博客,称其为最简单、最直观理解法。
彭:圣维南方程,学生不太好理解的是运动方程,它需要理解本地加速度与迁移加速度。
陈:运动方程,我还在学习与琢磨。连续方程貌似容易理解,但推敲起来,似乎也暗藏机关。芮老师书的讲解与张文华老师的讲解也不同。
陈:我对连续性方程用“水涨面高、水落面绌”来表达,公式仅用一行式子就可以直观解释。现在的连续方程,形状上就有点怪异。我是改用了等价、直观的表达。按标准连续方程直接解释,不太直观。
陈:目前的连续方程解释中,既然有dL,那么,dA的位置在断面的哪个部位?现有解释中说法不明确,现在我也不确定,难道是dL的各处皆可?
陈:连续性方程,貌似水量平衡,没有无源之水,但公式不太直观,学生怎么理解,我不太清楚,我感觉不易理解。我由于本科毕业不久即改行多年,水文学一些定量公式的学习与理解十分欠缺,现在因考虑做些水文史梳理,因此重新学习一些基础理论,彭老师一生以水文为业,没有我这样重操旧业、水文知识框架松散的困惑。
彭:dA 是在某个断面处的断面面积,而dL是满足单位长度的需要而除,或者理解为该断面左右的单位长的邻域。
彭:我现在花了很多时间在思考这些水力学方程的认知哲学基础。
彭:连续性方程用水库的入库来水与闸门放水理解,更直观。或者用每分钟从教室前门出去的学生数量减去同时段从后门进来的学生数量,就是这一分钟内教室学生数量的增加量来比喻。
陈:我用截图左上角的一行改写式,使连续方程直观易懂了。我不是说道理层面不易懂,而是说连续方程的式子不易懂,虽然表面上貌似简单。
彭:是的,将单位长与单位时间分别从分母位置改变为现在的样子,变微分为类似积分式,学生好理解一些。
微分方程要比差分方程抽象一些,用差分方程的模式讲,更直观。
陈:就我自己而言,这样一转换,感到亲切多了。
彭:从物理哲学意义上讲,严格地税,真正的可描述的物理过程,应该是差分式的,而不可能式微分式的。
陈:我这个转换,数学上不是很严谨,但便于理解。彭老师的差分方式,数学上比较严谨。呵呵,学生课时不通,要求不高的,能直观看懂就很好了。
彭:从物理哲学意义上讲,严格地税,真正的可描述的物理过程,应该是差分式的,而不可能是微分式的。
陈:我说的是普及,彭老师是提高。
彭:微分式是思维美学意义上的数学抽象。
陈:彭老师如觉得合适,我把这段对话放入科学网博客。放入前,我会发给彭老师过目一下。这至少反映了我对连续方程直观理解的思考。
彭:谢谢您。我的话,不是很严谨,但是意思大致就是如此。
附2:https://mp.weixin.qq.com/s?__biz=MjM5Mzg1MzEyOA==&mid=100000596&idx=1&sn=804aa469a873cb59175329c3548af887
水环境编Cheng长
圣维南方程组是对于明渠非恒定渐变流基本微分方程组,由连续性方程和运动方程式所构成,其广泛应用于洪水推流,描述了水流的运动,当然水质模型中,水动力是基础,所以其对水质模型也非常重要。其具体的推导方式就略过,这里仅仅描述结论。
连续性方程:
能量方程:
上述所述,这是每项的具体含义,简化为如下形式:
我们得到了Q和Ac的另一种关系,带入之前的曼宁公式,偏微分方程可以通过数值解进行求解。
附3:https://baike.baidu.com/item/%E5%9C%A3%E7%BB%B4%E5%8D%97%E6%96%B9%E7%A8%8B%E7%BB%84/4349325?fromtitle=%E5%9C%A3%E7%BB%B4%E5%8D%97%E6%96%B9%E7%A8%8B&fromid=5030851&fr=aladdin
1871年,法国科学家圣维南提出了圣维南方程组(英文: de Saint-Venant system of equations),该方程在气动、水力以及其他工程领域中有着广泛的应用。
含 义
描述水道渐变不恒定水流运动规律
1 简介
2 方程组的形式
3 基本假定
4 求解方法
5 参考书目
简介
编辑
描述水道和其他具有自由表面的浅水体中渐变不恒定水流运动规律的偏微分方程组。由反映质量守恒律的连续方程和反映动量守恒律的运动方程组成。1871年由法国科学家A.J.C.B.de圣维南提出,故名。
一百多年来,虽然为了考虑更多的因素和实际应用方便对它的基本假定作了某些简化或改进,产生出多种不同的表达形式,但其实质没有变化。主要进展表现在求解方法的改进和创新。1877年法国工程师克莱茨提出了瞬态法。1938年苏联С.А.赫里斯季安诺维奇提出另一类解法──特征线法。但均因计算量较大,不得不进行各种简化处理,使实际应用受到限制。自50年代以来,随着电子计算机的普及,研究和提出了一整套解法,并研究出若干个通用性较强的应用软件(即程序系统),促进了圣维南方程组在水文和其他工程领域中的应用。
方程组的形式
编辑
式中t为时间;s 为距水道某固定断面沿流程的距离;h、v、Z0分别为相应于s处过水断面的水深、断面平均流速和水底高程;Hf为由于摩阻损失而引起的能量比降;g为重力加速度;t和s为自变量;h和v为因变量;Z0、Hf可由s、h和v确定。(1)式为连续方程,反映了水道中的水量平衡,即蓄量的变化率(第一项)应等于沿程流量的变化率(第二项)。 (2)式为运动方程。其中第一项反映某固定点的局地加速度,第二项反映由于流速的空间分布不均匀所引起的对流加速度。以上两项称为惯性项。第三项反映由于底坡引起的重力作用,称为重力项。第四项反映了水深的影响,称为压力项。第三、四项可合并为一项,即水面比降。第五项为水流内部及边界的摩阻损失。该式表达了重力与压力的联合作用使水流克服惯性力和摩阻引起的能量损失而获得加速度。
圣维南方程组还有许多其他形式。例如:以断面流量代替流速,以面积代替水深作为因变量;也可考虑河道两侧的沿程入流、地转力和水面风力的影响;还可把垂线平均流速作为因变量,写出二维水体渐变不恒定明流的运动方程。
基本假定
编辑
①流速沿整个过水断面(一维情形)或垂线(二维情形)均匀分布,可用其平均值代替。不考虑水流垂直方向的交换和垂直加速度,从而可假设水压力呈静水压力分布,即与水深成正比;
②河床比降小,其倾角的正切与正弦值近似相等;
③水流为渐变流动,水面曲线近似水平。此外,在计算不恒定的摩阻损失Hf时,常假设可近似采用恒定流的有关公式,如曼宁公式(见河水运动)。
圣维南方程组描述的不恒定水流运动是一种浅水中的长波传播现象,通常称为动力波。因为水流运动的主要作用力是重力,属于重力波的范畴。如忽略运动方程中的惯性项和压力项,只考虑摩阻和底坡的影响,简化后方程组所描述的运动称为运动波。如只忽略惯性项的影响,所得到的波称为扩散波。运动波、扩散波及其他简化形式可以较好地近似某些情况的流动,同时简化计算便于实际应用。
求解方法
编辑
圣维南方程组在数学上属于一阶拟线性双曲型偏微分方程组。联解方程组并使其符合给定的初始条件和边界条件,就可得出不恒定水流的流速和水深(或其他因变量)随流程和时间的变化,即v=v(s,t)和h=h(s,t)。初始条件为某一起始时刻的水流状态,如水道沿程各断面的水深和流速。边界条件为所计算的水体的边界水流状态,如某一河段上、下游边界断面处的水位过程、流量过程或水位流量关系等。给定的初始条件和边界条件的数目和形式必须恰当,符合水流的性质,才能保证方程组的解存在和唯一,保证不致因数据的微小变化而使方程的解发生很大的变化。此时,问题称为是适定的,求解才有意义。
除特殊情况外,很难用解析方法求得圣维南方程组的解析解。一般只能通过数值计算获得个别情况的近似解。常用的数值计算方法主要有以下五类 [1] :①有限差分法。将所计算的水体按照一定的网格划分,每个网格点处的微分形式的圣维南方程组,用某种形式的差分方程组来逼近。边界条件也写成差分形成。然后逐时段地求解差分方程组,得出各网格点(如断面)处的水深及流速。根据所采用的差分计算方法的不同,对每一计算时段来说,或可逐个算出各网格点处的水力要素,或是必须联立求解各网点处的水力要素。前者称为显式差分法,后者称为隐式差分法。克莱茨提出的瞬态法就属于一种简化的显式差分法。②特征法。把圣维南方程组由偏微分方程组变换为在所谓“特征”上成立的常微分方程组,通常称为特征方程组。在空间为一维的情况下,“特征”的几何表示称为特征线,而在二维则为特征面。不恒定水流中的波动和干扰是沿“特征”传播的。用有限差分法联立求解表达“特征”几何位置的方程和特征方程组,即可求得所需的数值解。③有限单元法。把水体划分成几何形状简单的单元(如一维的直线段,二维的矩形、直边或曲边三角形等),在每一单元内,解用数学处理比较简单的内插函数来逼近。把圣维南方程组应用于每个单元,变换为积分形式,并根据某种准则(如逼近的残差最小)来确定内插函数中的待定系数便可定解。常用的是伽辽金半离散有限单元法。 ④有限元。常见的有限元计算方法有直接法、变分法、加权余量法及能量平衡法等。⑤有限分析法。在局部单元上线性化微分方程和插值近似边界的条件下,在局部单元上求微分方程的解析解,而构成整体的线性代数方程组。
除了求解完全或简化形式的圣维南方程组的上述解法,在水文学中多年来还对一维流动发展出许多简化计算方法。例如,把运动方程简化为计算时段内计算河段的蓄水量与出流量之间关系的方程,然后联立求解。同时,已对水文学中常用的方法与求解圣维南方程组的关系进行了研究。如应用广泛的马斯金格姆(曾译“马斯京根”)流量演算法,可列为扩散波中的特殊情形。水文学方法简单,而且能较好地适用于某些情况,今后仍将长期广泛地被应用。
对于非渐变的流动,水流通过激波把两部分渐变流连接起来。如通过水跃实现由急流(超临界流)到缓流(次临界流)的过渡。在涨潮和溃坝波中也常出现近乎垂直的波前。此时,两边的渐变流仍可用圣维南方程组来描述。只要补充激波处的跳跃条件和用以判别物理上是否许可的某种准则(如熵条件等)即可求解。
圣维南方程组所描述的具有自由表面的水体的渐变不恒定流动的计算具有重要的实际意义。洪流演进计算是洪水预报、堤防设计和防洪系统运用的重要依据。水电厂引水渠、下游尾水、灌溉和通航水道中的不恒定流计算,是确定堤岸、尾水管出口高程,论证发电、航运、给水等工程设施的安全和效益的根据。此外,潮汐河口的潮流计算、溃坝决堤造成的洪水灾害的估计等也都具有十分重要的经济意义。
参考书目
编辑
V.Yevjevich,K.Mahmood,ed.,Unsteady Flow in Open Channels,Vol.1,Water Resources Publ.,Fort Collins,Colorado,1975.