第4章 导热问题的数值解基础.ppt

第四章 导热问题的数值解基础 Basic for Numerical Methods of Heat Conduction § 4-1 导热问题数值求解的基本思想 及内节点离散方程的建立 数值求解的基本步骤 完 § 4-4 计算机解题介绍 为了有一定的精度，网格要细，节点往往很多。这样就要解很多个方程的线性方程组。这样的方程组如用手是无法求解的，只能用计算机。 * * §4-1 数值计算的基本思想 数值求解通常是对微分方程直接进行数值积分或者把微分方程转化为一组代数方程组再进行求解。这里要介绍的是后一种方法。 如何实现从微分方程到代数方程的转化又可以采用不同的数学方法，如有限差分法、有限元法和边界元法等。这里仅简要地介绍用有限差分析方法从微分方程确立代数方程的处理过程。 有限差分法的基本思想是把原来在时间和空间坐标中连续变化的物理量（如温度、压力、速度和热流等），用有限数目的离散点上的数值集合来近似表达。有限差分的数学基础是用差商代替微商（导数），而几何意义是用函数在某区域内的平均变化率代替函数的真实变化率。 导热问题一般数学描述为: 上述问题的解法有以下两种: 1. 理论解(analytical method): 通过对上述方程积分求得(有限情况)。 2. 数值解(numerical method): 用某种方式把微分方程化为关于各个离散点(节点)的代数方程,通过解代数方程获得问题近似解的方法。 连续——离散(任意情况) 数学描述 区域离散化 建立节点物理量的代数方程 设立迭代初场 求解代数方程组 解的分析 数学描述 导热问题一般为: 无限长棱柱（如图）导热、沿高度各截面的温度分布相同，可简化为二维问题。 2. 区域离散化 有限差分法原理 finite difference 有限元法 finite element 边界元法 boundary element 有限分析法 finite analysis 网格划分 grid 节点(node): 网格线交点. 控制容积(control volume): 节点代表的区域 ,其边界位于两点之间. 界面(interface): 控制容积的边界. 均分网格: 节点编号: 从小往大排 网格划分方法: practice A 先确定节点，后定界面 practice B 先确定界面，后定节点 3. 代数方程的建立 （1）Taylor 级数展开法 对点(i,j)作Taylor 展开: 代表二阶导数和更高阶导数之和，称为截断误差 一阶向前截差公式 一阶向后截差公式 3. 代数方程的建立 两式相减: （1）Taylor 级数展开法 对点(i,j)作Taylor 展开: 一阶导数的中心差分表达式，比一阶准确 代入微分方程得: 对于正方形网格 则有: 同理 两式相加: （2）热平衡法（热力学第一定律） 说明: 1 用此方法所得的边界方程是有O(?x2)精度 2解析解是温度（物理量）的连续函数 3数值解得出离散点上的数值 § 4-2 稳态导热的数值计算----内节点离散方程的建立 以常物性、无内热源的二维稳态导热为例 和周围节点温度有关，写出每个节点的离散方程，得到方程组 求解方程组，得到各节点的温度 只适用于内节点，对边界节点不适用，因为有边界条件 热平衡法的优点 对于边界节点要根据边界条件来确定。 1) 第一类边界条件, 将边界温度直接代入即可,方程封闭。 1. 边界上离散方程的建立 2) 第二类边界条件, 设边界的热流密度为qw，则边界的热平衡： 3)对于第三类边界条件, 对控制 体直接应用热力学第一定律 其他不同条件下的边界节点的离散方程：P89, 表4-14 P90，例题4-1 2. 节点离散方程的求解—P91 直接求解 矩阵求逆 消元法 迭代法 Gauss—Seidel迭代 点迭代 线迭代 块迭代 直接求解法 内存大 迭代法 使用较多 线性方程组: 可以写成: 上角标为计算序号,计算时先给出ti 的初值,然后用上式进行迭代。 终止计算的方法： Gauss---Seidel 迭代 3. 例题 针肋如右图所示，碳钢 ?=43.2W/(m.K)，求其温度分布及换热量。 解： 以上是精确解，现在我们用数值方法求解： 该问题的数学描述为 节点2： 同理得节点3 节点4 用热力学第一定律，导入的热量应等于