步长和差分格式对一维非稳态导热数值计算的影响.pdf

ISSN 1008- 9446 承 德 石 油 高 等 专 科 学 校 学 报 第 4 卷第 1 期, 2002 年 3 月 　CN 13- 1265T E 　　　　　　Jou rn a l o f Ch engde P et ro leum Co llege　　　　　　 V o l. 4,N o. 1,M arch 2002 步长和差分格式对一维非稳态导热数值计算的影响 司广树 刘金平　杨泽亮 (华南理工大学电力学院, 广东 广州　5 10640) 摘要: 以一个一维非稳态导热问题为例, 研究了不同的步长和差分格式对计算结果的影响。研究表明: 空 间步长不能选得过大, 否则会使得网格数较少, 计算误差增加; 空间步长一定时, 计算结果随傅立叶数变 化而波动, 且波动幅度因空间步长的增加而加剧; 空间步长较大时, 从显式格式、 格式到 C rank N ico lson 全隐格式的计算精确度逐渐变差。另外, 显式格式计算结果的合理性较强依赖于傅立叶数的取值; 格式只在空间步长较大时, 才显示出这种依赖性; 而全隐格式没有此类问题。 C rank N ico lson 关键词: 非稳态导热; 数值计算; 傅立叶数; 单元数; 差分格式 + ( ) 中图分类号: T K 11 2　　文献标识码: A 　　文章编号: 2002 0 100 1003 　　非稳态导热是工程实际中常见的现象。本文以一个一维非稳态导热问题为研究对象, 采用有限差分 [ 1 ] 法, 数值研究了网格划分中的空间步长和时间步长对计算结果的影响, 并探讨了显式格式 、全隐格 [2 ] [3 ] 式 和 格式 等三种不同差分格式对计算结果准确性和稳定性的影响。该问题可看作 C rank N ico lson 对实际工程中存在的一维非稳态导热问题的简化和改造。 1　研究对象及其数学模型 本文的研究对象是如图 1 所示的大平壁。壁厚为 = 1 、热扩散率为 = 1. 14 2, 壁内的电热器 L m a m h 使其保持为稳定状态。初始温度分布为: = 100+ 400 ( ) 。时间 = 0 时突然切断电流, 并将两侧壁 t sin x 温控制在 100 ℃不变。现要研究壁内最高温度下降到低于 150 ℃时所需的时间。显然这是一个一维非稳 态导热问题, 且壁内沿x 方向的温度分布是对称的, 温度最高点恰在大平 ( ) 壁的中间界面 即 = 0. 5 处 。由于此问题中满足几何对称和热对称, 所 x m ( ) 以可以只对大平壁的左半部分 0