03传热学第三章非稳态热传导

合集下载

第三章 非稳态导热传热学

第三章 非稳态导热传热学
基本思想: 基本思想:当所研究的问题非常复杂, 当所研究的问题非常复杂,涉及到的参数很多, 涉及到的参数很多, 为了减少问题所涉及的参数, 为了减少问题所涉及的参数,于是人们将这样一些参数组合 起来, 起来,使之能表征一类物理现象, 使之能表征一类物理现象,或物理过程的主要特征, 或物理过程的主要特征, 并且没有量纲。 并且没有量纲。因此, 因此,这样的无量纲数又被称为特征数, 这样的无量纲数又被称为特征数,或 者准则数。 者准则数。
§3.1 非稳态导热的基本概念
二、非稳态导热的研究内容
1. 研究内容
温度分布和热流量分布随时间和空间的变化规律
t = f ( x, y , z ,τ ) ;
2. 数学模型
Φ = f(τ )
∂t ∂ ∂t ∂ ∂t ∂ ∂t ɺ ρ c = ( λ ) + ( λ ) + ( λ )+Φ ∂τ ∂x ∂x ∂y ∂y ∂z ∂z 解的唯一性定律 初 始 条 件 边 界 条 件
τ4 τ3
τ2
t
1
τ1
t
0
τ0
第3章 非稳态热传导
§3.1 非稳态导热的基本概念
一、非稳态导热
6. 导热量的特点
Φ1
Φ2
由于物体各处本身温度的变化 要积聚或消耗热量, 要积聚或消耗热量,非稳态导热过 程中在与热流方向相垂直的不同截 面上热流量处处不等。 面上热流量处处不等。
第3章 非稳态热传导
Φ1--板左侧导入的热流量 --板左侧导入的热流量 Φ2--板右侧导出的热流量 --板右侧导出的热流量

t
tf,h x
q
rh
rh = 1 h
rλ = δ λ

3传热学-第三章

3传热学-第三章

2019/2/22
2
3 温度分布
(1)左侧壁面温度突然升高到t1,并保持不变

t
1
4 3
2
1
t
0
0
2019/2/22
3
(2)今有一无限大平板,突 然放入加热炉中加热,平板受 炉内烟气环境的加热作用,其 温度就会从平板表面向平板中 心随时间逐渐升高,其内能也 逐渐增加,同时伴随着热流向 平板中心的传递。
Φ f(x, y, z, )
(2) 非稳态导热的导热微分方程式:
t t t t c ( ) ( ) ( ) x x y y z z
(3) 求解方法: 分析解法、近似分析法、数值解法
2019/2/22
7
7
毕渥数
t
本章以第三类边界条件为重点。 已知:平板厚 2 、初温 t 0 、表面 传热系数 h 、平板导热系数 , 将其突然置于温度为 t 的流体中 冷却。
t c
hA(t t ) V

dt hA(t t ) - Vc d
2019/2/22
初始条件
0, t t0
13
令: t t — 过余温度,则有
控制方程 hA - Vc d d ( 0) t t 初始条件 0 0
16
2019/2/22
hA 即 与 1 的量纲相同,当 Vc 时,则 Vc hA
2019/2/22
10
Bi 准则对温度分布的影响
0
t t0 0
3
2
1
t t0 0
1 0
t t0
1 0
2 1 2 1

传热学第三章 非稳态导热

传热学第三章 非稳态导热
Bi hl ≤0.1
时、物体中最大与最小的过余温度之差小于5%,对于一 般工程计算,此时已经足然特确地可以认为整个物体温度 均匀。按照这样要求,由于l=V/A对圆柱有球分别是半轻 的1/2与1/3、因而如果以l作为Bi数的特征长度,则该Bi数 对平板、国柱与球应该分别小于0.1、0.05和0. 033。
方程中指数的量纲:
hA
W m2K
m2
w1
Vc
kg m3
Jkg K
[
m3
]
J
s
第三章 非稳态导热
9
即与 1 的量纲相同,当 Vc 时,则
hA
hA
1 Vc
此时,
e1 36.8%
0

Vc
hA
为时间常数,用 c 表示。
第三章 非稳态导热
10
如果导热体的热容量( Vc )小、换热条件好(h大),
有一直径为 5cm 的钢球,初始温度为 450 ℃,将其突然置 于温度为 30 ℃空气中。设钢球表面与周围环境间的总换热 系数为 24w/(m2 . K),试计算钢球冷却到 300 ℃所需的 时间。已知钢球的 c=0.48kJ/(kg·K ) , ρ =7753kg/m3 , λ =33w/(m. K ).
Fo
l2
a
换热时间 边界热扰动扩散到l 2面积上所需的时间
无量纲 热阻
Fo越大,热扰动就能越深入地传播到物体 内部,因而,物体各点地温度就越接近周
围介质的温度。
无量纲 时间
第三章 非稳态导热
12
对于平板、圆柱、球的一维非稳态第三类边界条件条件下 的导热问题,当按特征长度
l= 、厚度为2 的平板,
l=R、圆柱 l=R.球 定义的Bi数满足

传热学课件 第三章 非稳定导热

传热学课件 第三章 非稳定导热
2
0e
特征尺寸
V A

hA cV
0e FoBi
a Fo V / A2
Bi
h(V / A)
4. 热量计算

0
hAt t hA hA0e BiF hA0ehA cV f
累积传热量 0 时传热量
Biv

hV A


h d 2 l 4 dl 2 d 2 4




h dl 4 140 0.50 0.3 0.049 0.05 l d 2 33 0.3 0.025
1
可以采用集总参数法.又
hA h V h 4l d 2 cV c A c dl 140 4 0.325 0.326 10 2 s 1 0.48 10 3 7753 0.005 0.3
cV
10.36 5 60 1.89103 0.138103 13110 0.953103 2


exp Bi Fo exp1.07103 1.89103 exp 2.02 0.133 V V
0
即经5min后温度计读数的过余温度的确13.3%.也就是说,在这 段时间内温度计的读数上升了这次测量中温度跃升的86.7%
2. 数学描述
t t a x 2
2
引入过余温度
t t
0 t0 t

t x,0 t0
t x 0, 0 x x 0 t x , ht t x
t
t x x
方程简化为:
例题3-3 一直径为5cm,长为30cm的钢圆柱体,初始温度为300C,将其 放入炉温为12000C的加热炉中加热,升温到8000C方可取出.设钢圆 柱体与烟气间的复合换热表面传热系数为140W/(m2.K),钢的物性参 数取与例3-1中一样的值,问需多少时间才能达到要求。 解:首先检验是否可用集总参数法.为此计算Biv

3第三章 非稳态导热

3第三章 非稳态导热

Bi
n
2.一维非稳态导热的分析解
(2)总传热量
设从初始时刻至某一时刻τ所传递的热量为Q,则有:
分离变量积分并代入初始条件得:
hA
=e cV
0
思考:上述结果是对物体被冷却 的情况导出的,如果要用于被加 热的场合,该怎么办?
6.集总参数系统的分析解
hA hV cV A
A2 cV 2
h(V / A) a (V / A)2
BiV FoV
Bi hl l= 物体内部导热热阻 1 h 物体表面对流换热热阻
• 在某厂生产的测温元件说明书上,标明该元件的 时间常数为1s。你怎么看待这个值?
cV
c hA
——根据定义式,时间常数中物性参数ρ、c、V、A可 以看作是常数,但表面传热系数h却是与具体过程 有关的量。
——说明书上的标明的时间常数需要具体分析,不能 盲目相信。
【内容小结】
• 集总参数系统的分析 • 时间常数的导出和意义 • 时间常数对测温系统的指导
一个集总参数系统,其体积
为V、表面积为A、密度为、 比热为c、初始温度为t0,突 然放入温度为tf (设t0> tf )、 对流换热系数为h的环境中,
求系统温度变化。
A h, tf
ΔE
Qc
ρ, c, V, t0
——表面对流换热对其过程有着重要影响,如何处理?
4. 微分方程
-
t n
ht
t
f
集总参数系统内部没有温差, 不能用第三类边界条件。
不断减小,在其它各截面上,其
截面温度开始升高之前通过该截
面的热流量是零,温度开始升高
A
之后,热流量才开始增加。
BC D 3

《传热学》第3章_非稳态热传导分析

《传热学》第3章_非稳态热传导分析
9
第3章 非稳态热传导
3.2 零维问题的分析法——集中参数法
定义: 固体内部的导热热阻远小于其表面的换热热阻时,任何时刻 固体内部的温度都趋于一致,以致可以认为整个固体在同一 瞬间均处于同一温度下。(类比质心) 特点: 温度仅是时间τ的一元函数 Bi ,而与空间坐标无关 0 t f (。 ) 应用: 物体导热系数极大;几何尺寸极小;表面传热系数极低。 方法: 忽略物体内部导热热阻的简化方法——集中参数法
p
笛卡尔、圆柱和圆球坐标系下的导热微分方程可以统一表示为 t cP div grad t
初始条件 t x, y, z,0 f x, y, z
t x, y, z,0 t0
(第三类边界条件,n为换热表面外法 t 边界条件 ( n ) w h(t w t f ) t 线,h,tf已知,tw, 未知) n w


由此解得
570 s 0.158 h
17
第3章 非稳态热传导
3.3 典型一维物体非稳态导热的分析解
内容:介绍平板、圆柱与球的一维非稳态导热温度场的分析解。
重点:分析解的应用,了解求解过程,掌握公式中各部分表示的含义
3.3.1 三种几何形状物体的温度场分析解 1. 平板
厚度为2δ的无限大平板, 初始温度为t0,外部温度为t∞, 平板关于中心截面对称,因此 只研究x≥0的半块平板
导热微分方程变为
导热微分方程
零维问题完整数学描述
初始条件
d
hA d Vc
11
第3章 非稳态热传导
导热微分方程求解:
d
从0到τ积分 hA d Vc


0
d

第三章非稳态导热_传热学


tm m tf 9 8= 17c C
m 2 0.9064 exp 1.1347 2 0.22 0.9 0 1.1347 0.9064 0.4224
平壁表面处 x 的过余温度为:
w 2sin 1 cos 1 exp 12 Fo 0 1 sin 1 cos 1
(2)在垂直于热量传递的方向上,每个截面上热流量不相等; (3)温度随时间变化,热流也随时间变化。
3.讨论非稳态导热问题的目的:
(1)在加热和冷却时,物体内部某一点温度达到预定温度 时所需要的时间,以及该时间内物体吸收和放出的热量;
(2)对物体加热或冷却之后,物体内部温度分布以及物体 温度随时间的变化率
1 0
1 d a d
c1 exp a
c1 exp a
2
1 d2 X 2 2 X dx
X c2 cos x c3 sin x
x, X x
第三章 非稳态导热
§3-1 §3-2 §3-3 §3-4 §3-5 非稳态导热的基本概念 无限大平壁的瞬态导热 半无限大物体的瞬态导热 其他形状物体的瞬态导热 周期性非稳态导热(自学)
• • • •
1.加热冷却过程 2.动力机械中的开关车
应用背景
3.地球的气候变化
4.医疗中激光技术(控制温度范围)
x, 2sin 1 x cos 1 exp 12 Fo 0 1 sin 1 cos 1
查表3-1,当Bi=2.5时, 1 1.1347
180 sin 1 sin 1.1347 0.9064
x A sin 0 B cos 0 exp a 2 0

《传热学》第3章_非稳态热传导分析

《传热学》第3章_非稳态热传导分析非稳态热传导分析是传热学中一个重要的研究内容。

在真实的物理系统中,尤其是工程实际中,非稳态热传导过程往往更为常见。

非稳态热传导分析主要研究物体内部温度分布随时间的变化规律,以及热传导过程中的能量交换。

本文将重点介绍非稳态热传导分析的基本原理和方法。

非稳态热传导分析需要考虑时间因素以及物质的热传导性质。

在非稳态热传导过程中,物体内部的温度分布随时间的变化满足热传导方程。

传热方程的一般形式为:∂(ρcT)/∂t=k∇²T+Q其中ρ是物质密度,c是比热容,T是温度,k是热传导系数,∇²是拉普拉斯算子,Q是热源项,即热传导过程中的能量增减。

解决非稳态热传导分析的一般步骤如下:1.建立热传导方程。

根据实际情况,确定适当的坐标系,并根据系统的几何形状和边界条件,建立热传导方程。

2.确定边界条件。

边界条件包括物体表面的温度、热通量以及对流边界等。

根据具体情况,选择适当的边界条件。

3.选择合适的数值方法。

非稳态热传导问题通常需要借助数值方法进行求解。

有限差分法、有限元法、迭代法等都可以应用于非稳态热传导分析,具体选择哪种方法需要根据具体问题的特点进行判断。

4.数值求解。

根据使用的数值方法,将热传导方程离散化,并进行数值求解。

通常需要在计算过程中进行迭代,直到得到满足要求的结果。

5.结果分析和验证。

得到物体内部温度随时间的变化规律后,可以通过实验进行验证。

比较模拟结果与实验结果,判断模拟的准确性。

非稳态热传导分析的典型应用包括热处理过程中的温度变化分析、电子元器件的散热分析、建筑物内部温度分布分析等。

通过对非稳态热传导问题的分析,可以更好地理解和控制物体内部温度分布的变化规律,为实际工程提供指导。

然而,非稳态热传导分析也存在一些挑战和限制。

首先,非稳态热传导分析通常需要考虑物质性质的非线性以及边界条件的复杂性,这增加了问题的难度。

其次,非稳态热传导问题的求解往往需要较长的计算时间和大量的计算资源。

《传热学》第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章_非稳态热传导



0
hAexp

hA
cV

d

t0

t
cV
1

exp

hA
cV


2. 时间常数
导热量
过余温度随时间呈指数曲线关系,即:开始温度变化快,后来变化较慢
0

t t t0 t
exp
hA
Vc


cV /hA
热扰动就能越深入地传播到物体内部。
14
第3章 非稳态热传导
3.2.3 集中参数法的适用范围及应用举例
毕渥数越小,越适合应用集中参数法,小到什么程度?
如果特征长度
l ,
l R,
厚度为2的平板
圆柱
那么Bi hl

0.1最大最小
l R, 圆球
过余温度之差小于5%。
如果特征长度lc

第3章 非稳态热传导
3-1 非稳态导热基本概念 3-2 零维问题的分析法——集中参数法 3-3 典型一维非稳态导热的分析解 3-4 半无限大物体的非稳态导热 3-5 简单几何形状物体多维非稳态导热的分析解
1
第3章 非稳态热传导
3.1 非稳态导热的基本概念
定义:物体的温度随时间而变化的导热过程称为非稳态导热 分类:物体温度随时间推移逐渐趋近于恒定的值(金属热处理)
3
第3章 非稳态热传导
非稳态导热问题定性分析:
左金属、右保温,初始t0,左边升高至t1, 温度变化情况。
第一阶段:P-B-L
第二阶段:P-D-L
第三阶段:P-E-L
第四阶段:P-H-M (HM斜率大于PH斜率)
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

cV
dt
d
cV (t0
t )(
hA
cV
)
exp(
hA
cV
)
hA0
exp(
hA
cV
)
※0~ 时间内传给流体的总热量:
Q 0 d
0
h
A
0
e
xp(
hA
cV
)d
2021/1/14
0 cV
1
exp
hA
cV
15
(2) 时间常数

c
cV
hA
ቤተ መጻሕፍቲ ባይዱ
e c
0
※当 时
0 0
即t t
※当

c
与几何参数、物理性 质、换热条件有关
(, ) m ( )
cos(1)
f
( Bi , )
则平板中任意点过余温度比 m 0 m 0
2021/1/14
31
相当于第一 类边界条件
2021/1/14
32
2021/1/14
任意时刻平板 内温度均匀
33
书中的诺谟图仅适用一维平板第一类边界条件下的加热及冷却
过程以及具有恒温介质的第三类边界条件,并且Fo>0.2
Q0
cV (t0 t )
0
τ时刻的平均 过余温度
当Fo>0.2时,正规状况阶段温度场与导热量的计算式可统一表示为:
( , 0
)
A exp(
12 Fo)
f
( 1 )
Q Q0
1
A exp(12Fo)B
其中,A、B、f(μ1η)的表达示见表3-1。
2021/1/14
30
3、非稳态导热正规状况阶段的工程计算方法
2021/1/14
(, ) (0, )
(, ) m ( )
cos( 1 )
与时间无关
29
(3)一段时间间隔内所传导的热量计算式
从初始到平衡所传递 的热量(最大传热量)
考察热量的传递 Q0 cV (t0 t )
若令Q为从初始到某一时刻τ内所传递热量,则
Q cV [t0 t(x, )]dV 1
因此,这样的无量纲数又被称为特征数,或者准则数, 比如,毕渥数又称毕渥准则。以后会陆续遇到许多类似 的准则数。特征数涉及到的几何尺度称为特征长度,一
般用符号 lc 表示。
对于一个特征数,应该掌握其定义式+物理意义, 以及定义式中各个参数的意义。
2021/1/14
12
§3.2 零维问题的分析法—集中参数法
exp(1) 0.368
0
2021/1/14
16
※ hA hV A2 cV A cV 2
hV
A
a
V A2
Biv Fov
方程中指数的量纲:
hA
Vc
W m2K
m2
kg m3
Jkg K
[
m3
]
w J
1 s
物体中的温度
呈指数分布
2021/1/14
0
t t t0 t
exp(Biv Fov )
w
h
tw
tf
若某函数 t(x, y, z, )满足上述方程及初始、边界条
件,此函数为唯一解。
2021/1/14
9
3、第三类边界条件下Bi数对平板中温度分布的影响
物体内部导热热阻:
表面对流换热热阻:1 h
特性长度 特征数 (准则数)
2021/1/14
毕渥数:Biv
1
h
h
10
0
21 3
t t0 0
的表面传热系数为100W/(m2·K)。
鸡肉丝的 810kg / m3、c 3.35kJ /(kg K )、 1.1W /(m K )。
2021/1/14
20
解:∵ Biv
hl
100 2 103 1.1
/4
0.045
0.05
可以使用集中参数法

h
cl
100 810 3.35103 2 103 / 4
17
(3)毕渥数与傅立叶数的物理意义
Bi hlc lc 1h
固体内部单位面积导热热阻 固体表面单位面积对流换热热阻
无量纲 热阻
Fo lc2 a
换热时间 边界热扰动扩散到lc2面积上所需要的时间
无量纲 时间
lc
V A
特征 长度
Fo越大,热扰动就能越深入地传播到物体内部, 因而,物体各点温度就越接近周围介质的温度。
—图线法(诺谟图法)
以无限大平板为例,Fo>0.2 时,取其级数首项即可
e (,
)
0
1
2 sin 1 sin 1 cos
1
12F0 cos(1 ) f (Fo, Bi, )
三个变量,因此,需要分开来画
根据简化表达式先绘出
再根据任意点与中心处 过余温度之比绘出
m ( ) f (Fo, Bi ) 0
n 为方程 1 n cos(n ) Bi , n 1,2, 的根
2021/1/14
27
2、非稳态导热正规状况阶段分析解的简化
(1)非稳态导热正规阶段的物理概念与数学含义
三种几何形状分析解均有 n f (Bi)
当Bi一定,n↑, n↑迅速,
则分析解中时间影响部分
exp(
2 n
Fo)
随Fo增加而迅速衰减。
2021/1/14
35
§3.4 半无限大物体的非稳态导热
主要内容: ※半无限大 物体的定义 难点
( definition of semi-infinite body )
※半无限大物体在非稳态导热中的温度分布
( distribution of temperature in semi-infinite
0
th(
m
H
)
3
实际散热量
7. f 假设整个肋表面处于肋基温度下的散热量
8. 肋面总效率
0
Ar f Af
Ar Af
2021/1/14
4
第三章 非稳态热传导
2021/1/14
5
主要内容:
非稳态导热过程中温度场的变化规律及换热量的分析 求解方法。包括:
1. 零维非稳态导热的集中参数分析法;
Fo>0.2,只保留第1项与完整级数误差小于1%
相当于Cn=0(n≥2) 则三个分析解无穷级数的第1项即为正规状况阶段温度场的解。
t(, ) t f (Fo, Bi,)
0
t0 t
2021/1/14
28
(2)正规状况阶段三个分析解的简化表达式
平板
( , 0
)
1
2 sin 1 sin 1 cos
(, 0
)
n1
Cn
exp(
2 n
FO
)
cos(
n
)
其中,
FO
a 2
,
x
Cn
n
2 sin n cos n sin
n
2021/1/14
n 为方程
tan
n
Bi
n
,n
1,2,
的根
25
(2)圆柱 无量纲温度的分析解为
(, )
0
n1
Cn
exp( n2FO )J0 (n)
其中,
FO
a
R2
,
r R
Cn
1
cos(1
x
) exp( 12Fo)
圆柱
(, ) 0
2
1
J
2 0
(
J1
1 )
( 1 )
J12
(1
)
exp(
12
Fo)
J
0
(
1
)

( , 0
)
2
sin( 1) 1 cos(1) 1 sin( 1) cos(1)
exp( 12Fo)
sin( 1) 1
以平板为例,正规状况阶段的任意时刻,任意处( η )与 平板中心处( η=0.0)处的过余温度之比为
03传热学第三章非稳态热传导
3.边界条件 第一类边界条件:给定壁面温度 tw 第二类边界条件:给定壁面热流密度 qw
第三类边界条件:对流换热条件 h t f
4.无内热源平壁导热(单层和多层)
q (t1 t2 )
(单层)
q
tw1
tw(n1)
n i
(多层)
i1 i
2021/1/14
2
5.圆筒壁导热(单层和多层)
始状态)。
t
t1
t1
初温t0
上升到tw
被 测 t2 介 质
热探针
2021/1/14
满足半无限大物体要求
t2
t0
end
38
1、三种边界条件下半无限大物体温度场的分析解
c t
x
t x
y
t y
z
t

z
控制方程 t a 2t ,0 x
2021/1/14
2
n
J1(n )
J
2 0
(n
)
J12
(n
)
n 为方程
n
J1(n ) J0 (n )
Bi , n
1,2,
的根
2021/1/14
26
(3)球 无量纲温度的分析解为
(,
0
)
n1
Cn
exp(
2 n
FO
)
1
n
sin(
相关文档
最新文档