首页 > 其他分享 >圣维南方程组

圣维南方程组

时间:2024-05-06 10:37:51浏览次数:14  
标签:方程 求解 维南 方程组 计算 水流

圣维南微分方程组包括两个方程组, 分别是质量守恒方程和动量守恒方程。质量守恒方程描述了流体的质量守恒规律,即在没有外力作用下,流体的质量应该保持不变。动量守恒方程描述了流体的动量守恒规律,即在流体流动过程中,流体的动量应该保持不变。

>>圣维南方程和ns方程的关系:https://www.zhihu.com/question/36902274/answer/419364508

圣维南方程忽略了z和y?只考虑沿着河道流动,把三维的河道简化成一条线?

考虑因素:河道弯曲程度?河道的高程斜率?

一维单宽水流情况下,圣维南方程组的典型形式为:式中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]:①有限差分法。将所计算的水体按照一定的网格划分,每个网格点处的微分形式的圣维南方程组,用某种形式的差分方程组来逼近。边界条件也写成差分形成。然后逐时段地求解差分方程组,得出各网格点(如断面)处的水深及流速。根据所采用的差分计算方法的不同,对每一计算时段来说,或可逐个算出各网格点处的水力要素,或是必须联立求解各网点处的水力要素。前者称为显式差分法,后者称为隐式差分法。克莱茨提出的瞬态法就属于一种简化的显式差分法。②特征法。把圣维南方程组由偏微分方程组变换为在所谓“特征”上成立的常微分方程组,通常称为特征方程组。在空间为一维的情况下,“特征”的几何表示称为特征线,而在二维则为特征面。不恒定水流中的波动和干扰是沿“特征”传播的。用有限差分法联立求解表达“特征”几何位置的方程和特征方程组,即可求得所需的数值解。③有限单元法。把水体划分成几何形状简单的单元(如一维的直线段,二维的矩形、直边或曲边三角形等),在每一单元内,解用数学处理比较简单的内插函数来逼近。把圣维南方程组应用于每个单元,变换为积分形式,并根据某种准则(如逼近的残差最小)来确定内插函数中的待定系数便可定解。常用的是伽辽金半离散有限单元法。 ④有限元。常见的有限元计算方法有直接法、变分法、加权余量法及能量平衡法等。⑤有限分析法。在局部单元上线性化微分方程和插值近似边界的条件下,在局部单元上求微分方程的解析解,而构成整体的线性代数方程组。

除了求解完全或简化形式的圣维南方程组的上述解法,在水文学中多年来还对一维流动发展出许多简化计算方法。例如,把运动方程简化为计算时段内计算河段的蓄水量与出流量之间关系的方程,然后联立求解。同时,已对水文学中常用的方法与求解圣维南方程组的关系进行了研究。如应用广泛的马斯金格姆(曾译“马斯京根”)流量演算法,可列为扩散波中的特殊情形。水文学方法简单,而且能较好地适用于某些情况,今后仍将长期广泛地被应用。

对于非渐变的流动,水流通过激波把两部分渐变流连接起来。如通过水跃实现由急流(超临界流)到缓流(次临界流)的过渡。在涨潮和溃坝波中也常出现近乎垂直的波前。此时,两边的渐变流仍可用圣维南方程组来描述。只要补充激波处的跳跃条件和用以判别物理上是否许可的某种准则(如熵条件等)即可求解。

圣维南方程组所描述的具有自由表面的水体的渐变不恒定流动的计算具有重要的实际意义。洪流演进计算是洪水预报、堤防设计和防洪系统运用的重要依据。水电厂引水渠、下游尾水、灌溉和通航水道中的不恒定流计算,是确定堤岸、尾水管出口高程,论证发电、航运、给水等工程设施的安全和效益的根据。此外,潮汐河口的潮流计算、溃坝决堤造成的洪水灾害的估计等也都具有十分重要的经济意义。

标签:方程,求解,维南,方程组,计算,水流
From: https://www.cnblogs.com/2008nmj/p/18174399

相关文章

  • 扩展欧几里得 解二元一次方程组
    扩展欧几里得: 最大公约数-OIWiki(oi-wiki.org) intexgcd(inta,intb,int&x,int&y){if(!b){x=1,y=0;returna;}intd=exgcd(b,a%b,y,x);inttmp=x;x=y;y=tmp-(a/b)*y;returnd;}P5656【模板】二元一次不定......
  • 【习题】5.2 一阶常系数线性微分方程组
    [T050201]求方程组\(\begin{cases}x'+y'=y+z\\y'+z'=z+x\\z'+x'=x+y\end{cases}\)的通解.    解将方程组写成矩阵形式:\[\begin{aligned}&\begin{pmatrix}1&1&0\\0&1&1\\1&0&1\end{pmatrix}\begin{pmat......
  • 鄂维南院士:应用数学新时代的到来
    鄂维南院士:应用数学新时代的到来Kepler范式和Newton范式自从Newton(牛顿)时代以来,做科学研究存在着两种不同的范式:Kepler(开普勒)范式(paradigm)和Newton范式.在Kepler范式或称为数据驱动的(datadriven)方法中,人们通过对数据的分析提炼科学发现.经典的例子是行星运动的Ke......
  • 线性代数 方程组的几何解释
    本次,我们将讲述线性代数的基础——求线性方程组,下面从方程组讲起,它有\(n\)个未知数以及\(n\)个方程,方程数量与未知数个数相等,这是最普遍的状况。我们会了解到“行图像(\(RowPictrue\))”和“列图像(\(ColumnsPictrue\))”,行图像想必大家都见过就是两个方程的函数图像交于一......
  • 第一章 矩阵与线性方程组
    §1.1域§1.2矩阵的基本运算§1.3\(\text{Gauß}\)消元与矩阵的相抵标准型第一节域定义:一个群是指一个集合\(G(|G|\ge2)\)与二元运算\(*\),满足:封闭性:\(\foralla,b\inG,a*b\inG;\)结合律:\(\foralla,b,c\inG,(a*b)*c=a*(b*c);\)存在恒元:\(\existse\inG.......
  • 第一章 矩阵与线性方程组
    §1.1域§1.2矩阵的基本运算§1.3\(\text{Gauß}\)消元与矩阵的相抵标准型第一节域定义:一个群是指一个集合\(G(|G|\ge2)\)与二元运算\(*\),满足:封闭性:\(\foralla,b\inG,a*b\inG;\)结合律:\(\foralla,b,c\inG,(a*b)*c=a*(b*c);\)存在恒元:\(\existse\inG.\fo......
  • 一次线性方程组 高斯消元笔记
    高斯消元原理高斯消元用来解如下形式的方程组:\[\begin{cases}a_{1,1}x_1+a_{1,2}x_2+\cdots+a_{1,n}x_n=b_1\\a_{2,1}x_1+a_{2,2}x_2+\cdots+a_{2,n}x_n=b_2\\\cdots\\a_{n,1}x_1+a_{n,2}x_2+\cdots+a_{n,n}x_n=b_n\end{cases......
  • 克拉默法则求方程组
    例1:求二元一次方程组把方程写成矩阵的形式:第1个矩阵为系数矩阵(方阵),第2个矩阵为变量矩阵根据克拉默法则,xi=Di/D,Di表示第i列被最后那个列向量替换后的行列式,D为系数矩阵行列式 例2:三元一次方程组把方程写成矩阵形式:根据克拉默法则,x,y,z的解为3阶行列式可以用......
  • 麦克斯韦方程组
    麦克斯韦方程组\(\newcommand{\big}{\displaystyle}\)\(\newcommand{\d}{\text{d}}\)\(\newcommand{\e}{\epsilon}\)到目前为止,我们已经零碎地研究过麦克斯韦方程组。现在我们开始讨论完整地电磁场理论,对于可能以任何形式随时间变化的电磁场,我们将有完整而正确的描述。完整的麦......
  • 【数值分析】第6章-解线性方程组的迭代法
    第6章-解线性方程组的迭代法\[A\vec{x}=\vec{b}\Leftrightarrow\vec{x}=B\vec{x}+\vec{f}\]建立迭代\[\vec{x}^{(k+1)}=B\vec{x}^{(k)}+\vec{f}\]B称为迭代矩阵Jacobi迭代的矩阵形式\[\begin{align}A\vec{x}=\vec{b}&\Leftrightarrow(D+L+U)\vec{x}=......