第五章导热问题的数值解解析.ppt

高等传热学内容 第一章 导热理论和导热微分方程 第二章 稳态导热 第三章 非稳态导热 第四章 凝固和熔化时的导热 第五章 导热问题的数值解 第六章 对流换热基本方程 第七章 层流边界层的流动与换热 第八章 槽道内层流流动与换热 第九章 湍流流动与换热 第十章 自然对流 第十一章 热辐射基础 第十二章 辐射换热计算 第十三章 复合换热 第五章 导热问题的数值解 数值求解方法是以离散数学为基础，以计算机为工具的一种求解方法。与各种分析求解的方法相比，它在应用方面表现出很大的适用性，对于处理诸如非线性、复杂几何形状、复杂边界条件的问题以及藕合的偏微分方程组等都能较好地解决。尤其是随着计算机技术的迅速发展，数值计算的精度和速度都大大提高。因此，稍微复杂一些的导热问题现在主要依靠数值法求解。 目前用于求解偏微分方程的数值方法主要有有限差分法(FDM)、有限元法(FEM) [1,2]和边界元法(BEM)。有限元法起源于固体力学和结构分析。与有限差分法相比，有限元法在整个区域内单元的划分比较随意，更易于处理不规则几何形状的问题。有限差分法涉及的数学基础及其表达式比较简单，本章将主要介绍有限差分法及其在导热中的应用。 5-1 导数的有限差分近似表达式 有限差分的数学基础是用差商代替微商，即用有限差分代替导数。若f(x)是连续函数，则它的导数定义为 (5-1-1) 在这里，df /dx称为微商(导数)，△f/△x 称为有限差商。微商是有限差商当△x趋于零时的极限。在△x没有达到零以前，△f/△x只是df/dx的近似，而两者的差值就是用差商代替微商的偏差。 利用泰勒级数可以得到一个函数的导数的有限差分近似表达式，如图5-1所示 朝前差分 (5-l-7 ) 朝后差分 (5-1-8) 5-1 导数的有限差分近似表达式 式中的O(△x)是用差商代替微商的误差，即舍去泰勒级数的高阶项引起的误差，称为截断误差。这里O 表示数量级，O(△x)表示该截断误差是与△x 同量级的小量。 一阶导数的中心差分表达式 (5-1-9) 注意中心差分格式的截断误差是(△x)2量级的。 二阶导数的中心差分格式 (5-l-10) 前面介绍的一阶导数的差分格式都只用到两个点的函数值，为了得到更高精度的差分格式，可以利用更多点的函数值。 二阶导数的朝前差分： (5-1-14) 将以上方法推广，可以得到更高精度的各种有限差分表达式。 这里以二维稳态导热为例讨论有限差分法的应用。为了进行数值求解，首先要把求解的区域离散化。有限差分法要求来用正交的网格，在直角坐标系中就是矩形网格，如图5-2 所示。网线的交点称为节点，在区域内的节点称为内节点，落在边界上的节点为边界节点。网线之间的距离称为步长。x 方向的步长记为△x，y 方向的步长记为△y。为了便于计算，需要对节点编号。若节点P的坐标为(x, y)，x＝i△x，y＝i△y，则用(i, j )表示节点的位置，节点的温度t(x，y)则记为ti,j。 数值求解的第二步是建立节点方程，即对每一个节点写出一个代数方程。建立节点方程的方法之一是从微分形式出发，也就是在微分方程或边界条件中用差商近似微商，可以得到以节点温度为未知量的代数方程。