传热学第四章热传导问题的数值解法解读.pptVIP

* * * * § 4-4 非稳态导热问题的数值解法 非稳态导热数值解法的特点： （1）非稳态导热微分方程多了非稳态项，因此单值性条件中增加了初始条件； （2）除了对空间域进行离散外，还需要对时间域进行离散； （3）利用热平衡法导出节点温度方程时需要考虑控制容积的热力学能随时间的变化； （4）由于时间和空间同时离散，在有些情况下空间步长和时间步长不能任意选择，否则会带来节点温度方程求解的稳定性问题。 * * 第三类边界条件下，常物性、无内热源无限大平壁的一维非稳态导热问题为例。 1) 求解域的离散 2) 节点温度差分方程的建立 运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。 1. 一维非稳态导热的数值求解： * * 内部节点温度差分方程 内部节点n所代表的控制容积在i时刻的热平衡： 时间项一阶导数采用向前差分，扩散项采用当前（i）时层上的值来表示，热平衡方程式可写成 左侧导入热量 右侧导入热量 计算第i+1时层温度值时，用第i时层的已知值 ——内部节点温度方程的显式差分格式 * * 两点结论： (a) 任意一个内部节点n在(i+1)时刻的温度都可以由该节点及其相邻节点(n-1) 、(n+1)在i 时刻的温度由上式直接求出，不必联立求解方程组，这是显式差分格式的