2有限差分法及热传导数值计算课件.ppt

第二章 有限差分法及热传导的数值计算 下面先对稳态导热问题中位于计算区域内部的节点(简称内节点)介绍其离散方程的建立方法，而位于边界上的节点及非稳态导热中的非稳态项的离散将在以后讨论。 为讨论方便，把如图中的节点(m，n)及其邻点取出并放大，如图所示。 泰勒级数展开法 根据泰勒级数展开式，用节点(m,n)的温度tm,n 来表示节点(m+1,n)而温度tm+1,n 用节点(m,n)的温度tm,n来表示节点(m-1,n)的 温度tm-1,n 当计算区域中出现曲线边界或倾斜的边界时，常常用阶梯形的折线来模拟真实边界，然后再用上述方法建立起边界节点的离散方程。例如，如要用数值方法确定如图4-6a所示二维区域的形状因子，显然，根据对称性我们只要考虑四分之一的计算区域即可。图4-6a中的内圆边界可以来用图4-6b所示的阶梯形的折线边界来近似。只要网格取得足够密，这种近似处理方法仍能获得相当准确的结果。处理不规则边界的更好的方法要用到坐标变换，这里不做介绍。 2.4 非稳态导热问题的数值解法 非稳态导热与稳态导热的主要差别在于控制方程中多了一个非稳态项，而其中扩散项的离散方法与稳态导热是一样的。因此，本节讨论重点将放在非稳态项的离散以及扩散项离散时所取时间层的不同对计算带来的影响上。 1.泰勒展开法 首先以一维非稳态导热为例讨论时间—空间区域的离散化。如图4-8所示，x为空间坐标，我们将计算区域划分为(N-1)等份，得到N个空间节点；τ为时间坐标，我们将时间坐标上的计算区域划分为(I-1)等份，得到I个时间节点。从一个时间层到下一个时间层的间隔Δτ称为时间步长。空间网格线与时间网格线的交点，如(n，i)，代表了时间—空间区域中的一个节点的位置，相应的温度记为tn(i)。 非稳态项 的离散有三种不同的格式。如果将函数在节点(n，i+1)对点(n，i)作泰勒展开，可有 于是有 由式(b)可得在点(n，i)处一阶导数的一种差分表示式 , 的向前差分: 类似地，将t在点(n，i-1)对点(n，i)作泰勒展开，可得 的向后差分的表达式： 如果将t在点(n，i+1)及(n，i-1)处的展开式相加，则可得一阶导数的中心差分的表达式： 在非稳态导热问题的数值计算中，非稳态项的上述三种差分格式都有人采用，本书主要采用向前差分的格式，但也简单介绍了向后差分的格式。 至此，对于形如式(3-10)所示的一维非稳态导热方程，如扩散项取中心差分，非稳态项取向前差分，则有 此式可进一步改写为 求解非稳态导热方程就是从已知的初始温度分布出发，根据边界条件依次求得以后各个时间层上的温度值，式(4-14b)是对平板中各内点进行这种计算的公式。由该式可见，一旦i层上各节点的温度已知，可立即算出(i+1)时层上各内点的温度，而不必求解联立方程，因而式(4-14)所代表的计算格式称为显式差分格式。显式的优点是计算工作量小，缺点是对时间步长及空间步长有一定的限制，否则会出现不合理的结果。 如果把式(4-14a)中的扩散项也用(i+1)时层上的值来表示，则有 式中已知的是i时层的值tn(i)，而未知量有3个，因此不能直接由上式立即算出tn(i+1)之值，而必须求解(i+1)时层的一个联立方程才能得出(i十1)时层各节点的温度，因而式(4-15)称为隐式差分格式。从时—空坐标系中的节点(n，i+1)来看，式(4-15)的左端是非稳态项的一种向后差分。隐式格式的缺点是计算工作量大，但它对步长没有限制，不会出现解的振荡现象。 以上是将一维非稳态导热方程中的两个导数项用相应的差分表示式代替而建立差分方程的，这种方法称为泰勒展开法。 2.热平衡法 这种方法不受网格是否均分及物性是否为常数等限制，是更为一般的方法。 式中 是以Δx特征长度的傅里叶数，称为网格傅里叶数， 一项可作如下变化： 式中FoΔ及BiΔ分别为网格傅里叶数及网格毕渥数。于是式(4-16b)又可改写为 至此，我们可以把第三类边界条件下、厚度为2δ的无限大平板的数值计算问题作一归纳。由于问题的对称性，只要求解—一半厚度即可。设将计算区域等分为N-1