瞬态热传导问题的精细积分边界元法  

Precise Integration Boundary Element Method for Transient Heat Conduction Problems

作　　者：周枫林 袁小涵 钦宇 潘先云 余江鸿[1] ZHOU Feng-lin;YUAN Xiao-han;QIN Yu;PAN Xian-yun;YU Jiang-hong(College of Mechanical Engineering,Hunan University of Technology,Zhuzhou 412007,China)

机构地区：[1]湖南工业大学机械工程学院,株洲412007

出　　处：《科学技术与工程》2022年第20期8588-8596,共9页Science Technology and Engineering

基　　金：国家自然科学基金(11602082);湖南省自然科学基金(2021JJ30211);湖南省教育厅优秀青年项目(19B145)。

摘　　要：热传导为最普遍的工程问题之一,为解决含内部热源的三维瞬态热传导问题,通过双重互易边界元法(dual reciprocal boundary element method, DRBEM)与精细积分法(precise integration method, PIM)耦合进行求解。该方法根据含有内部热源的各向同性介质瞬态常系数热传导问题的控制方程,通过加权余量法推导出相应的边界积分方程,然后用双互易法(dual reciprocity method, DRM)处理得到的边界积分方程,将热源项和温度关于时间导数项引起的域积分通过径向基函数(radial basis function, RBF)逼近后转化为边界积分。之后将边界积分方程离散,得到与时间相关的一阶常系数微分方程组,最后,在获得解析解的过程中,通过PIM处理其中的矩阵指数函数(matrix exponential function, MEF)。最后通过3个数值算例进行验证。结果表明:在3种边界条件下,数值求解结果都具有很好的准确性。可见该方法的可行性和稳定性。Heat conduction is one of the most common engineering problems.In order to solve the three-dimensional transient heat conduction problem with internal heat source,it was solved by coupling the dual reciprocal boundary element method(DRBEM)and the precise integration method(PIM).According to the governing equation of transient constant coefficient heat conduction problem in isotropic medium with internal heat source,the corresponding boundary integral equation was derived by the weighted residual method.The dual reciprocity method(DRM)was applied to deal with the obtained boundary integral equation.The domain integral caused by the heat source term and the time derivative term of temperature was transformed into the boundary integral after being approximated by radial basis function(RBF).Then,the boundary integral equation was discretized to obtain the first-order constant coefficient differential equations related to time.Finally,in the process of obtaining the analytical solution,the matrix exponential function(MEF)was processed by PIM.Finally,three numerical examples were used to verify the results.The results show that the numerical results have good accuracy under the three boundary conditions.It can be seen that the method is feasible and stable.

关 键 词：边界元法 双互易方法 精细积分方法 瞬态热传导 

分 类 号：O241.82[理学—计算数学]