石河子大学《传热学》课件-第4章导热问题数值解法.pptVIP

非稳态导热与稳态导热的主要区别：控制方程中多一个非稳态项；温度随空间和时间变化. 能量平衡关系：网格单元不仅与相邻的网格单元之间有热量的导入或导出，网格单元本身的热力学能也随时间发生变化 温度对时间的一阶导数：向前差分、向后差分 差分方程：显式差分格式、隐式差分格式 §4-4 非稳态导热问题的数值解法 一、内部节点的显式差分格式 常物性、无内热源、一维、非稳态导热 内节点 i、在 k?? 时刻热平衡关系式： 可用泰勒级数展开法或热平衡法 内节点 i、 k?? 时刻： 用Taylor级数 展开法可以得 到同样的结果 式中 是以 为特征长度的傅里叶数，称为网格傅里叶数。 只要知道 k?? 时刻各节点温度，就可以利用上式计算出 (k+1) ?? 时刻各节点的温度 显式差分格式： (k+1) ?? 时刻各节点的温度直接利用前一个时刻（ k?? ）各节点的温度以显函数的形式表示 由初始时刻(τ=0) 的温度 算出1Δτ时刻的温度 由1Δτ的温度，算出2Δτ时刻的温度 依次进行，就可算出各时刻的温度—显式差分格式 ……………………………………………………….. 由k时刻的温度, 算出k+1时刻的温度 当计算从一个时层向下一个时层推进时，任何一个时层的解取决于前一个时层的计算结果，而每次计算时，四舍五入总是难免的 解不稳定的原因