第五章传热学
第五章
复习题
1、试用简明的语言说明热边界层的概念。
答：在壁面附近的一个薄层内，流体温度在壁面的法线方向上发生剧烈变化，而在此薄层之外，流体的温度梯度几乎为零，固体表面附近流体温度发生剧烈变化的这一薄层称为温度边界层或热边界层。 2、与完全的能量方程相比，边界层能量方程最重要的特点是什么？
答：与完全的能量方程相比，它忽略了主流方向温度的次变化率
σ
α2
2
x A
，因此仅适用于边界层内，不适用整个流体。
3、式（5—4）与导热问题的第三类边界条件式（2—17）有什么区别？
答：
=???-
=y
y
t t
h λ（5—4）
)
()(
f w t t h h
t -=??-λ （2—11）
式（5—4）中的h 是未知量，而式（2—17）中的h 是作为已知的边界条件给出，此外（2—17）中的λ为固体导热系数而此式为流体导热系数，式（5—4）将用来导出一个包括h 的无量纲数，只是局部表面传热系数，而整个换热表面的表面系数应该把牛顿冷却公式应用到整个表面而得出。
4、式（5—4）表面，在边界上垂直壁面的热量传递完全依靠导热，那么在对流换热中，流体的流动起什么作用？
答：固体表面所形成的边界层的厚度除了与流体的粘性有关外还与主流区的速度有关，流动速度越大，边界层越薄，因此导热的热阻也就越小，因此起到影响传热大小
5、对流换热问题完整的数字描述应包括什么内容？既然对大多数实际对流传热问题尚无法求得其精确解，那么建立对流换热问题的数字描述有什么意义？
答：对流换热问题完整的数字描述应包括：对流换热微分方程组及定解条件，定解条件包括，（1）初始条件 （2）边界条件 （速度、压力及温度）建立对流换热问题的数字描述目的在于找出影响对流换热中各物理量之间的相互制约关系，每一种关系都必须满足动量，能量和质量守恒关系，避免在研究遗漏某种物理因素。 基本概念与定性分析
5-1 、对于流体外标平板的流动，试用数量级分析的方法，从动
量方程引出边界层厚度的如下变化关系式：
x
x
Re
1~
δ
解：对于流体外标平板的流动，其动量方程为： 2
2
1xy
u v
dx
d y
u v
x
y u
?+-
=??+??ρ
ρ
根据数量级的关系，主流方的数量级为1，y 方线的数量
级为δ
则有
2
21
1
11
1
111δ
ρ
δ
δv
+?
-
=?
+?
从上式可以看出等式左侧的数量级为1级，那么，等式右侧
也是数量级为1级，
为使等式是数量级为1，则v 必须是2
δ量级。
x
δ 从量级看为1δ
级
1
~
1
1~
1
11~
1Re
12
δ
δ
δ
?=
∞v
x u x
量级
两量的数量级相同，所以x δ
与
x
Re
1成比例
5-2、对于油、空气及液态金属，分别有1>>r
P ，1?r P ，1< 试就外标等温平板的层流流动，画出三种流体边界层中速度分布和温度分布的大致图象（要能显示出x δδ与的相对大小）。 解：如下图： 5-3、已知：如图，流体在两平行平板间作层流充分发展对流换热。 求：画出下列三种情形下充分发展区域截面上的流体温度分布曲线：（1）2 1 w w q q =；（2）212w w q q =；（3）01=w q 。 解：如下图形： 5-4、已知：某一电子器件的外壳可以简化成如图所示形状。 c h t >t 。 求：定性地画出空腔截面上空气流动的图像。 解 ： 5-5、已知：输送大电流的导线称为母线，一种母线的截面形状如图所示，内管为导体，其中通以大电流，外管起保护导体的作用。设母线水平走向，内外管间充满空气。 求：分析内管中所产生的热量是怎样散失到周围环境的。并定性地画出截面上空气流动的图像。 解：散热方式：（1）环形空间中的空气自然对流 （2）内环与外环表面间的辐射换热。 5-6、已知：如图，高速飞行部件中广泛采用的钝体是一个轴对称的物体。 求：画出钝体表面上沿x方向的局部表面传热系数的大致图像，并分析滞止点s附近边界层流动的状态。（层流或湍流）。 解：在外掠钝体的对流换热中，滞止点处的换热强度是很高的。该处的流动几乎总处层流状态，对流换热的强烈程度随离开滞止点距离的增加而下降。 5-7．温度为80℃的平板置于来流温度为20℃的气流中．假设平板表面中某点在垂直于壁面方向的温度梯度为40mm ℃，试确定该处 的热流密度． 边界层概念及分析 5-8、已知：介质为25℃的空气、水及14号润滑油，外掠平板边界层的流动由层流转变为湍流的灵界雷诺数5 10 5Re ?=c ，s m u /1=∞。 求：以上三种介质达到c Re 时所需的平板长度。 解：（1）25℃的空气 v =15.536 10-?s m /2 v x u x ∞= Re =5 6 10 510 53.151?=??-x x=7.765m （2）25℃的水 s m v /10 9055.02 6 -?= x=0.45275m （3）14号润滑油 s m v /10 7.3132 6 -?= x=156.85m 5-9、已知：20℃的水以2m/s 的流速平行地流过一块平板，边界层内的流速为三次多项式分布。 求：计算离开平板前缘10cm 及20cm 处的流动边界层厚度及两截面上边界层内流体的质量流量（以垂直于流动方向的单位宽度计）。 解：20℃的水 s m v /10 006 .12 6 -?= s m u /2= （1）x=10cm=0.1m 6 10 00.101.02Re -∞??= = v x u x =19880.72 小于过渡雷诺数x Re . 按（5—22） m u vx 3 6 10 0406.12 1 .010 006.164 .464 .4--∞?=??==δ 设3 ) ( 2 12 3δ δ y y u u y ?-? = ∞ y y y y d y y u d u u u ud u u ud m ])( 2 12 3[ 3 δ δ ρρρ ρδ δ δ δ ?- ? === = ? ? ? ? ∞ ∞ ∞ ∞ = ] 8 43[)](8143 [034 2δδρδδρδ -=-?∞∞u y y u =998.2 ? 2 δ 8 5? =1.298 2 /m kg （2）x=20cm=0.2m 6 10 006.102.02Re -??= x =39761.43 （为尽流） 3 6 10 47.1202 .010 006.164 .464 .4--∞ ?=??==u vx δ m 834 .18 522.9980 =? ?== ? δρδ y x d u m 2 /m kg 5-10、已知：如图，两无限大平板之间的流体，由于上板运动而引起的层流粘性流动称 为库埃流。不计流体中由于粘性而引起的机械能向热能的转换。 求：流体的速度与温度分布。 解：（1）动量方程式简化为 22 =+-dy u d dx dp μ ，y=0, u=0, y=H, ()σ =y u ，σ为上板速度。平行平板间的流动 =dx dp 。积分两次并 代入边界条件得 ()σ ??? ??=H y y u 。 （2）不计及由于粘性而引起机械能向热能的转换，能量方 程为： t k y t x t u c 2?=??? ? ? ???+??ν ρρ,对于所研究的情形，0=ν，0 =??x t ， 因而得0 2 2 =dy t d ，y=0,1w t t =,y=H,2w t t =,由此得() 121w w w t t H y t t -??? ??+=。 5-11、已知：如图，外掠平板的边界层的动量方程式为： 22 y u y u v x u u ??=??+??ν 。 求：沿y 方向作积分（从y=0到δ≥y ）导出边界层的动量积分方程。 解：任一截面做y=0到∞ →y 的积分 dy y u v dy y u v dy x u u ???∞ ∞ ∞ ??= ??+??0 2 2 根据边界层概念y>∞≈u u ,δ 故在该处0 . 0, 02 2 ≈??≈??≈??y u y u x u 则有dy y u v dy y u v dy x u u ???∞ ??= ??+??0 2 2 δ δ (1) 其中dy y u u u v dy y u v ????- =??∞δ δ δ0 由连续行方程可得 dy x u v dy y u dy y v ? ? ???-=??-=??δ δδ δ ; 所以dy x u u dy x u u dy y u v ?????+ ??-=??∞ δ δ δ (2) 又因为 00 2 2 =?? ??? ????-=??y y u v dy y u v δ (3) （1）（2）代入（3）()????∞ -= ??+??-??δ δ δ δ dy u u u dx d dy x u u dy x u dy x u u 故边界层的动量积分方程为()00=∞?? ??? ????=-y y u v dy u u u dx d δ 5-12、已知：Pa 510013.1?、100℃的空气以v=100m/s 的速度流过一块平板，平板温 度为30℃。 求：离开平板前缘3cm 及6cm 处边界层上的法向速度、流动边界 层及热边界层厚度、局 部切应力和局部表面传热系数、平均阻力系数和平均表面传热系数。 解：定性温度65 2 30 100=+= m t ℃ ()K m W ?=/0293.0λ ,695 .0Pr =， s m /10 5.192 6 -?=ν， 3 /045,1m kg =ρ。 （1）cm x 3=处，5 6 10 538.110 5 .1910003.0Re ?=??= =∞ν x u x () s m v /2218.010538.187 .01002 15=??=δ 动量边界层厚度() mm 355.010538.103.064.42 15 =???=-δ mm t 398.0355.0695 .0Pr 3 13 1=?==--δδ () 2 5 2 2 61.810 538.1100 045.1323.0Re 323.0s m kg u x w ?=???= = ∞ ρτ () K m W x h x x ?=???? ==2 5 3 12 16.112695.010 538.103 .00293.0332.0Pr Re 332 .0λ 比拟理论 5-13．来流温度为20℃、速度为4m/s 空气沿着平板流动，在距离前沿点为2m 处的局部切应力为多大？如果平板温度为50℃，该处的对流传热表面传热系数是多少？ 5-14．实验测得一置于水中的平板某点的切应力为1.5Pa ．如果水温与平板温度分别为15℃与60℃，试计算当地的局部热流密度． 5-15．温度为160℃、流速为4m/s 的空气流过温度为30℃的平板．在离开前沿点为2m 处测得局部表面传热系数为149()℃2 m W ．试 计算该处的f x x c j St ,, ,Nu ,Re x 之值． 5-16、已知：将一块尺寸为m m 2.02.0?的薄平板平行地置于由风洞造成的均匀气体流场中。在气流速度s m u /40=∞ 的情况下用测力仪 测得，要使平板维持在气流中需对它施加0.075N 的力。此时气流温度20 =∞ t ℃，平板两平面的温度120 =w t ℃。气体压力为Pa 3 10 013.1?。 求：试据比拟理论确定平板两个表面的对流换热量。 解： Pa m N 9375.0/9375.02 .02.02/075.02 ==?= τ，边界层中空气定性温度为 70℃， 物性： ()694 .0Pr ,/10 02.20,//1009,/029.12 6 3 =?=?==-s m K kg J c m kg p νρ 利用Chilton-Colburn 比拟： 3 /24 2 3 /2Pr 2 ,10 69.52 /40029.19375.02 12 /212 Pr 2p f h f f b c u h c j u c St c j ∞-∞= = ?=?? = = ? == ρρτ ( ) K m W c u c h p f ?=?=?????== ∴---∞2 3 /24 3 /2/1.30276.16.23694 .010*******.11069.5Pr 2 ρ ()()W t t hA w 9.240201202.001.3222 =-???=-=Φ ∞。 这说明Chilton-Colburn 比拟对层流运动也是适用的，即适用于平均值也适用于局部值。 工程应用 5-17．一飞机在10000m 高空飞行，时速为600km/h ．该处温度为-40℃．把机翼当成一块平板，试确定离开机翼前沿点多远的位置上，空气的流动为充分发展的湍流？空气当作干空气处理． 5-18．将一条长度为原型1/4的潜水艇模型放在一闭式风洞中进行阻力试验．潜水艇水下的最大航速为16m/s ，风洞内气体的压力为 Pa 5 106?，模型长3m ，使确定试验时最大的风速应为多少？潜水艇在 水下工作，风洞中的阻力试验结果能否用于水下工作的潜水艇？ 5-19．一火车以25m/s 的速度前进，受到140N 的切应力．它由 1节机车及11节客车车厢组成．将每节车厢都看成是由四个平板所组成，车厢的尺寸为9m（长）m 3? ?（宽）．不计各节车厢间的间 m5.2 隙，车外空气温度为35℃，车厢外表面温度为20℃．试估算该火车所需的制冷负荷． 5-20．在一热处理工程中将一块尺寸为cm 70?平板置于30℃ cm70 的空气气流中，空气流速为 1.2m/s．作用在平板一侧的切应力为0.14N．试估计当该金属板的温度为200℃时平板的散热量． 小论文题目 5-21．夏天，常常将饮料容器置于冰水中来冷却饮料．为了加速冷却，有人提出了这样一个专利（见附图）：将饮料壳体（例如易拉罐）绕其轴线在冰水中做转动．如果能实现饮料瓶或易拉罐绕其轴线的纯转动，试从对流传热基本方程出发，分析这样的方法能否加速饮料的冷却？ 精确解： p=[1,5,10]; x=0:1/19:1; for i=1:1:3 for j=1:1:20 y(i,j)=(exp(p(1,i)*19*x(1,j))-1)/(exp(p(1,i)*19)-1); end plot(x,y(i,:)); hold on ; end 由题对中心差分、一阶迎风、混合格式进行模块编程： 他们之间可以通用，只需更改ae 关于p 的函数即可： 程序如下: (1)中心差分 p=[1,5,10]; for i=1:1:3 ae=1-0.5*p(1,i); x/L (Φ-ΦL )/(Φ0-ΦL ) 精确解图像 aw=p(1,i)+ae; ap=ae+aw; for i=1:1:18 for j=1:1:20 a(i,j)=0; end end for i=1:1:18 j=i; a(i,j)=aw; a(i,j+1)=-ap; a(i,j+2)=ae; end for i=1:1:17 n=i+1; for m=i:-1:1 b(1,1)=a(m,n); a(m,n)=-a(i+1,n)/a(i+1,n)*b(1,1)+a(m,n); a(m,n+1)=-a(i+1,n+1)/a(i+1,n)*b(1,1)+a(m,n+1); a(m,n+2)=-a(i+1,n+2)/a(i+1,n)*b(1,1)+a(m,n+2); end end F(1)=0; F(20)=1; F(19)=(-a(1,20)*F(20)-a(1,1)*F(1))/a(1,19); for i=2:1:18 F(i)=(-a(i,20)*F(20)-a(i,19)*F(19))/a(i,i); end x=0:1/19:1; y(1,:)=F; plot(x,y); hold on end 三、编程题 4.16 设计工程，已知圆的半径r，求圆面积S。 【解答】设圆半径为r，圆面积为S。根据数学知识，已知圆半径r，求圆面积S的公式为：2r Sπ =。 设计步骤如下。 （1）建立应用程序用户界面，如图4-1所示。 （2）设置对象属性： Label1的Caption属性为“已知圆半径r=”； Text1的Text属性为空； Command1的Caption属性为“圆面积为：”； Label2的Caption属性为空； Label2的BorderStyle属性为1-Fixed Single。 各控件的属性设置如图4-2所示。 图4-1 建立用户界面图4-2 设置各控件的属性（3）编写程序代码。 写出“圆面积为：”命令按钮Command1的Click事件代码为： Private Sub Command1_Click( ) Const pi = 3.14 Dim r As Single, S As Single r = V al(Text1.Text) S = pi * r ^ 2 Label2.Caption = S End Sub 运行程序时，在文本框输入圆半径的值，单击“圆面积为：”按钮后，输出结果如图4-3所示。 也可以不用文本框接收输入值，改用InputBox函数接收圆的半径r，求圆面积S，代码如下。 图4-3 程序运行结果 Private Sub Form_Load( ) Show Const pi = 3.1415926 Dim r As Single, S As Single r = V al(InputBox("输入半径：", "计算圆面积", "10")) FontSize = 18 S = pi * r ^ 2 Print "圆面积："; S End Sub 程序运行时，首先显示如图4-4所示的对话框，在该对话框的文本框中输入数字，按Enter 键或单击“确定”按钮后，才能显示窗体。 图4-4 输入对话框 用InputBox 函数输入文本虽然很方便，但是由于输入框弹出后将暂停程序的运行，直到用户响应，因此输入框不符合VB 自由环境的精神。输入框适合于像要求用户输入口令等这样不常见的输入方式。还可以用更好的用户输入方式，如文本框、选项按钮等。 4.17 已知平面坐标系中两点的坐标，求两点间的距离。 【解答】 由数学知识可知，已知两点坐标（x A , y A ）、（x B , y B ），求两点间距离的计算公式为 2 A B 2 A B )()(y y x x s -+-= 建立用户界面如图4-5所示。在该界面中用TextBox 控件输入数据，用Label 控件输出数据。为了形象地表示两点之间的距离，可用Picture 控件插入一幅图，该图用画图软件绘制。 命令按钮Command1的Click 事件代码为： Private Sub Command1_Click( ) Dim xa As Single, xb As Single Dim ya As Single, yb As Single Dim s As Single xa = Val(Text1.Text) ya = V al(Text2.Text) xb = V al(Text3.Text) yb = V al(Text4.Text) s = Sqr((xb - xa) ^ 2 + (yb - ya) ^ 2) Label6.Caption = s End Sub 程序运行结果如图4-6所示。 传热学试题 第一章概论 一、名词解释 1．热流量：单位时间所传递的热量 2．热流密度：单位传热面上的热流量 3．导热：当物体有温度差或两个不同温度的物体接触时，在物体各部分之间不发生相对位移的情况下，物质微粒(分子、原子或自由电子)的热运动传递了热量，这种现象被称为热传导，简称导热。 4．对流传热：流体流过固体壁时的热传递过程，就是热对流和导热联合用的热量传递过程，称为表面对流传热，简称对流传热。 5．辐射传热：物体不断向周围空间发出热辐射能，并被周围物体吸收。同时，物体也不断接收周围物体辐射给它的热能。这样，物体发出和接收过程的综合结果产生了物体间通过热辐射而进行的热量传递，称为表面辐射传热，简称辐射传热。 6．总传热过程：热量从温度较高的流体经过固体壁传递给另一侧温度较低流体的过程，称为总传热过程，简称传热过程。 7．对流传热系数：单位时间单位传热面当流体温度与壁面温度差为1K是的对流传热量，单位为W／(m2·K)。对流传热系数表示对流传热能力的大小。 8．辐射传热系数：单位时间单位传热面当流体温度与壁面温度差为1K是的辐射传热量，单位为W／(m2·K)。辐射传热系数表示辐射传热能力的大小。 9．复合传热系数：单位时间单位传热面当流体温度与壁面温度差为1K是的复合传热量，单位为W／(m2·K)。复合传热系数表示复合传热能力的大小。 10．总传热系数：总传热过程中热量传递能力的大小。数值上表示传热温差为1K时，单位传热面积在单位时间的传热量。 二、填空题 1．热量传递的三种基本方式为、、。 （热传导、热对流、热辐射） 2．热流量是指，单位是。热流密度是指，单位是。 （单位时间所传递的热量，W，单位传热面上的热流量，W/m2） 3．总传热过程是指，它的强烈程度用来衡量。 (热量从温度较高的流体经过固体壁传递给另一侧温度较低流体的过程，总传热系数) 4．总传热系数是指，单位是。 （传热温差为1K时，单位传热面积在单位时间的传热量，W／(m2·K)） 5．导热系数的单位是；对流传热系数的单位是；传热系数的单位是。 （W／(m·K)，W／(m2·K)，W／(m2·K)） 第五章 复习题 1、试用简明的语言说明热边界层的概念。 答：在壁面附近的一个薄层内，流体温度在壁面的法线方向上发生剧烈变化，而在此薄层之外，流体的温度梯度几乎为零，固体表面附近流体温度发生剧烈变化的这一薄层称为温度边界层或热边界层。 2、与完全的能量方程相比，边界层能量方程最重要的特点是什么 答：与完全的能量方程相比，它忽略了主流方向温度的次变化率，因此仅适用于边界层内，不适用整个流体。 3、式（5—4）与导热问题的第三类边界条件式（2—17）有什么区别 答：（5—4）（2—11） 式（5—4）中的h是未知量，而式（2—17）中的h是作为已知的边界条件给出，此外（2—17）中的为固体导热系数而此式为流体导热系数，式（5—4）将用来导出一个包括h的无量纲数，只是局部表面传热系数，而整个换热表面的表面系数应该把牛顿冷却公式应用到整个表面而得出。 4、式（5—4）表面，在边界上垂直壁面的热量传递完全依靠导热，那么在对流换热中，流体的流动起什么作用 答：固体表面所形成的边界层的厚度除了与流体的粘性有关外还与主流区的速度有关，流动速度越大，边界层越薄，因此导热的热阻也就越小，因此起到影响传热大小 5、对流换热问题完整的数字描述应包括什么内容既然对大多数实际对流传热问题尚无法求得其精确解，那么建立对流换热问题的数字描述有什么意义 答：对流换热问题完整的数字描述应包括：对流换热微分方程组及定解条件，定解条件包括，（1）初始条件（2）边界条件（速度、压力及温度）建立对流换热问题的数字描述目的在于找出影响对流换热中各物理量之间的相互制约关系，每一种关系都必须满足动量，能量和质量守恒关系，避免在研究遗漏某种物理因素。 基本概念与定性分析 5-1 、对于流体外标平板的流动，试用数量级分析的方法，从动量方程引出边界层厚度的如下变化关系式： 解：对于流体外标平板的流动，其动量方程为： 根据数量级的关系，主流方的数量级为1，y方线的数量级为 则有 从上式可以看出等式左侧的数量级为1级，那么，等式右侧也是数量级为1级，为使等式是数量级为1，则必须是量级。 取步长δx=0.02。已知x=0,Φ=0;x=1,Φ=1.令k=ρｕ/Γ计算结果图表： 程序及数据结果： 追赶法： #include a[i]=2+0.02*k; b[i]=4; c[i]=2-0.02*k; f[i]=0; } tdma(a,b,c,f,x); for(i=0;i 传热学习题集第一章 思考题 1．试用简练的语言说明导热、对流换热及辐射换热三种热传递方式之间的联系和区别。 答：导热和对流的区别在于：物体内部依靠微观粒子的热运动而产生的热量传递现象，称为导热；对流则是流体各部分之间发生宏观相对位移及冷热流体的相互掺混。联系是：在发生对流换热的同时必然伴生有导热。 导热、对流这两种热量传递方式，只有在物质存在的条件下才能实现，而辐射可以在真空中传播，辐射换热时不仅有能 量的转移还伴有能量形式的转换。 2．以热流密度表示的傅立叶定律、牛顿冷却公式及斯忒藩－玻耳兹曼定律是应当熟记的 传热学公式。试写出这三个公式并说明其中每一个符号及其意义。 答：① 傅立叶定律： dx dt q λ-=，其中，q －热流密度；λ－导热系数；dx dt －沿x 方向的温度变化率，“－”表示热量传递的方向是沿着温度降低的方向。 ② 牛顿冷却公式：)(f w t t h q -=，其中，q －热流密度；h －表面传热系数；w t －固体表面温度；f t －流体的温度。 ③ 斯忒藩－玻耳兹曼定律：4T q σ=，其中，q －热流密度；σ－斯忒藩－玻耳 兹曼常数；T －辐射物体的热力学温度。 3．导热系数、表面传热系数及传热系数的单位各是什么？哪些是物性参数，哪些与过程 有关？ 答：① 导热系数的单位是：W/(m.K)；② 表面传热系数的单位是：W/(m 2.K)；③ 传热系数的单位是：W/(m 2.K)。这三个参数中，只有导热系数是物性参数，其它均与过程有关。 4．当热量从壁面一侧的流体穿过壁面传给另一侧的流体时，冷、热流体之间的换热量可 以通过其中任何一个环节来计算（过程是稳态的），但本章中又引入了传热方程式，并说它是“换热器热工计算的基本公式”。试分析引入传热方程式的工程实用意义。答：因为在许多工业换热设备中，进行热量交换的冷、热流体也常处于固体壁面的两侧，是工程技术中经常遇到的一种典型热量传递过程。 5．用铝制的水壶烧开水时，尽管炉火很旺，但水壶仍然安然无恙。而一旦壶内的水烧干 后，水壶很快就烧坏。试从传热学的观点分析这一现象。 答：当壶内有水时，可以对壶底进行很好的冷却（水对壶底的对流换热系数大），壶底的热量被很快传走而不至于温度升得很高；当没有水时，和壶底发生对流换热的是气体，因为气体发生对流换热的表面换热系数小，壶底的热量不能很快被传走，故此壶底升温很快，容易被烧坏。 6．用一只手握住盛有热水的杯子，另一只手用筷子快速搅拌热水，握杯子的手会显著地 感到热。试分析其原因。 答：当没有搅拌时，杯内的水的流速几乎为零，杯内的水和杯壁之间为自然对流换热，自热对流换热的表面传热系数小，当快速搅拌时，杯内的水和杯壁之间为强制对流换热，表面传热系数大，热水有更多的热量被传递到杯壁的外侧，因此会显著地感觉到热。 7．什么是串联热阻叠加原则，它在什么前提下成立？以固体中的导热为例，试讨论有哪 些情况可能使热量传递方向上不同截面的热流量不相等。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各串联环节热阻的和。例如：三块无限大平板叠加构成的平壁。例如通过圆筒壁，对于各个传热环节的传热面积不相等，可能造成热量传递方向上不同截面的热流量不相等。 8.有两个外形相同的保温杯A 与B ，注入同样温度、同样体积的热水后不久，A 杯的外表面就可以感觉到热，而B 杯的外表面则感觉不到温度的变化，试问哪个保温杯的质量较好？ 4-1 解：采用区域离散方法A 时；网格划分如右图。内点采用中心差分 23278.87769.9 T T T === 22d T T=0dx - 有 i+1i 12 2+T 0i i T T T x ---=? 将2点,3点带入 32122 2+T 0T T T x --=? 即321 209T T -+= 432322+T 0T T T x --=?4321322+T 0T T T x --=? 即4 321 209 T T T -+-= 边界点4 （1）一阶截差 由x=1 1dT dx =，得 431 3 T T -= （2）二阶截差 11B M M q x x x T T S δδλλ -=++ 所以 434111. 1. 36311 T T T =++ 即 431 22293 T T -= 采用区域离散方法B 22d T T=0dx - 由控制容积法 0w e dT dT T x dT dT ????--?= ? ????? 所以代入2点4点有 322121011336 T T T T T ----= 即 239 028T T -= 544431011363 T T T T T ----= 即 34599 02828T T T -+= 对3点采用中心差分有 432 32 2+T 013T T T --=?? ??? 即 23499 01919 T T T -+= 对于点5 由x=1 1dT dx =，得 541 6 T T -= （1）精确解求左端点的热流密度 由 ()2 1 x x e T e e e -= -+ 所以有 ()22 20.64806911x x x x dT e e q e e dx e e λ -====- +=-=++ （2）由A 的一阶截差公式 21 0.247730.743113 x T T dT q dx λ =-=-= =?= （3）由B 的一阶截差公式 0 0.21640 0.649213 x dT q dx λ =-=-= = （4）由区域离散方法B 中的一阶截差公式： 210.108460.6504()B B T T dT dx x δ-?? ==?= ? ?? 通过对上述计算结果进行比较可得：区域离散B 有控制容积平衡法建立的离散方程与区域离散方程A 中具有二阶精度的格式精确度相当！ 4-3 解：将平板沿厚度方向3等分，如图 传热学习题_建工版V 0-14 一大平板，高3m ，宽2m ，厚0.2m ，导热系数为45W/(m.K), 两侧表面温度分别为w1t 150C =?及w1t 285C =? ，试求热流密度计热流量。 解：根据付立叶定律热流密度为： 2 w2w121t t 285150q gradt=-4530375(w/m )x x 0.2λλ??--??=-=-=- ? ?-???? 负号表示传热方向与x 轴的方向相反。 通过整个导热面的热流量为： q A 30375(32)182250(W)Φ=?=-??= 0-15 空气在一根内经50mm ，长2.5米的管子内流动并被加热，已知空气的平均温度为85℃，管壁对空气的h=73(W/m 2.k)，热流密度q=5110w/ m 2, 是确定管壁温度及热流量?。 解：热流量 qA=q(dl)=5110(3.140.05 2.5) =2005.675(W) πΦ=?? 又根据牛顿冷却公式 w f hA t=h A(t t )qA Φ=??-= 管内壁温度为： w f q 5110t t 85155(C)h 73 =+ =+=? 1-1．按20℃时，铜、碳钢（1.5%C ）、铝和黄铜导热系数的大小，排列它们的顺序；隔热保温材料导热系数的数值最大为多少？列举膨胀珍珠岩散料、矿渣棉和软泡沫塑料导热系数的数值。 解： (1)由附录7可知，在温度为20℃的情况下， λ 铜 =398 W/(m ·K)，λ 碳钢 ＝36W/(m ·K)， λ 铝 ＝237W/(m ·K)，λ 黄铜 ＝109W/(m ·K). 所以,按导热系数大小排列为: λ 铜 >λ 铝 >λ 黄铜 >λ钢 (2) 隔热保温材料定义为导热系数最大不超过0.12 W/(m ·K). (3) 由附录8得知，当材料的平均温度为20℃时的导热系数为: 膨胀珍珠岩散料:λ=0.0424+0.000137t W/(m ·K) =0.0424+0.000137×20=0.04514 W/(m ·K); 矿渣棉: λ=0.0674+0.000215t W/(m ·K) =0.0674+0.000215×20=0.0717 W/(m ·K); 第五章 复习题 1、试用简明的语言说明热边界层的概念。 答：在壁面附近的一个薄层内，流体温度在壁面的法线方向上发生剧烈变化，而在此薄层之外，流体的温度梯度几乎为零，固体表面附近流体温度发生剧烈变化的这一薄层称为温度边界层或热边界层。 2、与完全的能量方程相比，边界层能量方程最重要的特点是什么？ 答：与完全的能量方程相比，它忽略了主流方向温度的次变化率 σ α2 2 x A ，因此仅适用于边界层内，不适用整个流体。 3、式（5—4）与导热问题的第三类边界条件式（2—17）有什么区别？ 答： =???- =y y t t h λ（5—4） ) ()( f w t t h h t -=??-λ （2—11） 式（5—4）中的h 是未知量，而式（2—17）中的h 是作为已知的边界条件给出，此外（2—17）中的λ为固体导热系数而此式为流体导热系数，式（5—4）将用来导出一个包括h 的无量纲数，只是局部表面传热系数，而整个换热表面的表面系数应该把牛顿冷却公式应用到整个表面而得出。 4、式（5—4）表面，在边界上垂直壁面的热量传递完全依靠导热，那么在对流换热中，流体的流动起什么作用？ 答：固体表面所形成的边界层的厚度除了与流体的粘性有关外还与主流区的速度有关，流动速度越大，边界层越薄，因此导热的热阻也就越小，因此起到影响传热大小 5、对流换热问题完整的数字描述应包括什么内容？既然对大多数实际对流传热问题尚无法求得其精确解，那么建立对流换热问题的数字描述有什么意义？ 答：对流换热问题完整的数字描述应包括：对流换热微分方程组及定解条件，定解条件包括，（1）初始条件 （2）边界条件 （速度、压力及温度）建立对流换热问题的数字描述目的在于找出影响对流换热中各物理量之间的相互制约关系，每一种关系都必须满足动量，能量和质量守恒关系，避免在研究遗漏某种物理因素。 基本概念与定性分析 5-1 、对于流体外标平板的流动，试用数量级分析的方法，从动 量方程引出边界层厚度的如下变化关系式： x x Re 1~ δ 解：对于流体外标平板的流动，其动量方程为： 2 2 1xy u v dx d y u v x y u ?+- =??+??ρ ρ 根据数量级的关系，主流方的数量级为1，y 方线的数量 级为δ 则有 2 21 1 11 1 111δ ρ δ δv +? - =? +? 从上式可以看出等式左侧的数量级为1级，那么，等式右侧 也是数量级为1级， 为使等式是数量级为1，则v 必须是2 δ量级。 x δ 从量级看为1δ 级 1 ~ 1 1~ 1 11~ 1Re 12 δ δ δ ?= ∞v x u x 量级 两量的数量级相同，所以x δ 与 x Re 1成比例 5-2、对于油、空气及液态金属，分别有1>>r P ，1?r P ，1< 传热学习题答案 第一章 蓝色字体为注释部分 1-4、对于附图中所示的两种水平夹层，试分析冷、热表面间的热量交换方式有什么不同？如果要通过实验来测定夹层中流体的导热系数，应采用哪种布置？ 答：图（a）的热量交换方式为导热（热传导），图（b）的热量交换方式为 导热（热传导）及自然对流。应采用图(a)的方式来测定流体的导热系数。 解释：因为图(a)热面在上，由于密度不同，热流体朝上，冷流体朝下，冷 热流体通过直接接触来交换热量，即导热；而图（b）热面在下，热流体密度小，朝上运动，与冷流体进行自然对流，当然也有导热。 因为图（a）中只有导热，测定的传热系数即为导热系数；而图（b）有导热和自然对流方式，测定的传热系数为复合传热系数。 1-6、一宇宙飞船的外形如附图所示，其中外遮光罩是凸出于飞船船体之外的一个光学窗口，其表面的温度状态直接影响飞船的光学遥感器。船体表面各部 分的表面温度与遮光罩的表面温度不同。试分析：飞船在太空中飞行时与外遮光 罩表面发生热交换的对象可能有哪些？换热方式是什么？ 答：可能与外遮光罩表面发生热交换的对象有两个：一个是外遮光罩表面与 外太空进行辐射换热，另一个是外遮光罩表面与船体表面进行辐射换热。 解释：在太空中，只有可能发生热辐射，只要温度大于0K，两个物体就会发生辐射换热。 1-9、一砖墙的表面积为12m2,厚260mm，平均导热系数为1.5W/(m.K)，设面向室内的表面温度为25℃，外表面温度为-5℃，试确定此砖墙向外界散失的热 = 8.5 =0.7?5.67?10-8?(2504-0) 量。 Φ=A λ(t-t)δw1w2 解：=12?1.5 ? (25-(-5)) 0.26 =2076.92W 此砖墙向外界散失的热量为2076.92W。 1-12、在一次测定空气横向流过单根圆管的对流换热实验中，得到下列数据：管壁平均温度t w=69℃，空气温度t f=20℃，管子外径d=14mm，加热段长80mm，输入加热段的功率为8.5W。如果全部热量通过对流传热传给空气，试问此时的对流传热表面传热系数多大？ 解：此题为对流传热问题，换热面积为圆管外侧表面积，公式为： Φ=hA(t-t)=h?πdl?(t-t) w f w f h=Φ πdl?(t-t) w f ∴ 3.14?0.014?0.08?(69-20) =49.3325W (m2?K) 此时的对流传热表面传热系数49.3325W/(m2.K) 1-18、宇宙空间可近似地看成为0K的真空空间。一航天器在太空中飞行，其外表面平均温度为250K，表面发射率为0.7，试计算航天器单位表面上的换热量。 解：此题为辐射换热问题，公式为： q=εσ(T4-T4) 12 =155.04W m2 航天器单位表面上的换热量为155.04W/m2。 复习题 1、试用简明的语言说明热边界层的概念。 答：在壁面附近的一个薄层内，流体温度在壁面的法线方向上发生剧烈变化，而在此 薄层之外，流体的温度梯度几乎为零， 固体表面附近流体温度发生剧烈变化的这一薄层称为 温度边界层或热边界层。 2、与完全的能量方程相比，边界层能量方程最重要的特点是什么? 答：与完全的能量方程相比，它忽略了主流方向温度的次变化率 适用于边界层内，不适用整个流体。 3、式（5—4）与导热问题的第三类边界条件式（ 2 —17）有什么区另 一个包括h 的无量纲数，只是局部表面传热系数，而整个换热表面的表面系数应该把 牛顿冷却公式应用到整个表面而得出。 4、式（5—4）表面，在边界上垂直壁面的热量传递完全依靠导热，那么在对流换热中，流 体的流动起什么作用? 答：固体表面所形成的边界层的厚度除了与流体的粘性有关外还与主流区的速度有关, 流动速度越大，边界层越薄，因此导热的热阻也就越小，因此起到影响传热大小 5、对流换热问题完整的数字描述应包括什么内容？既然对大多数实际对流传热问题尚无法 求得其精确解，那么建立对流换热问题的数字描述有什么意义? 答：对流换热问题完整的数字描述应包括：对流换热微分方程组及定解条件，定解条件 包括，（1）初始条件 （2 ）边界条件 （速度、压力及温度）建立对流换热问题的数字描述 目的在于找出影响对流换热中各物理量之间的相互制约关系，每一种关系都必须满足动量, 能量和质量守恒关系，避免在研究遗漏某种物理因素。 基本概念与定性分析 5-1、对于流体外标平板的流动， 试用数量级分析的方法， 从动量方程引出边界层厚度 解：对于流体外标平板的流动，其动量方程为: 第五章 2 / 2 A / X ，因此仅 h 答： (5— 4) (丄)h(t w t f ) h (2—11) 式（5—4）中的 h 是未知量，而式（2 —17）中的h 是作为已知的边界条件给出, 此外（2 —17）中的 为固体导热系数而此式为流体导热系数，式（ 5— 4）将用来导出 的如下变化关系式: x 5-2 解：根据课本p158式（5—1a ）得一维稳态无源项的对流－扩散方程如下所示： 2 2x x u ??Γ =??φ φρ （取常物性） 边界条件如下： L L x x φφφφ====,; ,00 由（5—2）得方程的精确解为： 1 1)/(00--=--?Pe L x Pe L e e φφφφ Γ=/uL Pe ρ 将L 分成15等份，有：?=P Pe 15 对于中心差分、一阶迎风、混合格式和QUICK 格式分别分析如下： 1） (CD)中心差分 节点离散方程： 2 )5.01()5.01(1 1-?+?++-=i i i P P φφφ 10,2 =i 2） 一阶迎风 节点离散方程： ? -?++++=P P i i i 2)1(1 1φφφ 10,2 =i 3） 混合格式 当1=?P 时，节点离散方程：2 )5.01()5.01(1 1-?+?++-= i i i P P φφφ ，10,2 =i 当10,5=?P 时，节点离散方程： 1-=i i φφ ， 10,2 =i 4） QUICK 格式，节点离散方程： ??? ???--++++++= +-?? -??+?)336(8122121 1111i i i i i i P P P P P φφφφφφ， 2=i ?? ????---++++++= +--? ? -??+?)35(8122121 12111i i i i i i i P P P P P φφφφφφφ， 2≠i 用matlab 编程如下：（本程序在x/L=0-1范围内取16个节点进行离散计算，假设y(1)= 0φ=0,y(16)=L φ=1,程序中Pa 为?P ，x 为题中所提的x/L 。由于本程序假设 y(1)=0φ=0,y(16)=L φ=1，所以 y y y y y y L =--=--=--0 10 )1()16()1(00φφφφ） Pa=input('请输入Pa=') x=0:1/15:1 Pe=15*Pa; y=(exp(Pe*x)-1)/(exp(Pe)-1) plot(x,y,'-*k') %精确解 hold on y(1)=0,y(16)=1; for i=2:15 y(i)=((1+0.5*Pa)*y(i-1)+(1-0.5*Pa)*y(i+1))/2; end plot(x,y(1:16),'-or') %中心差分 hold on for i=2:15 y(i)=((1+Pa)*y(i-1)+y(i+1))/(2+Pa); end plot(x,y(1:16),'-.>g') %一阶迎风 hold on for i=2:15 if Pa==1 y(i)=((1+0.5*Pa)*y(i-1)+(1-0.5*Pa)*y(i+1))/2; else y(i)=y(i-1) end end plot(x,y(1:16),'-+y') %混合格式 hold on for i=2:15 if i==2 y(i)=y(i+1)/(2+Pa)+(1+Pa)*y(i-1)/(2+Pa)+(Pa/(2+Pa))*(6*y(i)-3*y(i-1)-3*y(i+1))/8 else y(i)=y(i+1)/(2+Pa)+(1+Pa)*y(i-1)/(2+Pa)+(Pa/(2+Pa))*(5*y(i)-y(i-1)-y(i-2)-3*y(i+1))/8 end end plot(x, y(1:16),'- 绪论 1.冰雹落体后溶化所需热量主要是由以下途径得到： Q λ——与地面的导热量 f Q——与空 气的对流换热热量 注：若直接暴露于阳光下可考虑辐射换热，否则可忽略不计。6．夏季：在维持20℃的室内，人体通过与空气的对流换热失去热量，但同时又与外界和内墙面通过辐射换热得到热量，最终的 总失热量减少。（T T? 外内 ） 冬季：在与夏季相似的条件下，一方面人体通过对流换热失去部分热量，另一方面又与外界和内墙通过辐射换热失去部分 热量，最终的总失热量增加。（T T? 外内 ）。挂上窗帘布阻断了与外界的辐射换热，减少了人体的失热量。 7．热对流不等于对流换热，对流换热 = 热对流 + 热传导热对流为基本传热方式，对流换热为非基本传热方式 8．门窗、墙壁、楼板等等。以热传导和热对流的方式。 9．因内、外两间为真空，故其间无导热和对流传热，热量仅能通过胆壁传到外界，但夹层两侧均镀锌，其间的系统辐射系数 降低，故能较长时间地保持热水的温度。 当真空被破坏掉后，1、2两侧将存在对流换热，使其保温性 能变得很差。 10．t R R A λλ = ? 1t R R A λ λ = = 221 8.331012 m --=? 11．q t λσ =? const λ=→直线 const λ≠ 而为λλ=（t ） 时→曲线 12. i R α 1 R λ 3 R λ 0 R α 1 f t ??→ q 首先通过对流换热使炉子内壁温度升高，炉子内壁通过热传导，使内壁温度生高，内壁与空气夹层通过对流换热继续传递热量，空气夹层与外壁间再通过热传导，这样使热量通过空气夹层。（空气夹层的厚度对壁炉的保温性能有影响，影响a α的大小。） 13．已知：360mm σ=、0.61()W m K λ=? 1 18f t =℃ 2187() W h m K =? 2 10f t =-℃ 22124() W h m K =? 墙高2.8m ，宽3m 求：q 、1 w t 、2 w t 、φ 解：12 11t q h h σλ?= ++＝ 18(10) 45.9210.361 870.61124 --=++2W m 传热学习题集 第一章 思考题 1． 试用简练的语言说明导热、对流换热及辐射换热三种热传递方式之间的联系和区别。 答：导热和对流的区别在于：物体内部依靠微观粒子的热运动而产生的热量传递现象，称为导热；对流则是流体各部分之间发生宏观相对位移及冷热流体的相互掺混。联系是：在发生对流换热的同时必然伴生有导热。 导热、对流这两种热量传递方式，只有在物质存在的条件下才能实现，而辐射可以在真空中传播，辐射换热时不仅有能 量的转移还伴有能量形式的转换。 2． 以热流密度表示的傅立叶定律、牛顿冷却公式及斯忒藩－玻耳兹曼定律是应当熟记的传 热学公式。试写出这三个公式并说明其中每一个符号及其意义。 答：① 傅立叶定律： dx dt q λ -=，其中，q －热流密度；λ－导热系数；dx dt －沿x 方向的温度变化率，“－”表示热量传递的方向是沿着温度降低的方向。 ② 牛顿冷却公式：) (f w t t h q -=，其中，q －热流密度；h －表面传热系数；w t －固体表面温度； f t －流体的温度。 ③ 斯忒藩－玻耳兹曼定律：4 T q σ=，其中，q －热流密度；σ－斯忒藩－玻耳 兹曼常数；T －辐射物体的热力学温度。 3． 导热系数、表面传热系数及传热系数的单位各是什么？哪些是物性参数，哪些与过程有 关？ 答：① 导热系数的单位是：W/(m.K)；② 表面传热系数的单位是：W/(m 2.K)；③ 传热系数的单位是：W/(m 2.K)。这三个参数中，只有导热系数是物性参数，其它均与过程有关。 4． 当热量从壁面一侧的流体穿过壁面传给另一侧的流体时，冷、热流体之间的换热量可以 通过其中任何一个环节来计算（过程是稳态的），但本章中又引入了传热方程式，并说它是“换热器热工计算的基本公式”。试分析引入传热方程式的工程实用意义。 答：因为在许多工业换热设备中，进行热量交换的冷、热流体也常处于固体壁面的两侧，是工程技术中经常遇到的一种典型热量传递过程。 5． 用铝制的水壶烧开水时，尽管炉火很旺，但水壶仍然安然无恙。而一旦壶内的水烧干后， 水壶很快就烧坏。试从传热学的观点分析这一现象。 答：当壶内有水时，可以对壶底进行很好的冷却（水对壶底的对流换热系数大），壶底的热量被很快传走而不至于温度升得很高；当没有水时，和壶底发生对流换热的是气体，因为气体发生对流换热的表面换热系数小，壶底的热量不能很快被传走，故此壶底升温很快，容易被烧坏。 6． 用一只手握住盛有热水的杯子，另一只手用筷子快速搅拌热水，握杯子的手会显著地感 到热。试分析其原因。 答：当没有搅拌时，杯内的水的流速几乎为零，杯内的水和杯壁之间为自然对流换热，自热对流换热的表面传热系数小，当快速搅拌时，杯内的水和杯壁之间为强制对流换热，表面传热系数大，热水有更多的热量被传递到杯壁的外侧，因此会显著地感觉到热。 7． 什么是串联热阻叠加原则，它在什么前提下成立？以固体中的导热为例，试讨论有哪些 情况可能使热量传递方向上不同截面的热流量不相等。 答：在一个串联的热量传递过程中，如果通过每个环节的热流量都相同，则各串联环节的总热阻等于各串联环节热阻的和。例如：三块无限大平板叠加构成的平壁。例如通过圆筒壁，对于各个传热环节的传热面积不相等，可能造成热量传递方向上不同截面的热流量不相等。 8.有两个外形相同的保温杯A 与B ，注入同样温度、同样体积的热水后不久，A 杯的外表面就可以感觉到热，而B 杯的外表面则感觉不到温度的变化，试问哪个保温杯的质量较好？ 第七章 凝结与沸腾换热 1．凝液量：m=(kg/s) 2.水平放置时，凝水量m=(kg/s) 3.壁温t w =1000 , h=12029 w/(m 2·k) 4. 5.此时管下端液膜内已出现紊流。 H=6730 w/(m 2·k) 6.竖壁高 h= mm 7.单管与管束平均表面传热系数之比：管束 单h h = 8．凝结水量 m=? (kg/s) 9.考虑过冷度时，m=?(kg/s) 相差： %39.0%10014 .512 .514.5=?- 10．管长 m L 1= ,管长减少量31 5 .115.1= - 11．凝结表面传热系数 h= w/(m 2·k) 凝液量：m=?(kg/s) 12. 管长能缩短 13．用于水时， h= w/(m 2·k) 与11题相比换热系数倍率 63.72 .7001 .5341= 15．氟利昂 12： φ=42143（W ） 氟利昂 22： φ=50810（W ） 差异：% 16．用电加热时，加热方式是控制表面的热流密度。而采用蒸汽加热则是壁面温度可控的情形。由大容器饱和沸腾曲线可知，当加热功率q 稍超过max q 值时，工况将沿max q 虚线跳至稳定膜态沸腾线，使壁面温度飞升，导致设备烧坏。总之，电加热等依靠控制热流来改变工况的设备，一旦热流密度超过峰值，工况超过热流密度峰值后，沸腾温差将剧烈上升到1000℃左右，壁温也急剧升高，发生器壁烧毁现象。 采用蒸气加热时，工况点沿沸腾曲线依次变化。不会发生壁面温度急剧上升情况。 18．由式（7）t T R s ?= υγρσ2min ,在一定的s T t ,,,,υργσ?五个量中，只有υ ρ随压强变化最大，P 增加时，υρ的增加值将超过T s 的增值和γ的减少，最终使R min 随P 的增加而减小。 19．h=? w/(m 2·k) 20. h=67140 w/(m 2·k) 21.温度降为183℃ h=1585 w/(m 2·k) 与自然对流相比较， 485.01585 769 == 沸腾 自然对然h h 22．Q= w/(m 2·k) ,t w =℃ 第四章导热问题的数值解法 1 、重点内容：①掌握导热问题数值解法的基本思路； ②利用热平衡法和泰勒级数展开法建立节点的离散方程。 2 、掌握内容：数值解法的实质。 3 、了解内容：了解非稳态导热问题的两种差分格式及其稳定性。 §4—1导热问题数值求解的基本思想及内节点方程的建立由前述 3 可知，求解导热问题实际上就是对导热微分方程在定解条件下的积分求解，从而获得分析解。但是，对于工程中几何形状及定解条件比较复杂的导热问题，从数学上目前无法得出其分析解。随着计算机技术的迅速发展，对物理问题进行离散求解的数值方法发展得十分迅速，并得到广泛应用，并形成为传热学的一个分支——计算传热学（数值传热学），这些数值解法主要有以下几种： （１）有限差分法（ 2 ）有限元方法（ 3 ）边界元方法 数值解法能解决的问题原则上是一切导热问题，特别是分析解方法无法解决的问题。如：几何形状、边界条件复杂、物性不均、多维导热问题。 一．分析解法与数值解法的异同点： ?相同点：根本目的是相同的，即确定① t=f(x ， y ， z) ；② 。 ?不同点：数值解法求解的是区域或时间空间坐标系中离散点的温度分布代替连续的温度场；分析解法求解的是连续的温度场的分布特征，而不是分散点的数值。 数值求解的基本思路及稳态导热内节点离散方程的建立 二．解法的基本概念 ?实质 对物理问题进行数值解法的基本思路可以概括为：把原来在时间、空间坐标系中连续的物理量的场，如导热物体的温度场等，用有限个离散点上的值的集合来代替，通过求解按一定方法建立起来的关于这些值的代数方程，来获得离散点上被求物理量的值。该方法称为数值解法。 这些离散点上被求物理量值的集合称为该物理量的数值解。 2 、基本思路：数值解法的求解过程可用框图 4-1 表示。 由此可见： 1 ）物理模型简化成数学模型是基础； 2 ）建立节点离散方程是关键； 3 ）一般情况微分方程中，某一变量在某一坐标方向所需边界条件的个数等于该变量在该坐标方向最高阶导数的阶数。 ?数值求解的步骤 如图 4-2 （ a ），二维矩形域内无内热源、稳态、常物性的导热问题采用数值解法的步骤如下：（1）建立控制方程及定解条件 控制方程：是指描写物理问题的微分方程 针对图示的导热问题，它的控制方程（即导热微分方程）为：（ a ） 边界条件： x=0 时， x=H 时， 当 y=0 时， 数 值 传 热 学 近代发展及数值方法 建环：屈锐 2011年10月5日 数值传热学的发展史及数值方法 一、计算传热学的发展史 首先，计算传热学（Numerical Heat Transfer）与计算流体动力学(Computational Fluid Dynamics)之间的关系密切，可以认为，他们的主要研究内容是一致的，因此，计算传热学的发展史很大程度上也就是计算流体动力学的发展史，但他们之间还有不少区别，流体动力学的一个主要研究内容是讨论无粘流动及跨、超音速流动数值计算中的一些特殊问题。应用计算机和数值方法求解流动及传热问题在全世界范围内逐渐形成规模而且得出有益的结果，大致始于60年代，故从60年代起，可以把数值传热学的发展过程分为3个阶段： 1、萌芽初创阶段 主要有以下重大事件： （1）交错网格的提出。初期的数值传热学出现的两大困难之一是，网格设置不当时会得出具有不合理的压力场的解。1965年美国科学家首先提出了交错网格的思想，有效解决了这一难题，促使了求解NS 方程的原始变量法的发展。 （2）对流项差分迎风格式的再次确认。初期发展遇到的另一难题是 对流项采用中心差分时，对流速较高的情况的计算会得出振荡的解，1966年，科学家撰稿介绍了迎风格式在求解可压缩流体及非稳态层流流动中的作用，使流动与对流换热问题的求解建立在一个健壮的数值方法上发展。 （3）世界上第一本介绍流体及计算传热学的杂志于1966年创刊。（4）求解抛物型流动的P-S方法出现。由于受到计算机资源的限制，边界层类型问题的数值计算得到更多的关注，如何把有限个节点数目都充分利用起来成为了一个重要的问题。 （5）1969年Spalding在英国帝国理工学院创建了CHAM，旨在把他们研究组的成果推广应用到工业界。 （6）1972年SIMPLE算法问世。所谓分离式的求解方法应运而生，这个算法的基本思路是，在流场迭代求解的任何一个层次上，速度场都必须满足质量守恒方程，这一思想被以后的大量数值计算实例证明，是保证流场迭代计算收敛的一个十分重要的原则。 1974年美国学者提出了采用微分方程来生成适体坐标的方法。由于有限元法对不规则区域有很强的适应性，有限差分法与有限容积法则对复杂区域的适应能力很差，但对于流动问题的数值处理则要比有限元法容易得多。TTM方法的提出，为有限差分法与有限容积法处理不规则边界问题提出了一条崭新的道路。 2、开始走向工业应用阶段陶文铨 数值传热学 第二版 第五章 5-2
第四章编程题
《传热学期末复习试题库》含参考答案
传热学第四版课后题答案第五章
计算传热学中国石油大学(华东)第四章大作业
(完整)传热学第一章答案第四版-杨世铭-陶文铨汇总,推荐文档
数值传热学陶文铨第四章作业
传热学第五版课后习题答案(1)汇编
第五章传热学
传热学-第一章习题答案
传热学第五章答案
数值传热学第五章作业
传热学答案+第五版+章熙民(完整版)
传热学第1章答案
传热学第5.7章答案
第四章导热题的数值解法
数值传热学报告