`Advances in Mechanical EngineeringVolume 2013 (2013), Article ID 173237, 4 pageshttp://dx.doi.org/10.1155/2013/173237`
Research Article

## Analytical Method for Heat Conduction Problem with Internal Heat Source in Irregular Domains

Department of Building Environment and Equipment Engineering, Beijing University of Technology, Beijing 100024, China

Received 20 June 2013; Accepted 30 September 2013

Copyright © 2013 Hongming Fan et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

In engineering field, analytical methods can only be applied to classic heat transfer problems with regular geometric boundaries. It is difficult to apply analytical methods in solving the mathematical and physical equations in nonorthogonal boundary of irregular domains. In this paper, we presented a method by conformal mapping from the solution for heat conduction problems in regular domains to solve problems in irregular domains.

#### 1. Introduction

In engineering field, all two-dimensional heat conduction problems can be solved by the partial differential equations in the initial and boundary conditions in mathematics. For problems in regular domains, the analytical methods can be applied, including the variable separation method, the integral transformation method, Green’s function method, and the conformal mapping method. However, it is hard to apply those analytical methods to solve problems in irregular domains. Therefore, numerical methods were extensively used to solve the practical problems. In fact, analytical methods have advantages in illustrating the basics of physical concepts, giving the concise conclusions, by which it is easier to analyze a different physical phenomenon. Besides, it is less time consuming compared with the numerical method.

Jiji introduced the method of singular perturbation to Slovenia two-dimensional heat conduction in irregular domains which retains the two-dimensional nature and gives accurate estimates of the heat flux in the boundary-layer region [1]. A simple finite-difference technique using the generalized finite-difference (GFD) discretization was presented by Chung for two-dimensional heat transfer problems of irregular geometry [2]. McCorquodale et al. presented an algorithm for solving the heat equation on irregular time-dependent domains. This leaded to a method that was second-order accurate in space and time [3]. Lately researchers also developed Weighted Least-Squares Collocation Method (WLSCM) and smoothed particle hydrodynamics method for Heat Conduction Problems in Irregular Domains [4, 5]. In general, the partial differential equations for heat conduction problems can be solved by the variable separation method, in which the infinite series solution can be constructed in the orthogonal curvilinear coordinate system with the corresponding functions. And the coefficients in the above infinite series can be determined by the Fourier-Liouville series expansion. Different coordinate systems can be chosen for different problems, according to the specific boundaries. For example, using spherical coordinates or cylindrical coordinates to separate the variables in the Laplace equations, wave equations, or transport equations, an associated Legendre equation, Bessel equation or Spherical Bessel equation will appear, respectively. The image method and conformal mapping technique are well known in solving two-dimensional problems [6]. For heat conduction problems with irregular domains, we can map these irregular domains to regular domains which are called the image spaces. The solutions of the problems in the image spaces are easy to be obtained. Then by mapping the solutions back to the initial space, we can get the solutions of the initial problems [711].

#### 2. Analytical Method

##### 2.1. Dirichlet Problem in the Eccentric Domain

For steady heat conduction issues with internal heat source in irregular domains, analytical methods solve the initial state and boundary condition of the Poisson equation. The fundamental equations and boundary conditions are In the equation, is temperature distribution, is heat flux, and is thermal conductivity.

Boundary description is shown in Figure 1. In the -plane domain, and are the outer and inner radius of the eccentric ring, and the boundary conditions are known.

Figure 1: Boundary of two numeration fields.
##### 2.2. Conformal Mapping Method

The traditional solution has been to map eccentric circles in the physical space -plane with nonorthogonal boundary to concentric circles in the -plane with orthogonal boundary using conformal mapping. Solving the problem in the -plane and inverse transformation to the -plane give the final analytical solution.

Then, using linear conformal mapping, we obtained Thus, the transformation of -plane eccentric circles to -plane concentric circles can be realized. Through conformal mapping, the equations and boundary conditions of the concentric ring in the -plane are and are the outer and inner radius of the concentric ring in the -plane after transformation: Obviously, the fundamental equation is still a Poisson equation in the -plane, but the strength of the “source’’ becomes times, and is not a constant, different point by point.

Through conformal mapping, the equation is transferred into Poisson equation with a concentric circles boundary. In polar coordinates, the equation can be described as (5). In the equation, is the polar radius, and is the polar angle:

##### 2.3. Particular Solution of the Poisson Equation in the Domain of Concentric Rings

Solving the Poisson equation in the concentric rings domain (5), we assumed the form of the solution to be where is the particular solution and is the general solution of the Laplace equation: In order to get a particular solution, a coordinate transformation is introduced: Then, the Poisson equation can be rewritten as Thus, we can construct Then Taking it into the equation above Undetermined coefficients and can be calculated Thus, the particular solution is According to the form of the particular solution , by trigonometric function can be expressed as Thus, the particular solution is

##### 2.4. General Solution of the Poisson Equation in the Domain of Concentric Rings

On the inner boundary and the outer boundary of the concentric rings, we get The theoretical solution of the Dirichlet problem in the Laplace equation of the concentric rings is According to the form of the boundary conditions, we can assume that By substitution of the boundary conditions on the inner boundary into (19) and the outer boundary into (15), then , , , , , , , and can be determined by comparing the factors of the two sides in (15) and (19), and we can obtain the analytical solution (6) of the Poisson equation (5):

##### 2.5. Calculations with Typical Positions of the Eccentric Rings with Different Radii

In this paper, for the Dirichlet problem, °C and °C. Further, the outer radius was set at  m, the inner radius  m, distance of the two centers , heat flux  kW/m3, and thermal conductivity of the pipe wall  W/(m·°C), and we determined the distribution of temperature in the eccentric circles with different radii.

We selected points of different radii on the typical position in the -plane and then conducted coordinate transformation and solved the Poisson equation of the concentric ring in the -plane, thus obtaining the solution to the original problem with nonorthogonal boundary. Table 1 lists the calculations of the typical position of three eccentric circles with different radii.

Table 1: Calculation and simulation results of three typical positions.

As seen from Table 1, the temperature distribution trend of the typical positions in the -plane was substantially the same; the temperature gradually reduced as the radius increased. When , the temperature along the x-axis direction gradually reduced as the radius increased. When , the temperature along the -axis gradually reduced as the radius increased. When , the temperature along the negative x-axis also had a downward trend as the radius increased. Comparing with , the temperature downtrend of the former was slower; the temperature gradient along the negative x-axis was larger than the one along the x-axis. Compared with the temperature distribution simulated by ANSYS, the analytical solution derived in the article is correct.

#### 3. Conclusions

Steady heat conduction issues with internal heat sources in irregular domains can be resolved by solving the Poisson equation with the initial and boundary conditions. In this study, the conformal mapping method is presented to solve this equation. Firstly, the eccentric rings in the physical space -plane with nonorthogonal boundary was mapped to the concentric rings with orthogonal boundary. And then the solution in the -plane is inverse transformed to the -plane. Finally we obtained the analytical solution to the initial problems. Also, the temperature field was analyzed to obtain the temperature distribution of typical eccentric ring positions. And the result is verified by the simulation.

#### References

1. L. M. Jiji, “Singular-perturbation solutions of conduction in irregular domains,” Quarterly Journal of Mechanics and Applied Mathematics, vol. 27, no. 1, pp. 45–55, 1974.
2. K. C. Chung, “A generalized finite-difference method for heat transfer problems of irregular geometries,” Numerical Heat Transfer, vol. 4, no. 3, pp. 345–357, 1981.
3. P. McCorquodale, P. Colella, and H. Johansen, “A Cartesian grid embedded boundary method for the heat equation on irregular domains,” Journal of Computational Physics, vol. 173, no. 2, pp. 620–635, 2001.
4. Y. J. Dai, X. H. Wu, and W. Q. Tao, “Weighted least-squares collocation method (WLSCM) for 2-D and 3-D heat conduction problems in irregular domains,” Numerical Heat Transfer B, vol. 59, no. 6, pp. 473–494, 2011.
5. V. Vishwakarma, A. K. Das, and P. K. Das, “Steady state conduction through 2D irregular bodies by smoothed particle hydrodynamics,” International Journal of Heat and Mass Transfer, vol. 54, no. 1–3, pp. 314–325, 2011.
6. R.-L. Lin, “Explicit full field analytic solutions for two-dimensional heat conduction problems with finite dimensions,” International Journal of Heat and Mass Transfer, vol. 53, no. 9-10, pp. 1882–1892, 2010.
7. C. S. Wu, Mathematical Methods for Physics, Press of Peking University, 2001.
8. D. R. Guo, Mathematical Methods for Physics, People's Education Press, Beijing, China, 1965.
9. G. Arfken, Mathematical Methods for Physics, Academic Press, New York, NY, USA, 4th edition, 1997.
10. H. M. Fan, Y. Jiang, and Z. Y. He, “Thermal analysis of a straight buried pipe,” Journal of Tsinghua University, vol. 42, no. 6, pp. 806–809, 2002 (Chinese).
11. Z. Y. He and H. M. Fan, “SH Wave scattering around the interface hole,” Mechanical Press, vol. 34, no. 1, pp. 68–76, 2002.