CN105319174B - 同时获取半透明材料温变导热系数及吸收系数的测量方法 - Google Patents
同时获取半透明材料温变导热系数及吸收系数的测量方法 Download PDFInfo
- Publication number
- CN105319174B CN105319174B CN201510907097.9A CN201510907097A CN105319174B CN 105319174 B CN105319174 B CN 105319174B CN 201510907097 A CN201510907097 A CN 201510907097A CN 105319174 B CN105319174 B CN 105319174B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- temperature
- msubsup
- 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.)
- Active
Links
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
同时获取半透明材料温变导热系数及吸收系数的测量方法,涉及同时获取半透明介质温度相关导热系数及吸收系数技术。测量过程中使用某波长连续激光照射待测样品,借助探测器测量待测样本的随时间变化的温度响应以及透射辐射强度,最后通过逆问题求解技术间接得到待测样品随温度变化的导热系数及吸收系数。本发明通过建立导热系数及吸收系数随温度变化的半透明介质导热辐射耦合换热的正、反问题模型,在介质其他参数已知的前提下,提出了采用微粒群优化算法同时反演得到半透明介质温度相关导热系数及吸收系数的方法。本发明适用于航天、国防和民用工业。
Description
技术领域
本发明涉及同时获取半透明介质温度相关导热系数及吸收系数技术,属于半透明介质物性测量技术领域。
背景技术
半透明介质辐射物性以及热物性参数是对半透明介质在其应用过程中进行分析、设计、优化所需的重要参数。近年来,随着航空航天、红外探测、目标与环境的红外特性、激光、电子器件、生物医学等现代高新技术的飞速发展,半透明介质在高温、多维等情况下的随温度变化物性参数变得尤为重要。进行参与性介质热辐射物性以及相关学科的研究对于军用和民用领域均具有重要意义。
在探测器光学窗口以及光致发光材料等研究领域,对于纯吸收介质的导热系数和吸收系数的研究显得尤为重要。深入理解此热物性参数并对其进行实验测量及理论分析在材料科学以及环境监测等领域也具有重要的应用价值。并且通常情况下,导热系数以及吸收系数是与材料温度相关的。因此,对于随温度变化的导热系数以及吸收系数的测量在实际应用过程中将具有重要意义。
由于实际测量过程中,实验设备存在一定的测量误差,某些情况下单独使用光或者热信息不能完成辐射热物性的测量或者获得的结果误差较大,并且对于温度相关热物性的反演需要更多的测量信息。
发明内容
本发明是为提高对半透明介质热物性测量的精度,从而提供一种同时获取半透明材料温变导热系数及吸收系数的测量方法。
同时获取半透明材料温变导热系数及吸收系数的测量方法,它由以下步骤实现:
步骤一、制作厚度为L的待测样品;
步骤二、利用波长为λ的连续激光沿着与厚度为L的待测样品表面垂直的方向入射到待测样品左侧表面,持续时间为t秒;采用探测器在待测样品的右侧表面分别测量其随时间变化的温度Tw(t)以及辐射强度R(t);
步骤三、利用逆问题算法假设出待测样品的对应波长随温度T变化的吸收系数κa(T)=a1+a2·T mm-1和随温度变化的导热系数λ(T)=b1+b2·TW/(m·K);式中a1,a2,b1和b2表示需要获取的系数。
然后通过对辐射传输方程以及导热微分方程的求解,获得计算域内的辐射强度场及温度场;同时得到待测样品右侧随时间变化的温度的预测值Tw,est(t);
步骤四、利用步骤三获得的辐射强度场结合下列公式:
获得右侧边界的辐射强度的预测值Rest(t);
式中:I0,λ是波长为λ的连续激光的强度;Iλ(L,θ)为θ方向上z=L处的右侧边界上散射光的辐射强度,θ为辐射方向角;Ic,λ(L,θc)为连续激光沿着入射方向θc衰减到样本右侧壁面时的辐射强度,θc为连续激光入射方向角,此处θc=0;
步骤五、利用步骤二获得的右侧边界处的随时间变化的温度Tw(t)以及辐射强度R(t)与步骤三中相应的预测值,结合公式:
获得逆问题算法中的目标函数F1,obj;式中:t1和t2为温度(或辐射强度)的测量时间。
步骤六、判断步骤五中的目标函数是否小于设定阈值ε1,若是,则将步骤三中假设的待测样品的导热系数λ(T)=b1+b2·T W/(m·K)作为结果输出,否则返回步骤三重新修正预测的导热系数及吸收系数;
步骤七、重复步骤三和四,其中导热系数使用步骤六输出的结果;
步骤八、利用步骤二中获得的右侧边界处的辐射强度R(t)与步骤四中相应的预测值,结合公式:
获得逆问题算法中的目标函数F2,obj;
步骤九、判断步骤八中的目标函数是否小于设定阈值ε2,若是,则将步骤七中获得的待测样品的吸收系数κa(T)=a1+a2·T mm-1作为结果输出,完成基于同时获取半透明介质温度相关导热系数及吸收系数的方法,否则,返回步骤七。
本发明提出的测量方法在逆问题求解的基础上引入了光热信息融合技术,能够大大提高对于半透明介质热物性测量的精度。本发明通过建立导热系数及吸收系数随温度变化的半透明介质导热辐射耦合换热的正问题和逆问题求解模型,解决半透明介质随温度变化的导热系数和吸收系数不能直接测量和测量结果不准确的问题,提出了一种同时获取半透明介质温度相关导热系数及吸收系数的方法。优点在于:采用连续激光,该激光器廉价购买方便,且模型简单,便于理论求解;采用量子微粒群优化算法,该算法求解优化问题时有简单、高效和灵敏度高等优点。该项发明为研究半透明介质随温度变化的导热系数和吸收系数提供一种快速准确的方法,对航天、国防和民用工业具有十分重要的意义。
附图说明
图1是具体实施方式一所述连续激光辐照下导热系数和吸收系数随温度变化的半透明介质辐射导热耦合模型示意图;图中左侧实心箭头为连续激光入射方向,左、右侧的空心箭头方向为辐射热流方向。
具体实施方式
具体实施方式一、结合图1说明本具体实施方式,同时获取半透明材料温变导热系数及吸收系数的测量方法,该方法的具体操作步骤为:
步骤一、制作厚度为L的待测样品;
步骤二、如图1所示,利用波长为λ的连续激光沿着与厚度为L的样本表面垂直的方向入射到待测样本左侧表面,持续时间为t秒;使用探测器在样本的右侧表面分别测量其随时间变化的温度Tw(t)以及辐射强度R(t);
步骤三、利用逆问题求解思路假设出待测样品的对应波长随温度变化的吸收系数κa(T)=a1+a2·T mm-1和随温度变化的导热系数λ(T)=b1+b2·T W/(m·K);然后通过对辐射传输方程以及导热微分方程的求解,获得计算域内的辐射强度场及温度场;同时可以得到样品右侧随时间变化的温度的预测值Tw,est(t);
步骤四、利用步骤三获得的辐射强度场结合下列公式:
获得右侧边界的辐射强度的预测值Rest(t)。式中I0,λ是波长为λ的连续激光的强度;θ为天顶角;Iλ(L,θ)为θ方向上z=L处的右侧边界上散射光的辐射强度,θ为辐射方向角;Ic,λ(L,θc)为连续激光沿着入射方向θc衰减到样本右侧壁面时的辐射强度,θc为连续激光入射方向角,此处θc=0;
步骤五、利用步骤二获得的右侧边界处的随时间变化的温度Tw(t)以及辐射强度R(t)与步骤三中相应的预测值,结合公式:
获得逆问题算法中的目标函数F1,obj。
步骤六、判断步骤五中的目标函数是否小于设定阈值ε1,若是,则将步骤三中假设的待测样品的导热系数λ(T)=b1+b2·T W/(m·K)作为结果输出,否则返回步骤三重新修正预测的导热系数及吸收系数。
步骤七、重复步骤三和四,其中导热系数不需要重复假设,而是使用步骤六输出的结果。
步骤八、利用步骤二中获得的右侧边界处的辐射强度R(t)与步骤四中相应的预测值,结合公式:
获得逆问题算法中的目标函数F2,obj;
步骤九、判断步骤八中的目标函数是否小于设定阈值ε2,若是,则将步骤七中获得的待测样品的吸收系数κa(T)=a1+a2·T mm-1作为结果输出,完成基于同时获取半透明介质温度相关导热系数及吸收系数的方法,否则,返回步骤七;
本实施方式首先设计导热系数和吸收系数随温度变化的半透明介质内瞬态辐射导热耦合物理模型,然后建立相应的数学模型和求解方法,通过测量得到待测样品的随时间变化的温度以及辐射强度,利用逆问题理论模型的重建出半透明介质的温度相关的吸收系数和导热系数。
具体实施方式二、本实施方式是对具体实施方式一所述的同时获取半透明介质温度相关导热系数及吸收系数的方法的进一步说明,步骤三获得计算域内的温度场的方法为:
利用导热微分方程:
T(t=0)=T0 (5)
实现,其中ρ和cp分别表示待测介质的密度及比热容,z表示待测样品厚度方向坐标,λ表示待测介质的导热系数,T和h分别表示温度和对流换热系数。qr表示热流密度,其中脚标w1和w2分别表示待测样品的左边界和右边界。
具体实施方式三、本实施方式是对具体实施方式一所述的同时获取半透明介质温度相关导热系数及吸收系数的方法的进一步说明,步骤三获得计算域内的辐射场强度的方法为:
利用辐射传输方程:
实现,式中κa表示待测介质的吸收系数,I表示介质内辐射强度,Ib表示相同温度下黑体的辐射强度。
具体实施方式四、本实施方式是对具体实施方式二所述的获得计算域内的温度场的方法的进一步说明,获取导热微分方程中的热流密度的方法为:利用方程
实现,式中ε1和ε2分别表示待测介质两侧壁面的发射率,σ表示斯特藩-玻尔兹曼常数。
Claims (4)
1.同时获取半透明材料温变导热系数及吸收系数的测量方法,其特征是:它由以下步骤实现:
步骤一、制作厚度为L的待测样品;
步骤二、利用波长为λ的连续激光沿着与厚度为L的待测样品表面垂直的方向入射到待测样品左侧表面,持续时间为t秒;采用探测器在待测样品的右侧表面分别测量其随时间变化的温度Tw(t)以及辐射强度R(t);
步骤三、利用逆问题算法假设出待测样品的对应波长随温度T变化的吸收系数κa(T)=a1+a2·T mm-1和随温度变化的导热系数λ(T)=b1+b2·T W/(m·K);式中a1,a2,b1和b2表示需要获取的系数;
然后通过对辐射传输方程以及导热微分方程的求解,获得计算域内的辐射强度场及温度场;同时得到待测样品右侧随时间变化的温度的预测值Tw,est(t);
步骤四、利用步骤三获得的辐射强度场结合下列公式:
<mrow>
<msub>
<mi>R</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>I</mi>
<mrow>
<mn>0</mn>
<mo>,</mo>
<mi>&lambda;</mi>
</mrow>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mn>2</mn>
<mi>&pi;</mi>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mrow>
<mi>&pi;</mi>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msubsup>
<msub>
<mi>I</mi>
<mi>&lambda;</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mi>d</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>c</mi>
<mo>,</mo>
<mi>&lambda;</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>,</mo>
<msub>
<mi>&theta;</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
获得右侧边界的辐射强度的预测值Rest(t);
式中:I0,λ是波长为λ的连续激光的强度;Iλ(L,θ)为θ方向上z=L处的右侧边界上散射光的辐射强度,z表示待测样品厚度方向坐标,θ为辐射方向角;Ic,λ(L,θc)为连续激光沿着入射方向θc衰减到样本右侧壁面时的辐射强度,θc为连续激光入射方向角,此处θc=0;
步骤五、利用步骤二获得的右侧边界处的随时间变化的温度Tw(t)以及辐射强度R(t)与步骤三中相应的预测值,结合公式:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mn>1</mn>
<mo>,</mo>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msubsup>
<mo>|</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<msub>
<mi>T</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1.0</mn>
<mo>|</mo>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
获得逆问题算法中的目标函数F1,obj;式中:t1和t2为温度或辐射强度的测量时间;
步骤六、判断步骤五中的目标函数是否小于设定阈值ε1,若是,则将步骤三中假设的待测样品的导热系数λ(T)=b1+b2·T W/(m·K)作为结果输出,否则返回步骤三重新修正预测的导热系数及吸收系数;
步骤七、重复步骤三和四,其中导热系数使用步骤六输出的结果;
步骤八、利用步骤二中获得的右侧边界处的辐射强度R(t)与步骤四中相应的预测值,结合公式:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mn>2</mn>
<mo>,</mo>
<mi>o</mi>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<msub>
<mi>t</mi>
<mn>1</mn>
</msub>
<msub>
<mi>t</mi>
<mn>2</mn>
</msub>
</msubsup>
<mo>|</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>e</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>R</mi>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mn>1.0</mn>
<mo>|</mo>
<mi>d</mi>
<mi>t</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
获得逆问题算法中的目标函数F2,obj;
步骤九、判断步骤八中的目标函数是否小于设定阈值ε2,若是,则将步骤七中获得的待测样品的吸收系数κa(T)=a1+a2·T mm-1作为结果输出,完成基于同时获取半透明介质温度相关导热系数及吸收系数的方法,否则,返回步骤七。
2.根据权利要求1所述的同时获取半透明材料温变导热系数及吸收系数的测量方法,其特征在于步骤三获得计算域内的温度场的方法为:
利用导热微分方程:
<mrow>
<msub>
<mi>&rho;c</mi>
<mi>p</mi>
</msub>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>t</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mi>&lambda;</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>)</mo>
</mrow>
<mfrac>
<mrow>
<msup>
<mo>&part;</mo>
<mn>2</mn>
</msup>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<msup>
<mi>z</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>q</mi>
<mi>r</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
T(t=0)=T0 (5)
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mn>1</mn>
</mrow>
<mi>r</mi>
</msubsup>
<mo>-</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mrow>
<mi>z</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mn>1</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>&infin;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mn>2</mn>
</mrow>
<mi>r</mi>
</msubsup>
<mo>-</mo>
<mi>&lambda;</mi>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>T</mi>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<msub>
<mo>|</mo>
<mrow>
<mi>z</mi>
<mo>=</mo>
<mi>L</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>h</mi>
<mrow>
<mi>w</mi>
<mn>2</mn>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>T</mi>
<mi>&infin;</mi>
</msub>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
实现,其中ρ和cp分别表示待测介质的密度及比热容,z为待测样品厚度方向坐标,λ表示待测介质的导热系数,T和h分别表示温度和对流换热系数;qr表示热流密度,其中脚标w1和w2分别表示待测样品的左边界和右边界。
3.根据权利要求1所述的同时获取半透明材料温变导热系数及吸收系数的测量方法,其特征在于,步骤三获得计算域内的辐射场强度的方法为:
利用辐射传输方程:
<mrow>
<mfrac>
<mrow>
<mi>d</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mi>d</mi>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&kappa;</mi>
<mi>a</mi>
</msub>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>&kappa;</mi>
<mi>a</mi>
</msub>
<msub>
<mi>I</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
实现,式中κa表示待测介质的吸收系数,I表示介质内辐射强度,Ib表示相同温度下黑体的辐射强度。
4.根据权利要求2所述的同时获取半透明材料温变导热系数及吸收系数的测量方法,其特征在于,获取导热微分方程中的热流密度的方法为:
利用方程:
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mn>1</mn>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mn>1</mn>
</msub>
<mo>&lsqb;</mo>
<msubsup>
<mi>&sigma;T</mi>
<mrow>
<mi>w</mi>
<mn>1</mn>
</mrow>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo><</mo>
<mn>0</mn>
</mrow>
</msub>
<mn>2</mn>
<mi>&pi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mn>0</mn>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mi>d</mi>
<mi>&theta;</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msubsup>
<mi>q</mi>
<mrow>
<mi>w</mi>
<mn>2</mn>
</mrow>
<mi>r</mi>
</msubsup>
<mo>=</mo>
<msub>
<mi>&epsiv;</mi>
<mn>2</mn>
</msub>
<mo>&lsqb;</mo>
<msubsup>
<mi>&sigma;T</mi>
<mrow>
<mi>w</mi>
<mn>2</mn>
</mrow>
<mn>4</mn>
</msubsup>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mrow>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
</msub>
<mn>2</mn>
<mi>&pi;</mi>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>L</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>|</mo>
<mi>c</mi>
<mi>o</mi>
<mi>s</mi>
<mi>&theta;</mi>
<mo>|</mo>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mi>d</mi>
<mi>&theta;</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfrac>
<mrow>
<mo>&part;</mo>
<msup>
<mi>q</mi>
<mi>r</mi>
</msup>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>z</mi>
</mrow>
</mfrac>
<mo>=</mo>
<mn>4</mn>
<mi>&pi;</mi>
<mo>&CenterDot;</mo>
<msub>
<mi>&kappa;</mi>
<mi>a</mi>
</msub>
<mo>&CenterDot;</mo>
<mo>&lsqb;</mo>
<msub>
<mi>I</mi>
<mi>b</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<msubsup>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>&pi;</mi>
</msubsup>
<mi>I</mi>
<mrow>
<mo>(</mo>
<mi>z</mi>
<mo>,</mo>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>&theta;</mi>
<mo>&rsqb;</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
实现,式中ε1和ε2分别表示待测介质两侧壁面的发射率,σ表示斯特藩-玻尔兹曼常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510907097.9A CN105319174B (zh) | 2015-12-09 | 2015-12-09 | 同时获取半透明材料温变导热系数及吸收系数的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510907097.9A CN105319174B (zh) | 2015-12-09 | 2015-12-09 | 同时获取半透明材料温变导热系数及吸收系数的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105319174A CN105319174A (zh) | 2016-02-10 |
CN105319174B true CN105319174B (zh) | 2018-01-30 |
Family
ID=55247037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510907097.9A Active CN105319174B (zh) | 2015-12-09 | 2015-12-09 | 同时获取半透明材料温变导热系数及吸收系数的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105319174B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105675646B (zh) * | 2016-03-23 | 2018-04-24 | 哈尔滨工业大学 | 基于本征光热信息同时测量高温半透明介质热导率及吸收系数的方法 |
CN106442347B (zh) * | 2016-09-09 | 2019-03-12 | 哈尔滨工业大学 | 高温熔融盐光谱吸收系数测量装置及测量方法 |
CN108333213B (zh) * | 2018-01-11 | 2020-04-24 | 哈尔滨工业大学 | 半透明多孔材料高温传导及辐射性质多参数同步测量方法 |
CN108362733B (zh) * | 2018-02-11 | 2020-06-09 | 哈尔滨工业大学 | 基于锁相热波与光学层析相结合的半透明材料光热特性分布测量方法 |
CN110657893B (zh) * | 2018-06-28 | 2020-12-08 | 上海梅山钢铁股份有限公司 | 一种水蒸气介质环境下辐射测温比色波长确定方法 |
CN110160964B (zh) * | 2019-05-24 | 2021-11-19 | 哈尔滨工业大学 | 基于Chirp调频激光辐照的半透明材料光热特性参数检测方法 |
CN113218872B (zh) * | 2021-04-08 | 2022-05-27 | 北华航天工业学院 | 一种高温半透明材料光学特性多参数同时辨识方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3926571B2 (ja) * | 2000-08-17 | 2007-06-06 | 独立行政法人科学技術振興機構 | 固体の温度伝導率等測定方法および装置 |
CN1584521A (zh) * | 2004-06-04 | 2005-02-23 | 北京科技大学 | 一种参与性介质遮蔽高温表面的辐射测温方法 |
CN102353478B (zh) * | 2011-10-10 | 2013-07-31 | 哈尔滨工业大学 | 半透明介质环境下非接触测温的校正方法 |
CN103389272B (zh) * | 2013-08-01 | 2015-06-17 | 哈尔滨工业大学 | 基于脉冲激光的半透明介质衰减系数和散射反照率的快速测量方法 |
CN103411905B (zh) * | 2013-08-23 | 2015-08-19 | 哈尔滨工业大学 | 一种基于短脉冲激光辐照及多信息逆问题求解算法的参与性介质辐射特性测量方法 |
CN103528978B (zh) * | 2013-11-01 | 2016-04-20 | 哈尔滨工业大学 | 利用脉冲激光加热产生的瞬态光热信号测量半透明材料热物性参数的方法 |
CN104915538B (zh) * | 2015-04-16 | 2017-09-12 | 河海大学常州校区 | 一种敞开式支架安装光伏逆变器元件温度预测方法 |
-
2015
- 2015-12-09 CN CN201510907097.9A patent/CN105319174B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN105319174A (zh) | 2016-02-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105319174B (zh) | 同时获取半透明材料温变导热系数及吸收系数的测量方法 | |
Dombrovsky et al. | Approximate analytical solution to normal emittance of semi-transparent layer of an absorbing, scattering, and refracting medium | |
CN104833695B (zh) | 基于红外热成像技术的金属薄板热导率测量方法 | |
Sfarra et al. | Ceramics and defects: Infrared thermography and numerical simulations—a wide-ranging view for quantitative analysis | |
CN105675646B (zh) | 基于本征光热信息同时测量高温半透明介质热导率及吸收系数的方法 | |
CN108333213A (zh) | 半透明多孔材料高温传导及辐射性质多参数同步测量方法 | |
CN107802243A (zh) | 基于点热源模型的生物体内部热源强度和温度分布的分析方法与装置 | |
CN105486711A (zh) | 基于红外热像仪的空间材料发射率测量系统及方法 | |
Wei et al. | Improved teaching-learning-based optimization for estimation of temperature-dependent radiative properties of semitransparent media | |
CN203798759U (zh) | 一种玻璃热物性测试仪 | |
Kessentini et al. | Three dimensional heat transfer analysis of combined conduction and radiation in honeycomb transparent insulation | |
Wang et al. | Non-uniform and volumetric effect on the hydrodynamic and thermal characteristic in a unit solar absorber | |
Niu et al. | Three-dimensional inhomogeneous temperature tomography of confined-space flame coupled with wall radiation effect by instantaneous light field | |
CN104764713A (zh) | 太赫兹瞬态热成像检测和层析成像系统及方法 | |
CN110160964A (zh) | 基于Chirp调频激光辐照的半透明材料光热特性参数检测方法 | |
CN103472036A (zh) | 基于脉冲激光辐照的半透明介质辐射特性测量方法 | |
CN106353361A (zh) | 一种测试材料涂覆涂层后激光吸收率的方法 | |
Sun et al. | Combined lock-in thermography and SQP algorithm for non-intrusive reconstruction of optical and thermal properties in semitransparent medium | |
CN103528963B (zh) | 采用多频调制激光加热与光热信息重建技术的半透明材料辐射特性测量方法 | |
CN109655483B (zh) | 一种基于深度学习算法的材料微观结构缺陷检测方法 | |
Wei et al. | Application of hybrid SPSO-SQP algorithm for simultaneous estimation of space-dependent absorption coefficient and scattering coefficient fields in participating media | |
D’Alessandro et al. | Modeling and measuring thermodynamic and transport thermophysical properties: A review | |
CN103411905B (zh) | 一种基于短脉冲激光辐照及多信息逆问题求解算法的参与性介质辐射特性测量方法 | |
Wei et al. | Inverse heat transfer analysis to determine the temperature or phase change-dependent refractive index of semitransparent materials | |
Hottel | Radiation as a diffusion process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |