4-19、有一砖墙厚为??0.3m，?＝0.85W/(m.K)，?c?1.05?10J/(m.K)室内温度为
63t1?20℃，h=6W/(m2.K)。起初该墙处于稳定状态，且内表面温度为15℃。后寒潮入侵，
室外温度下降为tf2??10℃，外墙表面传热系数h2?35W/(m.K)。如果认为内墙温度下降0.1℃是可感到外界温度起变化的一个定量判据，问寒潮入侵后多少时间内墙才感知到？
解：采用数值解法得t=7900s。
4-20、一冷柜，起初处于均匀的温度（20℃）。后开启压缩机，冷冻室及冷柜门的内表面温
2度以均匀速度18℃/h下降。柜门尺寸为1.2m?1.2m。保温材料厚8cm，?＝0.02W/(m.K)。冰箱外表面包裹层很薄，热阻可忽略而不计。柜门外受空气自然对流及与环境之间辐射的加热。自然对流可按下式计算： h?1.55??t/H?W/(m2.K)
其中H为门高。表面发射率??0.8。通过柜门的导热可看作为一维问题处理。试计算压缩
机起动后2h内的冷量损失。
4?c?1?104J/?m3?K?解：取保温材料的，用数值计算方法得冷量损失为5.97?10J。
1/4?＝0.81W/(m.K), ??1800kg/m,c?0.88J/?kg.K?。4-21、一砖砌墙壁，厚度为240mm，
3设冬天室外温度为24h内变化如下表所示。室内空气温度ti?15℃且保持不变；外墙表面传热系数为10W/(m.K)，内墙为6W/(m.K)。试用数值方法确定一天之内外墙，内墙及墙壁中心处温度随时间的变化。取???1h。设上述温度工况以24h为周期进行变化。
时刻/h 温度/C
0220：00 1：00 2：00 3：00 4：00 5：00 6：00 7：00 8：00 9：00
10：
00 11：00
-5.9 -6.2 -6.6 -6.7 -6.8 -6.9 -7.2 -7.7 -7.6 -7.0 -4.9 -2.3
时刻/h 温度/C
012：00 13：00 14：00 15：00 16：00 17：00 18：00 19：00 20：00 21：00 22：00 23：00
-1.0 2.4 1.8 1.8 1.6 0.5 -1.6 -2.8 -3.5 -4.3 -4.8 -5.3
解：采用数值解法得出的结果如下表所示。 时刻/h
0
1
2
3
4
5
6
7
8
环境温
度/0C
-5.9 外墙温度/0C -1.70 墙壁中心温度
3.65 /0C
内墙温度/0C
8.99 时刻/h 9 环境温度/0C -7 外墙温度/0C -3.58 墙壁中心温度
2.36 /0C 内墙温度/0C
8.31 -6.2 -6.6 -2.19 -2.44 3.32 3.15 8.82 8.73 10 11 -4.9 -2.3 -3.07 -1.34 2.70 3.87 8.49 9.11 -6.7 -6.8 -2.76 -2.85 2.92 2.87 8.61 8.58
12 13 -1 2.4 0.78 1.87 5.32 6.05 9.87 10.26
-6.9 -7.2 -2.93 -3.01 2.81 2.75 8.55 8.52 14 15 1.8 1.8 4.63 4.15 7.95 7.62 11.26 11.10 -7.7 -7.6
-3.26 -3.67
2.59 2.31
8.43 8.28
16 17
1.6 0.5
4.14 3.97
7.62 7.51
11.10 11.10
时刻/h 18 19 20 21 22 23
环境温度/C 外墙温度/C 墙壁中心温度
0/C 00-1.6 -2.8 -3.5 -4.3 -4.8 -5.3
3.06 1.34 0.36 -0.22 -0.87 -1.29
6.09 5.73 5.05 4.66 4.21 3.93
内墙温度/C
010.71 10.10 9.73 9.53 9.30 9.14
多维稳态导热问题
4-22、如附图所示，一矩形截面的空心电流母线的内外表面分别与温度为tf1,tf2的流体发生对流换热，表面传热系数分别为h1,h2，且各自沿周界是均匀的，电流通过壁内产生均匀热
?。今欲对母线中温度分布进行数值计算，试： 源?（1）划出计算区域
（2）对该区域内的温度分布列出微分方程式及边界条件；
（3）对于图中内角顶外角顶及任一内部节点列出离散方程式（?x??y），设母线的导热系数?为常数。
4-23、一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面的温度分布及每米长度上通过壁面的冷量损失： （1） 内外壁分别维持在10℃及30℃
2t?10h?20W/(m.K)，tf2?30℃，f11（2） 内外壁与流体发生对流换热，且有℃，
h2?4W/(m2.K)。
解：此题应采用计算机求解。如有墙角导热的热点模拟实验设备，则计算参数（如h，?t及网格等）可以取得与实验设备的参数相一致，以把计算结果与实测值作比较。
根据对称性，取1/4区域为计算区域。数值计算解出，对于给定壁温的情形，每米长通道的冷损失为39.84W，对于第三类边界条件为30.97W（取壁面导热系数
??0.53W/?m?K?）
。内外表面为给定壁温时等温线分布如下图所示。第三类边界条件的
结果定性上类似。
4-24、为了提高现代燃气透平的进口燃气温度以提高热效率，在燃气透平的叶片内部开设有冷却通道以使叶片金属材料的温度不超过允许值，为对叶片中的温度分布情况作一估算，把附图a所示的截片形状简化成为附图b所示的情形。已知
T0?1700K,h0?1000W/(m2.K)， Ti?400k,hi?250W/(m2.K)。试计算：（1）截
面中最高温度及其位置；（2）单位长度通道上的热量。
解：根据对称性选择1/4区域为计算区域，采用60?70网格，
??时得单位长度的传热量为987.8W，取壁面
等温线分布如图所示。截面中最高温度发生在左上角，该处
??15W/m?K0温度为1419.9C。
综合分析与分析、论述题
4-25、工业炉的炉墙以往常用红砖和耐火砖组成。由于该两种材料的导热系数较大，散热损失较严重，为了节省能量，近年来国内广泛采用在耐火砖上贴一层硅酸纤维毡，如附图所示。今用以下的非稳态导热简化模型来评价黏贴硅酸纤维毡的收益：
??0s时内壁设炉墙原来处于与环境平衡的状态，
表面突然上升到550℃并保持不变。这一非稳态导热过程一直进行到炉墙外表面的对流，辐射热损失
与通过墙壁的导热量相等为止。在炉墙升温过程中外表面的总表面传热系数由两部分组成，即自然对流引起的部分
1/3?hc?W/?m.K??1.12??tw?c??tf?c?
200及辐射部分
2h?4??T,Tm??Tw?Tf?/2 r0m
tf,Tf为内表面温度，tw,Tw为外表面温度，?1?240mm,?2?240mm,?3?40mm。其中：
为简化计算，设三种材料的导热系数分别为?1?1.6W/(m.K)，?2?0.8W/(m.K)，
?3?0.04W/(m.K)。试计算每平方炉墙每平方面积上由于粘贴了硅酸纤维毡而在炉子升温
过程中节省的能量。
解：采用数值计算方法，详细过程从略。
4-26、空气在附图所示的一长方形截面的送风管道中作充分发展的层流流动，其z方向的动量方程简化为
??2w?2w?dp????x2??y2???dz?0??
dp?而且u?v?0。上式可看成是源项为dz的一常物性导热
方程。试用数值方法求解这一方程并计算f,Re之值。f为阻
dp力系数，Re为特征长度为当量直径De。计算时可任取一个dz值，并按a/b＝0.5及1两种
情形计算。
解：假设壁温为常数，则不同a/b下换热充分发展时的fRe及Nu数的分析解为：
a/b 1 0.5
4-27、一家用烤箱处于稳定运行状态，箱内空气平均温度ti?155℃，气体与内壁间的表面
2h?40W/(m.K)。外壁面与20℃的周围环境间i传热系数
2h?10W/(m.K)。烤箱保温层厚0的表面传热系数
Nu 2.98 3.39
fRe 57 62
30mm,??0.03W/(m.K),保温层两侧的护板用金属制成且很
薄，分析中可不予考虑，然后，突然将烤箱调节器开大，风扇加速，内壁温度突然上升到185℃，设升温过程中烤箱外
h?c?tw?tf壁面与环境间的表面传热系数可用0?1/4计算，
环境温度tf仍保持为20℃，tw为烤箱外壁面温度，c之值与运行时一样。试确定烤箱内壁温度跃升后到达新的稳定状态所需时间。
解：需采用数值方法求解，过程从略。 小论文题目
4-28、一厚为2.54cm的钢管，初始温度为16℃。其后，温度为572℃的液态金属突然流过管内，并经历了10s。液态金属与内壁面间的表面传热系数h=2.84W/(m.K)。钢管可以按平壁处理，其外表面的散热由对流及辐射两条路径，并分别可按
2?hw?W/?m.K??1.2????t?c?1/3及hr03Tm??Tf?Tw?/2，?4??0Tm计算，周围环境温度tf＝
20℃。试用有限差分法确定在液态金属开始流入后的18s时截面上的温度分布。已知钢管的
3??41W/(m.K)，??7530kg/m，c=536J/(kg.K)。
解：在钢管壁厚方向上取27个点，以内壁为坐标原点，沿着壁厚方向为x正方向，数值计算结果如下。 位置/cm
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9