2.4 差分方程相容性、收敛性与稳定性.pptVIP

计 算 力 学 基 础; 一个微分方程采用不同的方法可以得到不同的差分方程。那么，我们要问，对于这些不同的差分方程是否都同样有效，同样可靠，而且能得到同样的计算结果呢？ 答案是否定的。事实上，不同的差分方程和原方程有完全不同的对应关系，它们具有各自不同的性质，因此，数值结果也完全不同。在这些差分方程中有些差分方程是有效的、可靠的；些差分方程只有在一定的条件下是有效的、可靠的；有些差分方程则是完全无效的、不可靠的。所以，如何判断和分析差分方程有效性和可靠性就成为非常必要和现实的问题了。 在这一节中我们首先对差分方程有效性的一些基本概念（如相容性、收敛性、稳定性）作简单介绍，为本章以后各节的分析讨论奠定基础。 ; 差分方程相容性是讨论当 时，差分方程逼近于微分方程的程度，因此，相容性是讨论差分方程和微分方程的关系。 定义：对于一足够光滑函数 ，若时间步长 ，空间步长 趋近于0时，差分方程截断误差 对于每一点 都趋近于0，则该差分方程 逼近微分方程 ，即差分方程与微分方程是相容的。 差分方程相容性可以通过Taylor展开方法来证明。例如，扩散方程的FTCS差分格式为： ;把 作为t的函数，在 邻域展开成Taylor级数，把 和 作为x的函数，在 邻域展开成Taylor级数： ; 当 时，上等式右侧所有项都趋近0，差分方程趋近于原微分方程，即FTCS差分方程和原方程是相容的。 关于差分方程相容性需要作以下说明： 相容性是对求解区域内任意一点差分方程逼近于微分方程的程度，相容性是有限差分算法（包括有限体积算法）首先必须满足的有效性条件。 ;相容性要求对于求解区域内任意点 ，在 同时趋近于0，截断误差 趋近于0。如果 不是同时趋近于0或并不趋近于0，而是趋近于某值，或结论并不是对每个点 都成立，则差分方程就不满足相容性条件，差分方程也就不逼近于微分方程。 相容性条件不仅要求差分方程截断误差 趋近于0，而且要求差分方程定解条件截断误差 也同时趋近于0。 差分格式有两种不同形式的相容性，即无条件相容和有条件相容。;2.4.2 收敛性（Convergence ）;差分方程收敛性有两种证明方法，直接证明法和数值试验法。;由（b）式可以看出离散化误差方程在形式上和差分方程是完全相同的，由此可以得到：;式中 表示在n层的所有节点上离散化误差 绝对值最大值，对于所有节点j有：;由此可得到：;二、数值试验法;关于差分方程收敛性需要作以下说明：; 粗看起来，差分方程相容性要求时，差分方程逼近于微分方程，似乎差分方程数值解也应该收敛于微分方程精确解。事实上，当我们在证明相容性时，已经假定了差分方程数值解就是微分方程精确解，在对微分方程进行展开时，截断误差中已经忽略了离散化误差的存在。因此，差分方程相容性并不能保证其收敛性。;2.4.3 稳定性（Stability ）;定义：在某一个时刻tn存在计算误差 ，若在 时刻满足：;下面我们用几个简单的例子来说明差分方程稳定性概念。 ;在（a）式中设在tn时刻xj的计算误差为 ，而计算到n+100时刻，（xj，tn+100）点的计算误差将发展到 ，假定只有在节点（xj，tn）上存在误差 ，其他各节点的计算误差为零，则若取r=0.8，则 。由此可以看出，这个计算误差必定会将差分方程精确解原来面目完全淹没了，所求得差分方程数值解已经没有任何意义了，因此，FTFS差分方程是不稳定的。 ;（2）对流方程FTBS差分格式的误差传播方程为：;由此可知，在n时刻的计算误差 是不会大于 ，因此，当a＞0， 时，FTBS差分格式是稳定的（见图a）。这是有条件的稳定，稳定的条件是a＞0， 。但是，对于不同的a，Δt，Δx，FTBS差分格式的稳定条件是不同的（见图b）。;当a=1，Δx=0.1，r=0.8，则有： ； 当a=1