Temperature distributions in a periodically stratified layer with a vertical dug well

The stationary problems of heat conduction in a periodically stratified layer with a vertical cylindrical hole is considered. The lateral surface of the hole is kept at a given piece-wise constant temperature. The boundary planes at zero temperature are taken into account, and the ideal thermal contact between the composite constituents is assumed. The problem within the homogenized model with microlocal parameters is solved.


Introduction
The knowledge of temperature distributions in layered bodies with periodic structures play an important role in many technological applications. Periodically stratified materials can be found in nature (varved clays, sandstonesandstone-slates, sandstone-shales, thin-layered limenstones) or are made by man. The thermal behavior of systems, which are made of a large number of repeated laminae, is described within the framework of the classical heat conduction approach by a partial differential equation with discontinuous, oscillating coefficients. Thus formulated problems are too complicated for an analytical and numerical analysis. For this reason, applications of some approximate models seem to be useful. Some effective approaches to the modeling of heat conduction problem in periodic composites are presented in papers of Aurialt [3], Bufler [4], Bufler and Meier [5]. One of the models is the homogenized model with microlocal parameters devised by Woźniak [2] and developed for microperiodic layered composites by Matysiak and Woźniak [1]. It is important that this model satisfies continuity conditions for temperature and heat fluxes on interfaces (the conditions of perfect thermal contact). The homogenized model with microlocal parameters was applied to solve many problems, which were partially reassumed in papers of Matysiak [6], Woźniak and Woźniak [7]. The problems of heat conduction in a periodically stratified layer with a vertical, cylindrical hole was considered by Matysiak and Perkowski [8]. In this paper, zero temperature or zero radial heat flux on the lateral surface of hole is assumed and the body is heated by given temperature on the upper boundary plane. The problem has been solved within the framework of the homogenized model with microlocal parameters.
In the present paper the stationary problem of heat conduction in a periodically stratified layer with vertical, cylindrical hole is considered. The upper and lower boundary planes with circular cut-out are assumed to be kept at zero temperature. The later surface of hole is heated by a given constant temperature on some part of the surface (lower), the remaining part is kept at zero temperature. Moreover, the ideal thermal contact between the laminae is assumed. The considered problem can be determined some approximated model of heat conduction in a periodically layered stratum with a vertical cylindrical dug well, partically filled by fluid (for instance water). This a reason to entitle of the presented contribution.
2 Formulation and solution of the problem Consider a rigid, periodically stratified layer (stratum) with a cylindrical hole (dug well) normal to the layering. The constituents of the nonhomogeneous body are assumed to be isotropic and homogeneous heat conductors, and denoted by 1 and 2; here and in the sequel, all quantities (e.g. thermal constants, heat fluxes) pertaining to these sublayers will be associated with the index a taking the values 1 and 2, respectively. The problem is investigated by using the cylindrical coordinate system (r, u, z) with the axis z being the symmetry axis of the hole, see Fig. 1. Let a be the radius of hole and l 1 , l 2 be the thicknesses of the subsequent layers, l = l 1 ? l 2 be thickness of the fundamental lamina (the repeated unit). Let K 1 , K 2 denote the coefficients of heat conductivities of the subsequent layers. Let the parts of planes z = h 1 , r C a, and z = -h 2 , r C a, where h 1 , h 2 are positive constants, be the boundaries of the nonhomogeneous body, and h = h 1 ? h 2 be thickness of the layer, see Fig. 1. Moreover it is assumed that h 1 = n 1 l, h 2 = n 2 l, n 1 , n 2 [ N and n = n 1 ? n 2 is a sufficiently large natural number. The upper and lower boundary planes with circular cutouts are assumed to be kept at zero temperature (zero temperature is taken as a surroundings temperature). Let the upper part of lateral surface of the hole r = a, -h 2 \ z \ 0, (see, Fig. 1) be kept at zero temperature, the reaming part of the lateral surface 0 \ z \ h 1 , r = a is heated by constant temperature 0 0 . The presented above problem is stationary and axially symmetric. The ideal contact between the layers being constituents of the periodically stratified stratum is taken into account. The continuity conditions on the interfaces lead to some difficulty in analytical and numerical approaches. For small thickness of the fundamental lamina l with a comparison to the thickness of the nonhomogeneous layer h, l ( h it is suitable to use the homogenized model with microlocal parameters (see Woźniak [2], Matysiak and Perkowski [8]) to the approximated formulation of the considered problem. This model has been derived by using the non-standarard analysis combined with some a priori postulated assumptions. An application of the homogenization procedure leads to equations given in terms of unknown macro-temperature and thermal microlocal parameters. The model permits to evaluate not only mean but also local values of heat fluxes in every material component of the composite. Moreover, the model satisfies the continuity conditions for temperature and heat flux vector on interfaces. Making reference to the microlocal modeling method given in ( where h(r, z) is an unknown function (called the macrotemperature and being the averaged temperature), c(r, z) stands for the unknown thermal microlocal parameter, and h(z) is a given l periodic function (called the shape function): hðzÞ ¼ z À 0:5l 1 ; for 0 z l 1 ; Àgz 1Àg À 0:5l 1 þ l 1 1Àg ; for l 1 z l; The shape function h (Á) satisfies the following conditions hðzÞ 2 0ðlÞ; h 0 ðzÞ 2 0ð1Þ; hhi 1 l Bearing in mind the relations (2), (4), the underlined term in (1) is small for small l and will be neglected, and the following approximate relations can be written where h 0 (z) is the derivative of function h(z) taking the value 1 if a = 1 and -g/(1 -g) if a = 2. The governing equations of the homogenized model with microlocal parameters in the stationary axisymmetric case take the form (Matysiak and Perkowski [8]): From Eq. (6) it follows that the macro-temperature h satisfies the following equation The heat flux vector q (j) (r, z), j = 1, 2 in the layer of j-th kind is given by q ðjÞ ðr; zÞ ¼ ½q ðjÞ r ðr; zÞ; q ðjÞ z ðr; zÞ; q ðjÞ r ðr; From Eqs. (8) and (10) it is seen that the continuity conditions on interfaces are satisfied by the homogenized model.
The considered problem is determined by the following boundary conditions.
and lim r!1 From the mathematical reason caused by boundary condition (13) we divide the investigated layer on the two sublayers as well as the continuity conditions To solve the boundary value problems in sublayers S I and S II the Weber-Orr integral transforms will be applied (Titchmarsh [15], Krajewski and Olesiak [16], Olesiak [17] leads to the following relation (Olesiak [17]): W 00 ½B 0 ðh I ðr; zÞ; h II ðr; zÞÞ; r ! n ¼ À 2 p ðh I ða; zÞ; h II ða; zÞÞ À n 2 " h I n; z ð Þ; " h II n; z ð Þ À Á : The inverse transforms take the form ðh I ðr; zÞ; h II ðr; zÞÞ Transforming the Eq. (8) for h I and h II in S I and S II , respectively and using boundary conditions (16) and (18), we obtain the following ordinary differential equations From Eqs. (26), (27) it follows that where A, B, C, D are unknowns, which should be determined from boundary conditions (15), (17), (19) and (20). Satisfying these boundary conditions by " h I ðn; zÞ given in (28) and " h II ðn; zÞ given in (29) we obtain four linear algebraic equations for unknowns A, B, C, D and the solution takes the form and " h II ðn; zÞ ¼ 2h 0 To obtain the distribution of temperature h I (r, z) in S I , and h II (r, z) in S II in the integral forms we will substitute Eqs.
(30), (31) and (32) into (25). Thus, we have and " h II ðr; zÞ ¼ 2h 0 p The obtained integrals will be calculated numerically, and the results will be presented in the form of graphs. Moreover, knowing the temperature we will determine the distributions of heat flux. By using the relation (Olesiak [17], Olesiak and Krajewski [16]): and Eqs.     for 0 z \ h 1 : ð39Þ

Numerical results
The obtained results given by Eqs. (33), (34), (36)-(39) determine the temperature and heat flux distributions in the periodically stratified layer with the vertical cylindrical hole in the form of integrals, which will be calculated numerically. For this aim, the following dimensionless variables and parameters are introduced  Figure 4a, b present in the same manner as in Fig. 3 radial flux component q ðjÞ r =ðh 0 K Ã Þ as a function of g, g = l 1 /l for ĥ 1 ¼ ĥ 2 ¼ 1:0; r ¼ 1:05 and K 1 /K 2 = 4.0; 8.0. Figure 4a shows q ðjÞ r =ðh 0 K Ã Þ for ẑ ¼ Àĥ 2 =2 and Fig. 4b for The dimensional component of heat flux q ðjÞ z =ðh 0 K Ã Þ normal to the layering as a function of r for two cases of depths ẑ ¼ Àĥ 2 =2 and ẑ ¼ ĥ 1 =2, and two cases of ratio K 1 /K 2 = 4.0; 8.0 is presented in Fig. 5a, b. This component is continuous on the interfaces. The black curves represent the case of ĥ 1 ¼ ĥ 2 ¼ 1:0 the grey curves are suitable for ĥ 1 ¼ 2:0; ĥ 2 ¼ 1:0. Figure 6a, b show the dimensional heat flux q ðjÞ z =ðh 0 K Ã Þ as a function of parameter g for two cases of (a) Fig. 2 The isothermal lines in the periodically stratified layer ẑ; ẑ ¼ Àĥ 2 =2 and ẑ ¼ ĥ 1 =2, and K 1 /K 2 = 4.0; 8.0, r ¼ 1:05. The results for g = 0 and g = 1 are adequate for a homogeneous layer.

Final remarks
The temperature and heat flux distributions in the periodically stratified layer with cylindrical vertical hole is obtained within the framework of the homogenized model with microlocal parameters. The obtained results stand for some approximated solution to the considered problem. The boundary conditions are satisfied exactly, the equation of heat conduction with the periodically pice-wice constant coefficients is replaced by the equation with constant coefficients. The analysis of accuracy of the model were investigated in papers (Kulchytsky-Zhyhailo and Matysiak [18,19], Matysiak and Perkowski [20]), where the satisfactory consistency of the results obtained within the framework of homogenized model with microlocal  1903-1911 1909 parameters and within the classical heat conduction are confirmed.
Open Access This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.