清华大学航天航空学院传热学试验报告

【实验（一）名称】瞬态热线法测量多孔介质的热导率【实验原理】图1.实验装置示意图白金丝图白金丝图2.物理模型实验装置如图1所示，将一根细长白金丝埋在初始温度均匀的待测材料中，白金丝同时充当加热器和温度传感器，通电加热后，测定白金丝温度随时间的变化，据此推出其周围介质的热导率。该实验的特点是测量时间短，对试样尺寸无特殊要求。物理模型如图2所示，单位长度上加热丝发出的热流为：q=12R/1=IU/1（1）式中，I和U为通过白金丝的电流与加载在白金丝上的电压，R是白金丝的电阻值。白金丝发热量较小，介质可视为无限体，导热微分方程、初始和边界条件：TOC\o"1-5"\h\zdT"2T1dT、 八pc——=九（+ ），r<r<8,t>0（2）pdt dr2rdr 0T=T，t=00一2兀r九 =q，r=r,t>0\o"CurrentDocument"08r， 0解得加热丝表面处待测介质温度：e/、e2q①281—exp(-atu2/r2),T(r,t)—T=\』 " 0dUu⑶0 0兀3九 U3A(U,3)

式中，3是试样与加热丝热容之比的2倍。A(u,3)=[uJ(u)-3J(u)]2+[uY(u)-uY(u)]2(4)式中，J0(u),J1(u)为第一类贝塞尔函数的零阶、一阶函数；Y0(u)、Y1(u)为第二类贝塞尔函数的零阶、一阶函数；u为积分变量。当t足够大：式(3)中指数积分可用级数展开近似，忽略小量，得到:q4兀九[1口史-C](6)q4兀九r2

0式中,欧拉常数C=0.5772,a为介质的热扩散率。令过余温度0:T1t)-T。,由式⑹可得：d0_qd1nt4兀九qd0IUd0 / - / 4兀d1nt4兀ldInt实验中白金丝长径比大于2000,可以忽略端部效应的影响，实验测得白金丝轴向平均温度T([,t)可视为以上各式中的T(r/t)，白金丝平均温度T([,t)与其电阻R的关系如下：R=R01+P(T(r,t)-T0](9)式中，R0是初始温度T0(取当时室温)时白金丝的零点(不通电加热)电阻；通入较大电流后，t时刻白金丝电阻和平均温度分别为Rt和T(r,t)；P为白金丝的电阻温度系数(0.0039K-1)。【实验器材】TOC\o"1-5"\h\z直流电源(AdvantestR6243) 1台多孔介质及样品槽 1套安捷伦数据采集器(主机34970A，模块34901A) 1台电压表 1台白金丝(直径100〃m,99.99%) 若干标准电阻 1个铜康铜热电偶 1支【实验流程】.将白金丝（长度、电阻约为0.2m、3Q）焊在支架上，连接电路、数据采集器和电压表线路后将其放入样品槽底部中间位置；.用待测多孔介质将白金丝埋好，多孔介质总高度约为230mm；.打开数据采集器及其软件，配置采集白金丝上加载电压、环境温度的通道（扫描间隔时间为1s）；.通入恒定小电流I0（白金丝基本不升温），测定白金丝上和标准电阻上（由电压表测定，U0）加载的电压，得到初始温度T0时白金丝的零点电阻R0，导出数据采集卡采集的数据并记录电压表读数；.只打开测定白金丝上加载电压的通道并将扫描间隔时间设为50ms，开始数据采集后通入较大的恒定电流I（白金丝平均温升5〜10K），待白金丝上加载的电压变化较小后停止实验（约采集4分钟），标准电阻上加载的电压U仍由电压表测定，导出数据采集卡采集的数据并记录电压表读数；.结束实验，断开加热丝的供电电路。【实验数据】表1基本参数记录白金丝长度白金丝直径标准电阻I0电压表读数u0I电压表读数U11通过对实验数据表”初始温度.xlsx”，“101（VDC）”一列取平均值得到。柏丝初始电阻R0（U/I）初始环境温度T22通过对实验数据表”初始温度.xlsx”，“103（C）”一列取平均值得到。0.2193m10-4m1Q5.04mA5.04mV5.04mA15.2mV3.02Q 0 19.06℃I1=300mA和口I2=330mA时的原始数据分别见“传热试验数据”文件夹“300mV.xlsx”和“330mV.xlsx”。2【数据整理】数据的整理步骤如下：.通入恒定小电流I0（白金丝基本不升温），测定白金丝上和标准电阻上（由电压表测定，U0）加载的电压，得到初始温度T0时白金丝的零点电阻R0；.通入较大的恒定电流I（白金丝平均温升5〜10K），待白金丝上加载的电压变化较小后停止实验（约采集4分钟），标准电阻上加载的电压U仍由电压表测定。根据R0、I『300mA和铂丝电压Ut，求得铂丝电阻Rt，进而由公式（9）求得铂丝过余温度0。.作“lnt-0曲线”，在直线段部分取间距较大的两点（即取点时舍弃前面一段弧线）,求得斜率k=（。—0）/（lnt—lnt）。2 1 2 14.再根据公式（8）推导出热导率入=IUf/（4klk）0第3步提到的“lnt-0曲线”图是通过matlab作出的。Ii=300mA和I2=330mA时，作图分

别见下:“r：：dhjs：■.■：：■二.•士触一K5.Z0E¥：7.*=0得k=5.714-5000=0.36814,947—3.0073,k=7,493-6,034=0.461别见下:“r：：dhjs：■.■：：■二.•士触一K5.Z0E¥：7.*=0得k=5.714-5000=0.36814,947—3.007IU0.3x0.9253U=1.0245V ,人=t1= =0.300W/(m.K)12 1 4兀lk 4兀x0.2x0,368九二公2 4九二公2 4兀lk20.33x1.02454兀x0.2x0.461=0.292W/(m.K)【误差分析】视介质为无限体，满足一维情形。实际上，介质是有尺寸的，由于与空气接触，对流传热和辐射广泛存在，使得介质表面温度低于理想情况下的数值，进而导致3在“300mA下的时间自然对数一过余温度曲线”图中取点的坐标是通过“DataCursor”实现的。300mA时方法与之相同。白金丝过余温度减小，介质的九偏大。实验推导过程中认为白金丝承受的电压为定值，很显然，电压随时间是有个增大的过程的，仅仅对电压取平均值，会对结果有一定的影响。仪器先测量零点电阻R0。在通较小电流时，白金丝的温度仍然会有较小的上升，导致测得的R0偏大，进而使得计算得到的过余温度偏小，九偏大。测量电压时，仪表是并在电阻两端的，导致一定的电流分流，如果这部分分流较大，对实验结果的影响较大。不过总体上，仪器的测量精确度较高，这方面引起的误差可以忽略不计。【分析与讨论】a.白金丝热容对温度随时间的变化曲线有什么影响？答：电流一定且电阻与温度关系一定的情况下，白金丝热容越大，升高相同的温度吸收的热量越多，对温度变化的滞后效应越明显，导致过余温度变化越慢，计算得到的九越大。b.如何减小端部效应的影响？答：理想情况下的热传导是沿径向的，所以可以在端部截面处放置隔热板，以减小轴向热流；增大长径比。c.根据你的理解，本实验还有什么问题，如何解决？答:本来是想先计算每两个相邻时间点测量得到的〜，再求平均值来得到精确的热导率。但实际上，测量仪器测量的电压脉动很大，导致计算得到的九没有实际价值。这一问题i一方面难以解决，毕竟仪器的测量误差(脉动)很难消除，另一方面也没有必要解决这一问题，本实验所用的方法已经相当合理了。【实验结果】不同热流密度下测量得到的介质热导率略有不同，热流密度越大，测得的热导率越低，但都接近于0.3W/(m.K)，即介质的热导率约为0.3W/(m.K)。【实验（二）名称】测量空气中细线的自然对流换热系数【实验目的】.掌握细线的大空间自然对流换热系数的测量方法；.掌握大空间自然对流换热实验关联式的确定方法；.加深对自然对流换热的理解。【实验原理】图1.实验装置示意图L图2.物理模型实验系统如图1所示，水平焊接在支架上的白金丝既是加热器也是温度传感器。初始时刻白金丝和两端热沉的温度均为环境温度T，通电加热白金丝产生的热量由对流换热、向0基底导热和向环境辐射带走，热沉温度保持不变。保持白金丝与环境间温差在小于20K，辐射换热量可以忽略。假设导热可以忽略，由热量平衡即可求出对流换热系数：h= ［（Lt ^（1）九LD7—Tptf式中，IPtU^为白金丝加热功率，D和L分别为白金丝直径和长度，h为平均对流换热系数，Tf（T0，因为实验中空气温度基本不变）为热边界层以外流体温度。白金丝的平均温度t与其电阻,，间关系表达式：RP=R01+P什-T0）（2）式中，R0是参考温度T0（环境温度）对应的白金丝零点（未通电加热）电阻，P为白金丝的电阻温度系数（0.0039K-1）。【确定准则方程】影响自然对流的换热系数的主要因素有：流体冷热部分的密度差产生的浮升力、流体流动状态、流体的热物性、换热壁面的热状态、换热壁面的几何因素，依据相似原理，其实验关联式可表示为：Nu=C（Gr•Pr＞（3）本实验介质为空气，在实验温度范围内其Pr数随温度变化较小（0.710~0.713），实验关联式可简化为：Nu=C（Gr）"（4）两边同时取自然对数：ln（Nu）=nln（Gr）+InC（5）改变白金丝加热功率，得到一组Nu与Gr数，采用最小二乘法计算出n和C，确定出实验关联式。拟合以ln（Gr）为横坐标，ln（Nu）为纵坐标的直线。努塞尔数Nu：Nu=hl/九（6）格拉晓夫数Gr：Gr=gaATl31V2（7）vAT=T-T（8）式（3）~（8）中，有关物理量为：h为平均对流换热系数，l为水平白金丝的特征长度，入为空气导热系数，av为空气膨胀系数，g为重力加速度，AT为白金丝平均温度T和热边界层以外空气温度T之差，v为空气运动粘度。f【实验器材】直流电源标准电阻空气槽及白金丝支架安捷伦数据采集器（主机34970A，模块34901A）白金丝（直径100〃m,99.99%） 若干铜康铜热电偶 1支【实验步骤】.将白金丝（长度、电阻约为0.2m、3Q）焊接在支架上，连接电路、数据采集器线路后放入样品槽底部中间位置，盖上空气槽盖板；.打开数据采集器，配置采集白金丝上加载电压、标准电阻上加载电压和空气温度的通道；.估算实验中通过白金丝的电流值，估算参数可以取为：h=120~140W/（m2.K），九尸二71.4W/（m-K）；.通入恒定小电流I°（白金丝基本不升温），测定加载在白金丝上、标准电阻上的电压和环境温度，计算出参考温度T°对应的零点电阻R°；.通入较大恒定电流I，使白金丝平均温升约为4、8、12、16、20K，测定加载在白金丝上、标准电阻上的电压和环境温度；.完成实验，断开白金丝供电电路。【实验记录及数据处理】表1数据记录与处理物理量符号或计算式12345白金丝长度L/m0.2072白金丝直径D/m10-4白金丝电阻温度系数P/K-10.0039白金丝零点电阻R=U/1/O14.4mV/5.1mA=2.824初始环境温度T/℃19.8标准电阻Rs/O1空气温度Tf/℃19.819.819.920.020.0标准电阻加载的电压Us/mV98.5130.3149.8175.9238.1通过白金丝的电流IP=U/R/mA98.5130.3149.8175.9238.1白金丝上加载的电压UPt/mV282.4377.4436.4517.8721.5白金丝电阻RP=UP/IP/O2.8672.8962.9132.9443.030白金丝加热功率P=IPUP/mW27.8249.1865.3791.08171.79

白金丝平均温升AT=(R-R)/(「R)/℃Pt 0 03.906.548.0810.9018.70空气膨胀系数(定压)a=-(ap/at)/p/K-143.41X10-33.41X10-33.41X10-33.41X10-33.41X10-3定性温度Tm=(Tpt+T)/2/℃21.6023.0723.9425.4529.35空气热导率5根据Tm查取/W/(m•k)2.6028X10-22.6146X10-22.6215X10-22.6336X10-22.6648X10-2空气运动粘度根据理想气体状态方程P=P根据理想气体状态方程P=PRT得a=-(ap/at)/p=1/T。5根据教材P559【附录5】利用插值法求得。 '6根据教材P559【附录5】利用插值法求得。根据Tm查取/Pa•s15.21X10-615.35X10-615.43X10-615.57X10-615.94X10-6换热系数h=P/W/(m2•k)九LDAT113.5119.7128.8133.0146.2努塞尔数Nu=hl/Ml=d)0.4360.4580.4910.5050.549格拉晓夫数Gr=gaATI3/v2v5.64X10-49.29X10-41.135X10-31.504X10-32.462X10-31)整理变加热功率实验结果，确定实验关联式对Nu与Gr采用最小二乘法计算出n和C,确定出实验关联式。拟合以ln(Gr)为横坐标，ln(Nu)为纵坐标的直线，将原曲线和拟合线对比如下见下：拟合得到

In(Nu)=0.1603ln(Gr)+0.3612n=0.1603,C=1.435即实验关联式为：Nu=1.435(Gr”16032)考虑白金丝向基底导热时对流换热系数的测定原理并求出对流换热系数由于白金丝长径比很大，可以认为径向温度均匀。列出微分方程为：、d2T九 十6=03X2其中小为等效内热源，分为两部分：电流产生的热量和空气的对流换热①IUdx;L-hlD(T-TJdx

0.25兀D2dx取过余温度9=T—TfTOC\o"1-5"\h\z326 4h6 4IU - Z-+ 3x2 D九 lD2九L所以测定温度沿白金丝轴向分布，作出6-x曲线，根据下式即可求出h।D九326 IUh= +463x2 lD6L3)选取其中一个加热功率，画出考虑白金丝导热时轴向温度分布图。微分方程323263x24h6 4PD九+lD2九L根据图2可知，边界条件为x=0,6=0；x=L,6=0.4h, 4P方程较为复杂，其通解为一个非齐次特解加上一个齐次通解。令a二灰，b二访迁。先求解齐次通解:326.a6=03x20=ce-ax+ceax再求非齐次特解:所以通解为0=ce-■ax+ceax代入边界条件得:bbc=—(e-al-1)/(eal-e-al)；c=---c

2a 1a2设定相关数据为加热功率为65.37mW下的数据。由于a=1.9653e+008,b=1.5875e+009,将cl、c2代入原方程忽略小量化简得：0二—（1—e-ax—e'ax-'"al）aMatlab代码为clearall;clc;h=128.8;d=10A-4;lambda=2.6215*10A-2;p=0.06537;l=0.2;a=4*h/(d*lambda)b=4*p/(pi*dA2*lambda*l)fori=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（'过余温度'）求得的图像为aD.02D.Ofl016Q.L60.1DJ2D.I40.16 0.18 0.2弦画］也遇所以轴向过余温度基本为定值8.08,只在两端有较大偏差。这表明将白金丝轴向温度视为均匀是可行的。【分析与讨论】1）观察铜棒顶部，分析为什么设置斜面？答：顶部设置为斜面，有利于迅速地将上升气流排开，防止气体滞留，保证能够实时测量上升气流温度。2）考虑辐射换热，对实验结果会有什么样的影响？答：测量的实际上是辐射与对流换热的总和，所以考虑到辐射换热，实际的对流换热系数会减小。3）将白金丝换成纳米线，可以研究纳米材料的自然对流换热吗？会遇到什么问题？答：可以，但计算方法需要改变。纳米线是微尺度的，必须考虑分子运动的影响，在传热学方面有专门的Kn数来衡量微观尺度的影响。在这种尺度下，傅里叶导热定律等并不适用，需要用其他计算方法。4）测量流体热导率时，白金丝应如何放置，为什么？答：流体具有流动性，所以在测量热导率时，应尽量减小对流换热的影响。为此，应该将流体层设置的足够薄。5）根据你的理解，本实验还哪些方面应改进？答：测量次数太少，导致得到的关系曲线相当扭曲。所以建议多测几次。【实验结果】本实验中，空气中细线的自然对流换热系数 h与格拉晓夫数Gr的实验关联式为h=1.435-（Gr）0」603【实验（三）名称】测量水中水平细线自然对流和过冷沸腾换热系数【实验目的】.掌握水的自然对流换热和过冷沸腾换热系数的测定方法；.深化对对流换热系数影响因素的理解；.观察过冷沸腾现象，建立对过冷沸腾的认识。【实验原理】过冷沸腾是指液体主体温度低于饱和温度，而加热壁面温度高于饱和温度时所发生的沸腾换热。随主流体温度的增加，过冷沸腾可分为单相对流区、过冷沸腾第一区（高过冷沸腾区）、过冷沸腾第二区（低过冷沸腾区）。在单相对流区，加热壁面温度TW未达到过冷沸腾起始点对应的壁面温度，加热壁面上不生成汽泡。在过冷沸腾起始点，由于产生汽泡吸收了潜热，扰动边界层，改善了传热，壁面温度会有一些下降。随着液体温度继续升高，出现沿壁面移动并在运动过程中长大或碰并长大而脱离壁面的气泡，该温度称为充分发展过冷沸腾起始温度。当主流体温度达到沸点时，开始饱和沸腾。实验系统如图1所示，由辅助加热器将去离子水加热升温至45〜60℃，浸没于水中的白金丝同时充当加热器和温度传感器。向白金丝通入较小电流I测定水的自然对流换热系数，1通入较大电流12测定水的过冷沸腾换热系数并观察过冷沸腾时白金丝表面气泡运动现象。由白金丝上加载的电压和通过电流得到其热负荷。和电阻R。图1.实验装置示意图由通电加热后白金丝电阻R与参考温度［（当时水温）下白金丝的电阻R0确定白金丝表面温度T和过热度AT：wR=R0（1+PAT）（1）AT=T-T0（2）忽略导热和辐射，由热平衡确定水的自然对流换热系数或过冷沸腾换热系数：h=IU/（AAT）（3）式（1）〜（3）中，P为白金丝的电阻温度系数（0.0039K-1），h为自然对流或过冷沸腾换热系数，IU为加热功率，A为白金丝的换热面积。【实验器材】TOC\o"1-5"\h\z直流电源 1台标准电阻（1Q和0.001Q） 2个水槽及白金丝支架 1套安捷伦数据采集器 1台辅助加热器 1个白金丝（直径100^m，99.99%） 若干铜康铜热电偶 1支【实验流程】.将白金丝（长度约为0.2m、环境温度下电阻约为3Q）焊接在支架上，连接好电路、数据采集器线路后放入水槽中，向水槽内注去离子水至水槽高度的2/3处；.估算实验中通过白金丝的电流值，估算参数可以取为：自然对流换热（白金丝温度低于100℃）系数h=104W/（m2K），过冷沸腾换热（白金丝温度为105〜110℃）系数取为h=2.2x104W/（m2・K）；.向白金丝通入小电流I（白金丝基本不升温，用1Q标准电阻），测定其在当时水温T0下电阻R；0.通入较大电流（（发生自然对流，用0.001Q标准电阻），测定加载在白金丝上和标准电阻上的电压及水温；.通入大电流12（发生过冷沸腾，用0.001Q标准电阻），记录观察到的沸腾现象，测定加载在白金丝上和标准电阻上的电压及水温；.完成实验，断开白金丝的供电电路。【实验数据】表1基本参数记录白金丝长度/m白金丝直径/mI0/mAR0/mAI