第4章 非稳态导热
简述非稳态导热的基本特点
非稳态导热的基本特点简介非稳态导热是指热量在物体中传播过程中,温度场随时间发生变化的现象。
相比稳态导热,非稳态导热的特点更加复杂和动态。
基本概念在进一步探讨非稳态导热的基本特点之前,我们先来了解一些相关的基本概念:1.热传导:指物质内部由高温区向低温区传递热量的现象,热能以分子碰撞的方式传导。
2.热扩散:是指物质内部温度差引起的分子的热运动,导致热量向周围传播的过程。
3.热传导方程:描述了非稳态导热的基本规律,形式为∂T(x,t)∂t =α∂2T(x,t)∂x2,其中T(x,t)为温度场,x为空间坐标,t为时间,α为热扩散系数。
非稳态导热的主要特点1. 温度场随时间变化在稳态导热中,物体的温度场是随空间坐标变化而稳定的;而在非稳态导热中,温度场不仅随空间坐标变化,还会随时间发生变化。
这意味着物体的温度分布在传播过程中会发生改变。
2. 传热时间依赖非稳态导热过程中,传热时间不再是一个常数,而是随着时间的推移而改变。
不同部位的温度差越大,传热时间会越短。
这是因为温度差越大,传热速率越快。
3. 初始和边界条件的影响非稳态导热过程中,初始和边界条件起着重要的作用。
初始条件是指导热开始时物体内部温度分布的初始状态,边界条件是指物体表面与外界的热交换情况。
不同的初始条件和边界条件会引起不同的非稳态导热行为。
4. 温度传播的迟滞效应相比稳态导热,在非稳态导热中温度的传播速度较慢。
这是由于热传导所需要的时间,以及物质自身的热容和热传导性质等因素共同作用的结果。
5. 温度波动的存在在非稳态导热过程中,温度场会出现波动现象。
这是由于能量在物体内部的传递是通过分子的热运动实现的,而分子的热运动是随机的,因此温度场会存在一定的波动性。
6. 热传导方程的求解非稳态导热过程可以通过求解热传导方程来描述和预测。
热传导方程是一个偏微分方程,可以通过数值方法或解析方法进行求解。
求解热传导方程可以得到物体中温度场随时间演变的规律。
第四章传热学
4. 非稳态导热4.1 知识结构1. 非稳态导热的特点;2. (恒温介质、第三类边界条件)一维分析解求解方法(分离变量,特解叠加)及解的形式(无穷级数求和);3. 解的准则方程形式,各准则(无量纲过余温度、无量纲尺度、傅里叶准则、毕渥准则)的定义式及其物理涵义; 4. 查诺谟图求解方法;5. 多维问题的解(几个一维问题解(无量纲过余温度)的乘积);6. 集总参数法应用的条件和解的形式;7. 半无限大物体的非稳态导热。
4.2 重点内容剖析4.2.1 概述在设备启动、停车、或间歇运行等过程中,温度场随时间发生变化,热流也随时间发生变化,这样的过程称为非稳态导热。
一.过程特点分类1. 周期性非稳态导热(比较复杂,本书不做研究) 如地球表面受日照的情况 (周期为24小时)对于内燃机气缸壁受燃气冲刷的情况,周期为几分之一秒,温度波动只在很浅的表层,一般作为稳态处理。
2. 非周期性非稳态导热:(趋于稳态的过程,非稳态 稳态) 例子:如图4-1,一个无限大平板,初始温度均匀,某一时刻左壁面突然受到一恒温热源的加热,分析平壁内非稳态温度场的变化过程: (1) 存在两个阶段初始阶段:温度变化到达右壁面之前(如曲线A-C-D ),右侧不参与换热,此时物体内分为两个区间,非稳态导热规律控制区A-C 和初始温度区C-D 。
正规状况阶段:温度变化到达右壁面之后,右侧参与换热,初始温度分布的tx1t 0t ABCDEF图4-1 非稳态导热过程的温度变化影响逐渐消失。
(2) 热流方向上热流量处处不等因为物体各处温度随时间变化而引起内能的变化,在热量传递路径中,一部分热量要用于(或来源于)这些内能,所以热流方向上的热流量处处不等。
二. 研究任务1. 确定物体内部某点达到预定温度所需时间以及该期间所需供给或取走的热量,以便合理拟定加热和冷却的工艺条件,正确选择传热工质;2. 计算某一时刻物体内的温度场及温度场随时间和空间的变化率,以便校核部件所承受的热应力,并根据它制定热工设备的快速启动与安全操作规程。
传热学-第4章-非稳态导热的计算与分析
10
4.2 对流边界条件下的一维非稳态导热
❖ 对几何形状简单、边界条件不太复杂的情形,仍然可 以通过数学分析的方法获得分析解
❖ 这里以(无限大)平壁被流体对称加热的非稳态导热 过程为例,说明非稳态导热的基本特征、分析方法和 过程
❖ 定性地、定量两个方面
11
4.2.1 平壁内非稳态过程的基本特征
问题描述: ❖ 厚为2δ、无内热源的常物性平壁 ❖ 初始时刻温度分布均匀,为t0 ❖ 某时刻突然投入到温度为t∞的高
conduction):物体内任意位置的温度随时间持续升高 (加热过程)或连续下降(冷却过程) 边界条件或内热源不变时,过程将最终逐渐趋于某个 新的稳定温度场
6
4.1 概述
研究目的:
❖ ——确定非稳态过程中的温度场:在此基础上确定物体中
某个部位到达某个预定温度所需经历的时间,或者在预定时间 内可以达到的温度,或者物体的温度对时间的变化速率。
8
4.1 概述
研究方法与过程:与稳态导热的完全相同 (1)简化假设给出物理模型 (2)给出数学模型(方程+定解条件) (3)采用适当的数学方法求解 (4)分析讨论
9
4.1 概述
❖ 非稳态导热的控制方程:
τ
ρct
x
λ
t x
y
λ
t y
z
λ
t z
Φ
❖ t=f(x,y,z,t)
❖ 控制方程:偏微分方程,数学求解难度很大
❖ 随着时间的延续,壁面加热的波及区域将继续向平壁中
心推进
16
4.2.1 平壁内非稳态过程的基本特征
17
4.2.1 平壁内非稳态过程的基本特征
❖ 当温度扰动刚刚传到平壁对称 面的那个时刻,称为穿透时间, 记作τc
第四章 非稳态导热(5)14
④ 某一时刻物体表面的热流量或从某一时刻起经一定时间后表面传递的总热量。
5
2)求解方法:主要有分析解法、数值解法、图解法和热电模拟法等。 本章仅介绍分析解法,而且只针对第三类B、C下一维非稳态导热的求解。
二、一维非稳态导热的分析解及诺谟图
工程上常见的非稳态导热问题分以下三种:
一维非稳态导热问题:
无限大平壁 无限长圆柱体
一、概 述
1.1 定义:非稳态导热是指发生在非稳态温度场内的导热过程。
其数学表达式为:t f (x, y, z, )
按照其过程进行的特点,可分为以下二种:
(1)周期性非稳态导热:导热物体内的温度随时间周期性地 变化。
(2)非周期性非稳态导热(瞬态导热):物体内的温度随时 间不断的 升高或降低。
2
1.2 非稳态导热过程的特点
大平壁非稳态导热分析
由左侧表面导入的热量到达右侧表面之前的一段时间。
② 正常情况阶段。
当左侧表面导入的热量到达右侧表面之后,使右侧壁温不断升高,直到它达
到新的平衡状态的这段时间。
4
B)大平壁两侧被加热过程
一初始温度均匀为t0的无限大平壁,突 然投入到温度为tf的热流体中对称加热。平 壁内发生了非稳态导热过程。平壁中的温
1.3 求解的目的和方法
1) 求解非稳态导热问题主要目的有四个:
① 物体的某一部分从初始温度上升或下降到某一确定温度所需的时间,或经某 一时间后物体各部分的温度是否上升或下降到某一指定值;
② 物体在非稳态导热过程中温度分布,为求材料热应力和热变形提供必要资料; ③ 物体在非稳态导热过程中的温升速率;
二维非稳态导热问题:短圆柱体、长的方柱体
三维非稳态导热问题:短方柱体、长方体
第四章 非稳态导热..
工程中:
机器启动、停机、变工况时部件的导热过程; 冶金、热加工、热处理工艺中工件的加热及冷却过程等; 石油工程中钻井、焖井、采油等过程中热量在地层内的扩散过程。
具有实际意义。
2
第三节
本节讨论:
非稳态导热
——基本概念和特点
——非稳态导热问题的求解及诺模图
——集总参数法
3
m (0, ) f1 ( Fo,Bi) 0
( x, ) ( x, ) m 0 m 0
意味着初始条件的影响已经消失, 这是正规状况阶段。
( x, ) x f 2 ( Bi, ) m (0, )
工程上常采用两种简化的计算方法:
诺模图方法——由海斯勒(Heisler)提出;
13
第四章 / 第三节 非稳态导热
(一)无限大平壁的分析解及诺模图
1、平壁内温度分布的求解
t 2t a 2 0 x , 0 x
初始条件: t | 0 t 0
0 x
边界条件: t | 0 (对称性) x 0
x
t |x h t |x t f x
物体内部各点在同一时刻的温度趋于一 致,温度场与空间位置无关,只是时间 的单值函数。
这样的物体称为集总热容系统。
工程中取Bi<0.1
25
第四章 / 第三节 非稳态导热
(一)无限大平壁的分析解及诺模图
2、Fo数和Bi数的物理意义以及对非稳态过程的影响
h Bi数的影响: Bi 1h
—出现在特征数中的几何尺度 —不同情况下,不同形状的物体特征长度是不同的。 Fo 数 、 Bi数称为特征数,习惯上又称准则数, 具有特定的物理意义。
第四章 非稳态导热(6)14
(b)
36 .8%
可以得出内部热阻可被忽略的非稳 态导热过程具有以下二个特点: (1)物体温度 随时间按指数函数关系下降,如 图所示,开始下降快,随后变化减慢。
0
, t t f 0,即t t f
Tτ
集总参数系统θ -τ曲线
τ
(2)物体温度随周围流体温度变化的快慢与该物体的时间常数Tτ有关。 什么是时间常数?式(b)中 ρcV/(hA) 具有时间的量纲,此外,对于常物 性物体,一旦几何尺寸确定( V/A 确定), ρcV/(hA) 的值也就确定了。 cV T 在以上二个意义上,把 ρcV/(hA) 称为时间常数,记为Tτ,即 。 hA
代人(a)式得 cV
Ah(t t f ) V
集总参数系统的微分方程
dt = Ah (t t f ) d
(2)根据能量守恒定律:物体内能(焓)的变化等于物体表面对外散去的热量:
cV
dt =Ah (t t f ) d
3
求解微分方程:
引入过余温
初始条件:
t tf
d = Ah , 上式变成 cV d
o
d
hAo (
cV
hA
o
)(e
hA cV
hA cV
o
1)
Φ的单位—W或kW; Qτ的单位—J或kJ。
cVo (1 e
)
请大家思考:瞬时的传热量Φ和总传热量Qτ的单位是什么?
7
三、集总参数法的适用条件
集总参数法比较简单,但应用它是有条件的,必须满足: Bi
1 R R 2l V BiV Bi 0.05 L 2 2 2Rl A
传热学第四章非稳态导热例题
(V / A)
3
85 K) 3.885 10 39.63 W/(m2· 0.025 / 3
BiV FoV 3.885 10 2.07945
3
535.25
2013-9-10
9
由式(4-6)计算换热量:
hA Q cV(t 0 t f)1 exp( ) cV
a 6
D(t 0 t f)1 e (
3
BiV FoV
)
85 3 2.07945 0.05 ( 60 (1 e 300 ) ) 5 6 2.95 10
=39.6 kJ
返回
2013-9-10 10
【例4-3】一根直径为1m,壁厚40mm 的钢管,初温为-20℃,后将温度为60℃的 热油泵入管中,油与管壁的换热系数为 500 W/(m2· K),管子外表面可近似认为是绝 热的。管壁的物性参数ρ=7823kg/m3, c=434J/(kg·K),λ=63.9 W/(m· K)。
1.882 10 8 60 Fo 2 5.646 2 0.04
a
5
2013-9-10
14
(2) 由于Bi>0.1, 故不能采用集总参
数法,需用线算图求解。
管子外表面, 1 3.195
Bi
查图4-7得
m 0.24 0
管子外表面温度为:
t m m t f 0.24 0 t f 0.24 20 60 60 40.8 ( )
V 准则中的特征尺寸是用 LV 确定的, A
而不是 R/2 ,所以,是否可采用集总参 数法的判别用BiV<0.1M。
传热学基础(第二版)第四章教学课件非稳态导热
23/250291/4/16
0~τ范围内积分,得凝固层厚度的表达式
2 b L t w c ttp 0tw K
此式称为平方根定律,即凝固层厚度与凝固时 间的平方根成正比。式中
K2 b L t w c ttp 0tw
ms12
K 称为 凝固系数
24/250291/4/16
几种材质在不同冷却条件下的K值
由于砂型的导热系数较小,型壁较厚,所以平面 砂型壁可按半无限大平壁处理。本节得到的公式 应用于铸造工艺,可以计算砂型中特定地点在τ 时刻达到的温度和0~τ时间内传入砂型的累积热量。 瞬时热流密度qw和累计热量Q w都与蓄热系数成正 比,所以选择不同造型材料,即改变蓄热系数, 就成为控制凝固进程和铸件质量的重要手段。
物性的这种组合可表成: a c
cb W /m (2Cs1/2)
a b称为蓄热系数。它完全由材料的热物性构 成,它综合地反映了材料的蓄热能力,也是个热 物性。
15/250291/4/16
铸铁和铸型蓄热系数b的参考值。
热物性 材料
铸铁
导热系数 比热容 密度 热扩散率 蓄热系数
λ
c
ρ
a
b
46.5 753.6 7000 8.82×10-6 15600
5 /59 2021/4/16
积蓄(或放出)热 量随时间而变化是过 程的又一个特点。于 是在工程计算中,确 定瞬时热流密度和累 计热量也是非稳态导 热问题求解的任务。 在图中,累计热量由 指定时间τ与纵坐标 间曲线下的面积表示。
6/59 2021/4/16
4-2 第一类边界条件下的一维非稳态导热
式:
qw ' Lctptw
d d
与式
第4章-非稳态导热的计算分析
是与物体几何形状
Biv
h( V
A)
1、非稳态导热的分类
周期性非稳态导热:物体的温度随时间而作周期 性的变化 非周期性非稳态导热(瞬态导热):物体的温度 随时间不断地升高(加热过程)或降低(冷却过 程),在经历相当长时间后,物体温度逐渐趋近 于周围介质温度,最终达到热平衡,物体的温度 随时间的推移逐渐趋近于恒定的值。
❖ 300℃的铁块在冷水中的冷却
x, 0,
cos
1
x
它表明:当Fo>0.2后,虽然θ(x,τ)与θ(0,τ)各自均与τ相关, 但它们的比值却与τ无关而仅取决于平壁的几何位置(x/δ) 和Bi数
这意味着初始条件的影响已经消失,这就是正规状况阶段
❖ 计算正规状况阶段的温度需要根据Bi数确定相应 的特征值,使用时不甚方便
❖ 工程上常采用两种简化的计算方法,由海斯勒 (Heisler)提出的诺模图(nomogram)方法和由 Campo提出的近似拟合公式
数时,即 τ=τr,
=e1 0.386 0
0.386 01
τ/τr
τ=4τr,
=e4.6 0.01 工程上认为 =4τr时导热
0
体已达到热平衡状态
瞬态热流量:
Φ( ) hA(t( ) t ) hA
总热量:
hA
hA0e Vc
W
导热体在时间 0~ 内传给流体的总热量:
Q
0
Φ(
)d
一、无限大平板加热(冷却)过程分析
厚度 2 的无限大平壁,、a 为已知常数;=0时温度为 t0;
突然把两侧介质温度降低为 t 并保持不变;壁表面与介质之 间的表面传热系数为h。 两侧冷却情况相同、温度分布 对称。中心为原点。
非稳态导热PPT课件
h
ctg
上式中:为书写方便,令=,只要求出即可得,另上
式左侧的分母即为毕渥准则数 Bi=h/ 上式可写为:
/Bi=ctg
<13>
式<13>为一特征方程(超越方程),其解理论上有无穷
多个,如图3-4所示,为 y1= /Bi 与 y2=ctg 间的交点, 其具体数值如P57表3-1所示。于是也有无穷多个解
p116图34见图34a由于表面换热热阻可以忽略一开始平板表面温度就被冷却到随着时间的推移平板内各点的温度逐渐下降而趋近见图34b由于平板导热热阻可以忽略任一时刻各点的温度一致即tf并随时间的推移整体下降逐渐趋近于当两种热阻的数值比较接近即bi为有限值时其温度分布见图34c
第一节 非稳态导热的基本概念
t
2t 2t 2t
( )
x2 y2 z2 t(x,y,z,0)t0
c
nt wh(twtf)
数学上可以证明其解t=f(x,y,z,τ)是唯一的。
.
6
第一节 非稳态导热的基本概念
4、非稳态导热的三种情形
流设体一中块冷厚却2δ,的表金面属换平热板系,数初为始h温,度平为板的t 0,导突热然系将数它为置λ。于根温据度平为板的t的f
1=1/、 2=2/、 3=3/. 、… n=n/ 即:=F(B1i)6
第三节 典型一维物体非稳态导热
一、无限大平板非稳态导热
2.解:③解微分方程
于是,在给定Bi条件下,将对应的1 、2、…、n代入式 <12>,即得一组温度分布:
1(x,)=D1cos(1x)exp(2(x,)=D2cos(2x)exp(-
b.随着时间的推移,a、b、c处的 温度分别自a、b、c时刻后 开始上升;
第四章_导热问题的数值方法
5 热传导问题的数值方法5.1一维稳态导热一维稳态导热在直角坐标系下的控制方程可表示为:0)(=+s dxdT k dx d (5-1) 式中k 为导热系数,T 是温度,s 是单位容积的热产生率。
首先选定控制体和网格,如图5.1所示,并对方程(5-1)在所选定的控制体进行积分,即得:0)()(=+-⎰dx s dxdTk dx dT k e w w e (5-2)图5.1 控制体和网格然后进行离散化。
如果用分线段性分布来计算方程(5-2)中的微商dxdT,那么最终的方程为:0)()()()(=∆+---x s x T T k x T T k wW P w e P E e δδ (5-3)假设源项s 在任一控制体中之值可以表示为温度的线性函数,即P P c T s s s +=,则导出的离散化方程为:b T a T a T a W W E E P P ++= (5-4)式中x s b xs a a a x k a x k a c P W E P w wW ee E ∆=∆-+=δ=δ=)()( (5-5) 式(5-4)就是一维稳态导热方程的离散形式,系数a E 和a W 分别代表了节点P 与E 间及W 与P 间导热阻力的倒数,它们的大小反映了节点W 和E 处的温度对P 点的影响程度。
式中的k e 和k w 是控制容积中的e 和w 界面上的当量导热系数。
进行计算时,物理参数值存储在节点的位置上。
为了确定k e 和k w ,还需规定由节点上的物理量来计算相应界面上的量的方法。
常用的方法由两种,即算术平均法与调和平均法。
1、算术平均法假定k 与x 呈线性关系,由P 与E 点的导数系数确定e k 的公式为:eeE e e P e x x k x x k k )()()()(δδ+δδ=-+ (5-6)2、调和平均法利用传热学的基本公式可以导出确定界面上当量导热系数的调和平均公式。
控制容积中P 和E 的导热系数不相等,但界面上热流密度应该连续,则由Fourier 定律可得:()()()()EePePE EeeE PePe e k x k x T T k x T T k x T T q +-+-δ+δ-=δ-=δ-=(5-7)而()Pe PE e k x T T q δ-=则()()()Ee Pe eek x k x k x +-+=δδδ (5-8)这就是确定界面上当量导热系数的调和平均公式,它反映了串联过程热阻的迭加原则。
第四章集总参数法
第四章 / 第三节 非稳态导热
(4)适用条件 Bi﹤0.1
Bi hl
厚度2 的无限大平壁
无限长圆柱体 球体
l
半厚 半径R 半径R
Bi
h / hR/ hR/
V/A
R/2 R/3
8
第四章 / 第三节 非稳态导热
几点说明
(1)以上分析结果既适用于物体被冷却的情况, 也适用于物体被加热的情况。
(2)以上计算公式针对第三类边界条件下导出。 在其他边界情况下的非稳态导热,只要物体内部
流换热条件hA。 5
第四章 / 第三节 非稳态导热
计算公式的应用:
ttf 0 t0tf
ex p(hcA V)
e xp BV (iFV o)
exp( c
)
(1)已知温度t,求时间; (2)已知时间,求温度t; (3)已知温度t和时间,求c或h。
6
第四章 / 第三节 非稳态导热
(3)换热量的计算
d t (2)以上计算公式针对第三类边界条件下导出。
h A(tt )cV f 1,应采用诺模图或其他方法重新计算。 d (3)利用集总参数法计算时,必须首先检验Bi﹤0.
求解:
❖ 物体温度分布t = f ( );
—使用的前ห้องสมุดไป่ตู้条件: Bi﹤0.
利用=集 4 总c参数稳法(态求第出h三或l;类边界条件)
后,温度变化dt。 这样热电偶越能迅速地反映被测流体的温度变化。
例如:小金属块炉内加热或空气中冷却;
为常数。
第四章 / 第三节 非稳态导热
热平衡方程: 第四章 / 第三节 非稳态导热 很大,或几何尺寸l很小,或h极低
假设满足Bi0.1的条件。
第四章 / 第三节 非稳态导热 —使用的前提条件: Bi﹤0.
传热学课件第四章非稳态导热
exp
hA
cV
hA
cV
h V
A
c
V
A2
hl
c
l2
hl
a
l2
BiV
FoV
0
e BiV FoV
exp
BiV FoV
下角标V表示以 l=V/A为特征长度
在0~ 时间内物体和周围环境之间交换的热量
升高到t1并保持不变,而右侧仍与温度为t0的 空气接触。这时紧挨高温表面那部分的温度
很快上升,而其余部分则仍保持初始温度t0, 如图中曲线HBD所示。随着时间的推移,经τ 1, τ 2,τ 3…平壁从左到右各部分的温度也依次 升高,从某一时刻开始平壁右侧表面温度逐
渐升高,图中曲线HCD、HE、HF示意性地表示
• 二、Bi数对导热体温度分布的影响
•
Bi hL L / 的大小对非稳态导热过程中导
热体内的 温1度/ h 分布有重要的影响。
• 厚为2δ的平壁突然置于流体中冷却时 ,Bi数 不同壁中温度场的变化会出现三种情形 。
思考题: 试说明毕渥数的物理意义。 毕渥数趋于
零和毕渥数趋于无穷各代表什么样的换热条件? 有人认为,毕渥数趋于零代表了绝热工况,你 是否赞同这一观点,为什么?
圆
球 Bi hR
Fo
a 2
BiV
h
FoV
a 2
Fo
a
R2
BiV
h(R / 2)
FoV
第四章 非稳态导热
f (0 , h, , , a, , x)
15
第四章 / 第三节 非稳态导热
(一)无限大平壁的分析解及诺模图
1、平壁内温度分布的求解
无穷级数
x , x 2 a 分析解: Cn exp n 2 cos n 0 n 1
2
2 a
分子—表示边界上发生热扰动时刻算起到计算时刻 为止的时间; 从过程开始到 时刻的时间 分母δ /温度变化波及到 a—表示热扰动经过一定厚度的固体层传播到 2面积所需的时间 面积δ2上所需要的时间。 Fo数看成是反映非稳态进程的无量纲时间。 Fo数越大,边界上的热扰动就能更深入地传播到 物体内部,非稳态过程进行得越充分。
物体内的温度分布受初始温度的影响很大,温度分布呈现部 壁面温差引起 分为非稳态导热规律和部分为初始温度区的混合分布。
B 正规状况阶段——整个物体参与变化
物体内的温度分布不再受初始温度的影响,而只受控于非 稳态导热的规律(边界条件、物性和几何因素的影响)。
热应力,会致 热变形!
(2)在非稳态导热热量传递的路径中,每一个与热流 方向垂直的截面上的热流量是处处不等的。
13
第四章 / 第三节 非稳态导热
(一)无限大平壁的分析解及诺模图
1、平壁内温度分布的求解
t 2t a 2 0 x , 0 x
初始条件: t | 0 t 0
0 x
边界条件: t | 0 (对称性) x 0
x
t |x h t |x t f x
(例如半无限大物体的导热)。
22
第四章 / 第三节 非稳态导热
第四章 非稳态导热(5)14
13
注
意
① 图4-4纵坐标为对数坐标,而图4-5和图4-6横坐标为对数坐标。三个图均为半 对数坐标系。 ② 图4-4中为直线关系,只在 Fo 0.2 时才是这样,即当过程进入正规状态阶段, 求解的无穷级数只取第一项( n 1)即满足精确要求。因此,成简单的指数函数 关系,它们在半对数坐数上为线性关系。否则,第二项以后的余项不能舍去,结果 就不是简单的指数函数关系,线图就不是图4-4的形式。 当 Fo 0.2 时会是什么样?例如:取数据: a 1.489105 m2 s, 100mm,求 得 134 s ,相对时间很短,一般工程上都不会加热或冷却这样短的时间,由图 4-4可见,数据集中在左上角很小的范围内,在整个图上占的份额很小。 如果确实需要计算 Fo 0.2 时,可用式(4-13)计算,即无穷级数解多取几级。 上述分析解的应用范围可以作三点推广:
2 1 2 2 1 2
m / e m 0 0 1 sin 1 cos 1
x 2 sin 1 cos(1 )
a
a
11
通过上述两个线算图分别查出 m ,
,利用 如果已知温度分布 t、x、Bi ,求τ,可以先计算出 ,再由图查出 0 m 1 m 公式 和 反查图得出 Fo 数,求出加热或冷却到此温度 m ,求出 o 0 m o Bi
(1)对无限大平板问题的分析是以平板被加热的情况为例的,上述结果对物体被冷 却的情况同样适用; (2)从无限大平板问题的数学描述式可以看出,分析解也适用于一侧绝热、另一侧 为第三类边界条件的厚为δ 的平板情形; (3) 当固体表面与流体间的表面传热系数趋于无穷大时,固体的表面温度就趋近于流 体温度,因而 Bi 时的上述分析解就是物体表面温度发生突然变化然后保 持不变时的解, 即第一类边界条件的解。
非稳态导热——精选推荐
非稳态导热一、基本概念本节基本概念主要包括:对物理问题进行分析,得出其数学描写(控制方程和定解条件);定性画出物体内的温度分布;集总参数法的定性分析;时间常数概念的运用;一维非稳态导热分析解的讨论;对海斯勒图(诺谟图)的理解;乘积解在多维非稳态导热中的应用;半无限大物体的基本概念。
1、由导热微分方程可知,非稳态导热只与热扩散率有关,而与导热系数无关。
你认为对吗?答:由于描述一个导热问题的完整数学描写不仅包括控制方程,还包括定解条件。
所以虽然非稳态导热的控制方程只与热扩散率有关,但边界条件中却有可能包括导热系数λ(如第二或第三类边界条件)。
因此上述观点不对。
2、无内热源,常物性二维导热物体在某一瞬时的温度分布为t=2y2cosx。
试说明该导热物体在x=0,y=1处的温度是随时间增加逐渐升高,还是逐渐降低。
答:由导热控制方程,得:当时,,故该点温度随时间增加而升高。
3、两块厚度为30mm的无限大平板,初始温度为20℃,分别用铜和钢制成。
平板两侧表面的温度突然上升到60℃,试计算使两板中心温度均上升到56℃时两板所需时间之比。
铜和钢的热扩散率分别为103×10-6m2/s,12.9×10-6m2/s。
答:一维非稳态无限大平板内的温度分布有如下函数形式:两块不同材料的无限大平板,均处于第一类边界条件(即Bi→∞)。
由题意,两种材料达到同样工况时,Bi 数和相同,要使温度分布相同,则只需Fo数相等,因此:,即,而δ在两种情况下相等,因此:4、东北地区春季,公路路面常出现“弹簧”,冒泥浆等“翻浆”病害。
试简要解释其原因。
为什么南方地区不出现此病害?东北地区的秋冬季节也不出现“翻浆”?答:此现象可以由半无限大物体(地面及地下)周期性非稳态导热现象的温度波衰减及温度波时间延迟特征来解释。
公路路面“弹簧”及“翻浆”病害产生的条件是:地面以下结冰,而地表面已解冻(表面水无法渗如地下)。
东北地区春季地表面温度已高于0℃,但由于温度波的时间延迟,地下仍低于0℃,从而产生了公路路面“弹簧”及“翻浆”等病害。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
2
材料
铸铁 砂型 金属型 46.5 0.314 61.64 753.6 963.0 544.3 7000 1350 7100
例题4-1 一大型平壁状铸铁件在砂型中凝固冷却。设砂 型内侧表面温度维持1200℃不变,砂型初始温度为 20℃,热扩散率������ = 2.41 × 10−7 m2 s,试求浇注后 1.5h砂型中离内侧表面50mm处的温度
4.3 伴有相变边界的一维非稳态导热
ⅆ������ ′ ������������ = −������������ ⅆ������ ������ ⟵ ������ + ������ ������������ − ������������
′ ������������
������ 砂 型 ������0 ⅆ������ ⅆ������
������
∞
1.温度场求解 常物性一维非稳态无内热源导热微分方程: ������������ ������ ������ 2 ������ ������ 2 ������ ������ 2 ������ ������ = + 2+ 2 + 2 ������������ ������������ ������������ ������������ ������������ ������������
������������ =0 边界条件: ������ = 0, ������������ ������������ ������ = ������, −������ = ℎ������ ������������
采用分离变量法求解:
������������ ������������
2 n
2sin n ( x, ) e 0 sin cos n 1 n n n
������ ������������
������
2
2������ = ������������ − ������0 ������ ������ = ������
������������������ = ������ J
m2 ⋅ ℃ ⋅ s 1
4.蓄热系数
热物性 导热系数 ������ W m℃ 比热容 ������ J m℃ 密度 ������ kg m3 热扩散率 ������ m2 s 8.82 × 10−6 2.41 × 10−7 1.58 × 10−5 蓄热系数 ������ J m2 ⋅ ℃ ⋅ s 1 15600 2030 15500
������
������������
������������
������3
������2 ������
������0
������1
������
积蓄(放出)热量随时间而变化是过程的又一特点。 瞬时热流密度+累计热量
第四章 非稳态导热
4.2 第一类边界条件下的一维非稳态导热
半无限大物体 一端面被一平面所限制,而另一端面延伸至无 +∞ 限远的物体。 ������
������������ ������ ������ 2 ������ ������ 2 ������ = = ������ 2 2 ������������ ������������ ������������ ������������
1.温度场求解 初始条件: ������ = 0 ⟶ ������ 边界条件: ������ > 0 ⟶ ������
第四章 界条件下一维非稳态导热问题
2.第三类边界条件下平板、圆柱、球体的导热 3.集总参数简化分析法 4.集中热源作用下的非稳态导热问题
4.1 非稳态导热基本概念
������ ������������
������3 ������2
������
������0
������1
������ erf 2 ������������
0.5
������ = 2 ������ ������ 2 = 2 ⇒ ������ = 16������ 2 ������������
0
0 1.0 ������ ������ = 2 ������������ 2.0
惰性时间
2.表面瞬时热流密度 ������������ ������������ = −������ ������������ ������ ������ = ������������ + ������0 − ������������ erf ������ ������������ ������0 − ������������ ������ 2 ������������ ������ ������ = ������0 − ������������ erf = exp − ������������ ������������ 2 ������������ 4������������ ������������������ ������������ ������������ = −������ ������������ = ������ ������������ − ������0
20℃ 砂型 1200℃
铸件
50
例题4-1
������ 50 × 10−3 m ������ = = 2 ������������ 2 2.41 × 10−7 m2 s × 1.5 × 3600s = 0.694
查表得:erf 0.694 = 0.6736
������ = ������������ + ������0 − ������������ erf ������ = 1200℃ + 20℃ − 1200℃ × 0.6736 = 405℃
铸件材质 灰口铸铁 可锻铸铁 铸钢
冷却条件
凝固系数 ������ × ������������������ / m s������
������
砂型
金属型 砂型
0.904
2.84 1.42
金属型
砂型 铸铁硬模
2.58
1.68 3.35
砂型
黄铜 铸铁硬模 水冷铜模
2.32
4.90 5.42
铝
铸铁硬模
4.0
������=0
1 ������������������
W m2
3.累计热量 ������������ = ������ ������������ − ������0
������
1 ������ = ������������ − ������0 ������������������ ������������
例题4-2
一个 140mm厚的平板铸钢件在砂型中冷却凝固, 其两侧释放热量给砂型。已知浇注温度������������ = 1540℃, 凝固潜热 ������ = 268 kJ kg ,钢的比热容 ������ =例题 873 J4-2 kg⋅℃ , 密度������ = 7200 kg m3 。砂型初始温度������0 = 10℃ ,蓄热 系数������ = 2030 J m2 ⋅ ℃ ⋅ s 1 2 。与铸件接触的砂型表 面温度近似等于钢的熔点������������ ≈ 1510℃ ,求全部凝固所 需时间。
例题4-2 ������
砂 型
������������ 液 凝 凝 态 固 固 金 层 属 层 ������ ������������ 砂 型
������0
������ =
2������ ������������ − ������0
������������ ������ + ������ ������������ − ������������
弗兰克著《传热和传质基本 原理》第六版,葛新石翻译
������=∈∀
= ������0 = ������1
������=0
= ������������ = ������2
erfx 2
������
������������ − ������ ������ = erf =erf ������ ������������ − ������0 2 ������������
ⅆ������ ⅆ������
ⅆ������ ������ ������������ − ������0 = ⅆ������ ������������������ ������ + ������ ������������ − ������������ 2������ ������������ − ������0
������������
������
������
������
������
������ = ������ − ������������ 过余温度
������������ ������ 2 ������ = ������ 2 导热微分方程: ������������ ������������
初始条件: ������ = 0, ������ = ������0 − ������������ = ������0
1 ������������
������ 0
������������ =
0
������������ ⅆ������ = ������ ������������ − ������0
= 2������ ������������ − ������0
ⅆ������ ������ J m2
������������ ∝ ������
初始条件: ������ = 0, ������ = ������0
������
������0
������������ 边界条件: ������ = 0, =0 ������������ ������������ ������ = ������, −������ = ℎ ������ − ������������ ������������