数值计算与计算机应用
       首页 |  期刊介绍 |  编委会 |  投稿指南 |  期刊订阅 |  下载中心 |  联系我们 |  重点论文 |  在线办公 | 
数值计算与计算机应用  2019, Vol. 40 Issue (3): 215-229    DOI:
论文 最新目录 | 下期目录 | 过刊浏览 | 高级检索 Previous Articles  |  Next Articles  
有限差分法求解声子热输运方程
陈建灵1,2, 冯仰德1
1 中国科学院计算机网络信息中心, 高性能计算技术与应用发展部, 北京 100190;
2 中国科学院大学, 北京 100049
FINITE DIFFERENCE METHOD FOR SIMULATING PHONON HEAT TRANSPORT PROCESS
Chen Jianling1,2, Feng Yangde1
1 Computer Network Information Center, Department of High Performance Computing Technology and Application Development, CAS, Beijing 100190, China;
2 University of Chinese Academy of Sciences, Beijing 100049, China
 全文: PDF (641 KB)   HTML (1 KB)   输出: BibTeX | EndNote (RIS)      背景资料
摘要 声子辐射输运方程是一个微分-积分方程,使用有限差分法数值求解具有边界条件的声子辐射输运方程,可借助Gauss-Seidel方法和有限差分离散化来保证得到稳定的数值解.通过对室温下Ge/Si/Ge薄膜的声子热输运过程进行一维数值模拟,以及对界面处使用漫射失配界面模型,可以得到温度沿着薄膜法向方向的变化曲线和Ge/Si薄膜厚度比例与温度在界面处跳跃值对材料整体结构的热导率的影响,以及随着薄膜厚度的增加薄膜热导率的变化趋势;对硅薄膜声子热输运过程进行二维数值模拟,可以得到温度沿着薄膜法向和面向方向的分布情况,以及当薄膜宽度和厚度比例不同时温度的变化和法向热导率的变化.
服务
把本文推荐给朋友
加入我的书架
加入引用管理器
E-mail Alert
RSS
作者相关文章
关键词声子辐射输运方程   有限差分法   离散化   温度   热导率     
Abstract: The Equation of Phonon Radiative Transfer is a differential-integral equation. The finite difference method is used to numerically solve the Equation of Phonon Radiative Transfer with boundary conditions. The Gauss-Seidel method and the finite difference discretization can be used to ensure stable numerical solutions. Through one-dimensional numerical simulation of the phonon heat transport process of Ge/Si/Ge film at room temperature and by using the diffuse mismatch interface model at the interface,it can evaluate the influence of the temperature along the normal direction of the film and the thickness ratio of the Ge/Si film and the temperature jump at the interface on the thermal conductivity of the overall structure of the material, and the change of the thermal conductivity of the film as the thickness of the film increases. The two-dimensional numerical simulation of the phonon heat transport process of the silicon film can obtain the temperature along the the normal and orientation-oriented distribution of the film, and the change in temperature and the thermal conductivity when the film width and thickness ratio are different.
Key wordsthe Equation of Phonon Radiative Transfer   the finite difference method   discretization   temperature   the thermal conductivity   
收稿日期: 2018-12-11; 出版日期: 2019-09-12
基金资助:

中国科学院科研信息化应用工程(XXH13506-407).

引用本文:   
. 有限差分法求解声子热输运方程[J]. 数值计算与计算机应用, 2019, 40(3): 215-229.
. FINITE DIFFERENCE METHOD FOR SIMULATING PHONON HEAT TRANSPORT PROCESS[J]. Journal on Numerical Methods and Computer Applicat, 2019, 40(3): 215-229.
 
[1] Majumdar A. Microscale Heat Conduction in Dielectric Thin Films[J]. Journal of Heat Transfer, 1993, 115(1):7-16.
[2] 韩亚芬. 纳米复合隔热材料导热机理与特性研究[D]. 哈尔滨工业大学, 2013:13-37.
[3] Venkatasubramanian R. Chapter 4——Phonon Blocking Electron Transmitting Superlattice Structures as Advanced Thin Film Thermoelectric Materials[J]. Semiconductors & Semimetals, 2001, 71:175-201.
[4] Chen G. Nanoscale energy transport and conversion[M]. Oxford University Press, 2005:348-366.
[5] Miyazaki K, Arashi T, Makino D, et al. Heat conduction in micro-structured materials[J]. IEEE Transactions on Components & Packaging Technologies, 2006, 29(2):247-253.
[6] Yilbas B S, Mansoor S B. Phonon Transport in Two-Dimensional Silicon-Diamond Film Pair[J]. Journal of Thermophysics & Heat Transfer, 2013, 27(3):465-473.
[7] Haider Ali. Heat transport characteristics in thin films[D]. King Fand University of Petroleum & Minerals, 2014:74-199.
[8] Chu A J H. Phonon transport in two-dimensional silicon thin film:influence of film width and boundary conditions on temperature distribution[J]. European Physical Journal B, 2012, 85(7):1-9.
[9] R.Yang, et al. Simulation of Nanoscale Multidimensional Transient Heat Conduction Problems Using Ballistic-Diffusive Equations and Phonon Boltzmann Equation[J]. Journal of Heat Transfer, 2005, 127(3):298-306.
[10] 裴志坚. 基于高斯-赛德尔迭代法及MATLAB软件的电路方程组求解方法[J]. 北京工业职业技术学院学报, 2017, 16(3):22-25.
[11] 吴专保. 线性方程组的几种数值方法的MATLAB程序[J]. 荆楚理工学院学报, 2007, 22(3):78-81.
[12] 杨洋. 界面热阻的晶格动力学模拟[D]. 东南大学, 2012:6-36.
[13] 张兴丽, 孙兆伟. 基于分子动力学对超晶格结构界面热阻的模拟研究[J]. 航空材料学报, 2011, 31(4):7-10.
[14] 赵兰萍, 徐烈, 李兆慈, 等. 固体界面间接触导热的机理和应用研究[J]. 低温工程, 2000, (4):29-34.
[15] M.Mohamadi, et al. Numerical Study of Thermal Conductivity Reduction in Nanolayered Si-Ge Structures by Solving the Equation of Phonon Radiative Transfer[J]. Journal of the Serbian Society for Computational Mechanics, 2017, 11(1):29-45.
[16] 冉鑫, 郭洋裕, 王沫然. 双层薄膜界面声子输运的Monte Carlo模拟[J]. 工程热物理学报, 2017, V38(4):841-846.
[17] Nield D A, Bejan A. Heat Transfer Through a Porous Medium[M]. Convection in Porous Media. Springer New York, 2017, 31-46.
[1] 邱俊, 胡晓, 王汉权. 数字图像修复的变分方法与实现过程[J]. 数值计算与计算机应用, 2016, 37(4): 273-286.
[2] 罗振东. 非饱和土壤水流问题的CN有限元格式[J]. 数值计算与计算机应用, 2014, 36(4): 355-362.
[3] 刘艳峰, 魏兵. 基于二维FDFD分析金属柱阵列的双站RCS[J]. 数值计算与计算机应用, 2014, 35(2): 81-91.
[4] 腾飞,罗振东. 非饱和土壤水流方程的CN广义差分法[J]. 数值计算与计算机应用, 2014, 36(2): 205-214.
[5] 吴正, 阳佳慧, 张依文. 具有巴黎期权特性的可转债定价问题研究[J]. 数值计算与计算机应用, 2013, 34(4): 295-304.
[6] 闵涛, 任菊成, 耿蓓. Chebyshev谱-Euler混合方法求解一类非线性Burgers方程[J]. 数值计算与计算机应用, 2013, 34(2): 81-88.
[7] 黄文艳, 魏剑英, 葛永斌. 三维非定常不可压涡量——速度Navier-Stokes方程组的有限差分法[J]. 数值计算与计算机应用, 2012, 33(4): 301-311.
[8] 方乐, 洪洁瑛. 弱非均匀 Voronoi 图细胞面积/体积的快速近似算法[J]. 数值计算与计算机应用, 2011, 32(2): 135-142.
[9] 周艳杰, 孙萍, 罗振东, 杨晓忠. 空气污染方程的全离散化混合元格式[J]. 数值计算与计算机应用, 2008, 30(4): 425-436.
[10] 吴泽艳,薛晓春. FDTD计算中一种UPML吸收边界与其内部计算区域的统一建模方法[J]. 数值计算与计算机应用, 2008, 29(1): 8-14.
[11] 王珏,张法勇,. 带有多项式非线性项的高维反应扩散方程有限差分格式的长时间行为[J]. 数值计算与计算机应用, 2007, 29(2): 177-188.
[12] 杜其奎,余德浩. 抛物型初边值问题的有限元与边界积分耦合的离散化及其误差分析[J]. 数值计算与计算机应用, 1999, 21(2): 199-208.
[13] 罗振东,王烈衡. 非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅱ):向后一步的Euler全离散化格式[J]. 数值计算与计算机应用, 1998, 20(4): 431-448.
[14] 罗振东,王烈衡. 非定常的热传导-对流问题的非线性Galerkin混合元法(Ⅰ):时间连续情形[J]. 数值计算与计算机应用, 1998, 20(3): 305-324.
Copyright © 2008 数值计算与计算机应用 版权所有
中国科学院数学与系统科学研究院 《数值计算与计算机应用》编辑部
北京2719信箱 (100190) Email: szjs@lsec.cc.ac.cn
Support by Beijing Magtech Co.ltd   E-mail:support@magtech.com.cn
京ICP备05002806号-10