导热问题数值解法基础幻灯片.ppt

§4-1 建立离散方程的方法 导热部分小结 改写为： 隐差分格式： (k+1) ?? 时刻各节点的温度不能直接利用前一个时刻（ k?? ）各节点的温度以显函数的形式表示（含有未知的 ） 隐差分格式无条件稳定 需将所有的内节点和边界节点方程式联立求解 三、边界节点温度差分方程 第一类边界条件： 问题简单；边界节点温度给定 第二类、第三类边界条件： 根据给定的具体条件，针对边界节点所在的网格单元，写出热平衡关系式，建立边界节点温度差分方程 显式或隐式差分格式 假设边界条件为第三类： 边界节点0， k?? 时刻热平衡： 1、边界节点的显式差分格式 稳定性条件： 综合考虑内节点与边界节点：对于第三类边界条件、应用显式差分格式时，其稳定性条件： 一维非稳态内节点显式差分格式的稳定性条件 综合考虑内节点与边界节点：对于第一、二类边界条件、 显式差分格式的 稳定性条件： 第三类边界条件、二维非稳态导热均匀网格的显式差分格式稳定性条件： 第三类边界条件、三维非稳态导热均匀网格的显式差分格式稳定性条件： 边界节点0， k?? 时刻热平衡： 2、边界节点的隐式差分格式 热力学能的变化项用向后差分 边界节点0， k?? 时刻热平衡： 隐式差分格式是无条件稳定的 4、节点方程组的求解 所有内节点差分方程 + 边界节点差分方程 （1）应用显式差分格式时：计算简单；根据初始条 件可依次计算出各个时刻、各节点的温度值 * 第四章 导热问题数值解法基础 分析解法的优点：求解过程中的数学分析较严谨；求解结果以函数形式表示，能清楚地显示各种因素对温度分布的影响 但是，工程技术中遇到的许多导热问题具有复杂的形状或边界条件，无法得出其分析解 Numerical Heat Transfer（A, B） — 数值传热学期刊 数值计算方法 — 有效解决复杂问题的方法；是具有一定精度的近似方法 数值解法：有限差分法（finite-difference）、有限元法（finite-element） 、边界元法(boundary-element)、分子动力学模拟（MD） 一、区域和时间的离散化 — 有限差分法的基本原理 把物体分割为有限个离散的单元体，用有限差商代替导数，从而将微分方程转化为差分方程。通过数值计算求取各网格单元节点的温度 二维导热问题；网格线；沿x、y方向的间距为?x、?y；网格单元 节点：网格线的交点；p(i,j) p(i,j) 边界节点 （Nodal point） 每个节点温度就代表了它所在网格单元的温度 此方法求得的温度场在空间上不连续 网格越细密、节点越多，结果越接近分析解 p(i,j) 对于非稳态导热：除了在空间上把物体分割成网格单元，时间上也要分割成许多间隔 ?? 网格越细密，计算所花时间越长 —— 在时间上也是不连续的 建立有限差分离散方程的常用方法： (1) Taylor（泰勒）级数展开法；(2) 热平衡法； (3) 多项式拟合法；(4) 控制容积积分法 （Energy balance method） p(i,j) Taylor（泰勒）级数展开法和多项式拟合法偏重于从数学的角度进行推导；而热平衡法和控制容积积分法则着重于从物理的观点来分析 建立差分方程的方法 (1) Taylor（泰勒） 级数展开法★ (2) 多项式拟合法 (3) 热平衡法★ (4) 控制容积积分法 前一类方法优点：便于对离散方程进行数学特性分析缺点：变步长网格的离散方程形式复杂、导出过程的物理概念不清晰、不能保证差分方程具有守恒特性 后一类方法优点：推导过程物理概念清晰、离散方程系数具有一定物理意义、保证差分方程具有守恒特性 缺点：不便于对离散方程进行数学特性分析 1、Taylor（泰勒）级数展开法 将函数 t(x, y, z) 在某点（i+1, j）对点（i, j）作Taylor展开 二、建立差分方程的方法 二阶导数和更高阶导数 项之和 —— 截断误差 二阶导数和更高阶导数 项之和 —— 截断误差 节点（i, j）一阶导数的向前差分表达式；一阶截差公式 将函数 t(x, y, z) 在均匀网格点（i-1, j）对点（i, j）作Taylor展开： 节点（i, j）一阶导数的向后差分表达式；一阶截差公式 节点（i, j）一阶导数的中心差分表达式；二阶截差公式 二方程相减，得： 节点（i, j）一阶导数的三种表达式中，中心差分表达式的截断误差最小；尽可能采用中心差分表达式 节点（i, j）二阶导数的中心差分表达式；二阶截差公式 二方程相加，得： 在表示温度对时间的一阶导数时一般采用向前或向后差分表达式——温度对时间的中心差分表达式不稳定 节点 P(i, j) 的温度离散方程 在写出导数差分表达式后，