4非稳态导热问题的数值解法

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