导热数值解法基础-zw.pptVIP

导热数值解法基础Numerical Solutions to Heat Conduction Problems 现代科学研究的三大基本方法及其关系 空调送风方案对比 §4-1 建立离散方程的方法 1、内节点离散方程 以常物性、无内热源的二维稳态导热为例，对于物体内任意点 P(i,j)，它的离散方程为 若网格划分均匀，那么 对于物体内每个节点，可逐个写出它们的离散方程，从而得到 一组离散方程。 2、边界节点离散方程 边界节点收边界条件的制约和影响，因此上述得到的方程对边 界节点不适用。下面分别讨论不同边界条件下离散方程的建立 a.第一类边界条件 边界节点的温度给定，可以直接以数值形式写入与边界节点相 邻内节点离散方程中。 b.第二类边界条件 2、边界节点离散方程 c.第三类边界条件 已知对流换热系数和周围流体温度，这时 带入上式得 边界节点离散方程的建立方法 边界节点离散方程的建立方法 边界节点离散方程的建立方法 边界节点离散方程的建立方法 §4-2 稳态导热的数值计算 1. 时间—空间区域的离散化 1. 时间—空间区域的离散化 1. 时间—空间区域的离散化 * 赵 薇 辽宁工业大学 土木建筑工程学院 2010 理论分析 Analytical 数值模拟 Numerical 实验研究 Experimental 导热数值解法基础 分析解精确 几何形状及边界条件简单情况 实验条件复杂 费用昂贵 弥补两者的缺点 分类 有限差分法 有限元法 边界元法 A B 导热数值解法基础 把原来在空间与时间坐标中连续的物理量的场（如速度场、温度场、浓度场等），用有限个离散点（称为节点，node）上的值的集合来代替；通过一定的方法建立起这些离散点上变量之间关系的代数方程（称为离散方程，discretization equation）；求解所建立起来的代数方程以获得所求解变量的近似解（也称为数值解）。 对于导热问题：求解对象为导热物体的温度场 1. 物理问题数值解的基本思想 导热数值解法基础 建立控制方程 及定解条件 确定节点 （区域离散化） 建立节点方程 （方程离散化） 设温度场的 迭代初值 求解离散方程 是否收敛？ 解的分析 改进初场 是 否 ?Taylor级数展开法；?热平衡法；?控制容积积分法等 ?直接求解；?迭代求解 建立控制方程 及定解条件 确定节点 （区域离散化） 建立节点方程 （方程离散化） 设温度场的 迭代初值 求解离散方程 是否收敛？ 解的分析 改进初场 是 否 ?热流量分析；?热应力分析等 2. 物理问题数值求解的过程 导热数值解法基础 节点 (node) —所求未知量的位置，包括内节点和边界节点 步长 (step) — 相邻节点的距离 控制容积 (control volume) /微元体—实施守恒定律的最小几何单位 界面 (interface) —控制容积的分界位置 网格线 (grid lines) —沿坐标方向相邻节点连接成的曲线簇 元素 3. 区域离散化 导热数值解法基础 International Journal of Numerical Methods in Fluids, 1998,28: 1371-1387. 网格越密，节点越多，不连续节点温度的集合越逼近分析解。但是解题花费时间越多。 当网格足够细密以至于再进一步加密网格已对数值计算结果基本上没有影响时，所得到的数值解称为网格独立解（grid-independent solutions） 导热数值解法基础 4. 建立离散方程的方法 基本思想为用差商代替导数，用线性代数方程代替导热微分方程。建立离散方程的方法基本上有两种：泰勒级数展开法和热平衡法。热平衡法是学习的重点。 泰勒级数展开法 应用泰勒级数展开式，把导热微分方程中的各阶导数用相应的差分表达式来代替。 用节点(i,j)的温度来表示节点(i+1,j)的温度时： (1) 导热数值解法基础 4. 建立离散方程的方法 泰勒级数展开法 合并上式右边第三项及以后各项，整理得 式中， 表示二阶导数和更高阶导数之和，称为截断误差 也表明未写出的级数余项中 的最低阶数为1。 当 趋近于零， 趋近于零，此时可用 代替 截断误差小于或等于 。 （2）式是节点(i,j)一阶导数的向前差分表达式。 (2) 导热数值解法基础 4. 建立离散方程的方法 类似的，用节点(i,j)的温度来表示节点(i-1,j)的温度时： 合并上式右边第三项及以后各项，整理得 （4）式是节点(i,j)一阶导数的向后差分表