扩散方程的差分解法 在研究热传导过程、扩散过程、边界层现象时，我们常常遇到抛物型方程，这类方程中最典型、最简单的就是热传导方程。热传导方程中的自变量中包括时间t，它是描述一种随时间变化的物理过程，即所谓不定常现象。这类问题的基本定解问题应是初值问题，即在初始时刻（t=0）时给定定解条件，求解t>0时的解。 本文主要运用有限差分法对一维扩散方程进行求解，并对差分解的适定性、相容性、收敛性及稳定性进行分析，同时与解析解进行对比。

1.扩散方程 一维扩散方程为： 22uutx





（1）

式中，u为因知量，为扩散系数，x为坐标，t为时间。 其定解条件如下： 初始条件： (,0)() 0xuxfxL （2） 边界条件：

12(0,)() , (,)()utftuLtft （3） 一般假定函数()fx，1()ft，2()ft满足连接条件，即1(0)(0) ff，2()(0) fLf。

2.有限差分法 有限差分法是数值计算解微分方程古老的方法之一，也是系统化地、数值地求解数学物理方法的方程。其控制方程中的导数用离散点上函数值的差商代替。 差分格式可以分为显格式和隐格式。所谓显格式是指在任一结点上因变量在新是时间层上的值可以通过之前的时间层上相邻结点变量的值显式解出来。由于这些层的变量值是已知的，当时间向前推进时，空间点上的新的变量值就只需逐点计算就行了，因此显格式计算起来比较省事。隐格式则是指任一结点上变量在新的时间层的值，不能通过之前的时间层上相邻结点的值显式解出来，它不仅与之前的时间层上的已知值有关，而且也与新时间层的相邻结点的变量值有关。因而一个差分方程常常包括几个相邻结点上的未知数，未知数的个数取决于格式的构成形式。为了解出这些未知数需要联立新的方程，而每引进一个新的方程往往又同时引进了新的未知数。因此，隐格式总是伴随着求解巨大的代数方程组。隐格式的主要缺点是计算工作量大，因而不如显格式计算得快，但这只是就时间步长一样的情况而言的。隐格式的主要优点是时间步长可以比显格式能够采用的最大步长大很多。显格式的时间步长受到稳定性条件的限制，而隐格式则几乎不受限制。

3.方程的离散 3.1 显格式 采用时间前差及第n时间层的空间中心差，得一维扩散方程的显格式解：

11122()nnnnnjjjjjuuuuutx



（4）

即 111(2)nnnnnjjjjjuuruuu （5） 式中，2()trx

3.2 全隐格式 采用时间前差及第（n+1）时间层的空间中心差，得一维扩散方程的全隐格式解：

11111122()nnnnnjjjjjuuuuutx



（6）

4.差分解的基本问题 差分解的基本问题包括：适定性、相容性、收敛性和稳定性四个方面。 4.1 适定性 在用差分方程作微分方程数值解时，首先，要求微分方程的问题是适定的。所谓适定性问题是指这一微分方程在一定的初始和边界条件下要有唯一解，并且在初始条件和边界条件稍有改变时，微分方程的解也只是稍有偏离。从数学的角度讲，若微分方程的解存在并且是唯一的，同时连续依赖于数据（初始条件、边界条件），则问题是适定的。 对于抛物线方程，这是初值问题，要求给出初始条件，并且在区域边界上给出边界条件。在本文的一维扩散问题中，除了给出t=0时的初值(,0)ux外，应在左、右两端边界，即0x及xL，各给定一个边界条件，第一类边界条件u的值，或给定第二类边界条件un的值或给定第三类边界条件u与un的组合。 由于抛物型方程具有扩散性质，它往往是随时间扩散的，u值将不断因扩散而衰减。 在本文的一维扩散方程中，给定了初始条件，同时在区域的左、右端边界给定了边界条件，满足适定性要求。 4.2 相容性 将一个偏微分方程用差分格式化为相应差分方程，当步长t和x趋近于零时，这个差分方程应当收敛于原微分方程，也就是说，相应的差分方程和微分方程之间的截断误差在任一时刻任一网格点上均应趋近于零，这样的差分方程和微分方程才是相容的。 对一维扩散方程： 220uuLutx

（7）

采用显格式差分格式，令t，xh，则 111220nnnnnjjjjjnhjuuuuuLuh

 （8）

用Taylor级数展开代入上述差分方程中，则有 222424224(),0212nnnnhjj

jj

uuuhuLuOhtxtx





（9）

当0，0h时，上式化为 22()0nnhjjuuLutx

（10）

可见，此差分格式所构成的差分方程与原来的微分方程是相容的，故该显格式为相容格式。 采用全隐格式差分格式，令t，xh，则 111111220nnnnnjjjjjnhjuuuuuLuh

 （11）

用Taylor级数展开代入上述差分方程中，则有 222(),0nnhjj

uuLuOhtx



（12）

当0，0h时，上式化为 22()0nnhjjuuLutx

（13）

可见，此差分格式所构成的差分方程与原来的微分方程是相容的，故该全隐格式为相容格式。 4.3 稳定性 差分法计算中所产生的误差（舍入误差，参数误差等）随时间衰减或不增大，则称离散格式是稳定的，反之，则是不稳定的。分析差分方程稳定性有不同的方法，如矩阵方法，谐波分析法等。下面用谐波分析法对以上两种离散格式的稳定性进行分析： （1）显格式

显格式的差分方程为 111220()nnnnnjjjjjnhjuuuuuLutx



（14）

即 111(12)()nnnnjjjjururuu （15）

其误差方程为 111(12)()nnnnjjjjrr （16）

任取一k次谐波分量 eejikxnikjhjkkAA （17） 则 1e(12)ee(122cos)nikjhikhikhnjkjArrrrkh

 （18）

误差放大因子为 122cos12(1cos)rrkhrkh （19）

要满足稳定性条件，则要求对所有的k值均有||1，须|14|1r，即102r。 因此，一维扩散方程显格式的稳定性条件为 2(),02xt

且 （20）

（2）全隐格式 全隐格式的差分方程为

111111220()nnnnnjjjjjnhjuuuuuLutx



（21）

即 11111(12)()nnnnjjjjruuruu （22）

其误差方程为 11111(12)()nnnnjjjjrr （23）

任取一k次谐波分量 1eejikxnikjhjkkAA （24）

则 111ennikhjj，111ennikhjj （25）

则误差方程为 11[(12)r(ee)][12(1cos)]ikhikhnnnjjjrrkh （26）

误差放大因子为 112(1cos)]rkh （27）

要满足稳定性条件，则要求对所有的k值均有||1。从（28）式中可以看出，当0r（即0）时，||1恒成立。因此，全隐格式是无条件稳定的。 4.4 收敛性 如果差分方程的解为ju，微分方程的解为u，若当0t，0x时，差分方程的解与微分方程的解之差

()0jjuux （28） 则称差分格式是收敛的。 拉克斯（Lax）等价定理指出，如果问题是适定的，并且差分格式满足相容性条件，那么差分格式的稳定性就是该格式收敛性的充分而必要的条件。 由4.1、4.2和4.3的分析可知，显格式和全隐格式都满足适定性、相容性及稳定性的条件，因而这两种格式满足收敛性要求。

5.差分方程的求解 5.1 初始条件及边界条件 由以上对一维扩散问题的分析，可知，求解一维扩散方程需给定初始条件及边界条件。 在本文计算中，取10Lm，1。 初始条件（0t时） 2/ 0/2(,0)() 2(1/) /2xLxLuxfxxLLxL





（29）

边界条件为 12

(0,)()0 0(,)()0utfttuLtft





（30）

其初始时刻（0t）时的u分布如图1所示，x=0m处u随时间变化情况如图2所示，x=10m处u随时间变化情况如图3所示。