§4-3非稳态导热问题的数值解法
由前可知:非稳态导热和稳态导热二者微分方程的区别在于控制方程中多了一个非稳态项,其中扩散项的离散方法与稳态导热一样。
本节重点讨论:( 1 )非稳态项离散的方法;
( 2 )扩散项离散时所取时间层的不同对计算带来的影响。
一、一维非稳态导热时间——空间区域的离散化
1 、基本概念
如图 4-8 所示, x 为空间坐标, τ 为时间坐标。
1 )时间步长
:指从一个时间层到下一个时间层的间隔
。
2 )节点( n, i )——表示空间网格线与时间网格线的交点,即表示了时间——空间区域中一个节点的位置,相应的记为:
。
2 、非稳态项的离散
非稳态项的离散有三种不同的格式:
1 )向前差分
2 )向后差分
3 )中心差分
1 )向前差分
将函数 t 在节点( n,i+1 )对点( n,i )作泰勒展开,则有:
其中 0 (
)截断误差表示余项中
的最低阶为一次。
由上式得:函数 t 在节点( n,i+1 )对点( n,i )处一阶导数的向前差分公式:
( 4-10 ) 2 )向后差分
将函数 t 在节点( n,i-1 )对点 (n,i) 作泰勒展开,可得
的向后差分公式:
( 4-11 )
3 )中心差分
的向前差分与向后差分之和,即得
的中心差分表达式:
二、一维非稳态导热微分方程的离散方法
1 、泰勒级数展开法
1 )一维非稳态导热微分方程中的扩散项离散与稳态导热微分方程中的方法相同,则
对一维非稳态导热微分方程中
的扩散项 → 中心差分;
非稳态项 → 向前差分。
( 1 )非稳态项:
采用向前差分为:
( 4-13 )
( 2 )稳态项:
采用中心差分则为:
( 4-14 )