传热学：第四章 导热问题数值解法.pptVIP

第四章 热传导问题的数值解法(Numerical Simulaton of Heat Conduction) 解不稳定的原因 当计算从一个时层向下一个时层推进时，任何一个时层的解取决于前一个时层的计算结果，而每次计算时，四舍五入总是难免的 能保证舍入误差不被放大的差分格式称为稳定的差分格式，否则为不稳定的差分格式。 举例说明 一无限大平板，厚 ， ，两侧表面温度突然升高至 ，试用数值法求平板内温度变化。 将平板按步长 分成10 份，并取 结果出现了 和 。从物理意义上分析，可知 。 因此，上述计算结果是荒谬的，原因是 不满足稳定性条件，要改变这种情况，只有 也就是说，在时间与空间离散时， 和 的选取是相互制约的。 二、 内部节点的隐式差分格式 有时为了提高计算精度，需要令 。但在显格式中，令 ，必须同时令 ，且下降速率大于 。这就使得计算量大大增加。如何解决？ 采用隐格式。注意，扩散项用 时层上的值来表示。 扩散项取中心差分，非稳态项取向后差分。 本章小结 1.理解数值解法的基本思想，熟悉数值解法的基本步骤； 2.掌握有限差分法及热平衡法的原理; 3.能够根据导热问题的特点，合理地进行求解域的离散。 作业 内节点 i、 k?? 时刻： 用Taylor级数 展开法可以得 到同样的结果 式中 是以 为特征长度的傅里叶数，称为网格傅里叶数。 只要知道 k?? 时刻各节点温度，就可以利用上式计算出 (k+1) ?? 时刻各节点的温度 显式差分格式： (k+1) ?? 时刻各节点的温度直接利用前一个时刻（ k?? ）各节点的温度以显函数的形式表示 由初始时刻(τ=0) 的温度 算出1Δτ时刻的温度 由1Δτ的温度，算出2Δτ时刻的温度 依次进行，就可算出各时刻的温度—显式差分格式 ……………………………………………………….. 由k时刻的温度, 算出k+1时刻的温度 如果舍入误差在以后各时层的计算中逐渐减小，或总是一个有限值，则可求得与真值相近的结果，否则计算结果将严重偏离真值。 注意：稳定性问题与微分方程无关，它只是差分方程本身内在性质，它反映了差分方程进行计算时，误差积累（或传播）的性质。 根据稳定性分析： 必须保证每项前面的系数 一维非稳态内节点显式差分格式的稳定性条件。 一维非稳态内节点显式差分格式的稳定性条件 二维和三维非稳态导热内节点显式格式的稳定性条件： 若违反稳定性条件，则计算值波动、出现违反热力学第二定律的现象 解：初始条件： 边界条件： 因边界温度已知，故只需用中间节点计算式，此时 k 隐差分格式无条件稳定 需将所有的内节点和边界节点方程式联立求解 三 边界节点温度差分方程 第一类边界条件： 问题简单；边界节点温度给定 第二类、第三类边界条件： 根据给定的具体条件，针对边界节点所在的网格单元，写出热平衡关系式，建立边界节点温度差分方程 显式或隐式差分格式 假设边界条件为第三类： 边界节点0， k?? 时刻热平衡： 1、边界节点的显式差分格式 稳定性条件： 综合考虑内节点与边界节点：对于第三类边界条件、应用显式差分格式时，其稳定性条件： 综合考虑内节点与边界节点：对于第一、二类边界条件、 显式差分格式的 稳定性条件： 第三类边界条件、二维非稳态导热均匀网格的显式差分格式稳定性条件： 第三类边界条件、三维非稳态导热均匀网格的显式差分格式稳定性条件： 边界节点0， k?? 时刻热平衡： 2、边界节点的隐式差分格式 热力学能的变化项用向后差分 边界节点0， k?? 时刻热平衡： 隐式差分格式是无条件稳定的。 4、节点方程组的求解 所有内节点差分方程 + 边界节点差分方程 （1）应用显式差分格式时：计算简单；根据初始条 件可依次计算出各个时刻、各节点的温度值 * 求解导热问题的三种基本方法：(1) 理论分析法；(2) 数值计算 法；(3) 实验法 三种方法的基本求解过程 1) 所谓理论分析方法，就是在理论分析的基础上，直接对微分方程在给定的定解条件下进行积分，这样获得的解称之为分析解，或叫理论解； 2) 数值计算法，把原来在时间和空间连续的物理量的场，用有限个离散点上的