计算水力学第三章.pptVIP

计算水力学 第三章 有限差分的基本理论 第一节 基本概念 一维对流方程 计算平面为ｘ－ｔ的上半平面。在平面上画出两族平行于坐标轴的直线,把求解域分成矩形的计算网格。网格线的交点称为节点,ｘ方向上网格线之间的距离Δｘ称为空间步长,ｔ轴方向上网格线之间的距离Δt称为时间步长 网格剖分使得每一空间步长、时间步长均相等，则称该网格为一均匀网格，否则称之为非均匀网格 数值解主要是求解节点上的末知变量的数值，利用有限的节点上的值来代替整个求解域内的连续函数值。 概念：离散、插值、误差 构造差分方程、分析数值误差 第二节 偏导数的差商近似 一、差分、差商的基本概念 解析函数 导数定义 差分 差商 向前差分 向后差分 中心差分 一阶导数，对应的差分称为一阶差分。对一阶差分再作一阶差分，所得到的称之为二阶差分。二阶向前差分： 任何阶差分都可以由其低一阶的差分得到： 函数的差分与自变量的差分之比，即为函数对自变量的差商 一阶向前差商 一阶向后差商 一阶中心差商 二阶中心差商 二、偏导数的差商近似 １．Ｔａｙｌｏｒ展开法 通过对差商近似点（ｉ，ｊ）的Ｔａｙｌｏｒ展开，可以分析差商对偏导近似的精度 一阶向前差商 一阶向后差商 二阶中心差商 边界处偏导数的差商近似 对点（０，ｊ）进行Ｔａｙｌｏｒ展开 构造一阶偏导数的二阶精度的差商近似必须有 解得 构造二阶偏导数的差商近似必须有 解得 构造二阶偏导数具有二阶精度的差商近似必须有 解得： ２．多项式插值法 用多项式插值法把待求函数表示成含待定系数的解析函数，由节点函数值确定该系数，然后对此函数求偏导数，得到逼近偏导数的差商表达式。 设函数ｕ可用抛物插值公式来近似： 设原点ｘ＝０在点ｉ的位置上，则有 解出待定系数 用高阶多项式插值可得到高阶差商表达式。 高阶多项式插值具有龙格不稳定性，使得插值对计算误差十分敏感。 多项式插值法在计算流体力学中多用于处理边界处的差商近似。 偏导数的差商近似还有其它多种方法，但最终均需用Ｔａｙｌｏｒ展开来计算其近似的误差，因此在实际计算中通常均用Ｔａｙｌｏｒ展开法来构造，因为此法在构造差商近似的同时还得出了其近似的误差精度。 第三节 差分方程 偏导数用其差商近似来代替 偏微分方程转变为相应的代数方程称之为差分方程。 对流方程 在点（ｉ，ｊ）成立 在点（ｉ，ｊ）的对流方程可以近似 Ｃｏｕｒａｎｔ数 差分方程 设Ｃ≥０，求解域（０≤ｘ≤ｌ）， 定解条件 离散表达 定解问题 ＦＴＢＳ格式差分方程 例定解问题 采用ＦＴＢＳ格式 ａ．取Δｘ＝１，Δｔ＝0.5 ， C=1.0 则 b．取Δｘ＝１，Δｔ＝2，C=1.0则 对同一定解问题的同一差分格式（ＦＴＢＳ）其不同的空间与时间步长，将得到不同的结果，如果作为原始定解问题的近似解，那一个解精度高呢？。 不稳定的解是不能作为原定解问题的近似解的。 偏导数的差商近似并非一种，同一偏微分方程的差分方程也并非一个，可以有若干个，对原始定解问题也相应有若干种差分格式。 ＦＴＣＳ格式 ＦＴＦＳ格式 蛙跳格式 显式格式：由第ｊ时间层上的值，可直接算出第ｊ＋１时间层上的值的格式。 隐式格式：不能直接从ｊ时间层上值直接解出，需联立求解ｊ＋１层上的值的格式。 对同一个定解问题，可以有多种差分格式，多种步长参数来近似，从而也得到若干个差分近似解。那么这些解是否可以都作为原定解问题的近似解？那些解精度高？为什么？ 相容性、稳定性及收敛性分析 第四节 截断误差和相容性 以ＦＴＢＳ格式为例 等价方程 截断误差 ＦＴＣＳ格式的截断误差 ＦＴＦＳ格式的截断误差 蛙跳格式的截断误差 对流方程的差分方程等价形式 差分方程和相应的微分方程相容 定解条件 差分算子 截断误差 定解问题相容 第五节 收敛性 相容性：是指当自变量的步长趋于零时,差分格式与微分问题的截误差的范数是否趋于零,从而可看出是否能用此差分格式来逼近微分问题。 收敛性：是指当自变量步长趋于零时,要求差分格式的解趋于微分方程定解问题的解。要求差分格式的解(数值解)与微分方程定解问题的解(精确解)是一致的。 差分格式的解 微分问题的解 离散误差 差分格式收敛 相容性是收敛性的必要条件, 相容性是形式上的逼近，收敛性是解的逼近，相容性不一定能保证收敛性 例微分方程定解问题 解析解为 将［０，ａ］等分为ｎ段 则步长Δｘ＝ａ／ｎ 差分解为 微分问题 ＦＴＢＳ格式 离散误差 由截断误差分析有 当Ｃ≥０和η≤１,即０≤η≤１,则 说明当０≤η≤１时,本