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 t T c d drd r q d q d q dr r q p r sin sin 2 2??=+??-??-??-? （1-2） 导热速率可根据傅里叶定律计算： ?θθd r rd t T k q r r sin ???-= ?θθθθd r dr T r k q sin ???-= （1-3） θ? θ? ?rd dr T r k q ???- =sin 将上述式子代入（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 对于各向异性材料，化简整理后可得到： t T c q T r k T r k r T r r r k p r ??=+??+????+?????ρ?θθθθθ?θ2 222222sin )(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,,θθθ，根据题意写出下列方程组 ????? ??? ?? ?=+??==??======??+??00 000212222θθ λθθθδθθθ θh y L y y y x x y x （2-1） 解上述方程可以把θ分解成两部分I θ和∏θ两部分分别求解，然后运用叠加原理∏+=θθθI 得出最终温度场，一下为分解的I θ和∏θ两部分：
高等传热学导热理论——相变导热（移动边界问题）讨论 第五讲：相变导热（移动边界问题）： 移动边界的导热问题有许多种，本讲只讲固液相变时的导热模型。 5.1 相变换热特点与分类： 特点： (1) 相变处存在一个界面把不同相的物质分成两个区间（实际不是一个面， 而是一个区）。 (2) 相变面随时间移动，移动规律时问题的一部分。 (3) 移动面可作为边界，决定了相变问题是非线性问题。 分类： (1) 半无限大体单区域问题（Stefan Question ） (2) 半无限大体双区域问题（Neumman Question ） (3) 有限双区域问题 5.2 相变导热的数学描述和解： 假定：固液两相内部只有导热，没有对流（适用于深空中相变）。 物性为常量。不考虑密度变化引起的体积变化。 控制方程： 对固相： 2 21s s s t t a x τ ??=?? 对液相： 2 2 1l l l t t a x τ ??= ?? 初值条件：0:s l t t t τ∞=== 边界条件： 0:::s l w l s l s x t ort t x t ort or x t ort t ∞ ===∞≠∞ =?= 在相变界面，热量守恒，温度连续，Q l 为相变潜热： ()():s l s l l l s l p t t d x Q and t t t x x d δτδτλλρτ ??==+==?? 5.2.1 半无限大体单区域问题（Stefan Question ）的简化解： 以融解过程为例： 忽略液相显热， 2 210l l l t t a x τ ??==??，方程解为一直线，由边界条件得： ()/l w p w t t t t x δ =+- 对固相，忽略温差：w p t t t ∞==，即固相温度恒等于相变温度等于初始温度。 由相变处得换热条件求δ的变化规律：
沈阳航空航天大学 预测燃气涡轮燃烧室出口温度场 沈阳航空航天大学 2013年6月28日
计算传热学 图1模型结构和尺寸图 1.传热过程简述 计算任务是用计算流体力学/计算传热学软件Fluent求解通有烟气的法兰弯管包括管内烟气流体和管壁固体在内的温度分布，其中管壁分别采用薄壁和实体壁两种方法处理。在进行分析时要同时考虑导热、对流、辐射三种传热方式。 (1) 直角弯管内外壁面间的热传导。注意：如果壁面按薄壁处理时，则不用考虑此项，因为此时管壁厚度忽略不计，内壁和外壁温度相差几乎为零。 (2) 管道外壁面与外界环境发生的自然对流换热。由于流体浮生力与粘性力对自然对流的影响，横管与竖管对流换热系数略有不同的。计算公式也不一样。同时，管道内壁面同烟气发生的强制对流换热。 (3) 管道外壁和大空间（环境）发生辐射换热 通过烟气温度和流量，我们可以推断出管道内烟气为湍流流动。这在随后的模
沈阳航空航天大学 拟计算中可以得到证实。 2.计算方案分析 2.1 控制方程及简化 2.1.1质量守恒方程： 任何流动问题都要满足质量守恒方程，即连续方程。其积分形式为： 0vol A dxdydz dA t ρρ?+=?????? 式中，vol 表示控制体；A 表示控制面。第一项表示控制体内部质量的增量，第二项表示通 过控制面的净通量。 直角坐标系中的微分形式如下： ()()()0u v w t x y z ρρρρ????+++=???? 上式表示单位时间内流体微元体中质量的增加，等于同一时间段内流入该微元体的净增量。 对于定常不可压缩流动，密度ρ为常数，该方程可简化为 0u v w x y z ???++=??? 2.1.2动量守恒方程： 动量守恒方程也是任何流动系数都必须满足的基本定律。数学式表示为： F m dv dt δδ= 流体的粘性本构方程得到直角坐标系下的动量守恒方程，即N-S 方程： ()()()u u p div Uu div gradu S t x ρρμ??+=+-?? ()()()v v p div Uv div gradv S t y ρρμ??+=+-?? ()()()w w p div Uw div gradw S t z ρρμ??+=+-?? 该方程是依据微元体中的流体的动量对时间的变化率等于外界作用在该微元体上的各种力之和。式中u S 、v S 、w S 是动量方程中的广义源项。和前面方程一样上式
第2章边界层方程 第一节Prandtl 边界层方程一.边界层简化的基本依据 外：粘性和换热可忽略 )(t δδ , l l t <<<<δδ或内：粘性和换热存在 )(t δδ特征尺寸 —l
二.普朗特边界层方程 常数性流体纵掠平板，层流的曲壁同样适用)。 δ v l u ∞∞ ∞u l v v l u δδ~~，可见，0=??+??y v x u )()((x x R δ>>曲率半径y x u v ∞ ∞T u ,w T ∞ ∞T u ,δ l
)(122 22 y u x u x p y u v x u u ??+??+??-=??+??νρδ δ ∞ ∞ u u l l u u ∞∞ 2 l u ∞ν2 δ ν ∞ u ) (2 l u ∞ 除以无因次化11 Re 12 ) )(Re 1 (δ l
因边界层那粘性项与惯性项均不能忽略，故 项可忽略，且说明只有Re>>1时，上述简化才适用。)(12 2 22y v x v y p y v v x v u ??+??+??-=??+??νρ1~))(Re 1(2 δ l l δ ；可见22 22 x u y u ??>>??δδ 1 ) (2 ∞u l l u l u /)(∞∞δ 2 /)(l u l ∞δ ν2 /)(δδ ν∞u l ： 除以l u 2 ∞ )(Re 1l δ))(Re 1(δ l l δ
可见，各项均比u 方程对应项小得多可简化为 于是u 方程压力梯度项可写为。 )(2 2 22y T x T a y T v x T u ??+??=??+??,0=??y p dx dp ρ1-),(l δ 乘了δθδ w u l )(∞l u w θ∞2 l a w θ除以： l u w θ∞Pe /12 )(/1δ l Pe 12δ θw a 1 ) (∞-=T T w w θPr) Re (?====∞∞贝克列数—导热量对流热量w w p l k u c a l u Pe θθρ
二维导热物体温度场的数值模拟
一、物理问题 有一个用砖砌成的长方形截面的冷空气通道， 于纸面方向上用冷空气及砖墙的温度变化很小， 可以近似地予以忽略。 在下列两种情况下试计算： 砖墙横截面上的温度分布；垂直于纸面方向的每 米长度上通过砖墙的导热量。 第一种情况：内外壁分别均匀维持在 0℃及 30℃； 第二种情况：内外壁均为第三类边界条 件， 且已知： t 1 30 C,h 1 10.35W / m 2 K 2 t 2 10 C, h 2 3.93W / m 2 K 砖墙导热系数 0.35/ m K 二、数学描写 由对称的界面必是绝热面， 态、无内热源的导热问题。 控制方程： 22 tt 22 xy 边界条件： 第一种情况： 由对称性知边界 1 绝热： 边界 2 为等温边界，满足第一类边界条件： t w 0 C ； 边界 3 为等温边界，满足第一类边界条件： t w 30 C 。 第一种情况： 由对称性知边界 1 绝热： q w 0； 边界 2 为对流边界，满足第三类边界条件： q w ( t )w h 2(t w 可取左上方的四分之一墙角为研究对象， 该问题为二维、 稳 图1-
t f )； n t 边界3 为对流边界，满足第三类边界条件：q w ( ) w h 2 (t w t f )。 w n w 2 w f
0,m 6,n 1~ 7;m 7 ~ 16,n 7 30,m 1,n 1~12;m 2 ~ 16,n 12 三、方程离散 用一系列与坐标轴平行的间隔 0.1m 的二维网格线 将温度区域划分为若干子区域，如图 1-3 所示。 采用热平衡法， 利用傅里叶导热定律和能量守恒定 律，按照以导入元体（ m,n ）方向的热流量为正，列写 每个节点代表的元体的代数方程， 第一种情况： 边界点： 1 边界 绝热边界） ： 边界 图1-3 t m ,1 t 16,n 等温内边界） ： 14 (2t m,2 1 4 (2t 15,n t m 1,1 t m 1,1)，m 2 ~ 5 t 16,n 1 t 16,n 1)， n 8 ~ 11 边界 等温外边界） ： 内节 点： 1 (t t t t ) 4 m 1,n m 1,n m ,n 1 m,n 1 m 2 ~ 5,n 2 ~11;m 6 ~ 15,n 8 ~ 11 t m,n 第二种情况 边界点： 边界 1（绝热边界） ： t m ,1 1 4 (2t m,2 t m 1,1 t m 1,1)，m 2 ~ 5 t 16,n 1 4 (2t 15,n t 16,n 1 t 16,n 1)， n 8 ~11 4 边界 2（内对流边界） ： t6,n 2t 5,n t 6,n 1 t 6,n 1 2Bi 1t 1 ，n 1~ 6 6,n 2(Bi 2) t m,n t m,n
高等传热学知识重点 1.什么是粒子的平均自由程，Knusen数的表达式和物理意义。 Knusen数的表达式和物理意义：（Λ即为λ，L为特征长度） 2.固体中的微观热载流子的种类，以及对金属/绝缘体材料中热流的贡献。 3.分子、声子和电子分别满足怎样的统计分布律，分别写出其分布函数的表达式 分子的统计分布：Maxwell-Boltzmann（麦克斯韦-玻尔兹曼）分布： 电子的统计分布：Fermi-Dirac（费米-狄拉克）分布： 声子的统计分布：Bose-Eisentein（波色-爱因斯坦）分布； 高温下，FD,BE均化为MB；
4.什么是光学声子和声学声子，其波矢或频谱分布各有特性？ 答：声子：晶格振动能量的量子化描述，是准粒子，有能量，无质量； 光学声子：与光子相互振动，发生散射，故称光学声子； 声学声子：类似机械波传动，故称声学声子； 5.影响声子和电子导热的散射效应有哪些？ 答：影响声子（和电子）导热的散射效应有（热阻形成的主要原因）： ①界面散射：由于不同材料的声子色散关系不一样，即使是完全结合的界面也是有热阻的； ②缺陷散射：除了晶格缺陷，最典型的是不纯物掺杂颗粒的散热，散射位相函数一般为Rayleigh散 射、Mie散射，这与光子非常相似； ③声子自身散射：声子本质上是晶格振动波，因此在传播过程中会与原子相互作用，会产生散射、 吸收和变频作用。
6.简述声子态密度（Density of State）及其物理意义,德拜模型和爱因斯坦模型的区别。答：声子态密度（DOS）[phonon.s/m3.rad]：声子在单位频率间隔内的状态数（振动模式数）Debye（德拜）模型： Einstein（爱因斯坦）模型： 7.分子动力学理论中，L-J势能函数的表达式及其意义。 答：Lennard-Jones 势能函数（兰纳-琼斯势能函数）,只适用于惰性气体、简单分子晶体，是一种合理的近似公式；式中第一项可认为是对应于两体在近距离时以互相排斥为主的作用，第二项对应两体在远距离以互相吸引（例如通过范德瓦耳斯力）为主的作用，而此六次方项也的确可以使用以电子-原子核的电偶极矩摄动展开得到。
硕士研究生《高等工程热力学与传热学》作业 查阅相关资料,回答以下问题： 1、一滴水滴到120度和400度的板上，哪个先干？试从传热学的角度分析？ 答：在大气压下发生沸腾换热时,上述两滴水的过热度分别是△ t=tw–ts=20℃和△t=300℃,由大容器饱和沸腾曲线,前者表面发生的是泡态沸腾,后者发生膜态沸腾。虽然前者传热温差小,但其表面传热系数大,从而表面热流反而大于后者。所以水滴滴在120℃的铁板上先被烧干。 2、锅铲、汤勺、漏勺、铝锅等炊具的柄用木料制成，为什么？ 答：是因为木料是热的不良导体,以便在烹任过程中不烫手。 3、滚烫的砂锅放在湿地上易破裂。为什么？ 答：这是因为砂锅是热的不良导体, 如果把烧得滚热的砂锅,突然放到潮湿或冷的地方,砂锅外壁的热就很快地被传掉，而内壁的热又一下子传不出来，外壁冷却很快的收缩，内壁却还很热，没什么收缩，加以陶瓷特别脆，所以往往裂开。 或者：烫砂锅放在湿地上时,砂锅外壁迅速放热收缩而内壁温度降低慢,砂锅内外收缩不均匀,故易破裂。 4、往保温瓶灌开水时，不灌满能更好地保温。为什么？ 答：因为未灌满时,瓶口有一层空气,是热的不良导体,能更好地防止热量散失。
5、煮熟后滚烫的鸡蛋放入冷水中浸一会儿，容易剥壳。为什么？ 答：因为滚烫的鸡蛋壳与蛋白遇冷会收缩,但它们收缩的程度不一样,从而使两者脱离。 6、用焊锡的铁壶烧水，壶烧不坏，若不装水，把它放在火上一会儿就烧坏了。为什么？ 答：这是因为水的沸点在1标准大气压下是100℃,锡的熔点是232℃,装水烧时,只要水不干,壶的温度不会明显超过100℃,达不到锡的熔点,更达不到铁的熔点,故壶烧不坏.若不装水在火上烧,不一会儿壶的温度就会达到锡的熔点,焊锡熔化,壶就烧坏了。 7、冬天水壶里的水烧开后，在离壶嘴一定距离才能看见“白气”，而紧靠壶嘴的地方看不见“白气”。这是因为紧靠壶嘴的地方温度高，壶嘴出来的水蒸气不能液化，而距壶嘴一定距离的地方温度低；壶嘴出来的水蒸气放热液化成小水滴，即“白气”。 答：这是因为紧靠壶嘴的地方温度高,壶嘴出来的水蒸气不能液化,而距壶嘴一定距离的地方温度低；壶嘴出来的水蒸气放热液化成小水滴,即“白气”。 8、某些表演者赤脚踩过炽热的木炭，从传热学角度解释为何不会烫伤？不会烫伤的基本条件是什么？ 答：因为热量的传递和温度的升高需要一个过程，而表演者赤脚接触炽热木炭的时间极短，因此在这个极短的时间内传递的温度有限，不足以达到令人烫伤的温度，所以不会烫伤。 基本条件：表演者接触炽热木炭的时间必须极短，以至于在这段时间内所传递的热量不至于达到灼伤人的温度
高等传热学复习题 1．简述求解导热问题的各种方法和傅立叶定律的适用条件。 不论如何，求解导热微分方程主要依靠三大方法： 理论法、试验法、综合理论和试验法 理论法：借助数学、逻辑等手段，根据物理规律，找出答案。它又分： 分析法；以数学分析为基础，通过符号和数值运算，得到结果。方法有：分离变量法，积分变换法(Laplace变换，Fourier变换)，热源函数法，Green函数法，变分法，积分方程法等等，数理方程中有介绍。 近似分析法：积分方程法，相似分析法，变分法等。 分析法的优点是理论严谨，结论可靠，省钱省力，结论通用性好，便于分析和应用。缺点是可求解的对象不多，大部分要求几何形状规则，边界条件简单，线性问题。有的解结构复杂，应用有难度，对人员专业水平要求高。 数值法：是当前发展的主流，发展了大量的商业软件。方法有：有限差分法，有限元法，边界元法，直接模拟法，离散化法，蒙特卡罗法，格子气法等，大大扩展了导热微分方程的实用范围，不受形状等限制，省钱省力，在依靠计算机条件下，计算速度和计算质量、范围不断提高，有无穷的发展潜力，能求解部分非线性问题。缺点是结果可靠性差，对使用人员要求高，有的结果不直观，所求结果通用性差。 比拟法：有热电模拟，光模拟等 试验法：在许多情况下，理论并不能解决问题，或不能完全解决问题，或不能完美解决问题，必须通过试验。试验的可靠性高，结果直观，问题的针对性强，可以发掘理论没有涉及的新规律。可以起到检验理论分析和数值计算结果的作用。理论越是高度发展，试验法的作用就越强。理论永远代替不了试验。但试验耗时费力，绝大多数要求较高的财力和投入，在理论可以解决问题的地方，应尽量用理论方法。试验法也有各种类型：如探索性试验，验证性试验，比拟性试验等等。 综合法：用理论指导试验，以试验促进理论，是科学研究常用的方法。如浙大提出计算机辅助试验法(CA T)就是其中之一。 傅里叶定律向量形式说明，热流密度方向与温度梯度方向相反。它可适用于稳态、非稳态，变导热系数，各向同性，多维空间，连续光滑介质，气、液、固三相的导热问题。 2．定性地分析固体导热系数和温度变化的关系 3．什么是直肋的最佳形状与已知形状后的最佳尺寸？ Schmidt假定：如要得到在给定传热量下要求具有最小体积或最小质量的肋的形状和尺寸，肋片任一导热截面的热流密度都应相等。 1928年，Schmidt等提出了一维肋片换热优化理论：设导热系数为常数，沿肋高的温度分布应为一条直线。Duffin应用变分法证明了Schmidt假定。Wikins[3]指出只有在导热系数和换热系数为常数时，肋片的温度分布才是线性的。Liu和Wikins[4]等人还得到了有内热源及辐射换热时优化解。长期以来肋片的优化问题受到理论和应用两方面的重视。 对称直肋最优型线和尺寸的无量纲表达式分析： 假定一维肋片，导热系数和换热系数为常数，我们有对称直肋微分方程（忽略曲 线弧度）： yd2θ/dx2+(dy/dx)dθ/dx-θh/λ=0 由Schmidt假定，对任意截面x： dθ/dx＝－q/λ＝const
高等传热学作业修订版 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 t T c d drd r q d q d q dr r q p r sin sin 2 2??=+??-??-??-? （1-2） 导热速率可根据傅里叶定律计算： ?θθ θθd r dr T r 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 对于各 向异性材料，化简整理后可得到： t T c q T r k T r k r T r r r k p r ??=+??+????+?????ρ?θθθθθ?θ2222222sin )(sin sin )( （1-6）
7-4，常物性流体在两无限大平行平板之间作稳态层流流动，下板静止不动，上板在外力作用下以恒定速度U 运动，试推导连续性方程和动量方程。 解：按照题意 0, 0=??=??=x v y v v 故连续性方程 0=??+??y v x u 可简化为 0=??x u 因流体是常物性，不可压缩的，N-S 方程为 x 方向： )(12222y u x u v y p F y u v x u u x ??+??+??-=??+??ρρ 可简化为 022=??+??-y v x p F x η y 方向 )(12222y v x v v y p F y v v x v u y ??+??+??-=??+??ρρ 可简化为 0=??= y p F y 8-3，试证明，流体外掠平壁层流边界层换热的局部努赛尔特数为 12121 Re Pr x Nu r = 证明：适用于外掠平板的层流边界层的能量方程
22t t t u v a x y y ???+=??? 常壁温边界条件为 0w y t t y ∞ ==→∞时，时，t=t 引入量纲一的温度w w t t t t ∞-Θ= - 则上述能量方程变为22u v a x y y ?Θ?Θ?Θ+=??? 引入相似变量1Re ()y y x x ηδ= == 有 11()(()22x x x ηη ηηη?Θ?Θ?''==Θ-=-Θ??? ()y y ηηη?Θ?Θ?'==???；22()U y x ηυ∞ ?Θ''= Θ? 将上三式和流函数表示的速度代入边界层能量方程，得到 1 Pr 02 f '''Θ+Θ= 当Pr 1时，速度边界层厚度远小于温度边界层厚度，可近似认为温度边界层内 速度为主流速度，即1,f f η'==，则由上式可得 Pr ()2d f d η''Θ'=-'Θ，求解可得 12 12 ()()Pr 2 Pr (0)()erf η ηπ Θ='Θ= 则1212 0.564Re Pr x x Nu = 8-4，求证，常物性不可压缩流体，对于层流边界层的二维滞止流动，其局部努
2011年《高等传热学》结课作业 ———放假前提交作业 一、【15分】无内热源物体内的稳态导热，材料为常物性。请选择合适的坐标系，写出其导 热微分方程及边界条件。 (1) 巨型薄板（0≤x≤L1，0≤y≤L2，0≤z≤L3），L3<
高等传热学作业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 t T c d drd r q d q d q dr r q p r sin sin 2 2??=+??-??-??-? （1-2） 导热速率可根据傅里叶定律计算： ?θθ θθd r dr T r 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 对于各向异性材料，化简整理后可得到： t T c q T r k T r k r T r r r k p r ??=+??+????+?????ρ?θθθθθ?θ2 222222sin )(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,,θθθ，根据题意写出下列方程组
东南大学能源与环境学院 课程作业报告 作业名称：传热学大作业——利用matlab程序解决热传导问题 院系：能源与环境学院 专业：建筑环境与设备工程 学号： 姓名： 2014年11月9日
一、题目及要求 1.原始题目及要求 2.各节点的离散化的代数方程 3.源程序 4.不同初值时的收敛快慢 5.上下边界的热流量（λ=1W/(m℃）） 6.计算结果的等温线图 7.计算小结 题目：已知条件如下图所示： 二、各节点的离散化的代数方程 各温度节点的代数方程 ta=(300+b+e)/4 ; tb=(200+a+c+f)/4; tc=(200+b+d+g)/4; td=(2*c+200+h)/4 te=(100+a+f+i)/4; tf=(b+e+g+j)/4; tg=(c+f+h+k)/4 ; th=(2*g+d+l)/4 ti=(100+e+m+j)/4; tj=(f+i+k+n)/4; tk=(g+j+l+o)/4; tl=(2*k+h+q)/4
tm=(2*i+300+n)/24; tn=(2*j+m+p+200)/24; to=(2*k+p+n+200)/24; tp=(l+o+100)/12 三、源程序 【G-S迭代程序】 【方法一】 函数文件为： function [y,n]=gauseidel(A,b,x0,eps) D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1); G=(D-L)\U; f=(D-L)\b; y=G*x0+f; n=1; while norm(y-x0)>=eps x0=y; y=G*x0+f; n=n+1; end 命令文件为： A=[4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0,0; -1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0,0; 0,-1,4,-1,0,0,-1,0,0,0,0,0,0,0,0,0;
高等传热学导热理论 参考书：高等传热学 贾力 方肇洪 钱兴华 ?S .K a k a c ,Y .Y e n e r , H e a t C o n d u c t i o n 1985, T K 124/Y K 3 ?G .E .M y e r s , A n a l y t i c a l M e t h o d s i n C o n d u c t i o n H e a t T r a n s f e r ,1971,T K 124/Y M 1 ?M .N .O z i s i k ,H e a t C o n d u c t i o n ,1980,(中译本)O 551.3/A 2 ?俞昌铭，热传导及数值分析，1981，清华大学出版社， O 551.3/Y 2 ?J .E .P a r r o t t ,A .D .S t u c k e s ,T h e r m a l C o n d u c t i o n o f S o l i d s ,1975, O 551.3/Y P 1 ?U .G r i g u l l ,H .S a n d n e r , ,H e a t C o n d u c t i o n ,1984,Y K 124/Y G 3 ?E c k e r t E .R .G ,A n a l y s i s o f H e a t a n d M a s s T r a n s f e r , O 551.3/Y E 1(英)， O 551.3/A 3，（中） ?V .C .A r p a c i ,C o n d u c t i o n H e a t T r a n s f e r ,1966, ?钱壬章等，传热分析与计算，高教出版社 ?林瑞泰，热传导理论与方法，天津大学出版社 ?屠传经等，热传导，浙江大学出版社 第一讲 导热规律及其数学描述 导热可发生在物体的各种状态：气态、固态和液态。描述传热规律最基本的规律是傅里叶导热定律： 1. F o u r i e r L a w : dx dt q λ-= 傅里叶定律适用于稳态、非稳态，变导热系数，各向同性，多维空间，连续光滑介质，气、液、固三相的导热问题，但其表现形式上为已知热流方向的一维问题。用起来不方便。在已知温度场的情况，我们把傅里叶定律推广成向量形式： n n t t q ??-=?-=λλ 其中?叫n a b l a 算子，作用于温度叫温度梯度。n 为温度梯度单位方向向量。在 不同的坐标系中，?有不同的表现形式，在直角坐标系中： k z j y i x ??+??+??=? 傅里叶定律向量形式说明，热流密度方向与温度梯度方向相反。它可适用于稳态、非稳态，变导热系数，各向同性，多维空间，连续光滑介质，气、液、固三相的导热问题。 2．各向异性材料，导热系数张量； 许多物体的导热能力与方向有关，如木材。正确描述物体中一点的导热系数需采用二阶张量形式：
高等传热学问题及答案 1. 简述三种基本传热方式的传热机理并用公式表达传热定律；传热问题的边界条件有哪两类？ 2. 有限元法求解传热问题的基本思想是什么？基本求解步骤有哪些？同有限差分方法相比其优点是什么？ 3. 什么是形函数？形函数的两个最基本特征是什么？ 4. 加权余量法是建立有限元代数方程的基本方法，请描述四种常见形式并用公式表达。 5. 特征伽辽金法（CG ）在处理对流换热问题时遇到什么困难？特征分离法（CBS ）处理对流换热问题的基本思想是什么？ 第一题： （1）热传导 传热传导模式是因为从一个分子到另一个分子的能量交换，没有分子的实际运动，如果自由电子存在，也可能因为自由电子的运动。因此，这种形式的热输送在很大程度上取决于介质的性质，如果存在温度差，热传导发生在固体，液体和气体。 书上补充： 当两个物体有温差，或者物体内部有温度差时，在物体各部分之间不发生相对位移的情况下，物体微粒（分子，原子或自由电子）的热运动传递了热量。 （2）热对流 ()a w T T h q -=（牛顿冷却定律） 存在于液体和气体中的分子具有运动的自由，它们随身携带的能量（热量），从热区域移动到冷区域。由于在液体或气体的宏观运动，热量传递从一个地区到另一个地方 ,加上流体内的热传导能量传递,称为对流换热。对流可能是自然对流、强制对流，或混合对流。 百度补充： 对流仅发生于流体中，它是指由于流体的宏观运动使流体各部分之间发生相对位移而导致的热量传递过程。由于流体间各部分是相互接触的，除了流体的整体运动所带来的热对流之外，还伴生有由于流体的微观粒子运动造成的热传导。在工程上，常见的是流体流经固体表面时的热量传递过程，称之为对流传热。 （3）辐射 4w T q εσ=（斯蒂藩-玻耳兹曼定律）
第九章 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 )= 2 1 (100+70) ℃=85℃。 查教材附录5，水的物性为：ρ=958.4kg/m 3;λ=0.683 W ／(m 2·℃);μ=282.5×10-6N·s/m 2 假设流态为层流： 4 1 3 2)(13.1? ? ? ???-=w s t t l r g h μλρ 41 6 3 3 2 )70100(2.1105.282102257683.081.94.95813.1?? ????-???????=- W ／(m 2 ·℃) =5677 W ／(m 2·℃) 3 6102257105.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 ?t 61068.330 60 3061312?=??=?Φ?= t mc p τJ/℃ 16．当液体在一定压力下做大容器饱和沸腾时，欲使表面传热系数增加10倍，沸腾温 差应增加几倍?如果同一液体在圆管内充分发展段做单相湍流换热，为使表面传热系数增加10倍，流速应增加多少倍?维持流体流动所消耗的功将增加多少倍?设物性为常数。 解 ①由米洛耶夫公式： { 5 .033.22 25.033.211122.0122.0p t h p t h ?=?= 10)(33.21 212=??=t t h h 所以 69.21033.211 2 ==??t t 即当h 增大10倍时，沸腾温差是原来的2.69倍。 ②如为单相流体对流换热，由D-B 公式可知8 .0m u h ∝，即
能源与机械类研究生用教材 高等传热学 Advanced Heat Transfer 编著：李菊香 南京工业大学能源学院 二○○四年八月
绪论 传热是一种最常见的自然物理现象。当物体内部或物体与物体之间存在温度差异时，就会发生热量从温度较高区域传输到温度较低区域的能量传递过程，这就是通常所说的“传热”，因而由于温差或温度梯度的存在而传递的能量就定义为“热”。 在所有的工程领域，几乎都有传热问题。如动力工程、冶金工程、化学工程、石油工程、核反应堆等等，都有热量传递的问题，另外，在机械、纺织、宇航、电子、农业、环保、生物等领域同样也有许多传热问题。学习传热的目的在于理解传热的基本概念，掌握传热的基本规律，了解传热的计算方法，以及探索传热学中的一些新领域，以解决工程技术中的各种传热问题，在各种传热工程中选择合理的传热方式，优化各种参数，使用合理的设备结构，以达到取得最大的经济效益。 传热问题归纳起来有两大类型：一类是着眼于传热速率的大小及其控制的问题，如采用强化传热以缩小传热设备的尺寸，或削弱传热以减少能量损失。另一类是着眼于温度分布及其控制的问题，如采用何种手段保持物体上的温度分布不变或怎么变化。这两类问题是相互联系的，从它们的本质上说，实际是同一类问题。 有必要将热力学和传热学之间的区别和联系加以说明，热力学研究的是平衡体系，应用热力学定律可以预计一个体系从一种平衡状态转变为另一种平衡状态有多少能量输出或输入，但不能指出这一变化过程需要多长时间，因为在变化过程中体系是不平衡的。而传热学可以在给定的具体条件下指出热量将以多大的速率传播。传热学能以确定能量传播速率的基本定律（经验），补充了热力学第一定律和第二定律。例如，一根灼热的钢棒在一桶水中冷却，热力学可以预算出钢棒与水这一体系最终的平衡温度，但不能告知需要多长时间才能达到这个平衡温度，或者在到达这个平衡状态之前的每一个瞬间钢棒的温度情况，而传热学就可以预计出钢棒和水的温度随着时间的变化关系。传热过程必然遵循热力学第一定律和热力学第二定律，因为这两个热力学定律是自然定律，但传热的过程是一个典型的不可逆过程，还有传热自身的规律。所以也可以说，传热学是对热力学的一个补充， 从学科的角度来看，传热学是工程热物理技术学科的一个分支，工程热物理学科是由热力学、传热传质学、气动热力学、流体力学与燃烧学等组成。 从传热学单独形成一门系统的科学至今，只有两个多世纪的历史。但是随着原子能、宇宙航行等尖端技术的发展，不断地出现新的传热问题，促使传热学得到了迅速的发展。电子计算机和测试技术的发展，更加丰富了研究传热学的手段，使得传热学的研究范围不断扩大，研究方法不断更新，理论分析也不断完善。目前各种传热学的分支学科如“计算传热学”、“纳米尺度传热学”、“微通道流体的传热”、“分子传热学”、“传热优化设计”等也在不断地发展。
Homework 5.2 1. 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: 根据能量方程， 22 1p p T T T T c u c v k r k r x r r r x ρρ???????+=+ ???????? (1) 对于充分发展的管内弹状流，u=0, 且由于 m dT T x dx ?=?= constant ，因此(1)式可化简为： 1m p dT T c V k r dx r r r ρ???? = ????? (2) 解之可得， 2 12(r)ln 4p m V c dT T r C r C k dx ρ= ??++ (3) 根据边界条件， 0, , T(r)=T w T r r r R ?==?= 代入(3)中，得到C 1=0, 2 24p m w V c dT C T R k dx ρ=- ? 平均温度201 ()2R m T T r rdr R ππ=? 28()m w m p dT k T T dx V c R ρ=- 因此，得到
222 22(r)(1)22()4m w m w w m dT VR r r T T T T T T dx R R α=+??-=-+- 22 88q''h()()()222p p m w m w m w m p mc c R V dT k k T 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 由于 m dT T x dx ?=?= constant ，(1)式可化为： 22p T T c u k x y ρ??=?? (2) 根据01()h m u u y dy h =?，且2 2()6()m u y y y u h h =-，代入得到 22226()m m p dT y y T u c k h h dx y ρ?-=? 解得温度分布为 3 4 122 ()()2p m m c u dT y y T y C y C k dx h h ρ= ?-++ (3)
