传热讲义2009-第4章.pptVIP

第4章 导热问题的数值解法 4－2 有限差分法的基本原理 4-3 非稳态导热问题的数值解法 一厚度为0.1m的无限大平壁，两侧均为对流换热边界条件，初始时 两侧流体温度与壁内温度一致，tf1=tf2=t0=5 ℃；已知两侧对流换热系数分别为h1=11 W/m2K、h2=23W/m2K, 壁的导热系数?=0.43W/mK，导温系数a=0.3437×10-6 m2/s。如果一侧的环境温度tf1突然升高为23℃并维持不变，计算在其它参数不变的条件下，平壁内温度分布及两侧壁面热流密度随时间的变化规律（用图形表示）。 同内部节点温度方程的显式差分格式的道理一样，上式必须满足显式差分格式的稳定性条件，即 因为 ，所以只要满足上式，自然满足内部节点温度方程显式差分格式的稳定性条件。因此上式是第三类边界条件下一维非稳态导热所有节点温度方程显式差分格式的稳定性条件。 对比 如果节点0的温度对时间的变化率采用向后差分，热平衡方程式可写成 引进网格付里叶数 和网格毕渥数 上式写成显函数的形式 边界节点温度方程的显式差分格式 边界节点温度方程的隐式差分格式 2) 节点温度差分方程的求解 （1）对显式差分格式，在满足稳定性条件的前提下，从初始温度开始,依次求出每个节点下一时刻的温度,不需要不必联立求解方程组; （2）对隐式差分格式，对每一个时刻的所有