第四章 导热问题的数值解.pdfVIP

第四章 导热问题的数值解 导热问题如遇到物体不规则的几何形状、非常物性、非线性边界条件等情况时，利用 分析求解就无能为力了，即使能得到结果，也往往会出现无穷级数等非常繁琐的形式，为了 克服导热问题解析求解过程中的困难，除了前面已提及的积分解等近似解法外，常用的方法 就是数值解法。数值求解方法中以有限差分法和有限单元法最为通用。有限差分法的原理是 采用离散化的处理方法，把导热问题的偏微分方程化为一个线性代数方程组。有限单元法是 采用对泛函数求极值的方法得到线性方程组。这两种方法都要借助于电子计算机计算，从而 获得令人满意的足够准确的近似解。其中用差分代替微分的有限差分法较简捷，也常被采用。 这里主要介绍利用有限差分法求解导热问题。 分析解法的优点在于是精确解，而数值解法则具有内在的误差，但这往往不是一个严 重的缺点，因为能够迅速地、较准确地计算出结果来，而且在许多的实际情况中，所考虑的 数学模型本身就只是物理模型的一种近似。数值解可以做到像工程数据所保证的那样准确。 随着高速电子计算机的进展，数值解法的应用范围也不断得到扩展。到目前为止，数 值解法几乎可以处理任何复杂程度的各种问题。虽然这种技术得到的是一定的时间间隔的各 离散点的值，但可利用插入功能和图形输出设备来弥补。 4— 1 离散化方法与差分格式 利用有限差分法求解导热问题的基本步骤是，首先将导热方程的定解问题进行离散化 处理，使其构成一个封闭的、线性的代数方程组，即构成差分格式，离散化处理包括对导热 物体几何区域的离散化和对导热过程进行的时间区段的离散化；然后对所构成的差分格式讨 论解的存在性和唯一性；最后，在满足存在性和唯一性的前提下求出各离散点上的温度值， 从而得到区域内的温度分布。 一、离散化处理 N 对导热物体几何上连续的区域用一系列的网格线划分成 j 若干的网域，具体的处理方法有两种： Δy i, j （1） 划分区域时的网格线的交点称为节点，用这些离散的 Δx i 节点群来代替原有的连续区域，如图 4—1 所示。 M （2 ） 把区域原有的区域用网格线划分成若干个小区域，这 图4—1 几何区域的离散化 些小区域称为“单元”或“控制容积”，单元内物理量 的数值常用内部点（也称节点）处的物理量的值来表示。这些单元与节点一一对应， 至于节点处于单元内什么位置为宜，视具体条件而定。 相邻两节点之间的距离称为步长，如对一个空间的二维区域，取直角坐标系(x , y ) ，x 方向的步长用Δx 表示，y 方向的步长用Δy 表示，节点的位置用(i,j ) 表示，习惯上用i 表 示x 方向的坐标位置，用j 表示y 方向的坐标位置，随x 的增大，i 由小到大，同样，随y 的增大，j 由小到大，T 就表示节点(i,j ) 处的温度。 i,j 对于非稳态导热问题，还必须对时间进行离散化，把连续的时间τ 用离散的τ 来代替， n 时间的步长用Δτ表示，T n 表示节点 (i, j ) 在τ 时刻的温度。 i j n , 1 在物理意义上，每个节点都集中着周围区域（尺度即为步长）的热容，或者说，区域内 连续分布的