Boundary element method analyses of transient heat conduction in an unbounded solid layer containing inclusions

The boundary element method (BEM) is used to compute the three-dimensional transient heat conduction through an unbounded solid layer that may contain heterogeneities, when a pointwise heat source placed at some point in the media is excited. Analytical solutions for the steady-state response of this solid layer when subjected to a spatially sinusoidal harmonic heat line source are presented when the solid layer has no inclusions. These solutions are incorporated into a BEM formulation as Green’s functions to avoid the discretization of flat media interfaces. The solution is obtained in the frequency domain, and time responses are computed by applying inverse (Fast) Fourier Transforms. Complex frequencies are used to prevent the aliasing phenomena. The results provided by the proposed Green’s functions and BEM formulation are implemented and compared with those computed by a BEM code that uses the Green’s functions for an unbounded media which requires the discretization of all solid interfaces with boundary elements. The proposed BEM model is then used to evaluate the temperature field evolution through an unbounded solid layer that contains cylindrical inclusions with different thermal properties, when illuminated by a plane heat source. In this model zero initial conditions are assumed. Different simulation analyses using this model are then performed to evaluate the importance of the thermal properties of the inclusions on transient heat conduction through the solid layer.

[1]  Song-Ping Zhu,et al.  Solving linear diffusion equations with the dual reciprocity method in Laplace space , 1994 .

[2]  António Tadeu,et al.  3D sound scattering by rigid barriers in the vicinity of tall buildings , 2001 .

[3]  Günther Kuhn,et al.  Dual reciprocity BEM without matrix inversion for transient heat conduction , 2002 .

[4]  A. Tadeu,et al.  Wave motion between two fluid-filled boreholes in an elastic medium , 2002 .

[5]  W. Garvin,et al.  Exact transient solution of the buried line source problem , 1956, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[6]  Richard Paul Shaw,et al.  An integral equation approach to diffusion , 1974 .

[7]  Gary F. Dargush,et al.  Application of the boundary element method to transient heat conduction , 1991 .

[8]  Eduardo Kausel,et al.  Frequency Domain Analysis of Undamped Systems , 1992 .

[9]  António Tadeu,et al.  2.5D Green's Functions for Elastodynamic Problems in Layered Acoustic and Elastic Formations , 2001 .

[10]  M. Bouchon.  Discrete wave number representation of elastic wave fields in three-space dimensions , 1979 .

[11]  Vasco Peixoto De Freitas,et al.  Moisture migration in building walls—Analysis of the interface phenomena , 1996 .

[12]  H. Power,et al.  An efficient direct BEM numerical scheme for heat transfer problems using Fourier series , 2000 .

[13]  E. Kausel,et al.  Green's Functions for Two-and-a-Half-Dimensional Elastodynamic Problems , 2000 .

[14]  Edward L. Wilson,et al.  Numerical methods in finite element analysis , 1976 .

[15]  Alok Sutradhar,et al.  Transient heat conduction in homogeneous and non-homogeneous materials by the Laplace transform Galerkin boundary element method , 2002 .

[16]  K. Davey,et al.  SOURCE‐WEIGHTED DOMAIN INTEGRAL APPROXIMATION FOR LINEAR TRANSIENT HEAT CONDUCTION , 1996 .

[17]  Frank J. Rizzo,et al.  A method of solution for certain problems of transient heat conduction , 1970 .

[18]  Robert A. Phinney,et al.  Theoretical calculation of the spectrum of first arrivals in layered elastic mediums , 1965 .

[19]  Paulo Santos,et al.  Closed-form integration of singular terms for constant, linear and quadratic boundary elements. Part 2. SV-P wave propagation , 1999 .

[20]  Carlos Alberto Brebbia,et al.  A formulation of the boundary element method for axisymmetric transient heat conduction , 1981 .