第二章 有限差分基础.docVIP

PAGE 1 PAGE 13 有限差分基础（finite difference method，FDM） 偏微分方程的一般形式 （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 1 ） 网格划分 一般有限差分采用结构化网格划分。即节点对应于当地坐标系统的原点。它的轴同网格线一致。即两个同一族的网格线不相交，且没对网格线对应不同的族。每一个节点可用唯一的一个坐标表示，如（?1, ?2）。网格线能用?1＝const, ?2＝const 表示。 i-1i i-1 i i+1 N 1 N Nj j+1 j-1 j 1 1 i-1 i i+1 (i,j) 2D 有限差分法就是要将节点上的偏微分方程用相邻点上的值表示，变成线性代数方程式。 为流体力学的微分方程的数值求解方法之代表。必要条件： 连续领域内的分配有限的网格 领域内的函数分布可用网格点上的值代表 PWEEENNSSSWW P W E EE NN S SS WW N P P E W B T S N P P E W S N 5点计算分子 15点计算分子 7点计算分子 这些节点又称为计算分子。方程的个数应与未知数相同，即每个节点有一个方程式。 Taylor展开 例如：一维时间变量?的理论解为?(t,x)，它在离散点上的值为投影（projection）的近似值为：, n: 时间的step数 i:空间的step数 为了求得此近似解，需对微分方程进行差分近似。利用Taylor展开可得到几个差分表示形式，仅考虑空间依存问题： 在?x很小时，位置j?x内的物理量?用?j来表示，则位置(j+1)?x上的值?j+1表示为： （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 2 ） (j-1)?x上的值?j+1表示为： （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 3 ） 基本差分格式 一阶导数（first derivative）的近似 （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 4 ） 向前差分(forward difference，FDS) 利用( 1 ) 式，可得到1阶微分的向前差分形式: （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 5 ） 向后差分(backward difference，BDS) 利用( 1 ) 式，可得到1阶微分的向后差分形式 （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 6 ） 中心差分（central difference，CDS） (1)-(2) 得1阶微分的2 次精度中心差分法： （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 7 ） 上风法、迎风法（upwind difference, UDS） 与速度有关的微分 （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 8 ） 二阶导数的近似 中心差分（central difference，CDS） 利用（j?1/2）?x的Tayor展开，可得过且1阶微分的2次精度的向前向后差分形式： （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 9 ） 将上二式相减，得2阶微分的差分方程式中心差分：(?x相当) （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 10 ） 其它还有各种形式。但最常用的是中心差分。对于扩散项常采用 （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 11 ） 向后差分 向后差分 向前差分 中心 差分 精确值 i+1 i i-1 ?xi ?xi+1 ? x 数值误差 截断误差（truncation error） 例如：热传递方程 （ STYLEREF 1 \s 2 SEQ （ \* ARABIC \s 1 12 ） 时间向前差分，空间中心差分得