下载文档原格式
高等传热复习题传热是指热量在物体间或物体内部的传递过程。
在高等传热学中，我们研究了传热的机制和方法，以及如何计算传热过程中的热量转移和温度变化。
本文将为大家提供一些高等传热学的复习题，以帮助大家巩固知识和提高解题能力。
1. 请简要解释传热方式中的传导、传动和辐射。
2. 一块厚度为5 cm的钢板，长度和宽度都是20 cm。
钢板顶部受到一个温度为100°C的热源作用，底部则接触着一个恒温为20°C的冷源。
假设钢板的导热系数为40 W/(m·K)，计算在稳态下，钢板内部的温度分布。
3. 一个圆柱形水管的半径为2 cm，长度为50 cm。
水管内部的温度分布满足稳态条件，圆柱体表面的温度为100°C，水管底部与一个恒温为20°C的冷源接触。
圆柱体的导热系数为0.5 W/(m·K)，求圆柱体内部的温度分布。
4. 辐射传热是通过辐射方式进行的热量传递。
请简要介绍辐射传热的基本原理。
5. 一块长宽均为20 cm的黑色平板，其表面辐射率为0.8，温度为300 K。
平板的一侧受到恒温为100 K的冷源作用，另一侧则环境温度为25°C。
计算平板从一侧向另一侧的净辐射热通量。
6. 在传热过程中，常用形状的传热表面有什么特点？请简要介绍。
7. 计算传热过程中的对流传热的换热系数是非常重要的。
请简要说明对流传热的基本原理，并介绍如何计算换热系数。
8. 一块长宽均为30 cm的金属板，表面具有一定的粗糙度。
当板的一侧受到恒温为100°C的热源作用时，板的另一侧与一个恒温为20°C的冷源接触。
试根据实验数据计算金属板的对流传热系数。
9. 请简要说明传热过程中的热阻和热导的概念，并介绍如何计算热阻。
10. 在传热过程中常用的传热模型有哪些？请简要介绍各个模型的适用范围。
11. 热管是一种特殊的传热装置，用于高效传输热量。
请简要说明热管的工作原理，并描述其应用领域。
第一章1-4、试写出各向异性介质在球坐标系)(ϕθ、、r 中的非稳态导热方程，已知坐标为导热系数主轴。
解：球坐标微元控制体如图所示：热流密度矢量和傅里叶定律通用表达式为：→→→∂∂+∂∂+∂∂-=∆-=k T r k j T r k i r T k T k q r ϕθθϕθsin 11'' （1-1）根据能量守恒：st out g in E E E E ••••=-+ϕθθρϕθθϕϕθθϕθd drd r tT c d drd r q d q d q dr r q p r sin sin 22∂∂=+∂∂-∂∂-∂∂-• （1-2） 导热速率可根据傅里叶定律计算：ϕθθθθd r dr Tr k q sin ⋅∂∂-= （1-3） 将上述式子代入（1-4-3）可得到)51(sin sin )sin ()sin (sin )(222-∂∂=+⋅⋅∂∂∂∂+⋅⋅∂∂∂∂+⋅⋅∂∂⋅∂∂⋅ϕθθρϕθθϕθϕθϕϕθθθθϕθθϕθd drd r t T c d drd r q d rd dr T r k rd d dr T r k d d dr r T r k r p r 对于各向异性材料，化简整理后可得到：tTc q T r k T r k r T r r r k pr ∂∂=+∂∂+∂∂∂∂+∂∂∂∂⋅ρϕθθθθθϕθ2222222sin )(sin sin )( （1-6）第二章2-3、一长方柱体的上下表面(x=0，x=δ)的温度分别保持为1t 和2t ，两侧面(L y ±=)向温度为1t 的周围介质散热，表面传热系数为h 。
试用分离变量法求解长方柱体中的稳态温度场。
解：根据题意画出示意图：（1）设f f f t t t t t t -=-=-=2211,,θθθ，根据题意写出下列方程组⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=+∂∂==∂∂======∂∂+∂∂00000212222θθλθθθδθθθθh y L y y y x x y x （2-1） 解上述方程可以把θ分解成两部分I θ和∏θ两部分分别求解，然后运用叠加原理∏+=θθθI 得出最终温度场，一下为分解的I θ和∏θ两部分： （2）首先求解温度场I θ用分离变量法假设所求的温度分布),(y x I θ可以表示成一个x 的函数和一个y 的函数的乘积，即)()(),(11y Y x X y x I =θ （2-2）将上式代入I θ的导热微分方程中，得到012121212=+X dyY d Y dx X d ，即21''11''1ε=-=Y Y X X ，上式等号左边是x 的函数，右边是y 的函数，只有他们都等于一个常数时才可能成立，记这个常数为2ε。
高等传热学作业修订版 IBMT standardization office【IBMT5AB-IBMT08-IBMT2C-ZZT18】第一章1-4、试写出各向异性介质在球坐标系)(ϕθ、、r 中的非稳态导热方程，已知坐标为导热系数主轴。
解：球坐标微元控制体如图所示：热流密度矢量和傅里叶定律通用表达式为：→→→∂∂+∂∂+∂∂-=∆-=k T r k j T r k i r T k T k q r ϕθθϕθsin 11'' （1-1）根据能量守恒：st out g in E E E E ••••=-+ϕθθρϕθθϕϕθθϕθd drd r tT c d drd r q d q d q dr r q p r sin sin 22∂∂=+∂∂-∂∂-∂∂-• （1-2） 导热速率可根据傅里叶定律计算：ϕθθθθd r dr Tr k q sin ⋅∂∂-= （1-3） 将上述式子代入（1-4-3）可得到)51(sin sin )sin ()sin (sin )(222-∂∂=+⋅⋅∂∂∂∂+⋅⋅∂∂∂∂+⋅⋅∂∂⋅∂∂⋅ϕθθρϕθθϕθϕθϕϕθθθθϕθθϕθd drd r t T c d drd r q d rd dr T r k rd d dr T r k d d dr r T r k r p r 对于各向异性材料，化简整理后可得到：tTc q T r k T r k r T r r r k pr ∂∂=+∂∂+∂∂∂∂+∂∂∂∂⋅ρϕθθθθθϕθ2222222sin )(sin sin )( （1-6）第二章2-3、一长方柱体的上下表面(x=0，x=δ)的温度分别保持为1t 和2t ，两侧面(L y ±=)向温度为1t 的周围介质散热，表面传热系数为h 。
试用分离变量法求解长方柱体中的稳态温度场。
解：根据题意画出示意图：（1）设f f f t t t t t t -=-=-=2211,,θθθ，根据题意写出下列方程组⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=+∂∂==∂∂======∂∂+∂∂00000212222θθλθθθδθθθθh y L y y y x x y x （2-1）解上述方程可以把θ分解成两部分I θ和∏θ两部分分别求解，然后运用叠加原理∏+=θθθI 得出最终温度场，一下为分解的I θ和∏θ两部分： （2）首先求解温度场I θ用分离变量法假设所求的温度分布),(y x I θ可以表示成一个x 的函数和一个y 的函数的乘积，即)()(),(11y Y x X y x I =θ （2-2）将上式代入I θ的导热微分方程中，得到012121212=+X dy Y d Y dx X d ，即21''11''1ε=-=Y Y X X ，上式等号左边是x 的函数，右边是y 的函数，只有他们都等于一个常数时才可能成立，记这个常数为2ε。
高等传热学问题及答案1. 简述三种基本传热方式的传热机理并用公式表达传热定律；传热问题的边界条件有哪两类？2. 有限元法求解传热问题的基本思想是什么？基本求解步骤有哪些？同有限差分方法相比其优点是什么？3. 什么是形函数？形函数的两个最基本特征是什么？4. 加权余量法是建立有限元代数方程的基本方法，请描述四种常见形式并用公式表达。
5. 特征伽辽金法（CG ）在处理对流换热问题时遇到什么困难？特征分离法（CBS ）处理对流换热问题的基本思想是什么？第一题：（1）热传导传热传导模式是因为从一个分子到另一个分子的能量交换，没有分子的实际运动，如果自由电子存在，也可能因为自由电子的运动。
因此，这种形式的热输送在很大程度上取决于介质的性质，如果存在温度差，热传导发生在固体，液体和气体。
书上补充：当两个物体有温差，或者物体内部有温度差时，在物体各部分之间不发生相对位移的情况下，物体微粒（分子，原子或自由电子）的热运动传递了热量。
（2）热对流()a w T T h q -=（牛顿冷却定律） 存在于液体和气体中的分子具有运动的自由，它们随身携带的能量（热量），从热区域移动到冷区域。
由于在液体或气体的宏观运动，热量传递从一个地区到另一个地方 ,加上流体内的热传导能量传递,称为对流换热。
对流可能是自然对流、强制对流，或混合对流。
百度补充：对流仅发生于流体中，它是指由于流体的宏观运动使流体各部分之间发生相对位移而导致的热量传递过程。
由于流体间各部分是相互接触的，除了流体的整体运动所带来的热对流之外，还伴生有由于流体的微观粒子运动造成的热传导。
在工程上，常见的是流体流经固体表面时的热量传递过程，称之为对流传热。
（3）辐射4w T q εσ= （ 斯蒂藩-玻耳兹曼定律）任何（所有）物体和任何（所有）温度都能产生热辐射。
（绝对零度以上）这是唯一一种发生热传递不需要介质的方式。
热辐射本质上是从物体的表面发射电磁波，由电磁波携带能量进行能量传输。
《高等传热学》课程自学及作业安排2014届硕士研究生适用本课程教学方式：以自学为主，教师指导为辅。
考核方法：开卷笔试（50%）+平时成绩（作业及讲课30%）+两次大作业（20%）一、教学资料1.教材孙德兴编.高等传热学—导热与对流的数理解析.北京：中国建筑工业出版社，2005（图书馆均可借到）2.主要参考书张靖周编.高等传热学.北京：科学出版社，2009*王瑞金等编.Fluent技术基础与应用实例.北京：清华大学出版社，20073.参考资料[1]杨强生，高等传热学.上海：上海交通大学出版社，1996[2][美]E.R.G.埃克特，R.M.德雷克著，航青译.传热与传质分析.北京：科学出版社，1983[3][美]M. N.奥齐西克,俞昌铭主译.热传导.北京：高等教育出版社，1983[4]杨强生.对流传热与传质.北京：高等教育出版社，1985[5]赵镇南译.对流传热与传质（第4版）.北京：高等教育出版社，2007*[6][美]E.M.斯帕罗，R.D.塞斯著，顾传保，张学学译.辐射传热.北京：高等教育出版社，1982*[7]陶文铨编著.数值传热学.西安:西安交通大学出版社，1988[8]周俊杰等编. FLUENT工程技术与实例分析.北京：中国水力水电出版社，2010（除*外，均提供电子版）4.课件、教案、FLUENT软件及其他提供光盘！二、自学、收集整理资料及讲课1.自学根据教案及课件提前查资料并自学相关内容。
如：2.收集整理资料及讲课每三位同学负责一至二次课内容，具体分工自行商量。
内容包括：(1)收集整理资料按照教案要求，收集、整理、加工相关教学资料，如“典型一维稳态导热现象（参考文献[1]PP27-40）”，形成电子版提交到qq群，供全班同学共享。
(2)讲课其中一位同学讲解该次课教学内容，时间为45分钟，重点讲解教案中提出的“重点需要理解的问题”；另一位同学讲解作业，时间15分钟，重点讲解分析思路。
高等传热学导热练习题1. 试求圆柱坐标),,(z r φ的拉梅系数。
圆柱坐标(,,)r z φ和直角坐标(,,)x y z 的 关系是：cos x r φ=，sin y r φ=，z z = 解：由题目条件得：2222221cos sin 1x y z a r r r φφ∂∂∂⎛⎫⎛⎫⎛⎫=++=+= ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭得：11a =()()22222222sin cos x y z a r r r φφφφφ⎛⎫⎛⎫⎛⎫∂∂∂=++=−+= ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭得：2a r =222231x y z a z z z ∂∂∂⎛⎫⎛⎫⎛⎫=++= ⎪ ⎪ ⎪∂∂∂⎝⎭⎝⎭⎝⎭得：31a =123a a a a r ==3. 一维无限大平板，0≤x ≤L ，初始温度为F(x)。
当时间0>τ时，x=0处与x=L 处的边界温度维持零度。
试求时间0>τ时，平板内温度),(τx t 的表达式。
并求当初始温度F(x)=t 0=常数这种特殊情况下的温度),(τx t 。
解：该导热问题的数学描写为：()()()()()()22,,1,0,00,0,0,0t x t x x L x t t L t x F x τττατττ⎧∂∂=<<>⎪∂∂⎪⎪=⎨⎪=⎪⎪=⎩ 分离变量：()()(),t x X x ττ=⋅Γ 代入温度微分方程得：()()()()22211d X x d const X x dx d τβαττΓ==−=Γ得时间函数：()2e αβττ−Γ=空间变量的特征值问题为：()()()()222000d X x X x dxX X L β⎧+=⎪⎨⎪==⎩查表得：()(),sin m m X x x ββ=，()12m N Lβ=，m β是()sin 0m L β=的正根 温度通解为：()()21,,m m m m t x c X x e αβττβ∞−==∑代入初始条件可得：()()(),Lm mm X x F x dxc N ββ=⎰将上式代入温度的通用级数解，可得：()()()()2012,sin sin m L m m m t x x F x dx x e Lαβττββ∞−='''=⋅⋅∑⎰ 对于()0F x const t ==的情形，可得：()()()2011cos 2,sin m m m m m L t t x x e Lαβτβτββ∞−=−=⋅∑4. 一维无限大平板，0≤x ≤L ，初始温度为F(x)。
上海理工大学高等传热学试题及答案(总2页)--本页仅作为文档封面，使用时请直接删除即可----内页可以根据需求调整合适字体及大小--1.试求出圆柱坐标系的尺度系数，并由此导出圆柱坐标系中的导热微分方程。
2 .一无限大平板，初始温度为T 0；τ>0时，在x = 0表面处绝热；在x =L 表面以对流方式向温度为t f 的流体换热。
试用分离变量法求出τ>0时平板的温度分布（常物性）。
（需求出特征函数、超越方程的具体形式，范数（模）可用积分形式表示）。
（15分）3.简述近似解析解——积分法中热层厚度δ的概念。
答：近似解析解：既有分析解的特征：得到的结果具有解析函数形式，又有近似解的特征：结果只能近似满足导热解问题。
在有限的时间内，边界温度的变化对于区域温度场的影响只是在某一有限的范围内，把这个有限的范围定义为热层厚度δ。
4.与单相固体导热相比，相变导热有什么特点答：相变导热包含了相变和导热两种物理过程。
相变导热的特点是1.固、液两相之间存在着 移动的交界面。
2.两相交界面有潜热的释放（或吸收）对流部分（所需量和符号自己设定）1 推导极坐标系下二维稳态导热微分方程。
2 已知绕流平板流动附面层微分方程为y uy u V x u u 22∂∂=∂∂+∂∂ν取相似变量为：x u y νη∞= x u f νψ∞=写出问题的数学模型并求问题的相似解。
3 已知绕流平板流动换热的附面层能量积分方程为：⎰=∞∂∂=-δ00)(y y ta dy t t u dx d当Pr<<1时，写出问题的数学模型并求问题的近似积分解及平均Nu （取三次多项式）。
4 写出常热流圆管内热充分发展流动和换热问题的数学模型并求出速度和温度分布及Nu x .O x辐射1．请推导出具有n个表面的净热流法壁面间辐射换热求解公式，并简要说明应用任一种数值方法的求解过程。
2．试推导介质辐射传递方程的微分形式和积分形式，要求表述出各个步骤和结果中各个相关量的含义。
传热学实验试题及答案高中一、选择题（每题2分，共10分）1. 传热的基本方式有哪些？A. 导热、对流、辐射B. 导热、对流、蒸发C. 导热、对流、扩散D. 蒸发、对流、辐射2. 以下哪个不是影响导热速率的因素？A. 材料的导热系数B. 温度差C. 材料的厚度D. 材料的颜色3. 对流换热的特点是什么？A. 只发生在固体内部B. 只发生在液体和气体中C. 需要流体的流动D. 与流体的流动无关4. 辐射换热的特点是？A. 需要介质B. 与温度无关C. 可以在真空中进行D. 只发生在固体表面5. 以下哪个公式用于计算稳态导热的热流量？A. Q = h * A * ΔTB. Q = k * A * ΔT / dC. Q = c * ρ * ΔTD. Q = m * L * ΔT / t答案：1-5 A B C C B二、填空题（每空1分，共10分）1. 导热的基本定律是________，表达式为Q = k * A * ΔT / d。
2. 当流体的流动状态发生变化时，换热方式可能由________转变为对流换热。
3. 辐射换热不需要介质，可以在________中进行。
4. 影响对流换热的因素包括流体的________、________和流速。
5. 稳态导热条件下，物体内部某点的温度随时间________。
答案：1. 傅里叶定律 2. 导热 3. 真空 4. 温度、压力 5. 不变三、简答题（每题5分，共20分）1. 简述导热、对流和辐射三种传热方式的区别。
2. 描述稳态导热和非稳态导热的特点。
3. 解释为什么说辐射换热是所有物体都具有的自然现象。
4. 举例说明在日常生活中如何应用对流换热的原理。
答案：1. 导热是热量通过物体内部分子振动和自由电子运动传递的过程，主要发生在固体中；对流是流体中温度不同的各部分之间通过宏观位移传递热量的过程，主要发生在液体和气体中；辐射是物体因温度产生的电磁波向外传递热量的过程，可以在真空中进行，所有物体都具有辐射能力。
高等传热学复习题1. 太空飞行物伸出的细长散热棒，以辐射方式与外部进行换热，棒长L 、截面积A 、截面周长U 、导热系数λ、发射率ε、棒根部温度t 0 ，外部空间为绝对黑体，写出该问题的完整数学描述。
2. 半径为R 的实心球，初时温度为t 0，突然放入t f 冷水中，已知球的物性λ、c 、ρ及表面传热系数h ，写出球冷却的完整数学描述。
3. 直径为d 、单位长度电阻为R 、发射率为ε的金属棒，初始时与温度为T ∞的环境处于热平衡状态，后通过电流I ，已知棒与环境的表面传热系数为h 。
试导出通电流期间金属棒温度随时间变化的规律，并写出处于新的热平衡状态的条件。
（不用求解）4. 大平板：δ,Φ1） 已知两侧为对称第三类边界条件,h ,f t 求t 的分布；2） 一侧为第三类边界条件,h ,f t 另一侧绝热, 求t 的分布。
3） 一侧为第一类边界条件，另一侧为绝热，,求t 的分布。
4） 两侧为相同的第一类边界条件,求t 的分布。
5） 两侧为不同的第一类边界条件,求t 的分布。
5. 厚为L 、导热系数λ =1.5W/(m K)的浇注混凝土墙，两边保持温度为20℃，由于混凝土的固化，单位体积释放100W/m 2的化学热能。
若要求浇注时墙内任意处每米墙厚的温度梯度不大于50℃，墙的最大厚度是多少？6. 敷设肋片就一定能强化传热? 增加散热量满足的条件?解：敷设肋片时: ()()0sh()ch()ch()sh()mH h m mH ΦmH h m mH λλ+=+ 不敷设肋片时:0nf ΦhA θ= ()()00sh()ch()ch()sh()nf mH h m mH A mmH h m mH ΦΦhA λλθλθ++= ()()sh()ch()ch()sh()nf mH h m mH Φm Φh mH h m mH λλλ+=+ th()11th()nf m mH Φhh ΦmH m λλ+=+ 1>=< >1 增强换热；=1 不增强不减弱；<1 减弱换热。
高等传热学复习题答案1. 试述傅里叶定律的物理意义及其数学表达式。
傅里叶定律描述了在稳态条件下，热量通过材料的传导过程。
其物理意义是热量的传递速率与温度梯度的负值成正比，且与材料的热导率有关。
数学表达式为：\( q = -k \frac{dT}{dx} \)，其中 \( q \) 表示热量传递速率，\( k \) 表示材料的热导率，\( \frac{dT}{dx} \) 表示温度梯度。
2. 什么是热对流？请简述热对流的两种主要类型。
热对流是指流体中热量的传递过程，它依赖于流体的宏观运动。
热对流的两种主要类型为自然对流和强制对流。
自然对流是由流体内部密度差异引起的，而强制对流则是由外部力（如风扇或泵）驱动的流体运动。
3. 简述辐射换热的基本原理。
辐射换热是指物体之间通过电磁波传递能量的过程。
它不需要任何介质，可以在真空中进行。
辐射换热的基本原理是物体根据其温度和表面特性发射和吸收辐射能。
斯特藩-玻尔兹曼定律和普朗克定律是描述辐射换热的基本定律。
4. 试分析在不同边界条件下，热传导问题的解法。
在不同的边界条件下，热传导问题的解法会有所不同。
例如，在狄利克雷边界条件下，物体表面的温度是已知的；在诺伊曼边界条件下，物体表面的热流密度是已知的；而在罗宾边界条件下，物体表面的热流密度与温度的函数关系是已知的。
对于这些不同的边界条件，可以采用分离变量法、有限差分法或有限元法等方法求解。
5. 描述在不同工况下，流体流动的类型及其特点。
流体流动的类型通常根据流动的雷诺数（Re）来分类。
当Re小于2300时，流动为层流，特点是流线平行，无涡旋；当Re大于4000时，流动为湍流，特点是流线混乱，存在涡旋。
在过渡流区域（2300 < Re < 4000），流动状态不稳定，可能同时存在层流和湍流的特点。
6. 试解释热辐射中的黑体、灰体和选择性辐射体的概念。
黑体是指能够吸收所有入射辐射的物体，其辐射能力与温度有关，遵循斯特藩-玻尔兹曼定律。
1、 两种粗糙材料相互接触时会产生接触热阻，造成温度分布在界面上不再连续。
假设两种材料的接触热阻为c R ，导热系数分别是1λ和2λ，请写出稳态导热时两种材料接触面上的边界条件？dxdt t t Rc dx dt 222111)(λλ-=-=- 2、 对一维非稳态导热问题利用显示差分进行数值求解，请写出对流边界节点的差分方程式，并从物理意义上分析其稳定性条件是什么？[传热学书第四版176页] 3、 平板上流动的摩擦系数通常表示为m x fx C c -⋅=Re ，式中C 和m 为常数。
试以fx c 的关系式为基础，给出L 长度上的平均努塞尔数的一般表达式？[传热学书第四版212页]答案：3/1Pr Re m l l C Nu =4、 对外掠平板层流边界层流动和换热进行相似求解，Blasius 最早提出可以通过引入一个相似变量x U y νη∞=，将偏微分控制方程组变换为常微分方程。
请根据对边界层厚度的数量级分析的结果说明η的物理意义是什么？5、 一根直径为5mm 的不锈钢电缆，通过600A 的电流。
电缆的单位长度电阻率为m Ω⨯-4106，电缆外面包裹有一层导热系数为0.5W/(m ⋅℃)的绝缘层。
该电缆置于25℃的大气环境中，外表面与周围环境之间的表面传热系数为25W/(m 2⋅℃)，问绝缘层多厚时其内表面温度最低，并求此时内表面温度？[与传热学书第四版92页2-15题目类似，题目2-15答案如下]6、 2dA 为圆筒壁上距离为微元dx 的两个截面3A 和4A 之间的筒壁内侧面积。
已知1A 和3A 之间的角系数为222223,1242rr x x r x +-+=ϕ，求1A 和2dA 之间的角系数？［这个题目有没有问题］ 7、 单层大平壁的两个表面分别维持均匀的温度t 1和t 2(t 1＞t 2)，平壁内带有均匀分布的内热源，其体积发热率为v q ，求平壁内的温度分布，并定性画出平壁内的温度分布曲线？试分析在什么条件下，最大的温度可能出现在平壁内部？[传热学书第四版71页]8、 对一维非稳态导热问题利用显示差分进行数值求解。
第九章4．一工厂中采用0.1MPa 的饱和水蒸气在—金属竖直薄壁上凝结，对置于壁面另一侧的物体进行加热处理。
已知竖壁与蒸汽接触的表面的平均壁温为70 ℃，壁高1.2m ，宽300 mm 。
在此条件下，一被加热物体的平均温度可以在半小时内升高30℃，试确定这一物体的平均热容量(不考虑散热损失)。
解:本题应注意热平衡过程，水蒸气的凝结放热量应等于被加热物体的吸热量。
P=0.1Mpa=105Pa,t s =100℃,r=2257.1kJ/kg, t m =21( t s + t w )= 21(100+70) ℃=85℃。
查教材附录5，水的物性为：ρ=958.4kg/m 3;λ=0.683 W ／(m 2·℃);µ=282.5×10-6N·s/m 2假设流态为层流：4132)(13.1⎥⎦⎤⎢⎣⎡-=w s t t l r g h μλρ416332)70100(2.1105.282102257683.081.94.95813.1⎥⎦⎤⎢⎣⎡-⨯⨯⨯⨯⨯⨯⨯=- W ／(m 2·℃) =5677 W ／(m 2·℃)36102257105.2822.13056774)(4Re ⨯⨯⨯⨯⨯⨯=-=-r t t hl w s c μ=1282<1800流态为层流，假设层流正确Φ=ωl t t h w s )(-=5677×（100−70）×1.2×0.3W=61312W凝结换热量=物体吸热量Φ∆τ=mc p ∆t61068.330603061312⨯=⨯⨯=∆Φ∆=t mc p τJ/℃ 16．当液体在一定压力下做大容器饱和沸腾时，欲使表面传热系数增加10倍，沸腾温差应增加几倍?如果同一液体在圆管内充分发展段做单相湍流换热，为使表面传热系数增加10倍，流速应增加多少倍?维持流体流动所消耗的功将增加多少倍?设物性为常数。
1-5 椭球坐标系（),,ϕθη由η=常数的椭球面，θ=常数的双曲线面和ϕ=常数的平面组成。
如果椭球坐标系与直角坐标系的关系为：θηϕθηϕθηcos sin sin cos sin Ach z Ash y Ash x === 试证明该椭球坐标系的拉梅系数为：1H =ηH =2222cos sin θηθηsh ch A + θH H =12222cos sin θηθηsh ch A += ψH H =1θηsin Ash =并证明椭球坐标系中拉普拉斯算子的表达式为：⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂+∂∂+=∇θθθηηηθηθηt t t cth t sh ch A t cot )cos sin (12222222222ψθη22222sin 1∂∂+tsh A 解：（1）由式1 -2-18知 2221)()()(ηηηη∂∂+∂∂+∂∂==zy x H H 22222222222sin sin cos cos cos θηψθηψθηch A sh A sh A ++=θηsin Ash =（2）由式1 -2 - 25知ψθηθθθηηηθηθηψθηθθηθηθηθηθηηθηθηψθθθηηηψψψθθθηηηψψηθηη22222222222222222222222222222223332233222122321321223322333212322212312sin 1cot )cos sin (1sin 1sin cos sin sin )cos sin (11]11[1])([1)( )()([1)3,2,1()(1∂∂+⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂+∂∂+=∂∂+⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂+∂∂+=∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂=∂∂+∂∂∂∂+∂∂+∂∂∂∂+∂∂+∂∂∂∂=∂∂∂∂+∂∂∂∂+∂∂∂∂==∂∂∂∂=∇∑=tsh A t t t cth t sh ch A tsh A t Ash Ach t t Ash Ach t sh ch A t H t H H t t H H t H H t H H H t H H H t H t H t H t H H H H tH H t H H t H H H i x t H H x Ht ii i i3-2 大平壁的初始温度均匀为0t ，从某一时刻起其两表面的温度突然降为w t 并保持不变，试求:（1）写出该导热问题的数学描述； （2）用分离变量法求解平壁中的温度场。
高等传热学复习题1．简述求解导热问题的各种方法和傅立叶定律的适用条件。
答：导热问题的分类及求解方法：按照不同的导热现象和类型，有不同的求解方法。
求解导热问题，主要应用于工程之中，一般以方便，实用为原则，能简化尽量简化。
直接求解导热微分方程是很复杂的，按考虑系统的空间维数分，有0维，1维，2维和3维导热问题。
一般维数越低，求解越简单。
常见把高维问题转化为低维问题求解。
有稳态导热和非稳态导热，非稳态导热比稳态导热多一个时间维，求解难度增加。
有时在稳态解的基础上分析非稳态稳态，称之为准静态解，可有效地降低求解难度。
根据研究对象的几何形状，又可建立不同坐标系，分平壁，球，柱，管等问题，以适应不同的对象。
不论如何，求解导热微分方程主要依靠三大方法：甲．理论法乙．试验法丙．综合理论和试验法理论法：借助数学、逻辑等手段，根据物理规律，找出答案。
它又分：分析法；以数学分析为基础，通过符号和数值运算，得到结果。
方法有：分离变量法，积分变换法(L a p l a c e变换，F o u r i e r变换)，热源函数法，G r e e n函数法，变分法，积分方程法等等，数理方程中有介绍。
近似分析法：积分方程法，相似分析法，变分法等。
分析法的优点是理论严谨，结论可靠，省钱省力，结论通用性好，便于分析和应用。
缺点是可求解的对象不多，大部分要求几何形状规则，边界条件简单，线性问题。
有的解结构复杂，应用有难度，对人员专业水平要求高。
数值法：是当前发展的主流，发展了大量的商业软件。
方法有：有限差分法，有限元法，边界元法，直接模拟法，离散化法，蒙特卡罗法，格子气法等，大大扩展了导热微分方程的实用范围，不受形状等限制，省钱省力，在依靠计算机条件下，计算速度和计算质量、范围不断提高，有无穷的发展潜力，能求解部分非线性问题。
缺点是结果可靠性差，对使用人员要求高，有的结果不直观，所求结果通用性差。
比拟法：有热电模拟，光模拟等试验法：在许多情况下，理论并不能解决问题，或不能完全解决问题，或不能完美解决问题，必须通过试验。
Homework 5.21. Slug flow in a tube(u(r) = V) with a fully developed temperature profile:- Constant heat flux - What is T(r)?·Remember dT m /dx = constant·Use dT/dr = 0 at r = 0 and T w at r = R to get ingration constants - What is Nu? A: 根据能量方程，221p p T T T T c u c v k rk r x r r rx ρρ∂∂∂∂∂⎛⎫+=+ ⎪∂∂∂∂∂⎝⎭(1) 对于充分发展的管内弹状流，u=0, 且由于mdT T x dx∂=∂= constant ，因此(1)式可化简为： 1m p dT Tc Vk r dx r r rρ∂∂⎛⎫= ⎪∂∂⎝⎭(2) 解之可得，212(r)ln 4p m V c dT T r C r C k dxρ=⋅⋅++ (3) 根据边界条件，0,, T(r)=T w Tr rr R ∂==∂= 代入(3)中，得到C 1=0, 224p m w V c dT C T R k dxρ=-⋅平均温度201()2Rm T T r rdr R ππ=⎰ 28()m w m p dT kT T dx V c Rρ=- 因此，得到22222(r)(1)22()4m w m w w m dT VR r r T T T T T T dx R Rα=+⋅⋅-=-+-2288q''h()()()222p p m w m w m w m p mc c R V dT k kT T T T T T R dx R V c R Rρπππρ=-=⋅=⋅-=- 28h d h R Nu k k⋅⋅===2. Fully developed Poiseuille flow between parallel plates- Constant heat flux- Top and bottom at the same temperature - Neglect viscous dissipation - What is T(y)?·Remember dT m /dx = constant·Use T w at y = 0 and y = h to get integration constants - What is Nu? A: 根据能量方程，2222()p p T T T T c u c v k x y x yρρ∂∂∂∂+=+∂∂∂∂ (1)对于x 方向充分发展的Poiseuille 流，v=w=0 由于mdT T x dx∂=∂= constant ，(1)式可化为： 22p T T c u k x yρ∂∂=∂∂ (2)根据01()hm u u y dy h =⎰，且22()6()m u y y y u h h=-，代入得到22226()m m p dT y y Tu c k h h dx yρ∂-=∂解得温度分布为34122()()2p m m c u dT y y T y C y C k dx h h ρ=⋅-++ (3)根据Poiseuille 流边界条件：0,w y T T == ,w y h T T ==带入(3)式中解得，12p m mc u dT C h kdxρ=-, 2w C T =因此，得到温度分布为342()()22p m p m m m w c u c u dT dT y y T y y T k dx h h k dxρρ=⋅--⋅+ 根据342220011()()6()22h hp m p m m m m w m m m c u c u dT dT y y y y T T y udy y T u dy hu hu k dx h h k dx h h ρρ⎡⎤⎡⎤==⋅--+-⎢⎥⎢⎥⎣⎦⎣⎦⎰⎰∴2140()17m w m m p dT kT T dx u c hρ=- ∴ 140q''h()()17p m w m w m mc dT kT T T T L dx L =-=⋅=- Nu 14017=Homework 5.31. Find the equation for the boundary layer thickness δ and for C f for:23()u y a b c d Uηηη=+++ Compare to the exact values and the fourth order equation solutions. A: 对于23()u y a b c d Uηηη=+++，根据边界条件， 0, u=0y =，∴ a = 0,y u U δ==，∴ 1 = b + c + d 对于平板，有220y uy =∂=∂，∴ c = 0根据0y u yδ=∂=∂，∴ b + 3d = 0由此可得，a = 0, b = 1.5, c = 0, d = -0.5.3()3122u y U ηη=- 2200333()222w y y u y UU yμτμμδδδ==∂=-=-⋅-=-∂1039(1)280i u u y d U U δδδ=-=⎰ 对于U = constant ，23-2i w UU x δδμρτδδ∂==∂积分后得到，xδ=与精确解相差7.2%，与4级近似解相差20.5%02223=1122y f uy U C U U U μτμδρρρ=∂∂==，且x δ=f C ==与精确解相差2.6%，与4级近似解相差5.5%2. Derive the ordinary differential equation and the boundary conditions for the Blasius solution energy equation for flow over a flat plate. (See pages 29-30) A: 根据能量方程，22()k p T T T c u v x y yρ∂∂∂+=∂∂∂w w T T T T θ∞-=-，22()k p c uv x y y θθθρ∂∂∂+=∂∂∂ 常壁温边界条件：0,0,1y y θθ===∞=根据Blasius solution δη≈===3/22222221122x x x x y y U y y y x θθηθθηηηηθθηθηηθθηθθηηνη∂∂∂∂-∂⎛⎫⎛⎫===⋅-⋅ ⎪ ⎪∂∂∂∂∂⎝⎭⎝⎭∂∂∂∂==∂∂∂∂⎛∂∂∂∂∂∂==⋅=⋅ ∂∂∂∂∂∂⎝代入能量方程中得到，)1'''''2p p U c f U c f f k x vx ηρθρηθθ⎛⎫-+-= ⎪⎝⎭即 2''Pr '0f θθ+= 边界条件为00;1ηθηθ===∞=，，3. Start from the local friction coefficient for flow over a flat plate, C f , on slide 31 andderive the average friction coefficient over the entire plate, C L . Show your work. A:2(x)12f C U τρ===01LL C L ==⎰4. A projectile in the form of a bluff-ended cylinder 20 cm in diameter and 60 cm long, moves through the air in the direction of its long axis at a velocity U of 100 m/s.212Dfrontal F C A U ρ=The drag coefficient, C, is equal to 1.0 for this object. Frontal area=πD 2/4, F D = total dragAssuming that the boundary layer thickness over the cylindrical surface of the projectile at a distance x from the leading edge is given by:xδ=and that the momentum thickness of the boundary layer is given by:772i δδ=Find what proportion of the total drag on the projectile is attributable to skin friction over the curved surface, assuming no pressure gradient in the boundary layer in the streamline direction.Data. Air kinematic viscosity is 0.15 cm 2/s Air density is 1.15 Kg/m 3Assumptions: Assume this is like a flat plate, then use the equation on slide 15.A: 22*()i w U dU U x dxδρρδτ∂+=∂ 若将圆柱表面看成平板，U = constant,0dUdx=，则22()772i w U U x xδδτρρ∂∂==∂∂因为xδ=∴4155x x x δ-⎛⎫⎛⎫∂∂∂=== ⎪ ⎪∂∂∂⎝⎭125772w U τρ-=1125500070.131372LLL w f dx U dx x dx τρ--====⎰⎰⎰12.16N22142D F D U πρ=⋅=180.64N DfF =6.7%Homework 5.41. Derive the equation relating q” to the temperature difference for natural convection driven flow in a round tube (like on slide 19 for parallel plates). - water properties at 25 C - tube length = 100 cm - tube diameter = 1 mm - constant heat flux- fully developed flow in a tube (from the first homework)a) Plot ∆T and Q (m 3/s) versus heat flux for fluxes of 500 to 5000 W/m 2 b) Is the assumption of fully developed flow valid? - hint: is L/(D Re) > 0.05c) Is the Boussinesq approximation valid here? A: 圆管内充分发展的流体，其速度分布为 22214R p r u z R μ⎛⎫∂=- ⎪∂⎝⎭平均速度2222201(1)2r 48RR p r R p u dr R z R zππμμ∂∂=⋅⋅-=-⋅∂∂⎰体积流量为44288R p R pQ R u z Lπππμμ∂∆===∂()()()22L LL L Lp gdy g T T dy g T T g T ρρρβρβρβ∞∞∞∞∞∞∆=-=-=-=∆⎰⎰ 根据能量平衡，''2p Qc T q RL ρπ∆=即4''28p R pc T q RL L πρπμ∆∆=将p ∆代入，得到()232''32p c gR q T Lρβμ∞=∆Pr ''''32R cond Gr q q =，''cond T q R λ∆= a) 根据以上推导过程可以得到，T ∆=Q =分别作出当2''500~5000/q W m =时T ∆与Q 随q’’变化的曲线，如下图所示。
7-4，常物性流体在两无限大平行平板之间作稳态层流流动，下板静止不动，上板在外力作用下以恒定速度U 运动，试推导连续性方程和动量方程。
解：按照题意0,0=∂∂=∂∂=xv y v v 故连续性方程0=∂∂+∂∂yv x u 可简化为0=∂∂xu因流体是常物性，不可压缩的，N-S 方程为 x 方向：)(12222yu x u v y p F y u v x u u x ∂∂+∂∂+∂∂-=∂∂+∂∂ρρ 可简化为022=∂∂+∂∂-yv x p F x ηy 方向)(12222yv x v v y p F y v v x v u y ∂∂+∂∂+∂∂-=∂∂+∂∂ρρ 可简化为0=∂∂=ypF y8-3，试证明，流体外掠平壁层流边界层换热的局部努赛尔特数为12121Re Prx Nu r =证明：适用于外掠平板的层流边界层的能量方程22t t t u v a x y y∂∂∂+=∂∂∂ 常壁温边界条件为0w y t t y ∞==→∞时，时，t=t引入量纲一的温度wwt t t t ∞-Θ=-则上述能量方程变为22u v a x y y∂Θ∂Θ∂Θ+=∂∂∂引入相似变量12Re ()y yx x ηδ===有11()(()22x x xηηηηη∂Θ∂Θ∂''==Θ-=-Θ∂∂∂()y y ηηη∂Θ∂Θ∂'==∂∂∂；22()U y x ηυ∞∂Θ''=Θ∂ 将上三式和流函数表示的速度代入边界层能量方程，得到1Pr 02f '''Θ+Θ=当Pr1时，速度边界层厚度远小于温度边界层厚度，可近似认为温度边界层内速度为主流速度，即1,f f η'==，则由上式可得Pr ()2d f d η''Θ'=-'Θ，求解可得 11()()Pr 2Pr(0)()erf ηηπΘ='Θ=则12120.564RePrx xNu =8-4，求证，常物性不可压缩流体，对于层流边界层的二维滞止流动，其局部努赛尔特数满足10.4220.57Re Pr x Nu =⋅证明：对于题中所给情况，能量方程可表示为22u v x y yθθθα∂∂∂+=∂∂∂其中，,,()u v y x ψψψθθηθ∂∂==-===∂∂ 故上式可转化为Pr02θζθ'''+⋅⋅= 经两次积分，得到0000Pr [exp()]2()Pr [exp()]2d d d d ημμζηηθμζηη∞-=-⎰⎰⎰⎰ 定义表面传热系数s x s q h T T ∞=-，则(0)q '= 进一步，进行无量纲化处理，引入局部努赛尔特数12(0)Re x x x h x Nu k ⋅'===其中1200Re (0)Pr [exp()]2xd d μθζηη∞'=-⎰⎰ 针对层流边界层的条件，查由埃克特给出的计算表如下：不同Pr 数下，常物性层流边界层，12Re x Nu -⋅的值故可看出，12Re x Nu -⋅=常数，进而，12()=x h xu k υ-∞⋅=1常数C ，由1m u C x ∞=⋅，得11212m C kh xυ-=⋅对于二维滞止流，m=1，则h 也为常数，从x=0到x 处的平均热导率h m 定义为1xm h hdx x =⎰故11112212120121m m x m C k C k h x dx x x m υυ--=⋅=⋅⋅+⎰， 则21m h h m =+，由此可看出， 在m=1时，努赛尔特数的近似解可以很好的表示为10.4220.57Re Pr x Nu =⋅ 同样的，我们也可以得到三维滞止流的近似解10.4220.76Re Pr x Nu =⋅9-1，试证明：圆管内充分发展流动的体积流量可表示为： ()0408p p Lr V i -=μπ9-2，常物性不可压缩流体在两平行平板间作层流流动，下板静止，上板以匀速U 运动，板间距为2b ，试证明充分发展流动的速度分布为⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛--=b y b y dx dp b b y U u 2222μ 证：二维流体质量、动量方程0=∂∂+∂∂yvx u ① ⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂-=∂∂+∂∂2222)(y u xu x py u v x u u μρ ②⎪⎪⎭⎫ ⎝⎛∂∂+∂∂+∂∂-=∂∂+∂∂2222)(y v xv y py v v x v u μρ ③ 在充分发展区，截面上只有沿流动方向的速度u 在断面上变化，法向速度v 可以忽略，因此可由方程①得:0=v ,0=∂∂xu④ 将式④代入③得到，0=∂∂yp，表明压力P 只是流动方向x 的函数，即流道断面上压力是均匀一致的进一步由式②得，t cons y udx dp tan 22=∂∂=μ ⑤相应的边界条件：Uu b y u y ====,20,0对⑤积分得：11C y dx dpyu +=∂∂μμ21221C y C y dxdp U ++=μ ddp b b u C μ-=21，02=C ⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛--=⇒b y b y dx dp b b y U u 2222μ1. 强迫流动换热如何受热物性影响？答：强迫对流换热与Re 和Pr 有关；加热与对流的粘性系数发生变化。
高等传热学作业Revised on November 25, 2020第一章1-4、试写出各向异性介质在球坐标系)(ϕθ、、r 中的非稳态导热方程，已知坐标为导热系数主轴。
解：球坐标微元控制体如图所示：热流密度矢量和傅里叶定律通用表达式为：→→→∂∂+∂∂+∂∂-=∆-=k T r k j T r k i r T k T k q r ϕθθϕθsin 11'' （1-1）根据能量守恒：st out g in E E E E ••••=-+ϕθθρϕθθϕϕθθϕθd drd r tT c d drd r q d q d q dr r q p r sin sin 22∂∂=+∂∂-∂∂-∂∂-• （1-2） 导热速率可根据傅里叶定律计算：ϕθθθθd r dr Tr k q sin ⋅∂∂-= （1-3） 将上述式子代入（1-4-3）可得到)51(sin sin )sin ()sin (sin )(222-∂∂=+⋅⋅∂∂∂∂+⋅⋅∂∂∂∂+⋅⋅∂∂⋅∂∂⋅ϕθθρϕθθϕθϕθϕϕθθθθϕθθϕθd drd r t T c d drd r q d rd dr Tr k rd d dr T r k d d dr r T r k r p r 对于各向异性材料，化简整理后可得到：tTc q T r k T r k r T r r r k p r ∂∂=+∂∂+∂∂∂∂+∂∂∂∂⋅ρϕθθθθθϕθ2222222sin )(sin sin )( （1-6） 第二章2-3、一长方柱体的上下表面(x=0，x=δ)的温度分别保持为1t 和2t ，两侧面(L y ±=)向温度为1t 的周围介质散热，表面传热系数为h 。
试用分离变量法求解长方柱体中的稳态温度场。
解：根据题意画出示意图：（1）设f f f t t t t t t -=-=-=2211,,θθθ，根据题意写出下列方程组⎪⎪⎪⎪⎩⎪⎪⎪⎪⎨⎧=+∂∂==∂∂======∂∂+∂∂00000212222θθλθθθδθθθθh y L y y y x x y x （2-1） 解上述方程可以把θ分解成两部分I θ和∏θ两部分分别求解，然后运用叠加原理∏+=θθθI 得出最终温度场，一下为分解的I θ和∏θ两部分： （2）首先求解温度场I θ用分离变量法假设所求的温度分布),(y x I θ可以表示成一个x 的函数和一个y 的函数的乘积，即)()(),(11y Y x X y x I =θ （2-2）将上式代入I θ的导热微分方程中，得到012121212=+X dyY d Y dx X d ，即21''11''1ε=-=Y Y X X ，上式等号左边是x 的函数，右边是y 的函数，只有他们都等于一个常数时才可能成立，记这个常数为2ε。
由此得到一个待定常数的两个常微分方程001221212212=+=-Y dyY d X dxX d εε （2-3） 解得)()()(1x Bsh x Ach x X εε+= （2-4） )sin()cos()(1y D y C y Y εε+= （2-5） 把边界条件0,0=∂∂=yy Iθ代入（2-3-4）得到A=0，所以有)()(1x Bsh x X ε= （2-6） 把边界条件0,=∂∂=yL y Iθ代入（2-3-5）得到D=0，所以有 )cos()(1y C y Y ε= （2-7） 把边界条件0,=+∂∂=I Ih yL y θθλ联立（2-3-7）得到 λεε/)cot(hL LL =（2-8）设Bi hL L ==λβε/,，则有i B /)cot(ββ=，这个方程有无穷多个解，即常数β有无穷多个值，即)3,2,1( =n n β，所以对应无穷多个ε，即)3,2,1( =n n ε，所以有)cos()(1y C y Y n n ε= （2-9） 联立（2-3-6）可得∑∞==1)()cos(),(n n n n I x sh y K y x εεθ （2-10）把边界条件2,θθδ==I x 代入上式可得⎰⎰=Ln n n Ln dy y sh K dy y 0202)(cos )()cos(εδεεθ （2-11）解得])cos())[sin(/()sin(22n n n n n n L sh K βββδββθ+= （2-12）其中L n n εβ= )()cos(])cos())[sin(/()sin(2),(12x L sh y L L sh y x n n n n n n n n I βββββδββθθ∑∞=+= （2-13）（3）求解温度场∏θ与解I θ一样用分离变量法，假设所求温度分布),(y x ∏θ可以表示成一个x 的函数和一个y 的函数的乘积)()(),(22x Y x X y x =∏θ （2-14）将该式子代入∏θ的导热微分方程中得到022222222=+X dyY d Y dx X d ，即22''22''2ε=-=Y Y X X ，由此可得到两个常微分方程 02222=-X dx X d ε （2-15） 022222=+Y dy Y d ε （2-16）解式（2-3-15）时根据x 的边界条件可以把解的形式写为 )]([)]([)(2x Bsh x Ach x X -+-=δεδε （2-17） 把边界条件0,==∏θδx 代入上式，得到A=0，所以有 )]([)(2x Bsh x X -=δε （2-18） 其中i n n n n B L /)cot(,βββε==)]([)cos(),(1x sh y k y x n n n n I -=∑∞=δεεθ （2-19）把边界条件1,0θθ==∏x 代入上式可得⎰⎰-=LLn n nn dy y x sh K dy y 02'1)(cos )]([)cos(εδεεθ （2-20） ])cos())[sin(/()sin(21'n n n n n nL sh K βββδββθ+= （2-21）)]([)cos(])cos())[sin(/()sin(2),(11x L sh y L L sh y x n n n n n n n n -+=∑∞=∏δβββββδββθθ （2-22）（4）最终求得稳态温度场2-5、地热换热器是管中流动的流体与周围土地之间的换热，可应用于热能的储存、地源热泵等工程实际。
一种布置方式是把管子埋设在垂直于地面的钻孔中。
由于管子的长度远大于钻孔的直径，可把管子的散热简化为一个有限长度的线热源。
当运行的时间足够长以后，系统可以达到基本稳定的状态。
设土地是均匀的半无限大介质，线热源单位长度的发热量为ql ，地表面的温度均匀，维持为t0。
使用虚拟热源法求解土地中的稳态温度场。
解：根据题意画出示意图如下：设有限长热源长度为H ，单位长度热源发热量为l q ，电源强度为)(0w dz q l ⨯，设地面温度维持恒定温度00,t t t -=θ。
（1）求解点热源dz0产生的温度场有限长线热源在某点产生的温度可以看做是许多点源在该点产生的温度场的叠加，因此我们先来看下无限大介质中点源产生的温度场，这是一个球坐标系中的无内热源的稳态导热问题，其导热微分方程为：0)(122=drd r dr d r θ（3-1） 解微分方程可得rc c 12-=θ （3-2） 把边界条件0,=∞→θr 代入上式得到02=c ，所以有 rc 1-=θ （3-3） 在球坐标系点热源0dz 单位时间内的发热量等于它在任意球面上产生的热流量Q ，即01244dz q c r drd Q l =-=-=λππθλ （3-4） 所以得到014dz q c lπλ-=由此可得到球坐标系中点热源0dz 产生的温度场为 0*14dz rq l πλθ=（3-5） （2）分别求出两个线热源产生的温度场线热源产生的温度场可以看作是点热源产生的温度场的叠加，因此有 地下有限长线热源产生的温度场0014dz rl ⎰=πλθ （3-6） 对称的虚拟热源产生的温度场为 00214dz rq Hl ⎰--=πλθ （3-7） （3）虚拟热源法求解的地热换热器产生的温度场⎥⎥⎦⎤⎢⎢⎣⎡-++++++++-+-=⎥⎥⎦⎤⎢⎢⎣⎡++--+=-+=⎰⎰⎰-z z z z z H z H z H z H q dz z z z z q dz rq dz r q lH l H lHl 22222222002022020000)()(ln 4)(1)(141414ρρρρπλρρπλπλπλθ （3-8） 第三章3-1、用热电偶测量呈简谐波周期变化的气流温度，热电偶的感温节点可看作直径为1mm 的圆球，其材料的密度为8900kg/m3，比热容为390J/(KgK)，测温记录最高和最低温度分别为130℃和124℃，周期为20s 。
若已知气流与热电偶间的对流换热的表面传热系数为20W/(m2K)，试确定气流的真实温度变化范围。
解：气流温度按简谐波变化时，热电偶的温度响应为 )cos(*ϕτθ+=w B （4-1） 式中)arctan(122r rf w w A B τϕτ-=+=按题目要求102022πππ===T w ，s hA cv r 925.2820610139089003=⨯⨯⨯⨯==-ρτ，)/(202k m w h ⋅=，根据题目提供的热电偶测量的最高温度、最低温度，求出热电偶测量的温度变化的振幅如下式32124130122=-=+r f w A τ （4-2） 把r w τ,的数据代入上式中得到气流温度变化的振幅4.27=f A ，所以真实气体温度变化的最大值、最小值为 C t 0max 4.1544.272124130=++=（4-3）C t 0min 6.994.272=-=（4-4） 3-6、已知初始温度均匀的无限大介质中由连续恒定发热的线热源所引起的温度场由式子)4(4),(2τπλτa r E q r t i l --=确定。
若线热源的加热不是连续的而是间歇的，即从0=τ的时刻起，线热源进行周期性的间歇加热，周期为T ，其中加热的时段为T1，其余的T-T1时间不加热。
试利用线性叠加原理确定介质中的温度响应。
解：无限大介质连续恒定发热的线热源引起的温度场：)4(4),(2τπλτa r E q t r t i i --=-∞ （5-1） 其中：du ue z E zui ⎰∞-=)( 对于随时间变化的热流可以用一系列连续的矩形脉冲热流来近似如图所示：由叠加原理得到τ时刻的温度变化为：)0(,])(4[)(4101121=----=--=∞-∑l i i l l ni q a r E q q t t i i ττπλ （5-2）对于间歇性的脉冲，令T T C l /=为运行份额，如果在整个运行期间的平均热负荷为l q ，则脉冲加热的强度为C q l /，具体见下图：由叠加原理得到：∑∑∑∞=∞=∞=∞⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡---⎥⎦⎤⎢⎣⎡---=--------=-0222020)(4(44])(4[4])(4[4n i l i l l i n l i n l nT a r E T nT a r E C q T nT a r E q nT a r E q t t ττπλτπλτπλ （5-3）即温度响应为∑∞=∞⎭⎬⎫⎩⎨⎧⎥⎦⎤⎢⎣⎡---⎥⎦⎤⎢⎣⎡---=-=022)(4(41)(4n i l i l nT a r E T nT a r E C q t t ττπλθ （5-4）第四章4-1、处在x>0的半无限大空间内的一固体，初始温度为溶解温度t m 。
