《传热学》课件_第四章.ppt

* * * * * * * * * * * * * * * * * * * 第四章 导热问题的数值解法 Numerical Methods of Heat Conduction 导热问题一般为: 上述问题的解法有以下两种: 1. 理论解(analytical method): 通过对上述方程积分求得(有限情况)。 2. 数值解(numerical method): 用某种方式把微分方程化为关于各个离散点(节点)的代数方程,通过解代数方程获得问题近似解的方法。 连续——离散(任意情况) § 4-1 导热问题数值求解的基本思想 及内节点离散方程的建立 数学描述 区域离散化 建立节点物理量的代数方程 设立迭代初场 求解代数方程组 解的分析 数值求解的基本步骤 数学描述 导热问题一般为: 无限长棱柱（如图）导热、沿高度各截面的温度分布相同，可简化为二维问题。 2. 区域离散化 有限差分法原理 finite difference 有限元法 finite element 边界元法 boundary element 有限分析法 finite analysis 网格划分 grid 节点(node): 网格线交点. 控制容积(control volume): 节点代表的区域 ,其边界位于两点之间. 界面(interface): 控制容积的边界. 均分网格: 节点编号: 从小往大排 网格划分方法: practice A 先确定节点，后定界面 practice B 先确定界面，后定节点 3. 代数方程的建立 两式相加得: （1）Taylor 级数展开法 对点(m,n)作Taylor 展开: 代入微分方程得: 对于正方形网格 则有: 同理 （2）热平衡法（热力学第一定律） ?w ?s ?n ?e 说明: 1 用此方法所得的边界方程是有O(?x2)精度 2解析解是温度（物理量）的连续函数 3数值解得出离散点上的数值 对于边界节点要根据边界条件来确定。 1) 第一类边界条件, y=b 处将边界温度直接代入即可,方程封闭。 2)对于第三类边界条件, y=0对控制体直接应用热力学第一定律 当 时,上式为: § 4-2 边界节点离散方程的建立及代数方程的求解 1. 边界上离散方程的建立 3) 对于第二类边界条件 4）不规则边界的处理 折线法 坐标变换 2. 代数方程的求解 直接求解 矩阵求逆 消元法 迭代法 Gauss—Seidel迭代 点迭代 线迭代 块迭代 直接求解法 内存大 迭代法 使用较多 线性方程组: 可以写成: 上角标为计算序号,计算时先给出ti 的初值,然后用上式进行迭代。 终止计算的方法： Gauss---Seidel 迭代 3. 例题 针肋如右图所示，碳钢 ?=43.2W/(m.K)，求其温度分布及换热量。 解： 以上是精确解，现在我们用数值方法求解： 该问题的数学描述为 节点2： 同理得节点3 节点4 用热力学第一定律，导入的热量应等于对流散出的热量，固有： 网格划分如右图： 得 三种情况的计算结果如下 温度分布 热量计算： 误差 精确解 ?=15.06 W 四节点 ? =11.94 W 21% 三节点 ? =10.52 W 30% 如取5 节点, 则? 的误差为 19% X 0 10 15 20 30 175 139.5 127.9 119.7 113.4 175 139.8 120.13 113.8 175 128.13 114.29 § 4-3 非稳