计算水力学第三章教学文案.pptVIP

计算水力学;第一节 基本概念;ｘ－ｔ平面、计算网格 ;网格节点 ;网格剖分使得每一空间步长、时间步长均相等，则称该网格为一均匀网格，否则称之为非均匀网格 数值解主要是求解节点上的末知变量的数值，利用有限的节点上的值来代替整个求解域内的连续函数值。 概念：离散、插值、误差 构造差分方程、分析数值误差;第二节 偏导数的差商近似;向前差分 向后差分 中心差分 ;一阶导数，对应的差分称为一阶差分。对一阶差分再作一阶差分，所得到的称之为二阶差分。二阶向前差分：;任何阶差分都可以由其低一阶的差分得到： 函数的差分与自变量的差分之比，即为函数对自变量的差商 一阶向前差商 ;一阶向后差商 一阶中心差商 二阶中心差商;二、偏导数的差商近似;一阶向前差商 一阶向后差商 ;二阶中心差商;边界处偏导数的差商近似; 构造一阶偏导数的二阶精度的差商近似必须有 解得 ; 构造二阶偏导数的差商近似必须有 解得 ;;构造二阶偏导数具有二阶精度的差商近似必须有 解得：;２．多项式插值法;设原点ｘ＝０在点ｉ的位置上，则有 解出待定系数 ;;用高阶多项式插值可得到高阶差商表达式。 高阶多项式插值具有龙格不稳定性，使得插值对计算误差十分敏感。 多项式插值法在计算流体力学中多用于处理边界处的差商近似。 偏导数的差商近似还有其它多种方法，但最终均需用Ｔａｙｌｏｒ展开来计算其近似的误差，因此在实际计算中通常均用Ｔａｙｌｏｒ展开法来构造，因为此法在构造差商近似的同时还得出了其近似的误差精度。 ;第三节 差分方程;;在点（ｉ，ｊ）的对流方程可以近似 Ｃｏｕｒａｎｔ数 差分方程 ;设Ｃ≥０，求解域（０≤ｘ≤ｌ）， 定解条件 离散表达;定解问题 ;ＦＴＢＳ格式差分方程 ;例定解问题 ;采用ＦＴＢＳ格式 ;ａ．取Δｘ＝１，Δｔ＝0.5 ， C=1.0 则 b．取Δｘ＝１，Δｔ＝2，C=1.0则 ;;;对同一定解问题的同一差分格式（ＦＴＢＳ）其不同的空间与时间步长，将得到不同的结果，如果作为原始定解问题的近似解，那一个解精度高呢？。 不稳定的解是不能作为原定解问题的近似解的。 偏导数的差商近似并非一种，同一偏微分方程的差分方程也并非一个，可以有若干个，对原始定解问题也相应有若干种差分格式。;ＦＴＣＳ格式 ;ＦＴＦＳ格式 ;蛙跳格式 ;显式格式：由第ｊ时间层上的值，可直接算出第ｊ＋１时间层上的值的格式。 隐式格式：不能直接从ｊ时间层上值直接解出，需联立求解ｊ＋１层上的值的格式。 对同一个定解问题，可以有多种差分格式，多种步长参数来近似，从而也得到若干个差分近似解。那么这些解是否可以都作为原定解问题的近似解？那些解精度高？为什么？ 相容性、稳定性及收敛性分析;第四节 截断误差和相容性;等价方程 截断误差 ;ＦＴＣＳ格式的截断误差 ＦＴＦＳ格式的截断误差 蛙跳格式的截断误差 ;对流方程的差分方程等价形式 差分方程和相应的微分方程相容 ;定解条件 差分算子 截断误差 定解问题相容 ;第五节 收敛性;差分格式的解 微分问题的解 离散误差 差分格式收敛 相容性是收敛性的必要条件, 相容性是形式上的逼近，收敛性是解的逼近，相容性不一定能保证收敛性 ;例微分方程定解问题 解析解为 将［０，ａ］等分为ｎ段 则步长Δｘ＝ａ／ｎ;差分解为;;;微分问题 ＦＴＢＳ格式 ;离散误差 由截断误差分析有 ;当Ｃ≥０和η≤１,即０≤η≤１,则 ;; 说明当０≤η≤１时,本问题的ＦＴＢＳ格式收敛。这种离散化误差的最大绝对值趋于零的收敛性情况称为一致性收敛。;第六节 稳定性;依赖区间AB 和决定域pAB;;影响区域 ;;同一微分问题,当采用不同差分格式时,其依赖区间、决定区域和影响区域可以是不一致的。 依赖区间、决定区域和影响区域是由差分格式本身的构造所决定的并与步长比 有关 ;定解问题;假设在第j层上的第i点,由于计算误差，得到。 设i＝０,ｊ＝０， ε＝１，相应于FTBS格式 ;;;;算例表明了当 值不同时计算误差所产生的影响。 误差逐渐衰减传播 误差无衰减传播 误差震荡放大传播 数值误差有不同的传播方式,格式使误差逐渐衰减传播称为差分格式稳定，否则称为不稳定。 ;单增长型的不稳定称为静力不稳定性;过冲型振荡的不稳定称为动力不稳定 ;von Neumann稳定性分析方法;误差展开成傅氏级数 代入误差传播方程;对任意的ｋ有 Ｇ为放大因子 ; FTBS格式稳定条件 ;FTBS格式稳定条件 ＦＴＣＳ格式为一不稳定格式 ＦＴＦＳ格式稳定条件 蛙跳格式稳定条件 ;V