Steady-state nonlinear heat conduction in composite materials using the method of fundamental solutions

https://doi.org/10.1016/j.cma.2008.02.011Get rights and content

Abstract

The steady-state heat conduction in composite (layered) heat conductors with temperature dependent thermal conductivity and mixed boundary conditions involving convection and radiation is investigated using the method of fundamental solutions with domain decomposition. The locations of the singularities outside the solution domain are optimally determined using a non-linear least-squares procedure. Numerical results for non-linear bimaterials are presented and discussed.

Introduction

In many heat transfer problems the assumption of constant thermal conductivity, i.e. that the heat conductors are homogeneous within the whole temperature variation interval, may lead to unacceptable errors in high-temperature environments or if large temperature differences are present, see [31]. In the steady-state situation, the nonlinearity associated with the temperature dependence of the thermal conductivity can be removed by employing the Kirchhoff transformation, which replaces the original nonlinear partial differential equation in divergence form by the Laplace equation in the transformed space, see [11]. Boundary conditions of the Dirichlet (first kind) or Neumann (second kind) types pose no problem for the transformation, but the Robin convective (third kind) boundary conditions become non-linear. Although this non-linearity is not strong, convergence problems may arise if radiative heat transfer (fourth kind) boundary conditions are also present, see [12]. Since all the non-linearities are transferred to the boundary conditions, the Kirchhoff transformation approach is very well-suited for applying the boundary element method (BEM) [7], [20], or, more recently, the method of fundamental solutions (MFS) [22], [23]. In the same manner these techniques can be extended to composite bodies through the subregion technique. In it, each region is dealt with separately and then the whole body is linked together by applying compatibility and equilibrium conditions along the interfaces between subregions.

Two-dimensional boundary value problems of heat conduction in nonlinear composite materials have been the subject of several studies using the BEM [2], [6], [8]. However, the implementation of the BEM becomes quite tedious, especially in three-dimensional irregular domains. Moreover, the evaluation of the gradient of the temperature solution on the boundary requires the use of finite differences or the evaluation of hypersingular integrals. In order to alleviate some of these difficulties, we propose the use of the MFS. The merits and drawbacks of the MFS over the BEM for solving elliptic boundary value problems are thoroughly discussed in Refs. [10], [13], [14], [19], [25], [32]. Recently, the MFS has been made applicable to inhomogeneous elliptic equations [1] and inverse problems [33].

Prior to this study, the MFS was used for the solution of problems of heat conduction in linear layered materials with linear boundary conditions [3]. It is the purpose of this paper to extend this analysis to nonlinear materials with nonlinear boundary conditions.

The mathematical formulation of the problem is given in Section 2. The MFS and its implementation are described in Section 3. Numerical results are presented in Section 4 and in Section 5 some conclusions and ideas for future applications are given.

Section snippets

Mathematical formulation

Consider a bounded domain ΩRd, d2, with piecewise smooth boundary Ω, formed from two (or more) subregions Ω1 and Ω2 separated by the interfacial surface Γ12=Ω1Ω2. The material of subregion Ω1 has a temperature dependent thermal conductivity k1>0 and the material of subregion Ω2 has a different thermal conductivity k2>0. The governing steady-state heat conduction equations are·(ki(Ti)Ti)=0,inΩi,i=1,2,where Ti is the temperature solution in domain Ωi, i=1,2, and, for the sake of

The method of fundamental solutions (MFS)

As the sources of nonlinearity are associated with the boundary conditions (2.14), (2.15), (2.17), (2.18) only, the boundary value problem (2.11), (2.12), (2.13), (2.14), (2.15), (2.16), (2.17), (2.18) for each subregion can, following [23], be converted into a minimization problem, or equivalently an algebraic system of nonlinear equations, using the MFS.

From [9], [26], the MFS approximations for the solutions Ψ1 and Ψ2 of the Laplace Eq. (2.11) have the formΨNi(ci,ξi;x)=k=1NckiGd(ξki,x),xΩ¯i

Numerical results and discussion

In this section, we present numerical results obtained from the application of the MFS described in the previous section.

Conclusions

In this paper, the application of the MFS to steady-state non-linear heat conduction problems in composite heat conductors has been investigated. The MFS is used in conjunction with a domain decomposition technique and the method recasts the problem as a non-linear minimization problem. The numerical results obtained are in good agreement with the available analytical solutions showing high accuracy and stable convergence, and probably with the BEM results of [8] if these were to be corrected.

Acknowledgement

The authors would like to thank the UK Royal Society for supporting this research.

References (33)

  • J.R. Berger et al.

    The method of fundamental solutions for heat conduction in layered materials

    Int. J. Numer. Meth. Engrg.

    (1999)
  • J.R. Berger et al.

    Stress intensity factor computation using the method of fundamental solutions: mixed mode problems

    Int. J. Numer. Meth. Engrg.

    (2007)
  • R. Bialecki et al.

    Solving nonlinear steady-state potential problems in inhomogeneous bodies using the boundary-element method

    Numer. Heat Transfer, Part B

    (1989)
  • R. Bialecki et al.

    Boundary element solution of heat conduction problems in multizone bodies of non-linear material

    Int. J. Numer. Meth. Engrg.

    (1993)
  • A. Bogomolny

    Fundamental solutions method for elliptic boundary value problems

    SIAM J. Numer. Anal.

    (1985)
  • H.S. Carslaw et al.

    Conduction of Heat in Solids

    (1959)
  • Cited by (0)

    View full text