传热学报告
传热学第4章习题报告
姓名：潘清泉
学号：2011302690013
4-16一厚为2.54cm 的钢板，初始温度为650℃，后置于水中淬火，其表面温度突然下降为93.5℃并保持不变。试用数值方法计算中心温度下降到450℃所需的时间。已知a 1.16 10 m / s = ×?5 2 。建议将平板8 等分，取9 个节点，并把数值计算的结果与按海斯勒计算的结果作比较。
1.区域离散化（划分网格）：将平板16等分，取17个节点，由于对称性，计算时只考虑一侧的10个节点：t(-1) t(0) t(1) t(2) t(3) t(4) t(5) t(6) t(7) t(8)
其中t(0)代表中心温度,
t(8)=93.5℃,
t(-1)≡t1;
2，实际计算如下；
取Fo=0.4602889时
3由图像就行求解：
当t=3.8S时中心温度降到450℃
若取Fo=0.0460289 可得t=3.85S. 可见实际所需时间可约等于3.85S。
4-23 、一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面的温度分布及每米长度上通过壁面的冷量损失：
（1）内外壁分别维持在10℃及30℃
（2）内外壁与流体发生对流换热，且有10 1 = f t ℃，20 1 h = W /(m2 .K)，30 2 = f t ℃，
4 2 h = W /(m2 .K)。
1，区域离散化（划分网格）：将原图划分为30×22小正方形，相应的取31、23个节点；
2.取λ=1W/(M^2·K) 第一问：
第二问：
24.8%，在可以接受的范围内。
传热实验实验报告
传热实验 一、实验目的 1、了解换热器的结结构及用途。 2、学习换热器的操作方法。 3、了解传热系数的测定方法。 4、测定所给换热器的传热系数K。 5、学习应用传热学的概念和原理去分析和强化传热过程，并实验之。 二、实验原理 根据传热方程Q=KA△tm，只要测得传热速率Q，冷热流体进出口温度和传热面积A，即可算出传热系数K。在该实验中，利用加热空气和自来水通过列管式换热器来测定K,只要测出空气的进出口温度、自来水进出口温度以及水和空气的流量即可。 在工作过程中，如不考虑热量损失，则加热空气释放出的热量Q1与自来水得到的热量Q2应相等，但实际上因热损失的存在，此两热量不等，实验中以Q2为准。 三、实验流程和设备 实验装置由列管换热器、风机、空气电加热器、管路、转子流量计、温度计等组成。空气走管程，水走壳程。列管式换热器的传热面积由管径、管数和管长进行计算。 实验流程图： 四、实验步骤及操作要领 1、熟悉设备流程，掌握各阀门、转子流量计和温度计的作用。 2、实验开始时，先开水路，再开气路，最后再开加热器。 3、控制所需的气体和水的流量。 4、待系统稳定后，记录水的流量、进出口温度，记录空气的流量和进出口温度，记录设备的
有关参数。重复一次。 5、保持空气的流量不变，改变自来水的流量，重复第四步。 6、保持第4步水的流量，改变空气的流量，重复第四步。 7、实验结束后，关闭加热器、风机和自来水阀门。 五、实验数据记录和整理 1、设备参数和有关常数 换热流型错流；换热面积㎡
六、实验结果及讨论 1、求出换热器在不同操作条件下的传热系数。 计算数据如上表，以第一次记录数据序号1为例计算说明： 2、对比不同操作条件下的传热系数，分析数值，你可得出什么结论？ 答：比较一、二、三组可知当空气流量不变，水的流量改变时，传热系数变化不大，比较四、五组可知空气流量改变而水的流量不改变时，传热系数有很大变化，且空气流量越大，传热系数越大，传热效果越好；综上可知，K值总是接近热阻大的流体侧的α值，实验中，提高空气侧的α值以提高K值。。 3、转子流量计在使用时应注意什么问题？应如何校正读数？ 答：转子流量计不能用于流量过大的流体测量，使用时流量计必须安装在垂直走向的管段上，流体介质自下而上地通过转子流量计。 读数时应读转子的最大截面与玻璃管刻线相交处的数值，可以读初始值和最终值，取两者之差来校正读数。 4、针对该系统，如何强化传热过程才能更有效，为什么？ 答：该系统传热效果主要取决于热流体，所以可以通过增加空气流量，提高其所占比例来强化传热效果；减小水的流量；内管加入填充物或采用螺纹管，加热面在上，制冷面在下。因为由实验可知提高热阻大的流体的传热系数可以更有效的强化传热过程。 5、逆流换热和并流换热有什么区别？你能用实验装置加以验证吗？ 答:①逆流换热时热流体是冷热流体流动方向相反；而并流传热时，其冷热流体流动方向相同；②在相同操作条件下，逆流换热器比并流换热器所需传热面积小。可以改变冷热流体进出口方向，测得在相同传热效果下，逆并流所需传热面积大小，从而加以验证。 6、传热过程中，哪些工程因素可以调动？ t ；④换热过程的流型（并流，逆答：①增大传热面积S；②提高传热系数α；③提高平均温差 m 流，错流）。 7、该实验的稳定性受哪些因素的影响？ 答：①冷凝水流通不畅，不能及时排走；②空气成分不稳定，导致被冷凝效果不稳定；③冷热流体流量不稳定；④传热器管表面的相对粗糙度。 8、你能否对此实验装置作些改进，使之能够用于空气一侧对流传热系数的测定？ 答：让空气走壳程，水走管程，根据流体在管外的强制对流公式，可提出空气一侧的对流传热系数α值。
传热学数值计算大作业2014011673
数值计算大作业 一、用数值方法求解尺度为100mm×100mm 的二维矩形物体的稳态导热问题。物体的导热系数λ为1.0w/m·K。边界条件分别为： 1、上壁恒热流q=1000w/m2; 2、下壁温度t1=100℃； 3、右侧壁温度t2=0℃; 4、左侧壁与流体对流换热，流体温度tf=0℃，表面传热系数 h 分别为1w/m2·K、10 w/m2·K、100w/m2·K 和1000 w/m2·K; 要求： 1、写出问题的数学描述； 2、写出内部节点和边界节点的差分方程； 3、给出求解方法； 4、编写计算程序(自选程序语言)； 5、画出4个工况下的温度分布图及左、右、下三个边界的热流密度分布图； 6、就一个工况下（自选）对不同网格数下的计算结果进行讨论； 7、就一个工况下（自选）分别采用高斯迭代、高斯——赛德尔迭代及松弛法（亚松弛和超松弛）求解的收敛性（cpu 时间，迭代次数）进行讨论； 8、对4个不同表面传热系数的计算结果进行分析和讨论。 9、自选一种商业软件（fluent 、ansys 等）对问题进行分析，并与自己编程计算结果进行比较验证（一个工况）。（自选项） 1、写出问题的数学描述 设H=0.1m 微分方程 22220t t x y ??+=?? x=0，0 y=H ，0 《传热学C》课程教学大纲 课程编号：80210127 课程名称：传热学C 英文名称：Heat Transfer C 总学时：24 学分：1.5 适用对象：机械工程及其自动化专业，测控技术及仪器专业 先修课程：高等数学，流体力学 一、课程性质、目的和任务 传热学C是机械工程及其自动化专业和测控技术及仪器专业的一门专业选修课程。其目的在于使学生掌握有关热量传递的基本理论知识，具备一定的传热学分析计算能力。它不仅为以后专业课的学习提供必要的理论基础，也是培养提高学生综合分析能力和解决工程实际问题能力的重要环节之一。 二、教学内容、方法及基本要求 教学内容 1．绪论 了解传热学与工程热力学在研究内容和方法上的异同。认清传热学的研究对象及其在工程和科学技术中的应用。掌握热量传递的基本方式：导热、对流和热辐射的概念和所传递热量的计算公式。了解复合换热过程的计算方法，了解辐射换热表面传热系数的概念。认识到工程实际问题的热量传递过程往往不是单一的方式而是多种形式的组合，以加深传热过程的概念及传热方程的理解。初步理解热阻在分析传热问题中的重要地位。 2．导热基本定律及稳态导热 掌握傅里叶定律的意义和应用方法，了解常见材料导热系数的大致范围。理解推导导热微分方程的理论依据和思路，以及导热微分方程中各项的物理意义，能够正确书写导热问题的初始条件和三类边界条件。能应用傅里叶定律或导热微分方程对常物性、无内热源的一维稳态导热问题（平壁、圆筒壁）进行分析求解，得出温度场及导热量的计算公式。了解肋片在工程中的应用场合。加深理解热阻概念及其在分析导热问题时的重要性。 3．非稳态导热 了解非稳态导热过程的特点。掌握集总参数法的分析求解方法，了解其限制条件。 4．对流换热 牛顿冷却公式是对流换热计算的基础，要求重点掌握。理解影响对流换热的因素。掌握流动边界层和温度边界层的概念。理解相似原理在指导对流换热实验中的作用，准则方程的导出。掌握实验数据的整理方法。掌握管内换热入口段与充分发展段的概念。掌握定型尺寸和定性温度的概念。能正确和熟练地运用准则方程（实验关联式）计算简单的对流换热问题。了解有限空间自然对流换热的概念。掌握强化单相流体对流换热的途径。 5．凝结与沸腾换热 一、 实验名称： 传热实验 二、实验目的： 1.熟悉套管换热器的结构； 2.测定出K 、α，整理出e R N -u 的关系式，求出m A 、. 三、实验原理： 本实验有套管换热器4套，列管式换热器4套，首先介绍套管换热器。 套管换热器管间进饱和蒸汽，冷凝放热以加热管内的空气，实验设备如图2-2-5-1（1）所示。 传热方式为：冷凝—传导—对流 1、传热系数可用下式计算： ]/[2m k m W t A q K m ???= (1) 传热实验 图2-2-5-1（1） 套管换热器示意图 式中：q ——传热速率[W] A ——传热面积[m 2] △t m —传热平均温差[K] ○ 1传热速率q 用下式计算： ])[(12W t t C V q p S -=ρ （2） 式中：3600/h S V V =——空气流量[m 3/s] V h ——空气流量[m 3/h] ρ——空气密度[kg/m 3 ]，以下式计算： ]/)[273(4645.031 m kg t R p P a ++=ρ （3） Pa ——大气压[mmHg] Rp ——空气流量计前表压[mmHg] t 1——空气进换热器前的温度[℃] Cp ——空气比热[K kg J ?/]，查表或用下式计算： ]/[04.01009K kg J t C m p ?+= （4） t m =(t 1+t 2)/2——空气进出换热器温度的平均值（℃） t 2——空气出口温度[℃] ②传热平均面积A m ： ][2m L d A m m π= （5） 式中：d m =传热管平均直径[m] L —传热管有效长度[m ] ③传热平均温度差△t m 用逆流对数平均温差计算： T ←——T t 1——→t 2 )(),(2211t T t t T t -=?-=? 2 1 2 1ln t t t t t m ???-?= ? (6) 式中：T ——蒸汽温度[℃] 2、传热膜系数（给热系数）及其关联式 空气在圆形直管内作强制湍流时的传热膜系数可用下面准数关联式表示： n r m e P AR Nu = （7） 式中：N u ——努塞尔特准数 R e ——雷诺准数 P r ——普兰特准数 A ——系数，经验值为0.023 传热学大作业报告二维稳态计算 院系：能源与环境学院 专业：核工程与核技术 姓名：杨予琪 学号：03311507 一、原始题目及要求 计算要求： 1. 写出各未知温度节点的代数方程 2. 分别给出G-S 迭代和Jacobi 迭代程序 3. 程序中给出两种自动判定收敛的方法 4. 考察三种不同初值时的收敛快慢 5. 上下边界的热流量（λ=1W/(m ℃）） 6. 绘出最终结果的等值线 报告要求： 1. 原始题目及要求 2. 各节点的离散化的代数方程 3. 源程序 4. 不同初值时的收敛快慢 5. 上下边界的热流量（λ=1W/(m ℃）） 6. 计算结果的等温线图 7. 计算小结 二、各节点的离散化的代数方程 左上角节点 )(21 1,22,11,1t t t += 右上角节点 )(2 15,24,15,1t t t += 左下角节点 C t ?=1001,5 右下角节点 )2(211,24,55,5λ λ x h t t x h t ?++?+= 左边界节点 C t i ?=1001,，42≤≤i 上边界节点 C t j ?=200,1，42≤≤j 右边界节点 )2(415,15,14,5,+-++= i i i i t t t t ，42≤≤i 下边界节点 )42()2(211,51,5,4,5∞+-?+++?+=t x h t t t x h t j j j j λλ ，42≤≤j 内部节点 )(2 1,1,11,1,,j i j i j i j i j i t t t t t +-+-+++= ，4,2≤≤j i 三、源程序 1、G-S 迭代法 t=zeros(5,5); t0=zeros(5,5); dteps=0.0001; for i=2:5 %左边界节点 t(i,1)=100; end for j=2:4 %上边界节点 t(1,j)=200; end t(1,1)=(t(1,2)+t(2,1))/2; t for k=1:100 for i=2:4 %内部节点 for j=2:4 t(i,j)=(t(i-1,j)+t(i+1,j)+t(i,j-1)+t(i,j+1))/4; end end t(1,5)=(t(1,4)+t(2,5))/2;%右上角节点 for i=2:4;%右边界节点 t(i,5)=(2*t(i,4)+t(i-1,5)+t(i+1,5))/4; end for j=2:4; %下边界节点 一维稳态导热的数值计算 1.1物理问题 一个等截面直肋，处于温度t ∞=80 的流体中。肋表面与流体之间的对流换热系数为 h =45W/(m 2?℃),肋基处温度t w =300℃，肋端绝热。肋片由铝合金制成，其导热系数为λ=110W/(m ?℃)，肋片厚度为δ=0.01m ，高度为H=0.1m 。试计算肋内的温度分布及肋的总换热量。 1.2数学描述及其解析解 引入无量纲过余温度θ = t?t ∞t w ?t ∞ ，则无量纲温度描述的肋片导热微分方程及其边界条件： 22 20d m dx θθ-= x=0,θ=θw =1 x=H, 0x θ?=? 其中m = 上述数学模型的解析解为：[()] ()() w ch m x H t t t t ch mH ∞∞--=-? ()()w hp t t th mH m ∞?= - 1.3数值离散 1.3.1区域离散 计算区域总节点数取N 。 1.3.2微分方程的离散 对任一借点i 有：22 2 0i d m dx θ θ??-= ??? 用θ在节点i 的二阶差分代替θ在节点i 的二阶导数，得：211 2 20i i i i m x θθθθ+--+-= 整理成迭代形式：()1122 1 2i i i m x θθθ+-=++ (i=2,3……,N-1) 1.3.3边界条件离散 补充方程为：11w θθ== 右边界为第二类边界条件，边界节点N 的向后差分得：1 0N N x θθ--= ，将此式整理为 迭代形式，得：N 1N θθ-= 1.3.4最终离散格式 11w θθ== ()1122 1 2i i i m x θθθ+-= ++ (i=2,3……,N-1) N 1N θθ-= 1.3.5代数方程组的求解及其程序 假定一个温度场的初始发布，给出各节点的温度初值：01θ，02θ，….，0 N θ。将这些初值代 入离散格式方程组进行迭代计算，直至收敛。假设第K 步迭代完成，则K+1次迭代计算式为： K 11w θθ+= () 11 11 2212i i K K K i m x θθθ+-++= ++ (i=2,3……,N-1) 1 11N K K N θθ-++= #include 《传热学》三级项目报告书 钢生锈表面黑度测量 姓名： 课程名称：传热学 指导教师：金昕 2013年6月 固体表面黑度的测定 摘要： 通过课上老师讲解，我们了解到实际物体的辐射力与同温度下黑体的辐射力之比称为实际物体的黑度。为了测量固体表面黑度，首先设计一个已知外表面积的工件，将其放入不存在吸收热辐射介质的空腔内（实验中为真空中的系统），测出换热量，通 过调节试件外表面的温度，研究试件表面黑度随温度的变化，再通过对实验数据的分析整理，计算出黑度，通过图像处理，得到固体表面黑度随温度的变化规律，进一步加深了对黑度的认识。 前言： 本项目目的在于熟悉物体表面黑度的测试原理；熟练应用三维设计软件对实验试件的设计；分析影响物体表面黑度的因素；了解实验原理，并对实验设备进行拆装，学会使用各种测试仪表进行测试；提高学生的动手能力、理论联系实际能力和团队的协作能力；得到查阅文献、阅读相关技术资料和调查研究能力的训炼；通过研究报告的撰写使学生在科技文献写作方面获得训练。 本项目主要内容是测定试件表面黑度以及分析黑度随温度变化关系。黑度是辐射换热的重要特性，黑度取决于物体的性质，物体的温度，表面状态，波长，方向，通过本实验的学习黑度概念和黑度测量，自己动手测绘试件、测量数据、分析结果，增强了动手实验能力，培养灵活运用知识的能力和创新思维。 本实验采用真空辐射法测定固体黑度的实验方法，根据公式及日常经验预测试件表面黑度随温度升高呈下降趋势。 实验目的 1、巩固辐射换热理论. 2、掌握用真空辐射法测定固体表面黑度的试验方法. 3、分析固体表面黑度随温度的变化规律. 实验设备介绍及实验原理分析： 通过查阅相关资料，我们了解到研究黑度的方法主要有：辐射法、量热计法和正规热工况法。辐射法是建立在以被测物体的辐射和绝对黑体或其他辐射系数已知的辐射系数为已知的物体辐射，相比较的基础之上的。本实验采用辐射法。 实验原理 当一物体放在另一物体的空腔内，且空腔内不存在吸收热辐射的介质时（如空气），彼此以辐射换热方式进行热交换，其辐射换热量为： 化工原理实验报告 实验三 传热膜系数测定实验 实验日期：2015年12月30日 班级： 学生姓名： 学号： 同组人： 报告摘要 本实验选用牛顿冷却定律作为对流传热实验的测试原理,通过建立不同体系的传热系统,即水蒸汽—空气传热系统、分别对普通管换热器和强化管换热器进行了强制对流传热实验研究。确定了在相应条件下冷流体对流传热膜系数的关联式。此实验方法可以测出蒸汽冷凝膜系数和管内对流传热系数。采用由风机、孔板流量计、蒸汽发生器等组成的自动化程度较高的装置，让空气走内管，蒸汽走环隙，用计算机在线采集与控制系统测量了孔板压降、进出口温度和两个壁温，计算了传热膜系数α，并通过作图确定了传热膜系数准数关系式中的系数A 和指数m （n 取0.4），得到了半经验关联式。实验还通过在内管中加入混合器的办法强化了传热，并重新测定了α、A 和m 。 二、 目的及任务 1.掌握传热膜系数α及传热系数K 的测定方法； 2.通过实验掌握确定传热膜系数准数关系式中的系数A 和指数m 的方法； 3.了解工程上强化传热的措施。 三、基本原理 对流传热的核心问题是求算传热膜系数α，当流体无相变时对流传热准数关 系式的一般形式为：p n m Gr A Nu Pr Re 对于强制湍流而言。Gr 数可忽略，即 n m A Nu Pr Re = 本实验中，可用图解法和最小二乘法计算上述准数关系式中的指数m 、n 和系数A 。 用图解法对多变量方程进行关联时，要对不同变量Re 和Pr 分别回归。本实验可简化上式，即取n=0.4(流体被加热)。这样，上式即变为单变量方程，在两边取对数，得到直线方程为 Re lg lg Pr lg 4.0m A Nu += 在双对数坐标中作图，求出直线斜率，即为方程的指数m 。在直线上任取一点函数值带入方程中，则可得系数A ，即 m Nu A Re Pr 4.0= 用图解法，根据实验点确定直线位置有一定人为性。而用最小二乘法回归，可得到最佳关联结果。应用计算机辅助手段，对多变量方程进行一次回归，就能的道道A 、m 、n 。 对于方程的关联，首先要有Nu 、Re 、Pr 的数据组。其特征数定义式分别为 μρ du = Re ， λμ Cp = Pr ， λαd Nu = 实验中改变空气的流量，以改变Re 值。根据定性温度(空气进、出口温度的算数平均值)计算对应的Pr 值。同时，由牛顿冷却定律，求出不同流速下的传热膜系数值，进而求得Nu 值。 牛顿冷却定律为 Q=αA △t m 式中α——传热膜系数，W/(m 2.℃)； 《传热学》上机大作业 二维导热物体温度场的数值模拟 学校：西安交通大学 姓名：张晓璐 学号:10031133 班级：能动A06 一．问题（4-23） 有一个用砖砌成的长方形截面的冷空气通道，形状和截面尺寸如下图所示，假设在垂直纸面方向冷空气和砖墙的温度变化很小，差别可以近似的予以忽略。在下列两种情况下计算：砖墙横截面上的温度分布；垂直于纸面方向上的每米长度上通过墙砖上的导热量。 第一种情况：内外壁分别维持在10C ?和30C ? 第二种情况：内外壁与流体发生对流传热，且有C t f ?=101， )/(2021k m W h ?=，C t f ?=302，)/(422k m W h ?=，K m W ?=/53.0λ 二．问题分析 1.控制方程 0222 2=??+??y t x t 2.边界条件 所研究物体关于横轴和纵轴对称，所以只研究四分之一即可，如下图： 对上图所示各边界： 边界1：由对称性可知：此边界绝热，0=w q 。 边界2：情况一：第一类边界条件 C t w ?=10 情况二：第三类边界条件 )()( 11f w w w t t h n t q -=??-=λ 边界3：情况一：第一类边界条件 C t w ?=30 情况二：第三类边界条件 )()( 22f w w w t t h n t q -=??-=λ 三：区域离散化及公式推导 如下图所示，用一系列和坐标抽平行的相互间隔cm 10的网格线将所示区域离散化，每个交点可以看做节点，该节点的温 度近似看做节点所在区域的平均温度。利用热平衡法列出各个节点温度的代数方程。 第一种情况： 内部角点： 11 ~8,15~611 ~2,5~2) (4 1 1,1,,1,1,====++++=+-+-n m n m t t t t t n m n m n m n m n m 平直边界1： 11~8),2(4 1 5~2),2(4 1 1,161,16,15,161,11,12,1,=++==++=+-+-n t t t t m t t t t n n n n m m m m 平直边界2： 传热学上机实验 班级: 学号： 姓名： 一：实验问题 一个长方形截面的冷空气通道的尺寸如附图所示。假设在垂直于纸面的方向上冷空气及通道墙壁的温度变化很小，可以忽略。试用数值方法计算下列两种情况下通道壁面中的温度分布及每米长度上通过壁面的冷量损失： (1)内、外壁面分别维持在10℃及30℃； (2)内、外壁面与流体发生对流传热，且有λ=0.53W/（m·K），t f1=10°C、h1=20W/（m2·K）, t f2=30°C、h2=4W/（m2·K）。 二：问题分析与求解 本题采用数值解法，将长方形截面离散成31×23个点，用有限个离散点的值的集合来代替整个截面上温度的分布，通过求解按傅里叶导热定律、牛顿冷却公式及热平衡法建立的代数方程，来获得整个长方形截面的温度分布，进而求出其通过壁面的冷量损失。 1. 建立控制方程及定解条件 对于第一问，其给出了边界上的温度，属于第一类边界条件。 ????? ??? ??=?==??+??C C y t x t 301002222外壁温内壁温 对于第二问，其给出了边界上的边界上物体与周围流体间的表面传热系数h 及周围流体的温度 t f ，属于第三类边界条件。 ()?????? ?-=??? ????-=??+??f w w t t h n t y t x t λ02222 2. 确定节点（区域离散化） 用一系列与坐标轴平行的网格线把长方形截面划分为31×23个节点。则步长为0.1m ，记为△x=△y=0.1m 。 3. 建立节点物理量的代数方程 对于第一问有如下离散方程： ()()()()()()()()()()? ??? ???? ? ????? ???+++==?==?==?==?==?==?==?==?=+-+-代表内部点,，点41 26~6,1018,26~6,106,18~6,10,2618~6,10,631~1,3023,31~1,301,23~1,30,3123~1,30,11,1,,1,1,n m t t t t t n C m t n C m t n C n t n C n t n C m t n C m t n C n t n C n t n m n m n m n m n m 对于第二问有如下离散方程： 对于外部角点（1,1）、（1,23）、（31,1）、（31,，23）有： ()()02 222,1,,22,,1,22 =??-+-?+??-+-?±±x y t t t t x h y x t t t t y h n m n m n m f n m n m n m f λλ 得到： ()()()()????? ??? ?? ? ++ =++=++=++=22,3123,3023,312,311,301,3122 ,123,223,12,11,21,11865331400186533140018653314001865331400t t t t t t t t t t t t 同理可得： 对于内部角点（6,6）（6,18）（26,6）（26,18） ，有 ()() ()()()()()()????? ??? ??? ++++ =++++ =++++=++++=7,2618,2518,2719,2618,267,266,256,275,266,2618 ,717,619,618,518,67,66,75,66,56,671853359533592000718533595335920007185335953359200071853359533592000t t t t t t t t t t t t t t t t t t t t 传热膜系数测定实验（第四组） 一、实验目的 1、了解套管换热器的结构和壁温的测量方法 2、了解影响给热系数的因素和强化传热的途径 3、体会计算机采集与控制软件对提高实验效率的作用 4、学会给热系数的实验测定和数据处理方法 二、实验内容 1、测定空气在圆管内作强制湍流时的给热系数α1 2、测定加入静态混合器后空气的强制湍流给热系数α1’ 3、回归α1和α1’联式4 .0Pr Re ??=a A Nu 中的参数A 、a *4、测定两个条件下铜管内空气的能量损失 二、实验原理 间壁式传热过程是由热流体对固体壁面的对流传热，固体壁面的热传导和固体壁面对冷流体的对流传热三个传热过程所组成。由于过程复杂，影响因素多，机理不清楚，所以采用量纲分析法来确定给热系数。 1）寻找影响因素 物性：ρ，μ ，λ，c p 设备特征尺寸：l 操作：u ，βgΔT 则：α=f （ρ，μ，λ，c p ，l ，u ，βgΔT ） 2）量纲分析 ρ[ML -3]，μ[ML -1 T -1]，λ[ML T -3 Q -1]，c p [L 2 T -2 Q -1]，l [L] ，u [LT -1]， βg ΔT [L T -2]， α[MT -3 Q -1]] 3）选基本变量（独立，含M ，L ，T ，Q-热力学温度） ρ，l ，μ, λ 4）无量纲化非基本变量 α：Nu ＝αl/λ u: Re ＝ρlu/μ c p : Pr ＝c p μ/λ βgΔT : Gr ＝βgΔT l 3ρ2/μ2 5）原函数无量纲化 ??? ? ???=223,,μρβλμμρλαtl g c lu F l p 6）实验 Nu ＝ARe a Pr b Gr c 强制对流圆管内表面加热：Nu ＝ARe a 圆管传热基本方程： m t A K t T t T t T t T A K Q ???=-----?=111 22112211 1ln ) ()( 热量衡算方程： )()(12322111t t c q T T c q Q p m p m -=-= 圆管传热牛顿冷却定律： 2 2112211 22211221121 1ln ) ()(ln )()(w w w w w w w w T T T T T T T T A t t t t t t t t A Q -----?=-----?=αα 圆筒壁传导热流量：)] /()ln[)()()/ln(11221122121 2w w w w w w w w t T t T t T t T A A A A Q -----?-?=δλ 空气流量由孔板流量测量：54 .02.26P q v ??= [m 3h -1，kPa] 空气的定性温度：t=(t 1+t 2)/2 [℃] 传热学题目 传热学 1．热流密度 q 与热流量的关系为（以下式子 A 为传热面积，λ为导热系数，h 为对流传热系数）：( ) (A)q＝φA (B)q＝φ/A (C)q=λφ (D)q＝hφ 2．如果在水冷壁的管子里结了一层水垢，其他条件不变，管壁温度与无水垢时相比将：( ) (A)不变 (B)提高 (C)降低 (D)随机改变 3. 当采用加肋片的方法增强传热时，最有效的办法是将肋片加在哪一侧? ( ) (A)传热系数较大的一侧（B）传热系数较小的一侧 (C)流体温度较高的一侧（D）流体温度较低的一侧 4. 导温系数的物理意义是什么? ( ) （A）表明材料导热能力的强弱 （B）反映了材料的储热能力 （C）反映材料传播温度变化的能力 （D）表明导热系数大的材料一定是导温系数大的材料 5. 温度梯度表示温度场内的某一点等温面上什么方向的温度变化率? ( ) (A)切线方向 (B)法线方向 (C)任意方向 (D)温度降低方向 6. 接触热阻的存在使相接触的两个导热壁面之间产生什么影响? ( ) (A)出现温差 (B)出现临界热流 (C)促进传热 (D)没有影响 7. 金属含有较多的杂质，则其导热系数将如何变化? ( ) （A）变大（B）变小 （C）不变（D）可能变大，也可能变小 8. 物体之间发生热传导的动力是什么? ( ) (A)温度场 (B)温差 (C)等温面 (D)微观粒子运动 9. 通过大平壁导热时，大平壁内的温度分布规律是下述哪一种？( ) (A)直线 (B)双曲线 (C)抛物线 (D)对数曲线 10. 已知某一导热平壁的两侧壁面温差是 30℃，材料的导热系数是 22W/(m. K)，通过的热流密度是 300W/m2，则该平壁的壁厚是多少? ( ) (A) 220m (B)22m (C)2.2m (D)0.22m 11. 第二类边界条件是什么? ( ) （A）已知物体边界上的温度分布。 （B）已知物体表面与周围介质之间的传热情况。 （C）已知物体边界上的热流密度。 （D）已知物体边界上流体的温度与流速。 《传热学》 实验指导书工程热物理教研室编 华北电力大学（北京） 二00六年九月 前言 1．实验总体目标 通过本实验，加深学生对传热学基本原理的理解，掌握相关的测量方法，熟练使用相关的测量仪表，培养学生分析问题、解决问题的能力。 ⒉适用专业 热能与动力工程、建筑环境与设备工程、核科学与核工程 ⒊先修课程 高等数学、大学物理 ⒋实验课时分配 ⒌实验环境（对实验室、机房、服务器、打印机、投影机、网络设备等配置及数量要求） 实验用器材及水电等齐全，布局合理，实验台满足2～3人一组，能够满足实验的要求。 在醒目的地方有实验原理的说明，便于教师讲解及学生熟悉实验的基本原理和方法。⒍实验总体要求 每一学期前两周下达实验教学任务，实验教师按照实验任务准备相应的实验设备，实验期间要求学生严格遵守实验室的各项规章制度。学生要提前预习实验内容，完成实验后按照规定格式写好实验报告，交给实验指导教师，实验指导教师根据实验课上的表现和实验报告给出成绩评定结果。 ⒎本实验的重点、难点及教学方法建议 传热学实验的重点是非稳态（准稳态）法测量材料导热性能实验、强迫对流单管外放热系数测定试验、热管换热器实验，这三个实验都是综合性实验，涉及到传热学的导热基本理论，单相对流换热，热边界层理论，相变换热理论，相似原理等方面的内容，实验仪器涉及到热电耦的使用及补偿方法，比托管的使用，电位差计的使用等等，需要学生对传热学的基本内容和基本概念有清楚的了解。 难点是准稳态法测量材料的热性能实验时的测量时间不好把握，单管外放热系数实验中实验数据的整理中存在着一些技巧，另外如何调节热管的加热功率使之得到更好的热管性能曲线也存在一些技巧。 教学方法：学生提前预习，做实验之前老师提问；学生仔细观察指导教师的演示；实验室对学生开放，一次没有做成功，或者想更好地掌握实验技巧的学生，可以跟指导教师预约时间另做。 二维导热物体温度场的数值模拟 一、物理问题 有一个用砖砌成的长方形截面的冷空气通道，其截面尺寸如下图1-1所示，假设在垂直于纸面方向上用冷空气及砖墙的温度变化很小，可以近似地予以忽略。在下列两种情况下试计算： 砖墙横截面上的温度分布；垂直于纸面方向的每米长度上通过砖墙的导热量。 第一种情况：内外壁分别均匀维持在0℃及30℃； 第二种情况：内外壁均为第三类边界条件，且已知： K m K m W h C t K m W h C t ?=?=?=?=?=∞∞/35.0/93.3,10/35.10,302 22211λ砖墙导热系数 二、数学描写 由对称的界面必是绝热面，可取左上方的四分之一墙角为研究对象，该问题为二维、稳态、无内热源的导热问题。 控制方程： 02 222=??+??y t x t 边界条件： 第一种情况： 由对称性知边界1绝热： 0=w q ； 边界2为等温边界，满足第一类边界条件： C t w ?=0； 边界3为等温边界，满足第一类边界条件： C t w ?=30。 第一种情况： 由对称性知边界1绝热： 0=w q ； 边界2为对流边界，满足第三类边界条件： )()( 2f w w w t t h n t q -=??-=λ； 边界3为对流边界，满足第三类边界条件： )()(2f w w w t t h n t q -=??-=λ。 1 -1图2 -1图 三、方程离散 用一系列与坐标轴平行的间隔0.1m 的二维网格线将温度区域划分为若干子区域，如图1-3所示。 采用热平衡法，利用傅里叶导热定律和能量守恒定律，按照以导入元体（m,n ）方向的热流量为正，列写每个节点代表的元体的代数方程， 第一种情况： 边界点： 边界1（绝热边界）： 5~2)2(4 1 1,11,12,1,m =++= +-m t t t t m m m ， 11~8)2(4 1 1,161,16,15,16=++=+-n t t t t n n n n ， 边界2（等温内边界）： 7,16~7;7~1,6,0,=====n m n m t n m 边界3（等温外边界）： 12,16~2;12~1,1,30,=====n m n m t n m 内节点： 11 ~8,15~6;11~2,5~2)(41 1,1,,1,1,====+++= -+-+n m n m t t t t t n m n m n m n m n m 第二种情况 边界点： 边界1（绝热边界）： 5~2)2(4 1 1,11,12,1 ,m =++=+-m t t t t m m m ， 11~8)2(4 1 1,161,16,15,16=++=+-n t t t t n n n n ， 边界2（内对流边界）： 6~1) 2(2221 11,61,6,5,6=++++= ??-+n Bi t Bi t t t t n n n n ， 3 -1图 哈工大-传热学虚拟仿真实验报告 Harbin Institute of Technology 传热学虚拟仿真实验报告 院系：能源科学与工程学院 班级：设计者： 学号： 指导教师：董士奎 设计时间：2016.11.7 传热学虚拟仿真实验报告 1 应用背景 数值热分析在核工业、铁道、石油化工、航空航天、机械制造、能源、汽车交通、国防军工、电子、土木工程、造船、生物医学、轻工、地矿、水利、以及日用家电等各个领域都有广泛的应用。 2 二维导热温度场的数值模拟 2.1 二维稳态导热实例 假设一用砖砌成的长方形截面的冷空气通道，其截面如图2.1所示，假设在垂直于纸面方向上冷空气及砖墙的温度变化很小，可以近似地予以忽略。 图2.1一用砖砌成的长方形截面的冷空气通道截面 2.2二维数值模拟 基于模型的对称性，简化为如图所示的四分之一模 型。 图2.2 二维数值模拟 2.3 建立离散方程 此时对于内部节点，如图2.3： ,1,,1,,,1,,1=? ? - +??-+??-+??--++-x y t t x y t t y x t t y x t t j t j i j t j i j t j i j t j i λ λ λ λ 对于平直边界上的节点，如图2.4： 2 22,,1,,1,,,1=?+Φ??+??-+??-+??-? -+-w j i j t j i j t j i j t j i yq y x x y t t x y t t y x t t λλλ 对于外部和内部角点，如图2.5： 2 43220 2422,,,1,1,,1,,,1,,1,,,1=?+?+Φ??+??-+??-+??-+??-=?+?+Φ??+??-+??-?+-+-?--w n m n m n m n m n m n m n m n m n m w n m n m n m n m n m q y x y x y x t t x y t t x y t t y x t t q y x y x x y t t y x t t λλλλλλ 东南大学能源与环境学院 课程作业报告 作业名称：传热学大作业——利用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; 传热学计算机实习指导书 本指导书是为配合本科生传热学课中计算机应用方面的教学而编写的。 应用计算机解决工程实际问题，是现代工程技术人员所必备的技能。在传热学课程中引入计算机实习的目的，是使学生初步掌握用计算机求解传热问题的技能，从而提高学生应用计算机解决工程实际问题的能力。 大量的传热问题能够用计算机求解。研究如何用计算机求解传热问题的专门知识数值传热学（或称计算传热学）已经发展成了传热学的一个分支学科。传热学课中所涉及的只是数值传热学的初步知识。因此，本次计算机实习也仅仅是作为数值传热学的入门。 本指导书给出了三个练习题及相应的算法。这三个练习题分别涉及了一维稳态导热、二维稳态导热和一维非稳态导热。要求学生在掌握问题的数值计算方法的基础上，独立编写计算机程序并用所编的程序计算出这三个练习题的数值结果。 1 练习题一：一维稳态导热的数值计算 1．1 物理问题 图1示出了一个等截面直肋，处于温度t ∞=80℃的流体中。肋表面与流休之间的对流换热系 数为h=45W/m 2. ℃，肋基处温度t w =300℃，肋端绝热。肋片由铝合金制成，其导热系数为λ=110W/m ℃，肋片厚度为δ=0.01m ，高度为H=0.1m 。试计算肋内的温度分布及肋的总换热量。 1． 2数学描述及其解析解 引入无量纲过余温度∞ ∞ --=t t t t w θ，则以无量纲温度θ描述的肋片导热微分方程及其边界条 件为: 0,1,002 2 2=??====-x H x x m dx d θθθθ (1-2) (1-1) (1-3) 其中A hp m λ= (其中符号含义与教科书杨世铭陶文铨编著《传热学》相同，以下同)。 上述数学模型的解析解为: ()()[]() ()()m H th t t m hp m H ch H x m ch t t t t w w ∞∞∞-= -? -=-φ (1-4) 按式(1-4)计算得到的在肋内各点的温度由表1给出。 1． 3 数值离散 1．3．1 区域离散 在对方程(1-1)~(1-3)进行数值离散之前,应首先进行计算区域的离散。计算区域的离散如图1所示，总节点数取N 。 1．3．2 微分方程的离散 由于方程(1-1)在计算区域内部处处成立，因而对图1所示的各离散点亦成立。对任一节点i 有： 0222=-???? ??i i m dx d θθ 用θ在节点i 的二阶差分代替θ在节点i 的二阶导数,得: 0222 1 1=-?+--+i i i i m x θθθθ 整理上式成迭代形式: ()112 221 -++?+= i i i x m θθθ (i=2,3,…,N-1) (1-5) 1．3．3 边界条件离散 上面得到的离散方程式(1-5)，对所有内部节点都成立，因此每个内部节点都可得出一个类似的方程。事实上，式(1-5)表达的是一个代数方程组。但这个方程组的个数少于未知数i θ (i=1,2, ……,N)的个数。因此，还需要根据边界条件补充进两个方程后代数方程组才封闭。左边界(x=0)为第一类边界条件，温度为已知，因此可以根据式(1-2)直接补充一个方程为： 11==-=∞w w t t θθ 右边界为第三类边界条件，由图1中边界节点N 的向后差分来代替式(1-3)中的导数，得： 01 =?--x N N θθ80210127传热学C
传热实验实验报告
传热学大作业报告 二维稳态导热
传热学上机C程序源答案之一维稳态导热的数值计算
传热学报告
传热膜系数实验报告
西安交通大学传热学大课后复习
传热学上机实验
化工原理实验传热实验报告
传热学题目培训讲学
传热学实验指导书
西安交通大学传热学大作业---二维温度场热电比拟实验
哈工大-传热学虚拟仿真实验报告
传热学MATLAB温度分布大作业完整版
传热学课程上机实习