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

第 29 卷 第 2 期（上） 赤 峰 学 院 学 报（自 然 科 学 版 ） Vol. 29 No.2 2013 年 2 月 Journal of Chifeng University （Natural Science Edition ） Feb. 2013 基于 Crank-Nicolson 格式的输运方程的数值求解方法 魏丙涛 （文山学院 数理系， 云南 文山 663000） 摘 要： 本文将介绍一种求解输运方程的数值方法 Crank-Nicolson 格式，通过数值求解的方法得到方程的数值解并与 该方程的解析解对比，给出数值解法的稳定性以及数值方法精度. 关键词： 加权隐格式；输运方程；数值方法 中图分类号：0175.2 文献标识码：A 文章编号 ：1673-260X（2013 ）02-0012-02 输运方程是数学物理方程中比较重要的方程之一，能 在（4 ）式的等式左右两边同乘 兹 ，在（3 ）式的两边同乘 够描述一些输运的物理过程，比如：制作半导体器件就常 [4] (1-兹) ，并将两式相加得到 ： 用到扩散法， 把含有所需杂质的物质涂敷在硅片表面，把 n n-1n n n n n-1 n-1 n-1 u -u u -2u +u u -2u +u j jj -a j+1 j j-1 j+1 j j-1 =0 （5 ） 鄣 鄣 硅片放在扩散炉里，杂质就向硅片里面扩散，扩散运动的 子 兹 h2 +(1-兹) h2 [1] 姿 表 方向基本上是垂直于硅片表面而指向硅片深处 ，这种只 上式为加权隐格式的差分格式 ， 在上面的式子中， 沿某一方向进行的扩散就可以用一个一维的输运方程来 示时间步长，h 表示空间步长，兹 是一个大于 0 小于 1 的数， 描述这一物理过程 .输运方程的解法很多，一般的输运方 当 兹=0 时上式成为向前差分格式 ，