林贵平《传热学》非稳态导热问题的数值解法.pptVIP

第四章　导热问题的数值解法 上海交通大学 第四章　导热问题的数值解法 网格划分如图所示 整理上式： 其中 作业： 计算所用公式 温度节点 变量标志符如下： 节点的坐标变量 　节点　　　的温度 　所算出的节点温度 输入数据 输入数据 壁表面的对流换热系数 控制打印各打印各节点温度的时间间隔数 终止计算的时间 节点的编号 K时刻的节点温度 K+1时刻节点温度 时间间隔k变量 变量标志符 * * * * §4-3非稳态导热问题的数值解法 非稳态导热与稳态导热的主要区别：温度不仅随空间变化，还随时间变化，控制方程中多一个非稳态项； 非稳态项 热源项 ？ 能量平衡关系：网格单元不仅与相邻的网格单元之间有热量的导入或导出，网格单元本身的热力学能也随时间发生变化 下面我们直接用一维、有内热源、常物性的非稳态导热问题为例给出非稳态项的处理方法 一维、有内热源、常物性的非稳态导热问题 离散方程的建立过程——空间和时间的离散化 时间步长：从一个时层到下一个时层的间隔 称为时间步长 ? x ? m-1, m, m+1 M i+1 ? 0 i i-1 (m+1,i) (m-1,i) (m,i) (m,i+1) (m,i-1) ? x m-1, m, m+1 M 0 m+1 m m-1 ? 表示形式 一维、有内热源、常物性的非稳态导热问题离散方程 的建立过程——热平衡法建立内部节点的离散方程 考察(m,i)点，则扩散项可直接写出： 源项： 非稳态项 ：不同的型线导致了不同的差分格式 假设温度随空间分段线性分布 t ? 显式 隐式 C-N格式 t ? 向前差分 向后差分 中心差分 t ? 显式 隐式 (1) 向前差分 以网格尺寸 为特征尺度的Fourier 数 t ? 显式 隐式 (2) 向后差分 第3种形式(中心差分格式)感兴趣的自己推导 可以看出，对于第一个公式，一旦i层时间上个节点的温度已知，则可以立即算出(i+1)时层上个内部节点的温度，而不必求解联立方程组，因而，称之为显式差分格式； 第二个公式则必须求解第i时层的一个联立方程才能得出(i)时层各节点的温度(迭代求解)，因此，称之为隐式差分格式。 两种差分格式的优缺点: (1) 显式差分格式计算速度快，但对时间步长和空间步长有限制，如果 和 取得不好，很有可能导致计算结果发散 稳定性条件： 的系数必须大于或等于零 即： 同理，对于二维不稳定，均匀网格的显式差分格式，稳定性条件为： 若温度对时间的一阶导数采用向后差分，则 式等价写为： 将上式移项整理： 此式为隐式差分格式 隐式差分式 显式差分式 优缺点 或 比较 用热平衡法建立边界节点的节点方程 边界节点也有 显式格式 隐式格式 考察一无限大平板，其左侧面为第三类边界条件，针对边界节点，其节点方程 边界的热容项 网格毕渥数，从为特征长度 移项整理 上式与内节点稳定性条件相比更为严格，在第三类边界条件下，应采用上式作为稳定性条件在第一类边界条件下，只需采用内节点方程定性条件 此式为 的显式差分表达式： 稳定性条件 即 对于绝热边界条件,可令边界上的对流换热量为零 即 在第三类边界条件下，二维不稳定态导热均匀网格的显式差分格式，其稳定性条件为： 即令 针对图中的边界节点1应用热平衡法 也可以写出其隐式差分格式即 整理上式，得 令 移项整理 隐式差分格式无条件稳定 同样，若是绝热边界条件可令 4-10，4-15 说明： 4-15：只列出1，2，4三个节点的离散方程即可，无需化简，也不用求解 §4-4导热问题数值计算实例 稳态导热问题： 例：一矩形薄板，几何尺寸及节点布置。如图所示。薄板左侧边界给定温度为200℃，其他三个界面给定温度为50℃，求各节点的温度 沿X方向和Y方向网格划分数 左侧边界温度 右侧边界温度 顶部边界 底部边界温度 节点温度的初始假定值 控制迭代过程终止的误差 允许的最大迭代次数 开始 输入M，N，EPS，K，TTB，TLB，TRB，TBB 迭代次数IT＝0 YES NO YES NO 打印“IT ” 打印“不收敛” 打印Ti,j 停机 非稳态导热问题 例：一厚度为0.06m的无限大屏壁，初始温度为20℃，给定壁两侧的对流换热边界条件，流体温度为150℃壁表面对流换热系数 已知壁的导热系数 ，试计算2分钟后，无限大屏壁内各节点的温度。 选定壁的半厚度作为计算对象，将半壁厚度等分为10层即N=10 节点1　　 Q 热边界面 节点11　　对流边界面 令 *