[2018年最新整理]导热数值解法.pptVIP

第四章 导热问题的数值解法基础 1. 有限差分法的基本原理 2. 节点温度差分方程组的求解方法 3. 非稳态导热问题的数值解法 * * 数值解法的基本思想: 用导热问题所涉及的空间和时间区域内有限个离散点(称为节点)的温度近似值来代替物体内实际连续的温度分布，将连续温度分布函数的求解问题转化为各节点温度值的求解问题，将导热微分方程的求解问题转化为节点温度代数方程的求解问题。 数值解法的基本内容与步骤: （1）对实际导热问题的几何、物理性质进行分析，做必要的、合理的简化，建立符合实际的物理模型。 （2）根据物理模型建立完整的数学模型，即给出导热微分方程和单值性条件。 第(1)、(2)步是导热问题所有求解方法的基础。 （3）求解域离散化：用与坐标轴平行的网络线将所涉及的空间和时间区域划分成有限个子区域，将网络线的交点作为节点, 每个节点就代表以它为中心的子区域（控制容积），节点温度就代表子区域的温度。 （4）建立节点温度代数方程组； （5）求解节点温度代数方程组，得到所有节点的温度值； （6）对计算结果进行分析，若不符合实际情况，则修正上述步骤，重复进行计算，直到结果满意为止。 目前求解导热问题常用的数值解法主要有：有限差分法、有限元法等。其中有限差分法比较成熟，应用广泛。 下面主要介绍有限差分法的基本原理。 以常物性、无内热源的二维稳态导热为例: 用有限差分近似微分，用有限差商近似微商（导数） 例如:?x?dx， ,进而将导热偏微分方程转化为节 点温度差分代数方程。 (1) 求解域的离散化 1)子区域的划分 选择网格宽度?x、?y（步长），划分子区域。步长大小根据问题的需要而定。 2)节点的选择 选择网线和网线及网线与物体边界的交点作为节点，标定节点位置，如( i, j )、( i+1,j ) 等。 (2) 节点温度差分方程的建立 两种方法：泰勒级数展开法与控制容积热平衡法。只介绍控制容积热平衡法。 控制容积热平衡法：根据节点所代表的控制容积在导热过程中的能量守恒来建立节点温度差分方程。 1) 内部节点温度差分方程 内部节点( i, j )所代表的控制容积在导热过程中的热平衡 对于垂直于纸面方向单位宽度 选择?x=?y 上式可整理为 可见，物体内每一个节点温度都等于相邻4个节点温度的算术平均值。 2) 边界节点温度差分方程 对于具有第三类边界条件的边界节点 ( i,j )所代表的控制容积，根据其热平衡 选择步长?x=?y ，将上式简化 令 称为网格毕渥数。 上式可整理为 第三类边界条件下的外拐角边界节点： 第三类边界条件下的内拐角边界节点： 绝热边界节点： 运用有限差分方法可以建立导热物体所有内部节点和边界节点温度的差分方程。求解这些差分方程构成一个线性代数方程组就可以得节点温度的数值。 线性代数方程组的求解方法有消元法、矩阵求逆法、迭代法等，这里仅简单介绍在导热的数值计算中常用的迭代法中的两种： (1) 简单迭代法 (2) 高斯-塞德尔迭代法 (1) 简单迭代法 其中aij、bi为常数，且aij?0。改写为显函数形式： 假设 (2)高斯-塞德尔迭代法 高斯-塞德尔迭代法是在简单迭代法的基础上加以改进的迭代运算方法。它与简单迭代法的主要区别是在迭代运算过程中总使用最新算出的数据。 高斯-塞德尔迭代法比简单迭代法收敛速度快 。 非稳态导热数值解法的特点： （1）非稳态导热微分方程多了非稳态项，因此单值性条件中增加了初始条件； （2）除了对空间域进行离散外，还需要对时间域进行离散； （3）利用热平衡法导出节点温度方程时需要考虑控制容积的热力学能随时间的变化； （4）由于时间和空间同时离散，在有些情况下空间步长和时间步长不能任意选择，否则会带来节点温度方程求解的稳定性问题。 以第三类边界条件下常物性、无内热源的无限大平壁的一维非稳态导热问题为例。 (1) 求解域的离散 (2) 节点温度差分方程的建立 运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。 空间步长为?x，时间步长为 ??, ?x、??大小的选择需要保证节点温度方程求解的稳定性。 表示空间节点i在k??时刻（简称k时刻）的节点温度。 1) 内部节点温度差分方程 内部节点i所代表的控制容积在k时刻的热平衡： 如果节点i的温度对时间的变化率采用向前差分，热平衡方程式可写成 令 网格傅里叶数 内部节