有限差分方法24差分方程相容性收敛性和稳定性.pptVIP

计算力学基础 第二章有限差分方法 24差分方程的相容性、收敛性和稳定性 个微分方程采用不同的方法可以得到不同的差分方程。那么, 我们要问,对于这些不同的差分方程是否都同样有效,同样可靠,而 且能得到同样的计算结果呢? 答案是否定的。事实上,不同的差分方程和原方程有完仝不同的 对应关系,它们具有各自不同的性质,因此,数值结果也完全不同。 在这些差分方程中有些差分方程是有效的、可靠的;些差分方程只有 在一定的条件下是有效的、可靠的:有些差分方程则是完全无效的 不可靠的。所以,如何判断和分析差分方程有效性和可靠性就成为非 常必要和现实的问题了 在这一节中我们首先对差分方程有效性的一些基本概念(如相容 性、收敛性、稳定性)作简单介绍,为本章以后各节的分析讨论奠定 基础 24.1相容性( Consistency) 差分方程相容性是讨论当№,∧x→0时,差分方程逼近于微分方 程的程度,因此,相容性是讨论差分方程和微分方程的关系 定义:对于一足够光滑函数l,若时间步长At,空间步长Ax 趋近于0时,差分方程截断误差R对于每一点(x,)都趋近于 0,则该差分方程L=0逼近微分方程LM=0,即差分方程与 微分方程是相容的 差分方程相容性可以通过 Taylor展开方法来证明。例如,扩散方 程的FTCS差分格式为: 2u1+1-1=0 T 把+1作为的函数,在n邻域展开成Tayo级数,把t1 和l1作为x的函数,在x邻域展开成 Taylor级数 a=/ouM+1(ana2+可ar+o(△r) au of U,= △x2+ △x3+O(△x4) △x3+O(△x2) 2 ax 将L 小1和l11代入FCS格式中,即可得到: i a-u A2-O(43)-O(△x2) 2 at 6 at 当Δ,Δx→o时,上等式石侧所有项都趋近0,差分方程趋近 于原微分方程,即FTCS差分方程和原方程是相容的。 关于差分方程相容性需要作以下说明: ①相容性是对求解区域内任意一点差分方程逼近于微分方程的 程度,相容性是有限差分算法(包括有限体积算法)首先必 须满足的有效性条件