不同层显式格式及在微尺度热传导中的应用.pdfVIP

第 21 卷 第 4 期 计 　　算 　　物 　　理 Vol . 21 ,No . 4 2004 年 7 月 　CHINESE JOURNAL OF COMPUTATIONAL PHYSICS 　 Jul . , 2004 　 ( ) [文章编号] 　1001246X 2004 不同层显式格式及在微尺度热传导中的应用 1 ,2 1 1 1 1 1 唐 　晨 , 　张桂敏 , 　闫海青 , 　李文润 , 　张 　　 , 　刘 　铭 ( 1 天津大学应用物理系 , 天津 　300072 ; 2 中国科学院力学研究所非线性力学国家重点实验室 , 北京 　100080) [摘 　要] 　基于指数拟合和牛顿插值多项式构造了用于数值计算偏微分方程的不同层显式格式 , 应用不同层显 式格式可获得不同精度的数值计算结果 , 并将该算法应用于微结构热传导方程和薄膜强瞬态热传导方程中. [ 关键词] 　多层显式格式 ; 指数拟合 ; 微尺度热传导方程 ; 薄膜强瞬态热传导方程 [ 中图分类号] 　O241 [文献标识码] 　A 0 　 引言 近几年来自然科学和工程技术发展的一个重要趋势是朝微型化发展 ,其中微电子机械系统尤其取得了 巨大成功 ,其影响遍及相当广泛的领域. 由于能量传输和交换的普遍性 ,现在人们已经普遍地将注意力集中 到一些小尺度和快速的热现象中. 构造用于数值计算微尺度热传导方程新的高精度方法 ,对了解微电子材料 的热物性 、热行为有着重要的意义. ( ) 指数拟合法 exponential fitting EF 是近期提出的求解常微分方程方法 ,这种方法放弃了传统的差分而改 为积分 ,所以是高精度的方法[ 1～3 ] . 文[ 1]应用 Hermite 插值函数构造了2 阶和4 阶指数拟合算法 ,但推导出的 公式较为复杂. 显式指数拟合 RungeKutta 法由文[2 ]提出 ,但显式指数拟合 RungeKutta 法的右函数计算次数 要多于显式指数拟合多步法. 我们知道插值节点的个数不同 , 即不同次数插值多项式构造不同阶次的算法. 当需要增加插值节点时牛顿插值多项式只需在原插值多项式的后面再添加新项即可 ,所以本文应用牛顿插 值多项式来重新构造具有不同阶次指数拟合算法 ,并将该算法应用于微结构热传导方程和薄膜强瞬态热传 导方程中 ,取得了较好的效果. 1 　基于指数拟合不同层显式格式的构造 以微尺度热传导方程为例 ,说明新算法构造的原理 ,微尺度热传导方程为[4 ,5 ] 2 2 3 1 T 3 T T 2 T + 2 2 = 2 + 2 2 , 　0 ≤x ≤1 , t 0 , ( 1) 2 t 2π t x π t x 2 ( ) - π t (π) 其解析解是 T x , t = e sin x , 初边值条件是 T