求解热传导方程的Crank_Nicolson方法.pdfVIP

20 12 10 Oct． 20 12 年 月 枣庄学院学报 29 5 JOURNAL OF ZAOZHUANG UNIVERSITY Vol． 29 NO． 5 第 卷 第 期 求解热传导方程的Crank － Nicolson 方法 陶燕燕 ( ， 266061) 青岛科技大学 数理学院 山东 青岛 ［ ］ Crank － Nicolson ， 2 2 摘 要 给出了数值求解热传导方程的一种 格式 其截断误差为O ( + h )， τ 并且分析了该差分格式 的稳定性． ， ， 在最后的数值例子中 验证了该格式求解出的数值解可以很好的逼近精确解 以及当空间步长和时间步长 1 1 ， ． 同时缩小 倍时 最大误差约缩小为原来的 2 4 ［ ］ ; ;Crank － Nicolson ; 关键词 热传导方程 差分方法 格式 无条件稳定 ［ ］O175． 26 ［ ］A ［ ］1004 － 7077 (2012)05 － 0004 － 05 中图分类号 文献标识码 文章编号 0 引言 、 Crank － Nicolson 求解热传导方程的差分方法有古典显式格式 古典隐式格式和 格 ， ， ， ， 式 单纯地从并行计算层面来看 古典显式格式计算量小 便于编程 但是这种格式是条 1 (r ， r ); Crank － Nicolson 件稳定的 ≤ 其中 是步长比 而古典隐式格式和 格式是绝对稳定 2 ， ， ． 的 但是需要解整体的线性代数方程组 不能直接实现并行运算 有限差分区域分解算法 综合了显格式和隐格式的优点， ． 是一种高效实用的方法 目前已有许多关于区域分解方 法的结果， ［2］ H = mh ， 文献