第三章非稳导热分析解法文档格式.docx
- 文档编号:21504956
- 上传时间:2023-01-30
- 格式:DOCX
- 页数:26
- 大小:181.74KB
第三章非稳导热分析解法文档格式.docx
《第三章非稳导热分析解法文档格式.docx》由会员分享,可在线阅读,更多相关《第三章非稳导热分析解法文档格式.docx(26页珍藏版)》请在冰豆网上搜索。
前者存在着有区别的两个不同阶段,而后者不存在。
3、特点;
非稳态导热过程中,在与热流量方向相垂直的不同截面上热流量不相等,这是非稳态导热区别于稳态导热的一个特点。
原因:
由于在热量传递的路径上,物体各处温度的变化要积聚或消耗能量,所以,在热流量传递的方向上
。
二、非稳态导热的数学模型
1、数学模型
非稳态导热问题的求解
规定的{初始条件,边界条件}下,求解导热微分方程。
2、讨论物体处于恒温介质中的第三类边界条件问题
在第三类边界条件下,确定非稳态导热物体中的温度变化特征与边界条件参数的关系。
已知:
平板厚2
、初温to、表面传热系数h、平板导热系数
,将其突然置于温度为
的流体中冷却。
试分析在以下三种情况:
<
1/h、
>
=1/h时,平板中温度场的变化。
1)1/h<
因为1/h可忽略,当平板突然被冷却时,其表面温度就被冷却到
,随着时间的延长,平板内各点
t→
如图3-3(a)。
2)1/h>
因为
忽略不计,即平板内导热的流量接近于无穷大,所以任意时刻平板中各点温度接近均匀,随着时间的延长,平板内各点t→
,而且整体温度下降如图3-3(b)。
3)1/h=
平板中的温度分布介于二者之间,如图3-3(c)。
由此可见,表面对流换热热阻1/h与导热热阻
的相对大小对物体中非稳态导热的温度场的分布有重要影响,因此,引入表征二者比值的无量纲数,毕渥数。
3、毕渥数
1)定义式:
(3-1)
毕渥数属特征数(准则数)。
2)Bi物理意义:
Bi的大小反映了物体在非稳态条件下内部温度场的分布规律。
3)特征数(准则数):
表征某一物理现象或过程特征的无量纲数。
4)特征长度:
是指特征数定义式中的几何尺度。
3—2集总参数法的简化分析
一、集总参数法
当固体内的
时,固体内的温度趋于一致,此时可认为整个固体在同一瞬间均处于同一温度下,这时需求解的温度仅是时间的一元函数,而与坐标无关,好象该固体原来连续分布的质量与热容量汇总到一点上,而只有一个温度值那样。
这种忽略物体内部导热热阻的简化分析方法称为集总参数法。
2、集总参数法的计算
有一任意形状的物体,其体积为V,面积为A,初始温度为t0,在初始时刻,突然将其置于温度恒为
的流体中,且to>
固体与流体间的表面传热系数h,固体的物性参数均保持常数。
试根据集总参数法确定物体温度随时间的依变关系。
解:
①建立非稳态导热数学模型
方法一:
椐非稳态有内热源的导热微分方程:
∵物体内部导热热阻很小,忽略不计。
∴物体温度在同一瞬间各点温度基本相等,即t仅是τ的一元函数,二与坐标x、y、z无关,即
=0
则:
(a)
∵
可视为广义热源,而且热交换的边界不是计算边界(零维无任何边界)。
∴界面上交换的热量应折算成整个物体的体积热源,即:
(b)
∵t>
物体被冷却,∴
应为负值
由(a),(b)式得:
(3-2)
这就是瞬时时刻导热微分方程式。
方法二:
根据能量守恒原理,建立物体的热平衡方程,即
物体与环境的对流散热量=物体内能的减少量
则有:
②物体温度随时间的依变关系
引入过余温度:
则上式表示成:
其初始条件为:
将
分离变量求解微分方程,
对时间
从0
积分,则:
In
即:
(
)(3-3)
其中:
V/A是具有长度的量纲,记为
;
毕渥数;
傅立叶数;
而V说明Fov、Biv中的特征长度为V/A
故得:
(3-4)
由此可见,采用集总参数法分析时,物体内的过余温度随时间成指数曲线关系变化。
而且开始变化较快,随后逐渐变慢。
指数函数中的
的量纲与
的量纲相同,如果时间
,则
称时间常数,记为
的物理意义:
表示物体对外界温度变化的响应程度。
当时间
时,物体的过余温度已是初始过余温度值的36.8%。
③确定从初始时刻到某一瞬间这段时间内,物体与流体所交换的热流量
首先求得瞬时热流量:
带入瞬时热流量的定义式得:
=
(3-5)
式中负号是为了使Φ恒取正值而引入的。
若
(物体被加热),则用
代替
即可。
然后求得从时间
0到
时刻间的总热流量:
(3—6)
3、集总参数法的判别条件
对形如平板、圆柱和球这一类的物体,如果毕渥数满足以下条件:
=h(V/A)/
0.1M(3-7)
则物体中各点间过余温度的偏差小于5%。
其中M是与物体几何形状有关的无量纲数。
无限大平板:
M=1
无限长圆柱:
M=1/2
球:
M=1/3
毕渥数的特征长度为V/A,不同几何形状,其值不同,对于:
厚度为2
的平板:
半径为R的圆柱:
半径为R的球:
由此可见,对平板:
=Bi
圆柱:
=Bi/2
球体:
=Bi/3
二、毕渥数
与傅立叶数
的物理意义
1、
1)定义:
表征固体内部单位导热面积上的导热热阻与单位面积上的换热热阻(即外部热阻)之比。
越小,表示内热阻越小,外部热阻越大。
此时采用集总参数法求解更为合适。
2)物理意义:
的大小反映了物体在非稳态导热条件下,物体内温度场的分布规律。
2、
表征两个时间间隔相比所得的无量纲时间。
分子τ是从边界上开始发生热扰动的时刻起到所计时刻为止的时间间隔。
分母可视为边界上发生的有限大小的热扰动穿过一定厚度的固体层扩散到
的面积上所需的时间。
表示非稳态导热过程进行的程度,
越大,热扰动就越深入地传播到物体内部,因而物体内各点的温度越接近周围介质的温度。
3—3一维非稳态导热的分析解
本节介绍第三类边界条件下:
无限大平板、无限长圆柱、球的分析解及应
用。
如何理解无限大物体,如:
当一块平板的长度、宽度>
厚度时,平板
的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少,以至于可
以把平板内各点的温度看作仅是厚度的函数时,该平板就是一块“无限大”
平板。
若平板的长度、宽度、厚度相差较小,但平板四周绝热良好,则热量
交换仅发生在平板两侧面,从传热的角度分析,可简化成一维导热问题。
一、无限大平板的分析解
厚度
的无限大平板,初温
,初始瞬间将其放于温度为
的
流体中,而且
,流体与板面间的表面传热系数为一常数。
试确定在非稳态过程中板内的温度分布。
如图3-5所示,平板两面对称受热,所以其内温度分布以其中心截面为对称面。
对于x
0的半块平板,其导热微分方程:
(0<
x<
)(3-8)
定解条件:
t(x,0)=
(0
x
)
(边界条件)
则
(0<
)(3-9)
(x,0)=
)(初始条件)
(边界条件)
对偏微分方程
分离变量求解得:
(3-10)
其中离散值
是下列超越方程的根,称为特征值。
……(3-11)
其中Bi是以特征长度为
的毕渥数。
由此可见:
平板中的无量纲过余温度
与三个无量纲数有关:
以平板厚度一半
为特征长度的傅立叶数、毕渥数及
(3-12)
二、非稳态导热的正规状况阶段
1、平板中任一点的过余温度与平板中心的过余温度的关系
前述得到的分析解是一个无穷级数,计算工作量大,但对比计算表明,当Fo>
0.2时,采用该级数的第一项与采用完整的级数计算平板中心温度的误差小于1%,因此,当Fo>
0.2时,采用以下简化结果:
(3-13)
其中特征值
之值与Bi有关。
由上式(3-13)可知:
Fo>
0.2以后平板中任一点的过余温度
(x,τ)与平板中心的过余温度
(0,τ)=
(τ)之比为:
(3-14)
此式反映了非稳态导热过程中一种很重要的物理现象:
即当Fo>
0.2以后,虽然
(x,τ)与
(τ)各自均与τ有关,但其比值则与τ无关,而仅取决于几何位置(
)及边界条件(Bi)。
也就是说,初始条件的影响已经消失,无论初始条件分布如何,只要Fo>
0.2,
之值是一个常数,也就是无量纲的温度分布是一样的。
由此可见,当Fo>
0.2时,非稳态导热过程进入正规状况阶段。
2、在一个时间间隔内非稳态导热过程中传递的热量
1)从物体初始时刻平板与周围介质处于热平衡,这一过程中传递的热量:
(3-15)
此值为非稳态导热过程中传递的最大热量。
2)从初始时刻到某一时间τ,这段时间内所传递的热量
:
(3-16)
3)
之比:
(3-17)
是时刻τ物体的平均过余温度,
对于无限大平板,当Fo>
0.2,将式(3-13)代入
的定义式,可得:
(3-18)
对圆柱体、球体
0.2时,无穷级数的解也可用第一项近似代替,并且
及
可表示为:
(3-19)
(3-30)
η为无量纲几何位置,对平板
,对柱体及球体
,R为外表面半径,系数A、B及函数
的表达式取决于几何形状,见教材表3-2所示。
三、正规阶段状况的实用计算方法
当Fo>
0.2时,可采用上述计算公式求得非稳态导热物体的温度场及交换的热量,也可采用简化的拟合公式和诺模图求得。
1、诺模图:
工程技术中,为便于计算,采用按分析解的级数第一项绘制的一些图线,叫诺模图。
2、海斯勒图:
诺模图中用以确定温度分布的图线,称海斯勒图。
首先根据(3—13)式给出
随Fo及Bi变化的曲线(此时x/δ=0),然后根据(3—14)式确定
的值,于是平板中任意一点的
值便为:
(3-21)
同样,从初始时刻到时刻τ物体与环境间所交换的热量,可采用(3—15)、(3—17)作出
曲线。
3、诺模图法评述
优点:
简洁方便。
缺点:
准确度有限,误差较大。
目前,随着计算技术的发展,直接应用分析解及简化拟合公式计算的方法受到重视。
四、分析解应用范围的推广及讨论
1、推广范围
1)对物体被冷却的情况也适用;
2)也适于一侧绝热,另一侧为第三类边界条件的厚为δ的平板;
3)当固体表面与流体间的表面传热系数h
时,即表面换热热阻
0时,所以
时分析解就是固体表面温度发生一突然变化然后保持不变时的解,即第一类边界条件的解。
2、讨论Bi与Fo对温度场的影响:
1)傅立叶数Fo:
由(3-10)、(3-13)式及诺模图可知:
物体中各点的过余温度随时间τ的增加而减小;
而Fo与成正比,所以物体中各点过余温度亦随Fo的增大而减小。
2)毕渥数Bi
Bi对温度的影响从以下两方面分析:
一方面,从教材图3—6可知,Fo相同时,Bi越大,
越小。
因为,Bi越大,意味着固体表面的换热条件越强,导致物体的中心温度越迅速地接近周围介质的温度;
当Bi
时,意味着在过程开始瞬间物体表面温度就达到介质温度,物体中心温度变化最快,所以在诺模图中1/Bi=0时的线就是壁面温度保持恒定的第一类边界条件的解。
另一方面Bi的大小决定于物体内部温度的扯平程度。
对于平板,从诺模图3—7中可知:
当
10(即Bi<
0.1)时,截面上的过余温度差小于5%
当Bi下限一直推到0.01时,其分析解与集总参数法的解相差极微。
综上可得如下结论:
介质温度恒定的第三类边界条件下的分析解;
时,转化为第一类边界条件下的解,Bi
0时,则与集总参数法的解相同。
3—4二维及三维非稳态导热问题的求解
一、求解方法
对于典型的几何形状的物体,可利用一维非稳态导热问题分析解的组合求得。
如图3-9所示:
无限长方柱体的非稳态导热问题,属二维导热问题。
截面
尺寸为:
的方柱体可视为两块厚度分别为
的无限大平板
垂直相交所截出的物体。
讨论的目的:
找出二维温度场与两块无限大平板的温度场之间的关系。
方柱体初温为to,初始时放于t
流体中,表面传热系数为h.
试求:
温度场分布。
如图3-9所示,建立坐标系,由于其对称性,只研究其
截面的温度分布,截面上的温度分布由下列导热微分方程和定解条件确定:
(3-22)
(I)
(II)
(III)
(IV)
式中:
为无量纲过余温度。
如果无量纲过余温度
与
分别是处于与方柱体同样定解条件下的厚度分别为
的无限大平板的分析解,则它们必须满足各自的导热微分方程及定解条件,即:
(3-23)
(V)
(VI)
(VII)
及
(3-24)
只要证明:
两块无限大平板分析解的乘积就是上述无限长方柱体的分析解,即:
(3-25)
证明:
首先证明式(3-25)满足导热微分方程(3-22),为此将式(3-25)代入式(3-22)的左右两端得:
左端:
右端:
左端减去右端得:
-
+
∴证明
满足微分方程
其次证明:
满足初始条件
根据
的初始条件
和
得:
最后证明:
满足边界条件
将式(3-25)代入边界条件(I),并注意到式(VII)的关系得
[
]=
*0=0
同样可以证明它也满足式(3-34)。
再将式(3-25)代入边界条件(III),并注意到式(VI)的关系得
=0*
同理可证明它也满足式(IV)。
综上可知:
是上述无限长方柱体导热微分方程的解。
结论:
此方法是多维非稳态导热的求乘积解法,此法适用于第一类边界条件,且
时。
同理,对长圆柱体,矩六柱体等二维,三维非稳态导热问题,可以用相应的二个或三个一维问题的解的乘积来表示其温度分布。
——表示无限大平板的解
——表示无限长圆柱体的解
(a):
(b):
(c):
二、乘积解法的适用条件
•
初始温度为常数,
const;
第一类边界条件,
第三类边界条件,
const、h=const;
线性微分方程,且定解条件均为齐次,即乘积解中温度必须以过余温度或无量纲过余温度的形式表示。
说明:
对于形状复杂或边界条件复杂,分析解法无能为力,应借助其它的求解的方法,
如①数值解法;
②实验模拟法。
3—5半无限大物体的非稳态导热
一、半无限大物体的概念
几何上是指从x=0的界面开始可以向正的x方向及其他两个坐标(x,y)方向无限延伸的物体,称半无限大物体。
实际中不存在该物体,但研究物体中非稳态导热的初始阶段,可把实物看为该物体处理。
有限厚度的平板,起初有均匀温度,后其侧表面突然受到热扰动,如
壁温突然升高到一定值并保持不变;
壁面突然受到恒定的热流量密度加热;
壁面受到温度恒定的流体的加热或冷却。
当扰动的影响只局限在表面附近,而尚未进入平板内部时,就可视该平板为,“半无限大”物体。
二、第一类边界条件下半无限大物体非稳态导热温度场的分析解
如图3-11所示:
已知半无限大物体初始温度均匀为(to),当
=0
时,x=0侧表面温度突然升高到
,并保持不变,试确定物体内温度随时间
的变化和在时间间隔[0,τ]内的热流量。
1、物体内的温度分布
根据半无限大物体的定义,得出其导热微分方程:
(3-26)
初始条件为:
τ=0时,
边界条件为:
x=0时,
x→
时,
(3-27)
τ=0时,
x=0时,
将微分方程
分离变量并求解得分析解为:
=erf(
)=erfη(3-28)
无量纲变量η=
erfη称为误差函数,它随η的变化而变化,由附录表可知:
当η=2时,
=0.9953,就是说当η
2即
2时,该处x的温度仍认为等于to(无量纲过余温度的变化小于5%),由此得到以下两个重要参数:
①从几何位置上说,若
,则时刻τ时x处的温度可认为未发生变化。
所以,对
且厚为2δ的平板,当其一侧温度突然变化到另一恒定温度时,若δ
,则在τ时刻之前该平板中瞬时温度场的计算可采用半无限大物体模型处理。
②从时间上看,如果
,则此时x处的温度可认为完全不变,所以把
视为惰性时间,即当
时x处的温度可认为仍等于
2、表面上的瞬时热流密度及在[0,τ]时间间隔内放出或吸收的热量:
物体中任意一点的热流密度:
(3-29)
则,表面上的热流密度为:
在时间[0,τ]内,流过面积A的总热流流量:
(3-31)
①半无限大物体在第一类边界条件影响下被加热或冷却时,界面上的瞬时热流量与时间的平方根成反比;
②在时间[0,
]内交换的总热量则正比于
及时间的平方根。
称为吸热系数,表示物体与其接触的高温物体吸热的能力。
三、半无限大物体概念的适用范围
只适于物体非稳态导热的初始阶段,当物体表面上的热扰动已深入传递到物体内部时,就不再适用,则应采用前述分析方法。
补充:
一、非稳态导热问题求解思路
解一维
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 第三 章非稳 导热 分析 解法