Научная статья на тему 'THE PERIODIC TEMPERATURE OSCILLATIONS IN A CYLINDRICAL PROFILE WITH A LARGE THICKNESS'

THE PERIODIC TEMPERATURE OSCILLATIONS IN A CYLINDRICAL PROFILE WITH A LARGE THICKNESS Текст научной статьи по специальности «Физика»

CC BY
60
13
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Magazine of Civil Engineering
Scopus
ВАК
RSCI
ESCI
Область наук
Ключевые слова
ПОЛЫЙ ЦИЛИНДР / УРАВНЕНИЕ ТЕПЛОПРОВОДНОСТИ / ФУНКЦИЯ ХАНКЕЛЯ / КОНЕЧНО-РАЗНОСТНАЯ СХЕМА / ЦИЛИНДРИЧЕСКАЯ СИММЕТРИЯ / КОЭФФИЦИЕНТ ЗАТУХАНИЯ / ГРАЖДАНСКОЕ СТРОИТЕЛЬСТВО / ЗДАНИЕ / СТРОИТЕЛЬНАЯ ИНДУСТРИЯ / HOLLOW CYLINDER / THERMAL CONDUCTIVITY EQUATION / HANKEL FUNCTION / FINITE-DIFFERENCE SCHEME / TEMPERATURE WAVE / CYLINDRICAL SYMMETRY / DAMPING COEFFICIENT / CIVIL ENGINEERING / BUILDING / CONSTRUCTION INDUSTRY

Аннотация научной статьи по физике, автор научной работы — Samarin O.D.

A hollow cylinder with thick walls is one of the most complex objects to calculate the unsteady temperature field, so this field is the least studied. However, such objects are found in many modern constructions of systems of generation and distribution of heat. In the proposed work it deals with the study of propagation of temperature waves in the wall of the hollow cylinder with harmonic temperature change of external environment arising from its diurnal fluctuations. The approximate analytical solution is presented by separation of variables in the complex domain with the use of cylindrical functions. The algorithm of calculation of temperature fields numerically is shown using an explicit finite-difference scheme of high accuracy in conditions of cylindrical symmetry with boundary conditions of the first kind. The results of calculations according to the considered algorithm, depending on the time since the start of heat exposure and their comparison with the analytic solution are given for its implementation. Calculated radial profiles of the temperature in the cylindrical wall within the temperature waves and the analytical approximation relations for the description of its damping coefficient are presented. The results are compared with the existing analytical solution in rectangular coordinates and it is marked that they have some differences but the common results are found regardless of the material and geometry of the cylinder, as well as of temperatures of inner and outer environment. Presented dependences are invited to apply for the analytical evaluation of the temperature amplitude on the inner surface of the heated cylindrical structures that will allow the use of engineering methods to verify compliance with industrial safety requirements.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

ПЕРИОДИЧЕСКИЕ ТЕМПЕРАТУРНЫЕ КОЛЕБАНИЯ В ЦИЛИНДРИЧЕСКОМ СЛОЕ ПРИ БОЛЬШОЙ ТОЛЩИНЕ СТЕНКИ

Полый цилиндр с толстыми стенками является одним из наиболее сложных объектов для расчета нестационарного температурного поля, поэтому такое поле является наименее изученным. Вместе с тем подобные объекты встречаются во многих современных конструкциях систем генерации и распределения теплоты. В предлагаемой работе рассматривается исследование распространения температурной волны в стенке полого цилиндра при гармоническом изменении температуры наружной среды, возникающем при ее суточных колебаниях. Представлено приближенное аналитическое решение задачи методом разделения переменных в комплексной области с использованием цилиндрических функций. Показан алгоритм расчета температурного поля численным методом с помощью явной конечно-разностной схемы повышенной точности в условиях цилиндрической симметрии при граничных условиях первого рода. Приведены результаты вычислений по рассмотренному алгоритму в зависимости от времени с момента начала теплового воздействия и их сопоставление с аналитическим решением для осуществления его идентификации. Представлены рассчитанные радиальные профили температуры в стенке цилиндра в пределах температурной волны и предложены аналитические зависимости для коэффициента ее затухания. Полученные результаты сопоставлены с имеющимся аналитическим решением в прямоугольных координатах и отмечены их различия, а также общность найденных результатов независимо от материала и геометрии цилиндра, а также температур внутренней и наружной среды. Представленные зависимости предложено применять для аналитической оценки амплитуды колебаний температуры на внутренней поверхности цилиндрических нагреваемых конструкций, что позволит использовать инженерные методы проверки выполнения требований промышленной безопасности.

Текст научной работы на тему «THE PERIODIC TEMPERATURE OSCILLATIONS IN A CYLINDRICAL PROFILE WITH A LARGE THICKNESS»

Magazine of Civil Engineering. 2019. 85(1). Pp. 51-58 Инженерно-строительный журнал. 2019. № 1(85). С. 51-58

Magazine of Civil Engineering ISSN

° 2071-0305

journal homepage: http://engstroy.spbstu.ru/

DOI: 10.18720/MCE.85.5

The periodic temperature oscillations in a cylindrical profile with a large thickness

O.D. Samarin,

National Research Moscow State University of Civil Engineering, Moscow, Russia E-mail: [email protected]

Keywords: hollow cylinder; thermal conductivity equation; Hankel function; finite-difference scheme; temperature wave; cylindrical symmetry; damping coefficient; civil engineering; building; construction industry

Abstract. A hollow cylinder with thick walls is one of the most complex objects to calculate the unsteady temperature field, so this field is the least studied. However, such objects are found in many modern constructions of systems of generation and distribution of heat. In the proposed work it deals with the study of propagation of temperature waves in the wall of the hollow cylinder with harmonic temperature change of external environment arising from its diurnal fluctuations. The approximate analytical solution is presented by separation of variables in the complex domain with the use of cylindrical functions. The algorithm of calculation of temperature fields numerically is shown using an explicit finite-difference scheme of high accuracy in conditions of cylindrical symmetry with boundary conditions of the first kind. The results of calculations according to the considered algorithm, depending on the time since the start of heat exposure and their comparison with the analytic solution are given for its implementation. Calculated radial profiles of the temperature in the cylindrical wall within the temperature waves and the analytical approximation relations for the description of its damping coefficient are presented. The results are compared with the existing analytical solution in rectangular coordinates and it is marked that they have some differences but the common results are found regardless of the material and geometry of the cylinder, as well as of temperatures of inner and outer environment. Presented dependences are invited to apply for the analytical evaluation of the temperature amplitude on the inner surface of the heated cylindrical structures that will allow the use of engineering methods to verify compliance with industrial safety requirements.

1. Introduction

In this article, an unsteady temperature profile of a hollow thick wall cylinder under periodic thermal influence is studied.

The problem of unsteady thermal conductivity in bodies of various geometric forms is studied for a long time. Moreover, currently, due to achievements of computational equipment, numerical methods of solving the problem are more attractive. Nevertheless, in most cases, they are either related to a single dimension case or consideration is done in rectangular coordinates. This rather corresponds to the most really existing problems in both unsteady and steady modes [1-2]. Under the condition of cylindrical symmetry, heating and cooling of solid hollow cylinders are the most developed issues [1].

Over the last years a number of researches, where such issues are studied in both analytic and numerical ways, is emerged. However, the results obtained by the authors are either very difficult for use in engineering practice [3, 4, 9, 14, 17] or too rough, contrary [11]. The other solution variants are related to specific types of constructions used in limited areas and functioning in super critical modes [8], in nuclear power industry [13] or composite material production [5-7, 10], as well as in case of phase transitions [12], [15, 16, 18] or fuel combustion [19], or for underground pipes [20].

Samarin, O.D. The periodic temperature oscillations in a cylindrical profile with a large thickness. Magazine of Civil Engineering. 2019. 85(1). Pp. 51-58. DOI: 10.18720/MCE.85.5.

Самарин О.Д. Периодические температурные колебания в цилиндрическом слое при большой толщине стенки // Инженерно-строительный журнал. 2019. № 1(85). С. 51-58. DOI: 10.18720/MCE.85.5

("ОЕ^^^И This open access article is licensed under CC BY 4.0 (https://creativec0mm0ns.0rg/licenses/bv/4.0Л

In the same time, the calculation of an unsteady temperature profile of thick walled hollow cylinders, including a periodic mode, is interesting. For instance, this can be useful for solving the problem of condensate formation either on the internal surface of fume stacks or on the external surface of thermal pipe insulation with daily fluctuation of ambient air temperature and heat carrier temperature. This is especially interesting for open-laid pipelines of heating systems with a significant moisture content of the outside air, which happens at a temperature close to zero and above, when it may be necessary to assess the conditions of freezing and thawing of atmospheric moisture and snow. In addition, such a calculation is needed to confirm the absence of condensation on the inner surface of the air ducts of exhaust ventilation systems with their outer gasket, which in some cases is allowed in case of structural necessity with the condition of their thermal insulation. The same applies to the outer surface of the air intake ducts, when they are laid as an exception inside the building, also in a heat-insulated form. Finally, the considered calculation may be required to assess the temperature fluctuations on the outer surface of cylindrical furnaces with changes in their operation mode in order to solve the question of the admissibility of this temperature for sanitary and hygienic requirements. The author of the published paper [21] has obtained a simple solution for thick walled cylinder cooling. However, the solution is related only to aperiodic mode conditions, whereas an analytic solution, but for a small diameter linear source is given in the paper [22]. Therefore, the relevance of the proposed research is in the necessity to search the precise and physically well-reasoned (in the same time acceptable for engineering usage) dependences of temperature alteration in periodically heated/ cooled cylindrical structures. Obtained results may be acceptable for a large number of energetic facilities of such structure.

The calculation of a temperature profile at harmonic fluctuations of ambient temperature near the external surface of a cylinder is the research purpose. Research tasks are as follows:

- building up an analytical solution for an equation of thermal conductivity in the cylinder wall for a regular mode;

- development of an algorithm, which implements finite-difference approximation of the equation;

- obtaining analytic dependences for temperature across cylinder cross section and temperature wave amplitude by results of correlating the theoretical results with software generation data.

Figure 1 shows a diagram of studied cylinder and some conventional symbols used further. The main of them are ro and ri, i.e. the external and internal radii of a cylinder, m, respectively.

Figure 1. A diagram of a hollow cylinder and unsteady temperature profile in its wall.

Then a differential equation of unsteady heat conductivity for the cylindrical symmetry case may be written in the following way [1-2]:

2. Methods

r.

~7

dt _ a dt + d2t

(1)

where a = — - coefficient of cylinder body material heat conductivity, m2/s; cp

X is its thermal conductivity, W/(m K);

c and p are specific heat capacity, J/(kgK), and density, kg/m3, respectively. The solution in this case will show a current temperature profile in the temperature wave distribution zone.

The boundary condition on the cylinder external surface for a periodic mode in the complex view of the harmonic oscillation simplest case may be written as follows:

t = Ate exp( -/rax), (2)

where Ate is oscillation amplitude of ambient air temperature te, K;

ro = 2n/T is wave circular frequency of te, s-1. Here Tis oscillation period, s. As the initial condition we may adopt t(r) = 0, considering an average period value te.av = 0, and, if it is not, we shall consider excessive temperature 0 = t(r) - te.av. Therefore, we adopt conditions of 1-st type, and in case of known heat exchange intensity on the surfaces, in first approximation it is possible to introduce additional layers with thickness equal to krcond = X/a, where a - full heat exchange coefficient, W/(m2K), at a corresponding side.

Then we may search for solution of equation (1) by the method of dividing the variables as the production of cofactors depending respectively only on t and r:

t = exp(-irax)©(r') , (3)

where r' = ro - r is radial coordinate, with respect to the fact that a temperature wave distributes from outsides (see Figure 1). After substituting in (1) we obtained a regular differential equation for 9:

, d2© dm irar' _ r'—^ + + = (4)

dr 12 dr a

Such equation is Bessel's equation in the complex area. The equation may be solved by Hankel cylindrical function of zero order first type H01 [1] under applied boundary and initial conditions. The final shape of the function t may be found by substitution into (3):

t = Ae exp(-/rax)H01 (r, (5)

Through selecting a real component and using the known asymptotic approximation of the function H0(1) at sufficiently large values of the argument, which may be used, since we are interested in the established oscillation process, we finally find:

Re(') = ^ •■ cos- r^ + n) ^(^] . <6>

Therefore, unlike temperature waves in a flat wall, for a cylindrical wall it is more likely to have rapid decrease of amplitude as depth increases: not only due to the exponential factor, but also in a proportional

way lA/7. In this case, there is an additional shift of a phase in relation to oscillations te expressed by cosine

argument summand equal to n/n. The theoretical value of the parameter n makes 1/8, however, taking into account the approximated character of the dependency (6), it is reasonable to clarify it through comparing with numeric calculation data based on the finite-differential approximation of the equation (1).

Since we need the obtained result precision, whereas used memory volume and quantity of made operations are not significant due to high computational resources of modern computers, a finite-differential diagram will be the most favorable due to its programming simplicity. Then a temperature value in i-n grid node at j+1-n time moment may be calculated according to the following expression:

f f \ \ IV 1 /1 \ IV 1

ti, j+1 = Fo A

2i -1

v2-2 V--A / -- - y

ti+1, j +

1 .-2

A y

Fo

2i - 3 t;,j + 2-2 t/-1,j

(7)

Here ti,j, ti i j and ti+i,j are temperature values at j-n time moment in i-n node and two adjoining nodes on the right and left (node numbering from the cylinder axis to the external surface side);

аАт

FoA =-- is dimensionless local Fourier criterion, where Ax, s and Лг, m is respectively, time and

(Ar )2

coordinate steps representing parameters of finite-differential diagram. As it is known [1-2], the adopted diagram converges at Foa ^ 1/2. In this case, the value Foa = 1/6, providing high precision of 4-th order approximation by a space coordinate and 2-nd order approximation by a time coordinate, was chosen.

3. Results and discussion

Calculation results according to the formula (6) and author provided software in Fortran language according to the diagram (7) for Ate = 1 K, ro = 1 m, a = 5.1-10-7 m2/s and ю = 2п/86400 = 7.2710-5 s-1, i.e. with daily oscillation period are given in Figure 2. On the internal radius ri = 0.5 m, temperature was kept constant.

t, к

3

ч NN \>/ У\

V^2 \ 4

\\ \ \

\ \ 1 \

Figure 2. Temperature alteration along the cylinder radius (continuous lines - numeric calculation, dotted lines - according to the formula (6)); 1 - wx = 23.7; 2 - wx = 47.4; 3 - wx = 71.1.

In this case, for maximal matching it is necessary to adopt the value n = 32, i.e. real additional phase shift at cylindrical symmetry is very insignificant. It is clear, that in general, correctness of the relation (6) is

confirmed, except r' < 0.05. It is true that the function approximation H01 used in the conclusion (6) poorly works in this area, however for the considered problem it is meaningless, since in most cases a value of the relation Arcond/ro, showing relative radius of the additional external cylindrical layer with considered external heat exchange is more than 0.05 or, at least, equal to 0.05.

Since in practical problems, temperature fluctuation amplitude is interesting, it is possible to introduce its attenuation coefficient Vr = Ate/At, which may be associated with thermal inertia of a cylindrical layer Dr in the same way as in the single dimension problem. From (6) the following expression may be obtained for it:

VnDf expDrA/2) n - fc

vr = -----, where Dr = r — . (8)

2 \ a

It is not hard to see that the obtained dependency differs from the single dimension case [2] only by the

JnDr ,

additional factor — ■ It is clear that at Dr > 4/n = 1.26 this factor is more than 1, therefore damping in a cylindrical layer is really more significant. In the considered example for the entire cylinder

Dr =(1 - 0.5). I7-27 '10- = 5.965, then Vr = 147, or 2.16 times more than for a flat wall with the same heat r v 5.1-10

inertia. For illustrative purpose, the ratio between damping coefficients for a flat and cylindrical layers shown in Figure 3.

Since the ratio (8) is written in a dimensionless form, it is common, and having no dependency on actual air temperature values, cylinder radius, cylinder wall material, and the ratio ri/ro, at least at ri/ro > 0.5, for which numeric calculations are done.

Damping coefficient

160

140 120 100 80 60 40 2D 0

0 1 2 3 d 4 5 6 7

r

Figure 3. The dependence of the heat wave damping coefficient on heat inertia for a flat wall (dotted line) and cylindrical profile.

The results are principally matched with the data given in [2] for a flat wall with precision down to a factor

JnDr

—. Besides, a general form of calculated temperature profiles corresponds to results of some other

authors, for instance, [15], whereas their analytical description discovers similarities in theoretical solutions of other sources, in particular, [13], [14].

To illustrate the practical use of the ratio (8), we calculate the temperature fluctuations in the insulation layer of an open air duct with the parameters ri = 0.25 m, ro = 0.35 m, i.e. with the insulation thickness of 100 mm, with the amplitude of the outdoor air temperature fluctuations Ate = 2.7 °C for the average conditions of the heating season. In this case the condition ri/ro > 0.5 is satisfied. For thermal insulation of mineral wool X = 0.044 W/(mK), c = 840 J/(kgK), p = 50 kg/m3, so a = 1.0510-6 m2/s. Since we are primarily interested in the temperature on the inner surface of the air duct t(ri) from the point of view of assessing the possibility of condensation of water vapor, then r' = ro - ri = 0.1 m; and when the daily period of oscillations according to the formula (8) Dr = 0.83, where Vr = 1.45, and therefore the oscillation amplitude of the fluctuations of t(ri) is equal to Ate/Vr = 2.7/1.45 = 1.86 °C. This means that when checking the adequacy of the adopted thickness of the insulation will need to take into account that the minimum during the day the value of t(ri) will be below the average of 1.86 degrees.

4. Conclusion

- It is noted that temperature wave distribution is controlled by the same laws as in the single dimension case, but with a slightly different attenuation coefficient and phase shift.

- It is proved that unlike temperature waves in a flat wall, for a cylindrical wall it is more likely to have rapid decrease of amplitude as depth increases: not only due to the exponential factor, but also in reverse proportion to square root of a radial coordinate.

- It is shown that a real additional phase shift of temperature fluctuation at cylindrical symmetry in comparison with the single dimension case is very insignificant and equals to 1/32 of a period.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

- It is shown that the ratio of the internal and external radii of a hollow cylinder does not explicitly influences on the temperature profile character in the temperature wave zone, at least, at ro/R > 0.5.

- It is shown that temperature field distribution in the zone of temperature wave penetration into the cylinder, according to results of analytic and numeric solutions, coincides within engineering calculation limits, which means that the obtained dependences are correct.

-It is proposed to apply the ratios obtained in the research for analytic evaluation of temperature fluctuation amplitude on cylindrical surfaces of heated and cooled structures. It is necessary, first of all, to solve the problem of condensate formation on the internal surface of fume stacks in case of boiler unit load

fluctuation or on the external surface of heat line insulation as well as on the surfaces of the open-laid heat-insulated air ducts in case of daily ambient temperature alteration and on the outer surface of cylindrical furnaces in variable operating modes. It will allow to use not only software methods, but also engineering methods to check meeting the industrial safety requirements.

References

1. Carslaw, H.S., Jaeger, J.C. Conduction of heat in solids. Oxford University Press, USA, 1986. 520 p.

2. Fokin, K.F. Stroitel'naya teplotekhnika ograzhdayushchikh chastey zdaniya [Thermal physics of building envelopes]. Under ed. of Yu.A. Tabunshchikov and V.G. Gagarin. 5th ed., revised. Moscow : AVOK-PRESS Publishers, 2006. 256 p. (rus)

3. Shokouhmand, H., Ghaffari, S. Investigation of cooling process of a high-temperature hollow cylinder in moving induction heat treatment. Journal of Thermal Science and Engineering Applications. 2012. Vol. 4. No. 1. Pp. 011002.

4. Kalynyak, B.M. Fredholm equations of the second kind for radial stresses aimed at the determination of the thermoelastic state of an inhomogeneous hollow long cylinder. Journal of Mathematical Sciences. 2015. Vol. 205. No. 5. Pp. 659-666.

5. Norouzi, M., Birjandi, A.K., Rahmani, H., Joneidi, A.A. A general exact analytical solution for anisotropic non-axisymmetric heat conduction in composite cylindrical shells. International Journal of Heat and Mass Transfer. 2016. Vol. 93. Pp. 41-56.

6. Norouzi, M., Shariati, M., Rezaei Niya, S.M., Karimi Demneh, M., Naghavi, M.S., Kayhani, M.H. Exact solution of unsteady conductive heat transfer in cylindrical composite laminates. Journal of Heat Transfer. 2012. Vol. 134. No. 10. Pp. 101301.

7. Shevchuk, V.A. Nonstationary one-dimensional problem of heat conduction for a cylinder with a thin multi-layer coating. Journal of Mathematical Sciences. 2012. Vol. 184. No. 2. Pp. 215-223.

8. Brunner, G. Heat transfer. Supercritical fluid science and technology. 2014. No. 5. Pp. 228-263.

9. Horikiri, K., Yao, Y., Yao, J. Modelling conjugate flow and heat transfer in a ventilated room for indoor thermal comfort assessment. Building and Environment. 2014. No. 77. Pp. 135-147.

10. Yun, T.S., Jeong, Y.J., Han, T.-S., Youm, K.-S. Evaluation of thermal conductivity for thermally insulated concretes. Energy and Buildings. 2013. No. 61. Pp. 125-132.

11. Gabrielaitiene, I. Numerical simulation of a district heating system with emphases on transient temperature behavior. Proceedings of the 8th International Conference «Environmental Engineering». Vilnius: VGTU Publishers. 2011. Vol. 2. Pp. 747-754.

12. Vasilyev, G.P., Gornov, V.F., Konstantinov, P.I., Kolesova, M.V., Korneva, I.A. Analysis of ground temperature variations, on the basis of years-long measurements. Magazine of Civil Engineering. 2017. 72(4). Pp. 62-72. DOI: 10.18720/MCE.72.8.

13. Loginov, V.S., Parpiev, A.T. Nestacionarnye temperaturnye rezhimy pologo cilindricheskogo aktivnogo elementa i paroperegrevatelya v parogeneratore [Non-stationary thermal regimes of the hollow cylindrical active element and the vapor superheater in the vapor generator]. Herald of the Sholom-Aleihem Amur State Univ. 2013. No. 2(13). Pp. 90-97 [Online]. URL: https://elibrary.ru/download/elibrary_21298568_86404092.pdf (date of treatment: 23.11.18). (rus)

14. Eliseev, V.N., Tovstonog, V.A., Borovkova, T.V. Algoritm resheniya obobshchennoy zadachi nestacionar-noy teploprovodnosti [Solution Algorithm of Generalized Non-Stationary Heat Conduction Problem in the Bodies of Simple Geometric Shapes]. Herald of the Bauman Moscow State Tech. Univ., Mech. Eng. 2017. No. 1. Pp. 112-128 [Online]. URL: https://elibrary.ru/download/ elibrary_28415996_48648202.pdf (date of treatment: 23.11.18). (rus)

15. Vasilyev, G.P., Gornov, V.F., Peskov, N.V., Popov, M.P., Kolesova, M.V., Yurchenko, V.A. Ground moisture phase transitions: Accounting in BHE'S design. Magazine of Civil Engineering. 2017. 74(6). Pp. 102-117. DOI: 10.18720/MCE.74.9

16. Parfentyeva, N., Valancius, K., Samarin, O., Paulauskaite, S., Ciuprinskine J. Solving the problem of pipeline freezing with respect to external heat exchange. Mechanika. 2015. Vol. 21 (5). Pp. 393-396.

17. Liu, C.-S. An integral equation method to recover non-additive and non-separable heat source without initial temperature. International Journal of Heat and Mass Transfer. 2016. Vol. 97. Pp. 943-953.

18. Bahadori, A., Vuthaluru, H.B. Simple method for estimation of unsteady state conduction heat flow with variable surface temperature in slabs and spheres. International Journal of Heat and Mass Transfer. 2010. Vol. 53. Pp. 4536-4542.

19. Li, Y., Kong, S.-C. Coupling conjugate heat transfer with in-cylinder combustion modeling for engine simulation. International Journal of Heat and Mass Transfer. 2011. Vol. add.

20. Bouhacina, B., Saim, R., Oztop, H.F. Numerical investigation of a novel tube design for the geothermal borehole heat exchanger. Applied Thermal Engineering. 2015. No. 79. Pp. 153-162.

21. Samarin, O.D. The temperature waves motion in hollow thick-walled cylinder. Magazine of Civil Engineering. 2018. 78(2). Pp. 161-168. DOI: 10.18720/MCE.78.13

22. Podledneva, N.A., Krasnov, V.A., Magomadov, R.S. Opredeleniye koefficientov teploprovodnosti i tem-peraturoprovodnosti za odin opyt metodom lineynogo istochnika teploty postoyannoy moshchnosti [De-termination of coefficients of heat conductivity and heat diffusivity for one experience by a method of a linear source of heat of constant power]. Herald of Astrakhan State Technical University. 2013. No 2. Pp. 50-55 [Online]. URL: https://elibrary.ru/download/elibrary_20859905_11088693.pdf (date of treatment: 23.11.18). (rus)

Contacts:

Oleg Samarin, +7(916)107-77-01; [email protected]

© Samarin, O.D., 2019

Magazine of Civil Engineering

journal homepage: http://engstroy.spbstu.ru/

ISSN

2071-0305

РСН: 10.18720/МСЕ.85.5

Периодические температурные колебания в цилиндрическом слое при большой толщине стенки

0.Д. Самарин,

Федеральное государственное бюджетное образовательное учреждение высшего образования «Национальный исследовательский Московский государственный строительный университет», Москва, Россия Е-mail: [email protected]

Ключевые слова: полый цилиндр; уравнение теплопроводности; функция Ханкеля; конечно-разностная схема; цилиндрическая симметрия; коэффициент затухания; гражданское строительство; здание; строительная индустрия

Аннотация. Полый цилиндр с толстыми стенками является одним из наиболее сложных объектов для расчета нестационарного температурного поля, поэтому такое поле является наименее изученным. Вместе с тем подобные объекты встречаются во многих современных конструкциях систем генерации и распределения теплоты. В предлагаемой работе рассматривается исследование распространения температурной волны в стенке полого цилиндра при гармоническом изменении температуры наружной среды, возникающем при ее суточных колебаниях. Представлено приближенное аналитическое решение задачи методом разделения переменных в комплексной области с использованием цилиндрических функций. Показан алгоритм расчета температурного поля численным методом с помощью явной конечно-разностной схемы повышенной точности в условиях цилиндрической симметрии при граничных условиях первого рода. Приведены результаты вычислений по рассмотренному алгоритму в зависимости от времени с момента начала теплового воздействия и их сопоставление с аналитическим решением для осуществления его идентификации. Представлены рассчитанные радиальные профили температуры в стенке цилиндра в пределах температурной волны и предложены аналитические зависимости для коэффициента ее затухания. Полученные результаты сопоставлены с имеющимся аналитическим решением в прямоугольных координатах и отмечены их различия, а также общность найденных результатов независимо от материала и геометрии цилиндра, а также температур внутренней и наружной среды. Представленные зависимости предложено применять для аналитической оценки амплитуды колебаний температуры на внутренней поверхности цилиндрических нагреваемых конструкций, что позволит использовать инженерные методы проверки выполнения требований промышленной безопасности.

Литература

1. Carslaw H.S., Jaeger J.C. Conduction of heat in solids. Oxford University Press, USA, 1986. 520 p.

2. Фокин К. Ф. Строительная теплотехника ограждающих частей здания. / Под ред. Ю.А. Табунщикова, В.Г. Гагарина. 5-е изд., пересмотр. М.: АВОК-ПРЕСС, 2006. 256 с.

3. Shokouhmand H., Ghaffari S. Investigation of cooling process of a high-temperature hollow cylinder in moving induction heat treatment // Journal of Thermal Science and Engineering Applications. 2012. Vol. 4. No. 1. Pp. 011002.

4. Kalynyak B.M. Fredholm equations of the second kind for radial stresses aimed at the determination of the thermoelastic state of an inhomogeneous hollow long cylinder // Journal of Mathematical Sciences. 2015. Vol. 205. No. 5. Pp. 659-666.

5. Norouzi M., Birjandi A.K., Rahmani H., Joneidi A.A. A general exact analytical solution for anisotropic non-axisymetric heat conduction in composite cylindrical shells // International Journal of Heat and Mass Transfer. 2016. Vol. 93. Pp. 41-56.

6. Norouzi M., Shariati M., Rezaei Niya S.M., Karimi Demneh M., Naghavi M.S., Kayhani M.H. Exact solution of unsteady conductive heat transfer in cylindrical composite laminates // Journal of Heat Transfer. 2012. Vol. 134. No. 10. Pp. 101301.

7. Shevchuk V.A. Nonstationary one-dimensional problem of heat conduction for a cylinder with a thin multi-layer coating // Journal of Mathematical Sciences. 2012. Vol. 184. No. 2. Pp. 215-223.

8. Brunner G. Heat transfer // Supercritical fluid science and technology. 2014. No. 5. Pp. 228-263.

9. Horikiri K., Yao Y., Yao J. Modelling conjugate flow and heat transfer in a ventilated room for indoor thermal comfort assessment // Building and Environment. 2014. No. 77. Pp. 135-147.

10. Yun T.S., Jeong Y.J., Han T.-S., Youm K.-S. Evaluation of thermal conductivity for thermally insulated concretes // Energy and Buildings. 2013. No. 61. Pp. 125-132.

11. Gabrielaitiene I. Numerical simulation of a district heating system with emphases on transient temperature behavior // Proceedings of the 8th International Conference «Environmental Engineering». Vilnius: VGTU Publishers. 2011. Vol. 2. Pp. 747-754.

12. Васильев Г.П., Горнов В.Ф., Константинов П.И., Колесова М.В., Корнева И.А. Анализ изменения температуры грунта на основе многолетних измерений // Инженерно-строительный журнал. 2017. № 4(72). С. 62-72. DOI: 10.18720/MCE.72.8.

13. Логинов В.С., Парпиев А.Т. Нестационарные температурные режимы полого цилиндрического активного элемента и пароперегревателя в парогенераторе // Вестник Приамурского государственного университета им. Шолом-Алейхема. 2013. № 2(13). С. 90-97 [Электронный ресурс]. URL: https://elibrary.ru/download/elibrary_21298568_86404092.pdf (дата обращения: 23.11.18)

14. Елисеев В.Н., Товстоног В.А., Боровкова Т.В. Алгоритм решения обобщенной задачи нестационарной теплопроводности в телах простой геометрической формы // Вестник МГТУ им. Н.Э. Баумана. Сер. Машиностроение. 2017. № 1. С. 112-128 [Электронный ресурс]. URL: https://elibrary.ru/download/elibrary_28415996_48648202.pdf (дата обращения: 23.11.18)

15. Васильев Г.П., Горнов В.Ф., Песков Н.В., Попов М.И., Колесова М.В., Юрченко В.А. Фазовые переходы влаги в грунте. Учет при проектировании грунтовых теплообменников // Инженерно-строительный журнал. 2017. № 6(74). С. 102-117.

DOI: 10.18720/MCE.74.9

16. Parfentyeva N., Valancius K., Samarin O., Paulauskaite S., Ciuprinskine J. Solving the problem of pipeline freezing with respect to external heat exchange // Mechanika. 2015. Vol. 21 (5). Pp. 393-396.

17. Liu C.-S. An integral equation method to recover non-additive and non-separable heat source without initial temperature // International Journal of Heat and Mass Transfer. 2016. Vol. 97. Pp. 943-953.

18. Bahadori A., Vuthaluru H. B. Simple method for estimation of unsteady state conduction heat flow with variable surface temperature in slabs and spheres // International Journal of Heat and Mass Transfer. 2010. Vol. 53. Pp. 4536-4542.

19. Li Y., Kong S.-C. Coupling conjugate heat transfer with in-cylinder combustion modeling for engine simulation // International Journal of Heat and Mass Transfer. 2011. Vol. add.

20. Bouhacina B., Saim R., Oztop H.F. Numerical investigation of a novel tube design for the geothermal borehole heat exchanger // Applied Thermal Engineering. 2015. No. 79. Pp. 153-162.

21. Самарин О.Д. Распространение температурных волн в пустотелом толстостенном цилиндре // Инженерно-строительный журнал. 2018. № 2(78). С. 161-168. DOI: 10.18720/MCE.78.13

22. Подледнева Н.А., Краснов В.А., Магомадов Р.С. Определение коэффициентов теплопроводности и температуропроводности за один опыт методом линейного источника теплоты постоянной мощности // Вестник Астраханского государственного технического университета. 2013. № 2. С. 50-55 [Электронный ресурс]. URL: https://elibrary.ru/download/elibrary_20859905_11088693.pdf (дата обращения: 23.11.18).

Контактные данные:

Олег Дмитриевич Самарин, +7(916)107-77-01; эл. почта: [email protected]

© Самарин О.Д., 2019

i Надоели баннеры? Вы всегда можете отключить рекламу.