32 边界节点温度差分方程.ppt

第四章 导热问题数值解法基础 主要内容 建立导热问题的温度节点方程式 节点方程式的求解方法 第一节 建立离散方程的方法 1.2 建立离散方程的方法 泰勒级数展开法 二维稳态导热节点方程式 热平衡法：(优点：适应性强，物理概念明确) 第二节 稳态导热问题的数值计算 能量守恒关系式 2.3 节点方程组的求解 第三节 非稳态导热问题的数值计算 显式格式 3.2 边界节点温度差分方程 （采用热平衡法列节点方程式） 节点方程式 作业：4-3，4-7，4-9 * * 子区域 节点 微元体 边界节点 1.1 区域和时间的离散化 区域离散：用网格线将连续的区域离散为子区域，确定网格节点及微元体。 时间离散：将时间分割成许多间隔，Δτ，2 Δτ ，3 Δτ，… 网格间距称为空间步长，时间间隔称为时间步长 以二维常物性无内热源稳态导热问题为例，介绍节点方程式建立的两种基本方法。导热微分方程式为 下面将探讨如何建立均匀网格上的节点方程式 基本泰勒展开式 二阶导数中心差分 向前差分 向后差分 中心差分 一阶导数 同理，y方向二阶导数中心差分 二阶导数中心差分 控制方程 + 差分方程 P T B L R 邻节点导入P节点的热量 P节点能量守恒关系式 P节点方程式 + 2.1 内节点离散方程 (以常物性，无热源的二维稳态导热为例) 若网格的划分是均匀的， ? x = ?y 2.2 边界节点离散方程的建立??热平衡法 P T B L 邻节点传入的热量 1 , , 1 ′ - = - y x t t Q j i j i LP D D l 或 时： — 称为对流边界网格毕渥数 节点方程式 或 方法：逆矩阵法，消元法，迭代法等。 该问题涉及的方程组系数主要是稀疏矩阵，应采用迭代法。 简单迭代法的基本形式 用迭代法求解节点的离散方程组一定要在计算机上进行。为加速迭代速度，可以用高斯－赛德尔迭代法，即等式的右边及时使用新值。 内部节点方程式 P T B L R 邻节点导入P节点的热量 P节点能量守恒关系式 P节点的能量增加量 或 稳定性条件 或 时： 隐式格式 无稳定性条件 网格傅利叶数 或 时： P T B L 邻节点传入的热量 P节点能量守恒关系式 P节点的能量增加量 P节点方程式 或 时：