第4章­热传导问题的数值解法.pptVIP

4.2 内节点离散方程的建立方法 4.2.1 Taylor 级数展开法 包括Taylor级数展开法和热平衡法 4.2.2 热平衡法（热力学第一定律） 返回 两式相加得: 4.2.1 Taylor 级数展开法 对节点(m+1,n)和节点(m-1,n)分别写出t对节点(m,n)的Taylor级数展开: 首先推导温度在x方向二阶导数的代数表达式 略去截断误差，得到温度在x方向二阶导数的中心差分表达式: 同理,得温度在y方向二阶导数的中心差分表达式: 将差分表达式代入控制方程 如果 得: 则有: 在均分网格中，一、二阶导数常见的差分表达式如下表所示： 返回 4.2.2 热平衡法（热力学第一定律） ?n ?e ?w ?s 热平衡法不是在控制方程的基础上进行离散，而是直接对元体应用热力学第一定律和傅里叶定律，从而得到该节点温度的离散方程。 二维稳态常导热系数无内热源的稳态导热问题，对元体(m,n)列出能量守恒方程： 从元体西界面导入的热量为: 从元体东界面导入的热量为: 从元体南界面导入的热量为: 从元体北界面导入的热量为: ?n ?e ?w ?s 将各表达式代入对元体(m,n) 能量守恒方程得： 整理得: 所得结果与Taylor级数法结果相同 采用热平衡法建立节点的离散方程，物理概念清晰，推导过程简单，并且对于建立边界节点的离散方程也能适用，需要很好的掌握。 返回 4.3 边界节点离散方程的建立及代数方程的求解 4.3.1 边界节点离散方程的建立 4.3.2 处理不规则区域的阶梯型逼近法（不要求） 4.3.3 求解代数方程的迭代法 返回 边界节点的离散方程的形式与边界条件的类型有关 一、第一类边界条件情形 4.3.1 边界节点离散方程的建立 如果所有边界均为第一类边界条件类型，由于此时边界温度值为已知，所有内节点的离散方程组成了一个封闭的代数方程组，可以封闭求解。因此这种情形边界节点不需要离散方程。 此时边界温度值未知，需建立边界节点温度的离散方程。 设边界热流密度为qw，并且导热体内有内热源，下面采用元体能量平衡法来建立边界节点温度的离散方程。 如 二、第二类边界条件情形 1、平直边界上的节点 边界节点(m,n)代表的区域为半个普通大小元体。对该半个元体应用能量平衡（稳态情形） ，则有： 2、边界上的外部角点 边界节点D代表的区域为1/4个普通元体大小的面积。对该外部节点元体应用能量平衡 如 ，则有： 3、边界上的内部角点 边界节点F代表的区域为3/4个普通元体大小的面积。对该外部节点元体应用能量平衡 如 ，则有： 三、第三类边界条件情形 将该热流密度的表达式代入第二类边界条件中，可得第三类边界条件下边界节点的离散方程。 对于Δx=Δy的情形，有: 平直边界节点： 内部角点： 外部角点： 返回 4.3.2 处理不规则区域的阶梯型逼近法 不要求掌握 返回 4.3.3 求解代数方程的迭代法 代数方程组的求解方法分为直接解法（高斯消元法等）和迭代法。本书仅介绍迭代法中的高斯-赛德尔迭代法。 下面以简单的三元方程组为例说明该方法的步骤： （1）将方程组改写成关于t1,t2,t3的显式形式，即迭代方程。 1. 高斯-赛德尔迭代法 （2）假设一组解（即迭代初场），记为t1(0)、t2(0)、t3(0),由迭代公式逐一计算出改进值t1(1)、t2(1)、t3(1)。每次计算均用t的最新值代入。 （3）以计算所得之值作为初场，重复上述计算，直到相邻两次迭代值之差小于允许值，此时称为已经达到迭代收敛，迭代计算终止。 2. 迭代过程是否已经收敛的判据 判断迭代是否已经收敛的判据常用的有三种： 允许的相对偏差ε之值一般在10-3—10-6之间，视具体情况而定 3. 迭代过程能否收敛的判据 对于常物性导热问题所组成的差分方程组，迭代公式的选择应使每一个迭代变量的系数总是大于或等于该式中其他变量系数绝对值之和，此时用迭代法求解代数方程一定收敛。 该条件在数学上称为主对角线占优。 返回 4.4 非稳态导热问题的数值解法 4.4.1 时间-空间区域的离散化 相比稳态问题，在非稳态导热微分方程中多了非稳态项，因此需要学习非稳态项的离散方法。扩散项的离散方法与前者相同。 4.4.2 非稳态导热问题离散的显式格式 4.4.3 非稳态导热问题离散的隐式格式 4.4.4 边界节点的离散方程 4.4.5 非稳态导热显式格式离散方程组及稳定性分析 返回 以一维非稳态导热问题为例，介绍非稳态导热问题的离散化。 在空间-时间坐标系中对所研究的空间区域和时间区域进行离散 时间坐标 ? 分成I-1等份，?? 为时间步长 空间网格线与时间网格线的交点（n,i）代表了时间-空间区域中的一个节点位置tn(i) 4.4.1 时间-空间区域的离散化 空间坐标 x 分成 N-1等份， ?x 为