求解二维扩散方程的一种高精度紧致差分格式.PDFVIP

2015 5 ( ) May ，2015 年 月 四川师范大学学报 自然科学版 38 3 Journal of Sichuan Normal University (Natural Science) Vol． 38 ，No． 3 第 卷 第 期 求解二维扩散方程的一种高精度紧致差分格式 * ， ， 尹治丹 陈建华 葛永斌 ( ， 750021) 宁夏大学数学计算机学院 宁夏银川 : ， 摘要 扩散方程通常用来描述扩散现象中的物质密度的变化或者与扩散相类似的现象 针对二维扩散 ， ， (2 ， 方程提出了一种高精度紧致差分格式 该格式基于四次样条函数对空间变量进行离散 对时间导数采用 2)Padé ， ． ． 逼近 从而得到了时间和空间均为四阶精度的紧致差分格式 然后证明了该格式是无条件稳定的 ， ． 最后通过数值实验 验证方法的精确性和稳定性 : ; ;Padé ; ; 关键词 扩散方程 四次样条函数 逼近 高精度紧致格式 无条件稳定 中图分类号:O24 1． 82 文献标志码:A 文章编号:1001 － 8395 (2015)03 － 0365 － 06 doi :10． 3969 /j ． issn． 1001 － 8395． 2015． 03． 011 ， ， 扩散方程是一类重要的偏微分方程 通常用来 由于多项式样条函数具有一定光滑性 各相邻 ， 描述扩散现象中的物质密度的变化或者与扩散相 段上的多项式之间既具有某种连续性性质 又在各 ， ． 类似的现象 例如在群体遗传学中的等位基因在群 段之间保持了相对独立的局部性质 因此样条函数 、 、 、 作为一种很好的数值逼近工具已被不少学者用来 体中的扩散 气体的扩散 液体的渗透 半导体材料 ［5 － 7］ ［5］ 中的杂质或河流混合污染物浓度的变化等问题． 由 离散空间导数 ． 如H． W． Liu 等 提出了求解 于问题的精确解往往很难得到，因此寻找精确稳定