基于Crank-Nicolson格式的输运方程的数值求解方法

魏丙涛

（文山学院 数理系，云南 文山 663000）

输运方程是数学物理方程中比较重要的方程之一，能够描述一些输运的物理过程，比如：制作半导体器件就常用到扩散法，把含有所需杂质的物质涂敷在硅片表面，把硅片放在扩散炉里，杂质就向硅片里面扩散，扩散运动的方向基本上是垂直于硅片表面而指向硅片深处[1]，这种只沿某一方向进行的扩散就可以用一个一维的输运方程来描述这一物理过程.输运方程的解法很多，一般的输运方程是可以用解析的解法求解的，对于复杂的方程、边界条件及初始条件，解析解就变的比较困难了，于是需要寻求数值解法.对于偏微分方程，数值解法一般采用有限差分法，在此介绍一种基于Crank-Nicolson格式的差分方法的数值求解.

本文第一部分将对Crank-Nicolson格式做个简单的介绍，第二部分将给出数值结果，第三部分则是本文的结论.

1 Crank-Nicolson格式简介

考虑如下的定解问题：

上述方程是一个输运方程，该方程的最基本差分格式为向前差分和向后差分格式，向前差分格式可以写为[2]：

可以将这个差分格式改写为：

向后差分的的格式可以写为[3]：

在（4）式的等式左右两边同乘θ，在（3）式的两边同乘(1-θ)，并将两式相加得到[4]：

上式为加权隐格式的差分格式，在上面的式子中，λ表示时间步长，h表示空间步长，θ是一个大于0小于1的数，当θ=0时上式成为向前差分格式，当θ=1时上式变成向后差分格式.

这个格式一般称作Crank-Nicolson格式，这个差分格式是无条件稳定的，这……