非稳态导热习题
第三章 非稳态导热习题
例一腾空置于室内地板上的平板电热器，加在其上的电功率以对流换热和辐射换热的方式全部损失于室内。电热器表面和周围空气的平均对流换热系数为h ，且为常数，室内的空气温度和四壁、天花板及地板的温度相同，均为t f 。电热器假定为均质的固体，密度为ρ，比热为c ，体积为V , 表面积为A ，表面假定为黑体，因其导热系数足够大，内部温度均布。通电时其温度为t 0。试写出该电热器断电后温度随时间变化的数学描述。
[解] 根据题意，电热器内部温度均布，因此可用集中参数分析法处理。 电热器以辐射换热方式散失的热量为：
44r f ()A T T σΦ=- （1）
以对流换热方式的热量为：
c f ()hA T T Φ=- （2）
电热器断电后无内热源，根据能量守恒定律，散失的热量应等于电热器能量的减少。若只考虑电热器的热力学能
r c d d T
cV
ρτ
-Φ-Φ= （3） 因此，相应的微分方程式为：
44f f d ()()d T
A T T hA T T cV
σρτ
-+-=- （4） 初始条件为：
τ=0， t =t 0 （5）
上述两式即为该电热器断电后温度随时间变化的数学描述。
例 电路中所用的保险丝因其导热系数很大而直径很小可视为温度均布的细长圆柱体，电流的热效应可视为均匀的内热源。如果仅考虑由于对流换热的散热量，保险丝表面和温度为t f 的周围空气之间的平均对流换热系数为h ，且为常数。试求该保险丝通电后温度随时间的变化规律。
[解] 根据题意，保险丝内部温度均布，因此可用集中参数分析法处理。 保险丝表面以对流换热方式散失的热量为：
c f ()hA T T Φ=- （1）
保险丝的内热源为：
Q 0=IR 2 （2）
式中：I ——保险丝通过的电流，（A ）； R ——保险丝的电阻，Ω。
根据能量守恒，散失的热量与内热源所转变成的热量的和应等于保险丝能量的变化。若只考虑保险丝的热力学能
c 0d d T
Q cV
ρτ
-Φ+= （3）
因此，相应的微分方程式为：
2f d ()d T
hA T T I R cV
ρτ
--=- （4） 初始条件为：
τ=0， t =t f （5）
上述两式即为该保险丝通电后温度随时间变化的数学描述。
令2f I R
t t hA
θ=--，则上述微分方程改写为
d d ρθ
θτ
=-
cV hA （6）
该微分方程的解为
θθτθρ=-
=-00exp()exp(BiFo)hA
cV
（7） 以温度t 表示该解
τρ--=---22f 0f (）exp()I R I R hA
t t t t hA hA cV
（8）
由初始条件τ=0， t 0= t f ，该式可写为
τρ-=--2f [1exp()]I R hA
t t hA cV
（9）
上式即为该保险丝通电后温度随时间的变化规律，从中可以看出内热源对保险丝的温度
变化的作用。
例 一块厚10 mm 的纯铝板置于温度为10 ℃的空气中，铝板和空气之间的平均对流换
热系数h =10 W/(m 2
·K)，且为常数。
求该铝板从100 ℃降到20 ℃所需时间及当时的热流密度。 [解] 求解瞬态导热问题，应先计算比渥准则Bi 的数值，确定是否能采用简单的集总参
数法。查取铝的物性参数，密度ρ=2702 kg/m 3
，比热容c =903 J/kg ，导热系数λ=237 W/(m·K)。
λ-??=
==??4100.01'
Bi 2.11102372'
hV A A A （1） Bi<, 可用集总参数法计算。
0exp()hA
cV
θθτρ=-
（2） 102'
2010(10010)exp()27029030.01'
A A τ?-=--
??? （3）
τ=2680 s （4）
铝板从100 ℃降到20 ℃时，铝板的表面温度，空气温度，铝板和空气之间的平均对流换热系数h 均为已知，因此热流密度可用牛顿冷却公式计算。
q =h ( t- t f )=10×(20-10)=100 W/m 2
（5）
例 用球形热电偶接点作动态温度测量时，对热电偶的响应速度有一定要求。现要求一个初温为t 0的球形热电偶与温度为t f 的被测流体接触后，在 1 s 内所指示的过余温度比
f 00f
95%t t t t θθ-==-。现有一铜-康铜球形热电偶接点，它与被测流体之间的对流换热表面传热系数h =50 W/(m 2
·K)，且为常数。试求该球形热电偶接点的最大允许半径r 0。
[解] 求解瞬态导热问题，应先计算比渥准则Bi 。查取铜的物性参数，密度ρ=8954 kg/m 3
，比热容c =384 J/kg ，导热系数λ=398 W/(m·K)；查取康铜的物性参数，密度ρ=8922
kg/m 3
，比热容c =410 J/kg ，导热系数λ=22 W/(m·K)。球形热电偶接点是这两种材料的熔
化物，因此取平均值，密度ρ=8938 kg/m 3
，比热容c =397 J/kg ，导热系数λ=210 W/(m·K)。
λ=
=00
50Bi 210
hr r （1）
因半径r 0未知，比渥准则Bi 的数值无法计算。但可假定Bi<，先用集总参数法计算， 然后进行较核。
0exp()hA
cV
θθτρ=-
（2）
ππ-?-?=-
=??
2
5
030
504 4.206100.95exp(1)exp(
)489383973
r r r （3）
r 0=×10-4 m (4)
校核r 0，λ==00
50Bi 0.1210
hr r （5）
因为比渥准则Bi<<，上述分析计算合理。
讨论：求解瞬态导热问题，应先计算比渥准则Bi ，一旦Bi<，就可以用简单的集总参数法计算，但是Bi 数值的确定需要先知道定型尺寸的数值。本题中定型尺寸的数值是所求对
象，因此只能先假定Bi<，能用集总参数法计算，计算完后需要根据算出的定型尺寸校核集总参数法的应用条件Bi<是满足的。
例某种电路中所用的保险丝的直径为0.5 mm,长20 mm ，导热系数λ=20 W/(m·K)，热
扩散率a =5×10-5 m 2
/s,电阻为Ω,熔点为900 ℃。如果仅考虑由于对流换热的散热量，保险
丝表面和温度为20 ℃的周围空气之间的平均对流换热系数为10 W/(m 2
·K)，且为常数。试确定该保险丝通过2 A 的电流后多少时间会熔断。
[解] 该保险丝因其导热系数较大而直径很小可视为温度均布的细长圆柱体，电流的热效应可视为均匀的内热源。瞬态导热问题先计算比渥准则Bi 的数值。
πλπ?=≈
(1)
根据解析题,
2f [1exp()]I R hA t t hA cV
τρ-=-- (2)
代入具体数值
22
5
20.8
1090020[1exp()]2010/4510
dl
dl
d l πτππ-??-≈--
?? (3)
880=×104
[1-exp τ)] (4) τ= s (5)
例 将直径为30 mm 、初温为20 ℃的生红肠放入温度为180 ℃的烘箱中烤熟。假定生
红肠的密度ρ=960 kg/m 3
，比热容c =5000 J/kg ，导热系数λ= W/(m·K)，仅考虑由于对流
换热的加热量，红肠和烘箱中空气之间的平均对流换热系数为30 W/(m 2
·K)，且为常数。试求生红肠放入烘箱中10 min 时红肠的中心温度。
[解] 红肠可视为细长圆柱体。瞬态导热问题先计算比渥准则Bi 的数值。
πλπ???=≈>???2300.03/4Bi 0.10.90.03hV l A l
(1)
因此本题不能用集总参数法计算，只能用查计算线图（海斯勒图）的方法。
λ?=
==300.015
Bi 0.50.9hR (2) λτρ?=
==??22
0.9600
Fo 0.596050000.015cR (3)
查计算线图（海斯勒图）得
m m 01800.5620180
t θθ-=≈- (4) 解得生红肠放入烘箱中10 min 时的中心温度 t m =90.4 ℃ (5)
例 直径为400 mm 、初温为20 ℃的钢棒放入温度为600 ℃的炉中加热。钢棒的密度ρ=7833 kg/m 3，比热容c =465 J/kg ，导热系数λ=54 W/(m·K)，仅考虑由于对流换热的加
热量，钢棒与炉中气体之间的平均对流换热系数为130 W/(m 2
·K)，且为常数。
试求钢棒中心温度达到400 ℃时所需的时间，并确定此时钢棒的表面温度。 [解] 钢棒可视为细长圆柱体。瞬态导热问题先计算比渥准则Bi 的数值。
πλπ???=≈>???21300.4/40.1540.4hV l Bi A l
(1)
因此本题不能用集总参数法计算，只能用查计算线图（海斯勒图）的方法。
λ?=
==1300.2
Bi 0.481454
hR (2) 钢棒中心温度达到400 ℃时，
m 04006000.344820600
θθ-=≈- (3) 查计算线图（海斯勒图）得Fo 准则的数值为
λτττρ-?==
==???4
22
54Fo 1.3 3.7061078334650.2
cR (4) 解得钢棒中心温度达到400 ℃时所需的时间 τ=3508 s (5)
R R m 600
0.79400600
t θθ-=≈- (6) 解得此时钢棒的表面温度 t R =442 ℃
例 截面为1 m×1 m 的耐火砖方形长柱体，初温为20 ℃，与600 ℃的高温烟气接触，
仅考虑由于对流换热的加热量，柱体与燃气之间的平均对流换热系数为20 W/(m 2
·K)，且为
常数。耐火砖的密度ρ=2000 kg/m 3
，比热容c =960 J/kg ，导热系数λ= W/(m·K)，
试求耐火砖柱体与烟气接触120小时时方柱体的中心温度。 [解] 瞬态导热问题先计算比渥准则Bi 的数值。
λ???=
≈>???2011Bi 0.11.0741hV l
A l
(1) 因此本题不能用集总参数法计算，只能用查计算线图（海斯勒图）的方法。方形长柱体
的导热是二维导热问题，可用两个壁厚相同的无限大平壁的解的乘积求得。
δλ?=
==200.5
Bi 9.3461.07h (2) λτρδ?=
==??22
1.07432000
Fo 0.96320009600.5
c (3) 再查计算线图（海斯勒图）得
m
0.18θθ≈ (4) 因此，方形长柱体中心的过余温度比
m f m m m 0f 00
6000.180.1820600t t t t t θθ
θθ--==?=?-- (5)
最后解得120小时时方形长柱体的中心温度t m =581.2 ℃。
例 一块厚300 mm 的无限大钢板[密度ρ=7753 kg/m 3
，比热容c =486 J/kg ，导热系数λ=36 W/(m·K)]初温为900℃，突然置于20 ℃的空气气流中，空气和钢板间的对流换热系
数为400 W/( m 2
·K)。求钢板表面温度达到200 ℃时所需时间。
[解] 瞬态导热问题先计算比渥准则Bi 的数值。
λ??=
≈>?4000.3'
Bi 0.1362'
hV A A A (1) 因此不能用集总参数法，只能用查计算线图（海斯勒图）的方法。
由于中心温度未知，不能用含有傅立叶数Fo 的第一张图直接查时间，只能先用第二张图确定中心温度。
δλ?=
≈=4000.15
Bi 1.66736
h (2) δm m 200200.5420
t θθ-==- (3) 解得中心温度 t m =353.3 ℃ (4)
接下来用第一张图，Bi -1
=，中心处的过余温度比m 0353.320
0.378890020
θθ-==-，查得傅立叶数Fo≈
λττ
ρδ?=??22
36Fo=
77534860.15c (5)
最后解得钢板表面温度达到200℃时所需时间
τ=2255 s
例 一钢锭加热到400 ℃，浸在100 ℃的水中冷却。钢锭的密度ρ=7753 kg/m 3
，比热
容c =486 J/kg ，导热系数λ=36 W/(m·K)。
试求3 min 后深度为5 cm 处的温度。
[解]水沸腾时对流换热系数很大，若忽略钢锭表面的对流换热热阻，本题可视为半无限大物体在恒温边界条件下的温度分布问题。相应的温度表达式为
τ-=-w 0w (,)erf(t x t x
t t (1)
热扩散率a =λ/ρc =36/(7753×486)=×10-6
。
(,)100
erf(
erf(0.6029)=0.6061400100
t x τ-==- (2)
最后解得3min 后深度为5cm 处的温度 t = 281.8 ℃.
例 若将一以余弦函数形式的周期性温度变化加在一块很大的钢筋混凝土表面，使其温
度由35 ℃变化到90 ℃的一个完整循环需要15 min 。钢筋混凝土的密度ρ=2400 kg/m 3
，比热容c =840 J/kg ，导热系数λ= W/(m·K)
试求温度波动2小时后，距表面5cm 处的温度。
[解]本题可视为半无限大物体在周期性变化边界条件下的温度分布问题。根据已知条件，
温度波动幅度 A =(90-35)/2=27.5 ℃ 温度波动周期 T =60×15=900 s 相应的温度表达式为
w 2(,)exp()cos(x A T πθττ=--
3.379=== (1)
250.24 3.37946.81T
π
τ-=-== (2) 温度波动2小时后，距表面5cm 处的温度
(0.05,7200)27.5exp( 3.379)cos(2682.02)0.8914θ=-=- (3)
θ=t -t m (4) t m =(t max +t min )/2=(90+35)/2=62.5 ℃ (5) 因此解得温度波动2小时后，距表面5cm 处的温度 t =61.61 ℃
例 根据半无限大物体在周期性变化边界条件下的温度分布公式，试确定在表面和深度为x 处温度最大值的滞后时间。
[解] 半无限大物体在周期性变化边界条件下的温度分布公式如下：
w 2(,)exp()cos(x A T πθττ=-- (1) 深度为x
处温度出现最大值时，2cos(
1T
π
τ-=，即
20T
π
τ-= (2) 在x =0处，τ=0 (3)
在x =x 处， 2/2x T πτ=
=
(完整版)传热学期末考试试题
传热学（一） 第一部分选择题 ?单项选择题（本大题共 10 小题，每小题 2 分，共 20 分）在每小题列出的四个选项中只有一个选项是符合题目要求的，请将正确选项前的字母填在题后的括号内。 1. 在稳态导热中 , 决定物体内温度分布的是 ( B) A. 导温系数 B. 导热系数 C. 传热系数 D. 密度 2. 下列哪个准则数反映了流体物性对对流换热的影响 ?(C ) A. 雷诺数 B. 雷利数 C. 普朗特数 D. 努谢尔特数 3. 单位面积的导热热阻单位为 ( B)
A. B. C. D. 4. 绝大多数情况下强制对流时的对流换热系数 (C ) 自然对流。 A. 小于 B. 等于 C. 大于 D. 无法比较 5. 对流换热系数为 100 、温度为 20 ℃的空气流经 50 ℃的壁面，其对流换热的热流密度为（D ） A. B. C. D. 6. 流体分别在较长的粗管和细管内作强制紊流对流换热，如果流速等条件相同，则（ C） A. 粗管和细管的相同 B. 粗管内的大 C. 细管内的大 D. 无法比较 7. 在相同的进出口温度条件下，逆流和顺流的平均温差的关系为（ A） A. 逆流大于顺流 B. 顺流大于逆流 C. 两者相等 D. 无法比较
8. 单位时间内离开单位表面积的总辐射能为该表面的（A ） A. 有效辐射 B. 辐射力 C. 反射辐射 D. 黑度 9. （D ）是在相同温度条件下辐射能力最强的物体。 A. 灰体 B. 磨光玻璃 C. 涂料 D. 黑体 10. 削弱辐射换热的有效方法是加遮热板，而遮热板表面的黑度应（B ） A. 大一点好 B. 小一点好 C. 大、小都一样 D. 无法判断 第二部分非选择题 ?填空题（本大题共 10 小题，每小题 2 分，共 20 分） 11. 如果温度场随时间变化，则为。非稳态温度场
非稳态导热习题
第三章 非稳态导热习题 例3.1一腾空置于室内地板上的平板电热器，加在其上的电功率以对流换热和辐射换热的方式全部损失于室内。电热器表面和周围空气的平均对流换热系数为h ，且为常数，室内的空气温度和四壁、天花板及地板的温度相同，均为t f 。电热器假定为均质的固体，密度为ρ，比热为c ，体积为V , 表面积为A ，表面假定为黑体，因其导热系数足够大，内部温度均布。通电时其温度为t 0。试写出该电热器断电后温度随时间变化的数学描述。 [解] 根据题意，电热器内部温度均布，因此可用集中参数分析法处理。 电热器以辐射换热方式散失的热量为： 44r f ()A T T σΦ=- （1） 以对流换热方式的热量为： c f ()hA T T Φ=- （2） 电热器断电后无内热源，根据能量守恒定律，散失的热量应等于电热器能量的减少。若只考虑电热器的热力学能 r c d d T cV ρτ -Φ-Φ= （3） 因此，相应的微分方程式为： 44f f d ()()d T A T T hA T T cV σρτ -+-=- （4） 初始条件为： τ=0， t =t 0 （5） 上述两式即为该电热器断电后温度随时间变化的数学描述。 例 3.2 电路中所用的保险丝因其导热系数很大而直径很小可视为温度均布的细长圆柱体，电流的热效应可视为均匀的内热源。如果仅考虑由于对流换热的散热量，保险丝表面和温度为t f 的周围空气之间的平均对流换热系数为h ，且为常数。试求该保险丝通电后温度随时间的变化规律。 [解] 根据题意，保险丝内部温度均布，因此可用集中参数分析法处理。 保险丝表面以对流换热方式散失的热量为： c f ()hA T T Φ=- （1） 保险丝的内热源为： Q 0=IR 2 （2） 式中：I ——保险丝通过的电流，（A ）； R ——保险丝的电阻，Ω。 根据能量守恒，散失的热量与内热源所转变成的热量的和应等于保险丝能量的变化。若只考虑保险丝的热力学能 c 0d d T Q cV ρτ -Φ+= （3）
传热学传热学--第三章 第三节 一维非稳态导热问题
传热学--第三章第三节一维非稳态导热问题 §3 — 3 一维非稳态导热的分析解 本节介绍第三类边界条件下：无限大平板、无限长圆柱、球的分析解及应用。如何理解无限大物体，如：当一块平板的长度、宽度>> 厚度时，平板的长度和宽度的边缘向四周的散热对平板内的温度分布影响很少，以至于可以把平板内各点的温度看作仅是厚度的函数时，该平板就是一块“无限大”平板。若平板的长度、宽度、厚度相差较小，但平板四周绝热良好，则热量交换仅发生在平板两侧面，从传热的角度分析，可简化成一维导热问题。 一、无限大平板的分析解 已知：厚度的无限大平板，初温t0，初始瞬间将其放于温度为的流体中，而且> t0，流体与板面间的表面传热系数为一常数。 试确定在非稳态过程中板内的温度分布。 解：如图3-5 所示，平板两面对称受热，所以其内温度分布以其中心截面为对称面。对 于x 0 的半块平板，其导热微分方程：(0 （边界条件） （边界条件） 对偏微分方程分离变量求解得： （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 ）。也就是说，初始条件的影响已经消失，无论初始条件分布如何，只要 一维非稳态导热的数值计算 一、实验名称 一维非稳态导热的数值计算 二、实验内容 一块无限大平板（如图3所示），其一半厚度为L=0.1m ，初始温度T 0=1000℃，突然将其插入温度T ∞=20℃的流体介质中。平板的导热系数λ=34.89W/m ℃，密度ρ=7800 kg/m 3，比热c=0.712310 J/kg ℃，平板与介质的对流换热系数为h=233W/m 2.℃，求平板内各点的温度分布。 三、实验编程 #include for(k=0; k<=1000; k++) {for(i=1; i<=L-1; i++) for(j=1; j<=L-1; j++) {T[i][j] = T[i][j] + (a/4)*(T[i+1][j] + T[i][j+1] + T[i-1][j] + T[i][j-1] - 4*T[i][j]); } } printf(" a = %lf\n", a); printf("T[x][y] = ...\n"); for(i=0; i<=L; i++) for(j=0; j<=L; j++) {printf("%.1lf\t", T[i][j]); if(j == L) putchar(10); } return 0; } 四、运行结果 1、试分析室内暖气片的散热过程，各个环节有哪些热量传递方式？以暖气片管内走热水为例。 答：有以下换热环节及传热方式： (1) 由热水到暖气片管道内壁，热传递方式为强制对流换热； (2) 由暖气片管道内壁到外壁，热传递方式为固体导热； (3) 由暖气片管道外壁到室内空气，热传递方式有自然对流换热和辐射换热。 2、写出Nu 、Re 、Pr 、Bi 、Fo 数的表达式，并说明其物理意义。 答：（1）努谢尔特数，，它表示壁面法向无量纲过余温度梯度的大小，或反映对流换热的强弱。 （2）雷诺数，，它表示流体流动时惯性力与粘性力的相对大小。 （3）普朗特数，，它反映了流体的动量传递能力与热量 传递能力的相对大小。 （4）毕渥数，，它表示物体内部导热热阻与物体表面对 流换热热阻的比值。 （5）傅立叶数，，是非稳态导热过程的无量纲时间,它表示非稳态导热过程进行的深度。 3、热扩散系数是表征什么的物理量？它与导热系数的区别是什么？ 答：热扩散率，与导热系数一样都是物性参数， 它是表征物体被加热或冷却时，物体内各部分温度趋向均匀一致的能力。热扩散率取决于导热系数和的综合影响；而导热系数λ是反映物体的导热能力大小的物性参数。 一般情况下，稳态导热的温度分布取决于物体的导热系数，但非稳态导热的温度分布不仅取决于物体的导热系数，还取决于物体的导温系数。 4、集总参数法的适用条件是什么？满足集总参数法的物体，其内部温度分布有何特点？ 答：集总参数法的适用条件是Bi<0.1，其特点是当物体内部导热热阻远小于外部对流换热热阻时，物体内部在同一时刻均处于同一温度，物体内部的温度仅是时间的函数，而与位置无关 8、影响强制对流换热的表面换热系数的因素有哪些？ 答：影响强制对流换热的表面换热系数的因素有流态、流体的物性、换热表面的几何因素等，用函数表示为： 9、沸腾换热的临界热流密度的含义是什么？ 答：在泡态沸腾阶段时，液体温度与壁面温度之差若进 一步增大，汽泡在表面上生成、长大，随后引因浮力作用而离开表面。沸腾的液体主体温度这时有一定的过热 度，故汽泡通过液体层时还会继续被加热、膨胀，直到 逸出液面，由于气泡的大量迅速生成和它的剧烈运动，换热强度剧增，热流密度随的提高而急剧增大，直到达到热流密度的峰值，此时的热流密度称为临界热流密度。当进一步增大时，热流密度又开始下降。 14、不凝结气体对表面凝结换热强弱有何影响？ 答：不凝结气体的存在，一方面使凝结表面附近蒸汽的分压力降低，从而蒸汽饱和温度降低，使得传热驱动力即温差减小；另一方面，凝结蒸汽穿过不凝结气体层到达壁面依靠的是扩散，从而增加了阻力。因此，上述两方面原因导致凝结换热时的表面传热系数降低。 15、空气横掠垂直管束时，沿流动方向管排数越多，换热越强，而蒸汽在水平管束外凝结时，沿液膜流动方向管排数越多，换热强度降低，为什么？ 答：空气横掠垂直管束时，沿流动方向管排数越多，气流扰动越强，换热越强，而蒸汽在水平管束外凝结时，沿液膜流动方向管排数越多，凝结液膜越厚，凝结换热热阻越大，换热强度降低。 16、写出时间常数的表达式，时间常数是从什么导热问题中定义出来的？它与哪些因素有关？ 答：时间常数的表达式为， 是从非稳态导热问题中定义出来的，它不仅取决于几何参数V/A,物性参数PC ，还取决于换热条件h 。 8、其它条件相同时，同一根管子横向冲刷与纵向冲刷相比，哪个的表面换热系数大？为什么？ 答：同一根管子横向冲刷比纵向冲刷相比的表面换热系数大。因为纵向冲刷时相当于外掠平板的流动，热边界层较厚，热阻较大；而横向冲刷时热边界层较薄且在边 界层由于分离而产生的旋涡，增加了流体扰动，因而换 热增强。 23、在寒冷的北方地区，建房用砖采用实心砖还是多孔的空心砖好？为什么？ 答：采用空心砖较好，因为空心砖内部充满着空气，而空气的导热系数相对较小，热阻较大，空心砖导热性较之实心砖差，同一条件下空心砖的房间的散热量小保温性好。 25、北方深秋季节的清晨，树叶叶面上常常结霜。试问树叶上、下表面的哪一面上容易结霜？为什么？ 答：霜会容易结在树叶的上表面，因为树叶上表面朝向太空，而太空表面的温度会低于摄氏零度；下表面朝向地面，而地球表面的温度一般在零度以上。相对于下表面来说，树叶上表面向外辐射热量较多，温度下降的快，一旦低于零度时便会结霜。 27、窗玻璃对红外线几乎是不透过的，但为什么隔着玻璃晒太阳却使人感到暖和？ 答：窗玻璃对红外线几乎不透过，但对可见光则是可透过的，当隔着玻璃晒太阳时，太阳光可以穿过玻璃进入室内，而室内物体发出的红外线却被阻隔在室内，因房 间内温度越来越高，从而感到暖和。 29、用热电偶监测气流温度随时间变化规律时，应如何选择热电偶节点的大小？ 答：在其它条件相同时，热电偶节点越大，它的温度变化一定幅度所需要吸收（或放出）的热量越多，此时虽然节点换热表面积也有所增大，但其增大的幅度小于体 λ l h Nu = ν ul =Re λ l h Bi =a ν=Pr 2δ τa Fo =c a ρλ=λc ρ),,,,,,,,(l c t t u f h p f w μαρλ=t ?t ?t ?hA cV ρτ= 稳态导习题 1 固体内的一维导热问题 例1 具有均匀内热源强度q v 的无限大平壁处于稳态导热，其厚度为2δ，导热系数λ为常数，两侧壁温各自均布，分别为 t w1和t w2，试求该平壁内的温度分布表达式。 解: 根据题意，x 坐标的原点取平壁的中心线，描述该平壁内稳态导热现象的微分方程式为： 2v 2d 0d q t x λ += （1） 边界条件： x= -δ: t=t w1 x= δ: t=t w2 （2） 移项后积分该微分方程式两次可得其通解 v 1d d q t x C x λ =?+ 2v 122q t x C x C λ =?++ （3） 代入边界条件 2v w112()()2q t C C δδλ=??+?+ （4） 2v w2 122q t C C δδλ=?++ （5） 式（4）+式（5） 2 w1w2v 22δλ+= +t t q C （6） 式（4）－式（5） w2w1 12t t C δ ?= （7） C 1和C 2代入微分方程式的通解式（3）后得到壁内的温度表达式 22v w2w1w2w1(2)222 δλδ?+= ?++q t t t t t x x （8） 例2具有均匀内热源q v 的无限大平壁处于稳态导热，其厚度为2δ，导热系数λ为常数，两侧壁温各自均布且相同，均为t w ，试求该平壁内的温度分布表达式。 解: 根据题意，导热微分方程式同上题。由于两侧壁温相同，是一种对称情况，因此只需求解一半的求解域即可，x 坐标的原点取平壁的中心线。描述该平壁内稳态温度场的微分方程式为： 2v 2d 0d q t x λ += （1） 边界条件：x=0: d 0d t x = 第三章非稳态导热分析解法 本章主要要求： １、重点内容： ① 非稳态导热的基本概念及特点； ② 集总参数法的基本原理及应用； ③ 一维及二维非稳态导热问题。 2 、掌握内容： ① 确定瞬时温度场的方法； ② 确定在一时间间隔内物体所传导热量的计算方法。 3 、了解内容：无限大物体非稳态导热的基本特点。 许多工程问题需要确定：物体内部温度场随时间的变化，或确定其内部温度达某一极限值所需的时间。如：机器启动、变动工况时，急剧的温度变化会使部件因热应力而破坏。因此，应确定其内部的瞬时温度场。钢制工件的热处理是一个典型的非稳态导热过程，掌握工件中温度变化的速率是控制工件热处理质量的重要因素；金属在加热炉内加热时，要确定它在炉内停留的时间，以保证达到规定的中心温度。 §3—1 非稳态导热的基本概念 一、非稳态导热 1 、定义：物体的温度随时间而变化的导热过程称非稳态导热。 2 、分类：根据物体内温度随时间而变化的特征不同分： 1 2 ）物体的温度随时间而作周期性变化 如图 3-1 所示，设一平壁，初值温度 t 0 ，令其左侧的表面温 度突然升高到 并保持不变，而右侧仍与温度为 的空气接触，试分 析物体的温度场的变化过程。 首先，物体与高温表面靠近部分的温度很快上升，而其余部分仍 保持原来的 t 0 。 如图中曲线 HBD ，随时间的推移，由于物体导热温度变化波及范 围扩大，到某一时间后，右侧表面温度也逐渐升高，如图中曲线 HCD 、 HE 、 HF 。 最后，当时间达到一定值后，温度分布保持恒定，如图中曲线 HG （若 λ=const ，则 HG 是直线）。 由此可见，上述非稳态导热过程中，存在着右侧面参与换热与不参 与换热的两个不同阶段。 （ 1 ）第一阶段（右侧面不参与换热） 温度分布显现出部分为非稳态导热规律控制区和部分为初始温度区的混合分布，即：在此阶段物体温度分布受 t 分布的影响较大，此阶段称非正规状况阶段。 （ 2 ）第二阶段，（右侧面参与换热） 当右侧面参与换热以后，物体中的温度分布不受 to 影响，主要取决于边界条件及物性，此时，非稳态导热过程进入到正规状况阶段。正规状况阶段的温度变化规律是本章讨论的重点。 2 ）二类非稳态导热的区别：前者存在着有区别的两个不同阶段，而后者不存在。 3 、特点； 非稳态导热过程中，在与热流量方向相垂直的不同截面上热流量不相等，这是非稳态导热区别于稳态导热的一个特点。 计算传热学程序报告 题目：一维非稳态导热问题的数值解 姓名： 学号：学院：能源与动力工程学院专业：工程热物理日期：2014年5月25 日 一维非稳态导热问题数值解 求解下列热传导问题: 1. 方程离散化 对方程进行控制体积分得到: 非稳态项：选取T随x阶梯式变化，有扩散项：选取一阶导数随时间做显示变化，有进一步取T随x呈分段线性变化，有 T ()e x T E T P (T \ T P T W ，( )w (x)x ( X) 整理可以得到总的离散方程为: 2. 计算空间和时间步长 取空间步长为: h=L/N 网格Fourier数为： F。一^ —V (小于时稳定) x x 时间步长为： 3. 建立温度矩阵与边界条件 T=o nes(N+1,M+1) T(:,1)=Ti(初始条件温度都为0) T(1,:)=To(边界条件x=0处温度为1) T(N+1,:)=Te(边界条件x=L处温度为0) 4. 差分法求解温度 由离散方程可得到： 转化为相应的温度矩阵形式： 5. 输入界面 考虑到方程的变量，采用inputdlg函数设置5个输入变量，对这5个变量设置了默认值，如图1所示。在计算中可以改变不同的数值，得到不同的结果，特别注意稳定条件的临界值是。根据设置的默认值，得到的计算结果如图2所示。图1matlab变量输入界面 图2 默认值的计算结果 6. 结果分析根据上面的分析，给出了程序的输入界面，以及默认值状态下的数值解。可以通过改变不同的输入值，得到需要的分析结果，总结出了下面4 点结论： （1）取F o=，得到一维非稳态导热结果如下图所示 图2F。二时一维非稳态导热 从图中可以看出，对于长度L=1 的细杆，初始时刻t=0 时温度为0，边界条件 x=0时，T=1,边界条件x=1时，T=0。随着时间的增加，温度从x=0通过导热的形式传递到x=1，不同时刻不同位置杆的温度都不同，并且随着时间的增加，杆的温度也逐渐增加。 （2）取F o=，可以得到不同位置的温度响应曲线，如下图所示 图3F o=时不同X位置处的温度响应 图中红色曲线代表x=位置的温度瞬态响应，黑色曲线代表x=位置的温度瞬态响应，蓝色曲线代表X=位置的温度瞬态响应。从图中可以看出，随着X的增加，曲线与X 轴的交点值越大，温度开始传递到该位置的所需的时间越长。随着x 的增加，温度响应曲线的变化速率越慢，最终的达到的温度也越低。 （3）取F o=，得到不同位置的温度响应曲线如下图所示 图4F o=时不同X位置处的温度响应 图中三条曲线分别是X=，X=，x=位置的温度瞬态响应。与图3的F o=进行对比，两种情况下的F o值不同，F o值越大表明热扩散系数的值越大。从图中可以 看出热扩散系数对于导热的影响，尸0=时，与F o=相比较，各位置开始响应时所需的时间较长，而且各位置响应曲线的变化速率较小，最终的达到的温度也较低，说明了热扩散系数越小，热传导越慢，传递效率越低。 （4）取F o=，得到非稳定的数值解如图所示 图5F o二时一维非稳态导热 图6F o=时不同X位置处的温度响应 从图中可以看出，对于显示格式的离散方程，并不是所有的F o值都能得到有意义的解，必须要求F o＜时才能得到稳定的数值解，当F o＞时，会出现物理上不真实的解。 一、名词解释 稳态温度场：物体内各点温度不随时间变化的温度场。 等温面 ：温度场中同一瞬间温度相同点组成的面。 热扩散率（或导温系数）：表征物体内部温度趋于一致的能力，为c ρλα= 肋效率：肋片的实际散热量与假设整个肋片表面处于肋基温度下的散热量之比。 二、解答题和分析题 1．写出傅里叶定律的一般形式的数学表达式，并说明其中各个符号的意义。 答：傅里叶的一般表达式为：n n t gradt q ??-=-=λλ。 其中：q 是热流密度矢量；λ为导热系数，它表示物质导热本领的大小；gradt 是空间某点的温度梯度；n 是通过该点的等温线上的法向单位矢量，指向温度升高的方向，“-”号表示热量沿温度降低的方向传递。 2、写出傅里叶定律的文字表达式。 答：在导热中，单位时间内通过给定截面面积的导热量，正比于垂直该截面方向上的温度变化率和截面面积，而热量传递的方向则与温度升高的方向相反。 3、等温面与等温线的特点，不同温度的等温面(线)能相交不？ 答：1) 温度不同的等温面或等温线彼此不能相交； 2) 在连续的温度场中，等温面或等温线不会中断，它们或者是物体中完全封闭的曲面（曲线），或者就终止与物体的边界上； 3）物体的温度场通常用等温面或等温线表示，若每条等温线间的温度间隔相等时，等温线越密反映出该区域导热热流密度的越大。 不同温度的等温面(线)不能相交 4．得出导热微分方程所依据的是什么基本定律? 答：傅里叶定律和能量守恒定律。 5．解释材料的导热系数λ和导温系数α之间的区别和联系? （或热扩散率α的定义及物理意义。） 答：从表达式看，导温系数c ρλ/a =与导热系数成正比关系，但导温系数不但与材料的导热系数有关，还与材料的热容量(或储热能力)也有关；从物理意义看，导热系数表征材料导热能力的强弱，导温系数表征材料传播温度变化的能力的大小，两者都是物性参数。 “稳态导热”例题 例题1：某加热炉炉墙由厚460mm 的硅砖、厚230mm 的轻质粘土砖和厚5mm 的钢板组成，炉墙内表面温度为1600℃，外表面温度为80℃，三层材料的导系数分别为 1.85 W/(m ? K)、0.45 W/(m ? K)和40 W/(m ? K)。已知轻质粘土砖最高使用温度为1300℃，求该炉墙散热的热流密度？并确定轻质粘土砖是否安全？ 解： (1) 3 322114131 4 1λδλδλδλδ++-= -=Φ= ∑=w w i i i w w t t t t A q 2 W/m 200040 /05.045.0/23.085.1/46.080 1600=++-= (2) 1 12 1λδw w t t q -= 1300110285 .146.0200016001112<=?-=-=?λδq t t w w ℃ 因此，轻质粘土砖是安全的。 例题2：某炉壁由厚度=1 δ250mm 的耐火粘土制品层和厚度=2 δ500mm 的红砖层组成。内壁温度=w1t 1000℃，外壁温度=w3 t 50℃。已知耐火粘土制品的导热系数可表示为 t 000233.028.01+=λ，红砖的导热系数近似为 7.02 =λW/(m ? K)。试求稳定运行时，该炉壁单位面积上的散热损失和层间接触界面的温度。 解：由于接触界面温度w2 t 未知，因此无法计算耐火粘土制品层的平均温度，进而无法求得该层的导热系数。现用工程计算中广泛应用的试算法求解。 假设接触界面温度600w2 =t ℃，则耐火粘土制品层的导热系数为 ) K W/(m 466.0 ]2/)6001000[(000233.028.0 ] 2/)[(000233.028.0000233.028.0w2w11?=+?+=++=+=t t t λ 两层炉壁单位面积的散热损失为 2 3 32 2 11w3 w1 W/m 760 7 .010500466.01025050 1000=?+ ?-= +-= --λδλδt t q 校核所假设接触界面的温度w2 t ，得 1 1w2 w1/λδt t q '-= ℃ 593 466 .01025076010003 11w1w2 =?? -=-='-λδq t t 593w2 ='t ℃与假设600 w2 =t ℃相差不大，可认为上述 传热学C 程序源 二维稳态导热的数值计算 2.1物理问题 一矩形区域，其边长L=W=1，假设区域内无内热源，导热系数为常数，三个边温度为T1=0，一个边温度为T2=1，求该矩形区域内的温度分布。 2.2 数学描述 对上述问题的微分方程及其边界条件为：2222T T 0x y ??+=?? x=0，T=T 1=0 x=1，T=T 1=0 y=0，T=T 1=0 y=1，T=T 2=1 该问题的解析解：112121(1)sin n n n sh y T T n L x n T T n L sh W L ππππ∞=??? ?---????=? ?-????? ??? ∑ 2.3数值离散 2.3.1区域离散 区域离散x 方向总节点数为N ，y 方向总节点数为M ，区域内任一节点用I,j 表示。 2.3.2方程的离散 对于图中所有的内部节点方程可写为：2222,,0i j i j t t x y ??????+= ? ??????? 用I,j 节点的二阶中心差分代替上式中的二阶导数，得： +1,,-1,,+1,,-1222+2+0i j i j i j i j i j i j T T T T T T x y --+= 上式整理成迭代形式：()()22 ,1,-1,,1,-12222+2() 2()i j i j i j i j i j y x T T T T T x y x y ++=++++ (i=2,3……,N-1)，(j=2,3……,M-1) 补充四个边界上的第一类边界条件得：1,1j T T = (j=1,2,3……,M) ,1N j T T = (j=1,2,3……,M) ,1i j T T = (i=1,2,3……,N) !*************************************************************** ! 一维非稳态导热通用程序(不变部分) ! This is a general purpose program to solve 1-D diffusion ! problem in the form of: ! ρcdt/dz=1/a(x)d/dx(a(x)Γdt/dx)+s !******************2003.7 revised******************************** !...................Define Variables......................... MODULE VARIABLES INTEGER,PARAMETER::L1=130 REAL,DIMENSION(L1):: X,XF,XM,XP,R,RF,AP REAL,DIMENSION(L1)::AE,AW,CN,T,TA REAL,DIMENSION(L1)::TG,GM,RC INTEGER:: K=1,KM=1,KP=1,OM=1 INTEGER:: JB,JE,KE,KI,KF,KN,KR,KT,LS,MD,M1,M2,NF REAL:: AEC,AI,BE,BI,DF,DS,DT,EP,EX,PW,TU,TM,XE,XI REAL:: A1,A2,T2,TC,SC,SP,RO,TE,DN,LM END !...............................................Main Program........................................... PROGRAM MAIN USE V ARIABLES IMPLICIT NONE INTEGER I OPEN(1,FILE="q.dat") OPEN(2,FILE="temp.dat") NF=1！（求解变量指标） KN=1！（非稳态问题输出局部变量，输出一次加1） TU=0！（当前时间） 150 KT=1！（非线性问题迭代次数） CALL Speci !First to specify the problem CALL Grid!Set up grid points 200 CALL Difsor!Specify the diff-coeff and source term 220 CALL InterOutput!Output intermediate results CALL Coeff!Set up coefficients of discretization equation CALL TDMA!Solve the algebraic equation by TDMA IF(LS.EQ.2.OR.LS.EQ.4) THEN ！（对非线性问题） IF(DF.GT.EP) THEN！（如果最大偏差大于允许值） DO I=1,M1 TA(I)=TA(I)+OM*(T(I)-TA(I))！（采用亚松弛方式将当前解付给上一次迭代值） END DO DF=0！（最大偏差置零） 一维非稳态热传导热源反问题研究 摘要 本文是关于热传导的正反问题的研究，即利用偏微分方程中典型热传导方程 t时刻温度分布与热源位置。 求解含有内热源的金属细杆 本文从解偏微分方程出发，由已知条件最终得出温度分布函数及热源位置函数并建立了两个数学模型。 模型一：利用偏微分方程及初始温度分布函数建立了一段时间后的温度分布与热源强度、位置之间的数学模型，最终解出一段时间后长杆上的温度分布。 模型二：通过一类抛物型偏微分方程模型，解决已知初始温度分布函数、一段时候后的温度分布函数及热源强度的确定热源位置和中间任意时刻的温度分布函数。 u x t,即t时刻的温度根据模型一建立偏微分方程组，用分离变量法求解(,) 分布函数，并通过Matlab中的PDE(偏微分方程)工具箱求解偏微分方程组，且使解可视化。 u x T,结合抛物型方程，运用根据模型二依然建立偏微分方程组，通过测得(,) 离散正则法，确定热源位置，并通过论证说明问题的唯一性和确定性，给出反问题的数值解法。最后再简单介绍差分法解决热传导在非稳态导热问题中的应用。 最后是结论部分，主要总结本文的结果并提出一些尚待进一步研究的问题，以及研究该反问题的应用前景。 相同t不同x的温度变化曲线相同x不同t的温度变化曲线 一维非稳态热传导热源反问题研究 一、问题的提出 在金属细秆的传热过程中，温度差是导致其发生必要条件，有无热源决定传导效率的高低。从一维非稳态传导问题的数学模型和初始条件出发，经过对有内热源问题的进一步分析，在初始温度分布已知的情况下，对分布函数的处理显得很关键。对热源反问题的处理中，我们的问题是如何寻找某种合理的附件条件，通过已知方程来解决方程右端的热源的具体位置并使其具有唯一性。本文利用微分方程并建立了满足温度分布的数学物理模型，从理论上导出了温度分布函数和热源位置的求解，并借助计算机软件画出了温度分布图。 二、问题的分析 对于热传导问题，为了使函数解决起来更容易，对于细秆的初始温度分布() g x我们可以设它在区间[0,L]连续，那么() g x可以展成正弦或余弦级数，对于有内热源的处理，由于细秆边界条件是齐次的，我们采用叠加原理把一根金属细秆的导热问题分解为有热源的具有其次边界条件的稳态导热问题和一个非稳态 其次问题，则原问题的解为 (,)1(,)2() u x t u x t u x =+。 对于源反问题的解决有如下3个问题: 1、反问题的唯一性:附加条件给得是否合理，也就是说，这个附加条件是否可以唯一确定热源的具体位置。 2、反问题的稳定性:反演所得到的热源的具体位置，该热源是否是连续地依赖于测量数据() h t? 3、反问题的数值解法:如何用可行的数值方法反演该热源的具体位置。用离散正则法将温度分布离散化，由已知初始温度分布再利用计算机软件得出热源位置 三、模型假设 1、金属细杆边界与外界无热量交换，即与外界绝缘 “非稳态导热”例题 例题1：一温度为20℃的圆钢，长度为0.3m ，直径为60mm ，在一温度为1250℃的加热炉 内被加热。已知圆钢的导热系数为35 W/(m ?K)，密度为7800kg/m 3，比热容为0.460kJ/(kg ?K)， 加热炉长为6m ，圆钢在其中匀速通过，其表面和炉内烟气间的表面传热系数为100 W/(m 2?K)。现欲将该圆钢加热到850℃，试求该圆钢在加热炉内的通过速度。 解 特征尺寸A V /为 m 0136.0)1060(14.34 13.0)1060(14.33.0)1060(14.3414124133322=???+???????=?+=---d dL L d A V πππ 则毕渥数v Bi 为 05.02 11.01.0039.0350136.0100)/(v =?=<=?==M A V h Bi λ 因此可以采用集总参数法求解。 θθρτ0ln hA cV = 即 s 548.14 1250 850125020ln 100)10460.0(78003=--??=τ 则该圆钢在加热炉内的通过速度为 m /s 0109.014 .5486===τL v 例题2：两块厚度均为30mm 的无限大平板，初始温度为20℃，分别用铜和钢制成。平板 两侧表面的温度突然上升至60℃，计算使两板中心温度均达到56℃时两板所需时间之比。 已知铜和钢的热扩散率分别为610103-?m 2/s 和6 109.12-?m 2/s 。 (125.0==铜 钢钢铜a a ττ) 例题3：无内热源、常物性的二维导热物体在某一瞬时的温度分布为x y t cos 22=。试说明 该导热物体在x =0，y =1处的温度是随时间增加而逐渐升高，还是逐渐降低？ 例题4：一初始温度为20℃的钢板，厚度为10cm ，密度为为7800kg/m 3，比热容为460.5 J/(kg ?K)，导热系数为53.5W/(m ?K)，放置到温度为1200℃的加热炉中加热，钢板与烟气间 的表面传热系数为407 W/(m 2?K)。试求单面加热30min 时该钢板的中心温度以及两面加热 到相同的中心温度需要的时间。 解：(1) 考虑单面加热时，特征尺寸为1m .0cm 10==δ，则毕渥数Bi 为 1.076.05 .531.0407>=?==λδ h Bi 因此不能采用集总参数法求解，可采用图解分析法。钢板中心处无量纲尺寸η为 5.01.01052 =?==-δηx 30min 时的傅里叶数Fo 为 68.21.0)6030()]5.4607800/(5.53[)/(2 22=???= ==δρλδτc a Fo 而毕渥数的倒数1-Bi 为 31.176.011==-Bi 查诺模图可得 93.0 ,21.0m 0m ==θθθθ 则钢板中心的无量纲过余温度0/θθ为 195.093.021.0m 0m f 0f 0=?==--=θθθθθθt t t t 因此钢板中心温度t 为 970)120020(195.01200)(f 00 f =-?+=-+=t t t t θθ℃ (2) 考虑两面加热时，特征尺寸为0.05m cm 2/102/==δ，则毕渥数Bi 为 1.038.05 .5305.0407>=?==λδ h Bi 因此仍不能采用集总参数法求解，可应用图解分析法。此时钢板中心的无量纲过余温度为 第二章部分答案-稳态导热 2-46. 一厚度为7cm 的大平壁，一侧绝热，另一侧暴露于温度为30℃的流体中，其内热源热量为5103?W/m 3。已知该平壁材料的导热系数为18K)W/(m ?，平壁与流体间的对流表面传热系数为450)K W/(m 2?，试确定该平壁中的最高温度位置及其温度值？ 解： (1) 该题为具有内热源的一维平壁稳态导热问题，导热微分方程式为： 022=Φ+λ&dx t d 边界条件为：0=x ，0=dx dt ； δ=x ，∞+Φ==t t t w λ 2&（根据热平衡求得：δΦ=-∞&)(t t h w ) 解方程，并代入边界条件得温度场为： ∞+Φ+-Φ=t h x t δδλ&&)(222 (2) 该平壁中最高温度在0=x 处(即 0=dx dt )： 117.5 30450)107()103()107(182103225252=+???+????=+Φ+Φ=--∞t h t δδλ&&℃ 2-47 核反应堆的辐射防护壁因受γ射线的照射而发热，这相当于防护壁内有ax e -Φ=Φ0&&的内热源，其中0Φ&是X=0的表面上的发热率，a 为已知常数。已知x=0处t=t1,x=δ处t=2t ,试导出该防护壁中温度分布的表达式及最高温度的所在位置。导热系数λ为常数。 解：由题意导热微分方程 0022=Φ+-ax e dx t d &λ 又x=0处t=t1,x=δ处t=2t 积分并结合边界条件可得 λδλλλδ2012020210a t x a e a t t a e t a ax Φ++Φ-Φ+--Φ=--&&&& 令0=dx dt 可得：当()??????-+Φ--=-δδλδa e t t a a x a 1ln 1021时，t 最大。 2-48 核反应堆中一个压力容器的器壁可以按厚为δ的大平壁处理。内表面（x=0处）绝热，外表面维持在恒定温度2t 。γ射线对该容器的加热条件作用可以用一个当量热源Φ&来表示，且ax e -Φ=Φ 0&&，a 为常数，x 是从加热表面起算的距离。在稳态条件下，试： 导出器壁中温度分布的表达式。 确定x=0处的温度。 确定x=δ处的热流密度。 解： 022=Φ+λ&dx t d （1） 边界条件 r=0,0=dx dt (2) 00,t t r r == (3) 三式联立得 ()()20201t x a e e a t ax a +-Φ+-Φ=--δλδλδ x=0时；()202011t a e a t a +Φ+-Φ=-λδλδ 当x=δ时，2t t = 所以 ()110-Φ-=-=-ax e a dx dt q λ 2-49 一半径为1r 的长导线具有均匀内热源Φ&，导热系数为1λ。导线外包有一层 绝缘材料，其外半径为2r ,导热系数为2λ。绝缘材料与周围环境间的表面传热系数为h ，环境温度为∞t 。过程是稳态的，试： 列出导线与绝缘层中温度分布的微分方程及边界条件。 求解导线与绝缘材料中温度分布。 提示：在导线与绝缘材料的界面上，热流密度及温度都是连续的。 解：导线中温度场的控制方程为：0111=???? ??Φ+λ&dr dt r dr d r ； 环形绝缘层中温度场的控制方程为：012=??? ??dr dt r dr d r 。 边界条件：对为有限；时，，110t r t = dr dt dr dt t t r r 2211211,λλ-=-==时，。 第三章 非稳态导热习题 例一腾空置于室内地板上的平板电热器，加在其上的电功率以对流换热和辐射换热的方式全部损失于室内。电热器表面和周围空气的平均对流换热系数为h ，且为常数，室内的空气温度和四壁、天花板及地板的温度相同，均为t f 。电热器假定为均质的固体，密度为ρ，比热为c ，体积为V , 表面积为A ，表面假定为黑体，因其导热系数足够大，内部温度均布。通电时其温度为t 0。试写出该电热器断电后温度随时间变化的数学描述。 [解] 根据题意，电热器内部温度均布，因此可用集中参数分析法处理。 电热器以辐射换热方式散失的热量为： 44r f ()A T T σΦ=- （1） 以对流换热方式的热量为： c f ()hA T T Φ=- （2） 电热器断电后无内热源，根据能量守恒定律，散失的热量应等于电热器能量的减少。若只考虑电热器的热力学能 r c d d T cV ρτ -Φ-Φ= （3） 因此，相应的微分方程式为： 44f f d ()()d T A T T hA T T cV σρτ -+-=- （4） 初始条件为： τ=0， t =t 0 （5） 上述两式即为该电热器断电后温度随时间变化的数学描述。 例 电路中所用的保险丝因其导热系数很大而直径很小可视为温度均布的细长圆柱体，电流的热效应可视为均匀的内热源。如果仅考虑由于对流换热的散热量，保险丝表面和温度为t f 的周围空气之间的平均对流换热系数为h ，且为常数。试求该保险丝通电后温度随时间的变化规律。 [解] 根据题意，保险丝内部温度均布，因此可用集中参数分析法处理。 保险丝表面以对流换热方式散失的热量为： c f ()hA T T Φ=- （1） 保险丝的内热源为： Q 0=IR 2 （2） 式中：I ——保险丝通过的电流，（A ）； R ——保险丝的电阻，Ω。 根据能量守恒，散失的热量与内热源所转变成的热量的和应等于保险丝能量的变化。若只考虑保险丝的热力学能 c 0d d T Q cV ρτ -Φ+= （3） 稳态导热习题解析 1 固体内的一维导热问题 例1 具有均匀内热源强度q v 的无限大平壁处于稳态导热，其厚度为2δ，导热系数λ为常数，两侧壁温各自均布，分别为 t w1和t w2。 试求该平壁内的温度分布表达式。 解: 根据题意，这是导热系数为常数但有均匀内热源强度的稳态导热现象，x 坐标的原点取平壁的中心线。描述该平壁内稳态导热现象的微分方程式为 2v 2d 0d q t x λ += （1） 边界条件： x= -δ: t=t w1 x= δ: t=t w2 （2） 移项后积分该微分方程式两次可得其通解 v 1d d q t x C x λ =-+ 2v 122q t x C x C λ =-++ （3） 代入边界条件 2v w112()()2q t C C δδλ=--+-+ （4） 2v w2 122q t C C δδλ=-++ （5） 式（4）+式（5） 2 w1w2v 22δλ+= +t t q C （6） 式（4）－式（5） w2w1 12t t C δ -= （7） C 1和C 2代入微分方程式的通解式（3）后得到壁内的温度表达式 22v w2w1w2w1(2)222 δλδ-+= -++q t t t t t x x （8） 例2具有均匀内热源q v 的无限大平壁处于稳态导热，其厚度为2δ，导热系数λ为常数， 两侧壁温各自均布且相同，均为t w 。 试求该平壁内的温度分布表达式。 解: 根据题意，所用的导热微分方程式同上。由于两侧壁温相同，是一种对称情况，因此只需求解一半的求解域即可，x 坐标的原点取平壁的中心线。描述该平壁内稳态温度场的微分方程式为 2v 2d 0d q t x λ += （1）一维非稳态导热的数值计算
传热学问答题
稳态导热习题(2020年整理).pdf
第三章非稳态导热分析解法
一维非稳态导热问题的数值解
传热学2章稳态导热总结问答题及答案
稳态导热例题
一维非稳态导热的数值计算
一维非稳态导热通用程序
一维非稳态热传导热源反问题研究
非稳态导热例题
稳态导热问题
非稳态导热习题
稳态导热习题解析