非稳态传热计算方法及举例
- 格式:docx
- 大小:88.73 KB
- 文档页数:8
传热学--第三章第三节一维非稳态导热问题§3 — 3 一维非稳态导热的分析解本节介绍第三类边界条件下:无限大平板、无限长圆柱、球的分析解及应用。
如何理解无限大物体,如:当一块平板的长度、宽度>> 厚度时,平板的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少,以至于可以把平板内各点的温度看作仅是厚度的函数时,该平板就是一块“无限大”平板。
若平板的长度、宽度、厚度相差较小,但平板四周绝热良好,则热量交换仅发生在平板两侧面,从传热的角度分析,可简化成一维导热问题。
一、无限大平板的分析解已知:厚度的无限大平板,初温t0,初始瞬间将其放于温度为的流体中,而且>t0,流体与板面间的表面传热系数为一常数。
试确定在非稳态过程中板内的温度分布。
解:如图3-5 所示,平板两面对称受热,所以其内温度分布以其中心截面为对称面。
对于x 0 的半块平板,其导热微分方程:(0<x< , )定解条件:t(x,0)= t0(0 x )(边界条件)(边界条件)引入过余温度:则(0<x< , )(3-9)(x,0)= (0 x ) (初始条件)(边界条件)(边界条件)对偏微分方程分离变量求解得:(3-10 )其中离散值是下列超越方程的根,称为特征值。
其中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 ,τ) 与(τ )各自均与τ 有关,但其比值则与τ 无关,而仅取决于几何位置(x/ )及边界条件(Bi )。
第三章 非稳态导热分析解法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 分布的影响较大,此阶段称非正规状况阶段。
1 第三章 非稳态导热的分析计算 3-1 非稳态导热过程分析 一、非稳态导热过程及其特点 导热系统(物体)内温度场随时间变化的导热过程为非稳态导热过程。在过程的进行中系统内各处的温度是随时间变化的,热流量也是变化的。这反映了传热过程中系统内的能量随时间的改变。我们研究非稳态导热过程的意义在于,工程上和自然界存在着大量的非稳态导热过程,如房屋墙壁内的温度变化、炉墙在加热(冷却)过程中的温度变化、物体在炉内的加热或在环境中冷却等。归纳起来,非稳态导热过程可分为两大类型,其一是周期性的非稳态导热过程,其二是非周期性的非稳态导热过程,通常指物体(或系统)的加热或冷却过程。这里主要介绍非周期性的非稳态导热过程。下面以一维非稳态导热为例来分析其过程的主要特征。 今有一无限大平板,突然放入加热炉中加热,平板受炉内烟气环境的加热作用,其温度就会从平板表面向平板中心随时间逐渐升高,其内能也逐渐增加,同时伴随着热流向平板中心的传递。图3-1显示了大平板加热过程的温度变化的情况。 从图中可见,当0时平板处于均匀的温度0tt下,随着时间的增加平板温度开始变化,并向板中心发展,而后中心温度也逐步升高。当时平板温度将与环境温度拉平,非稳态导热过程结束。图中温度分布曲线是用相同的来描绘的。总之,在非稳态导热过程中物体内的温度和热流都是在不断的变化,而且都是一个不断地从非稳态到稳态的导热过程,也是一个能量从不平衡到平衡的过程。 二、加热或冷却过程的两个重要阶段
从图3-1中也可以看出,在平板加热过程的初期,初始温度分布0tt仍然在影响物体整个的温度分布。只有物体中心的温度开始变化之后(如图中τ>τ2之后),初始温度分布0tt的影响才会消失,其后的温度分布就是一条光滑连续的曲线。据此,我们可以把非稳态导热过程分为两个不同的阶段,即: 初始状况阶段――环境的热影响不断向物体内部扩展的过程,也就是物体(或系统)仍然有部分区域受初始温度分布控制的阶段; 正规状况阶段――环境对物体的热影响已经扩展到整个物体内部,且仍然继续作用于物体的过程,也就是物体(或系统)的温度分布不再受初始温度分布影响的阶段。 由于初始状况阶段存在初始温度分布的影响而使物体内的整体温度分布必须用无穷级数来加以描述,而在正规状况阶段,由于初始温度影响的消失,温度分布曲线变为光滑连续的曲线,因而可以用初等函数加以描述,此时只要无穷级数的首项来表示物体内的温度分布。 3 边界条件对导热系统温度分布的影响 从上面的分析不难看出,环境(边界条件)对系统温度分布的影响是很显著的,且在整个过程中都一直在起作用。因此,分析一下非稳态导热过程的边界条件是十分重要的,
题目
一厚度为0.1m的无限大平壁,两侧均为对流换热边界条件,初始时 两侧流体温
度与壁内温度一致,tf1=tf2=t0=5 ℃;已知两侧对流换热系数分别为h1=11
W/m2K、h2=23W/m2K, 壁材料的导热系数=0.43W/mK,导温系数a=0.3437×
10-6 m2/s。如果一侧的环境温度tf1突然升高为50℃并维持不变,计算在其它参
数不变的条件下,平壁内温度分布及两侧壁面热流密度随时间的变化规律(用图
形表示)。
问题分析
此题为两侧受恒温流体作用,并求其从非稳态传热过程温度场到接近稳态传热的
温度场,并算出其热流密度随时间的变化规律。
解法
建立离散方程及求解
将平板分割成如下网格:共计10个网格,11个节点,以恒温流体1处为节点1,
恒温流体2处为节点11。
列写节点方程,边界条件皆为恒温流体传热,初始条件为5摄氏度。以此对每个
单独时刻进行求解,解出该时刻各节点的温度,并在此解的基础上进一步解出之
后各时段的温度解,进行迭代计算,直到满足时间要求为止。
非稳态传热计数器
计算过程使用Excel实现,具体做法是利用Excel进行解方程,并求出温度解。
因使用10个网格,故方程类型为10元1次方程组,也就是说每个时刻都有10
个方程必须联立求解,使用Excel的行列式计算能很容易地用克拉姆方法解出该
方程。之后用该组温度解进行下一次迭代运算,如此反复,直到满足题设要求。
具体的温度求解请查阅 非稳态传热计算器.xlsx 文件,为了要求计算器的整洁
美观,繁琐的计算过程使用Hide功能隐藏,若需查阅解除Hide指令即可。
使用计算器时仅需输入相关系数,并输入合适的时间步长即可,计算器将按给定
的参数计算出平板在之后各个时刻各节点上的温度值。
计算器将列出各节点的温度值随时间变化的计算表格,同时输出三种图形:平板
内各节点温度随时间变化规律,平板内各节点温度在某一时刻的变化规律及平板
壁面热流密度随时间变化规律。
网格1 网格2 网格3 网格4 网格5 网格6 网格10 网格9 网格8 网格7 恒温流体
Tf2=5摄氏度
H2=23 恒温流体 Tf1=50摄氏度 H1=11
计算器使用实例(按作业题目要求)
两侧受恒温流作用无限大平板非稳态传热
计算器(沿轴向分为10个网格)
请输入相关参数
平板厚度 0.10
平板热扩散率 0.0000003437
平板导热系数 0.43
传热系数1 11.00
传热系数2 23.00
平板初温 5.00
流体1温度 50.00
流体2温度 5.00
时间步长 30.00
校核结果 校核结果正确
网络傅里叶数 0.412440
网络毕渥数1 0.127907
网络毕渥数2 0.267442
平板内各节点温度随时间变化规律
节点 0秒 30秒 60秒 90秒 120秒 150秒
1
5.00 7.74 9.51 10.78 11.77 12.58
2
5.00 5.65 6.48 7.31 8.08 8.79
3
5.00 5.16 5.45 5.83 6.26 6.71
4
5.00 5.04 5.13 5.28 5.48 5.71
5
5.00 5.01 5.04 5.09 5.17 5.28
6
5.00 5.00 5.01 5.03 5.06 5.10
7
5.00 5.00 5.00 5.01 5.02 5.04
8
5.00 5.00 5.00 5.00 5.01 5.01
9
5.00 5.00 5.00 5.00 5.00 5.00
10
5.00 5.00 5.00 5.00 5.00 5.00
11
5.00 5.00 5.00 5.00 5.00 5.00
平板壁面热流密度随时间变化规律
节点 0秒 30秒 60秒 90秒 120秒 150秒
1 0.50 0.46 0.45 0.43 0.42 0.41
11 0.00 0.00 0.00 0.00 0.00 0.00
0.00
5.00
10.00
15.00
20.00
25.00
30.00
35.00
0100020003000400050006000
温
度
(
摄
氏
度
)
时间(秒)
平板内各节点温度随时间变化规律
1234567891011
0.00
5.00
10.00
15.00
20.00
25.00
30.00
35.00
024681012
温
度
(
摄
氏
度
)
节点
平板内各节点温度在某一时刻的变化规律
30秒
210秒
600秒
1620秒
4440秒
5310秒
热流密度方向为向右流动为正值,向左流动为负值。
关于计算程式的具体做法请直接打开计算器文件。
0.00
0.10
0.20
0.30
0.40
0.50
0.60
0100020003000400050006000
热
流
密
度
(
千
瓦
每
平
方
米
)
时间(秒)
平板壁面热流密度随时间变化规律
壁面1
壁面2
总结与分析
随着时间增加,通过壁面1及壁面2的流体密度趋近于相等,最终结果为热流从高流体
那一侧流向低温流体。
壁面与流体温差越大,热流密度也就越大,也就是说随着热流充分发展,热流密度的变
化也趋近于平滑。
随着时间的增加,在某一时刻平板内部各点的温度值逐渐趋近于一条直线,最终成为稳
态导热形式,在此之前,温度值的变化规律符合热传导微分方程,并按照指数exp(-t)
的形式规律来变化。