工程热力学与传热学：9-4 导热问题的数值解法基础.pptVIP

工程热力学与传热学 传热学 第九章 导热(数值计算） * * 9—4 导热问题的数值求解方法 数值解法 分析解法 导热问题的求解方法： 对导热微分方程式在规定的边界条 件和初始条件下积分求解。 求解结果能清楚显示各种因素对温度 分布的影响，但仅适用于简单的导热 问题，同时解的形式复杂。 实验方法 数值解法的基本思想 理论基础：离散数学； 基本思想：用空间或空间与时间区域内有 限个离散点（节点）上的温度近似值， 近似代替物体内原来是连续分布的温度 分布，按照一定方法建立关于这些节点 温度的代数方程。所有节点上温度值的 集合就是导热问题的数值解。 基本方法： 有限差分法 (finite-difference) 有限元法 (finite-element) 边界元法 (boundary-element) 9-4-1. 有限差分法的基本原理 将实际的物理过程离散化，用有限 差分近似微分，用有限差商近似微商。 将导热微分方程转化为 节点温度的差分方程。 主要内容 举例 以二维，稳态，无内热源的导热为例。 1. 求解域的离散化 (discretization) 沿x方向和y方向分别以Δ x， Δ y为间隔把 物体分搁成很多个小的子区域。 步长： 网格的宽度Δ x, Δ y 步长大小的选择 依问题需要 . 控制容积 (control volum)： 以节点为中心, 边长等于Δ x, Δ y的小区域 . 节点 (nodes)： 网格线(边界线)的交点. 节点 节点温度代表以它为 中心控制容积的平均温度 . 对非稳态导热问题,还需对时间域进行离散 . 2. 节点温度差分方程的建立 数值计算过程的核心内容 . 两种方法: 泰勒级数展开法; 控制容积热平衡法 . 根据傅立叶定律: 内部节点温度差分方程 垂直纸面单位宽度 设步长Δ x= Δ y —— 二维稳态导热均匀步长时 内部节点温度差分方程。 说明 此时，物体内每一个节点的温度都等于 它周围相邻 4 个节点温度的算术平均值。 稳态时，流向任何节点的热量的总和 必等于零。 第一类边界条件 — 边界温度已知 第二类边界条件 第三类边界条件 边界节点温度 差分方程 边界节点温度差分方程 举例 如图示边界条件。 建立控制容积的热平衡方程： 设步长Δ x= Δ y 网格毕渥数： ——第三类边界条件下边界节点温度差分方程。 思考题 试确定下列边界节点温度的差分方程： 9-4-2. 节点温度差分方程组的求解方法 常用方法： 消元法; 矩阵求逆法： 迭代法。 简单迭代法； 高斯-赛得尔迭代法 . . . . . . 其中： 均为常数，且 。 举例 高斯-赛得尔迭代法 设节点温度差分方程的形式为： 将该方程组改为 t1 ，t2，… tn 显函数形式 . . . . . . 求解原则：首先假设一组变量初始值，在依 次求解上述方程组时不断采用最新的当前值更 新变量的估计值。 假设初始值： . . . . . . 节点温度的 k 次近似值： 确定： 终止条件： 连续两次迭代得出 温度ti 的最大差值： 或： 迭代运算收敛。 9-4-3. 传导的热流量