求解热传导方程的一个高精度格式.pdfVIP

４０ １ （） ２０１５ １ Vol. ４０ No. １ Journal of Southwest China Normal University （Natural Science Edition） Jan. ２０１５ DOI ：１０.１３７１８／ .cnki.xsxb.２０１５.０１.００４ j ① 华南理工大学广州学院计算机工程学院，广州 ５１０８００ ３ ６ ： 用待定系数法构造了求解热传导方程的一个高精度隐式格式．格式的截断误差达到 O （τ ＋ h ）．证明了当 ３ １８５ １ － ≤ r＜ 时，差分格式是稳定的．通过数值试验， 比较了差分格式的解和精确解的区别，说明了差分格 ８ ４０ ２ 式的有效性． ： 热传导方程；隐式差分格式；截断误差 ： O241.82 ： A ： １０００ ５４７１（２０１５）１ ００１８ ０５ 在渗流、 扩散等领域中经常会遇到求解热传导方程的问题．在一维的情形，其模型为初边值问题 ２ 抄u 抄 u ＝ a ２ ，０ ≤ x ≤ L ，０ ≤ t ≤ T ， a ＞ ０ 抄t 抄x （１） u（x ，０）＝ φ（x ），０ ≤ x ≤ L u（０ ， t）＝ u（L ， t）＝ ０ ，０ ≤ t ≤ T ［１ ２］ 对问题（１） 的求解， 有限差分法是解决此类问题的常用方法． 常见的差分格式 － 诸如古典隐格式、 Crank-Nicolson 格式和Dufort-Frankel 格式等，虽都是绝对稳定的，但它们的截断误差较低．前两者分别是 ２ ２ ２ ２ ２ ２ τ O（τ＋ h ），O（τ ＋ h ）；后者为O τ ＋ h ＋ ，当τ＝ h 时还失去了相容性．文献［３ － ９］对上述问题