清华大学热工基础ppt课件工程热力学加传热学10第九章 导热、稳态导热、非稳态、数值解法.ppt

《清华大学热工基础ppt课件工程热力学加传热学10第九章 导热、稳态导热、非稳态、数值解法.ppt》由会员分享，可在线阅读，更多相关《清华大学热工基础ppt课件工程热力学加传热学10第九章 导热、稳态导热、非稳态、数值解法.ppt（96页珍藏版）》请在三一办公上搜索。

1、1,第九章导 热,主要内容,本章首先阐述导热的基本概念、基本定律、导热问题的数学描述方法，为进一步求解导热问题奠定必要的理论基础，然后讨论几种简单的稳态导热、非稳态导热的分析解法，最后简要介绍导热问题的数值解法。,2,研究方法,从连续介质的假设出发，从宏观的角度来讨论导热热流量与物体温度分布及其他影响因素之间的关系。,连续介质,一般情况下，绝大多数固体、液体及气体都可以看作连续介质。但是当分子的平均自由行程与物体的宏观尺寸相比不能忽略时，如压力降低到一定程度的稀薄气体，就不能认为是连续介质。,3,9-1 导热理论基础,主要内容,（1）与导热有关的基本概念；,（2）导热基本定律；,（3）导热现象

2、的数学描述方法。,为进一步求解导热问题奠定必要的理论基础。,1.导热的基本概念,（1）温度场(temperature field),在 时刻，物体内所有各点的温度分布称为该物体在该时刻的温度场。,4,一般温度场是空间坐标和时间的函数，在直角坐标系中，温度场可表示为,非稳态温度场,：温度随时间变化的温度场，其中的导热称为非稳态导热。,稳态温度场,：温度不随时间变化的温度场，其中的导热称为稳态导热。,一维温度场,二维温度场,三维温度场,5,（2）等温面与等温线,在同一时刻，温度场中温度相同的点连成的线或面称为等温线或等温面。,等温面上任何一条线都是等温线。如果用一个平面和一组等温面相交,就会得到一

3、组等温线。温度场可以用一组等温面或等温线表示。,等温面与等温线的特征：,同一时刻，物体中温度不同的等温面或等温线不能相交；在连续介质的假设条件下，等温面（或等温线）或者在物体中构成封闭的曲面（或曲线），或者终止于物体的边界，不可能在物体中中断。,6,（3）温度梯度(temperature gradient),在温度场中，温度沿x方向的变化率(即偏导数),很明显,等温面法线方向的温度变化率最大，温度变化最剧烈。,温度梯度：等温面法线方向的温度变化率矢量：,n等温面法线方向的单位矢量，指向温度增加的方向。,温度梯度是矢量，指向温度增加的方向。,7,在直角坐标系中，温度梯度可表示为,分别为x、y、z

4、 方向的偏导数;i、j、k 分别为x、y、z 方向的单位矢量。,（4）热流密度(heat flux),热流密度的大小和方向可以用热流密度矢量q 表示,热流密度矢量的方向指向温度降低的方向。,8,在直角坐标系中，热流密度矢量可表示为,qx、qy、qz分别表示q在三个坐标方向的分量的大小。,2.导热的基本定律,傅里叶（Fourier）于1822年提出了著名的导热基本定律，即傅里叶定律，指出了导热热流密度矢量与温度梯度之间的关系。,对于各向同性物体,傅里叶定律表达式为,傅里叶定律表明,导热热流密度的大小与温度梯度的绝对值成正比，其方向与温度梯度的方向相反。,9,标量形式的傅里叶定律表达式为,对于各向

5、同性材料,各方向上的热导率相等,由傅里叶定律可知,要计算导热热流量,需要知道材料的热导率,还必须知道温度场。所以,求解温度场是导热分析的主要任务。,10,傅里叶定律的适用条件：,（1）傅里叶定律只适用于各向同性物体。对于各向异性物体，热流密度矢量的方向不仅与温度梯度有关,还与热导率的方向性有关,因此热流密度矢量与温度梯度不一定在同一条直线上。,（2）傅里叶定律适用于工程技术中的一般稳态和非稳态导热问题，对于极低温（接近于0K）的导热问题和极短时间产生极大热流密度的瞬态导热过程,如大功率、短脉冲(脉冲宽度可达10-1210-15s)激光瞬态加热等,傅里叶定律不再适用。,11,3.热导率（导热系数

6、）,热导率表明物质导热能力的大小。根据傅里叶定律表达式,绝大多数材料的热导率值都可以通过实验测得。,12,物质的热导率在数值上具有下述特点:,(1)对于同一种物质,固态的热导率值最大,气态的热导率值最小；,(2)一般金属的热导率大于非金属的热导率；,(3)导电性能好的金属,其导热性能也好；,(4)纯金属的热导率大于它的合金；,(5)对于各向异性物体,热导率的数值与方向有关；,(6)对于同一种物质,晶体的热导率要大于非定形态物体的热导率。,热导率数值的影响因素较多,主要取决于物质的种类、物质结构与物理状态,此外温度、密度、湿度等因素对热导率也有较大的影响。其中温度对热导率的影响尤为重要。,13,

7、温度对热导率的影响：,一般地说,所有物质的热导率都是温度的函数,不同物质的热导率随温度的变化规律不同。,纯金属的热导率随温度的升高而减小。,一般合金和非金属的热导率随温度的升高而增大。,大多数液体（水和甘油除外）的热导率随温度的升高而减小。,纯金属的热导率随温度的升高而减小。,14,在工业和日常生活中常见的温度范围内,绝大多数材料的热导率可近似地认为随温度线性变化,并表示为,0为按上式计算的0下的热导率值，并非热导率的真实值,如图所示。b为由实验确定的常数，其数值与物质的种类有关。,15,多孔材料的热导率,绝大多数建筑材料和保温材料（或称绝热材料）都具有多孔或纤维结构(如砖、混凝土、石棉、炉渣

8、等),不是均匀介质,统称多孔材料。,多孔材料的热导率是指它的表观热导率,或称作折算热导率。,用于保温或隔热的材料。国家标准规定，温度低于350时热导率小于0.12 W/(mK)的材料称为保温材料。,保温材料（或称绝热材料）：,16,多孔材料的热导率随温度的升高而增大。,多孔材料的热导率与密度和湿度有关。一般情况下密度和湿度愈大，热导率愈大。,典型材料热导率的数值范围,纯金属 50415 W/(mK)合金 12120 W/(mK)非金属固体 140 W/(mK)液体(非金属)0.170.7 W/(mK)绝热材料 0.030.12 W/(mK)气体 0.0070.17 W/(mK),17,4.导热

9、问题的数学描述（数学模型）,（1）导热微分方程式的导出,导热微分方程式+单值性条件,建立数学模型的目的：求解温度场,依据：能量守恒和傅里叶定律。,假设：,1）物体由各向同性的连续介质组成；,2）有内热源，强度为，表示单位时间、单位体积内的生成热，单位为W/m3。,1）根据物体的形状选择坐标系,选取物体中的微元体作为研究对象；,导热数学模型的组成：,步骤：,2）根据能量守恒,建立微元体的热平衡方程式；,3）根据傅里叶定律及已知条件,对热平衡方程式进行归纳、整理，最后得出导热微分方程式。,18,导热过程中微元体的热平衡：,单位时间内，净导入微元体的热流量d与微元体内热源的生成热dV之和等于微元体热

10、力学能的增加dU,即,d+dV=dU,d=dx+dy+dz,dx=dx dx+dx,=qx dydz qx+dx dydz,19,同理可得从y和z方向净导入微元体的热流量分别为,于是,在单位时间内净导入微元体的热流量为,单位时间内微元体内热源的生成热：,单位时间内微元热力学能的增加：,根据微元体的热平衡表达式 d+dV=dU 可得,20,导热微分方程式建立了导热过程中物体的温度随时间和空间变化的函数关系。,当热导率为常数时,导热微分方程式可简化为,式中2是拉普拉斯算子,在直角坐标系中,或写成,称为热扩散率,也称导温系数,单位为m2/s。,其大小反映物体被瞬态加热或冷却时温度变化的快慢。,木材a

11、=1.510-7 紫铜a=5.3310-5,21,导热微分方程式的简化,(1)物体无内热源：,(2)稳态导热：,(3)稳态导热、无内热源：,2t=0，即,22,圆柱坐标系下的导热微分方程式,如果为常数：,23,球坐标系下的导热微分方程式,为常数时,24,（2）单值性条件 导热微分方程式推导过程中没有涉及导热过程的具体特点,适用于无穷多个导热过程,也就是说有无穷多个解。为完整地描写某个具体的导热过程，必须说明导热过程的具体特点,即给出导热微分方程的单值性条件（或称定解条件），使导热微分方程式具有唯一解。导热微分方程式与单值性条件一起构成具体导热过程完整的数学描述。单值性条件一般包括：几何条件、物

12、理条件、时间条件、边界条件。,25,1)几何条件,说明参与导热物体的几何形状及尺寸。几何条件决定温度场的空间分布特点和分析时所采用的坐标系。,2)物理条件,说明导热物体的物理性质,例如物体有无内热源以及内热源的分布规律，给出热物性参数(、c、a等)的数值及其特点等。,3)时间条件,说明导热过程时间上的特点,是稳态导热还是非稳态导热。对于非稳态导热,应该给出过程开始时物体内部的温度分布规律（称为初始条件）：,26,4)边界条件,说明导热物体边界上的热状态以及与周围环境之间的相互作用。例如,边界上的温度、热流密度分布以及边界与周围环境之间的热量交换情况等。,常见的边界条件分为以下三类:,(a)第一

13、类边界条件,给出边界上的温度分布及其随时间的变化规律：,(b)第二类边界条件,给出边界上的热流密度分布及其随时间的变化规律：,27,(c)第三类边界条件,给出了与物体表面进行对流换热的流体的温度tf及表面传热系数h。,用电热片加热物体表面可实现第二类边界条件。,如果物体的某一表面是绝热的,即qw=0,则,物体内部的等温面或等温线与该绝热表面垂直相交。,根据边界面的热平衡，由傅里叶定律和牛顿冷却公式可得,第三类边界条件建立了物体内部温度在边界处的变化率与边界处对流换热之间的关系，也称为对流换热边界条件。,28,上式描述的第三类边界条件是线性的,所以也称为线性边界条件，反映了导热问题的大部分实际情

14、况。,如果导热物体的边界处除了对流换热还存在与周围环境之间的辐射换热,则边界面的热平衡表达式为,qr 为物体边界面与周围环境之间的净辐射换热热流密度，它与物体边界和周围环境的温度和辐射特性有关,是温度的复杂函数。这种对流换热与辐射换热叠加的复合换热边界条件是非线性的边界条件。,本书只限于讨论具有线性边界条件的导热问题。,29,对一个具体导热过程完整的数学描述（即导热数学模型）应该包括,建立合理的数学模型,是求解导热问题的第一步,也是最重要的一步。,目前应用最广泛的求解导热问题的方法有:(1)分析解法、(2)数值解法、(3)实验方法。这也是求解所有传热学问题的三种基本方法。,(1)导热微分方程式

15、;,(2)单值性条件。,对数学模型进行求解,就可以得到物体的温度场,进而根据傅里叶定律就可以确定相应的热流分布。,本章主要介绍导热问题的分析解法和数值解法。,30,9-2 稳态导热,稳态导热是指温度场不随时间变化的导热过程.,下面分别讨论日常生活和工程上常见的平壁、圆筒壁、球壁及肋壁的一维稳态导热问题。,1.平壁的稳态导热,当平壁的两表面分别维持均匀恒定的温度时，平壁的导热为一维稳态导热。,假设：,表面面积为A、厚度为、为常数、无内热源，两侧表面分别维持均匀恒定的温度tw1、tw2，且tw1 tw2。,（1）单层平壁的稳态导热,选取坐标轴x与壁面垂直,如图所示。,31,数学模型：,x=0,t=

16、tw1,x=,t=tw2,求解结果：,可见，当为常数时,平壁内温度分布曲线为直线，其斜率为,由傅里叶定律可得热流密度,通过整个平壁的热流量为,上式与绪论中给出的公式完全相同。,32,变热导率问题：,当平壁材料的热导率是温度的函数时,平壁一维稳态导热的数学模型为,当温度变化范围不大时,可以近似地认为材料的热导率随温度作线性变化,即,可见,当平壁材料的热导率随温度作线性变化时,平壁内的温度分布为二次曲线。,求解数学模型可得平壁内的温度分布为,33,根据傅里叶定律表达式,（1）当tw1 tw2 时，热流方向与x轴同向，q为正值,而热导率数值永远为正,所以由上式可见,温度变化率为负值。,（2）如果b0

17、,沿x方向随温度的降低而减小,温度曲线斜率的绝对值增大,曲线向上弯曲（上凸）；,（3）如果b0,温度曲线向下弯曲。,平壁内的温度分布二次曲线形状的讨论：,b0,b0,34,根据傅里叶定律，可由温度分布求得平壁的热流密度,为平壁的算术平均温度；,为平壁的算术平均温度下的热导率。,上式说明，当热导率随温度线性变化时,通过平壁的热流量可用热导率为常数时的计算公式来计算,但公式中的热导率为平壁算术平均温度下的热导率m。,式中,对照,35,（2）多层平壁的稳态导热,多层平壁由多层不同材料组成，当两表面分别维持均匀恒定的温度时，其导热也是一维稳态导热。,以三层平壁为例，假设,（1）各层厚度分别为1、2、3

18、，各层材料的热导率分别为1、2、3,且分别为常数；,（2）各层之间接触紧密,相互接触的表面具有相同的温度；,（3）平壁两侧外表面分别保持均匀恒定的温度tw1、tw4。,显然，通过此三层平壁的导热为稳态导热,各层的热流量相同。,36,三层平壁稳态导热的总导热热阻为各层导热热阻之和，由单层平壁稳态导热的计算公式可得,三层平壁稳态导热可以由三个相互串联的热阻网络表示。,由此类推,对于n层平壁的稳态导热,利用热阻的概念,可以很容易求得通过多层平壁稳态导热的热流量,进而求出各层间接触面的温度。,37,2.圆筒壁的稳态导热,（1）单层圆筒壁的稳态导热,主要讨论圆筒壁稳态导热过程中的壁内温度分布及导热热流量

19、。,假设：,内、外半径分别为 r1、r2，长度为l，为常数，无内热源，内外壁温度tw1、tw2均匀恒定。,按上述条件，壁内温度只沿径向变化，如果采用圆柱坐标,则圆筒壁内的导热为一维稳态导热，,数学模型,r=r1:t=tw1,r=r2:t=tw1,38,对导热微分方程式进行两次积分,可得通解为,圆筒壁内的温度分布为对数曲线。代入边界条件,可得,温度沿r 方向的变化率为,其绝对值沿r 方向逐渐减小。,根据傅里叶定律,沿圆筒壁r 方向的热流密度为,热流密度是r的函数。,39,但对于稳态导热,通过整个圆筒壁的热流量是不变的,R为整个圆筒壁的导热热阻,单位是K/W。,单位长度圆筒壁的热流量为,Rl为单位

20、长度圆筒壁的导热热阻,单位是mK/W。,实际上，由于l为常数,根据傅里叶定律，,将该式分离变量积分，同样可求得上面的公式。,40,（2）多层圆筒壁的稳态导热,运用热阻的概念，很容易分析多层圆筒壁的稳态导热问题。,以三层圆筒壁为例，无内热源，各层的热导率1、2、3均为常数，内、外壁面维持均匀恒定的温度tw1、tw2。这显然也是一维稳态导热问题。通过各层圆筒壁的热流量相等，总导热热阻等于各层导热热阻之和:,41,对于n层不同材料组成的多层圆筒壁的稳态导热,单位长度的热流量为,42,4.肋片的稳态导热,加装肋片的目的：,增加换热面积A，强化对流换热。,对流换热的牛顿冷却公式：=A h(twtf),4

21、3,几种常见的肋片,44,（1）等截面直肋的稳态导热分析,以矩形肋为例进行分析:,矩形肋的高度为H、厚度为、宽度为l，与高度方向垂直的横截面积为A,横截面的周长为U。,假设：,1）肋片材料热导率为常数；,2）肋片根部与肋基接触良好，温度一致；,3）肋片厚度方向的导热热阻/与表面的对流换热热阻1/h相比很小，可以忽略,肋片温度只沿高度方向发生变化,肋片导热可以近似地认为是一维的；,4）肋片表面各处对流换热的表面传热系数h都相同；,5）忽略肋片端面的散热，即认为肋端面是绝热的。,45,（2）将肋片导热看作是具有负的内热源的一维稳态导热。,数学模型：,内热源强度的确定：,对于图中所示的微元段，,肋片

22、导热微分方程的两种导出方法：,（1）由肋片微元段的热平衡导出；,46,代入导热微分方程式，得,令,称为过余温度。,数学模型变为,双曲余弦函数,47,肋片的过余温度从肋根开始沿高度方向按双曲余弦函数的规律变化，,肋片的过余温度沿高度方向逐渐降低，mH较小时，温度降低缓慢；mH较大时，温度降低较快。,一般取0.7 mH 2,48,肋端，x=H，肋端的过余温度,肋端过余温度随mH增加而降低。,在稳态情况下,肋片散热量应该等于从肋根导入的热量,随着mH增大，散热量增加，开始增加迅速，后来越来越缓慢，逐渐趋于一渐近值。(增加肋高的经济性）,49,（2）肋片效率,肋片效率定义：,肋片的实际散热量与假设整个

23、肋片都具有肋基温度时的理想散热量0之比,式中tm、m分别为肋面的平均温度和平均过余温度，t0、0分别为肋基温度与肋基过余温度。,由于m 0，所以肋片效率f 小于1。,因为假设肋表面各处h都相等，所以等截面直肋的平均过余温度可按下式计算：,50,矩形和三角形肋片效率随mH的变化规律如下图所示。,肋片效率的影响因素：,（1）肋片材料的热导率、,（2）肋片高度H，,（3）肋片厚度，,（4）肋片与周围流体间对流换热的表面传热系数h，,可见，mH愈大，肋片效率愈低。,51,（1）上述分析结果同样适用于其它形状的等截面直肋，如圆柱、圆管形肋的一维稳态导热问题；,（2）如果必须考虑肋端面的散热，可以将肋端面

24、面积折算到侧面上去，相当于肋加高为H+H，其中,对于矩形肋,几点说明：,（4）对于肋片厚度方向的导热热阻/与表面的对流换热热阻1/h相比不可忽略的情况，肋片的导热不能认为是一维的，上述公式不再适用；,（5）上述推导没有考虑辐射换热的影响，对一些温差较大的场合，必须加以考虑。,（3）上述分析结果既适用于肋片被加热的情况，也适用于肋片被冷却的情况；,52,肋片导热对热电偶测温精度的影响,例题94,热电偶测量的是测温套管端部的温度tH。,只考虑套管的导热。在稳态情况下，套管端部温度不等于空气的温度，测温误差就是套管端部的过余温度。,忽略测温套管横截面上的温度变化，并认为端部绝热，则套管导热可以看成是

25、等截面直肋的一维稳态导热问题。,如何减小测温误差？,53,6.接触热阻,接触热阻的定义：,由于固体表面之间不能完全接触而对两个固体间的导热过程产生的热阻，用Rc表示。,由于存在接触热阻，使两个接触表面之间出现温差,接触热阻的主要影响因素：,(1)相互接触的物体表面的表面粗糙度；,(2)相互接触的物体表面的硬度；,(3)相互接触的物体表面之间的压力等。,减小接触热阻的措施：抛光、加压、添加薄膜等。,54,9-3 非稳态导热,非稳态导热,：温度场随时间变化的导热过程。,非稳态导热的类型：,（1）周期性非稳态导热：,（2）非周期性非稳态导热：,在周期性变化边界条件下发生的导热过程，如内燃机汽缸壁的导

26、热、一年四季大地土壤的导热等。,在瞬间变化的边界条件下发生的导热过程，例如热处理工件的加热或冷却等。,讨论一维非周期性非稳态导热的分析解法及求解特殊非稳态导热问题的集总参数法。,了解和掌握非稳态导热过程中温度场的变化规律及换热量的计算方法。,本节主要内容：,主要目的：,55,1.一维非稳态导热问题的分析解,第三类边界条件下大平壁、长圆柱及球体的加热或冷却是工程上常见的一维非稳态导热问题。,（1）无限大平壁冷却或加热问题的分析解简介,假设：,厚度为、热导率、热扩散率a为常数，无内热源，初始温度与两侧的流体相同并为t0。两侧流体温度突然降低为t，并保持不变，平壁表面与流体间对流换热表面传热系数h为

27、常数。,考虑温度场的对称性，选取坐标系如图，仅需讨论半个平壁的导热问题。,这是一维的非稳态导热问题。,56,1)数学模型：,（对称性）,引进无量纲过余温度、,无量纲坐标，,Fo是无量纲特征数，称为傅里叶数,称为毕渥数,57,傅里叶数的物理意义：,Fo为两个时间之比，是非稳态导热过程的无量纲时间。,毕渥数的物理意义：,Bi为物体内部的导热热阻与边界处的对流换热热阻之比。,由无量纲数学模型可知，是Fo、Bi、X三个无量纲参数的函数,确定此函数关系是求解该非稳态导热问题的主要任务。,2）求解结果：,58,解的函数形式为无穷级数，式中1,2,，n 是下面超越方程的根,根有无穷多个，是Bi的函数。无论B

28、i取任何值，1,2,，n都是正的递增数列，的解是一个快速收敛的无穷级数。,由解的函数形式可以看出，确实是Fo、Bi、X三个无量纲特征数的函数,59,（2）分析解的讨论,1）傅里叶数 Fo 对温度分布的影响,分析解的计算结果表明，当Fo 0.2时，可近似取级数的第一项，对工程计算已足够精确，即,因为，所以将上式左、右两边取对数，可得,，m为一与时间、地点无关的常数，只取决于第三类边界条件、平壁的物性与几何尺寸。,式中,式右边的第二项只与Bi、x/有关，与时间 无关。,60,上式可改写为,该式说明，当Fo0.2时，即 时，平壁内所有各点过余温度的对数都随时间线性变化，并且变化曲线的斜率都相等，这一

29、温度变化阶段称为非稳态导热的正规状况阶段。,上式两边求导，可得,m的物理意义是过余温度对时间的相对变化率，单位是s1，称为冷却率（或加热率）。,上式说明，当Fo 0.2，进入正规状况阶段后，所有各点的冷却率都相同，且不随时间而变化，其大小取决于物体的物性、几何形状与尺寸及表面传热系数。,61,对于平壁中心，,上面两式之比,可见，当Fo 0.2，非稳态导热进入正规状况阶段以后，虽然与m都随时间变化，但它们的比值与时间无关，只取决于毕渥数Bi与几何位置x/。,认识正规状况阶段的温度变化规律具有重要的实际意义，因为工程技术中的非稳态导热过程绝大部分时间都处于正规状况阶段。,62,2）毕渥数Bi对温度

30、分布的影响,平壁非稳态导热第三类边界条件表达式,上式的几何意义：在整个非稳态导热过程中平壁内过余温度分布曲线在边界处的切线都通点,即，该点称为第三类边界条件的定向点。,63,毕渥数 Bi 对温度分布的影响分析,(a)Bi0:,平壁的导热热阻趋于零，平壁内部各点温度在任一时刻都趋于一致，只随时间而变化，变化的快慢取决于平壁表面的对流换热强度。定向点在无穷远处。,工程上只要Bi0.1，就可以近似地按这种情况处理，用集总参数法进行计算。,64,对流换热热阻趋于零，非稳态导热一开始平壁表面温度就立即变为流体温度，相当于给定了壁面温度，即给定了第一类边界条件，平壁内部的温度变化完全取决于平壁的导热热阻。

31、定向点位于平壁表面上。当Bi100时可近似按此处理。,(b)Bi:,(c)0Bi100，按一般情况处理。,65,3）平壁与周围流体之间交换的热量,在0时间内，微元薄层dx单位面积放出的热量等于其热力学能的变化，,在0时间内，单位面积平壁放出的热量,将Fo 0.2时无量纲过余温度的近似解代入上式，得,66,（3）诺谟图,1）,67,2）,68,3）,69,几点说明：,(1)上述分析是针对平壁被冷却的情况进行的，但分析结果对平壁被加热的情况同样适用；,(2)由于平壁温度场是对称的，所以分析时只取半个平壁作为研究对象，这相当于一侧（中心面）绝热、另一侧具有第三类边界条件的情况，因此分析结果也适用于同

32、样条件的平壁；,(3)线算图只适用于Fo 0.2的情况;,70,(4)对于圆柱体和球体在第三类边界条件下的一维非稳态导热问题，分别在柱坐标系和球坐标系下进行分析，也可以求得温度分布的分析解，解的形式也是快速收敛的无穷级数，并且是Bi、Fo和r/R的函数,(5)当Fo 0.2时，圆柱和球体的一维非稳态导热过程也都进入正规状况阶段，分析解可以近似地取无穷级数的第一项，近似结果也被绘成了线算图。（P 221）,71,2.特殊多维非稳态导热问题的简易求解方法,(1)无限长方柱,(2)短圆柱,(3)垂直六面体,72,3.集总参数法（Bi0.1）,当Bi0.1时，物体内部的导热热阻远小于其表面的对流换热热

33、阻，可以忽略，物体内部各点的温度在任一时刻都近似于均匀，物体的温度只是时间的函数。对于这种情况，只需求解物体温度随时间的变化规律以及物体放出或吸收的热量。,集总参数法:忽略物体内部导热热阻的简化分析方法。,73,假设：一个任意形状的物体，体积为V，表面面积为A，密度、比热容c及热导率为常数，无内热源，初始温度为t0。突然将该物体放入温度t恒定的流体中，物体表面和流体之间对流换热的表面传热系数h为常数。假设该问题满足Bi0.1的条件。,74,根据能量守恒，单位时间物体热力学能的变化量应该等于物体表面与流体之间的对流换热量，,下角标V表示以 l=V/A为特征长度,75,几点说明：,（1）集总参数法

34、中的毕渥数BiV与傅里叶数FoV以l=V/A为特征长度，不同于分析解中的Bi与Fo，,（2）对于形状如平板、柱体或球的物体，只要满足Bi0.1，就可以使用集总参数法计算，偏差小于5。,76,令,，c称为时间常数。,当c时，,即物体的过余温度达到初始过余温度的36.8%。说明时间常数反映物体对周围环境温度变化响应的快慢，时间常数越小，物体的温度变化越快。,热电偶的时间常数越小越好，因为越小，热电偶越能迅速地反映被测流体的温度变化。,影响时间常数大小的主要因素是,由 可见，,物体的热容量cV和物体表面的对流换热条件hA。,77,在0 时间内物体和周围环境之间交换的热量,令,，表示物体温度从t0 变

35、化到周围流体温度t所放出或吸收的总热量，则上式改写为,上述分析结果既适用于物体被加热的情况，也适用于物体被冷却的情况。,78,9-4 导热问题的数值解法基础,数值解法的基本思想:,用导热问题所涉及的空间和时间区域内有限个离散点(称为节点)的温度近似值来代替物体内实际连续的温度分布，将连续温度分布函数的求解问题转化为各节点温度值的求解问题，将导热微分方程的求解问题转化为节点温度代数方程的求解问题。,79,数值解法的基本内容与步骤:,（1）对实际导热问题的几何、物理性质进行分析，做必要的、合理的简化，建立符合实际的物理模型。,（2）根据物理模型建立完整的数学模型，即给出导热微分方程和单值性条件。,

36、第(1)、(2)步是导热问题所有求解方法的基础。,（3）求解域离散化：用与坐标轴平行的网络线将所涉及的空间和时间区域划分成有限个子区域，将网络线的交点作为节点,每个节点就代表以它为中心的子区域（控制容积），节点温度就代表子区域的温度。,80,（4）建立节点温度代数方程组；,（5）求解节点温度代数方程组，得到所有节点的温度值；,（6）对计算结果进行分析，若不符合实际情况，则修正上述步骤，重复进行计算，直到结果满意为止。,目前求解导热问题常用的数值解法主要有：有限差分法、有限元法等。其中有限差分法比较成熟，应用广泛。下面主要介绍有限差分法的基本原理。,81,1.有限差分法的基本原理,以常物性、无内

37、热源的二维稳态导热为例:,(1)求解域的离散化,1)子区域的划分,选择网格宽度x、y（步长），划分子区域。步长大小根据问题的需要而定。,2)节点的选择,选择网线和网线及网线与物体边界的交点作为节点，标定节点位置，如(i,j)、(i+1,j)等。,82,(2)节点温度差分方程的建立,两种方法：泰勒级数展开法与控制容积热平衡法。只介绍控制容积热平衡法。,控制容积热平衡法：根据节点所代表的控制容积在导热过程中的能量守恒来建立节点温度差分方程。,1)内部节点温度差分方程,内部节点(i,j)所代表的控制容积在导热过程中的热平衡,对于垂直于纸面方向单位宽度,选择x=y,83,上式可整理为,可见，物体内每一

38、个节点温度都等于相邻4个节点温度的算术平均值。,2)边界节点温度差分方程,对于具有第三类边界条件的边界节点(i,j)所代表的控制容积，根据其热平衡,选择步长x=y，将上式简化,84,令,称为网格毕渥数。,上式可整理为,第三类边界条件下的外拐角边界节点：,第三类边界条件下的内拐角边界节点：,85,绝热边界节点：,运用有限差分方法可以建立导热物体所有内部节点和边界节点温度的差分方程。求解这些差分方程构成一个线性代数方程组就可以得节点温度的数值。,86,2.节点温度差分方程组的求解方法,线性代数方程组的求解方法有消元法、矩阵求逆法、迭代法等，这里仅简单介绍在导热的数值计算中常用的迭代法中的两种：,(

39、1)简单迭代法,(2)高斯-塞德尔迭代法,87,(1)简单迭代法,其中aij、bi为常数，且aij0。改写为显函数形式：,假设,88,(2)高斯-塞德尔迭代法,高斯-塞德尔迭代法是在简单迭代法的基础上加以改进的迭代运算方法。它与简单迭代法的主要区别是在迭代运算过程中总使用最新算出的数据。,高斯-塞德尔迭代法比简单迭代法收敛速度快。,89,3.非稳态导热问题的数值解法,非稳态导热数值解法的特点：,（1）非稳态导热微分方程多了非稳态项，因此单值性条件中增加了初始条件；,（2）除了对空间域进行离散外，还需要对时间域进行离散；,（3）利用热平衡法导出节点温度方程时需要考虑控制容积的热力学能随时间的变化

40、；,（4）由于时间和空间同时离散，在有些情况下空间步长和时间步长不能任意选择，否则会带来节点温度方程求解的稳定性问题。,90,以第三类边界条件下常物性、无内热源的无限大平壁的一维非稳态导热问题为例。,(1)求解域的离散,(2)节点温度差分方程的建立,运用热平衡法可以建立非稳态导热物体内部节点和边界节点温度差分方程。,空间步长为x，时间步长为,x、大小的选择需要保证节点温度方程求解的稳定性。,表示空间节点i在k时刻（简称k时刻）的节点温度。,91,1)内部节点温度差分方程,内部节点i所代表的控制容积在k时刻的热平衡：,如果节点i的温度对时间的变化率采用向前差分，热平衡方程式可写成,令,内部节点温

41、度方程的显式差分格式,92,两点结论：,(a)任意一个内部节点i在(k+1)时刻的温度都可以由该节点及其相邻节点在k 时刻的温度由上式直接求出，不必联立求解方程组，这是显式差分格式的优点。这样就可以从初始温度出发依次求出各时刻的节点温度。,(b)必须满足显式差分格式的稳定性条件，即,物理意义？,稳定性条件说明，一旦空间步长x或时间步长的数值确定之后，另一个步长的数值的就不能任意选择，必须满足稳定性条件。,93,1)边界节点温度差分方程,边界节点0所代表的控制容积在k 时刻的热平衡：,如果节点0的温度对时间的变化率采用向前差分，热平衡方程式可写成,引进网格傅里叶数,和网格毕渥数,上式写成显函数的

42、形式,边界节点温度方程的显式差分格式,94,同内部节点温度方程的显式差分格式的道理一样，上式必须满足显式差分格式的稳定性条件，即,因为，所以只要满足上式，自然满足内部节点温度方程显式差分格式的稳定性条件。因此上式是第三类边界条件下一维非稳态导热所有节点温度方程显式差分格式的稳定性条件。,对比,95,第九章小结,（1）傅里叶定律表达式及其适用条件;,重点掌握以下内容：,（3）导热问题的数学描述（数学模型）;,（2）物体导热系数的数值范围及特点;,（4）平壁、圆筒壁、球壁及肋壁的一维稳态导热 问题的分析求解;,96,（5）非稳态导热的特点、无限大平壁冷却或加热问题的分析求解方法、分析解的特点及其影响因素(Bi,Fo);,（6）求解特殊非稳态导热问题的集总参数法；,（7）了解导热问题数值解法的基本思想与步骤及有限差分法的基本原理。,