对流扩散方程有差分方法.docVIP

对流扩散方程有限差分方法 求解对流扩散方程的差分格式有很多种，在本节中将介绍以下3种有限差分格式：中心差分格式、Samarskii格式、Crank-Nicolson型隐式差分格式。 3.1 中心差分格式 时间导数用向前差商、空间导数用中心差商来逼近，那么就得到了（1）式的中心差分格式 （3） 若令 ，，则（3）式可改写为 （4） 从上式我们看到，在新的时间层上只包含了一个未知量，它可以由时间层上的值，，直接计算出来。因此，中心差分格式是求解对流扩散方程的显示格式。 假定是定解问题的充分光滑的解，将，，分别在处进行Taylor展开： 代入(4)式，有 显然，当，时，，即中心差分格式与定解问题是相容的。由以上的讨论也可得知，对流扩散方程的中心差分格式的截断误差为。 对于我们上面构造的差分格式，是否可以直接用于实际计算呢？也就是说，如果初始值有误差，在计算过程中误差会不会扩大传播呢？这就是接下来我们要讨论的是差分方程的稳定性问题。下面用Fourier方法来分析中心差分格式的稳定性。 令，代入到（4）式 整理得 所以该差分格式的增长因子为： 其模的平方为 由于，所以（即差分