第三章 非稳态导热-108

1、第3章 非 稳 态 导 热 是指温度场随时间变化的导热过程。是指温度场随时间变化的导热过程。 绝大多数非稳态导热过程都是由边界条件的变化引起。 根据温度场随时间的变化规律不同，非稳态导热分为:1.周期性非稳态导热周期性非稳态导热在周期性变化边界条件下发生的导热过程.如内燃机气缸的气体温度随热力循环发生周期性变化，汽缸壁的导热就是周期性非稳态导热。 2.非周期性非稳态导热非周期性非稳态导热 通常是在瞬间变化的边界条件下发生的导热过程.例如热处理工件的加热或冷却等。 了解和掌握1. 非稳态导热过程中温度场的变化规律；2. 换热量的计算方法. 对解决诸如热处理工艺中加热或冷却过程的优化控制等工程实际

2、问题具有重要意义 .主要介绍: :1.1. 一维非周期性非稳态导热的解析解法及求解结果;2. 求解非稳态导热问题的集总参数法。 1. 一维非稳态导热问题的解析解 第三类边界条件下大平壁、长圆柱及球体的加热或冷却是工程上常见的一维非稳态导热问题，下面重点讨论大平壁。 1) 无限大平壁冷却或加热问题的分析解简介无限大平壁冷却或加热问题的分析解简介 1. 厚度无限大平壁；2. 材料热导率、热扩散率为常数；3. 无内热源；4. 初始温度与两侧的流体相同，为 。5. 突然将两侧流体温度降低为 ，并保持不变；6. 假设平壁表面与流体间对流换热的表面传热系数h为常数。 0t t 考虑到温度场的对称性，选取坐

3、标系如图，x坐标原点位于平壁中心，因此仅需讨论半个平壁的导热问题。 很显然，这是一个一维的非稳态导热问题，其导热微分方程式为: 22xtat 初始条件：边界条件：（对称性） 0, 0tt 0,0 xtxtthxtx,以上导热微分方程式及单值性条件组成了该非稳态导热问题的数学模型。 引进过余温度:于是导热微分方程式和单值性条件变为: tt22xa 初始条件： 边界条件： tt00, 00, 0 xxhxx, 再引进无量纲温度:无量纲坐标:可将上式及单值性条件无量纲化为: 即222Xa0 xX 222Xa初始条件：边界条件： 1, 00XX00,hXX , 1v通过量纲分析可以发现，参数组均为无量

4、纲数，称为特征数特征数，习惯上也称为准则数准则数，具有特定的物理意义。 ha、2 傅里叶数傅里叶数 aaFo22 分子为从非稳态导热过程开始到 时刻的时间，分母也具有时间的量纲，可理解为温度变化波及到面积 所需要的时间。所以Fo为两个时间之比，是非稳态导热过程的无量纲时无量纲时间间。 2毕渥数毕渥数为物体内部的导热热阻与边界处的对流换热热阻之比。 hhBi1 由前式和单值性条件可知， 是 三个参数的函数， 可表示为：Xha、2 f BiFo X,o确定上式所表达的函数关系，是求解该非稳态导热问题的主要任务。 XFoBif,0采用分离变量方法可求得解析解，这里只给出求解结果： 210,2sinc

5、ossincosnFonnnnnnxxe解的函数形式为无穷级数，式中 是超越方程 、21 Bitan的根，有无穷多个，是毕渥数的函数。 无论Bi取任何值，超越方程式的根都是正的递增数列，所以从函数形式可以看出，上式是一个快速收敛的无穷级数。 由上式可以看出，无量纲过余温度确实是三个无量纲参数Bi、Fo、 的函数，与前面由无量纲导热微分方程式分析得出的结果相一致。 xX 2) 关于解析解的讨论关于解析解的讨论 (1) 傅里叶数傅里叶数Fo对温度分布的影响对温度分布的影响 计算结果表明，当傅里叶数Fo 0.2时，取级数的第一项来近似整个级数产生的误差很小，对工程计算已足够精确。 Foexx2111

6、1110coscossinsin2,将上式左、右两边取对数：式中： xm111110coscossinsin2lnln221am 为超越方程的第一个根，只与Bi有关，即只取决于：1.第三类边界条件；2.平壁的物性；3.几何尺寸；所以当平壁及其边界条件给定之后，m为一常数，与时间、地点无关。而式右边的第二项只与Bi、 有关，与时间无关。 于是上式可改为：1x ,BiCmlnx上式说明，当 ，即 时，平壁内所有各点过余温度的对数都随时间线性变化，并且变化曲线的斜率都相等，如图所示。 2 . 0 Foa22 . 0正规状况阶段正规状况阶段在此之前的非稳态导热阶段称为非正规状况阶段。 v在正规状况阶段

7、，各点的温度都按前式的规律变化，这是非稳态导热正规状况阶段的特点之一。 将前式两边对时间求导，可得： 2211am由上式可见，m的物理意义是过余温度对时间的相对变化率，单位是 ，称为冷却率冷却率（或加热加热率率）。 s1v当Fo 0.2，物体的非稳态导热进入正规状况阶段后，所有各点的冷却率或加热率m都相同，且不随时间而变化，m的数值取决于物体的物性参数、几何形状与尺寸大小以及表面传热系数，这是非稳态导热正规状况阶段的特点之二。 如果用 表示平壁中心的过余温度，则由原式可得： mFoBifeFo,cossinsin22111110m 由原式与上式之比可得：xBifx,cos10m0mv从上式可见

8、，当Fo 0.2，非稳态导热进入正规状况阶段以后，虽然都随时间而变化，但它们的比值与时间无关，只取决于毕渥数与几何位置，这是正规状况阶段的另一重要特点。 认识正规状况阶段的温度变化规律对工程计算具有重要的实际意义，因为工程技术中的非稳态导热过程绝大部分时间都处于正规状况阶段 。 有关文献已证明，当Fo 0.2时，其它形状物体的非稳态导热也进入正规状况阶段，表现出上述特点，具有前面几式所表示的温度变化规律，只是m的数值不同而已。 非稳态导热的正规状况 对无限大平板 当 取级数的首项，板中心温度， 误差小于1%1% 20aF2 . 0F0eFxx021)cos(cossinsin2),(11111

9、0eFm021111100cossinsin2)(),0()cos()(),(1xxm),()cos(cossinsin2),(111110021xBiFofxxeF3 正规热状况的实用计算方法线算图法诺谟图诺谟图三个变量，因此，需要分开来画三个变量，因此，需要分开来画以无限大平板为例，以无限大平板为例，F00.2 时，取其级数首项即可时，取其级数首项即可(1)先画先画),(0BiFofm(2) 再再 绘制其线算图绘制其线算图),()cos()(),(1xBifxxm(3) 于是，平板中任一点的温度为于是，平板中任一点的温度为00mm同理，非稳态换热过程所交换的热量也可以绘制出。同理，非稳态换

10、热过程所交换的热量也可以绘制出。 无限大平壁与周围流体之间交换的热量 在平壁内x处平行于壁面取一厚度为dx的微元薄层，在时间 内，单位面积微元薄层放出的热量等于其热力学能的变化，即: 0dxcdxttcdQ00单位面积平壁所放出的热量为: dxcdxcQ000012 令为单位面积平壁从温度冷却到稳态所放出的热量，于是： 上式也同样被绘制成线算图。 002 cQ Fo,BifecossinsinQQFo211112112021将式代入上式，得: dxexcQFo011111021coscossinsin2124 4 分析解应用范围的推广和讨论分析解应用范围的推广和讨论（1）分析解应用范围的三点推

11、广。 对物体冷却也可适用。 对一侧绝热、另一侧为第三类边界条件的平板也可适用。 当对流换热系数趋于无穷大时，固体的表面温度就趋近于流体的温度，因而Bi时就是物体表面温度突然变化后保持不变的第一类边界条件的解。（2）Fo及Bi对温度场的影响 随着Fo（）数的增加，物体中各点的过余温度减少。 Bi 数的影响可以从两个方面说明： 在相同的Fo数条件下，Bi 数越大，m/0的值越小。Bi 数越大，意味着表面上的换热条件越强，导致物体中心温度越迅速接近周围介质的温度。在极限情况下，Bi，这相当于在过程开始瞬间物体表面就达到了周围介质的温度，物体中心温度的变化当然也最迅速。相当于壁温保持恒定的第一类边界条

12、件。 另一方面，Bi 数的大小还决定物体中温度的扯平程度。当Bi 0.1时，截面上的过余温度差值已小于5，可忽略内阻，用集总参数法。 由此可见：介质温度恒定的第三类边界条件下的分析解，在Bi时转化为第一类边界条件下的解，而在Bi0时与集总参数法的解相同。 对于:1.温度仅沿半径方向变化的圆柱体（如可近似按无限长圆柱处理的长圆柱或两端绝热的圆柱体）；2.球体在第三类边界条件下的一维非稳态导热问题； 分别在柱坐标系和球坐标系下进行分析，可求得温度分布的分析解，也是快速收敛的无穷级数，并且是Bi、Fo和 的函数：RrFoBif,0r R式中R为圆柱或球体的半径， 为圆柱或球体的初始过余温度。BihR

13、FoaR、20v当 时，无限长圆柱和球体的非稳态导热过程也都进入正规状况阶段，分析解可近似取无穷级数的第一项，近似结果也被绘制成了线算图。 Fo 02 .3. 集集 总总 参参 数数 法法 当 时，物体内部的导热热阻远小于其表面的对流换热热阻，可以忽略，物体内部各点的温度在任一时刻都趋于均匀，物体的温度只是时间时间的函数，与坐标无关。 对于这种情况下的非稳态导热问题，只须求出温度随时间的变化规律以及在温度变化过程中物体放出或吸收的热量。 Bi 01 .这种忽略物体内部导热热阻的简化分析方法称为集总参数集总参数法法。 Vc,Ath,一个任意形状的物体，如图所示：1.体积为V，表面面积为A；2.密

14、度、比热容c及热导率为常数；3.无内热源，初始温度为 ；4.突然将该物体放入温度恒定为 的流体之中，5.物体表面和流体之间对流换热的表面传热系数h为常数； 需要确定该物体在冷却过程中温度随时间的变化规律以及放出的热量。 tt00tt 假设该问题满足 的条件， 根据能量守恒，单位时间物体热力学能的变化量应该等于物体表面与流体之间的对流换热量，即 Bi 01 .tthAddtcV引进过余温度 ，上式可改写为 （*）tthAddcV初始条件为通过分离变量，(*)式可改写为： tt,000dcVhAd 将上式积分可得 : 00dcVhAdcVhAexpecVhA0式中:令 具有长度的量纲，称为物体的特

15、征长度特征长度 2AVcAVhcVhAV AlVVFoBilahllchlcVhA22得 注意，式中毕渥数与傅里叶数的下角标V 表示以 为特征长度： 1.对于厚度为 的无限大平壁，2.对于半径为R的圆柱，3.对于半径为R的圆球， VVFoBiFoBiexpeVV0lV A2llR12lR13v前面介绍的分析解及诺谟图中：1.厚度为 的无限大平壁的特征长度为 ， 与集总参数法分析结果中的相同，2.但圆柱和圆球的特征长度都为半径R，即 与 、 不同。 2hBi 2aFo B ihR2RaFoVBiVFo分析结果表明，对于形状如平板、柱体或球这样的物体，只要满足： 物体内各点过余温度之间的偏差小于5

16、%，就可以使用集总参数法计算。M是与物体形状有关的无量纲数。1.对于无限大平板， 2.对于无限长圆柱， 3.对于球， MBiV1 . 01M2/1M3/1M 当 时，物体的过余温度按指数函数规律下降，随着温差的减小，下降的速度越来越缓慢。 式中指数部分中的 具有时间的量纲，令 称为时间常数时间常数，单位是s。当 时:即物体的过余温度达到初始过余温度的36.8% 。 Bi 01 .hAcVhAcVcc%.e836368010v这说明，时间常数反映物体对周围环境温度变化响应的快慢，时间常数越小，物体的温度变化越快，越迅速地接近周围流体的温度。 00 . 1368. 0c由式可见，影响时间常数的主要

17、因素是：1.物体的热容量；2.物体表面的对流换热条件。物体的热容量愈小，表面的对流换热愈强，物体的时间常数愈小。v利用热电偶测量流体温度，时间常数越小，热电偶越能迅速地反映被测流体的温度变化，所以，热电偶端部的接点总是做得很小。 hAcVc如果1.几种不同形状的物体用同一种材料制作;2.和周围流体之间的表面传热系数也都相同;3.都满足 的条件. 单位体积的表面面积 越大的物体，时间常数越小. 在初始温度相同的情况下放在温度相同的流体中被冷却（或加热）的速度越快。 A VBi 01 .v用同一种材料制成的体积相同的圆球、长度等于直径的圆柱与正方体，三者的表面面积之比为：A圆球 : A圆柱 : A

18、正方体 = 1 : 1.146 : 1.242正方体的表面面积最大，时间常数最小 v直径为2R的球体、长度等于直径2R的圆柱体与边长为2R的正方体相比，三者单位体积的表面积相同，时间常数相同，在相同条件下的冷却（或加热）速度也相同。 时间内物体和周围环境之间交换的热量 000cVttcVQVVFoBiecVcV11000令 ，表示物体温度从 变化到周围流体温度 所放出或吸收的总热量，上式可改写成无量纲形式： 00cVQ t0tVVFoBieQQ10既适用于物体被加热的情况，也适用于物体被冷却的情况。 作业：一块厚200mm的大钢板，钢材的密度为 kg/m3，比热容为 J/(kgK)，导热系数为

19、43.2W/(mK)，钢板的初始温度为20，放入1000的加热炉中加热，表面传热系数为 W/(m2K)。试求加热40分钟时距离钢板中心50mm处的温度。7790470pc 300h导热问题的数值解法基础 分析解法的主要优点是求解过程所依据的数学分析比较严谨, 物理概念和逻辑推理比较清晰, 求解结果以函数的形式表示, 能清楚地显示各种因素对温度分布的影响。但是只有少数几何形状和边界条件都比较简单的导热问题才能精确地分析求解，对于工程上绝大多数稍复杂些的导热问题，分析解法无能为力。 数值解法的基本思想是: 用导热问题所涉及的空间和时间区域内有限个离散点(称为节点)的温度近似值来代替物体内实际连续的

20、温度分布，将连续温度分布函数的求解问题转化为各节点温度值的求解问题，将导热微分方程的求解问题转化为节点温度代数方程的求解问题。因此，求解域的离散化、节点温度代数方程组的建立与求解是数值解法的主要内容。数值解法求解导热问题的基本步骤如下： (1) 对实际导热问题的几何、物理性质进行分析，做必要的、合理的简化，建立符合实际的物理模型； (2) 根据物理模型建立完整的数学模型，即给出导热微分方程（即导热控制方程）和单值性条件； (3) 求解域离散化：将导热问题所涉及的空间和时间区域按一定的要求划分成有限个子区域，将子区域的顶点作为需要确定其温度值的空间点或时间点（即节点）, 每个节点就代表以它为中心

21、的子区域，节点温度就代表子区域的温度； (4) 建立节点温度代数方程组； (5) 求解节点温度代数方程组，得到所有节点的温度值； (6) 对计算结果进行分析，若计算结果不符合实际情况，则检查上述计算步骤，修正不合理之处，重复进行计算，直到结果满意为止。有限差分法有限差分法的基本原理就是用有限差分近似微分，用有限差商近似微商（导数），例如 ，进而将导热偏微分方程转化为节点温度差分代数方程。tdtxdx 以二维方程为例，中间和边界处差分代数方程的建立节点温度差分方程组的求解方法简单迭代法：设节点温度差分方程的形式为：其中 为常数，且 。11112 211121122 2222112 2j jn nj jn nnnnj jnn nna ta ta ta tba ta ta tt tba ta ta ta tb ijia b0iia 将该方程组改写为 的显函数的形式：先合理地假设一组节点温度的初始值 ，代入上方程组，求得一组节点温度值 ；再将 代入上方程组，又求得一组新的温度值 ；以此类推。 当max 或 结束迭代。