化工传递工程;第七章-热传导-精品课件.ppt

1、化工传递过程基础 1化工传递过程基础 2 第七章第七章 热热 传传 导导 热传导热传导( (导热导热) )是介质内无宏观运动时的传热现象，是介质内无宏观运动时的传热现象，导热在固体、液体和气体中均可发生，但严格而言，导热在固体、液体和气体中均可发生，但严格而言，只有在固体中才是单纯的导热，而流体即使处于静止只有在固体中才是单纯的导热，而流体即使处于静止状态，其中也会由于温度梯度所造成的密度差而产生状态，其中也会由于温度梯度所造成的密度差而产生自然对流，因此在流体中对流与导热同时发生。鉴于自然对流，因此在流体中对流与导热同时发生。鉴于此，本章将针对固体中的热传导问题进行讨论，重点此，本章将针对固

2、体中的热传导问题进行讨论，重点研究某些情况下热传导方程的求解方法，并结合实际研究某些情况下热传导方程的求解方法，并结合实际情况，探讨一些导热理论在工程实际中的应用。情况，探讨一些导热理论在工程实际中的应用。 描述导热的基本微分方程已在第六章中导出如描述导热的基本微分方程已在第六章中导出如式式(6-17a)(6-17a)所示，即所示，即kqtt21化工传递过程基础 3式式(6-17a)(6-17a)在不同坐标系的一般形式为在不同坐标系的一般形式为直角坐标直角坐标柱坐标柱坐标球坐标球坐标kqztytxtt2222221kqzttrrtrrrt22222111kqtrtrrtrrrt222222si

3、n1sinsin111(7-1)(7-1)(7-3)(7-3)(7-2)(7-2)化工传递过程基础 4 求解热传导的规律问题，即解出上述微分方程，获求解热传导的规律问题，即解出上述微分方程，获得温度得温度t t与时间与时间 及位置及位置(z(z，y y，z)z)的函数关系，即不的函数关系，即不同时刻温度在空间的分布同时刻温度在空间的分布( (温度场温度场) )，所得的解为，所得的解为t=t=f(xf(x，y y，z)z)，它不但要满足式，它不但要满足式(7-1)(7-1)或式或式(7-2)(7-2)、式、式(7-3)(7-3)，而且要满足每一问题的初始条件与边界条件。，而且要满足每一问题的初始

4、条件与边界条件。 上述热传导方程的求解方法是相当复杂的，除了几上述热传导方程的求解方法是相当复杂的，除了几种简单的典型问题可以采用数学分析方法求解外，绝种简单的典型问题可以采用数学分析方法求解外，绝大部分问题常常需要采用特殊的方法，例如数值计算大部分问题常常需要采用特殊的方法，例如数值计算等方法进行求解。等方法进行求解。本章将主要针对以直角坐标系和柱本章将主要针对以直角坐标系和柱坐标系表达的某些简单的工程实际导热问题的求解方坐标系表达的某些简单的工程实际导热问题的求解方法进行研究。法进行研究。化工传递过程基础 5 第一节第一节 稳态热传导稳态热传导 一、无内热源的一维稳态热传导一、无内热源的一

5、维稳态热传导 对于无内热源的一维稳态热传导，由于温度与时间对于无内热源的一维稳态热传导，由于温度与时间无关，无关， ，且无内热源，且无内热源， 。又设沿。又设沿x x或或r r方向进行一维导热，则热传导方程方向进行一维导热，则热传导方程(7-1)(7-1)、(7-2)(7-2)、(7-3)(7-3)可简化为一可简化为一维维的拉普拉斯方程，即的拉普拉斯方程，即直角坐标直角坐标柱坐标柱坐标0t0q 022dxtd0drdtrdrd)47( )57( kqztytxtt2222221化工传递过程基础 6球坐标球坐标 工程上一维工程上一维( (沿沿x x或或r r方向方向) )稳态热传导的例子很多，稳

6、态热传导的例子很多，如方形燃烧炉的炉壁、蒸汽管的管壁、列管式换热器如方形燃烧炉的炉壁、蒸汽管的管壁、列管式换热器的管壁以及球形压力容器的器壁等。的管壁以及球形压力容器的器壁等。21,)2(, 0) 1 (ttbxttx02drdtrdrd ( (一一) )单层平壁一维稳态热传导单层平壁一维稳态热传导 单层平壁单层平壁( (如方形燃烧炉的炉壁如方形燃烧炉的炉壁) )沿一个方向的导沿一个方向的导热问题是最简单的热传导问题，当导热系数热问题是最简单的热传导问题，当导热系数k k 为常数为常数时，式时，式(7-4)(7-4)即为描述该导热过程的微分方程，即即为描述该导热过程的微分方程，即设边界条件为设

7、边界条件为022dxtd化工传递过程基础 7将式将式(7-7)(7-7)积分两次，可得积分两次，可得)87( 21CxCt21,CCxbtttt211式中，式中， 为积分常数，代人边界条件为积分常数，代人边界条件(1)(1)，可求，可求出出 ；代入边界条件；代入边界条件(2),(2),可求出可求出 。将将 代入式代入式(7-7)(7-7)，即可得到此情况下的温度，即可得到此情况下的温度分布方程为分布方程为1C2C12tC bttC)(121 由式由式(7-8)(7-8)可知可知, ,平壁稳态热传导过程的温度分平壁稳态热传导过程的温度分布为一条直线。该式也可由傅立叶定律导出布为一条直线。该式也可

8、由傅立叶定律导出 求出温度分布之后，便可进一步求出沿求出温度分布之后，便可进一步求出沿x x方向通过方向通过平壁的导热通量。根据傅立叶定律平壁的导热通量。根据傅立叶定律, ,通过某通过某x x处的导热处的导热通量通量q/A q/A 可表示为可表示为21,)2(, 0) 1 (ttbxttx化工传递过程基础 8dxdtkAq将式将式(7-8)(7-8)对对x x求导后代入上式，得求导后代入上式，得)(21ttbkAqxbtttt211化工传递过程基础 9(二二) )单层单层筒筒壁的稳态热传导壁的稳态热传导 化工生产中，经常遇到筒壁的热传导问题，求解化工生产中，经常遇到筒壁的热传导问题，求解筒壁的

9、径向热传导问题时，应用柱坐标系比较方便，筒壁的径向热传导问题时，应用柱坐标系比较方便，若筒壁的长度很长，若筒壁的长度很长，L Lr r，则沿轴向的导热可略去不，则沿轴向的导热可略去不计，于是可认为温度仅沿径向变化，在此情况下，描计，于是可认为温度仅沿径向变化，在此情况下，描述无内热源的一维稳态热传导方程为式述无内热源的一维稳态热传导方程为式(7-5)(7-5)，即，即 设边界条件为设边界条件为 0drdtrdrd2211,)2(,) 1 (ttrrttrrkqzttrrtrrrt22222111化工传递过程基础 10将式将式(7-5)(7-5)积分两次，可得积分两次，可得式中，式中，C C1

10、1、C C2 2为积分常数，经向该式代人边界条件为积分常数，经向该式代人边界条件(1)(1)和和(2)(2)后，可得后，可得)ln(12121rrttC21lnCrCt1121212ln)ln(rrrtttC 将将C C1 1、C C2 2代人式代人式(7-11)(7-11) ，即可得到沿筒壁径向一维稳态即可得到沿筒壁径向一维稳态导热时的温度分布方程为导热时的温度分布方程为112211ln)ln(rrrrtttt2211,)2(,) 1 (ttrrttrr(7-12)(7-12)化工传递过程基础 11上式表明，通过筒壁进行径向一维稳态热传导时，温上式表明，通过筒壁进行径向一维稳态热传导时，温度

11、分布是度分布是r r的对数函数。的对数函数。 通过半径为通过半径为r r的筒壁处的传热速率或热通量，可由的筒壁处的传热速率或热通量，可由柱坐标系的傅立叶定律导出即柱坐标系的傅立叶定律导出即drdtkAqrdrdtkAqr 式中，式中，q q和和q qArAr分别为半径分别为半径r r处的导热速率和热通量；处的导热速率和热通量；ArAr为该处的导热面积，为该处的导热面积，ArAr=2=2丌丌rLrL，其中，其中L L为筒壁的长度；为筒壁的长度； 为该处的温度梯度。为该处的温度梯度。将式将式(7-12)(7-12)对对r r求导并代入式求导并代入式(7-13)(7-13)和式和式(7-13a)(7

12、-13a)可得可得)ln(21221rrttkLq)ln(1221rrttrkAqrdrdt(7-13)7-13)112211ln)ln(rrrrtttt化工传递过程基础 12 式式(7-14)(7-14)即为单层筒壁的导热速率方程。以上诸即为单层筒壁的导热速率方程。以上诸式表明，式表明，尽管壁温、筒壁的传热面积和热通最均随半尽管壁温、筒壁的传热面积和热通最均随半径径r r而变，但传热速率在稳态时依然是常量，即而变，但传热速率在稳态时依然是常量，即常量LrAqLrAqLrAqq3322112)(2)(2)(常量或332211)()()(rAqrAqrAq式式(7-14)(7-14)亦可写成与平

13、壁导热速率方程相类似的形式，即亦可写成与平壁导热速率方程相类似的形式，即1221rrttkAqm将式将式(7-17)(7-17)与式与式(7-14)(7-14)对比，可知对比，可知LrLrrrrAmm2)ln(21212(7-15)(7-15)ln(21221rrttkLq(7-17(7-17)化工传递过程基础 13mmAr式中式中 简壁的对数平均半径；简壁的对数平均半径； 筒壁的对数平均面积。筒壁的对数平均面积。 应予指出，当应予指出，当 2 2时，上述各式中的对数平时，上述各式中的对数平均值可用算术平均值代替。均值可用算术平均值代替。 通常，筒壁的导热速率采用单位筒长来表示，则通常，筒壁的

14、导热速率采用单位筒长来表示，则由式由式(7-14)(7-14)可得可得12221212ln22ln22AAAALrLrLrLrAm或12rr)ln(21221rrttkLq化工传递过程基础 14 以上均假定导热系数以上均假定导热系数k k为与温度无关的常数。当为与温度无关的常数。当k k为为温度温度t t的线性函数时，上述各式中的导热系数的线性函数时，上述各式中的导热系数k k可采可采用用 、 算术平均温度下的值算术平均温度下的值 来代替。来代替。1tmk2t 二、有内热源的一维稳态热传导二、有内热源的一维稳态热传导 有内热源的导热设备，以柱体最为典型，例如核反有内热源的导热设备，以柱体最为典

15、型，例如核反应堆的铀棒、管式固定床反应器和电热棒等。应堆的铀棒、管式固定床反应器和电热棒等。若柱体若柱体很长，且温度分布沿轴向对称，在此情况下的稳态导很长，且温度分布沿轴向对称，在此情况下的稳态导热问题，可视为沿径向的一维稳态热传导，热问题，可视为沿径向的一维稳态热传导，此时，柱此时，柱坐标下的能量方程式坐标下的能量方程式(7-2)(7-2)可化为可化为01kqrtrrrkqzttrrtrrrt22222111化工传递过程基础 15式式(7(7- -19)19)系用来描述具有内热源、沿径向作一维稳态系用来描述具有内热源、沿径向作一维稳态热传导时的微分方程。若内热源均匀，则热传导时的微分方程。若

16、内热源均匀，则 为常数为常数。 结合具体的边界条件可求出柱体内的温度分布，为结合具体的边界条件可求出柱体内的温度分布，为此，将式此，将式(7-19)(7-19)进行第一次积分，进行第一次积分，得得再积分一次，又得再积分一次，又得 式中，式中， 为积分常数，可根据两个边界条件确定，为积分常数，可根据两个边界条件确定，具体方法参见例具体方法参见例7-17-1和例和例7-27-2。q )217( 212ln4CrCrkqt21,CCrCrkqdrdt12)207( 01kqrtrrr化工传递过程基础 16例例7-1 7-1 有一半径为有一半径为R R，长度为，长度为L L的实心圆柱体，其发热的实心圆

17、柱体，其发热 速率为速率为 ，圆柱体的表面温度为，圆柱体的表面温度为 ，L LR R，温度仅，温度仅为径向距离的函数，设热传导是稳态的，圆柱体的导为径向距离的函数，设热传导是稳态的，圆柱体的导热系数热系数k k为常数，试求圆柱体内的温度分布及最高温度为常数，试求圆柱体内的温度分布及最高温度处的温度值。处的温度值。 解：柱体内一维径向稳态热传导时的温度分布方程为解：柱体内一维径向稳态热传导时的温度分布方程为Wtq 212ln4CrCrkqtdrdtRLkLRqRrttRrW2,)2(,) 1 (2依题意，设边界条件为依题意，设边界条件为化工传递过程基础 17边界条件边界条件(2)(2)表示稳态热

18、传导时，圆柱体内的发热速表示稳态热传导时，圆柱体内的发热速率必等于表面热损失速率。率必等于表面热损失速率。 由由边界边界条件条件(2)(2)可得可得将上式代入式将上式代入式(7-20)(7-20)中，并取中，并取r=Rr=R，得，得故故 将将C Cl l=0=0及边界条件及边界条件(1)(1)代人式代人式(7-21)(7-21)中，得中，得故故 01CRCkRqkRq122rCrkqdrdt12kRqdrdtRr2|224CkRqtWkRqtCW422drdtRLkLRqRr2,)2(2212ln4CrCrkqt化工传递过程基础 18最后解出温度分布为最后解出温度分布为由于圆柱体向外导热，显然

19、最高温度在圆柱体中心处，由于圆柱体向外导热，显然最高温度在圆柱体中心处，即即201RrttttWW)(422rRkqttWkRqttttWr4|200max上式亦可写成无因次形式，即上式亦可写成无因次形式，即(上两式相比)化工传递过程基础 19【例例7-27-2】 有一外径为有一外径为4cm4cm、内径为、内径为1.5 cm1.5 cm、载有电、载有电流密度流密度I I为为5000 A5000 A 的内冷钢制导体。导体单位时的内冷钢制导体。导体单位时间发出的热量等于流体同时带走的热量，导体内壁面间发出的热量等于流体同时带走的热量，导体内壁面的温度为的温度为7070。假定外壁面完全绝热。试确定导

20、体内。假定外壁面完全绝热。试确定导体内部的温度分布并求算导体内部最高温度处的温度值。部的温度分布并求算导体内部最高温度处的温度值。 已知钢的热传导系数已知钢的热传导系数k=380 wk=380 w( (m mK K) ) 电阻率电阻率2cm1Cm8102解：由式解：由式(7-21)(7-21)出发，求出导体内部的温度分布。为出发，求出导体内部的温度分布。为此首先求出此首先求出 、 、 各值。各值。q 2C72482105)105000)(102(Iq3/mW根据题意，可知本题的两个边界条件为根据题意，可知本题的两个边界条件为化工传递过程基础 20212ln4CrCrkqt3 .261C,2)2

21、(70,75. 0) 1 (211cmrtcmr00drdtdrdtk即将边界条件将边界条件(2)(2)代人式代人式(720)(720)中，得中，得002. 002. 0380210517Cdrdt由此得由此得再将边界条再将边界条件件(1)(1)代入式代入式(7-21)(7-21)中，得中，得2227)1075. 0ln(3 .2610075. 0380210570C解之得解之得6 .2002CrCrkqdrdt12化工传递过程基础 21将将C C1 1、C C2 2代入式代入式(7-21)(7-21)中，即可求出导体内部的温度中，即可求出导体内部的温度分布方程为分布方程为或或最高温度发生在外

22、壁面处，该处最高温度发生在外壁面处，该处r r2 2=2cm =0.02 m=2cm =0.02 m，故，故6 .200ln3 .26328952rrt6 .200ln3 .26380210527rrt6 .846 .20002. 0ln3 .26)02. 0(328952maxt化工传递过程基础 22 三、二维稳态热传导三、二维稳态热传导 上面讨论的稳态热传导问题，其温度均可以用一个上面讨论的稳态热传导问题，其温度均可以用一个空间坐标的函数来表示，但空间坐标的函数来表示，但工程中还常遇到二维或三维工程中还常遇到二维或三维稳态热传导问题。对于这类问题，仅当边界条件比较简稳态热传导问题。对于这类

23、问题，仅当边界条件比较简单时，才有可能应用分析解法，单时，才有可能应用分析解法，但求解过程相当麻烦，但求解过程相当麻烦，结果也很复杂，不便于工程应用。而对于边界条件比较结果也很复杂，不便于工程应用。而对于边界条件比较复杂的热传导问题，其分析求解就更加困难，有的甚至复杂的热传导问题，其分析求解就更加困难，有的甚至根本不能得到分析解，此时，根本不能得到分析解，此时，解决问题最有效的方法是解决问题最有效的方法是数值计算法，数值计算法，这种方法有许多优越性，特别是计算机的这种方法有许多优越性，特别是计算机的迅速发展，使得人们能够对以前认为不能求解的许多问迅速发展，使得人们能够对以前认为不能求解的许多问

24、题得到数值解。题得到数值解。下面以无内热源的二维稳态热传导为下面以无内热源的二维稳态热传导为例，说明数值计算法的应用。例，说明数值计算法的应用。 化工传递过程基础 23 ( (一一) )物体内部的结点温度方程物体内部的结点温度方程 无内热源的二维稳态热传导，在直角坐标系中以二无内热源的二维稳态热传导，在直角坐标系中以二维的拉普拉斯方程描述，即维的拉普拉斯方程描述，即根据上式求出的温度分布根据上式求出的温度分布t=t=f(zf(z，y)y)为一连续曲面，为一连续曲面，数值计算法的基础是将上述连续变化的偏微分方程用数值计算法的基础是将上述连续变化的偏微分方程用差分方程近似表达，从而求出温度分布。差

25、分方程近似表达，从而求出温度分布。 如图如图7-17-1所示，将物体分割成若干个由所示，将物体分割成若干个由x x、y y组组成的小方格，分割线的交点称为结点，成的小方格，分割线的交点称为结点，x x及及y y的长的长度视对计算精度的要求选取，度视对计算精度的要求选取，x x或或y y越小，所得结越小，所得结果就越接近于真实温度分布，当然相应的计算量也就果就越接近于真实温度分布，当然相应的计算量也就越大。越大。 02222ytxtkqztytxtt2222221化工传递过程基础 24化工传递过程基础 25温度梯度可以写为温度梯度可以写为xttxtjijiji, 1,21|xttxtjijiji

26、, 1,21|yttytjijiji,1,21,|yttytjijiji1,21,|化工传递过程基础 262, 1, 1,21,21,222|xtttxxtxtxtjijijijijiji2,1,1,21,21,222|ytttyytytytjijijijijiji由此，可将式由此，可将式(7-22)(7-22)近似地写成差分形式，即近似地写成差分形式，即令令x=x=y y，上式化为，上式化为0222,1,1,2, 1, 1ytttxtttjijijijijiji04,1, 1, 11jijiijijitttttj02222ytxt(7-23)(7-23)化工传递过程基础 27 式式(7-23

27、)(7-23)称为物体内部的结点温度分布方程，它表称为物体内部的结点温度分布方程，它表达任一结点达任一结点(i(i，j)j)的温度的温度t ti,ji,j与邻近四个结点温度之间与邻近四个结点温度之间的关系，即在无内热源的二维稳态温度场中，其内部的关系，即在无内热源的二维稳态温度场中，其内部某结点的温某结点的温度度可用邻近四个结点温度的算术平均值表可用邻近四个结点温度的算术平均值表示。示。显然，若将所有内部结点的温度分别与其相邻显然，若将所有内部结点的温度分别与其相邻的的四个结点的温度按式四个结点的温度按式(7-23)(7-23)的形式联系起来，便可建的形式联系起来，便可建立物体内部的结点温度方

28、程组。立物体内部的结点温度方程组。 ( (二二) )物体边界上的结点温度方程物体边界上的结点温度方程 处于物体表面的结点，由于外界的影响，其温度处于物体表面的结点，由于外界的影响，其温度就不能应用式就不能应用式(7-23)(7-23)来表达，而要根据具体情况来来表达，而要根据具体情况来建立。建立。 化工传递过程基础 28 简单的边界情况如图简单的边界情况如图7-2(a)7-2(a)(d)(d)所示，图所示，图7-2(a)7-2(a)为绝热边界，其余三种为对流边界，下面分别建立此为绝热边界，其余三种为对流边界，下面分别建立此四种边界情况下的结点温度方程，为简便计，推导时，四种边界情况下的结点温度

29、方程，为简便计，推导时，均取垂直纸面的距离为单位长度。均取垂直纸面的距离为单位长度。 1 1绝热边界绝热边界化工传递过程基础 29化工传递过程基础 30如图如图7-2(a)7-2(a)所示，对虚线包围的微元体作热量衡算所示，对虚线包围的微元体作热量衡算, ,得得022,1,1, 1xyttkxyttkyxttkjijijijijiji042, 11,1,jiijijittttj)(22,1,1, 1,bjijijijijijijittyhyttxkyttxkxttyk令令x=x=y y，则上式化为，则上式化为2 2对流边界对流边界如图如图7-2(b)7-2(b)所示，设周围流体主体温度为所示，

30、设周围流体主体温度为t tb b且维持且维持不变，微元体表面与流体之间的对流传热系数为不变，微元体表面与流体之间的对流传热系数为h h，亦维持不变，对虚线包围的微元体作热量衡算，可亦维持不变，对虚线包围的微元体作热量衡算，可得得化工传递过程基础 31bjijijiitkxhtkxhtttj,1,1, 1)2()2(21bjijiitkxhtkxhttj2) 1(2,1, 1bjijijijiitkxhtkxhttttj2)3(222,1, 11, 1令令x=x=y y，则上式化为，则上式化为 同理可求得图同理可求得图7-2(c)7-2(c)中对流边界上的外角结点中对流边界上的外角结点( (i,

31、ji,j) )的结点温度方程为的结点温度方程为 及图及图 7-2(d)7-2(d)中对流边界上的内角结点中对流边界上的内角结点( (i,ji,j) )的结点的结点温度方程为温度方程为化工传递过程基础 32 ( (三三) )二维稳态温度场的结点温度方程组二维稳态温度场的结点温度方程组 式式(7-23)(7-23)、(7-24)(7-24)表达了无内热源二维稳态温度场表达了无内热源二维稳态温度场中结点温度之间的关系，各式均为线性代数方程。中结点温度之间的关系，各式均为线性代数方程。求求解温度场时，可根据物体内部及边界情况，并考虑精解温度场时，可根据物体内部及边界情况，并考虑精度要求，将物体分割成若

32、干个等边的小方格，将分割度要求，将物体分割成若干个等边的小方格，将分割线的交点统一编号，线的交点统一编号，i=li=l，2 2，3 3，n n，然后根据每，然后根据每个结点所在的位置，分别写出相应的结点温度方程，个结点所在的位置，分别写出相应的结点温度方程，从而得到整个温度场的结点温度方程组，即从而得到整个温度场的结点温度方程组，即 nnnnnnnnnnbtatatabtatatabtatata22112222212111212111化工传递过程基础 33 式中，式中， 和和 ( (i,ji,j=1=1，2 2，n)n)均为常数；均为常数； (i=1 (i=1，2 2，n)n)为未知温度，式为

33、未知温度，式(7-25)(7-25)为线性方为线性方程组，共有程组，共有n n个方程个方程, ,未知温度亦为未知温度亦为n n个，求解此方程个，求解此方程组即可得出组即可得出 的数值，于是整个温度场即的数值，于是整个温度场即可解出。可解出。ib2mitnttt,21jia, 【例例7-37-3】如附图所示，某一边长为如附图所示，某一边长为1m1m的正方形物的正方形物体，体，左侧面恒温为左侧面恒温为100100，顶面恒温为，顶面恒温为500500，其余，其余两侧面暴露在对流环境中，环境温度为两侧面暴露在对流环境中，环境温度为100100。已知。已知物体导热系数为物体导热系数为l0(m)l0(m)

34、，物体与环境的对流传热，物体与环境的对流传热系数为系数为10 W/( )10 W/( )，试建立，试建立1 19 9各结点的温度方各结点的温度方程组并求出各点的温度值。程组并求出各点的温度值。 求解上述结点温度方程组可采用求逆矩阵法求解上述结点温度方程组可采用求逆矩阵法,迭代法迭代法 和高斯消去法等。和高斯消去法等。 化工传递过程基础 34(1)(1)建立结点温度方程组建立结点温度方程组 由于内部和边界上的结点温由于内部和边界上的结点温度方程不同，今以内部结点度方程不同，今以内部结点1 1及边界上的结点及边界上的结点3 3、9 9为代为代表建立各结点温度方程。表建立各结点温度方程。 对于结点对

35、于结点l l，应用式，应用式(7-24a)(7-24a)，得，得31kxh04500100142ttt)/(10),/(10,100313CmWhCmWkC，tx：b已知解6004421ttt结点结点3 3为一般对流边界上的点，应用式为一般对流边界上的点，应用式(7-24b)(7-24b)得得btkxhtkxhtt362)2()5002(2104,1, 1, 11jijiijijitttttjbjijijiitkxhtkxhtttj,1,1, 1)2()2(21化工传递过程基础 35化工传递过程基础 36代人数据，得代人数据，得56767. 42632ttt结点结点9 9为对流边界外角上的点，

36、应用式为对流边界外角上的点，应用式(7-24c)(7-24c)得得btkxhtkxhtt2) 1(2986代人数据，得代人数据，得7 .6667. 2986ttt 其余各结点的温度方程可用相应的方程建立，最后其余各结点的温度方程可用相应的方程建立，最后得得l l9 9各点的结点温度方程组为各点的结点温度方程组为bjijiitkxhtkxhttj2) 1(2,1, 1化工传递过程基础 376004421ttt56767. 42632ttt7 .6667. 2986ttt0486542ttttt50045321tttt10047541tttt7 .6667. 429653tttt16767. 42

37、874ttt7 .6667. 429875tttt化工传递过程基础 38(2)(2)各点的温度数值的计算结果各点的温度数值的计算结果 采用求逆矩阵法求解上述方程组，可得采用求逆矩阵法求解上述方程组，可得() () 1904t3073t1567t2146t2791t2275t1828t1739t3272t化工传递过程基础 39 第二节第二节 不稳态热传导不稳态热传导 物体内任一点的温度均随时间而变的导热称为不稳物体内任一点的温度均随时间而变的导热称为不稳态导热。态导热。在工程实际中，经常遇到不稳态导热问题，在工程实际中，经常遇到不稳态导热问题，例如燃烧炉的点火升温过程和熄火降温过程、金属的例如燃

38、烧炉的点火升温过程和熄火降温过程、金属的熔化、淬火等热加工处理均为不稳态导热。熔化、淬火等热加工处理均为不稳态导热。此外有些此外有些稳态导热问题，在其初始阶段也常存在不稳态导热过稳态导热问题，在其初始阶段也常存在不稳态导热过程，如燃烧炉的点火阶段即是如此。程，如燃烧炉的点火阶段即是如此。 由于不稳态导热过程中的温度既与时间有关又与位置由于不稳态导热过程中的温度既与时间有关又与位置有关，故其求解要较稳态导热问题复杂得多。有关，故其求解要较稳态导热问题复杂得多。通常求通常求解不稳态导热问题时，需应用热传导方程解不稳态导热问题时，需应用热传导方程(7-1)(7-1)、(7-2)(7-2)或或(7-3

39、)(7-3)，并需满足具体的初始条件及边界条件。通过，并需满足具体的初始条件及边界条件。通过求解满足这些定解条件的偏微分方程，求得温度分布求解满足这些定解条件的偏微分方程，求得温度分布随时随时间间的变化关系，从而求得特定时刻的传热速率。的变化关系，从而求得特定时刻的传热速率。化工传递过程基础 40 初始条件是指在导热过程开始的瞬时，物体内部的初始条件是指在导热过程开始的瞬时，物体内部的温度分布情况。温度分布情况。边界条件视具体情况一般可分为三类：边界条件视具体情况一般可分为三类：第一类边界条件是给出任何时刻物体端面的温度分布；第一类边界条件是给出任何时刻物体端面的温度分布；第二类边界条件是给出

40、所有时刻物体端面处的导热通第二类边界条件是给出所有时刻物体端面处的导热通量；第三类边界条件是物体端面与周围流体介质进行量；第三类边界条件是物体端面与周围流体介质进行热交换，端而处的导热速率等于端面与流体之间对流热交换，端而处的导热速率等于端面与流体之间对流传热速率。传热速率。 不稳态导热过程中的传热速率取决于介质内部热阻不稳态导热过程中的传热速率取决于介质内部热阻和表面热阻，根据它们的相对大小和表面热阻，根据它们的相对大小不稳态导热过不稳态导热过程程可可以分为三种情况，即忽略内部热阻、忽略表面热阻和以分为三种情况，即忽略内部热阻、忽略表面热阻和两种热阻都不能忽略。两种热阻都不能忽略。化工传递过

41、程基础 41 一、忽略内部热阻的不稳态导热与集总热容法一、忽略内部热阻的不稳态导热与集总热容法 内部热阻可忽略的不稳态导热问题是一种最简单的内部热阻可忽略的不稳态导热问题是一种最简单的不稳态导热问题，不稳态导热问题，若固体的导热系数很大或内热阻很若固体的导热系数很大或内热阻很小，而环境流体与该固体表面之间的对流传热热阻又小，而环境流体与该固体表面之间的对流传热热阻又比较大时，便可忽略内热阻比较大时，便可忽略内热阻，即认为在任一时刻固体，即认为在任一时刻固体内部各处的温度均匀一致。这种假设物体内部热阻与内部各处的温度均匀一致。这种假设物体内部热阻与外部热阻相比，可忽略不计的一种分析方法称为集总外

42、部热阻相比，可忽略不计的一种分析方法称为集总热容法。热容法。 例如有一个热的金属小球，被浸泡在冷的油类或其例如有一个热的金属小球，被浸泡在冷的油类或其他流体中他流体中( (参见图参见图7-3)7-3)，显然小球的温度分布除与其显然小球的温度分布除与其材质的导热系数有关外，还和小球表面与周围流体的材质的导热系数有关外，还和小球表面与周围流体的对流传热系数有关。假定对流传热系数有关。假定小球导热良好，其导热热阻小球导热良好，其导热热阻比表面对流热阻小得多比表面对流热阻小得多，则主要的温度梯度将产生于，则主要的温度梯度将产生于小球表面的流体层内，而小球本身的温度在任一瞬时小球表面的流体层内，而小球本

43、身的温度在任一瞬时化工传递过程基础 42化工传递过程基础 43均可认为是均匀一致的。均可认为是均匀一致的。 设金属球的密度为设金属球的密度为 ，比热容为，比热容为c c，体积为，体积为V V，表面，表面积为积为A A，初始温度均匀，为，初始温度均匀，为 ，环境流体的主体温度恒，环境流体的主体温度恒定，为定，为 。流体与金属球表面的对流传热系数为。流体与金属球表面的对流传热系数为h h，且，且不随时间而变不随时间而变。又设在又设在 时间内，金属球的温度变化时间内，金属球的温度变化为为dtdt。根据热量衡算，整个金属球的放热速率应等于其。根据热量衡算，整个金属球的放热速率应等于其表面与流体间的对流

44、传热速率，即表面与流体间的对流传热速率，即bt)(btthAddtVc0td0, 0tt 初始条件为初始条件为 化工传递过程基础 44 式中，式中，t t为任一瞬时金属球表面的温度，由于金属球为任一瞬时金属球表面的温度，由于金属球的内热阻可以忽略不计，故金属球各点的温度均为的内热阻可以忽略不计，故金属球各点的温度均为t t0 0式中的式中的“负负”号表示球内的温度随时间而降低。号表示球内的温度随时间而降低。由于由于物体的温度仅随时间改变而与位置无关，因此不存在物体的温度仅随时间改变而与位置无关，因此不存在边界条件。边界条件。 令令 ，则式，则式(7-26)(7-26)可化为可化为dVchAdb

45、tt 0, 0初始条件为初始条件为 积分式积分式(7-27)(7-27)得得dVchAd00)(btthAddtVc化工传递过程基础 45222)/()(VAkAVhcVkAkAhVVchAVchA0lnVchAbbetttt00或或式式(7-28)(7-28)即为忽略物体内热阻情况下，物体温度与时间即为忽略物体内热阻情况下，物体温度与时间的定量关系式。的定量关系式。 式式(7-28)(7-28)中右侧指数中的量还可以写成如下形式中右侧指数中的量还可以写成如下形式式式(7-29)(7-29)右侧的两个数群都是无因次的，现在对该二右侧的两个数群都是无因次的，现在对该二数群的物理意义作进一步分析。

46、数群的物理意义作进一步分析。 (7-28)(7-28)化工传递过程基础 46第一个数群第一个数群 称为毕渥数，记为称为毕渥数，记为BiBi，即，即 lkAVh)/(khlkAVhBi)/(由于毕渥数中的由于毕渥数中的V VA A具有长度因次具有长度因次( (以以 表示表示) )，故毕，故毕渥数的物理意义为渥数的物理意义为 对流传热热阻导热热阻长度导热系数对流传热系数长度)()(Bi即毕渥数表示了物体内部的导热热阻与表面对流热阻即毕渥数表示了物体内部的导热热阻与表面对流热阻之比。之比。BiBi值大时，表示传热过程中物体内部的导热热值大时，表示传热过程中物体内部的导热热阻起控制作用，阻起控制作用，

47、物体内部存在较大的温度梯度，此时物体内部存在较大的温度梯度，此时系统的传热不能采用集总热容法处理；反之系统的传热不能采用集总热容法处理；反之BiBi值小时，值小时，则表示物体内部的热阻很小，表面对流传热的热则表示物体内部的热阻很小，表面对流传热的热(7-30)(7-30)化工传递过程基础 47阻起控制作用，物体内部的温度梯度很小，在同一瞬阻起控制作用，物体内部的温度梯度很小，在同一瞬时各处温度较为均匀。时各处温度较为均匀。研究表明，当研究表明，当Bi0.1Bi0.1时，系统时，系统的传热可采用集总热容法处理，的传热可采用集总热容法处理，此时用式此时用式(7-28)(7-28)计算计算物体温度与

48、时间的关系，其结果与实际比较，误差不物体温度与时间的关系，其结果与实际比较，误差不超过超过5 5。因此，求解不稳态传热问题时，首先要计。因此，求解不稳态传热问题时，首先要计算算BiBi的值，视其是否小于的值，视其是否小于0.10.1，以便确定该传热问题，以便确定该传热问题能否采用集总热容法处理。能否采用集总热容法处理。2)/(AV第二个数群第二个数群 称为称为“傅立叶数傅立叶数”，记为，记为F F0 0，即，即220)/(lAVF傅立叶数的物理意义表示时间之比，即无因次时间。傅立叶数的物理意义表示时间之比，即无因次时间。 将式将式(7-30)(7-30)、(7-31)(7-31)代入式代入式(

49、7-28)(7-28)中，得中，得(7-31)(7-31)化工传递过程基础 48【例例7-47-4】 有一半径有一半径r r0 0为为25 mm25 mm的钢球，初始温度均匀，的钢球，初始温度均匀，为为700 K700 K，突然将此球放入某流体介质中，介质的温，突然将此球放入某流体介质中，介质的温度恒定，为度恒定，为400 K400 K。假定钢球表面与流体之间的对流。假定钢球表面与流体之间的对流传热系数为传热系数为h=11.36 wh=11.36 w ，且不随温度而变。，且不随温度而变。钢球的物性值为：导热系数钢球的物性值为：导热系数k=43.3 wk=43.3 w( (mKmK) )密密度度

50、 ，比热容，比热容c=0.46 kJc=0.46 kJ( (kgKkgK) )。试计算试计算1 1小时后钢球的温度。小时后钢球的温度。)(2Km 3/7849mkg00FBbbietttt解：由于解：由于h h值较小，值较小，k k值较大，估计可以采用集总热值较大，估计可以采用集总热容法，为此首先计算容法，为此首先计算BiBi数。数。化工传递过程基础 49化工传递过程基础 50Kt477mrrrAV33020301033. 83102534)34(1 . 000219. 03 .43)1033. 8)(36.11()/(3khlkAVhBi故故故可用式故可用式(7-28)(7-28)计算计算l