Three-dimensional transient heat conduction analysis by Laplace transformation and multiple reciprocity boundary face method

In this paper, a new multiple reciprocity formulation is developed to solve the transient heat conduction problem. The time dependence of the problem is removed temporarily from the equations by the Laplace transform. The new formulation is derived from the modified Helmholtz equation in Laplace space (LS), in which the higher order fundamental solutions of this equation are firstly derived and used in multiple reciprocity method (MRM). Using the new formulation, the domain integrals can be converted into boundary integrals and several non-integral terms. Thus the main advantage of the boundary integral equations (BIE) method, avoiding the domain discretization, is fully preserved. The convergence speed of these higher order fundamental solutions is high, thus the infinite series of boundary integrals can be truncated by a small number of terms. To get accurate results in the real space with better efficiency, the Gaver-Wynn-Rho method is employed. And to integrate the geometrical modeling and the thermal analysis into a uniform platform, our method is implemented based on the framework of the boundary face method (BFM). Numerical examples show that our method is very efficient for transient heat conduction computation. The obtained results are accurate at both internal and boundary points. Our method outperforms most existing methods, especially concerning the results at early time steps.