山西能源学院教案 授课班级能动1701-1704 授课时间 计2学时 课题（章节 及内容) 4.4非稳态导热问题的数值解法 理解几种差分形式： 教学目的 了解非稳态导热的显示和隐式格式： 和要求 了解扩散项离散时所取时间层的不同对计算带来的影响。 重点 时间空间区域的离散化： 难点 非稳态项离散的方法。 教学进程 (含课堂 教学内容：时间空间区域的离散化：非稳态导热方程的隐式格式： 教学内容、 边界节点的离散方程。 教学方法、 教学方法：讲授与练习、启发讨论、诱导式、归纳总结法。 辅助手段等) 作业布置 主要 《传热学》第四版，杨世铭，陶文铨， 参考资料 高等教育出版社，2006年8月 课后自我 本章内容以分析推导为主，要求学生具有扎实的数学功底，课堂 总结分析 上部分学生因高数基础不好而无法理解推导过程（如泰勒展开）， 老师应尽可能对相关数学知识做复习补充
山西能源学院教案 授课班级 能动 1701-1704 授课时间 计 2 学时 课题（章节 及内容） 4.4 非稳态导热问题的数值解法 教学目的 和要求 理解几种差分形式； 了解非稳态导热的显示和隐式格式； 了解扩散项离散时所取时间层的不同对计算带来的影响。 重 点 难 点 时间-空间区域的离散化； 非稳态项离散的方法。 教学进程 （含课堂 教学内容、 教学方法、 辅助手段等） 教学内容：时间-空间区域的离散化；非稳态导热方程的隐式格式； 边界节点的离散方程。 教学方法：讲授与练习、启发讨论、诱导式、归纳总结法。 作业布置 主 要 参考资料 《传热学》第四版，杨世铭，陶文铨， 高等教育出版社，2006 年 8 月 课后自我 总结分析 本章内容以分析推导为主，要求学生具有扎实的数学功底，课堂 上部分学生因高数基础不好而无法理解推导过程（如泰勒展开）， 老师应尽可能对相关数学知识做复习补充
山西能源学院教案 4-4非稳态导热问题的数值解法 由前可知：非稳态导热和稳态导热二者微分方程的区别在于控制方程中多了 一个非稳态项，其中扩散项的离散方法与稳态导热一样。 一、一维非稳态导热时间一一空间区域的离散化 1、基本概念 如图4-8所示，x为空间坐标，τ为时间坐标。 1)时间步长△t:指从一个时间层到下一个时间层的间隔A:。 2)节点(，1)一一表示空间网格线与时间网格线的交点，即 表示了时间一一空间区域中一个节点的位置，相应的记为： 2、非稳态项的离散 非稳态项的离散有三种不同的格式： 1)向前差分 Ar n-1,in,i n+l,i 2)向后差分 3)中心差分 1)向前差分 (本书主要采用) 图48一维非稳态导热时间空 间区城的离散化 将函数t在节点(n,i+1)对点(n,i)作泰勒展开，则有： +=+△x △2a21 +- -+0(B) 于是有 △x 其中0（△x)截断误差表示余项中△t的最低阶为一次。 由上式得：函数t在节点(n,i+1)对点(n,i)处一阶导数的向前差分 公式： +-4 △T (4-10) 2)向后差分 将函数t在节点(n,i-1)对点（血，i)作泰勒展开，可得的向后差分公
山西能源学院教案 4-4 非稳态导热问题的数值解法 由前可知：非稳态导热和稳态导热二者微分方程的区别在于控制方程中多了 一个非稳态项，其中扩散项的离散方法与稳态导热一样。 一、一维非稳态导热时间——空间区域的离散化 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) 作泰勒展开，可得 的向后差分公
a -- 式： 8t △T (4-11) 3)中心差分 a 的向前差分与向后差分之和，即得 的中心差分表达式： -- a 2△T (4-12) 二、一维非稳态导热微分方程的离散方法 1、泰勒级数展开法 1)一维非稳态导热微分方程中的扩散项离散与稳态导热微分方程中的方法相 同，则 a 821 对一维非稳态导热微分方程中 a 的扩散项→中心差分： 非稳态项→向前差分。 d (1)非稳态项： r采用向前差分为： +)- 8t △T (4-13) 8t a (2)稳态项： x2采用中心差分则为： 62t 鼎-249+t9 4r2 (4-14) 由此可得： 02.岛-29+绵 4x2 变形得： g+D=09+91)+1-20)9 x2 (4-15) 由此可见，只要i时层上各节点的温度已知，那么i+1时层上各节点的温 度即可算出，且不需设立方程组求解。此关系式即为显式差分格式
式： （ 4-11 ） 3 ）中心差分 的向前差分与向后差分之和，即得 的中心差分表达式： （ 4-12 ） 二、一维非稳态导热微分方程的离散方法 1 、泰勒级数展开法 1 ）一维非稳态导热微分方程中的扩散项离散与稳态导热微分方程中的方法相 同，则 对一维非稳态导热微分方程中 的扩散项 → 中心差分； 非稳态项 → 向前差分。 （ 1 ）非稳态项： 采用向前差分为： （ 4-13 ） （ 2 ）稳态项： 采用中心差分则为： （ 4-14 ） 由此可得： 变形得： （ 4-15 ） 由此可见，只要 i 时层上各节点的温度已知，那么 i+1 时层上各节点的温 度即可算出，且不需设立方程组求解。此关系式即为显式差分格式
2)显示差分与隐式差分格式 求解非稳态导热微分方程，是从己知的初始温度分布出发，根据边界条件依 次求得以后各个时间层上的温度值。 (1)显示差分格式 定义：就是指若已知1时层上各节点的温度值，根据该差分格式即可算出 (+1)时层上各内点的温度，而不必求解联立方程。即+少是前一时刻(1) n节点及相邻两节点温度的显函数。 优点：计算工作量小： 缺点：受时间及空间步长的限制。 (2)隐式差分格式 a 82: 对一维非稳态导热微分方程元“x中的扩散项在(+1)时层上采用中 心差分，非稳态项将t在节点(n,i+1)处对节点(n,i)采用向前差分， g+-9 -2g+0+8 =a2+1 得： △x △x2 (4-16) 式中，已知的是i时层上的值9，而未知量有3个，无法求解+”。 定义：就是指己知ⅰ时层上各节点的温度值9，根据差分格式不能直接 算出(1+1)时层上各节点的温度，而必须求解(i+1)时层上的一个联立方程 组，才能算出(1+1)时层各节点的温度，此种差分格式称隐式差分格式。 优点：不受时间及空间的步长影响： 缺点：计算工作量大。 综上可知：①非稳态导热微分方程中，扩散项采用中心差分，非稳态项采 用向前差分得到显式差分格式： ②非稳态导热微分方程中，扩散项采用中心差分，非稳态项采用向后差分得到 隐式差分格式。 2、热平衡法 1)优点：(1)不受网格是否均匀限制： 图49边界节点离 散方程的建立
2 ）显示差分与隐式差分格式 求解非稳态导热微分方程，是从已知的初始温度分布出发，根据边界条件依 次求得以后各个时间层上的温度值。 ⑴ 显示差分格式 定义：就是指若已知 i 时层上各节点的温度值，根据该差分格式即可算出 （ i+1 ）时层上各内点的温度，而不必求解联立方程。即 是前一时刻（ i ） n 节点及相邻两节点温度的显函数。 优点：计算工作量小； 缺点：受时间及空间步长的限制。 ⑵ 隐式差分格式 对一维非稳态导热微分方程 中的扩散项在 (i+1) 时层上采用中 心差分，非稳态项将 t 在节点（ n,i+1 ）处对节点（ n,i ）采用向前差分， 得： ( 4-16 ) 式中，已知的是 i 时层上的值 ，而未知量有 3 个，无法求解 。 定义：就是指已知 i 时层上各节点的温度值 ，根据差分格式不能直接 算出 (i+1) 时层上各节点的温度，而必须求解（ i+1 ）时层上的一个联立方程 组，才能算出 (i+1) 时层各节点的温度，此种差分格式称隐式差分格式。 优点：不受时间及空间的步长影响； 缺点：计算工作量大。 综上可知： ① 非稳态导热微分方程中，扩散项采用中心差分，非稳态项采 用向前差分得到显式差分格式； ② 非稳态导热微分方程中，扩散项采用中心差分，非稳态项采用向后差分得到 隐式差分格式。 2 、热平衡法 1 ）优点：（ 1 ）不受网格是否均匀限制；
(2)不受物体是否为常数限制。 2)求解方法 如图4-9所示，一无限大平板，右侧面受周围流体的冷却，表面 △x 传热导数为h,对于边界节点N代表了宽为2的元体。 对于该元体，根据傅立叶定律和能量守恒定律得： 型+能，x络拉 △x 2△T 变形为 -0-a+偶+ 2hAT (4-17) a△T 其中△x -一是以△x为特征长度的傅里叶数，称网络傅里叶数，记为： Fob。 h△x 一一是以△x为特征长度的毕渥数，称网络毕渥数，记为： A。 hAt Oc△x 项变形如下： AT-是△TAx=aAAx=而，B c△xc△x2A△x22 所以袋D=9（1-2D△B△-2RA)+20a鼎+2 FoA BiAI (4-18) 补充：4-17的推导过程 Ax 对于2的元体： (I)根据傅立叶定律，在i时层上，从节点N一1传导给节点N的热流量， 即从N一1传给元体 单位面积的热流量为：热。-堡型 △x (a) (2)根据牛顿冷却公式，平板右侧被冷却时，在1时层上其单位面积损失的热流 量为： g=h心-0) (b) 9=c |△x (3)在1 时层上元体宁热力学能的增量： 其中 x
（ 2 ）不受物体是否为常数限制。 2 ）求解方法 如图 4-9 所示，一无限大平板，右侧面受周围流体的冷却，表面 传热导数为 h ，对于边界节点 N 代表了宽为 的元体。 对于该元体，根据傅立叶定律和能量守恒定律得： = 变形为 （ 4-17 ） 其中 ——是以 为特征长度的傅里叶数，称网络傅里叶数，记为： 。 ——是以 为特征长度的毕渥数，称网络毕渥数，记为： 。 一项变形如下： 所以 = （ 2 2 ） +2 +2 （ 4-18 ） 补充： 4-17 的推导过程 对于 的元体： ⑴ 根据傅立叶定律，在 i 时层上，从节点 N — 1 传导给节点 N 的热流量， 即 从 N — 1 传 给 元 体 单 位 面 积 的 热 流 量 为 ： （ a ） ⑵ 根据牛顿冷却公式，平板右侧被冷却时，在 i 时层上其单位面积损失的热流 量为： （ b ） ⑶ 在 i 时层上元体 热力学能的增量： 其中