CN104215660B - 一种可同时测固体材料导热系数及热扩散率的方法及系统 - Google Patents
一种可同时测固体材料导热系数及热扩散率的方法及系统 Download PDFInfo
- Publication number
- CN104215660B CN104215660B CN201410415796.7A CN201410415796A CN104215660B CN 104215660 B CN104215660 B CN 104215660B CN 201410415796 A CN201410415796 A CN 201410415796A CN 104215660 B CN104215660 B CN 104215660B
- Authority
- CN
- China
- Prior art keywords
- thermal diffusivity
- sample
- heat conductivity
- temperature
- constant temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Abstract
本发明涉及一种可同时测固体材料导热系数及热扩散率的方法及系统,即利用超级恒温水浴作为热源,与试样盒底部的黄铜板形成恒温边界,温度均匀性和稳定性好,易于实现;而且直接基于恒温边界一维非稳态传热模型,利用参数估计法测固体材料导热系数及热扩散率,避免了对数学模型过度简化造成的误差。此外,根据比热容参数估计灵敏度较低,估计误差较大的情况,采取了先估计热扩散率,然后对比热容估计值进行校正的策略,实现了对材料导热系数、比热容和热扩散率的同时准确估计。
Description
技术领域
本发明属于固体材料热物性参数测试技术领域,具体涉及一种可同时测固体材料导热系数及热扩散率的方法及系统。
背景技术
随着我国对节能及能源利用效率的日益重视,材料热物性及其测试方法方面的研究方兴未艾。固体材料热物性测试总体上分为稳态法和非稳态法,稳态法测试由于所需时间长、测试环境要求高,而不能适应现代科技迅速发展的需要,非稳态法应运而生并因为其具有快速、准确、可实现多参数测量的优势而得到了更为广泛的应用。然而,非稳态法在实际应用中也存在一些问题,例如目前多采用加热板、加热丝、加热棒等固态电热器作为加热热源,存在着热源热流密度的准确给定、加热探头热容及时间延迟性、热丝电阻温度依变性和功率波动等问题,这些问题在一定程度上影响了热物性参数测试的精度。近年来,发展了一种陶瓷基底刻蚀双螺旋铂金属丝加热膜,热稳定性好且厚度极薄,可有效提高材料热物性测试精度,但成本不菲且制作困难,难以推广应用。因此,目前尚无一种精度理想的可同时测固体材料导热系数及热扩散率的方法及系统。
发明内容
本发明要解决的技术问题是提供一种可同时测固体材料导热系数及热扩散率的方法及系统,以克服目前不能同时测试固体材料导热系数和热扩散率及其测试精度不理想的缺陷。
为解决上述技术问题,本发明采用以下技术方案:
一种可同时测固体材料导热系数及热扩散率的方法,具体实现步骤为:
步骤1:通过超级恒温水浴产生具有一定温度的恒温水介质;
步骤2:将试样放入底部为黄铜板材质的试样盒,试样盒底部与恒温水介质接触,形成恒温边界下的一维非稳态传热过程,试样内温度响应解析解为:
式中θ(x,τ)为过余温度,θ=t-tw,t为试样温度,tw为黄铜板恒温温度;热扩散率a=λ/(ρCp),λ为导热系数,ρ为试样密度、Cp为比热。
式(1)为无穷级数,无法直接求解;而由于恒温边界条件等价于毕渥数Bi→∞时的热对流边界条件,从而式(1)可改写为:
式中βn为超越方程ctgβ=β/Bi的系列解,对应不同的Bi数,βn值也不同(具体可查表),而在Bi→∞时β1=1.5708,βn=β1+(n-1)×π(n=1,2,…),显然在试样密度ρ、导热系数λ、比热容Cp及厚度δ物理参数已知的情况下,可计算得到试样内任意位置x任意时刻τ的温度分布情况,此即恒温边界条件下的非稳态传热数学模型;
步骤3:在该数学模型的基础上,利用参数估计法得到被测材料的热扩散率、比热容和导热系数。
优选的,所述参数估计法采用先估计热扩散率,然后对比热容估计结果进行修正,实现被测试样导热系数、热扩散率和比热容多个参数同时准确获取。
本发明还提供一种可同时测固体材料导热系数及热扩散率方法的系统,包括超级恒温水浴、试样盒和温度采集与分析系统,所述超级恒温水浴,用于形成恒温边界条件;所述试样盒,用于盛放试样,与恒温水介质接触形成恒温边界,可测定除板状成型材料外的松散物料;所述温度采集与分析系统,用于采集试样绝热侧表面温度变化,以基于数学模型计算得到被测试样的导热系数及热扩散率。
优选,所述超级恒温水浴生成具有一定温度的恒温水介质,通过与试样盒底部接触形成恒温边界。
优选的,所述试样盒的内腔尺寸满足长与厚的比值大于8~10,宽与厚度的比值大于8~10。
本发明与现有技术相比所取得的有益效果为:利用超级恒温水浴作为热源,与试样盒底部的黄铜板形成恒温边界,温度均匀性和稳定性好,且易于实现;而且直接基于恒温边界一维非稳态传热模型,利用参数估计法测固体材料导热系数及热扩散率,避免了对数学模型过度简化造成的误差。此外,根据比热容参数估计灵敏度较低,估计误差较大的情况,采取了先估计热扩散率,然后对比热容估计值进行校正的策略,实现了对材料导热系数、比热容和热扩散率的同时准确估计。
附图说明
1.图1为本发明的硬件结构原理图;
2.图2为本发明计算原理模型;
图1和图2中,1、水介质2、加热元件3、温度传感器4、温控单元5、搅拌器6、绝热水浴箱7、底部黄铜板8、侧绝热板9、上绝热板10、试样11、热电偶信12、信号采集与传输单元13、微型计算机。
3.图3和图4为本发明待估参数灵敏度系数计算曲线图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步详细的说明。
如图1所示,水介质1温度由超级恒温水浴系统控制,超级恒温水浴系统包括加热元件2、温度传感器3、温控单元4、搅拌器5、绝热水浴箱6等组成,其中加热元件2用于加热水介质1,温度传感器3用于检测水介质1的实时温度,并反馈给温控单元4,控制加热元件2的启闭,以保证水介质1维持恒定温度。搅拌器5用以搅拌水介质1,以使其温度均匀,绝热层水浴箱6用于盛放水介质1,并使之与周围环境绝热,降低热量散发因素的影响。试样盒由底部黄铜板7、侧绝热板8和上绝热板9组成,当其底部黄铜板7与恒温的水介质1接触时,将形成恒温边界。如在某一时刻τ突然将厚度为δ,初始温度为t0的试样10与黄铜板7接触,试样10中将被加热,而产生恒温边界条件下的一维热传导过程。当试样盒的内腔尺寸满足长和宽大于厚度的8-10倍时,能更好地形成一维非稳态传热条件,使测得的结果精度更高。试样10被加热过程中,上绝热面处的温度变化通过热电偶11检测,并通过信号采集与传输单元12传送给微型计算机13,进而利用数学模型计算得到试样热物性参数。
如图2所示,设试样导热系数为λ,密度ρ、比热Cp等均为常数,其初始温度为t0,坐标原点取于试样上绝热面,则试样10与恒温T=tw的黄铜板7接触后,其内的非稳态导热微分方程及其初始、边界条件如下:
τ=0,t=t0(2a)
x=0,
x=δ,t=tw(2c)
式(1)的解析解如下:
式中θ(x,τ)为过余温度,θ=t-tw;t为试样温度,tw为黄铜板恒温温度;热扩散率a=λ/(ρCp),λ为导热系数,ρ为试样密度、Cp为比热。
式(3)为无穷级数,无法直接求解;而由于恒温边界条件等价于毕渥数Bi→∞时的热对流边界条件,从而式(3)可改写为:
式中βn为超越方程ctgβ=β/Bi的系列解,对应不同的Bi数,βn值也不同(具体可查表),而在Bi→∞时β1=1.5708,βn=β1+(n-1)×π(n=1,2,…),试算表明,针对固体材料热物性测试需求取前15项的计算误差即可控制在1%以内,显然在试样10的密度ρ、导热系数λ、比热容Cp及厚度δ等物理参数已知的情况下,可计算得到试样10内任意位置x任意时刻τ的温度分布情况,此即恒温边界条件下的非稳态传热数学模型。
在式(4)的基础上,且试样热物性参数未知的情况下,由测得的x=0处的温度数据进行反演推算,也称参数估计。如前所述,从某一时刻开始试样下表面被温度为tw的热板加热,该温度扰动为输入信号,记为u(τ),试样内任意一点处的温升记为Y(τ,η),其不但与试样导热系数、比热容等热物性参数有关,还取决于热板温度tw、时间、空间位置等,所有这些参数构成一个向量:
η=(η1,η2,η3,......,ηm)T(5)
其中m为参数个数。
此处试样导热系数λ、比热容Cp未知,密度ρ可实测获得,首先获得x=0表面处k(k>2)个离散测点处的温度测量值Yi(τ,η)(i=1…k),然后基于估计的参数值,通过式(4)计算各离散测点的温度值θi(τ,η)(i=1…k),并进行对比:
由于测量误差以及参数估计误差的存在,ε(η)>0,估计试样导热系数及热扩散率的反问题可归结为对式(4)进行迭代计算,最终使目标函数ε(η)→min。
参数向量η中的j个参数要能同时以足够的精度估计出来,必须满足在最小二乘估计意义下,测量时间范围内(参数估计的时间区间内),参数的灵敏度线性无关。灵敏度系数反映了参数βj的微小扰动对输出的影响。
式(4)中待估计参数为:
η1=λ,η2=Cp
为检验对不同导热系数大小的材料参数估计可行性,分别选取导热系数较大的玻璃和导热系数较小的松散煤体进行参数灵敏度分析。设玻璃物理参数为:厚度5mm,密度2530kg/m3,导热系数0.701W/m.℃,比热容750J/kg.℃。松散煤体物理参数为:密度1010kg/m3,试样厚度60mm,粒度范围0-6mm,导热系数0.14W/m.℃,比热容840J/kg.℃。
灵敏度系数计算结果分别如图3和图4所示。从图可以看出,玻璃和松散煤体的导热系数和比热容参数均灵敏度线性无关,但比热容Cp灵敏度极小(仅10e-3数量级),意味着直接对其进行估计误差将比较大。实际上在测得温度θ(x,τ)的情况下,式(4)中的λ/(ρCp)项即热扩散率a,将只有此一个未知数,可首先将热扩散率a求解出来,然后利用其对估计的比热容Cp进行修正。也即,本系统首先采用随机共轭梯度法估计热扩散率,然后再对比热容估计结果进行修正,计算过程如下:
1)估计热扩散率a
(1)根据先验信息确定待反演参数的初始猜测值η0(此时参数向量η=(a)T)、设定一较小的大于零的数ζ作为退出控制误差。
(2)修改模型中的参数向量ηj:
式中ηj new和ηj old分别为修改前后的第j个参数;dηj表示第j个参数的基本偏离值,通过手工设定,用于控制模型参数的修改量;sign为符号函数,表示在相同概率下(0.5的概率)取值为1和-1;Rand()产生[0,1]范围内的随机数。
(3)将参数修改前后的目标函数分别记为εold(η)和εnew(η),其差值记为:△ε=εold(η)-εnew(η),当△ε<0时,p=1;△ε>0时,p=0;
(4)当p=1时,接受修改后的模型,计算共轭系数、梯度、搜索方向及搜索步长等参量;共轭梯度法收敛后,保存反演估计值;
(5)充分步骤(2)、(3)和(4),比较每次接受的反演估计值,直到收敛到允许误差范围。
2)将待估参数向量更改为η=(λ,Cp)T,然后按照步骤(1)至(5)估计导热系数λ及比热容Cp。
3)根据热扩散率a及导热系数λ估计值修正比热容Cp估计值。
理论上每个特定时刻τ均可得到一组数值λ和Cp,对计算出的一系列值求加权平均作为最终参数估计值,同时可消除测试过程中因干扰因素造成的个别数据波动影响。
以上实例所假设的参数和尺寸都是帮助理解本发明的方法和核心思想,针对不同的固体材料应该查相关手册来确定其参数以及根据实际情况来测量其尺寸大小。综上所述,本说明书参数和尺寸不应理解为对本发明的限制,在不脱离本发明的精神和范围情况下,还可以做出各种变化,因此所有等同的技术方案也属于本发明的范畴,本发明的权利保护范围应由权利要求限定。
Claims (5)
1.一种可同时测固体材料导热系数及热扩散率的方法,其特征在于:具体实现步骤为:
步骤1:通过超级恒温水浴产生具有一定温度的恒温水介质;
步骤2:将试样放入底部为黄铜板、上表面和侧面用绝热材料密封的试样盒,试样盒底部与恒温水介质接触,形成恒温边界下的一维非稳态传热过程,试样内温度响应解析解为:
式中θ(x,τ)为过余温度,θ=t-tw,t为试样温度,tw为黄铜板恒温温度;热扩散率a=λ/(ρCp),λ为导热系数,ρ为试样密度、Cp为比热容;
式(1)为无穷级数,无法直接求解;而由于恒温边界条件等价于毕渥数Bi→∞时的热对流边界条件,从而式(1)可改写为:
式中βn为超越方程ctgβ=β/Bi的系列解,对应不同的Bi数,βn值具体可查表,而在Bi→∞时β1=1.5708,βn=β1+(n-1)×π,n=1,2,…,显然在试样密度ρ、导热系数λ、比热容Cp及厚度δ物理参数已知的情况下,可计算得到试样内任意位置x任意时刻τ的温度分布情况,此即恒温边界条件下的非稳态传热数学模型;
步骤3:在该数学模型的基础上,利用参数估计法得到被测材料的热扩散率、比热容和导热系数。
2.根据权利要求1所述的可同时测固体材料导热系数及热扩散率的方法,其特征在于:所述参数估计法采用先估计热扩散率,然后对比热容估计结果进行修正,实现被测试样导热系数、热扩散率和比热容多个参数同时准确获取。
3.一种根据权利要求1~2之一所述的可同时测固体材料导热系数及热扩散率方法的系统,其特征在于:包括超级恒温水浴、试样盒和温度采集与分析系统,所述超级恒温水浴,用于形成恒温边界条件;所述试样盒,用于测定除板状成型材料外的松散物料,与恒温水介质接触形成恒温边界;所述温度采集与分析系统,用于采集试样绝热侧表面温度变化,以基于数学模型计算得到被测试样的导热系数及热扩散率。
4.根据权利要求3中所述的系统,其特征在于:所述超级恒温水浴生成具有一定温度的恒温水介质,通过与试样盒底部接触形成恒温边界。
5.根据权利要求3所述的系统,其特征在于:所述试样盒的内腔尺寸满足长与厚的比值大于8~10,宽与厚度的比值大于8~10。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410415796.7A CN104215660B (zh) | 2014-08-21 | 2014-08-21 | 一种可同时测固体材料导热系数及热扩散率的方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410415796.7A CN104215660B (zh) | 2014-08-21 | 2014-08-21 | 一种可同时测固体材料导热系数及热扩散率的方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104215660A CN104215660A (zh) | 2014-12-17 |
CN104215660B true CN104215660B (zh) | 2016-04-13 |
Family
ID=52097378
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410415796.7A Expired - Fee Related CN104215660B (zh) | 2014-08-21 | 2014-08-21 | 一种可同时测固体材料导热系数及热扩散率的方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104215660B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203593A (zh) * | 2015-10-14 | 2015-12-30 | 东南大学 | 基于混合料细观结构特性的沥青和集料导热系数反演方法 |
CN109000828A (zh) * | 2018-08-23 | 2018-12-14 | 山东建筑大学 | 一种用于非均质材料蓄热量的检测装置及其检测计算方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106053527B (zh) * | 2016-07-04 | 2019-02-26 | 中国第一汽车股份有限公司 | 一种同时测试动力电池各向异性热扩散系数的方法 |
CN108051475B (zh) * | 2017-12-05 | 2020-08-14 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种对流换热系数的快速测量方法 |
CN113138207B (zh) * | 2021-04-22 | 2022-04-19 | 安徽理工大学 | 一种正交各向异性固体材料热扩散系数测试系统及方法 |
CN114199935B (zh) * | 2021-12-28 | 2024-02-09 | 浙江大学 | 一种海洋沉积物导热系数测量方法及装置 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4195935B2 (ja) * | 2004-03-01 | 2008-12-17 | 独立行政法人産業技術総合研究所 | 熱物性測定方法及び装置 |
CN101339149B (zh) * | 2008-07-02 | 2010-09-15 | 河海大学 | 测定早期混凝土导热系数、导温系数的测试装置及测试方法 |
CN101639455A (zh) * | 2009-09-04 | 2010-02-03 | 陈清华 | 松散煤体热物性测试方法及系统 |
CN202189027U (zh) * | 2011-08-29 | 2012-04-11 | 舟山市博远科技开发有限公司 | 混凝土热扩散率测定装置 |
CN103411996B (zh) * | 2013-08-05 | 2016-03-02 | 电子科技大学 | 固体材料导热系数测量装置及测量方法 |
CN103954650B (zh) * | 2014-04-30 | 2016-06-15 | 安徽理工大学 | 一种固体材料热扩散系数测试方法与系统 |
-
2014
- 2014-08-21 CN CN201410415796.7A patent/CN104215660B/zh not_active Expired - Fee Related
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105203593A (zh) * | 2015-10-14 | 2015-12-30 | 东南大学 | 基于混合料细观结构特性的沥青和集料导热系数反演方法 |
CN105203593B (zh) * | 2015-10-14 | 2018-02-02 | 东南大学 | 一种基于混合料细观结构特性的沥青和集料导热系数反演方法 |
CN109000828A (zh) * | 2018-08-23 | 2018-12-14 | 山东建筑大学 | 一种用于非均质材料蓄热量的检测装置及其检测计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104215660A (zh) | 2014-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104215660B (zh) | 一种可同时测固体材料导热系数及热扩散率的方法及系统 | |
CN101126729B (zh) | 双热流计稳态法材料热导率测量方法 | |
Tisserand et al. | Comparison between rough and smooth plates within the same Rayleigh–Bénard cell | |
EP2126533A1 (en) | Apparatuses and methods for measuring and controlling thermal insulation | |
CN101113963A (zh) | 一种测量液体导热系数的方法及其装置 | |
CN109752113B (zh) | 薄片温度传感器及其应用中位置确定方法和电路设计方法 | |
CN103940847B (zh) | 一种基于热流传感器各向异性薄膜热导率测试方法及装置 | |
CN106338323B (zh) | 一种流量计现场检定方法及检定系统 | |
CN109443601A (zh) | 基于热电阻的高精度多点温度测量系统及温度校准方法 | |
CN103954650B (zh) | 一种固体材料热扩散系数测试方法与系统 | |
CN109100051B (zh) | 温度传感器的动态响应的温度修正方法及装置 | |
CN102879421A (zh) | 导热系数的测量装置及其测量方法 | |
CN109738092B (zh) | 一种双向热电堆式薄膜热流计及热流测量方法 | |
CN204086184U (zh) | 一种可同时测固体材料导热系数及热扩散率的系统 | |
CN103257154B (zh) | 大温差样品的半球向全发射率与导热系数的测量方法 | |
CN108051475A (zh) | 一种对流换热系数的快速测量方法 | |
CN105572161B (zh) | 一种非恒态导热系数测试的方法及所用测试装置 | |
Jamroz | Relationship between dynamic coefficients of two temperature sensors under nonstationary flow conditions | |
CN107966472B (zh) | 一种高温接触热阻的无损快速测量方法 | |
CN108956688B (zh) | 一种建筑材料表面综合辐射吸收系数的测量系统及方法 | |
Cortellessa et al. | Experimental and numerical analysis in heat flow sensors calibration | |
CN103245692B (zh) | 基于稳态分析的半球向全发射率与导热系数的测量方法 | |
CN108776152A (zh) | 一种建筑材料表面对流热交换系数的测量系统及方法 | |
CN109506806B (zh) | 一种瞬态条件下高温结构内部温度及厚度的同时测量方法 | |
CN202421102U (zh) | 测定气体导热系数的装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160413 Termination date: 20170821 |