传热学Chap培训教案.ppt

* * 重庆大学动力工程学院 本科生专业基础课程 传 热 学 Ⅱ 李 友 荣 E-mail: liyourong@cqu.edu.cn; Tel主要内容 概 述 节点离散方程的建立 代数方程组的求解 非稳态导热问题的数值解法 第二章 导热问题的数值解法 2.1 概 述 为什么需要数值求解？ (a) 导热问题复杂，不能得到分析解 ； (b) 计算机逐渐普及，计算速度越来越快。 什么是数值求解方法？ (a) 基本思想：空间与时间区域的离散； (b) 导热微分方程 代数方程组。 常见导热问题的数值解法： (a)有限差分法; (b)有限元法; (c)边界元法; (d)有限分析法 数值求解导热问题的步骤 (1) 建立物理数学模型 控制方程 边界条件 2.1 概 述 (2) 区域离散化 用一组与坐标轴平行的网格线把求解区域划分成许多子区域 节点（结点）：网格线的交点，需要确定温度值的空间位置。节点位置标号用(m，n)表示。 空间步长：相邻两节点间的距离，记为Dx和Dy。 元体（控制容积）：以节点为中心的一个小区域，其由相邻两节点连线的中分线构成。 内节点，边界节点 2.1 概 述 (3) 建立节点物理量的代数方程 物理量的控制方程 离散方程或节点方程 对内节点（m, n)，当Dx=Dy时， (4) 设立迭代初场 求解方法 直接解法 迭代法 初 场 (5) 求解代数方程组 线性问题：系数为常数 非线性问题 ：系数需在迭代过程中不断更新 2.1 概 述 是 建立控制方程及定解条件 是否收敛？ 区 域 离 散 化 建立节点物理量代数方程 设 立 迭 代 初 场 求 解 代 数 方 程 组 改进初场 解 的 分 析 否 (6) 解的分析 计算热流量、热应力或热变形 2.1 概 述 2.2 节点离散方程的建立 建立节点离散方程方法 泰勒级数展开法 热平衡法 1. 泰勒级数展开法 泰勒级数展开式 应用于节点（m, n） 上述两式相加得 同理可得 节点（m，n）的离散方程 如果Dx=Dy，则 2.2 节点离散方程的建立 i点的中心差分 O(Dx2) i点的中心差分 O(Dx2) i点的向后差分 O(Dx) i点的向前差分 O(Dx) 备 注 截断误差 差分表达式 导数 一阶和二阶导数常用离散表达式 2.2 节点离散方程的建立 2. 热平衡法 对象：元体（控制体） 基础：能量守恒原理 通过界面导热量 ：傅里叶定律+线性温度分布 （1）内部节点 有内热源 内热源项 2.2 节点离散方程的建立 （2）边界节点 平直边界节点 Dx=Dy 2.2 节点离散方程的建立 外部角点 内部角点 Dx=Dy Dx=Dy 2.2 节点离散方程的建立 对第三类边界条件，qw=h(tf-tm,n)，则 平直边界 外部角点 内部角点 网格Biot数，常用BiD表示 2.2 节点离散方程的建立 例题1 对某扇形区域中的二维、稳态导热，边界条件如图，试用热平衡法建立节点1和7的离散方程。 解：对节点1 对节点7 2.2 节点离散方程的建立 例题2 如图为垂直放置的短圆柱体，上部绝热，下部维持恒定温度tw，侧壁与温度为tf的流体发生对流换热，试用热平衡法建立内节点（3, 3）和上部角点（4, 5）的离散方程。 解：对内节点（3, 3） Dr=Dz 导热面积变化的影响 2.2 节点离散方程的建立 对上部角点（4, 5） 2.2 节点离散方程的建立 2.3 代数方程组的求解 数学角度 解的存在和唯一性； 线性代数方程组求解方法。 物理意义 物体具有稳态温度分布的条件：全部边界上流出的热量=物体内部的发热量，满足能量守恒 线性代数方程组数值求解方法 直接法 迭代法 1. 直接解法 经有限次运算得到方程组精确解的方法； 常用方法：高斯消元法、系数矩阵求逆法； 适用性：阶数不高的方程组； 限 制：高阶代数方程组； 实 例：一维导热问题求解的“追赶法”。 2.3 代数方程组的求解 求解一维导热问题的“追赶法” 离散方程特点：每一个内节点的方程中只包含3个未知温度，边界节点的方程中只包含2个未知温度，因此，方程系数矩阵中只有主对角线及相邻的两条对角线上有非零元素。 “顺追赶” “逆追赶” 2.3 代数方程组的求解 2. 迭代法 基本思想： 构造一个无限序列，逐步逼近精确解； 经有限步迭代得到方程组的近似解； 特别适用于稀疏矩阵。 2.3 代