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

对流扩散方程有限差分方法 求解对流扩散方程的差分格式有很多种，在本节中将介绍以下3种有限差分格式：中心差分格式、Samarskii格式、Crank-Nicolson型隐式差分格式。 3.1 中心差分格式 时间导数用向前差商、空间导数用中心差商来逼近，那么就得到了（1）式的中心差分格式 （3） 若令 ，，则（3）式可改写为 （4） 从上式我们看到，在新的时间层上只包含了一个未知量，它可以由时间层上的值，，直接计算出来。因此，中心差分格式是求解对流扩散方程的显示格式。 假定是定解问题的充分光滑的解，将，，分别在处进行Taylor展开： 代入(4)式，有 显然，当，时，，即中心差分格式与定解问题是相容的。由以上的讨论也可得知，对流扩散方程的中心差分格式的截断误差为。 对于我们上面构造的差分格式，是否可以直接用于实际计算呢？也就是说，如果初始值有误差，在计算过程中误差会不会扩大传播呢？这就是接下来我们要讨论的是差分方程的稳定性问题。下面用Fourier方法来分析中心差分格式的稳定性。 令，代入到（4）式 整理得 所以该差分格式的增长因子为： 其模的平方为 由于，所以（即差分格式稳定）的充分条件为 上式可以改写为 注意到，所以上面不等式满足的条件为 ， 。 由此得到差分格式（3）的稳定性限制为 ， 。 故有结论：对流扩散方程的中心差分格式是条件稳定的。根据Lax等价定理，我们可以知道，对流扩散方程的中心差分格式是条件收敛的。 3.2 Samarskii格式 设0，先对方程(1)作扰动，得到另一个对流扩散方程 （5） 其中 ，当时，（5）式化为（1）式 对于（5）式，构造迎风格式 （6） 差分格式（6）称为逼近对流扩散方程的Samarskii格式。 首先推导（6）的截断误差。设是对流扩散方程（1）式的充分光滑的解 令 用 Taylor级数展开有 再令 用 Taylor级数展开有 由于 所以 当，时，，所以Samarskii格式与定解问题是相容的，并且其截断误差为。 现在看看Samarskii格式的稳定性。将（6）式两边同时加上，把（6）式化为 令 ，则上式即为： 根据中心显示格式稳定性的讨论，可以得到（6）式的稳定性条件为 ， 即 ， 稳定性的第二个条件等价于 而 利用不等式 所以 利用稳定性的第一个条件，有 ，从而可知稳定性条件的第二个条件可由第一个条件推出，因此差分格式的稳定性条件为 ， 即 。 由Lax等价定理可知，Samarskii格式也是条件收敛的。 3.3 Crank-Nicolson型隐式差分格式 前面讨论了求解对流扩散方程的两种显示格式，它们都是条件稳定的，为了放松稳定性条件，可以采用隐式格式进行求解。现在考虑Crank-Nicolson型隐式差分格式 （7） 令，，则（7）式可化为 （8） 把（8）式用矩阵的形式 = + （9） 设 ， ， ， 则有 下面讨论Crank-Nicolson型格式的截断误差和精度。该格式涉及到时间层和时间层上的，，处六个点。设是