有限差分方法基础.ppt

材料计算机数值模拟讲义 有限差分法 (1)请推导出y=f(x)的二阶向前差分，向后差分以及中心差分。 (2)请推导出对流方程 ? ? FTCS格式，并指出其截断误差的精度。 ? (3)什么是差分问题的收敛性和稳定性？二者有何关系？ 第二节 差分方程、截断误差和相容性/相容性（1/3） 一般说来，若微分方程为 其中D是微分算子，f是已知函数，而对应的差分方程为 其中 是差分算子，则截断误差为 这里 为定义域上某一足够光滑的函数，当然也可以取微分方程的解 。 （2-14） （2-15） （2-16） 如果当 、 时，差分方程的截断误差的某种范数 也趋近于零，即 则表明从截断误差的角度来看，此差分方程是能用来逼近微分方程的，通常称这样的差分方程和相应的微分方程相容（一致）。 如果当 、 时，截断误差的范数不趋于零，则称为不相容（不一致），这样的差分方程不能用来逼近微分方程。 （2-17） 第二节 差分方程、截断误差和相容性/相容性（2/3） 若微分问题的定解条件为 其中B是微分算子，g是已知函数，而对应的差分问题的定解条件为 其中 是差分算子，则截断误差为 （2-18） （2-19） （2-20） 第二节 差分方程、截断误差和相容性/相容性（3/3） 只有方程相容，定解条件也相容，即 和 整个问题才相容。 （2-21） 无条件相容 条件相容 以上3种格式都属于一阶精度、二层、相容、显式格式。 第三节 收敛性与稳定性/收敛性（1/6） ，也是微分问题定解区域上的一固定点，设差分格式在此点 的解为 , 相应的微分问题的解为 ，二者之差为 称为离散化误差。如果当 时，离散化误差的某种范数 趋近于零，即 则说明此差分格式是收敛的，即此差分格式的解收敛于相应微分问题的解， 否则不收敛。与相容性类似，收敛又分为有条件收敛和无条件收敛。 （3-1） 、 （3-2） 第三节 收敛性与稳定性/收敛性（3/6） 相容性不一定能保证收敛性，那么对于一定的差分格式，其解能否收敛到相应微分问题的解？答案是差分格式的解 收敛于微分问题的解是可能的。至于某给定格式是否收敛，则要按具体问题予以证明。下面以一个差分格式为例， 讨论其收敛性： 微分问题 的FTBS格式为 在某结点（xi , tn）微分问题的解为 ，差分格式的解为 ，则离散化误差为 （3-6） （3-5） （3-4） 第三节 收敛性与稳定性/收敛性（4/6） 按照截断误差的分析知道 以FTBS格式中的第一个方程减去上式得 或写成 若条件 和 成立，即 ，则 式中 表示在第n层所有结点上 的最大值。 （3-7） （3-8） （3-9） （3-10） 第三节 收敛性与稳定性/收敛性（5/6） 由上式知，对一切i有 故有 于是 综合得 （3-11） （3-13） （3-12） （3-14） 第三节 收敛性与稳定性/收敛性（6/6） 由于初始条件给定函数 的初值，初始离散化误差 。并且 是一有限量，因而 可见本问题FTBS格式的离散化误差与截断误差具有相同的量级。最后得到 这样就证明了，当 时，本问题的RTBS格式收敛。这种离散化误差的最大绝对值趋于零的 收敛情况称为一致收敛。 。 （3-15） （3-16） 此例介绍了一种证明差分格式收敛的方法，同时表明了相容性与收敛性的关系：相容性是收敛性的必要条件， 但不一定是充分条件，还可能要求其他条件，如本例就是要求 。 第三节 收敛性与稳定性/稳定性（1/8） 首先介绍一下差分格式的依赖区间、决定区域和影响区域。还是以初值问题 （3-17） (a) FTCS (b) FTFS (c) FTBS 图3-1 差分格式的依赖区间 第三节 收敛性与稳定性/稳定性（2/8） FTCS格式 (b) FTFS格式 (c) FTBS格式 图3-2 差分格式的影响区域 第三节 收敛性与稳定性/稳定性（3/8） 其解为零，即 若用FTBS格式计算，且计算中不产生任何误差，则结果也是零，即 当采用不同差分格式时，其依赖区间、决定区域和影响区域可以是不一样的。依赖区间、决定区域和影响区域是由差分格式本身的构造所决定的，并与步长比 有关。 （3-18） （3-19） （3-20） 假设在第k层上的第j点，由于计算误差得到 不妨设k=0, j=0, ，即