计算流体CFD2004-3.pptVIP

第三章 有限差分方程式 3.2 一般方法 推导： 推广： 3.3 高阶导数 中心差分二阶精度 3.4多维有限差分公式 多维 3.5 有限差分解的精度 3.5.1、差分离散的基本原则 1、相容性(Consistency) 2、收敛性(Convergence) 3、稳定性(Stability) 4、与物理问题的要求相适应 5、 Lax等价定理 3.5.2、差分离散的其它性质 2、输运性： 例如 3、耗散与频散 4、守恒性： * 华中科技大学能源与动力工程学院 基本思想 以离散网格点上的变量的差分形式代替微分，结果得到一个包括所有节点上未知量的联立代数方程组。 3.1 简单的方法 ——台劳级数展开 在i+1和i-1上的u ，对i点展开的Taylor级数 前差 后差 中心差 二阶中心差分 差分格式的精度(截断误差)容易得到 任意阶的精度的有限差分方程 利用相关的三个点推导出一阶导数的二阶精度差分方程 单边公式 二阶导数的一阶精度差分方程(单边三点公式 ) (3.2.1) (3.2.2a) (3.2.2b) aui+b*(3.2.2a)+c* (3.2.2b)得： 与(3.2.1)比较得： a=3/2 b= –2 c=1/2 引入位移算子E和微分算子D： 前向差分算子 一阶精度 二阶精度 (3.2.9) 三阶精度？ 由级数展开可知 后向差分算子 中心差分算子 二阶精度 (第一项) 另一个中心差分算子 二阶精度 四阶精度 中心差分——(整数+1/2) 的点 二阶精度 四阶精度 四阶精度 表3.3.1 不同阶的有限差分公式 (b) 后向差分,O( ) (e) 后向差分, O( ) (c) 中心差分,O( ) (f) 中心差分, O( ) 1 -2 1 -1 16 -30 16 -1 -1 0 1 1 -8 0 8 -1 (i-1,j+1) (i,j+1) (i+1,j+1) ? ? (i-1,j) (i+1,j) ? ? ? (i-1,j-1) (i,j-1) (i+1,j-1) 图3.4.1 二维空间网格 (i,j ) 利用中心差分和平均中心差分算子 利用中心差分算子 a b a b 1 1 1 1 1 1 〇 □ 〇 □ 〇 〇 〇 〇 〇 〇 -4b -4b □ 〇 □ 〇-4a 1〇 〇 〇1 1〇 〇 〇1 -4a -8 -20 〇 □ 〇 □ 〇 〇 〇 〇 〇 〇 a b a b 1 1 1 1 1 1 图3.4.3五点格式的奇偶震荡. (a) b =2/3. (b) b =1/3 图3.4.4 九点式公式 有限差分格式以及随后在边界值问题上的应用，都必须保证在描述模拟问题的物理关系上的精度。这个精度取决于下面定义的相容性，稳定性和收敛性。 该格式1940年提出，1951年被证明是错误的 实例：对流现象的方程和边界条件： 初始条件 t=0, u(x,0)=f(x) 边界条件 x=0, u(0,t)=f1(t) x=L, u(L,t)=f2(t) 0 f(x) L x t T f1(t) f2(t) 时间前差、空间前差 FTFS 时间前差、空间中心差 FTCS