传热学第四章导热问题数值解法基础.doc

第四章 导热问题数值解法基础 第一节 建立离散方程的方法 一、区域和时间的离散化 图4-1 二维物体中的网格 二、建立离散方程的方法 1、泰勒级数展开法 … （1） … （2） （4-2） （4-3） （4-4） （4-5） （4-6） 2、热平衡法 图4-2 二维网格单元的能量平衡 以常物性，无热源二维稳态导热为例，对节点P所代表的微元体写热平衡式， （4-7） 第二节 稳态导热问题的数值计算 一、内节点离散方程的建立 （4-8） 二、边界节点离散方程的建立 图4-3 第三类边界条件的边界节点 （4-9a） （4－9b） 节 点 方 程 式 表4-1 序号 节点特征 节点方程式 1.内部节点 2.对流边界节点 3. 对流边界外部拐角节点 4. 对流边界内部拐角节点 5.绝热边界节点 6. 曲面边界节点 7.对流换热边界条件下曲面边界上的节点2 同上图 三、节点离散方程组的求解 方程组可以写为下列形式 … … … （4-10） … 1，2，… 允许误差为止，即 ≤ ≤ （4-11） … … … 高斯-赛德尔迭代的计算过程是按一定程序循环进行的，所以它是计算机常用的计算方法之一。 第三节 非稳态导热问题的数值计算 一、显式差分格式 图4-7 一维非稳态导热问题的空间和时间划分 （1） （2） （3） （4-12） ≤ Fo≤ （4-13） （4） 今按顺序号将各节点温度编号，壁面温度为，依次的内节点温度为，，，… 。因为两侧是对称的，只计算一侧即可，按式（4）计算结果列表如下： 0 500 500 500 500 100 500 100 900 100 100 500 －300 100 100 100 100 100 100 ≥ 0 Fo ≤