当前位置：文档之家 › 第四章 非稳态导热的分析计算

第四章 非稳态导热的分析计算

  • 格式：ppt
  • 大小：130.00 KB
  • 文档页数：14

下载文档原格式

  / 14

合集下载

传热学-第4章-非稳态导热的计算与分析

传热学-第4章-非稳态导热的计算与分析
3
‹# ›
第4章 非稳态导热的计算与分析
本章着重讨论非稳态导热问题 ——非稳态导热的基本概念 ——对称加热的无限大平壁的非稳态导热过程 ——最简单的非稳态导热问题－集总热容系统
4
‹# ›
第4章 非稳态导热的计算与分析
4.1 概述
非稳态导热的分类： ——周期性的非稳态导热（periodic unsteady heat conduction）：由于边界条件（或内热源）随时间呈周 期性变化，使物体内的温度场也随时间按周期性规律变 化，这种状况通常称为准稳态
19
‹# ›
4.2.1 平壁内非稳态过程的基本特征
整个瞬态导热过程可以分为两个阶段： 初始阶段（initial regime）：也称为非正规状况阶段，
指在穿透时刻之前阶段，此时平壁内的温度分布主要受 初始温度分布t0的影响。
正规状况阶段（regular regime）：穿透时刻之后，非稳态 过程进行到一定的程度，平壁初始温度分布的影响逐渐消失，此 后不同时刻的温度分布主要受热边界条件的影响。这个阶段的非 稳态导热称为正规状况阶段。
第4章 非稳本态节导内热容的结计束 算与分析
1
• 稳态导热是一种理想化的情况 • 受环境温度变化的影响，生活和工程中真正意义上的稳 态导热是不存在的 • 只是对工程中的某些问题，忽略温度随时间变化所造成 的影响、误差不大，而将其简化为稳态导热
2
‹# ›
• 生活和工程中还存在着大量的不能简化为稳态导热的现 象和问题，其中物体内的温度明显随时间而变化 • ——冷冻食品的解冻过程 • ——烘烤食品（花生米、蛋糕等点心） • ——热处理工艺中金属在高温火炉内的加热以及加热后 在水或空气中的冷却过程等 • ——焖井过程热量在地层内的扩散过程

第4章 导热问题的数值解法(含控制容积)

第4章 导热问题的数值解法(含控制容积)
（1）建立符合实际的物理模型 对实际导热问题的几何、物理性质进行分析，做必要的、 对实际导热问题的几何、物理性质进行分析，做必要的、合 理的简化，建立符合实际的物理模型； 理的简化，建立符合实际的物理模型； （2）建立控制方程及定解条件 根据物理模型建立完整的数学模型， 根据物理模型建立完整的数学模型，即给出导热微分方程和 单值性条件； 单值性条件； 步是导热问题所有求解方法的基础。 第(1)、(2)步是导热问题所有求解方法的基础。 2012-5-9 4
ti−1, j
二方程相加， 二方程相加，得：
ti+1, j − 2ti, j +ti−1, j ∂2t 2 = + 0(∆x2 ) ∂x ∆x2 i, j ti, j+1 − 2ti, j +ti, j−1 ∂2t 2 = + 0(∆y2 ) ∂y ∆y2 i, j
ti, j −ti−1, j ∂2t ∆x ∂3t ∆x2 ∂t + 2 − 3 +...... = ∂x 2 ∂x ∆x ! ∂x i, j i, j ! i, j 3 = ti, j −ti−1, j ∆x + 0(∆x)
2012-5-9 2
§4-1 导热问题数值求解的基本思想 及内节点离散方程的建立
2012-5-9
3
一、数值解法的基本思想 用导热问题所涉及的空间和时间区域内有限个离散 称为节点 节点) 点(称为节点)的温度近似值来代替物体内实际连续的温 度分布， 度分布 ， 将连续温度分布函数的求解问题转化为各节 点温度值的求解问题， 点温度值的求解问题 ， 将导热微分方程的求解问题转 化为节点温度代数方程的求解问题。 化为节点温度代数方程的求解问题。 数值解法的基本内容与步骤: 数值解法的基本内容与步骤:

传热学第四章

传热学第四章

第四章 非稳态导热
第一节 概 述
a)温度分布；b)两侧表面上导热量随时间的变化
图4-1
第四章 非稳态导热
第一节 概 述
（1）温度场：【如图4-1a）所示】 ①首先，紧挨高温表面部分的温度很快上升， 而其余部分仍保持原来的温度t0，如图中曲线FBC所示； ②其次，随着时间的推移，温度变化波及的范围不断扩大， 以致在一定时间以后，右侧表面的温度也逐渐升高， 如图中曲线FC、FD所示； ③最后，达到一个新的稳态导热时，温度分布保持恒定， 如图中曲线FE所示。（λ为常数时，FE 为直线。）
t f ( x, y, z, )
dt （3）物体在非稳态导热过程中的温升速率： d
（4）某一时刻物体表面的热流量Φ(W) 或从某一时刻起经过一定时间后表面传递的总热量Q(J)。 要解决以上问题，必须首先求出： 物体在非稳态导热过程中的温度场。
第四章 非稳态导热
第一节 概 述
※求解非稳态导热过程中物体的温度场，通常可采用
第四章 非稳态导热
第一节 概 述 一、基本概念
非稳态导热即指温度场随时间而变化的导热过程 1、定义（P53）
t f ( x, y, z, )
※在自然界和工程中有许多非稳态导热问题。 例如，锅炉、蒸汽轮机和内燃机等动力机械在起动、停机和变 工况运行时的导热； 又如，在冶金、热处理和热加工等过程中，工件被加热或冷却 时的导热； 再有，大地和房屋等白天被太阳加热、夜晚被冷却时的导热。 ※由此可见，研究非稳态导热具有很大的实际意义。
l
—— 导热物体的某一尺寸，详见后述。
第四章 非稳态导热
第一节 概 述
1、毕渥数Bi （P55）
有时用引用尺寸l
e
l ——导热物体的某一尺寸

传热学-第4章-非稳态导热的计算与分析

传热学-第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)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

非稳态导热分析解法课件

复杂边界条件和几何形状
非稳态导热问题常常涉及到复杂的边界条件和几何形状，给分析带来很大挑战。未来发展需要研究更高效的数值方法 ，以处理更复杂的导热问题。
多物理场耦合
许多实际导热问题涉及到多物理场的耦合，如热-力、热-流体等。未来发展需要研究多物理场耦合的非稳态导热问题 ，以提高对复杂系统的理解和预测能力。
高效能材料和新能源技术
随着高效能材料和新能源技术的发展，非稳态导热问题将更加复杂和多样化。未来发展需要加强与相关 领域的交叉融合，以应对不断出现的新的挑战和机遇。
核能利用
在核能利用中，非稳态导热分析可用于研究反应堆的冷却系统、核废料的处理和存储等。 通过优化导热性能，可以提高核能系统的安全性和稳定性。
风能利用
在风能利用中，非稳态导热分析可用于研究风力发电机的散热性能和风能转换效率。通过 改进导热设计，可以提高风能发电的经济性和可靠性。
非稳态导热面临的挑战和未来发展方向
物理模拟实验
物理模拟实验是通过模拟实际系统的物理过程来研究其行为的方法。
在非稳态导热分析中，物理模拟实验通常采用加热棒、散热片等模拟导热过程，通 过测量温度场、热流密度等参数来研究非稳态导热规律。
物理模拟实验具有直观、可重复性高等优点，但实验条件和操作难度较高，且难以 模拟复杂实际系统的非稳态导热过程。
有限体积法
有限体积法是一种将连续的求解域离散化为 有限个小的体积，通过求解每个体积的近似 解来逼近原问题的数值解法。
有限体积法的基本思想是将导热问题分解为 若干个小的体积，每个体积具有简单的几何 形状和边界条件，然后通过求解每个体积的 近似解来逼近原问题的解。这种方法在处理 复杂的几何形状和边界条件时具有较高的精
度和可靠性。
CHAPTER

4-3非稳态导热的数值计算


t k 1 i 1
)
tik
已知k时层的温度值，求k+1 时层的温度值要联立求解方 程组，即求解复杂，但无条件稳定（、x的取值不受 限制）。
三、边界节点的离散方程
t
1. 第三类边界条件：已知tf、h
tf、h
L
节点的热平衡：
N-1点导入 对流换热传 N节点内
+
=
N点的热量 入N点的热量 能的增量 0 1 N 1 N x
2t ( x2 )i,k
tk i 1
tk i 1
2tik
(x)2
节点 时层
( t
)i,k
t k 1 i
tik
空间用中心差分格式 时间用向前差分格式
将上面两式代入微分方程:
t k 1 i
tik
a
tk i 1
tk i 1
(x)2
2tik
有缘学习更多＋谓ｙｇｄ３０７６考证资料或关注桃报：奉献教育（店铺）
tik 1
k N
稳定性条件: 1 2Fo Bi 2Fo 0
或:
Fo 1 2Bi 2
2. 第二类边界条件：已知qw
建立节点0的差分方程（显式格式） t
t k 1 0
2Fo t1k
(1
2Fo)t0k
tf、
L
稳定性条件:
1 2Fo 0
0 1 N 1 N x 绝热
THANKS
非稳态导热 的数值计算
讨论： 一维、无内热源、常物性、非稳态导热
t f (x, )
t
a
2t x 2
一、显式差分格式
1. 内节点
k 1
1) 离散化: t f (x, )
k

第4章-非稳态导热的分析与计算-简化


h

第4章 非稳态导热的分析与计算
4.2 对流边界条件下的一维非稳态导热／分析解
x, x 2 a Cn exp n 2 cos n 0 n 1
4sin n Cn 2n sin 2n
t |x h t |x t x
第4章 非稳态导热的分析与计算
4.2 对流边界条件下的一维非稳态导热／数学模型

完整的数学模型：
t 2t 控制方程： a 2 x
0 x , 0
初始条件： cV t t0
非稳态导热过程所传递的最大热量
第4章 非稳态导热的分析与计算
4.2 对流边界条件下的一维非稳态导热／吸热量
从初始时刻开始的某时间段的吸热量：
Q c t x, t0 dV cA t x, t0 dx V V
第4章 非稳态导热的分析与计算
4.2 对流边界条件下的一维非稳态导热／分析解
x, 0
x x 2 f Fo, Bi, Cn exp n Fo cos n n1
当Fo>0.2时，取Cn= Fo>0.2：

当Fo>0.2后：
x, x cos 1 0,
——θ(x,τ)与θ(0,τ) 的比值却与τ无关，仅取决于平壁的几
何位置（x/δ）和Bi数 ——初始条件的影响已经消失：正规状况阶段
第4章 非稳态导热的分析与计算
4.2 对流边界条件下的一维非稳态导热／分析解


几何上：平壁
物理上：沿高度和宽度方向的换 热均匀一致
  1. 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。


由式(4-1)可得
dt d ' d ( )(e cV
'


A cV
)
所以导热体在单位时间内传递给流体的热量为
dt Q cV 'Ae d

A cV
W
(4-2)
因导热体被冷却，故dt／dτ＜0，加负号以使Q 为正值。

利用上式，可得导热体在τ＝0到τ＝τ时 间内传入流体的总热量：
Q Qd cV (1 e
' 0


A cV
)J
(4-3)
二、计算判断

毕渥数的定义：
L Bi 1
L

毕渥数的物理意义：固体的导热热阻与对 流换热热阻之间的对比关系。
V A 0.1 M Biv
* 负号是由于dt为负值 令过余温度θ＝t－tf ，则dt＝dθ，代入上式得：
d
A d cV
当V、A、α、ρ、c等为已知定值时，对上式积 分得： d A d ' 0 d A cV d cV A ln cV '
第四章 非稳态导热的分析计算
§4-1 概述
1.定义：温度场随时间变化
2.分类：* 周期性非稳态导热
* 非周期性非稳态导热(瞬态导热)
3.目的：* 在加热或冷却时，确定物体内部某一 位置达到预定温度所需要的时间，以及在该时间 内物体吸收或放出的热量；
* 对物体加热或冷却一定时间后，确定 物体内部的温度分布和温度场随时间的变化率


'

t tf 上式是采用集总参数法求解非稳态导热问题的 基本公式，可用于已知温度求时间，或反之。

ttf
'

e
A cV
(4-1)

当时间τ＝ρcV／(αA)

t tf ' e ' t tf

A cV
e 0.368 36.8%
1
即导热在此时的过余温度θ已下降到初始过余 温度θ′的36.8%
与流体接触 的面积

大平壁: M=1 长圆柱(正方形长柱体）: M= １／２ 球（正立方体）: M=１／３

作业：4-1；4-2；4-3

一、基本计算公式
初始温度为t' 被周围温度为tf的流体冷却 换热系数α为定值 导热体的平均温度t

经dτ时间后，由于散热，温度下降dt。 由能量平衡，散热量＝△导热体本身能量，即：
dt A(t t f ) cV (散热) d
dt A(t f t ) cV (吸热) d

一温度均匀的物体，两侧被具有 恒定温度tf的高温介质所包围
§4-2 集总参数分析法
当所需求解的温度仅为时间 τ的函数而与坐标无关， 即 t＝f(τ) 条件： 导热热阻<<对流换热热阻 集 总 参 数 分 析 法 (Lumped Parameter Analysis)：壁内各 处温度相差不大，温度梯度 极小，可以把整个导热系数 看作一个处于平均温度下的 物体。
ρcV／（αA)称为时间常数τc 如果导热体的热容量(ρcV)小，换热条件好(αA ) 大，则单位时间所传递的热量大，导热体的温 度变化快，将使导热体的温度迅速接近流体温 度。

当τ＝4τc= 4ρcV／（αA)时，则：


'
e
A cV
e
4
0.0183 1.83%
工程上习惯认为，τ＝4τc时导热体已达到热平 衡状态。 时间常数关系到测温仪表的响应时间。