第4章-导热数值解法基础.pptVIP

举例：二维、稳态、常物性、无qv 导热问题 A: 首先，写出通用微分表达式： 2. 热平衡法 * 设有一矩形平板，如图示。 a=2b。在边界x=0和y=0是 绝热的，在x=a处给出第三 类边界条件，即给定h和tf 。 而y=b处，是第一类边界条 件，即温度为已知，t=c11，c12，··· C15 。试写出各节点的离散方程。 例1 解：采用均匀网格Δy =Δx=b/2，给各未知节点编号t1，t2， ··· t10 。 按节点所在位置和题目所示边界条件，写出各节点的离散方程。 * * × = 系数矩阵 变量矩阵 常数矩阵 * 2.节点方程组的求解 前已写出所有内节点和边界节点的温度差分方程 n个未知节点温度t1、t2、t3….tn，n个代数方程式： 代数方程组的求解方法：直接解法、迭代解法 * 直接解法：通过有限次运算获得代数方程精确解; 矩阵求逆、高斯消元法 迭代解法：先对要计算的场作出假设、在迭代计算过程中不断予以改进、直到计算结果与假定值的结果相差小于允许值。称迭代计算已经收敛。 缺点：所需内存较大、方程数目多时不便、不适用于非线性问题（若物性为温度的函数，节点温度差分方程中的系数不再是常数，而是温度的函数。这些系数在计算过程中不断更新） 迭代解法有多种：简单迭代（Jacobi 雅克比迭代）、 高斯-赛德尔迭代、块迭代、交替方向迭代等 * 解：二维稳态导热问题。 对内部节点，应用P89表4-1序号1的公式 一各向同性材料的方形物体，其λ为常量。已知各边界的温度如图所示，且Δx=Δy。试用高斯-赛德尔迭代求其内部网格节点1, 2, 3, 4的温度。 例2 1 2 3 4 t =240℃ t =60℃ t =60℃ t =60℃ * 开始，假设t1(0)= t2(0) =120℃; t3(0) =t4(0) =80℃。带入上式进行迭代，第1~5次迭代值汇总如下： 例2 1 2 3 4 t =240℃ t =60℃ t =60℃ t =60℃ 82.46 82.42 127.42 127.345 3 82.50 82.495 127.495 127.49 5 82.49 82.48 127.48 127.46 4 82.345 82.19 127.19 126.875 2 81.875 81.25 126.25 125 1 80 80 120 120 0 t4/℃ t3/℃ t2/℃ t1/℃ 迭代次数 * §4-3 非稳态导热的数值计算---1D 一、显式差分格式 2. 物体温度沿x方向变化，物体在x方向被分割为n段，所取步长Δx ，其中第i个节点坐标为iΔx ，简记为i。 3. 物体中温度随时间变化，所取时间间隔步长Δτ，其中第k个节点坐标为kΔτ，简记为k。 表示物体中kΔτ时刻，iΔx位置处温度： 1. 导热微分方程式 一维、非稳态、常物性、无qv * 4. 内节点（i, k）温度对x的二阶导数， 采用中心差分表达式应为： 5. 温度对时间的一阶导数，取向前差分： 6. 将4与5带入1的非稳态导热方程 7. 内节点在计算时刻节点温度的差分方程 代入上式得： 网格Fo数 * 由初始时刻(τ=0)的温度 算出Δτ时刻的温度 由Δτ的温度, 算出2Δτ时刻的温度 依次进行，就可算出各时刻的温度—显式差分格式 对应隐式差分格式 P97 公式（4-15） …………………………………………………………………….. 由k时刻的温度,算出k+1时刻的温度 8. 分析 * 9. 稳定性判据——为保证数值计算的稳定性，必须有: 若 的系数为负，不同时刻的温度计算出现波动，见P95所举实例。 故不能任意选择Δx、Δτ，必须考虑稳定性。 的系数为+ 二、隐式差分格式 ——向后差分（自阅） * 思考题： 1. 节点的概念。 2. 向前差分、先后差分、中心差分的概念. 3. 利用能量守恒定律和傅立叶定律, 推导内节点 和边界点离散方程的基本方法。 4. 两个导热系数不同的物体紧紧贴在一起, 不计 接触热阻, 如何推导接触面上节点离散方程。 5. 显式差分方程、稳定性判据。 上海交通大学 * numerical solution 数值解 第四章 导热数值解法基础 Principles of Numerical Methods of Conduction 主要内容 §4-0 概述 §4-1 建立离散方程的方法 §4-2 稳态导热的数值计算（2D） §4-3 非稳态导热的数值计算（1D）