第四章导热问题的数值求解

《第四章导热问题的数值求解》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、第四章导热问题的数值求解随着计算机的普及应用和性能的不断改善，以及相关的数值计算方法的发展和应用程序的开发，传热学数值计算方法作为数值求解传热问题的有效工具也得到了相应的发展，利用计算机求解传热学问题愈来愈受到人们的普遍重视，而且在计算复杂传热问题中显示出它的优越性，因而成为传热学的一个重要的分支。数值传热的相关内容也很自然地成为工程类学生学习传热学课程的不可缺少的部分。为了使学生能简要地掌握传热学数值计算的基本方法，在这里我们以导热问题为例对传热学数值计算方法做一个简单的介绍。4-1导热问题数值解概述在笫二章和笫三

2、章中我们对较为简单的导热问题，如一维、二维简单儿何形状利边界条件的稳态导热和非稳态导热、以及一些特殊导热问题，象通过肋片的导热和忽略内肌的集总导热系统，进行了分析求解。然而对于一些更为复杂的导热问题，如复杂的几何形状和边界条件以及物性变化较大的情况，分析求解往往很复杂或者根本不可能。此时求解问题的唯一途径是利用数值分析的办法获得数值解。数值求解通常是对微分方程克接进行数值积分或者把微分方程转化为一组代数方程组再求解。这里耍介绍的是后一种方法。如何实现从微分方程到代数方程的转化又可以采用不同的数学方法，如有限差分法、有

3、限元法和边界元法等。作为-•本入门的教材，这里仅向读者简要地介绍用有限差分析方法从微分方程确立代数方程的处理过程。有限差分法的基本思想是把原來在时间和空间坐标中连续变化的物理量（如温度、压力、速度和热流等），用有限个离散点上的数值集合来近似表示。有限差分的数学基础是用差商代替微商（导数），而几何意义是用函数在某区域内的平均变化率代替函数的真实变化率。在图4・1屮可以看出有限差分表示的温度场与真实温度场的区别。图中用几、T］、门…表示连续的温度场T；4丸为步长，它将区域的x方向划分为有限个数的区域，4宓JX】、力池…，

4、它们可以相等，也可以不相等。当4/相等时，力处的真实变化率a可以用平均变化率方、c或〃来表示，其中b、c和〃分别表示三种不同差分格式下的温度随吋间的变化率，即：b为向后差分格式如〜弘）-八2）;dxAxC九向前差分格式岂一©+A）—7'（UdxAxd为中心差分格式如〜八21（2）。dxAx这种差分格式也可以推广到高阶微商的情形。对于二阶微商的差分格式可以在一阶差分格式的基础上得出：d'T、一丁（州+心）一2丁（兀］）+八兀］—Ar）乔Q（AO5°采用这样的处理之后，反映温度场随时间、空间连续变化的微分方程就可以用反映

5、离散点间温度线性变化规律的代数方程来表示。当利用相应的数学办法求解这些代数方程组Z后，我们就能获得离散点上的温度值。这些温度值就可以近似表示温度场的连续的温度分布。从上面的分析不难看岀，当我们要对导热问题进行数值求解时一定耍采取三个大的步骤，即研究区域的离散化；离散点（节点）差分方程的建立；节点方程（代数方程）的求解。下面我们将导热问题的数值求解进行较为详细的讨论。4-2研究区域温度场的离散化导热问题的温度场是假设为时间和空问的连续函数，当进行数值求解时首先要做的事情是在所研究的时间和空间区域内把时间和空间分割成为有

6、限人小的小区域，尤如地球被人为地划分为不同的地域且冠以不同的名称，时间被年、月、口和时、分、秒分割。如果在所分割的每一个时间间隔和空间区域内均用同一个温度值來农示，那么原來连续变化的温度场就被一个离散的阶跃变化的温度分布所代替。这就是连续变化的温度场离散化处理的基本思路。这里我们以一个矩形长柱休的非稳态导热过程为例来讨论区域离散化问题。如果不考虑矩形长林体长度方向上的温度变化，那么它是一个二维非稳态导热问题，图4-2表示了长柱体矩形截而上区域离散化的情况。图中可见，対于给定的空间区域，在/方向上的步长为你,在），方向

7、上的步长为用它们作为空间尺度可以将矩形区域划分成纵横交错的网格，交点称为节点。然后以节点为中心，在两个节点的中心处划分界限，定出节点的控制面积，对于三维情况则为控制体积或控制容积，因而常在-•般意义上称之为节点的控制体。控制体的形状是随着坐标系的不同而改变的，这里的控制体是一个个的矩形面积。网格的步长在每一个方向上可以均匀划分，也可以不均匀的划分。因此，选用不同的步长和不同的划分方法，可以将同一区域划分出不同大小、不同数日的控制区域，以及不同数Id的节点数。显然，随着步长的不断减小，节点数目的不断增加，由节点温度表示

8、的离散的温度场就会更加接近连续的温度场，但计算工作量也会随之增加。在时间方向上离散化的步长常用/厂来表示，4厂的选取也是可大可小的，也可以随时间的进程而变化。显然，无限小的时间步反4厂亦会使得离散温度变化接近连续的温度改变，但随Z而来的是相应的计算工作量的增加。4-3温度场节点方程的建立为了得出所研究区域的节点温度，必须建立相应的节点方程。建立