第1章有限差分法基础.pptVIP

铸造过程计算机模拟讲义 复习-有限差分法 * * 主要内容 1、差分原理及逼近误差 2、差分方程，截断误差和相容性 3、收敛性与稳定性 4、Lax等价定理 第一节 差分原理及逼近误差/差分原理 向前差分 （1-2） 向后差分 （1-3） 中心差分 （1-4） 1．差分原理 第一节 差分原理及逼近误差/差分原理 对一阶差分再作一阶差分，所得到的称为二阶 差分，记为 。 以向前差分为例，有 （1-5） 第一节 差分原理及逼近误差/差分原理 函数的差分与自变量的差分之比，即为差商。 一阶向前差商为 一阶向后差商为 （1-7） （1-8） 第一节 差分原理及逼近误差/差分原理 一阶中心差商为 或 （1-9） （1-10） 第一节 差分原理及逼近误差/差分原理 二阶差商多取中心式，即 当然，在某些情况下也可取向前或向后的二阶差商。 （1-11） 第一节 差分原理及逼近误差/差分原理 多元函数f ( x,y,…)的差分与差商也可以类推。 如一阶向前差商为 （1-12） （1-13） 第一节 差分原理及逼近误差/逼近误差 差商与导数之间的误差表明差商逼近导数的程度，称为逼近误差。 由函数的Taylor展开，可以得到逼近误差相对于自变量差分（增量）的量级，称为用差商代替导数的精度，简称为差商的精度。 （1-15） 2．逼近误差 第一节 差分原理及逼近误差/逼近误差 一阶向后差商也具有一阶精度。 （1-16） 同理： 第一节 差分原理及逼近误差/逼近误差 将 与 的Taylor展开式相减可得 可见一阶中心差商具有二阶精度。 （1-17） 第一节 差分原理及逼近误差/逼近误差 这说明二阶中心差商的精度也为二阶。 （1-18） 将 与 的Taylor展开式相加可得 第二节 差分方程、截断误差和相容性/差分方程 如果将微分方程中的导数用相应的差商近似代替，就可得到有限形式的差分方程。 现以对流方程为例，列出对应的差分方程。 （2-1） 图2-1 差分网格 第二节 差分方程、截断误差和相容性/差分方程 若时间导数用一阶向前差商近似代替，即 空间导数用一阶中心差商近似代替，即 则在 点的对流方程就可近似地写作 （2-2） （2-3） （2-4） 第二节 差分方程、截断误差和相容性/差分方程 第二节 差分方程、截断误差和相容性/截断误差 按照前面关于逼近误差的分析知道，用时间向前差商代替时间导数时的误差为 , 用空间中心差商代替空间导数时的误差为 因而对流方程与对应的差分方程间也存在一个误差： （2-5） 第二节 差分方程、截断误差和相容性/截断误差 一个与时间相关的物理问题，还须给定初始条件， 从而形成完整的初值问题。对流方程的初值问题为： 相应差分格式： （2-7） （2-8） 第二节 差分方程、截断误差和相容性/截断误差 FTCS格式 （2-9） FTFS格式 （2-10） 第二节 差分方程、截断误差和相容性/截断误差 （2-11） FTBS格式 第二节 差分方程、截断误差和相容性/截断误差 (a) FTCS (b)FTFS (c)FTBS 图2-2 差分格式 第二节 差分方程、截断误差和相容性/截断误差 FTCS格式的截断误差为 FTFS和FTBS格式的截断误差为 （2-12） （2-13） 3种格式对 都有一阶精度。 第二节 差分方程、截断误差和相容性/相容性 如果当 时，差分方程的截断误差的某种范数 也趋近于零，即 则表明从截断误差的角度来看，此差分方程是能用来逼近 微分方程的，通常称这样的差分方程和相应的微分方程相容。 （2-17） 第二节 差分方程、截断误差和相容性/相容性 只有方程相容，定解条件也相容，即 和 整个问题才相容。 （2-21） 第三节 收敛性与稳定性/收敛性 当自变量步长趋于零时，要求差分格式的（精确）解趋于微分方程定解问题的（精确）解。我们称这种是否趋于微分方程定解问题的解的情况为差分格式的收敛性。 ，也是微分问题定解区域上的一固定点，设差分格式在此点 的解为 , 相应的微分问题的解为 ，二者之差为 （3-