第四章 非稳态导热的分析计算
合集下载
传热学-第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
‹# ›
第4章 非稳态导热的计算与分析
4.1 概述
非稳态导热的分类： ——周期性的非稳态导热（periodic unsteady heat conduction）：由于边界条件（或内热源）随时间呈周 期性变化，使物体内的温度场也随时间按周期性规律变 化，这种状况通常称为准稳态
19
‹# ›
4.2.1 平壁内非稳态过程的基本特征
整个瞬态导热过程可以分为两个阶段： 初始阶段（initial regime）：也称为非正规状况阶段，
指在穿透时刻之前阶段，此时平壁内的温度分布主要受 初始温度分布t0的影响。
正规状况阶段（regular regime）：穿透时刻之后，非稳态 过程进行到一定的程度，平壁初始温度分布的影响逐渐消失，此 后不同时刻的温度分布主要受热边界条件的影响。这个阶段的非 稳态导热称为正规状况阶段。
第4章 非稳本态节导内热容的结计束 算与分析
1
• 稳态导热是一种理想化的情况 • 受环境温度变化的影响，生活和工程中真正意义上的稳 态导热是不存在的 • 只是对工程中的某些问题，忽略温度随时间变化所造成 的影响、误差不大，而将其简化为稳态导热
2
‹# ›
• 生活和工程中还存在着大量的不能简化为稳态导热的现 象和问题，其中物体内的温度明显随时间而变化 • ——冷冻食品的解冻过程 • ——烘烤食品（花生米、蛋糕等点心） • ——热处理工艺中金属在高温火炉内的加热以及加热后 在水或空气中的冷却过程等 • ——焖井过程热量在地层内的扩散过程
第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
一、数值解法的基本思想 用导热问题所涉及的空间和时间区域内有限个离散 称为节点 节点) 点(称为节点)的温度近似值来代替物体内实际连续的温 度分布， 度分布 ， 将连续温度分布函数的求解问题转化为各节 点温度值的求解问题， 点温度值的求解问题 ， 将导热微分方程的求解问题转 化为节点温度代数方程的求解问题。 化为节点温度代数方程的求解问题。 数值解法的基本内容与步骤: 数值解法的基本内容与步骤:
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章-非稳态导热的计算与分析
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
非稳态导热分析解法课件
复杂边界条件和几何形状
非稳态导热问题常常涉及到复杂的边界条件和几何形状，给分析带来很大挑战。未来发展需要研究更高效的数值方法 ，以处理更复杂的导热问题。
多物理场耦合
许多实际导热问题涉及到多物理场的耦合，如热-力、热-流体等。未来发展需要研究多物理场耦合的非稳态导热问题 ，以提高对复杂系统的理解和预测能力。
高效能材料和新能源技术
随着高效能材料和新能源技术的发展，非稳态导热问题将更加复杂和多样化。未来发展需要加强与相关 领域的交叉融合，以应对不断出现的新的挑战和机遇。
核能利用
在核能利用中，非稳态导热分析可用于研究反应堆的冷却系统、核废料的处理和存储等。 通过优化导热性能，可以提高核能系统的安全性和稳定性。
风能利用
在风能利用中，非稳态导热分析可用于研究风力发电机的散热性能和风能转换效率。通过 改进导热设计，可以提高风能发电的经济性和可靠性。
非稳态导热面临的挑战和未来发展方向
物理模拟实验
物理模拟实验是通过模拟实际系统的物理过程来研究其行为的方法。
在非稳态导热分析中，物理模拟实验通常采用加热棒、散热片等模拟导热过程，通 过测量温度场、热流密度等参数来研究非稳态导热规律。
物理模拟实验具有直观、可重复性高等优点，但实验条件和操作难度较高，且难以 模拟复杂实际系统的非稳态导热过程。
有限体积法
有限体积法是一种将连续的求解域离散化为 有限个小的体积，通过求解每个体积的近似 解来逼近原问题的数值解法。
有限体积法的基本思想是将导热问题分解为 若干个小的体积，每个体积具有简单的几何 形状和边界条件，然后通过求解每个体积的 近似解来逼近原问题的解。这种方法在处理 复杂的几何形状和边界条件时具有较高的精
度和可靠性。
CHAPTER
非稳态导热问题常常涉及到复杂的边界条件和几何形状，给分析带来很大挑战。未来发展需要研究更高效的数值方法 ，以处理更复杂的导热问题。
多物理场耦合
许多实际导热问题涉及到多物理场的耦合，如热-力、热-流体等。未来发展需要研究多物理场耦合的非稳态导热问题 ，以提高对复杂系统的理解和预测能力。
高效能材料和新能源技术
随着高效能材料和新能源技术的发展，非稳态导热问题将更加复杂和多样化。未来发展需要加强与相关 领域的交叉融合，以应对不断出现的新的挑战和机遇。
核能利用
在核能利用中，非稳态导热分析可用于研究反应堆的冷却系统、核废料的处理和存储等。 通过优化导热性能，可以提高核能系统的安全性和稳定性。
风能利用
在风能利用中，非稳态导热分析可用于研究风力发电机的散热性能和风能转换效率。通过 改进导热设计，可以提高风能发电的经济性和可靠性。
非稳态导热面临的挑战和未来发展方向
物理模拟实验
物理模拟实验是通过模拟实际系统的物理过程来研究其行为的方法。
在非稳态导热分析中，物理模拟实验通常采用加热棒、散热片等模拟导热过程，通 过测量温度场、热流密度等参数来研究非稳态导热规律。
物理模拟实验具有直观、可重复性高等优点，但实验条件和操作难度较高，且难以 模拟复杂实际系统的非稳态导热过程。
有限体积法
有限体积法是一种将连续的求解域离散化为 有限个小的体积，通过求解每个体积的近似 解来逼近原问题的数值解法。
有限体积法的基本思想是将导热问题分解为 若干个小的体积，每个体积具有简单的几何 形状和边界条件，然后通过求解每个体积的 近似解来逼近原问题的解。这种方法在处理 复杂的几何形状和边界条件时具有较高的精
度和可靠性。
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 对流边界条件下的一维非稳态导热／分析解
几何上：平壁
物理上：沿高度和宽度方向的换 热均匀一致
传热学基础(第二版)第四章教学课件非稳态导热
Lctptw
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
与式
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
与式
无限大平板内含热源非稳态导热微分方程
无限大平板内含热源非稳态导热微分方程一、概述热传导是研究物体内部热量传输规律的重要领域之一，其数学描述主要通过热传导方程来完成。
在实际工程问题中，除了简单的稳态导热问题外，还存在着非稳态导热问题。
本文将针对无限大平板内含热源的非稳态导热微分方程进行解析和探讨。
二、问题描述考虑一块无限大的平板，在其内部存在一个热源。
热源和平板内的热量传导过程满足非稳态导热方程。
我们的目标是求解该问题下的温度分布，并分析其规律和特点。
三、问题分析1.建立问题的数学模型，首先需要考虑热传导方程的具体形式，以及边界条件和初始条件的设定。
2.根据问题的特点，选择合适的数学方法进行求解，可能涉及到偏微分方程的解析解或者数值解的方法。
3.对问题的解进行讨论和分析，探讨其物理意义和对实际问题的指导作用。
四、数学建模考虑平板在x、y、z三个方向上的温度分布为T(x,y,z,t)，其中t为时间。
根据热传导方程和热源的作用，可以列出平板内部温度分布的非稳态导热方程。
∂T/∂t = α(∂^2T/∂x^2+∂^2T/∂y^2+∂^2T/∂z^2)+ Q(x,y,z,t)其中α为热导率，Q(x,y,z,t)为热源的分布情况。
在这个方程中，除了热源的影响外，还有热传导导致的温度变化。
边界条件和初始条件可以根据实际情况进行设定，通常包括平板表面的温度、热源的功率和分布等信息。
五、数学求解对于这类非稳态导热问题，通常需要借助数值方法进行求解。
常见的方法包括有限差分法、有限元法等。
在进行数值求解时，需要考虑网格的划分和时间步长的选择，以保证数值解的稳定性和精度。
通过数值方法求解后，可以得到平板内部温度的分布情况，进而分析其变化规律。
六、数学分析根据所得的温度分布情况，可以对问题进行进一步的分析和讨论。
可以分析热源对平板温度分布的影响，或者研究热传导导致的温度变化规律。
通过分析结果，可以得出一些对实际问题有指导意义的结论，为工程设计和应用提供理论支持。
非稳态导热的分析计算(最全版)PTT文档
W (4-2) 当τ＝4τc= 4ρcV／（αA)时，则：
时间常数关系到测温仪表的响应时间。 被周围温度为tf的流体冷却
令过余温度θ＝t－tf ，则dt＝dθ，代入上式得： 由能量平衡，散热量＝△导热体本身内能，即： * 对物体加热或冷却一定时间后，确定物体内部的温度分布和温度场随时间的变化率 被周围温度为tf的流体冷却
由式(4-1)可得
A
A
dt d ' d (e cV ) '( A )(e cV )
d d d
cV
所以导热体在单位时间内传递给流体的热量为
A
Q cV dt 'Ae cV d
W (4-2)
因导热体被冷却，故dt／dτ＜0，加负号以使Q
为正值。
利用上式，可得导热体在τ＝0到τ＝τ时
d A d cV
当V、A、α、ρ、c等为已知定值时，对上式积
分得：
d A d cV
'
d
A cV
0 d
ln
'
A cV
'
ttf t' t f
A cV
e
上式是采用集总参数法求解非稳态导热问题的 基本公式，可用于已知温度求时间，或反之。
当时间τ＝ρcV／(αA)
t tf ' t' tf
一、基本计算公式
初始温度为t' 被周围温度为tf的流体冷却 换热系数α为定值 导热体的平均温度t
经dτ时间后，由于散热，温度下降dt。
由能量平衡，散热量＝△导热体本身内能，即：
A(t
θt
f
)
cV
dθt
d
(散热)
A(t f
t)
时间常数关系到测温仪表的响应时间。 被周围温度为tf的流体冷却
令过余温度θ＝t－tf ，则dt＝dθ，代入上式得： 由能量平衡，散热量＝△导热体本身内能，即： * 对物体加热或冷却一定时间后，确定物体内部的温度分布和温度场随时间的变化率 被周围温度为tf的流体冷却
由式(4-1)可得
A
A
dt d ' d (e cV ) '( A )(e cV )
d d d
cV
所以导热体在单位时间内传递给流体的热量为
A
Q cV dt 'Ae cV d
W (4-2)
因导热体被冷却，故dt／dτ＜0，加负号以使Q
为正值。
利用上式，可得导热体在τ＝0到τ＝τ时
d A d cV
当V、A、α、ρ、c等为已知定值时，对上式积
分得：
d A d cV
'
d
A cV
0 d
ln
'
A cV
'
ttf t' t f
A cV
e
上式是采用集总参数法求解非稳态导热问题的 基本公式，可用于已知温度求时间，或反之。
当时间τ＝ρcV／(αA)
t tf ' t' tf
一、基本计算公式
初始温度为t' 被周围温度为tf的流体冷却 换热系数α为定值 导热体的平均温度t
经dτ时间后，由于散热，温度下降dt。
由能量平衡，散热量＝△导热体本身内能，即：
A(t
θt
f
)
cV
dθt
d
(散热)
A(t f
t)
第四章 非稳态导热
边界条件 常用的三类
求解方法：
分析解法、数值解法、实验模拟法、图解法等
首先求出温度分布t f ( x, y, z, )； 然后由傅里叶定律算出 各点的瞬时热流量。
9
10
第四章 / 第三节 非稳态导热
二、非稳态导热问题的求解及诺模图
以第三类边界条件下的 无限大平壁的非稳态导热为例
11
第四章 / 第三节 非稳态导热
工程中：
机器启动、停机、变工况时部件的导热过程； 冶金、热加工、热处理工艺中工件的加热及冷却过程等； 石油工程中钻井、焖井、采油等过程中热量在地层内的扩散过程。
具有实际意义。
2
第三节
本节讨论:
非稳态导热
——基本概念和特点
——非稳态导热问题的求解及诺模图
——集总参数法
3
（一）无限大平壁的分析解及诺模图 1、平壁内温度分布的求解
物理模型：
常物性、无内热源、一维平壁
数学模型：
两侧受流体对称加热，中心面为对称面，
只需研究半厚的平壁。 将坐标原点置于平壁中心面，建立如图 直角坐标系。
t 2t a 2 x
定解条件
0 x , 0
原因：各处本身温度变化要积蓄(或放出)热量。
6
第四章 / 第三节 非稳态导热
一、 概述
3、研究目的
（1）确定非稳态过程中的温度场
——物体中某个部位到达某个预定温度所需的时间； ——在预定时间内物体可以达到的温度； ——物体的温度对时间的变化速率。 （2）确定非稳态过程的热流量
——物体在某一瞬间每一位置处的热流密度；
—出现在特征数中的几何尺度 —不同情况下，不同形状的物体特征长度是不同的。 Fo 数 、 Bi数称为特征数，习惯上又称准则数， 具有特定的物理意义。
第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。 两侧冷却情况相同、温度分布 对称。中心为原点。
传热学课件第四章非稳态导热
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
第四章 非稳态导热
| x 0 0 x | x h | x x
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
第四章 / 第三节 非稳态导热
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
第四章 / 第三节 非稳态导热
传热学 第4章-导热问题的数值解法
第四章 导热问题的数值解法1、重点内容： ① 掌握导热问题数值解法的基本思路；② 利用热平衡法和泰勒级数展开法建立节点的离散方程。
2、掌握内容：数值解法的实质。
3、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。
由前述3可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。
但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。
随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种：（1） 有限差分法 （2）有限元方法 （3）边界元方法数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。
如：几何形状、边界条件复杂、物性不均、多维导热问题。
分析解法与数值解法的异同点：1、 相同点：根本目的是相同的，即确定① t=f(x ，y ，z)；② ),,,(τz y x g Q =。
2、 不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。
§4—1 数值求解的基本思路及稳态导热内节点离散方程的建立一、 解法的基本概念1、 实质对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。
该方法称为数值解法。
这些离散点上被求物理量值的集合称为该物理量的数值解。
2、基本思路：数值解法的求解过程可用框图4-1表示。
由此可见：1）物理模型简化成数学模型是基础； 2）建立节点离散方程是关键；3）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 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时导热体已达到热平 衡状态。 时间常数关系到测温仪表的响应时间。