对流扩散方程有限差分方法.doc

对流扩散方程有限差分方法
求解对流扩散方程的差分格式有很多种，在本节中将介绍以下3种有限差分格式：中心差分格式、Samarskii格式、Crank-Nicolson型隐式差分格式。
3.1中心差分格式
时间导数用向前差商、空间导数用中心差商来逼近，那么就得到了（1）式的中心差分格式
（3）
若令，，则（3）式可改写为
（4）
从上式我们看到，在新的时间层上只包含了一个未知量，它可以由时间层上的值，，直接计算出来。因此，中心差分格式是求解对流扩散方程的显示格式。
假定是定解问题的充分光滑的解，将，，分别在处进行Taylor展开：

代入(4)式，有

显然，当，时，，即中心差分格式与定解问题是相容的。由以上的讨论也可得知，对流扩散方程的中心差分格式的截断误差为。
对于我们上面构造的差分格式，是否可以直接用于实际计算呢？也就是说，如果初始值有误差，在计算过程中误差会不会扩大传播呢？这就是接下来我们要讨论的是差分方程的稳定性问题。下面用Fourier方法来分析中心差分格式的稳定性。
令，代入到（4）式
整理得
所以该差分格式的增长因子为：

其模的平方为

由于，所以（即差分格式稳定）的充分条件为

上式可以改写为

注意到，所以上面不等式满足的条件为
，。
由此得到差分格式（3）的稳定性限制为
，。
故有结论：对流扩散方程的中心差分格式是条件稳定的。根据Lax等价定理，我们可以知道，对流扩散方程的中心差分格式是条件收敛的。
3.2Samarskii格式
设>0，先对方程(1)作扰动，得到另一个对流扩散方程
（5）
其中，当时，（5）式化为（1）式
对于（5）式，构造迎风格式
（6）
差分格式（6）称为逼近对流扩散方程的Samarskii格式。
首先推导（6）的截断误差。设是对流扩散方程（1）式的充分光滑的解
令

用Taylor级数展开有

再令

用Taylor级数展开有

由于
所以

当，时，，所以Samarskii格式与定解问题是相容的，并且其截断误差为。
现在看看Samarskii格式的稳定性。将（6）式两边同时加上，把（6）式化为

令，则上式即为：

根据中心显示格式稳定性的讨论，可以得到（6）式的稳定性条件为
，

即
，

稳定性的第二个条件等价于

而
利用不等式

所以

利用稳定性的第一个条件，有，从而可知稳定性条件的第二个条件可由第一个条件推出，因此差分格式的稳定性条件为
，
即。
由Lax等价定理可知，Samarskii格式也是条件收敛的。
3.3Crank-Nicolson型隐式差分格式
前面讨论了求解对流扩散方程的两种显示格式，它们都是条件稳定的，为了放松稳定性条件，可以采用隐式格式进行求解。现在考虑Crank-Nicolson型隐式差分格式

（7）
令，，则（7）式可化为
（8）

把（8）式用矩阵的形式

=
+（9）

设，
，
，则有

下面讨论Crank-Nicolson型格式的截断误差和精度。该格式涉及到时间层和时间层上的，，处六个点。设是定解问题的充分光滑的解，把(7)式中各的值用代替，然后将，，，，，分别在点处进行Taylor展开：

这里出现的的各阶偏导数假设都是存在而且连续的。于是（7）式的截断误差

显然，Crank-Nicolson型格式的精度是二阶的。再来看看该格式的稳定性情况，我们还是用Fourier方法来分析。
令，代入到（8）式

整理得所以Crank-Nicolson型格式的增长因子是

其模的平方

改写上式

由于及上式的分母为正，故

即，从而得出Crank-Nicolson型格式是无条件稳定的。根据Lax等价定理，Crank-Nicolson型格式也是无条件收敛的。

4、数值例子
给出如下对流扩散方程的初边值问题：

所讨论的对流扩散方程的精确解为

4.1三种差分格式的比较
在各种对流扩散问题中，有许多对流相对于扩散来说在问题中起主导作用。对流占有扩散问题的数值求解面临很多困难。因此，对流占有扩散问题的有效数值解法一直是计算数学中重要的研究内容。
取，，，，此时上面给出的就是一个对流占优扩散问题。那么，本文讨论的三种差分格式对对流占有扩散问题的求解效果是怎样的呢？现在我们就来看看这个问题。
首先，根据差分格式的稳定性条件，确定的取值范围。
（1）中心差分格式：根据稳定性条件，可知，要使中心差分格式稳定，的取值必须满足：
（2）Samarskii格式：根据稳定性条件可知，的取值必须满足：
（3）Crank-Nicolson格式：该差分格式是无条件稳定的，所以可以取任意值。
要使三种差分格式都是稳定的，不妨取。
首先，我们通过表格看看三种差分格式的数值解与准确解之间的相对误差。
表4.1时三种差分格式结果的比较

x中心