清华大学航天航空学院“传热学”实验报告

1、【实验（一）名称】 瞬态热线法测量多孔介质的热导率【实验原理】235647891. 直流电源6. 白金丝支架2. 标准电阻7. 支架底座3. 多孔介质槽8. 安捷伦数据采集器14. 接线柱9. 电压表5.白金丝图 1.实验装置示意图图 2.物理模型实验装置如图1 所示，将一根细长白金丝埋在初始温度均匀的待测材料中，白金丝同时充当加热器和温度传感器，通电加热后， 测定白金丝温度随时间的变化，据此推出其周围介质的热导率。该实验的特点是测量时间短，对试样尺寸无特殊要求。物理模型如图2 所示，单位长度上加热丝发出的热流为：q I 2 R / lIU / l （ 1）式中， I 和 U 为通过白金丝的电

2、流与加载在白金丝上的电压，R 是白金丝的电阻值。白金丝发热量较小，介质可视为无限体，导热微分方程、初始和边界条件：cp T(2T 1 T ) ， r0r,t 0 （ 2）tr 2rrTT0 ， t02 r0Tr0 ,t0q ， rr解得加热丝表面处待测介质温度：T (r0 ,t ) T02q21exp(tu2 / r02 )3u3(u,du （ 3）0)式中， 是试样与加热丝热容之比的2 倍。(u, ) uJ0 (u)J1 (u) 2uY0 (u)uY1(u) 2（ 4）式中， J0(u) ， J1(u)为第一类贝塞尔函数的零阶、一阶函数；Y0(u)、 Y1(u)为第二类贝塞尔函数的零阶、一阶

3、函数； u 为积分变量。当 t 足够大：r021（ 5）4t式（ 3）中指数积分可用级数展开近似，忽略小量，得到：T (r0 ,t )T0qln 4tC （6）4r02式中， 欧拉常数 C 0.5772，为介质的热扩散率。 令过余温度T (r0 ,t )T0，由式（ 6）可得：dq（ 7）d ln t4q /dIU/d（ 8）4d ln t4ld ln t实验中白金丝长径比大于2000，可以忽略端部效应的影响，实验测得白金丝轴向平均温度 T (r ,t ) 可视为以上各式中的T (r0 , t) ，白金丝平均温度T (r, t ) 与其电阻 Rt 的关系如下：00RtR01T (r0 ,t)-

4、T0（ 9）式中， R0 是初始温度 T0 （取当时室温）时白金丝的零点（不通电加热）电阻；通入较大电流后， t 时刻白金丝电阻和平均温度分别为Rt 和 T (r0 , t) ；为白金丝的电阻温度系数（0.0039K-1）。【实验器材】直流电源（ Advantest R6243）1 台多孔介质及样品槽1 套安捷伦数据采集器（主机 34970A，模块 34901A）1 台电压表1 台白金丝（直径 100m，99.99%）若干标准电阻1 个铜康铜热电偶1 支【实验流程】1. 将白金丝（长度、电阻约为 0.2m、3 ）焊在支架上，连接电路、数据采集器和电压表线路后将其放入样品槽底部中间位置；2.用待

5、测多孔介质将白金丝埋好，多孔介质总高度约为230mm ；3. 打开数据采集器及其软件， 配置采集白金丝上加载电压、 环境温度的通道 （扫描间隔时间为 1s）；4. 通入恒定小电流 I 0（白金丝基本不升温） ，测定白金丝上和标准电阻上（由电压表测定， U0）加载的电压，得到初始温度 T0 时白金丝的零点电阻 R0，导出数据采集卡采集的数据并记录电压表读数；5.只打开测定白金丝上加载电压的通道并将扫描间隔时间设为50ms，开始数据采集后通入较大的恒定电流I（白金丝平均温升510K ），待白金丝上加载的电压变化较小后停止实验（约采集4 分钟），标准电阻上加载的电压U 仍由电压表测定，导出数据采集卡

6、采集的数据并记录电压表读数；6. 结束实验，断开加热丝的供电电路。【实验数据】表 1 基本参数记录白金丝标准I0电压表读电压表铂丝初始电初始环境温度白金丝长度I读数 U12直径电阻数 U0阻 R0(U/ I)T00.2193m10-4m15.04mA5.04mV5.04mA15.2mV3.0219.06 CI1=300mA 和 I2=330mA 时的原始数据分别见“传热试验数据”文件夹“300mV.xlsx” 和“330mV.xlsx”。【数据整理】数据的整理步骤如下：1. 通入恒定小电流 I 0（白金丝基本不升温） ，测定白金丝上和标准电阻上（由电压表测定， U 0）加载的电压，得到初始温度

7、T0 时白金丝的零点电阻R0 ；2. 通入较大的恒定电流 I（白金丝平均温升 510K ），待白金丝上加载的电压变化较小后停止实验（约采集4 分钟），标准电阻上加载的电压U 仍由电压表测定。根据R0、I1=300mA 和铂丝电压 U t，求得铂丝电阻 Rt，进而由公式（ 9）求得铂丝过余温度 。3. 作“ lnt- 曲线”，在直线段部分取间距较大的两点（即取点时舍弃前面一段弧线），求得斜率 k ( 21）（/lnt 2lnt 1）。4. 再根据公式（ 8）推导出热导率IU t / (4lk ) 。第 3 步提到的“ lnt- 曲线”图是通过 matlab作出的。I1=300mA2时，作图分和

8、I =330mA1通过对实验数据表“初始温度.xlsx”，“ 101(VDC)”一列取平均值得到。2通过对实验数据表“初始温度.xlsx”，“ 103(C)”一列取平均值得到。别见下：得 k15.7145.0000.368 3 ， k2 = 7.4936.0340.461 ， U t 1 0.9253V，4.9473.0075.2022.035U t 21.0245V，1I1U t10.30.92530.300W / (m.K )，4 lk 14 0.2 0.3682I 2U t 20.33 1.02450.292W / ( m.K )4 lk 2 40.20.461【误差分析】1. 视介质为

9、无限体，满足一维情形。实际上，介质是有尺寸的，由于与空气接触，对流传热和辐射广泛存在，使得介质表面温度低于理想情况下的数值，进而导致3在“ 300mA 下的时间自然对数过余温度曲线”图中取点的坐标是通过“Data Cursor”实现的。 300mA 时方法与之相同。白金丝过余温度减小，介质的偏大。2.实验推导过程中认为白金丝承受的电压为定值，很显然，电压随时间是有个增大的过程的，仅仅对电压取平均值，会对结果有一定的影响。3. 仪器先测量零点电阻 R0。在通较小电流时，白金丝的温度仍然会有较小的上升，导致测得的R0 偏大，进而使得计算得到的过余温度偏小，偏大。4. 测量电压时，仪表是并在电阻两端

10、的，导致一定的电流分流，如果这部分分流较大，对实验结果的影响较大。不过总体上，仪器的测量精确度较高，这方面引起的误差可以忽略不计。【分析与讨论】a. 白金丝热容对温度随时间的变化曲线有什么影响？答：电流一定且电阻与温度关系一定的情况下，白金丝热容越大，升高相同的温度吸收的热量越多，对温度变化的滞后效应越明显，导致过余温度变化越慢，计算得到的越大。b. 如何减小端部效应的影响？答：理想情况下的热传导是沿径向的，所以可以在端部截面处放置隔热板，以减小轴向热流；增大长径比。c. 根据你的理解，本实验还有什么问题，如何解决？答：本来是想先计算每两个相邻时间点测量得到的i ，再求平均值来得到精确的热导率

11、。但实际上，测量仪器测量的电压脉动很大，导致计算得到的i 没有实际价值。这一问题一方面难以解决，毕竟仪器的测量误差（脉动）很难消除，另一方面也没有必要解决这一问题，本实验所用的方法已经相当合理了。【实验结果】不同热流密度下测量得到的介质热导率略有不同，热流密度越大，测得的热导率越低，但都接近于0.3W / (m.K ) ，即介质的热导率约为0.3W / ( m.K ) 。【实验（二）名称】 测量空气中细线的自然对流换热系数【实验目的】1. 掌握细线的大空间自然对流换热系数的测量方法；2. 掌握大空间自然对流换热实验关联式的确定方法；3. 加深对自然对流换热的理解。【实验原理】243576891

12、.直流电源6.接线柱2.标准电阻7.白金丝支架3.空气槽8.支架底座14.热电偶9.安捷伦数据采集器5.白金丝图 1.实验装置示意图LT0白金丝T0x图 2.物理模型实验系统如图 1 所示，水平焊接在支架上的白金丝既是加热器也是温度传感器。初始时刻白金丝和两端热沉的温度均为环境温度T0 ，通电加热白金丝产生的热量由对流换热、向基底导热和向环境辐射带走，热沉温度保持不变。保持白金丝与环境间温差在小于20K，辐射换热量可以忽略。假设导热可以忽略，由热量平衡即可求出对流换热系数：hI PtU Pt（ 1）LD TPtTf式中， I PtU Pt 为白金丝加热功率，D 和 L分别为白金丝直径和长度，h

13、 为平均对流换热系数， Tf （ T0 ，因为实验中空气温度基本不变）为热边界层以外流体温度。白金丝的平均温度TPt 与其电阻 RPt 间关系表达式：RPt R0 1TPt -T0 （ 2）式中， R0 是参考温度 T0 （环境温度）对应的白金丝零点（未通电加热）电阻，为白金丝的电阻温度系数（0.0039K-1）。【确定准则方程】影响自然对流的换热系数的主要因素有：流体冷热部分的密度差产生的浮升力、流体流动状态、流体的热物性、换热壁面的热状态、换热壁面的几何因素，依据相似原理，其实验关联式可表示为：NuC Gr Prn（ 3）本实验介质为空气，在实验温度范围内其Pr 数随温度变化较小（0.71

14、00.713 ），实验关联式可简化为：NuC Grn（ 4）两边同时取自然对数：ln Nun lnGrln C （ 5）改变白金丝加热功率，得到一组Nu 与 Gr 数，采用最小二乘法计算出n 和 C，确定出实验关联式。拟合以 ln Gr为横坐标， lnNu为纵坐标的直线。努塞尔数 Nu ：Nuhl /（ 6）格拉晓夫数 Gr ：GrgvTl 3 /2（7）TTPt Tf （ 8）式（ 3） （ 8）中，有关物理量为：h 为平均对流换热系数，l 为水平白金丝的特征长度， 为空气导热系数，v 为空气膨胀系数，g 为重力加速度，T 为白金丝平均温度 TPt和热边界层以外空气温度Tf 之差，为空气运动

15、粘度。【实验器材】直流电源1 台标准电阻1 个空气槽及白金丝支架1 套安捷伦数据采集器（主机34970A，模块 34901A）1 台白金丝（直径100m， 99.99%）若干铜康铜热电偶1 支【实验步骤】5. 将白金丝（长度、电阻约为 0.2m、3 ）焊接在支架上，连接电路、数据采集器线路后放入样品槽底部中间位置，盖上空气槽盖板；6. 打开数据采集器， 配置采集白金丝上加载电压、 标准电阻上加载电压和空气温度的通道；7. 估算实验中通过白金丝的电流值，估算参数可以取为：h 120140W / m2 K ，Pt 71.4W /(m K ) ；8. 通入恒定小电流 I 0（白金丝基本不升温） ，测

16、定加载在白金丝上、标准电阻上的电压和环境温度，计算出参考温度T0 对应的零点电阻R0 ；9. 通入较大恒定电流 I ，使白金丝平均温升约为 4、 8、12、16、20K ，测定加载在白金丝上、标准电阻上的电压和环境温度；10. 完成实验，断开白金丝供电电路。【实验记录及数据处理】表 1 数据记录与处理物理量白金丝长度白金丝直径白金丝电阻温度系数白金丝零点电阻初始环境温度标准电阻空气温度标准电阻加载的电压通过白金丝的电流白金丝上加载的电压白金丝电阻白金丝加热功率符号或计算式12345L/ m0.2072D/ m10-4/ K-10.0039R0U0/ I0/14.4mV/5.1mA=2.824T

17、0/C19.8Rs/ 1Tf/ C19.819.819.920.020.0Us/mV98.5130.3149.8175.9238.1I PtU s / Rs /mA98.5130.3149.8175.9238.1UPt/mV282.4377.4436.4517.8721.5RPtU Pt / I Pt / 2.8672.8962.9132.9443.030PI PtU Pt /mW27.8249.1865.3791.08171.79白金丝平均温升TRPtR0（/ C3.906.548.0810.9018.70R0）1 43.41×3.41×3.41×3.41

18、15;3.41×空气膨胀系数（定压）/TP / K10-3v10-310-3-3-31010定性温度空气热导率空气运动粘度换热系数努塞尔数格拉晓夫数Tm(TPtT f)/ 2/ C21.6023.0723.9425.4529.355根据 Tm 查取 /W/(m ·k)2.60282.61462.62152.63362.6648×10-2×10-2×10-2× 10-2-2× 1015.21×15.3515.43×15.57×15.94 ×6根据 Tm查取 /Pa· s10-6

19、10-6×10-610-610-6hP/W/(m 2· k)113.5119.7128.8133.0146.2LDTNuhl /(ld)0.4360.4580.4910.5050.5495.64×9.29×1.135×1.504×2.462 ×GrgTl3/210-4v10-410-310-310-31) 整理变加热功率实验结果，确定实验关联式对 Nu 与 Gr 采用最小二乘法计算出 n 和 C，确定出实验关联式。 拟合以 ln Gr 为横坐标， ln Nu 为纵坐标的直线，将原曲线和拟合线对比如下见下：拟合得到4根据理想气

20、体状态方程pRT 得 v/T P /1/ Tf 。5根据教材 P559【附录 5】利用插值法求得。6根据教材 P559【附录 5】利用插值法求得。ln Nu0.1603ln Gr +0.3612n0.1603, C1.435即实验关联式为：0.1603Nu1.435 Gr2) 考虑白金丝向基底导热时对流换热系数的测定原理并求出对流换热系数由于白金丝长径比很大，可以认为径向温度均匀。列出微分方程为：2Tx 20其中为等效内热源，分为两部分：电流产生的热量和空气的对流换热IUdx LhD(TTf )dx0.25D2dx取过余温度TTf24h4IU0x 2DD2L所以测定温度沿白金丝轴向分布，作出x

21、 曲线，根据下式即可求出 hD2IUh4x 2D L3) 选取其中一个加热功率，画出考虑白金丝导热时轴向温度分布图。微分方程244Ph0x 2DD2 L根据图 2 可知，边界条件为x0,0;xL,0.方程较为复杂，其通解为一个非齐次特解加上一个齐次通解。令a4h4P。， bD 2DL先求解齐次通解：2x 2a0c1eaxc2 e ax再求非齐次特解：ba所以通解为c1eaxc2e axba代入边界条件得：c2b (e al1) / (e ale al ); c1bc2aa设定相关数据为加热功率为65.37mW 下的数据。由于a = 1.9653e+008 ， b =1.5875e+009 ，将

22、 c1、 c2 代入原方程忽略小量化简得：b (1eaxe axal )aMatlab 代码为clear all ;clc;h=128.8;d=10-4;lambda=2.6215*10-2;p=0.06537;l=0.2;a=4*h/(d*lambda)b=4*p/(pi*d2*lambda*l)for i=1:101x(i)=0.002*(i-1);theta(i)=b/a*(1-exp(-sqrt(a)*x(i)-exp(sqrt(a)*(x(i)-l);endplot(x,theta,'r');xlabel( '弦向位置 ');ylabel( '

23、过余温度 ')求得的图像为所以轴向过余温度基本为定值 8.08 ，只在两端有较大偏差。这表明将白金丝轴向温度视为均匀是可行的。【分析与讨论】1） 观察铜棒顶部，分析为什么设置斜面？答：顶部设置为斜面，有利于迅速地将上升气流排开，防止气体滞留，保证能够实时测量上升气流温度。2） 考虑辐射换热，对实验结果会有什么样的影响？答：测量的实际上是辐射与对流换热的总和，所以考虑到辐射换热，实际的对流换热系数会减小。3） 将白金丝换成纳米线，可以研究纳米材料的自然对流换热吗？会遇到什么问题？答：可以，但计算方法需要改变。纳米线是微尺度的，必须考虑分子运动的影响，在传热学方面有专门的Kn 数来衡量微观

24、尺度的影响。在这种尺度下，傅里叶导热定律等并不适用，需要用其他计算方法。4） 测量流体热导率时，白金丝应如何放置，为什么？答：流体具有流动性，所以在测量热导率时，应尽量减小对流换热的影响。为此，应该将流体层设置的足够薄。5） 根据你的理解，本实验还哪些方面应改进？答：测量次数太少，导致得到的关系曲线相当扭曲。所以建议多测几次。【实验结果】本实验中，空气中细线的自然对流换热系数h 与格拉晓夫数Gr 的实验关联式为h 1.435Gr0.1603。d【实验（三）名称】 测量水中水平细线自然对流和过冷沸腾换热系数【实验目的】4. 掌握水的自然对流换热和过冷沸腾换热系数的测定方法；5. 深化对对流换热系

25、数影响因素的理解；6. 观察过冷沸腾现象，建立对过冷沸腾的认识。【实验原理】过冷沸腾是指液体主体温度低于饱和温度，而加热壁面温度高于饱和温度时所发生的沸腾换热。随主流体温度的增加，过冷沸腾可分为单相对流区、过冷沸腾第一区（高过冷沸腾区）、过冷沸腾第二区（低过冷沸腾区）。在单相对流区，加热壁面温度TW 未达到过冷沸腾起始点对应的壁面温度，加热壁面上不生成汽泡。在过冷沸腾起始点，由于产生汽泡吸收了潜热，扰动边界层，改善了传热，壁面温度会有一些下降。随着液体温度继续升高，出现沿壁面移动并在运动过程中长大或碰并长大而脱离壁面的气泡，该温度称为充分发展过冷沸腾起始温度。当主流体温度达到沸点时，开始饱和沸

26、腾。实验系统如图1 所示，由辅助加热器将去离子水加热升温至4560， 浸没于水中的白金丝同时充当加热器和温度传感器。向白金丝通入较小电流I 1 测定水的自然对流换热系数，通入较大电流I 2 测定水的过冷沸腾换热系数并观察过冷沸腾时白金丝表面气泡运动现象。由白金丝上加载的电压和通过电流得到其热负荷Q和电阻 R 。24985 637101. 直流电源6. 辅助加热器2. 标准电阻7.支架底座13. 接线柱8.热电偶4. 白金丝9.水槽5. 白金丝支架10. 安捷伦数据采集器图 1.实验装置示意图由通电加热后白金丝电阻R 与参考温度 T0 （当时水温） 下白金丝的电阻R0确定白金丝表面温度 Tw 和

27、过热度T ：R R01T （1）TTwT0 （ 2）忽略导热和辐射，由热平衡确定水的自然对流换热系数或过冷沸腾换热系数：h IU / A T （ 3）式（ 1）（3）中，为白金丝的电阻温度系数（0.0039K-1）， h 为自然对流或过冷沸腾换热系数， IU 为加热功率， A 为白金丝的换热面积。【实验器材】直流电源1 台标准电阻（ 1和 0.001 ）2 个水槽及白金丝支架1 套安捷伦数据采集器1 台辅助加热器1 个白金丝（直径 100 m， 99.99%）若干铜康铜热电偶1 支【实验流程】11. 将白金丝（长度约为 0.2m、环境温度下电阻约为 3）焊接在支架上，连接好电路、数据采集器线路

28、后放入水槽中，向水槽内注去离子水至水槽高度的2/3 处；12. 估算实验中通过白金丝的电流值， 估算参数可以取为： 自然对流换热 （白金丝温度低于 100）系数h=10 4W/(m 2·K) ，过冷沸腾换热（白金丝温度为105110）系数取4 2为 h=2.2 ×10 W/(m ·K) ；13. 向白金丝通入小电流 I 0（白金丝基本不升温， 用 1 标准电阻），测定其在当时水温 T0下电阻 R0 ；14. 通入较大电流 I1 （发生自然对流，用 0.001 标准电阻），测定加载在白金丝上和标准电阻上的电压及水温；15. 通入大电流 I 2 （发生过冷沸腾，用 0.001 标准电阻），记录观察到的沸腾现象，测定加载在白金丝上和标准电阻上的电压及水温；16. 完成实验，断开白金丝的供电电路。【实验数据】表 1基本参数记录白金丝长度白金丝直径I 0/mAR0 /mAI1/AI 2/A初始水温实验后水温/m/mT0/ CT / C0.