传热学上海交大学导热问题的数值解法.ppt

第四章　导热问题的数值解法 上海交通大学 §4-1 导热问题数值求解的基本思想 及内部节点离散方程的建立 物 理 问 题 的 数 值 求 解 过 程 (1) 泰勒级数展开法 根据泰勒级数展开式，用节点(m,n)的温度tm,n 来表示节点(m+1,n)的温度tm+1,n (1) 泰勒级数展开法 若取上面式右边的前三项，并将两式相加 移项整理即得二阶导数的中心差分： 对于二维稳态导热问题，在直角坐标中，其导热微分方程为： 作业： §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 ? 表示形式 一维、有内热源、常物性的非稳态导热问题离散方程的建立过程 Taylor 展开法 以(m,i)点为例，首先考察非稳态项： 向前差分 向后差分 中心差分 ? ?i ?i+1 ?i-1 热平衡法 1 假设温度分布线型 温度的阶梯型分布如右图所示，即温度的分布是跳跃的，并不是连续的。 t ? 阶梯分布 首先考察左下角的分布形式，则非稳态项 扩散项： t ? 阶梯分布 显示格式 x ? i+1 i m m+1 m-1 源项： 非稳态项： 扩散项： 守恒方程： 离散方程： 以网格尺寸 为特征尺度的Fourier 数 显式格式 t ? 阶梯分布 显示格式 非稳态项： 扩散项： 源项： 守恒方程： 离散方程： x ? i+1 i m m+1 m-1 t 阶梯分布 隐式格式 显式 隐式 非稳态导热差分方程的稳定性条件： (1) 主对角占优; (2) 的系数必须大于或等于零. 两种差分格式的特点： 方程数较少时，显式差分格式计算速度快，但对时间步长和空间步长有限制，如果 和 取得不好，很有可能导致计算结果发散，这是由于 隐式差分格式则没有这类稳定性的问题 用热平衡法建立边界节点的离散方程 边界节点也有显式格式和隐式格式，下面只针对显式格式，考察一无限大平板，其右侧面为第三类边界条件，针对边界节点，建立其离散方程 稳定性条件： 内节点的稳定性条件： 所以，边界节点的稳定性条件比内节点的要苛刻 对于绝热边界条件,可令边界上的对流换热量为零 即令 于是，上面方程变为： * * §4-0 引言 求解导热问题的三种基本方法： (1) 理论分析 (2)实验 (3)数值计算 三种方法的特点 三种方法的基本求解过程 第四章　导热问题的数值解法 区域离散化 建立离散方程 确定迭代方法 设定温度初场 求解代数方程组 是否收敛 解的分析 改进初场 是 否 实现过程：以二维矩形区域为例，所谓区域离散化，就是将研究区域分解成有限数量的小区域（单元），单元的顶点（或中心点）称作节点（结点），每个节点都有自己的控制区域，称作控制体（控制容积），控制体内所有特性都是均匀的，节点的温度代表每个控制体的 1 区域离散化 温度。节点之间的距离称为空间步长。节点之间的连线称为网格线，控制容积的分界面称为界面 节点的表示方法： （m，n） 横坐标节点编号 纵坐标节点编号 x y n m (m,n) M N x y n m (m,n) M N 基本概念：节点、控制体、网格线、界面线、步长 二维矩形区域的离散化 2 离散方程的建立 (1) Taylor级数展开法； (2) 多项式拟合法； (3) 控制容积积分法； (4) 控制容积平衡法(也称为热平衡法) 用节点(m,n)的温度tm,n来表示节点(m-1,n)的温度tm-1,n (m,n) (m,n+1) (m+1,n) (m,n-1) (m-1,n) 截断误差 同理可得： (m,n) (m,n+1) (m+1,n) (m,n-1) (m-1,n) 其节点方程为： (m,n) (m,n+1) (m+1,n