《传热学cai》4.7 非稳态导热数值解.pptVIP

第四章 4.7节 （8） * 4.7 非稳态导热数值解 对于不规则几何形状，非均匀边界条件，考虑物性变化等因素的影响时，数值解法是最得心应手的首选解题方法。 非稳态导热有限差分计算中对空间导数项的处理方式和稳态时一样，下面把讨论的重点放在对非稳态项的处理方法上。 第四章 4.7节 （8） * 4.7.1 差分方程：能量平衡方法 时间导数项的两种差分格式： 向前差分 向后差分 第四章 4.7节 （8） * 把非稳态项及扩散项的的差分格式代入 结论：节点(i, j)下一时刻的温度用四周相邻节点以及它自身的当前温度值来表示。 第四章 4.7节 （8） * 非稳态问题初始温度总是已知的，因此从初始温度出发，一个时层一个时层地向前推进，直到所要求的那个时刻为止，无需作任何迭代计算。称显式差分格式。 边界节点仍采用能量平衡方法： 第四章 4.7节 （8） * 取?x = ?y，整理成节点 (i, j) 下一时刻温度的显函数形式 4.7.2 显式格式的稳定性条件 显式格式无需作任何迭代计算，但在选择时间步长和空间步长时要受到严格的限制，称为稳定性条件。 对二维内节点： Fo?≤1/4 第四章 4.7节 （8） * 若违反该不等式，节点温度会出现违背能量守恒原则的增幅振荡现象。 问题：厚0.06m平壁，初始温度100℃，突然放入 0℃的流体中，平壁的导热系数等于40 W/(m?K)，表面传热系数 h =1 000 W/(m2?K)。若取 Fo? = 1, 计算结果如下： 第四章 4.7节 （8） * 4.7.2 隐式差分格式 对时间导数项采用向后差分格式，即得到隐式差分格式： 亦可以写作 第四章 4.7节 （8） * 与显式格式的基本区别：节点i 的下一时刻温度是用相邻节点的下一时刻（不再是当前时刻！）以及自身节点的当前时刻温度来表示的。于是，必须一次列出全部待求节点的差分方程并联立求解。仍然是从初始温度分布开始，逐个时层地向前推进，但每前进一个时层，就必需解一次联立方程组。 重新整理成： 下一章 * * * *