计算水力学--主要内容.ppt

1、本 课 主 要 内 容,水文09级计算水力学教学课件,绪论 明渠一维非恒定流的基本方程 有限差分的基本理论 河道水流计算 河网水流计算,本课主要内容,第二章 非恒定流基本方程,水文09级计算水力学教学课件,基本假设与定律 六条基本假设 两大基本定律 基本方程 连续方程 动量方程,本 章 要 点,方程的其他形式 Z，Q为因变量 Q，h为因变量 Q，A为因变量 U，h为因变量 U，Z为因变量 漫洪滩地的处理 动量校正系数 调蓄滩地宽度,1. 基本假设与定律,1、定床情况，即假设河床高程与时间无关。 2、断面代表水位，不考虑横比降。 3、浅水问题，满足静水压力分布规律。 4、水为不可压液体。 5、河

2、床底坡很小。 6、恒定流阻力公式仍然适用。,基本假设,1. 基本假设与定律,质量守恒定律： 单位时间内通过控制面流进控制体的净质量，等于同时段内控制体的质量增量。,基本定律：质量守恒定律、动量守恒定律,1. 基本假设与定律,动量守恒定律： 单位时间内通过控制面流进控制体的净动量与作用于控制体的外力之矢量和等于同时段内控制体的动量增量。,基本定律动量守恒定律,2. 基本方程,圣维南方程组：连续方程、动量方程,一般情况下：Vx=0,3. 方程的其他形式及问题讨论,方程的其他形式 Z，Q为因变量 Q，h为因变量 Q，A为因变量 U，h为因变量 U，Z为因变量 漫洪滩地的处理 动量校正系数 调蓄滩地宽

3、度,不同的研究目的， 对圣维南方程组的 变换形式,问题关键： 滩地、主槽一起计算,动量校正系数,漫洪滩地的处理,调蓄滩地宽度,反应河道断面流速分布均匀性系数,4. 定解条件及定解问题,定解问题：基本方程、定解条件 定解条件：初始条件、边界条件 定解问题的适定性：解的存在性、唯一性和稳定性。,定解条件少了： 定解条件多了： 定解条件适当：,欠定 过定 适定,变量 在 所确定的线上为 常数,定解条件,对流方程,拟线性双曲型偏微分方程组,定解条件,圣维南方程的定界条件,S-V方程的定解条件,急流,缓流,反向急流,圣维南方程的定界条件及组合情况,缓 流,急 流,洪水波的分类,描述河道一维水流运动的圣维

4、南方程组：,运动波 惯性波 扩散波 动力波,运动波,水位流量之间呈单一关系 下游扰动不影响上游断面 波形传播波峰不变无耗散,惯性波,忽略摩阻项 假定底坡水平 棱柱形河道,不计摩阻损失，波动在传播过程中只有能量的转换， 无能量损失。,扩散波,扩散项的存在所以洪水波的波峰会逐渐坦化。,Z,ZMax,绳套形水位流量关系,Q,QMax,涨水,落水,水位或流量在短期内有大幅度的变化时：,动力波,运动波、惯性波和扩散波是动力波的特殊情况,感潮河道中的水流运动 闸门启闭引起的水流波动,这种情况下，动量方程式中的各项均不能忽略， 这样一种波动称为动力波。,动力波是所有波动中最复杂的，只能用完全的圣维南方程组描

5、述。,洪水波的简化方法水库调洪演算,一般情况下f (Q)的函数关系为非线性，难于用显式表达，故常用图解法或试算法求解。,连续方程式严格满足，并写成差分形式,基本假定：,马斯京根法,假定河段槽蓄量V与出流量 Q 及入流量 I之间存在着线性关系,洪水波的简化方法马斯京根法,连续方程式严格满足,第三章 有限差分的基本理论,水文09级计算水力学教学课件,差分、差商、差分方程 截断误差和离散误差 相容性、收敛性和稳定性 LAX等价定理 差分方程的数值效应 逆风效应与逆风格式,本 章 要 点,1. 偏导数的差商近似, 差分、差商的基本概念,解析函数：,导数定义：,差 分：,差 商：,通过对差商近似点（i，

6、j）的Taylor展开，可以分析差商对偏导近似的精度,1. 偏导数的差商近似Taylor展开法,例 边界处偏导数的差商近似,对点(0, j)进行Taylor展开,1. 偏导数的差商近似Taylor展开法,构造一阶偏导数的二阶精度的差商近似必须有,2. 偏导数的差商近似Taylor展开法,解 得：,可 得：,解 得：,可 得：,偏导数用其差商近似来代替 偏微分方程转变为相应的代数方程称之为差分方程。,3. 差分方程,对流方程,定解条件,3. 差分方程定解问题,3. 差分方程定解问题,显式格式：由第j时间层上的值，可直接算出第j+1时间层上的值的格式。 隐式格式：不能直接从j时间层上值直接解出，需

7、联立求解j+1层上的值的格式。,3. 差分方程定解问题,对同一个定解问题，可以有多种差分格式，多种步长参数来近似，从而也得到若干个差分近似解。那么这些解是否可以都作为原定解问题的近似解？那些解精度高？为什么？,相容性、稳定性及收敛性分析,4. 截断误差和相容性, 以FTBS格式为例,等价方程,截断误差,相容性是收敛性的必要条件，稳定性与收敛性有一定的联系。 Lax等价定理阐述相容性、收敛性和稳定性三者之间关系。,7. Lax等价定理,6. 稳定性,von Neumann分析,稳定性,数值试验,6. 稳定性 von Neumann稳定性分析,定解问题,误差传播方程,误差展开成傅氏级数,G为放大因

8、子,FTBS格式稳定条件,6. 稳定性,蛙跳格式稳定条件,FTBS格式稳定条件,FTCS格式为不稳定 格式,FTFS格式稳定条件,马斯京根法稳定条件,微分方程是描述物理量在时间和空间上的连续变化的规律。 差分方程来描述离散化后物理量的变化规律。 离散误差使原系统的物理性质和规律遭到歪曲和破坏的作用称为数值效应或离散近似的伪物理效应。 必须对这些效应有明确的概念，从物理上来考虑数值格式的合理性，减少数值效应的影响。,8. 差分方程数值效应,8. 差分方程数值效应数值耗散与弥散,FTBS格式,FTFS格式,FTCS格式,蛙跳格式,马斯京根格式,物质的对流输运出现了与波速相反方向传播的不合理现象，称

9、为“逆风”效应，是一伪物理现象的数值效应。,8. 差分方程数值效应“逆风”效应,“逆风”效应,采用何种格式还与波速的方向有关，例如潮水河道，则可以采用“逆风”格式,FTCS格式：二阶精度，无条件不稳定的。,FTFS格式：一阶精度，当C为负时采用。,FTBS格式：一阶精度，当C为正时采用。,8. 差分方程数值效应“逆风”效应,“逆风”格式,第四章 河道水流计算,水文09级计算水力学教学课件,课 程 内 容,蛙跳格式 LAX-Wendroff格式 Abbott隐式格式 Presissmann隐式格式 四点线性隐格式 内边界处理,Preissmann 隐式格式,四点线性隐格式,四点线性隐格式离散,连

10、续方程,动量方程,四点线性隐格式离散,为书写方便，忽略上标n+1，任一河段差分方程写成：,四点线性隐格式离散,追赶法求解,1、水位边界条件的计算,2、流量边界条件的计算,3、水位流量关系边界条件的计算,追赶法求解,在河道水流计算中，除了外部边界条件外，还可能遇到内部边界条件。 内部边界条件是指：河道的几何形状的不连续或水力特性的不连续点。例如，集中入流，过水断面突然放大，堰闸过流等等。 内边界处，圣维南方程组不再适用，必须根据其水力特性作特殊处理。内部边界条件通常包含两个相容条件，即流量的连续性条件和能量守恒条件(或动量守恒条件)。,内边界的处理,内边界的处理,对于内边界的处理归结于特殊河段的追赶系数计算。计算依据的是特殊河段的相容方程(水量守恒与动量守恒，这些方程与非恒定流的基本方程无关)，对相容方程进行必要的处理，单独计算特殊河段的追赶系数，可同正常河段一样求解。,内边界的处理,第五章 河网水流计算,水文09级计算水力学教学课件,课 程 内 容,1. 河网的基本概念,树状河网,环状河网,河网的节点,有调蓄节点 无调蓄节点,1. 河网的基