高祎,寇琼月,张灏,唐亚曼
(1.航空工业北京长城计量测试技术研究所,北京 100095;2.中国人民解放军92228部队,北京 100072)
核电池因其使用寿命长、寿命周期内维护频次低、功率稳定等优点,在空间探测、核工业、反潜监听等领域发挥着重要作用[1-2]。由于核电池长期工作在难以进行人工维护的环境中,因此需要对其进行功率测量,以确保准确核算探测器寿命、极限功率载荷、核燃料纯度等关键参数。核电池功率通过量热计测量,量热计的标定与溯源需要标准热源提供参考标准值。
美国核安全管理局(NNSA)洛斯阿拉莫斯国家实验室应用惠斯通桥原理测量核电池热功率,其标准热源应用经认证的标准Pu-238(准确度优于0.02%)及电加热标准热源[3-4]。法国塞塔拉姆公司开发的CALVET LV量热仪,使用焦耳效应标定,量热准确度为±0.2%[5]。柳加成等[6]应用基于热电堆原理的量热计测量热功率,其标准热源选择同位素热源和电加热模拟热源两种,电加热模拟12.81 W热源时,短期稳定性相对标准偏差为0.44%。张萍等[7]设计了射频中功率量热计校准系统,其中标准热源负载选择标称阻值为50Ω的负载电阻。
相较同位素热源,电加热模拟热源成本较低且无放射性,因此更具有普适性。应用电加热模拟热源时,基于量热法测量核电池热功率的溯源链为:核电池热功率溯源到热功率测量装置,热功率测量装置溯源到标准热源,标准热源溯源到电学基准。
目前我国用于对核电池热功率测量装置进行计量的电加热标准热源存在准确度较低、稳定性较差的问题。针对此问题,本文研制了一套热功率标准装置,该装置基于电流的焦耳效应原理将电能转化为热能,在电流小于等于3 A的范围内,能够产生0.7~147 W的功率。本文对该装置的结构组成、工作原理进行详细介绍,并开展性能试验及不确定度分析以验证其可靠性与准确性,为相关热功率标准装置的国产化设计研发提供借鉴。
研制的标准热源装置由高精度恒流源、标准电阻、发热体、八位半数字多用表等部分组成,如图1所示。高精度恒流源为电路提供稳定的电流,发热体与标准电阻串联在电路中。高精度数字多用表用于测量发热体和标准电阻两端的电压,利用低热电势扫描开关实现测量位置的切换。
图1 标准热源装置图Fig.1 Schematic diagram of standard heat source
热功率标准热源装置基于电流的焦耳效应将电能转化为热能,通过电功率计算产生的热功率。由于发热体和标准电阻为串联结构,因此通过二者的电流相同,标准电阻两端电压U1与电阻R0的比值即为通路中的电流。功率P为发热体两端电压U2与电流的乘积,故
1.2.1 发热体骨架
发热体骨架的形状应与测量对象适配,本文研制的热功率标准热源主要用于对核电池热功率测量装置进行计量,目前大部分核电池热功率测量装置均为圆柱体,因此将发热体骨架计为圆筒状,这种设计方式还具有达到热平衡的时间较短的优点。
根据不同的功率,将发热体整体尺寸分为3种,功率与发热体尺寸的对应关系如表1所示。
表1 发热体尺寸Tab.1 Size of heating element
1.2.2 加热丝材料选择
根据发热体的工作特性要求选择加热丝材料。标准热源装置作为标准器具对核电池热功率进行校准,其长期稳定性至关重要,加热丝材料需要具有抗氧化、抗腐蚀、不易受环境影响、温度分布均匀的特点。对比常用加热丝材料的上述特性(如表2所示),选择铂丝作为加热丝材料。
表2 加热丝材料选型Tab.2 Material selection of heating wire
1.2.3 加热丝缠绕方式
加热丝缠绕方式可分为内绕式和外绕式,内绕式无法提供支撑应力,耐震性差;外绕式能够提供支撑应力,耐震性相对较好。本文采用双螺旋外绕式缠绕方式缠绕加热丝,以减小电流的磁效应。具体方法为:将加热丝对折,对折端固定在骨架底部,然后在骨架表面双线螺旋向上缠绕加热丝。发热体骨架外表面开有双螺旋状沟槽用以固定加热丝。为保证发热体表面均匀性,双螺旋沟槽中的两个螺旋线螺距一致,螺旋线方向一致。
1.2.4 加热丝直径和长度选择
加热丝的直径D、长度L和功率P的关系需要通过公式(2)和(3)计算获得
式中:R为电阻值,Ω;I为通路中电流,A。
式中:ρ为加热丝的电阻率,Ω·m;L为加热丝的有效长度,m;D为加热丝的直径,m。
由于高功率情况下电阻阻值较大,若采用直径过大的加热丝,会导致所需长度过长,缠绕时每匝间距过近,实际使用过程中容易出现短路。若采用直径过小的加热丝,会导致加热丝强度较低,缠绕过程中易折断。因此,加热丝直径的选择应综合考虑每匝缠绕间隔和加热丝强度。开展了不同功率下加热丝直径和长度选择试验(数据如表3所示),根据试验情况,最终设计:功率为1~10 W时,选择直径为0.2 mm的加热丝,长度约为314 mm;功率为10~100 W和100~150 W的情况下,选择直径为0.5 mm的加热丝,长度分别约为10239 mm和19130 mm。
表3 加热丝直径和长度Tab.3 Diameter and length of heating wire
发热体装置结构如图2所示。发热体的外壳、骨架、顶盖、固定柱均选择石英玻璃材料制造,具有绝缘性好、耐热性强、热膨胀系数小、化学热稳定性好等优点。发热体外壳和顶盖的厚度均为2 mm,主要起绝缘作用。固定柱伸出骨架内壁10 mm,便于固定加热丝。
图2 发热体装置图Fig.2 Schematic diagram of heating element
开展发热体电阻值及温度随时间变化试验、发热体产生功率随电流变化试验、发热体产生功率的短期稳定性试验,以验证本文研制的标准热源装置的实际性能。
标准热源装置如图3所示,高精度恒流源提供恒定电流,八位半数字多用表与寄生电势小于0.2μV的扫描切换开关相连。其中,标准电阻连接扫描开关的标准通道,发热体连接扫描开关的01通道,通过转换扫描开关的通道来测量发热体和标准电阻两端的电压值。
图3 标准热源装置图Fig.3 Image diagram of standard heat source
以下步骤为试验验证的通用步骤:
步骤一,在发热体表面贴一支铂电阻温度计,在标准电阻四周贴4支铂电阻温度计,用以监测发热体表面温度和标准电阻温度;
步骤二,打开高精度恒流源、扫描开关及数字多用表的电源;
步骤三,设置数字多用表DCV(直流电压档),扫描开关设置为标准通道,设置高精度恒流源的电流,设置完成后,发热体电阻值和温度随时间变化试验间隔为2 min,功率随电流的变化以及发热体产生功率的稳定性试验间隔为1 h,读取数字多用表的示数并记录;读数完成后,开关通道设置为01通道,读取数字电压表的示数并记录;
步骤四,更改恒流源电流,重复步骤三。
三个试验的区别在于发热体电阻值和温度随时间变化的试验中记录数据的时间间隔与其他试验不同。发热体电阻值和温度随时间变化试验不需要步骤四。
发热体由铂丝缠绕而成,随着通电时间的增加,铂丝温度不断升高,铂丝的电阻值随温度升高而增大。开展发热体电阻值及温度随时间变化试验以得出标准热源装置达到动态平衡(即发热体的电阻值、温度和功率几乎不随时间变化)的时间,以便确定后续试验中的时间间隔。试验结果如图4所示,可以看出随着时间的增长,热源温度和铂丝电阻值的变化率均越来越小,40 min后,温度变化范围在0.5℃以内,铂丝电阻值的变化范围在0.005Ω范围内,相对于10.651Ω的变化率为0.05%。
图4 发热体电阻值及温度随时间变化规律Fig.4 The resistance value and temperature change of heating element with time
后续发热体功率随电流变化试验中,设置恒流源的电流和电压后,需等待相同时间,以保证试验结果的重复性,即每个电流值条件下需要保持至少40 min再记录数字多用表的电压值。
本试验目的是验证功率随电流的变化规律及各型号发热体达到的功率上限。当本发热体作为标准器时,根据功率与电流的拟合曲线,找到指定功率下对应的电流值,试验结果如图5~图7所示。本试验中,每个电流值条件下保持1 h再记录数据;数据记录后设置下一电流值。
图5 0.7~10 W功率随电流变化Fig.5 0.7~10 W power varies with current
图7 100~147 W功率随电流变化Fig.7 100~147 W power varies with current
由图5和图6可知,在空气中自然对流的情况下,拟合曲线显示功率与电流呈幂函数关系。由图7可知,功率与电流呈线性关系。
图6 10~100 W功率随电流变化Fig.6 10~100 W power varies with current
发热体由铂丝缠绕而成,铂丝的阻值随温度的升高呈现幂函数关系[8],即
式中:Rt为实际电阻值,Ω;Rn为标称电阻,Ω;T为温度,℃。
联立式(2)和式(4),可以得到
发热体产生的热量一部分随自然对流换热而散失,另一部分使内能增加,则有
式中:P内为使内能增加的功率,W;P散为由于自然对流等散失的功率,W。
P内与温度之间的关系[9]可表示为
式中:t为时间,s;c为比热容,J/(kg·K);m为质量,kg;ΔT为温度变化量,K。
P散与温度之间的关系[10]可表示为
式中:h为对流换热系数,W/(m2·K);tw为壁面温度,K;tf为来流温度(室温温度),K;A为换热面积,m2。
h可根据努塞尔数Nu计算得出,Nu的表达式为
式中:λ为导热系数,W/(m·K);l为特征长度,m。
努塞尔数Nu与普朗特数Pr和格拉校夫数Gr的乘积相关。Pr的计算公式为
式中:ν为运动粘度系数,m2/s;a为热扩散系数,m2/s。
Gr的计算公式为
式中:aV为体膨胀系数,K-1;ΔT为壁面温度与无穷远处流体温度的温差(无穷远处流体温度为室温温度,是恒定的),K。
努塞尔数Nu与普朗特数Pr和格拉校夫数Gr的关系为
式中:m为定性温度,℃;C为常数;n为系数。C和n的值受流态和Gr影响。根据实际情况判断,自然对流流态为层流,Gr的范围为1.43×104~3×109,因此可得C=0.59,n=1/4[10]。
tm为定性温度采用边界层的算数平均温度,即
联立式(8)~式(13)可得
根据式(6)、式(7)和式(14),可以得到
联立式(5)和式(15),得到P与I的隐函数关系。根据试验数据拟合曲线,发现0.7~10 W和10~100 W的发热体产生的功率随电流呈现二次幂函数关系,100~147 W的发热体产生的功率随电流呈现线性关系。这一结果的产生与发热体的外形尺寸相关,产生0.7~10 W和10~100 W的发热体直径和高度相对较小,自然对流换热量相对较小;而产生100~147 W的发热体直径和高度相对较大,自然对流换热量相对较大。
综上可知,本文研制的发热体在电流小于等于3 A时,可产生0.7~147 W的功率。
本试验目的是验证功率的短期稳定性,共进行3次试验,测试结果如图8所示,功率稳定后的数据如表4所示,可知稳定后功率变化最大值为0.008 W,相对变化为0.08%。
图8 功率随时间变化重复试验Fig.8 Repeat test of power changes with time
表4 功率稳定后试验结果Tab.4 Test results after power stabilization
表5为电流间隔1 h更改一次的功率值,共进行3组试验,每组间隔24 h,结果分别为P11,P12,P13。由表5可知,相同电流下,最大功率差值为0.051 W,相对变化约为0.21%。
表5 功率与电流关系短期稳定性试验结果Tab.5 Short-term stability test results of the relationship between power and current
本装置的量程范围为0.7~147 W,其测量不确定度[11]来源主要包括:测量重复性引入的不确定度分量、标准电阻两端电压测量误差引入的不确定度分量、发热体两端电压测量误差引入的不确定度分量、扫描开关寄生电势引入的不确定度分量、标准电阻阻值引入的不确定度分量、标准电阻在空气中的阻值变化引入的不确定度分量、引线漏热引入的不确定度分量。
以2.98 W功率为例,计算各不确定度分量:
1)测量重复性引入的不确定度分量u1
进行4次重复试验,通过极差法评定不确定度,即
式中:R为极差,R=xmax-xmin,W;C为极差系数,查表得到,当n=4时,C=2.06。
由表6可知,功率P为2.98 W时,计算重复性引入的不确定度u1=s(xk)=0.0025 W。
表6 重复性Tab.6 Repeatability
2)标准电阻两端电压测量误差引入的不确定度分量u2
当功率为2.98 W时,标准电阻两端电压为100.47 mV。选择档位为1 V档位,根据校准证书可知,100 mV测量不确定度为5×10-6V(k=2),其引入的不确定度分量为
3)发热体两端电压测量误差引入的不确定度分量u3
当功率为2.98 W时,发热体两端的电压为2.973 V,使用数字电压表3458A测量时,应选择10 V档位,根据校准证书可知,3 V电压测量不确定度为9×10-6V(k=2),其引入的不确定度分量为
4)扫描开关寄生电势引入的不确定度分量u4
扫描开关的寄生电势小于0.2μV,按均匀分布计算,取半宽,扫描开关引入的的不确定度为
5)标准电阻阻值引入的不确定度分量u5
试验过程选择SRL-0.1Ω标准电阻,最大电流为3 A,标准电阻最大功率为0.3 W。以二等标准器使用,标准电阻误差为±1×10-5Ω,按照均匀分布计算,取半宽,标准电阻引入的不确定度分量为
6)标准电阻在空气中的阻值变化引入的不确定度分量u6
试验过程选择SRL-0.1Ω标准电阻,根据其说明书,当环境温度在18.0~28.0℃范围变化时,其阻值变化为0.10000044~0.09999945Ω,分散区间为0.1Ω±0.00000045Ω。在整个试验过程中标准电阻周围温度均在25.3℃左右,未超过28℃。因此标准电阻在空气中的阻值变化引入的不确定度分量为
7)引线漏热引入的不确定度分量u7
电流回路的引线和电压测量的引线均由发热体两端引出,因此加热丝产生的热功率会通过引线散失一部分。引线导热系数λ随温度的变化关系可以近似为λ=-0.0521t+399.53,根据傅里叶导热定律
式中:Φ为热流量,W;A为引线的截面积,m2;L为引线长度,m;t为温度,℃。
功率为2.98 W时,引线漏热引入的功率为1.26×10-3W,按照均匀分布,引线漏热引入的不确定度分量为
合成标准不确定度u计算公式为
相对扩展不确定度U为
功率为2.98 W时,对应的不确定度分量如表7所示。计算可得合成标准不确定度为2.7×10-3W,相对扩展不确定度为1.9×10-3(k=2)。
表7 2.98 W对应不确定度分量Tab.7 Uncertainty components of 2.98 W
标准热源装置0.7~147 W功率不确定度如表8所示,0.7~147 W对应的相对扩展不确定度不大于2×10-3(k=2)。
表8 不同功率对应不确定度Tab.8 Uncertainty under the condition of different power
目前我国缺乏准确度高、稳定性好的核电池热功率溯源用标准热源装置,为解决此问题,本文研制了一套标准热源装置。采用铂丝作为加热丝以提高热源稳定性,通过加热丝双螺旋外缠绕方式提高热功率准确度,利用石英玻璃材料制造发热体外壳、骨架、顶盖以保证热源的高可靠性。该装置能够在电流小于等于3 A的范围内产生0.7~147 W的功率,相对扩展不确定度不大于2×10-3(k=2),为保障核电池热功率测量数据的准确溯源提供了重要支撑,具有技术借鉴价值。