西交大《传热学》课件第4章 导热问题的数值解法 .pptVIP

求解n阶线性方程组，需计算n+1个行列式， 如果取n=20，N=9.7×1020 每个行列式是n！个乘积的和， 每个乘积是n个数相乘，需要做n-1次乘法， 共需乘法次数(n+1)n!(n-1)， 1亿次乘法/s的计算机，需计算30.8万年 ②Gauss消元法 共需乘法次数n3/3+n2-n/3， 如果取n=20，N=3060 适用变量个数小于50 直接解法缺点：所需内存较大,方程数目多时不 便、不适用于非线性问题。 2. 迭代法 假设初场 不断更新 收敛 关键：构造迭代格式 第k个迭代层次 迭代方程 Jacobi迭代：迭代代入值总为上轮得到的值 Gauss-Seidel迭代：将本轮最新值代入 计算机内存更省，收敛速度快 是否合理？ 3. 迭代是否收敛判据 主对角占优原则 迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值的代数和，此时，用迭代法求解代数方程，一定收敛 用每个方程采用导出该方程的中心节点的温度作为迭代变量，符合对角占优原则。 4. 迭代能否收敛判据 迭代次数 t1 t2 t3 0 1 0 3.625 0 5.675 0 3.769 2 1.735 4.545 4.996 3 1.864 4.038 5.058 4 1.983 3.980 5.013 5 2.003 3.994 5.000 6 2.0001 4.000 5.000 7 2.000 4.000 5.000 利用高斯赛德尔求解方程 迭代次数 t 0 1 2 3 4 t1 0 32 522 8722 143522 t2 0 -36 -396 -3996 -3996 t3 0 -155 -3355 -61755 -1068075 迭代发散 课程名称：传热学 Heat Transfer 内 容：4-4 非稳态导热的数值解法 §4-4 非稳态导热问题的数值解法 在非稳态导热问题中，不但需要对空间区域进行离散，还需要对时间变量进行离散，接下来以一个一维非稳态导热问题为例，重点介绍对非稳态项的离散方法，以及不同离散方法对计算带来的影响等。 一、物理问题和数学描述 厚度 2? 的无限大平壁，?、a为已知常数，?=0时温度为 t0，突然将其放置于两侧温度为 t?并保持不变的流体中，两侧表面与流体之间的表面传热系数为h。 2δ h, t∞ h, t∞ 2δ h, t∞ h, t∞ 二、时空坐标的离散 0 x t n i n,i n-1,i n+1,i n,i+1 n,i-1 把【0，?】分成（N-1）等分，有N个节点，网格步长为?x； 把要计算的时间分成（I-1）等分，有I个时层，时间步长为?? 。 任一时刻任一几何位置上的温度可以表示为： 0 x t n i n,i n-1,i n+1,i n,i+1 n,i-1 内节点 1. 热平衡法 ②右端边界节点 ③左端对称 显式代数方程组 向前差分（forward difference） 2.直接用差分代替微分 ②向后差分(backward difference) 中心差分（central difference） 内节点离散方程 三、稳定性分析（stability） 显式格式的数值稳定性 系数必须反映“正”影响！！！ ⅠⅡBC时只考虑对内点的限制； ⅢBC时必须考虑对边界点的限制。 隐式格式的数值稳定性 绝对稳定！！！ 计算复杂，工作量大！ * * * * * * * * * * * * 传热学 Heat Transfer 传热学 Heat Transfer 传热学 Heat Transfer 传热学 Heat Transfer 传热学 Heat Transfer 课程名称：传热学 Heat Transfer 内 容：4-1 导热问题数值求解基本思想 第四章 导热问题的数值解法 内容简介： 本章介绍了数值求解物理问题的基本思想，以二维稳态导热问题和一维非稳态导热为例，介绍了采用有限差分法数值求解导热问题的基本过程。 学习重点： 掌握从能量守恒定律出发建立离散方程的方法 代数方程的迭代求解方法 §4-1 数值求解的基本思想 4.1.1 基本思想 分析解：对导热微分方程在定解条件下的积分求解 数值解：用求解区域上空间、时间坐标系中的离散点的温度分布代替连续的温度场，用大量的代数方程代替微分方程 离散 连续 代数方程 微分方程 分析解和数值解的比较 ①分析解 (3)分析解具有普遍性，各种情况的影响清晰可见 (1)能获得研究问题的精确解，可以为实验和数值计 算提供比较依据 (2)局限性很大，对