第3章-非稳态导热分析解法
传热学 第3章-非稳态导热分析解法
第三章 非稳态导热分析解法1、 重点内容:① 非稳态导热的基本概念及特点;② 集总参数法的基本原理及应用;③一维及二维非稳态导热问题。
2、掌握内容:① 确定瞬时温度场的方法;② 确定在一时间间隔内物体所传导热量的计算方法。
3、了解内容:无限大物体非稳态导热的基本特点。
许多工程问题需要确定:物体内部温度场随时间的变化,或确定其内部温度达某一极限值所需的时间。
如:机器启动、变动工况时,急剧的温度变化会使部件因热应力而破坏。
因此,应确定其内部的瞬时温度场。
钢制工件的热处理是一个典型的非稳态导热过程,掌握工件中温度变化的速率是控制工件热处理质量的重要因素;金属在加热炉内加热时,要确定它在炉内停留的时间,以保证达到规定的中心温度。
§3—1 非稳态导热的基本概念一、非稳态导热1、定义:物体的温度随时间而变化的导热过程称非稳态导热。
2、分类:根据物体内温度随时间而变化的特征不同分:1)物体的温度随时间的推移逐渐趋于恒定值,即:const t =↑τ2)物体的温度随时间而作周期性变化1)物体的温度随时间而趋于恒定值如图3-1所示,设一平壁,初值温度t 0,令其左侧的表面温度突然升高到1t 并保持不变,而右侧仍与温度为0t 的空气接触,试分析物体的温度场的变化过程。
首先,物体与高温表面靠近部分的温度很快上升,而其余部分仍保持原来的t 0 。
如图中曲线HBD ,随时间的推移,由于物体导热温度变化波及范围扩大,到某一时间后,右侧表面温度也逐渐升高,如图中曲线HCD 、HE 、HF 。
最后,当时间达到一定值后,温度分布保持恒定,如图中曲线HG (若λ=const ,则HG 是直线)。
由此可见,上述非稳态导热过程中,存在着右侧面参与换热与不参与换热的两个不同阶段。
(1)第一阶段(右侧面不参与换热)温度分布显现出部分为非稳态导热规律控制区和部分为初始温度区的混合分布,即:在此阶段物体温度分布受t 分布的影响较大,此阶段称非正规状况阶段。
第三章非稳态导热
第三章⾮稳态导热第三章⾮稳态导热的分析计算 3-1 ⾮稳态导热过程分析⼀、⾮稳态导热过程及其特点导热系统(物体)内温度场随时间变化的导热过程为⾮稳态导热过程。
在过程的进⾏中系统内各处的温度是随时间变化的,热流量也是变化的。
这反映了传热过程中系统内的能量随时间的改变。
我们研究⾮稳态导热过程的意义在于,⼯程上和⾃然界存在着⼤量的⾮稳态导热过程,如房屋墙壁内的温度变化、炉墙在加热(冷却)过程中的温度变化、物体在炉内的加热或在环境中冷却等。
归纳起来,⾮稳态导热过程可分为两⼤类型,其⼀是周期性的⾮稳态导热过程,其⼆是⾮周期性的⾮稳态导热过程,通常指物体(或系统)的加热或冷却过程。
这⾥主要介绍⾮周期性的⾮稳态导热过程。
下⾯以⼀维⾮稳态导热为例来分析其过程的主要特征。
今有⼀⽆限⼤平板,突然放⼊加热炉中加热,平板受炉内烟⽓环境的加热作⽤,其温度就会从平板表⾯向平板中⼼随时间逐渐升⾼,其内能也逐渐增加,同时伴随着热流向平板中⼼的传递。
图3-1显⽰了⼤平板加热过程的温度变化的情况。
从图中可见,当0=τ时平板处于均匀的温度0t t =下,随着时间τ的增加平板温度开始变化,并向板中⼼发展,⽽后中⼼温度也逐步升⾼。
当∞→τ时平板温度将与环境温度拉平,⾮稳态导热过程结束。
图中温度分布曲线是⽤相同的?τ来描绘的。
总之,在⾮稳态导热过程中物体内的温度和热流都是在不断的变化,⽽且都是⼀个不断地从⾮稳态到稳态的导热过程,也是⼀个能量从不平衡到平衡的过程。
⼆、加热或冷却过程的两个重要阶段从图3-1中也可以看出,在平板加热过程的初期,初始温度分布0t t =仍然在影响物体整个的温度分布。
只有物体中⼼的温度开始变化之后(如图中τ>τ2之后),初始温度分布0t t =的影响才会消失,其后的温度分布就是⼀条光滑连续的曲线。
据此,我们可以把⾮稳态导热过程分为两个不同的阶段,即:初始状况阶段――环境的热影响不断向物体内部扩展的过程,也就是物体(或系统)仍然有部分区域受初始温度分布控制的阶段;正规状况阶段――环境对物体的热影响已经扩展到整个物体内部,且仍然继续作⽤于物体的过程,也就是物体(或系统)的温度分布不再受初始温度分布影响的阶段。
第三章 非稳态导热-108
x 0,
t 0 x
(对称性)
x ,
t x
ht
t
以上导热微分方程式及单值性条件组成了该非稳态导热 问题的数学模型。
• 引进过余温度: t t
于是导热微分方程式和单值性条件变为:
2
a x2
初始条件: 边界条件:
0,
x 0,
x ,
0 t0 t
0 x
h x
• 再引进无量纲温度: 0
需要确定该物体在冷却过程中温度随时间的变化规律 以及放出的热量。
A
,c,V h, t
• 假设该问题满足Bi 0.1的条件,
根据能量守恒,单位时间物体热力学能的变化量应该
等于物体表面与流体之间的对流换热量,即
cV
dt
d
hAt t
引进过余温度 t t ,上式可改写为
cV d hA d
初始条件为 0, 0 t0 t
e 12 Fo
f Bi,Fo
上式也同样被绘制成线算图。
4 分析解应用范围的推广和讨论
(1)分析解应用范围的三点推广。
①对物体冷却也可适用。 ②对一侧绝热、另一侧为第三类边界条件的平板也可适用。 ③当对流换热系数趋于无穷大时,固体的表面温度就趋近于流体的温度 ,因而Bi→∞时就是物体表面温度突然变化后保持不变的第一类边界条件 的解。
正规状况阶段,分析解可近似取无穷级数的第一项,近似结果也 被绘制成了线算图。
3. 集 总 参 数 法
• 当 Bi 0.1 时,物体内部的导热热阻远小于其表面的
对流换热热阻,可以忽略,物体内部各点的温度在任 一时刻都趋于均匀,物体的温度只是时间的函数,与 坐标无关。 • 对于这种情况下的非稳态导热问题,只须求出温度随 时间的变化规律以及在温度变化过程中物体放出或吸 收的热量。
第3章-非稳态导热分析解法
V c
hA
hA
V c
=1
0
e 1 36.8%
Vc
hA
上式表明:当传热时间等于
时,物体的过余
Vc
hA
温度已经达到了初始过余温度的36.8%。称 为时间常数,也称弛豫时间,用 c 表示。
V c
hA
时间常数反映了系统处于一定的环境中所表现出来
的传热动态特征,与其几何形状、密度及比热有关,
数时,即τ=τc,则
θ /θ
0
0
=e
1
0.386
1
τ=4τc,时
0
=e
4.6
0.386 0 1 τ /τ
0.01
s
工程上认为= 4τc时导热体已达到热平衡状态
3
Bi Biຫໍສະໝຸດ Fohl l 物理意义
= 物体内部导热热阻 物体表面对流换热热阻
无量纲 热阻
无量纲 时间
非正规状况阶段(右侧面不参与换热 ):温度
分布显现出部分为非稳态导热规律控制区和部分
为初始温度区的混合分布,即:在此阶段物体温
度分布受 t 分布的影响较大。 环境的热影响不断向物体内部扩展的过程, 即物体(或系统)有部分区域受到初始温度分布 控制的阶段。必须用无穷级数描述。
正规状况阶段(右侧面参与换热 ): 当右侧面参与换热以后,物体中的温度分布
求物体温度随时间变化的依 变关系
A
ΔΕ
ρ , c, V, t0
φc
h, t
建立数学模型-利用两种方法
根据导热微分方程的一般形式进行简化; 利用能量守恒 热平衡关系为:内热能随时间的变化率ΔΕ=通
传热学第3章-非稳态导热分析解法
传热学第3章-⾮稳态导热分析解法第三章⾮稳态导热分析解法1、重点内容:①⾮稳态导热的基本概念及特点;②集总参数法的基本原理及应⽤;③⼀维及⼆维⾮稳态导热问题。
2、掌握内容:①确定瞬时温度场的⽅法;②确定在⼀时间间隔内物体所传导热量的计算⽅法。
3、了解内容:⽆限⼤物体⾮稳态导热的基本特点。
许多⼯程问题需要确定:物体内部温度场随时间的变化,或确定其内部温度达某⼀极限值所需的时间。
如:机器启动、变动⼯况时,急剧的温度变化会使部件因热应⼒⽽破坏。
因此,应确定其内部的瞬时温度场。
钢制⼯件的热处理是⼀个典型的⾮稳态导热过程,掌握⼯件中温度变化的速率是控制⼯件热处理质量的重要因素;⾦属在加热炉内加热时,要确定它在炉内停留的时间,以保证达到规定的中⼼温度。
§3—1 ⾮稳态导热的基本概念⼀、⾮稳态导热1、定义:物体的温度随时间⽽变化的导热过程称⾮稳态导热。
2、分类:根据物体内温度随时间⽽变化的特征不同分:1)物体的温度随时间的推移逐渐趋于恒定值,即:const t =↑τ2)物体的温度随时间⽽作周期性变化1)物体的温度随时间⽽趋于恒定值如图3-1所⽰,设⼀平壁,初值温度t 0,令其左侧的表⾯温度突然升⾼到1t 并保持不变,⽽右侧仍与温度为0t 的空⽓接触,试分析物体的温度场的变化过程。
⾸先,物体与⾼温表⾯靠近部分的温度很快上升,⽽其余部分仍保持原来的t 0 。
如图中曲线HBD ,随时间的推移,由于物体导热温度变化波及范围扩⼤,到某⼀时间后,右侧表⾯温度也逐渐升⾼,如图中曲线HCD 、HE 、HF 。
最后,当时间达到⼀定值后,温度分布保持恒定,如图中曲线HG (若λ=const ,则HG 是直线)。
由此可见,上述⾮稳态导热过程中,存在着右侧⾯参与换热与不参与换热的两个不同阶段。
(1)第⼀阶段(右侧⾯不参与换热)温度分布显现出部分为⾮稳态导热规律控制区和部分为初始温度区的混合分布,即:在此阶段物体温度分布受t 分布的影响较⼤,此阶段称⾮正规状况阶段。
传热学第3章非稳态热传导
•
边长为b的立方体:l
b 6
例题
1.有一直径为5cm的钢球,初始温度为450℃,被突然置 于温度为30℃的空气中。设钢球表面与周围流体间的 总换热系数为24W/(m2.K),
试确定钢球冷却到300℃所需的时间。 (已知钢球的ρ=7753kg/m3, cp=0.48kJ/(kg.K), λ=33W/(m.K))
Φ1
➢ Ф1:从左侧面导入物体的热流量;
➢ Ф2;从右侧面导出的热流量。
0 τ0
Φ2 τ
平板非稳态导热过程中两侧
表面上导热量随时间的变化
3.1.2 导热微分方程解的唯一性定律
非稳态导热问题的数学描述
c p
t
(
2t x 2
2t y 2
2t z 2
)
•
(a)
I.C t(x, y, z,0) f (x, y, z)
分析问题
有一任意形状物体,体积V,表面积A,物性参数
ρ,λ,c为常数。初始温度t0,初始时刻突然置于
温度t∞(恒温)的流体中,
表面传热系数h为常数。 V A
ρ,λ,c
求解
t∞ h
t0
➢ 物体冷却过程中温度随时间的变化规律; ➢ 物体放出的热量。
1. 物体在冷却过程中温度随时间的变化规律
根据能量守恒:
t
用过余温度表示的导热微分方程: τ=0
I.C
B.C
a
2
x 2
(0 x , 0)
0 0 (0 x )
x0
x
0
x
h
x
τ1 τ2 τ3
t∞ -δ 0
用分离变量法求解,直接给出分析解:
ห้องสมุดไป่ตู้t0
第三章_非稳态导热问题的分析解
ρ C pV
初始条件为 令θ =
dT = q vV − σ XS (T 4 − T w4 ) dτ
(a) (b)
T = T0 τ = 0,
qv L σ XT 03 L aτ T V 4 +θw , Fo = 2 , M o = ,N = ,其中, L = 为 4 T0 λ S σ XT 0 L
dθ + M o (θ 4 − N 4 ) = 0 dFo
薄壁物体的温度响应在非稳态导热过程中如果物体内的温度始终是均匀一致的如导热系数很高的薄壁物体或者说当一个物体与周围环境进行热交换时若认为物体内部的温度分布并不重要而只是关心物体的总体温度随着时间的变化如用热电偶测量气流的温度我们常常只关心整个热电偶结点的温度随时间的变化而对于结点内部的温度分布并不重要
∞
r
r
∞
0
0
Bi =
αL λ
L
(3—2)
其中,α 是对流换热系数; L 是物体的特性尺寸,对于平板,即是厚度,对于圆柱体和球, 即是半径; λ 是物体的导热系数。实际上,Biot 数是物体的导热热阻( 换热热阻(
λ
)与表面的对流
1
α
)的比。一般情况下,当 Bi < 0.1 时,导热物体可近似为薄壁。
−
(e)
θ = C 1e
αS τ ρ C pV
(f)
取(d)的特解为 θ = 1 ,所以方程(d)的一般解为
θ = 1 + C 1e
−
αS τ ρ C pV
(g)
根据初始条件(c) ,求得 C1 = −1 ,因此,终解即热电偶结点的温度变化规律为
3
θ = 1 − exp( −
θ
第三章 非稳态导热
2 、无限长圆柱体或球体
(r , ) a hR r r f ( 2 , , ) f ( Fo, Bi, ) 0 R R R
=
m a 1 f1 ( , 2 ) f1 ( , Fo) 0 hR R Bi
r 1 r f2 ( , ) f2 ( , ) m hR R Bi R
习题:3-13,3-15
(5)
3-3 一维非稳态导热的分析解
一、无限大平板的分析解
一块厚为 2 的无限大平板为例,
t ( x, )
1、导热微分方程式及定解条件
t
t ( x, )
导热微分方程式,由式(2-8)得 t 2t 0 )(3-11) a 2 ,( 0 x , x 初始条件:(1)t ( x,0) t0 ,( 0 x )(3-12) 边界条件:(1)t ( x, ) 0 x x0 (3-13)
二、求解一维非稳态导热问题的图线法
诺谟图:(1)按分析解第一项计算绘制的图线 m (0, ) 中心位置温度随时间变化量(x=0时) 0 ( x,0)
( x, ) f ( 1 ) 任意位置与中心位置的温度比值 m (0, ) 式(3-23)与x无关
其解为: e 0
hA exp( ) cV
说明:1) V A 具有长度的量纲,记作 l ,则 hA hV A2 h(V / A) a BiV FoV (3-6) 2 2 cV A cV (V / A)
a 一般地: Bi , Fo 2 l hl
只有两边同为某一常数时,该式才成立
只与 x 有关
分析解为
n n 1
( x , ) t ( x , ) t 0 t0 t
《传热学》第3章_非稳态热传导分析
《传热学》第3章_非稳态热传导分析非稳态热传导分析是传热学中一个重要的研究内容。
在真实的物理系统中,尤其是工程实际中,非稳态热传导过程往往更为常见。
非稳态热传导分析主要研究物体内部温度分布随时间的变化规律,以及热传导过程中的能量交换。
本文将重点介绍非稳态热传导分析的基本原理和方法。
非稳态热传导分析需要考虑时间因素以及物质的热传导性质。
在非稳态热传导过程中,物体内部的温度分布随时间的变化满足热传导方程。
传热方程的一般形式为:∂(ρcT)/∂t=k∇²T+Q其中ρ是物质密度,c是比热容,T是温度,k是热传导系数,∇²是拉普拉斯算子,Q是热源项,即热传导过程中的能量增减。
解决非稳态热传导分析的一般步骤如下:1.建立热传导方程。
根据实际情况,确定适当的坐标系,并根据系统的几何形状和边界条件,建立热传导方程。
2.确定边界条件。
边界条件包括物体表面的温度、热通量以及对流边界等。
根据具体情况,选择适当的边界条件。
3.选择合适的数值方法。
非稳态热传导问题通常需要借助数值方法进行求解。
有限差分法、有限元法、迭代法等都可以应用于非稳态热传导分析,具体选择哪种方法需要根据具体问题的特点进行判断。
4.数值求解。
根据使用的数值方法,将热传导方程离散化,并进行数值求解。
通常需要在计算过程中进行迭代,直到得到满足要求的结果。
5.结果分析和验证。
得到物体内部温度随时间的变化规律后,可以通过实验进行验证。
比较模拟结果与实验结果,判断模拟的准确性。
非稳态热传导分析的典型应用包括热处理过程中的温度变化分析、电子元器件的散热分析、建筑物内部温度分布分析等。
通过对非稳态热传导问题的分析,可以更好地理解和控制物体内部温度分布的变化规律,为实际工程提供指导。
然而,非稳态热传导分析也存在一些挑战和限制。
首先,非稳态热传导分析通常需要考虑物质性质的非线性以及边界条件的复杂性,这增加了问题的难度。
其次,非稳态热传导问题的求解往往需要较长的计算时间和大量的计算资源。
第3章 非稳态导热
解之,得: 2 a 1 2 x, 2sin n x e cos n 0 n 1 n sin n cos n
式中离散值n是下列超越方程的根,称为特征值
tan n
hA d cV 0
hA cV
hA ln 0 cV
e 0
hA exp cV
l=V/A hA h V A hl cV c V A 2 c l 2
将微分方程分离变量并求解得分析解为 : t t0 1 2 u e u2 du erfc
物体内的温度分布 根据半无限大物体的定义,得出导热微 分方程为: 2 a x2 初始条件为: τ=0 时, ( x,0) t0 t0 0 边界条件为:x=0 时, t t
x0 w 0 w
x= ∞ 时,
x t0 t0 0
t 2t a 2 0 x , 0 x t x,0 t0 0 x t x, 0 x x 0 t x, h t x , t x x
对热量计算公式的说明
热量计算公式适用于物体被冷却时,温差取
热量计算公式适用于物体被加热时,温差取
t0 t t t0
物体内部导热热阻可以忽略时的加热或冷却,有时又称 为牛顿加热或牛顿冷却。
注意:由于用集总参数法求物体的温度分布时,认为物 体内没有温度梯度,温度只随时间而变化,所以不能 用傅立叶定律求热量。
中心点的温度
12
Fo
x cos 1
0, 2 sin 1 e 0 1 sin 1 cos 1
《传热学》第三章 非稳态导热
令:
—— 过余温度
使导热微分方程边界条件齐次化:
1.分离变量法求解导热微分方程:
对于此类偏微分方程,应采用分离变量法来进行求解: 假定:
代入导热微分方程,得出:
令:
并对两式分别求解
求解结果: 因φ 不可能是无限大或常数,所以只能有:μ <0,因而可令:
求解结果:
将两个求解结果合并,得到:
其中:
A c1c2 , B c1c3
集总热容体的温度分布:
其中:
L
V ——定型尺寸 A
cV
hA
——时间常数(表示物体温度接近流体温度的快慢)
集总热容体的温度分布亦可写成:
四、不同加热方式下的无限大平壁瞬态导热
t
qv
h, t f
h, t f
qw
qw
h, t f
h, t f
x
第三节 半无限大物体的瞬态导热
应用领域:大地 一、第一类边界条件
半无限大物体表面温度:
半无限大物体表热负荷:
——一定时间内将壁温提高至tw所需的热负荷
第四节 其他形状物体的瞬态导热
一、无限长圆柱体和球体——计算线图法 分无 布限 计长 算圆 步柱 骤温 度
计算Bi和Fo
由图3-13计算中心温度
由图3-14计算任意处温度 无限大平壁—— 半壁厚δ
定型尺寸
无限长圆柱体和球体—— 半径 R 其他不规则形状物体——V/A
或:
傅立叶准则——
二、正常情况阶段——Fo准则对温度分布的影响
对
进行收敛性分析: 随着β n的递增,级数中指数一项收敛很快,所以级数收敛很快,尤其当Fo较 大时,收敛性更加明显。 因此,当Fo>0.2时,仅用级数第一项来描述,已足够精确,即:
第三章--非稳态导热
Q 0 Φ ( )d Vc 0 (1 e
hA Vc
) J
当物体被加热时(t<t),计算式相同(为什 么?)
4 Biv Fov
hl
物理意义
l 物体内部导热热阻 Bi = 1 h 物体表面对流换热热阻
换热时间 Fo 2 l a 边界热扰动扩散到l 2面积上所需的时间
A Qc
ΔΕ
ρ , c, V, t0
α , t0
热平衡关系为:内热能随时间的变化率ΔΕ=通 过表面与外界交换的热流量Qc
当物体被冷却时(t >t),由能量守恒 可知 hA(t t ) - Vc dt
令: t t — 过余温度,则有 :
控制方程 hA - Vc d d ( 0) t t 初始条件 0 0
第三章 非稳态导热
第三章
• • • • • §3.1 §3.2 §3.3 §3.4 §3.5
非稳态导热
非稳态导热的基本概念 集总参数法的简化分析 一维非稳态导热的分析解 二维三维非稳态导热问题的求解 半无限大物体的非稳态导热
1、重点内容: ① 非稳态导热的基本概念及特点; ② 集总参数法的基本原理及应用; ③ 一维及二维非稳态导热问题。 2 、掌握内容: ① 确定瞬时温度场的方法; ② 确定在一时间间隔内物体所传导热量的计算方法。 3 、了解内容:无限大物体非稳态导热的基本特点。
对厚为2δ的
无限大平板 对半径为R的无 限长圆柱 对半径为R的 球
M 1 1 M 2 1 M 3
对于一个复杂形体的形状修正系数时,可以将 修正系数M取为1/3,即 Bi 0.0333
V
例题3-2 将一个初始温度为20℃、直径为100mm的 钢球投入1000℃的加热炉中加热,表面传热系数为 h=50W/(m2· K)。已知钢球的密度为7790kg/m3 , 比热容为470J/(kg· K),导热系数为43.2W/(m· K)。 试求钢球中心温度达到800℃所需要的时间。
《传热学》第3章-非稳态导热
特殊多维非稳态导热的简易求解方法
在第一类边界条件(初始温度均匀)或第三类边界条件(表面 传热系数h为常数)下的二维或三维的非稳态导热问题,在数学 上已经证明,它们的无量纲过余温度的解等于构成这些物体的 两个或三个物体在同样边界条件下一维非稳态导热问题解的连 乘。
特殊多维非稳态导热的简易求解方法
对于无限长方柱 θ (x, y,τ ) = θ (x,τ ) ⋅ θ (y,τ )
该问题的解可以由3块相应的无限大平板的 解得出。最低温度发生在钢锭的中心,即3 筷无限大平板中心截面的交点上,最高温度 发生在钢锭的顶角,即3块大平板表面的公 共点上。
4
例题3 θ
m/B则θi x0钢==锭hλδ(1θ中=m心3/ 4θ温840×0度).05x.2⋅5(θ=
2.14
m/θ 0
)
y
⋅ (θ
无限大平板的非稳态导热
当Fo ≥ 0.2时,可取
θ (x,τ )
θ0
=
β1
2 sin β1 + sin β1 cos β1
cos
β
1
x δ
e − β12 ⋅Fo
只与Bi、x/δ有关, 与时间无关
lnθ
=
−mτ
+ lnθ 0
β1
2sin β1 + sinτ β1 cos β1
cos
= 0.36
短圆柱的中心温度为
查图3-6得 θ
再讨论直径为
m2R/θ=600=0m0m.8的无θ限m长/ θ圆0柱=:0.13
×
0.8
=
0.104
Bi = hR = 232 × 0.3 = 1.72 λ 40.5
tm = 0.104θ0 + t∞ 查附=2图0.11得04θ×m(3/θ00−=103.0103) +1300
第3章-非稳态导热分析解法3
任一点的热流密度:
1 0 qx x a
0 a
e
x
2
4 a
令 x 0 即得边界面上的热流密度
qw
[0,]内累计传热量
q
0
q w d 2
c 0
吸热系数
导热理论分析方法的基本思路
1、根据几何条件、物理条件简化导 热微分方程式 几何条件:导热体的 几何形状和大小,
f (Bi,
x
), f (Bi, Fo)
Bi<0.1,集中参数法
( x, ) (Bi, ) m ( )
x
P130:图3-8
( x, ) ( x, ) m ( ) ; 0 m ( ) 0
f (Bi,
x
) f (Bi, Fo)
m ( ) (Bi, Fo) 0
Q f 3 Bi , Fo Q0
特征尺寸R为圆柱体或球体的半径, r为圆柱体或球体的径向方向。
§3-4 半无限大的物体
• 半无限大系统:一个半无限大 的空间,也就是一个从其表面可 以向其深度方向无限延展的物体 系统。
• 很多实际的物体在加热或冷却过程的初期都可以 视为是一个半无限大固体的非稳态导热过程。
Bi h Bi hR
F0 a F0 a
2
R2
此处的A,B及函数 f(μ1,η) 见P127表3-1
3、 非稳态导热正规状况的工程计算法 (1)拟合公式法
( , ) A exp( 12 Fo) f ( 1 ) 0
Q 1 A exp( 12 Fo) B Q0
x 2 4 a
erf ( 2 ) 0 . 9953 0
第3章-非稳态导热分析解法1
0 100 200 300 400 500
cV
hA
θ/θ0
时间常数τc 小 导热体的温度变化快; 单位时间所传递的热量 大
0.4 0.2 0
时间τ [s]
时间常数τc 小:cV 小、hA大
作业: 3-3 3-6 3-13
(3)求解方法:分析解法、近似分析法、数值解法 分析解法:分离变量法、积分变换、拉普拉斯变换 近似分析法:集中参数法、积分法、瑞利-里兹法 数值解法:有限差分法、蒙特卡洛法、有限元法、 分子动力学模拟
7 毕渥数 Bi
已知:平板厚2δ、初温t0、平板导热系 数λ,将其突然置于温度为 tf的流体中冷 却,表面传热系数为h。 本章以第三类边界条件为重点。 (1) 问题的分析 如图所示,存在两个换热环节: a 流体与物体表面的对流换热环节 b 物体内部的导热
t
x
5 热流量随时间变化
t
Φ1 Φ2 Φ1
Φ2
x
Φ1--平壁左侧导入的热流量 Φ2--平壁右侧导出的热流量
6 学习非稳态导热的目的:
(1) 温度分布和热流量分布随时间和空间的变 化规律 t f ( x, y , z , ) ; Φ f( ) (2) 非稳态导热的导热微分方程式:
t t t t c ( ) ( ) ( ) x x y y z z
§3-2 零维问题的分析法-集中参数法
集中参数法 定义:忽略物体内部导热热阻、认为物体温度 均匀一致的分析方法。此时, Bi 0 ,温度分 布只与时间有关,即 t f ( ) ,与空间位置无 关,因此,也称为零维问题。
适用条件:物体的导热系数相当大,或几何尺寸很 小,或者表面传热系数很低时,符合此类情况
传热学-第三章非稳态导热问题分析解
单位时间 0, t t0
物体内能 的减少(或 增加)
Φ hAt t
Φ cV dt d
当物体被冷却时(t 0 >t),由能量守恒可
知
hA(t t ) -Vc dt
d
令: t t — 过余温度,则有
hA
-Vc
d d
( 0) t0 t 0
控制方程 初始条件
方程式改写为:d hA d 分离变量法 Vc
由于表面对流换热热阻与导热热阻相对大小的不同, 平板中温度场的变化会出现以下三种情形:
(1) 1/ h / Bi
(2) / 1/ h Bi 0
(3) δ/ λ 与1/h 的数值比较接近 0 Bi
Bi 准则对温度分布的影响
1/ h /
/ 1/ h δ/ λ 与1/h的数值接近
是一种理想化模型; 物体内导热热阻忽略不计; 物体内温度梯度忽略不计,认为整个物体具有相
同的温度;
通过表面传递的热量立即使整个物体的温度同时 发生变化; 把一个有分布热容的物体看成是一个集中热容的物体;
只考虑与环境间的换热不考虑物体内的导热。
问题的提出:
2 温度分布 如图所示,任意形状的物体,参数均为已知。
0.049 0.05 可采用集总参数法。
F cp V
cp
dl 2d 2 d 2l 4
4
cp
4(l d dl
2)
140 4 (0.3 0.025) 480 7753 0.05 0.3
0.326102
t tf 800 1200 0.342
0 t0 tf 30 1200
由式(3-1)得:
???
§3-2 集总参数法
基本思想:对任意形状的物体,忽略物体内部的导热 热阻,认为物体温度均匀一致。
传热学课件-第三章非稳态热传导共66页文档
0
物体中的温度 呈指数分布
方程中指数的量纲:
hA
mW2Km2
w1
cV
kg m3
JKkg[m3]
J
s
即与 1 的量纲相同,当
时hVAc,则
hA 1
Vc
此时, e1 36.8% 0
上式表明:当传热时间等于 Vc时,物体的过余温度已经达
到了初始过余温度的36.8%。hA称 为Vc 时间常数,用 表示 。c
非周期性非稳态导热:物体的温度随时间的推移逐渐趋 近于恒定的值
非周期性非稳态导热实例(汽轮机外壳)
冷态启动前:tf1=tw1=tw2=tf2
进汽后 tf1
内壁 q1=h1(tf1-tw1) 到某一时刻 h1A1(tf1-tw1)=h2A2(tw2-tf2) 以后为稳态导热
3 温度分布:
问题描述:
3 、了解内容:
①无限大物体非稳态导热的基本特点。 ②二维非稳态导热问题。
§3-1 非稳态导热的基本概念
一、非稳态导热的特点及类型
1 非稳态导热的定义
物体的温度随时间而变化的导热过程称非稳态导热。
2 非稳态导热的分类
周期性非稳态导热:物体的温度随时间而作周期性的变化
例如太阳辐射的周期性变化引起的房屋的墙壁温度随时间的变化。
这 时 , 由 于 导 热 热 阻 δ/λ几乎可以 忽略,因而任一时刻平板中各点的 温度接近均匀,并随着时间的推移, 整体地下降,逐渐趋近于t∞ 。
(3) δ/λ与 1/h 的数值比较接近
这时,平板中不同时刻的温度分布介于上述两种极 端情况之间。
由此可见,上述两个热阻 的相对大小对于物体中非稳态 导热的温度场的变化具有重要 影响。为此,我们引入表征这 两个热阻比值的无量纲数毕渥 数: Bi h
第三章--非稳态导热
(3) 求解方法:
分析解法、近似分析法、数值解法 分析解法:分离变量法、积分变换、拉普 拉斯变换 近似分析法:集总参数法、积分法 数值解法:有限差分法、蒙特卡洛法、有 限元法、分子动力学模拟
二、讨论物体处于恒温介质中的第三类边界 条件问题
在第三类边界条件下,确定非稳态导热物体中的 温度变化特征与边界条件参数的关系。
e 0
e
Biv Fov
物体中的温度 呈指数分布
方程中指数的量纲:
W 2 m 2 hA w 1 m K kg Jkg Vc J s 3 [m ] 3 K m
即与
1
的量纲相同,当
可解得 Fov=83.6
83.6R 3
R 3
a
2
83.6
2 2
2
c p
83.6 0.05m 3m 1968s 32.8 min 43.2W/(m K) 7790kg/m 3 470J/(kg K)
3
例题3-2 一温度计的水银泡呈圆柱状,长20mm,内径为4mm,初始温 度为t0,今将其插入到温度较高的储气罐中测量气体温度.设水银泡同 气体间的对流换热表面传热系数h=11.63W/(m2.K),水银泡一层薄玻 璃的作用可忽略不计,试计算此条件下温度计的时间常数,并确定插 入5min后温度计读数的过余温度为初始温度的百分之几?水银的物 性参数如下: k g m3 c 0.138k J kg K 10.36W m K 13110 解:首先检验是否可用集总参数法.考虑到水银泡柱体的上端面 不直接受热,故 V R 2l 0.002 0.02 3 0 . 953 10 m 2 A 2Rl R 2 0.20 0.001
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
解的唯一性定理 数学上可以证明,如果某一函数t(x,y,z,τ)满足 数学上可以证明,如果某一函数t(x,y,z,τ)满足 t(x,y,z,τ) 方程( 1a)( 1b)以及一定的初始和边界条件, )(3 方程(3-1a)(3-1b)以及一定的初始和边界条件, 则此函数就是这一特定导热问题的唯一解。 则此函数就是这一特定导热问题的唯一解。 本章所介绍的各种分析法都被认为是满足特定问题 的唯一解。 的唯一解。
非稳态导热问题的求解实质: 非稳态导热问题的求解实质:在规定的初始条件 及边界条件下求解导热微分方程式, 及边界条件下求解导热微分方程式,是本章主要 任务。 任务。 三个不同坐标系下导热微分方程式,用矢量形 三个不同坐标系下导热微分方程式, 式统一表示为: 式统一表示为:
∂t & ρcp = λ div( grad t ) + φ (3-1a) ∂τ
温度的拉普拉斯算子
∇ 2t
& ∂t φ = a∇ 2t + ρcp ∂τ
(3-1b)
初始条件的一般形式 初始条件的一般形式
t ( x, y, z , 0) = f ( x, y, z )
简单特例 f(x,y,z)=t0
边界条件: 边界条件:着重讨论第三类边界条件
∂t −λ ( ) w = h(tw − t f ) ∂n
(1) 1/ h << δ / λ ) 这时, 这时,由于表面对流换 几乎可以忽略, 热热阻1 / h 几乎可以忽略, 因而过程一开始平板的表 t∞ 面温度就被冷却到 。 并随着时间的推移, 并随着时间的推移,整体 t 地下降, 地下降,逐渐趋近于∞ 。
t∞
δ λ δh Bi = = 1h λ
(2) )
方法一
椐非稳态有内热源的导热微分方程: 椐非稳态有内热源的导热微分方程:
λ ∂ 2t ∂ 2t ∂ 2t Φ ∂t + = + + ∂τ ρc ∂x 2 ∂y 2 ∂z 2 ρc
⋅
∵物体内部导热热阻很小,忽略不计。 物体内部导热热阻很小,忽略不计。 物体温度在同一瞬间各点温度基本相等,即t仅是τ 物体温度在同一瞬间各点温度基本相等, 仅是τ 的一元函数,与坐标x 无关, 的一元函数,与坐标x、y、z无关,即
δ / λ << 1/ h
这时, 这时,平板内部导 热热阻 δ / λ 几乎可以忽 略,因而任一时刻平板 中各点的温度接近均匀, 中各点的温度接近均匀, 并随着时间的推移, 并随着时间的推移,整 体地下降, 体地下降,逐渐趋近 于 。t∞
Bi =
δ λ δh = 1h λ
(3)δ / λ 与 1 / h 的数值比较接近 )
θ
0
dθ
⇒
θ hA ln =− τ θ0 ρVc
⇒
θ t −t∞ = =e θ0 t0 −t∞
−
hA τ ρVc
hA hV λ A2 其中的指数: 其中的指数: τ= ⋅ 2 τ ρ cV λ A V ρc h(V A) aτ = ⋅ = Bi ⋅ Fo 2 (V A) λ过Fra bibliotek温度 比 傅立叶数
温度呈指数 分布
∂ 2t ∂ 2t ∂ 2t 2 + 2 + 2 =0 ∂x ∂y ∂z
∂t Φ = ∂τ ρ c
⋅
φ可视为广义热源,而且热交换的边界不是计算边 可视为广义热源, 零维无任何边界) 界(零维无任何边界) 界面上交换的热量应折算成整个物体的体积热源,即: 界面上交换的热量应折算成整个物体的体积热源,
3.2 零维问题的分析法-集总参数法
定义:忽略物体内部导热热阻、 定义:忽略物体内部导热热阻、认为物体温度均 匀一致的分析方法。此时, 匀一致的分析方法。此时,Bi → 0 ,温度分布只 与时间有关,即 t = f (τ ) ,与空间位置无关, 与空间位置无关, 与时间有关, 因此,也称为零维问题。 零维问题 因此,也称为零维问题。
V 特征长度 lc = A
θ = exp(−Bi, Fo) θ0
θ θ0
Bi ⋅ Fo
应用集总参数法时, 应用集总参数法时,物体过余温度随时间的变化 关系是一条负自然指数曲线, 关系是一条负自然指数曲线,或者无因次温度的 对数与时间的关系是一条负斜率直线
导热量计算式、 3.2.2 导热量计算式、时间常数与傅立叶数 1、导热量计算
2
非稳态导热的分类
周期性非稳态导热: 周期性非稳态导热:物体的温度随时间而作周期 性的变化 非周期性非稳态导热(瞬态导热) 非周期性非稳态导热(瞬态导热):物体的温度 随时间不断地升高(加热过程)或降低( 随时间不断地升高(加热过程)或降低(冷却过 ),在经历相当长时间后 在经历相当长时间后, 程),在经历相当长时间后,物体温度逐渐趋近 于周围介质温度,最终达到热平衡。 于周围介质温度,最终达到热平衡。 物体的温度随时间的推移逐渐趋近于恒定的值. 物体的温度随时间的推移逐渐趋近于恒定的值. 着重讨论瞬态非稳态导热。 着重讨论瞬态非稳态导热。
第三类边界条件下Bi Bi数对平板中 3.1.3 第三类边界条件下Bi数对平板中 温度分布的影响
在第三类边界条件下, 在第三类边界条件下,确定非稳态导热物体中的 温度变化特征与边界条件参数的关系。 温度变化特征与边界条件参数的关系。 已知:平板厚 2δ 、初温 t 0 、表面传热系数 h 、 已知: 平板导热系数 ,将其突然置于温度 为t ∞ λ 的流体中冷却。 的流体中冷却。 平板中温度场的变化会出现以下三种情形: 平板中温度场的变化会出现以下三种情形:
这时平板中不同时刻的温度分 布介于上述两种极端情况之间。 布介于上述两种极端情况之间。 由此可见,上述两个热阻的 由此可见, 相对大小对于物体中非稳态导 热的温度场的变化具有重要影 为此, 响。为此,我们引入表征这两 个热阻比值的无量纲数毕渥数 毕渥数。 个热阻比值的无量纲数毕渥数。
Bi =
δ λ δh = 1h λ
适用于本问题的导 热微分方程式
令: = t − t∞ —过余温度 θ
控制方程 hAθ =-ρVc dθ dτ θ (τ = 0) = t − t =θ 初始条件 0 ∞ 0
方程式改写为: 方程式改写为:
dθ
hA =− dτ θ ρVc
dθ
hA =− dτ θ ρVc
积分 ⇒ ⇒
hA τ ∫θ θ = − ρVc ∫0 dτ
第三章 非稳态导热
第三章 非稳态热传导
3.1 非稳态导热的基本概念 3.2 零维问题的分析法-集总参数法 零维问题的分析法- 3.3 典型一维物体非稳态导热的分析解
1、重点内容: 重点内容: 非稳态导热的基本概念及特点; ① 非稳态导热的基本概念及特点; 集总参数法的基本原理及应用; ② 集总参数法的基本原理及应用; 一维非稳态导热问题。 ③ 一维非稳态导热问题。 2 、掌握内容: 掌握内容: 确定瞬时温度场的方法; ① 确定瞬时温度场的方法; ② 确定在一时间间隔内物体所传导热量的计算 方法。 方法。
3.1 非稳态导热的基本概念
3.1.1 非稳态导热过程及其特点
物体的温度随时间而变化的导热过程为非稳态导热。 物体的温度随时间而变化的导热过程为非稳态导热。 温度随时间而变化的导热过程为非稳态导热 自然界和工程上许多导热过程为非稳态,t= f(τ) f(τ 自然界和工程上许多导热过程为非稳态, 冶金、热处理与热加工中工件被加热或冷却; 例:冶金、热处理与热加工中工件被加热或冷却; 锅炉、内燃机等装置起动、停机、变工况; 锅炉、内燃机等装置起动、停机、变工况;自然环 境温度;供暖或停暖过程中墙内与室内空气温度。 境温度;供暖或停暖过程中墙内与室内空气温度。
A
ΔΕ
ρ, c, V, t0
φc
h, t∞
建立数学模型-利用两种方法 建立数学模型-利用两种方法 两种
根据导热微分方程的一般形式进行简化; 根据导热微分方程的一般形式进行简化; 导热微分方程的一般形式进行简化 利用能量守恒 热平衡关系为:内热能随时间的变化率ΔΕ= 热平衡关系为:内热能随时间的变化率ΔΕ=通 ΔΕ 过表面与外界交换的热流量φ 过表面与外界交换的热流量φc 。
− ΦV = Ah(t − t ∞ )
物体被冷却, 物体被冷却,∴φ应为负值
dt ρ cV = − Ah(t − t∞ ) dτ
适用于本问题的导 热微分方程式
方法二
当物体被冷却时( 当物体被冷却时(t >t∞),由能量守恒可知 物体与环境的对流散热量= 物体与环境的对流散热量=物体内能的减少量
dt hA(t − t∞ ) = -ρVc dτ
瞬态热流量:
dt hA Φ = −ρcV = −ρcV (t0 −t∞ ) e ρcV dτ = (t0 − t∞ )hAe
− hA τ ρcV − hA τ ρcV
[ W]
内传给流体的总热量: 导热体在时间 0-τ 内传给流体的总热量:
Q = ∫ Φ(τ )dτ = (t0 − t∞ )∫ hAe τ
5 热量变化
Φ1--板左侧导入的热流量 --板左侧导入的热流量 Φ2--板右侧导出的热流量 --板右侧导出的热流量
各阶段热流量的特征: 各阶段热流量的特征: 非正规状况阶段: 非正规状况阶段:Φ1急剧减小,Φ2保持不变; 非正规状况阶段: 非正规状况阶段: Φ1逐渐减小,Φ2逐渐增大。
3.1.2 导热微分方程解的唯一性定律
正规状况阶段(右侧面参与换热 ): 正规状况阶段( 当右侧面参与换热以后, 当右侧面参与换热以后,物体中的温度分布 不受初始温度的影响, 不受初始温度的影响,主要取决于边界条件及物 此时非稳态导热过程进入到正规状况阶段。 性,此时非稳态导热过程进入到正规状况阶段。 环境的热影响已经扩展到整个物体内部, 环境的热影响已经扩展到整个物体内部,即 物体(或系统) 物体(或系统)不再受到初始温度分布影响的阶 段。可以用初等函数描述。 可以用初等函数描述。 二类非稳态导热的区别: 二类非稳态导热的区别:瞬态导热存在着有区别 的两个不同阶段,而周期性导热不存在。 的两个不同阶段,而周期性导热不存在。