第4章-热传导问题的数值解法.pptVIP

第4章 热传导问题的数值解法 导热问题数值求解的基本原理 内节点离散方程的建立方法 边界节点离散方程的建立及代数方程的求解 非稳态导热问题的数值解法 4.1 导热问题数值求解基本思想 1. 基本思想 对物理问题进行数值求解的基本思想可以概括为： （2）区域离散化 如图 步长 网格线的交点为节（结）点 元体：把节点代表的小区域称为元体，又叫控制容积。 框图 4.2 内节点离散方程的建立方法 1. 区域离散化 区域的划分取决于几何上的方便、所要求的计算精度 节点↑，网格密，计算精度↑，但计算时间长 将上两式相加 4.3 边界节点离散方程的建立 及代数方程的求解 对于第一类边界条件的热传导问题，处理比较简单，因为已知边界的温度，可将其以数值的形式加入到内节点的离散方程中，组成封闭的代数方程组，直接求解。 对于第二类边界条件或第三类边界条件的热传导问题，就必须用热平衡的方法，建立边界节点的离散方程，边界节点与内节点的离散方程一起组成封闭的代数方程组，才能求解。 为了求解方便，这里将第二类边界条件及第三类边界条件合并起来考虑，用 表示边界上的热流密度或热流密度表达式。用 表示内热源强度。 迭代过程能够收敛的判据 有的差分格式的计算结果与真值十分相近。有的差分格式的计算结果严重偏离真值，甚至发生上下震荡，得不到结果。 当计算从一个时层向下一个时层推进时，任何一个时层的解取决于前一个时层的计算结果，而每次计算时，四舍五入总是难免的 能保证舍入误差不被放大的差分格式称为稳定的差分格式，否则为不稳定的差分格式。 上式表明：点n上i+1时刻的温度是在该点i时刻温度的基础上计及了左右两邻点温度的影响后得出的。 如果两邻点的影响保持不变，合理的情况是： 要满足这种合理性，必须要求 前的系数大于或等于零，即 举例说明 一无限大平板，厚 ， ，两侧表面温度突然升高至 ，试用数值法求平板内温度变化。 将平板按步长 分成10 份，并取 结果出现了 和 。从物理意义上分析，可知 。 因此，上述计算结果是荒谬的，原因是 有时为了提高计算精度，需要令 。但在显格式中，令 ，必须同时令 ，且下降速率大于 。这就使得计算量大大增加。如何解决？ 采用隐格式。注意，扩散项用 时层上的值来表示。 扩散项取中心差分，非稳态项取向后差分。 作业 习题4-9 直接解法：通过有限次运算获得代数方程精确解; 矩阵求逆、高斯消元法 迭代解法：先对要计算的场作出假设,在迭代计算过程中不断予以改进,直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 缺点：所需内存较大、方程数目多时不便、不适用于非线性问题。 迭代解法有多种：简单迭代（Jacobi迭代）、高斯-赛德尔迭代等 迭代解法求解步骤 （1）构建迭代方程 （2）迭代求解 假设迭代初场，记为t1(0)、 t2(0)、……、 tn(0)，由上式逐一计算出改进值t1(1)、 t2(1)、……、 tn(1)。 注意：每次迭代均用温度的最新值代入。 （3）以计算所得之值作为初场，重复上述计算，直到相邻两次迭代值之差小于允许值，此时称为已经达到迭代收敛，迭代计算终止。 判断迭代是否收敛的准则 k及k+1表示迭代次数； —第k次迭代得到的最大值 对于常物性导热问题所组成的差分方程组，迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值之和，此时用迭代法求解代数方程一定收敛。数学上称为主对角线占优。 在用热平衡法导出差分方程时，若每一个方程都选用导出该方程的中心节点的温度作为迭代变量，则上述条件必满足，迭代一定收敛。 4.4 非稳态导热问题的数值解法 1. 一维非稳态导热时间—空间区域的离散化 非稳态项的离散 向前差分： 向后差分： 中心差分： 2. 一维非稳态导热微分方程的离散方法 （1）泰勒级数展开法（内部节点） 中心差分 向前差分 显式差分格式 优点：计算工作量小； 缺点：受时间及空间步长的限制。 优点：①不受网格是否均匀限制； ???????? ?? ②不受物性是否为常数限制。 （2）热平衡法（边界节点） 3. 一维平板非稳态导热离散方程 4. 讨论一维导热问题显式差分格式稳定性 限制的物理意义 如果舍入误差在以后各时层的计算中逐渐减小，或总是一个有限值，则可求得与真值相近的结果，否则计算结果将严重偏离真值。 注意：稳定性问题与微分方程无关，它只是差分方程本身内在性