粗糙表面之间接触热阻反问题研究

余远锋; 李泽卫; 郑晓亚

doi:10.6052/0459-1879-18-076

粗糙表面之间接触热阻反问题研究

doi: 10.6052/0459-1879-18-076

1 西北工业大学航天学院,西安 710072
2 中国航天科技集团公司四院四十三所,西安 710025

基金项目: 国家自然科学基金资助项目(11502207).

详细信息

作者简介:
通讯作者:郑晓亚,副教授,主要研究方向：热力耦合和界面力学. E-mail:zhengxy_8042@nwpu.edu.cn

通讯作者:
郑晓亚

中图分类号: TU311;
计量
- 文章访问数: 1131
- HTML全文浏览量: 171
- PDF下载量: 522
- 被引次数: 0
出版历程
- 收稿日期: 2018-03-19
- 刊出日期: 2018-05-18

THE INVERSE PROBLEM OF THERMAL CONTACT RESISTANCE BETWEEN ROUGH SURFACES

1 School of Astronautics, Northwestern Polytechnical University, Xi’an 710072, China
2 The 43st Institute of the Forth Academy of CASC,Xi’an 710025, China

摘要

摘要: 当两个固体表面相互接触时,由于接触面粗糙度的影响,界面间就形成了非一致接触,这种接触导致热流收缩,进而产生接触热阻. 目前的理论研究主要集中在正问题研究,对反问题的研究相对较少. 接触热阻反问题研究是通过研究部分边界温度、热流和部分测量点的温度来反演得到界面上的接触热阻. 反问题研究在很多工程领域都有应用,如航空航天、机械制造、微电子等,是工程中确定接触热阻一种快速有效的方法. 本文采用边界元法和共轭梯度法研究了二维空间随坐标变化的接触热阻反问题. 为了验证方法的准确性和可行性,假定在已知部分测量点温度和真实接触热阻的情况下,反演计算得到界面的温度和热流,进而得到接触热阻,并与真实接触热阻进行比较. 结果表明采用边界元法和共轭梯度法在无测量误差的情况下,可以准确反演获得界面的真实接触热阻. 若存在测量误差,反演计算结果对测量误差极其敏感,反演结果误差会由于测量误差的引入而被放大. 为处理这种不适定性, 采用最小二乘法对反演计算结果进行校正,结果表明采用最小二乘法能够避开反问题中一些偏离实际值较大的测量点,显著提高反演计算结果的准确性.
- 接触热阻 /
- 反问题 /
- 边界元法 /
- 共轭梯度法 /
- 最小二乘法
Abstract: When two solid surfaces are in contact, it leads to non-uniform contact because of surfaces roughness. This causes constriction of heat flux and forms thermal contact resistance. The theoretical research is mainly focused on the positive problem, but there are few studies on the inverse problem. The inverse problem of thermal contact resistance is to obtain thermal contact resistance by a part of the boundary temperature, heat flux and some of the measured point temperature. The research has been applied in many fields, such as aerospace, mechanical manufacturing, microelectronics and other fields. It is a fast and effective method to determine thermal contact resistance in engineering field. In this paper, the inverse problem of thermal contact resistance with 2-D coordinate variation was solved by the boundary element method (BEM) and the conjugate gradient method (CGM). In order to verify the accuracy and feasibility of the method, according to the measured point temperature and the assumed thermal contact resistance, the temperature and the heat flux of the interface could be obtained, and then calculated and compared with the value of actual thermal contact resistance. The results show that the actual thermal contact resistance can be accurately obtained by using the BEM and CGM without the measurement error. But there exists the measurement error, the calculated result will be extremely sensitive to the measurement error, and the error of inversion result will be amplified due to the measurement error. In order to deal with this ill-posed problem, the least-squares method (LSM) was used to correct the calculated results. The results show that it can avoid some points deviating from the actual value in the inverse problem, and obviously improve the accuracy of calculations.
- thermal contact resistance /
- inverse problem /
- boundary element method /
- conjugate gradient method /
- least square method

HTML全文

参考文献(1)

[1]	Madhusudana CV. Thermal Contact Conductance.Switzerland: Springer International Publishing, 2014: 1-3
[2]	Yovanoovich MM.Four decades of research on thermal contact gap, and joint resistance in microelectronics. IEEE Trans Compon Packagechnol, 2005, 28(2): 182-206
[3]	张平,宣益民,李强. 界面接触热阻的研究进展.化工学报, 2012, 63(2): 335-349
[3]	(Zhang Ping, Xuan Yimin, Li Qiang.Development on thermal contact resistanc. CIESC Journal, 2012, 63(2): 335-349 (in Chinese))
[4]	Levy A, Heider D, Tierney J, et al.Inter-layer thermal contact resistance evolution with the degree of intimate contact in the processing of thermoplastic composite laminates. Journal of Composite Materials, 2014, 48(4): 491-503
[5]	Wang AL, Zhao JF.Review of prediction for thermal contact resistance. Sci China Tech Sci, 2010, 53: 1798-1808
[6]	Oclon P, Lopata S.Study of the effect of fin-and-tube heat exchanger fouling on its structural performance. Heat Transfer Engineering, 2017, 1(1): 1-17
[7]	方文振,任兴杰,陶文铨.基于多块格子Boltzmann方法模拟接触热阻.工程热物理学报, 2017, 38(3): 595-599
[7]	(Fang Wenzhen, Ren Xingjie, Tao Wenquan.Numerical predictions of thermal contact resistance based on the multi-block lattice Boltzmann method. Journal of Engineering Thermophysics, 2017, 38(3): 595-599 (in Chinese))
[8]	Zhu Z, Zhang LW, Wu QK, et al.An experimental investigation of thermal contact conductance of Hastelloy C-276 based on steady-state heat flux method. Int Commun Heat Mass Transfer, 2013, 41: 63-67
[9]	刘东欢,郑小平,黄拳章等. C/C复合材料与高温合金GH600之间高温接触热阻的试验研究.航空学报, 2010, 31(11): 2189-2194
[9]	(Liu Donghuan, Zheng Xiaoping, Huang Quanzhang, et al.Experimental investigation of high-temperature thermal contact resistance between C/C composite material and superalloy GH600. Acta Aeronautica et Astronautica Sinica, 2010, 31(11): 2189-2194 (in Chinese))
[10]	Gill J, Divo E, Kassab AJ.Estimating thermal contact resistance using sensitivity analysis//Proceedings of International Mechanical Engineering Congress and Exposition. California, 2004
[11]	Zhang X, Cong PZ, Fujiwara SJ, et al.A new method for numerical simulation of thermal contact resistance in cylindrical coordinates. International Journal of Heat and Mass Transfer, 2004, 47: 1091-1098
[12]	Ngo IL, Vattikuti SVP, Byon C.Effects of thermal contact resistance on the thermal conductivity of core-shell nanoparticle polymer composites. International Journal of Heat and Mass Transfer, 2016, 102: 713-722
[13]	Lopata S, Oclon P.Numerical study of the effect of fouling on local heat transfer conditions in a high-temperature fin-and-tube heat exchanger. Energy, 2015, 92: 100-116
[14]	Cui TF, Li Q, Xuan YM, et al.Multiscale simulation of thermal contact resistance in electronic packaging. International Journal of Thermal Sciences, 2014, 83: 16-24
[15]	Taler D, Oclon P.Thermal contact resistance in plate fin-and-tube heat exchangers, determined by experimental data and CFD simulations. International Journal of Thermal Sciences, 2014, 84: 309-322
[16]	Barber AD, Werner JH, Bole BA, et al.An analysis of the effect of thermal contact resistance in a sheet-stringer structure. Journal of The Aeronautical Sciences, 1957, 3: 232-240
[17]	Hay BA, Bourouga B, Dessain C.Thermal contact resistance estimation at the blank/tool interface: Experimental approach to simulate the blank cooling during the hot stamping process. Int J Mater Form, 2010, 3: 147-163
[18]	Mutlu M, Kiliç M.Three-dimensional conjugate numerical analysis of fin and tube heat exchangers with various fin thermal conductivity values and geometric parameters. J of Thermal Science and Technology, 2016, 36(1): 85-98
[19]	de Assis LHP, Rom{ao EC. Numerical simulation of 1d heat conduction in spherical and cylindrical coordinates by fourth-order finite difference method. International Journal of Mathematics Trends and Technology, 2017, 46(3): 125-128
[20]	沈军,马骏,刘伟强. 一种接触热阻的数值计算方法. 上海航天, 2002 (4): 33-36
[20]	(Shen Jun, Ma Jun, Liu Weiqiang.A numerical calculation method of thermal contact resistance. Aerospace Shanghai, 2002, (4): 33-36 (in Chinese))
[21]	李磊涛,郑晓亚. 一种瞬态接触热导数值模拟方法.机械工程学报,2016, 52(6):174-180
[21]	(Li Leitao,Zheng Xiaoya.A numerical simulation method of transient thermal contact conductance. Journal of Mechanical Engineering, 2016, 52(6): 174-180 (in Chinese))
[22]	王安良,马松阳.一种测量板间接触热阻的新方法.工程热物理学报, 2017, 38(11): 2393-2398
[22]	(Wang Anliang, Ma Songyang.A new method for measuring the thermal contact resistance between plates. Journal of Engineering Thermophysics, 2017, 38(11): 2393-2398 (in Chinese))
[23]	黄正兴,王志刚,张鑫. 高压下碳纳米管阵列与金属接触热阻研究.大连理工大学学报,2017, 57(6): 638-643
[23]	(Huang Zhengxing, Wang Zhigang, Zhang Xin.Study of thermal contact resistance between CNT array and metal under high pressure. Journal of Dalian University of Technology, 2017, 57(6): 638-643 (in Chinese))
[24]	Dou RF, Ge TR, Liu XL, Zhi Wen.Effect of contact pressure, interface temperature, and surface roughness on thermal contact conductance between stainless steel surfaces under atmosphere condition. Int J Heat Mass Transfer, 2016, 94: 156-163
[25]	Wilson WE, Angadi SV, Jackson RL.Surface separation and contact resistance considering sinusoidal elastic-plastic multi-scale rough surface contact. Wear, 2010, 268: 190-201
[26]	Ma C, Zhao L, Shi H, et al.A geometrical mechanical thermal predictive model for thermal contact conductance in vacuum environment. J Engineering Manufacture, 2016, 230(8): 1451-1464
[27]	邓二平,赵志斌,张朋等. 温度对压接型 IGBT 器件内部接触热阻的影响. 半导体技术, 2016, 41(12): 906-923
[27]	(Deng Erping, Zhao Zhibin, Zhang Peng,et al.Influence of the temperature on the thermal contact resistance within press pack IGBTs. Semiconductor Technology, 2016, 41(12): 906-923 (in Chinese))
[28]	Fieberg C, Kneer R.Determination of thermal contact resistance from transient temperature measurements. International Journal of Heat and Mass Transfer, 2008, 51: 1017-1023
[29]	薛齐文. 多宗量热传导反问题的数值求解. [博士论文].大连:大连理工大学,2005
[29]	(Xue Qiwen.Numerical solutions of inverse heat conduction problems with multi-variables. [PhD Thesis]. Dalian: Dalian University of Thchnology, 2005 (in Chinese))
[30]	李斌. 导热几何反问题的边界元法. [博士论文].哈尔滨:哈尔滨工业大学,2008
[30]	(Li Bin.Boundary element method for inverse geometry problem of heat Conduction. [PhD Thesis]. Harbin: Harbin Institute of Technology, 2008 (in Chinese))
[31]	Alifanov OM.Solution of an inverse problem of heat conduction by iteration methods. Inzhenemo-Fizicheskii Zhurnal, 1974, 26(4): 682-689
[32]	陈生. 基于边界元方法求解热传导正反问题. [博士论文]. 长沙:湖南大学,2007
[32]	(Chen Sheng.Direct and inverse heat conduction problems solving by the boundary element method. [PhD Thesis]. Changsha: Hunan University, 2007 (in Chinese))
[33]	姚振汉, 王海涛. 边界元法. 北京: 高等教育出版社, 2010
[33]	(Yao Zhenghan,Wang Haitao. Boundary Element Method.Beijing: Higher Education Press, 2010 (in Chinese))