第4章导热5(免费阅读).ppt

为什么要学习数值解法 数值解法适用性强：解析法求解过程中的数学分析较严谨；求解结果以函数形式表示，能清楚地显示各种因素对温度分布的影响，但是仅能针对特殊的几何形状、特殊的定解条件 数值解法的计算结果具有一定精度 (1)采用解析法需要做一些简化的假设，例如把一个鸡蛋假设成一个球、把人体假设成一个圆柱体；(2)采用工程上的计算方法（例如形状因子）很粗糙。 一个思考题 导热部分小结 改写为： 隐差分格式： (k+1) ?? 时刻各节点的温度不能直接利用前一个时刻（ k?? ）各节点的温度以显函数的形式表示（含有未知的 ） 隐差分格式无条件稳定 需将所有的内节点和边界节点方程式联立求解 三、边界节点温度差分方程 第一类边界条件： 问题简单；边界节点温度给定 第二类、第三类边界条件： 根据给定的具体条件，针对边界节点所在的网格单元，写出热平衡关系式，建立边界节点温度差分方程 显式或隐式差分格式 假设边界条件为第三类： 边界节点0， k?? 时刻热平衡： 1、边界节点的显式差分格式 稳定性条件： 综合考虑内节点与边界节点：对于第三类边界条件、应用显式差分格式时，其稳定性条件： 一维非稳态内节点显式差分格式的稳定性条件 综合考虑内节点与边界节点：对于第一、二类边界条件、 显式差分格式的 稳定性条件： 第三类边界条件、二维非稳态导热均匀网格的显式差分格式稳定性条件： 第三类边界条件、三维非稳态导热均匀网格的显式差分格式稳定性条件： 边界节点0， k?? 时刻热平衡： 2、边界节点的隐式差分格式 热力学能的变化项用向后差分 边界节点0， k?? 时刻热平衡： 隐式差分格式是无条件稳定的 4、节点方程组的求解 所有内节点差分方程 + 边界节点差分方程 （1）应用显式差分格式时：计算简单；根据初始条 件可依次计算出各个时刻、各节点的温度值 注意：稳定性条件； （2）应用隐式差分格式时：需采用迭代法 （1）应用显式差分格式时：计算简单；根据初始条 件可依次计算出各个时刻、各节点的温度值 （参见教科书及陶文铨著的《数值传热学》） §4-4 控制容积积分法简介 导出离散方程的控制容积积分法也叫有限容积法，是传热问题的数值计算中最常用的方法 主要步骤： 1、将控制方程（导热微分方 程）在所选取的控制容积 及时间间隔内对空间和时 间进行积分 2、选定未知函数及其导数对时间及空间的局部分布 曲线形式（型线或插值方法） 3、对各个项按选定的型线进行积分，并整理成关于 节点上位置值的代数方程 用控制容积积分法导出离散方程应用举例： 一维非稳态导热微分方程 1、将导热微分方程在选取的控制容积及时间间隔内 对空间和时间进行积分 2、选定未知函数（温度）及其导数对时间及空间的局部分布曲线形式（型线或插值方法） 分段线性 阶梯式 分段线性 阶梯式 3、对各个项按选定的型线进行积分，并整理成关于 节点上位置值的代数方程 分段线性 阶梯式 分段线性 阶梯式 （1）非稳态项： 取 t 随 x 变化的型线为阶梯式 即：同一控制容积中各处 t 值相同，为节点P上之值 分段线性 阶梯式 分段线性 阶梯式 （2）扩散项： 取 t 的一阶导数随时间作显式阶跃式变化 取 t 随 x 呈分段线性变化 扩散项： （3）源项： 假设源项对 x 和 ? 均呈阶梯式变化 扩散项： 源项： 非稳态项： 一维非稳态导热问题的显式差分格式 同理可导出隐式差分格式 导热微分方程式 + 单值性条件 + 求解方法 温度场 热力学第一定律 傅里叶定律 几何条件 物理条件 时间条件 边界条件 第一类 第二类 第三类 分析解法 近似解法 数值解法 分离变量法 集总参数法 有限差分法 导 热 稳态导 热 非稳态导 热 导 热 的 数值解法 例题一： 对如图所示的圆截面直肋的一维稳态、无内热源、常物性导热问题，试分别列出内节点 m 和端部节点 M 的离散方程式。已知圆截面直径为 d 。 m m-1 m+1 M-1 M-2 M ?x ?x/2 h, t? 例题二： 厚度为3mm，??=1.05 W/(m·K)的玻璃板，底部绝热，顶部暴露在温度为 30℃的空气中，表面传热系数为 h=100W/(m2·K)，宽为 1 mm，中心相继 30 mm的电热带紧附在玻璃上表面，每米长电热带发热为 20W。取图中阴影区为计算区域，不计电热带的厚度。试用数值方法确定玻璃中的温度分布。并确定相邻两电热带间自上边界每米玻璃板散失的热量。取x、y方向的网格宽度??x = ?y = 0.001m。 电热带 玻璃板 绝热 30mm 30mm 15mm 3mm h, t? i j o (1) 泰勒级数