Mathematical Problems in Engineering

Volume 2010, Article ID 784794, 13 pages

http://dx.doi.org/10.1155/2010/784794

## Calculation of Some Integrals Arising in Heat Transfer in Geothermics

^{1}Departamento de Ciencias Experimentales y Matemáticas, Universidad Católica de Valencia, 46002 Valencia, Spain^{2}Departamento de Matemáticas, Universidad de Pinar del Río, 20200 Pinar del Río, Cuba^{3}Instituto Universitario de Matemática Pura y Aplicada, Universidad Politécnica de Valencia, 46022 Valencia, Spain

Received 28 September 2010; Accepted 29 November 2010

Academic Editor: Miguel Cerrolaza

Copyright © 2010 J. L. G. Santander 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

We calculate some integrals involved in the temperature field evaluation of the ground, heated by a borehole heat exchanger. This calculation allows a faster computation of that component of the temperature field which involves the periodic oscillation of the ambient temperature or the ambient heat flux.

#### 1. Introduction

Ground-coupled heat pumps (GCHPs) are an attractive choice of system for heating and cooling buildings [1]. By comparison with standard technologies, these pumps offer competitive levels of comfort, reduced noise levels, lower greenhouse gas emissions, and reasonable environmental safety. Furthermore, their electrical consumption and maintenance requirements are lower than those required by conventional systems and, consequently, have lower annual operating cost [2].

A GCHP exchanges heat with the ground through a buried U tube loop. Since this exchange strongly depends on the thermal properties of the ground, it is very important to have a knowledge of these properties when designing GCHP air-conditioning systems. Methods to estimate them include literature searches, conducting laboratory experiments on soil/rock samples, and/or performing field tests. The *in situ* tests are based on studying the thermal response of the borehole heat exchanger to a constant heat injection or extraction. The outputs of the thermal response test are the inlet and outlet temperature of the heat-carrier fluid as a function of time. From these experimental data, and with an appropriate model describing the heat transfer between the fluid and the ground, the thermal conductivity of the surroundings is inferred.

The results presented in this study are based on a model of the underground heat transfer due to a borehole heat exchanger assuming that it works as a line heat source of finite depth inside the ground. We assume as well a constant power of length unit for the heat source considered and a periodic oscillation of the ambient temperature or heat flux on the ground surface.

This paper is organized as follows. Section 2 describes the equations that model the heat transfer inside the ground due to a borehole heat exchanger. Section 3 solves the integrals given in literature for periodic oscillations of the ambient temperature or heat flux, on the ground surface. Section 4 shows that our results agree with the literature in some particular cases, such as the solution on the ground surface and the solution of the quasistationary regime.

#### 2. The Equations

Let us consider the heat equation in cylindrical coordinates, in which is the temperature field, On the one hand, the thermal parameters of (2.1) are the ground thermal diffusivity , the heat flow per length unit , and the volumetric heat capacity of ground . On the other hand, the geometry of the heat source, depicted in Figure 1, has been taken into account introducing the Heaviside function for its finite length and the Dirac delta distribution for its infinitely thin width. Equation (2.1) is subject to to the initial condition where is the undisturbed ground temperature and is the geothermal gradient. We may consider that (2.1) is subject to a Dirichlet boundary condition, which represents a time-dependent temperature on the surface , or to a Neumann boundary condition, which represents a time-dependent heat flux on the surface , where is the ground thermal conductivity. We split the Cauchy problem posed in (2.1), (2.2), and (2.3) or (2.4) in three terms [3], so that and where the functions , and satisfy the Cauchy problems presented below.

##### 2.1. The Problem for

The partial differential equation that satisfies is subject to a homogeneous initial condition, and to a homogeneous Dirichlet boundary condition, or to a Neumann boundary condition, homogenous as well, The solution to (2.6), (2.7), and (2.8) is [4] and the solution to (2.6), (2.7), and (2.9) is that we will present in a future work.

##### 2.2. The Problem for

The partial differential equation that satisfies is subject to the initial condition and to a Dirichlet boundary condition, or to a Neumann homogeneous boundary condition The solution to (2.12), (2.13), and (2.14) is [5, Equation 2.4.13], and the solution to (2.12), (2.13), and (2.15) is [5, Equation 2.9.7] where the error function is defined as

##### 2.3. The Problem for

The partial differential equation that satisfies is subject to an homogeneous initial condition and to the time-dependent Dirichlet boundary condition or to the time-dependent Neumann boundary condition The solution to (2.19), (2.20) and (2.21) is [5, Equation 2.5.1.] and the solution to (2.19), (2.20) and (2.22) is [5, Equation 2.9.9]

#### 3. Harmonic Analysis of the Boundary Condition

We may consider that on the ground surface we have a periodic temperature or heat flux oscillation, due to daily or annual cycles. Therefore, is a periodic function of a certain period , continuous and bounded, that we may expand in Fourier series where the frequency is , and the coefficients and are given by where .

##### 3.1. Dirichlet Boundary Condition

Substituting (3.1) in (2.23), we obtain where is the complementary error function. The integrals given in (3.3) for may be rewritten as where we have defined In order to calculate (3.5), let us apply the integral [6, Equation 7.4.33] taking , , and , so that Finally, substituting (3.7) in (3.4), we obtain for , where we have defined As far as we know, the integrals given in (3.8) are not reported in the literature.

###### 3.1.1. Neumann Boundary Condition

Substituting (3.1) in (2.24), we obtain Integrating by parts and performing the substitution , considering , we may calculate the integral given in (3.10) In order to calculate the integrals given in (3.11) and (3.12), let us perform the substitution , considering , so that and similarly where we have defined Taking , , and in (3.6), we may rewrite (3.16) as Finally, substituting (3.17) in (3.14) and (3.15), for , we obtain As far as we know, the integrals given in (3.18) are not reported in the literature.

#### 4. Particular Cases

As a consistency check, let us verify that in certain particular cases the new integrals given in (3.8) and (3.18) are reduced to integrals reported in the literature. These particular cases have to do with the quasistationary regime and the solution on the ground surface.

##### 4.1. Quasistationary Regime

Since there is a periodic oscillation on the ground surface, the surroundings of the borehole heat exchanger never reach a stationary regime. However, according to [5, Section 2.6], we may define a quasistationary regime in which the ground temperature field is stabilized periodically. In this quasistationary regime, (3.8) becomes where we have defined Similarly, (3.18) becomes where we have defined

###### 4.1.1. Dirichlet Case

Taking into account that and , we may calculate the limit given in (4.3), so that (4.1) becomes where we have substitute . Thus, for , Similarly, we may obtain Notice that (4.7) and (4.8) are sinusoidal in , with an amplitude exponentially decreasing in depth . In the literature, we may find [7, Equation 3.928.1-2]. where and thus, taking , , and , Therefore, rewriting the left side of (4.7) as and applying (4.11), taking , , we eventually get the same result as (4.7). Similarly, we may check the result given in (4.8).

###### 4.1.2. Neumann Case

In a similar way as in the previous subsection, we may calculate (4.4) for , arriving at Once again, (4.13) and (4.7) are sinusoidal in , with an amplitude exponentially decreasing in depth . In the literature, we may find [7, Equation 3.957.1-2] where , and . Therefore, taking in (4.15),, , , , and knowing that the Macdonald function of order is [8, Equation 5.5.5], we eventually get the same results as (4.13) and (4.14).

##### 4.2. Solution on the Ground Surface

Notice that in the special case , (3.8) becomes trivial,
while (3.18) are reduced to
Since [8, Equation 2.4.1]
where
are the *Fresnel integrals*; then (4.18) may be expressed as
The results given in (4.21), agrees with [7, Equation 2.653.1-2]

#### Acknowledgments

The authors wish to thank the financial support received from Generalitat Valenciana under Grant no. GVA 3012/2009 and from Universidad Politécnica de Valencia under Grant no. PAID-06-09. This work has been partially supported by the Structural Funds of the European Regional Development Fund (ERDF).

#### References

- J. F. Urchueguía, M. Zacarés, J. M. Corberán, A. Montero, J. Martos, and H. Witte, “Comparison between the energy performance of a ground coupled water to water heat pump system and an air to water heat pump system for heating and cooling in typical conditions of the European Mediterranean coast,”
*Energy Conversion and Management*, vol. 49, no. 10, pp. 2917–2923, 2008. View at Publisher · View at Google Scholar - J. W. Lund, “Ground source (geothermal) heat pumps,” in
*Course on Heating with Geothermal Energy: Conventional and New Schemes. World Geothermal Congress*, P. J. Lineau, Ed., pp. 1–21, Kazuno, Japan, 2000, Short Courses. View at Publisher · View at Google Scholar - T. V. Bandos, A. Montero, and A. Montero, “Finite line-source model for borehole heat exchangers: effect of vertical temperature variations,”
*Geothermics*, vol. 38, no. 2, pp. 263–270, 2009. View at Publisher · View at Google Scholar - P. Eskilson,
*Thermal analysis of heat extraction boreholes*, Ph.D. thesis, Department of Mathematical Physics, University of Lund, Lund, Sweden, 1987. - H. S. Carslaw and J. C. Jaeger,
*Conduction of Heat in Solids*, Oxford Science Publications, The Clarendon Press, Oxford University Press, New York, NY, USA, 2nd edition, 1988. - M. Abramowitz and I. A. Stegun,
*Handbook of Mathematical Functions*, Dover Publications Paperback, New York, NY, USA, 1972. - I. S. Gradshteyn and I. M. Ryzhik,
*Table of Integrals, Series, and Products*, Elsevier/Academic Press, Amsterdam, Netherlands, 7th edition, 2007. - N. N. Lebedev,
*Special Functions and Their Applications*, Dover, New York, NY, USA, 1972.