导热数值解-2015春.ppt

2 求解步骤 (2) 用 ， 值解方程式(a)，得 (3) 用 ， 值解方程式(b)，得 (4) 用 ， 值解方程式(c)，得 (5) 用 ， 值解方程式(a)，得 (1) 设初始值 ， （带括号的上角码为迭代次数序号） 只要该方程的解确实存在，则采用高斯—赛德尔迭代法必能求得解。 照此逐次计算，对于 的第 次近似值，可表示如下： 计算过程的终止 (b) 相对差值： (a) 绝对差值： 4.3 非稳态导热数值计算 一、时间—空间区域离散化 对一维非稳态导热问题，可用 坐标表示时—空区域离散化后的节点。 空间区域离散化： ——空间间隔 过程经历时间域离散化：将 划分成许多 ( 个) 小的时间区间 ，则 二、离散化方程的建立 泰勒级数展开法 控制容积法 1 显格式 导热项用 时层上的值来表示。 式中 是以 为特征长度的傅里叶数，称为网格傅里叶数。 上式表明：点 上 时刻的温度是在该点 时刻温度的基础上计及了左右两邻点温度的影响后得出的。 如果两邻点的影响保持不变，合理的情况是： 时刻 高 时刻 也高 也低 低 要满足这种合理性，必须要求 前的系数大于或等于零 网格傅里叶数 若用判别式表示，则为必须保证 2、隐格式 导热项用 时层上的值来表示。 * * * * * * * * * * 第 4 章 导热数值解法基础 航空航天热物理研究所 2015年3月 引言 工程中常遇到几何形状或边界条件复杂的导热问题，由于数学上的困难，目前还无法得出其分析解。 近三十年来，随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法得到日益广泛的应用。 导热问题的数值计算方法包括：有限差分法、有限元法、边界元法。 本章介绍物理概念明确、实施方法简便的有限差分法。 4.1 建立离散方程 一、基本思想 对物理问题进行数值求解的基本思想可以概括为： 把原来在时间、空间坐标系中连续的物理量场（如温度场），用有限个离散点上的值的集合来代替 通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值 二、求解步骤 结合二维矩形域内稳态、无内热源、常物性的导热问题，说明求解过程的各个步骤 1. 建立控制方程及定解条件 四个边界条件为第一、第二、第三类。 2. 区域离散化 步长 网格线的交点为节 (结) 点 节点所代表的小区域称为控制体 (容积) 3. 建立节点物理量的代数方程—离散方程 4. 设立迭代初场。代数方程的求解方法有两大类：直接解法、迭代解法。 采用迭代解法时，需要对被求的物理场预先假定一个解，称为初场，并在求解的过程中不断改进。 5. 求解代数方程 6. 解的分析 框图 改进初场 否 若： 、 ，则求 收敛否 解的分析 是 建立控制方程及定解条件 确定节点（区域离散化） 建立节点物理量的代数方程 设立温度场的迭代初值 求解代数方程 4.2 稳态导热的数值计算 一、区域离散化 区域的划分取决于几何上的方便、所要求的计算精度 节点↑，网格密，计算精度↑，但计算时间增加 称为节点，代表 这一单元（温度、物性）。 二、内部节点离散方程的建立 1. 泰勒级数展开法（见教课本 p88） 2. 热平衡法（控制容积法） 在它的四周有 4个控制体，用节点 和 表示 节点 实际上代表 一个小控制体 对于无内热源、稳态导热，由热平衡原理 4 个相邻控制体向内部节点 代表的控制体的导热量 = 0 如果 ，则： 三、边界节点离散方程的建立 边界上温度已知，即第一类边界条件 边界上温度未