CC BY
41
BeemHunJBry^T/Proceedings of VSUET, № 2, 20Î6L
Оригинальная статья/Original article_
УДК 517.951
DOI: http://doi.org/10.20914/2310-1202-2016-2-65-68
Нестационарное температурное поле в параллелепипеде в режиме теплопроводности при граничных условиях первого рода
Виталий К. Битюков, 1 president@vsuet.ru Анатолий А. Хвостов, 2 khvtol1974@yandex.ru Александра В. Сумина 3 vsum@rambler.ru
1 кафедра информ. и управл. систем, Воронеж. гос. ун т. инж. технол., пр-т Революции, 19, г. Воронеж, Россия
2 кафедра математики, ВУНЦ ВВС «ВВА им. профессора Н. Е. Жуковского и Ю. А. Гагарина», ул. Старых Большевиков, 54, г. Воронеж,
3 факультет прикладной математики, информатики и механики, Воронеж. гос. ун-т., Университетская пл. 1, г. Воронеж_
Реферат. Аналитическое изучение процессов теплопроводности является одним из основных разделов современных инженерных исследований в машиностроительной, энергетической, атомной промышленности, в технологических процессах химической, строительной, текстильной, пищевой, геологической и других отраслях промышленности. Достаточно указать, что практически все процессы в той или иной степени связаны с изменением температурного состояния и переносом теплоты. Следует также отметить, что инженерные исследования кинетики множества физических и химико-технологических процессов аналогичны задачам стационарной и нестационарной теплопроводности. К ним можно отнести процессы диффузий, седиментации, вязкого течения, замедления нейтронов, течения жидкостей через пористую среду, электрические колебания, сорбции, сушки, горения и др. Существуют различные методы решения классических краевых задач нестационарной теплопроводности и задач обобщённого типа: метод разделения переменных (метод Фурье); метод продолжений; метод произведения решений; метод Дюамеля; метод интегральных преобразований; операционный метод; метод функции Грина (для нестационарной и стационарной теплопроводности); метод отражения (метод источников). В данной работе на основе последовательного применения преобразования Лапласа по безразмерному времени в и конечного интегрального синус-преобразования по пространственным координатам X и Y решается задача нестационарного распределения температуры по механизму теплопроводности в параллелепипеде при граничных условиях первого рода. В результате получено аналитическое решение задачи распределения температуры в параллелепипеде для кондуктивного режима свободной конвекции, когда одна из боковых граней параллелепипеда
поддерживается при постоянной температуре, а остальные при другой одинаковой постоянной температуре._
Ключевые слова: аналитическое решение, конечные интегральные преобразования, теплопроводность, граничные условия
Nonstationary thermal field in the parallelepiped in the mode of heat _conduction under boundary conditions of first kind_
Vitaly K. Bityukov, 1 president@vsuet.ru Anatoly A. Khvostov, 2 khvtol1974@yandex.m Alexandra V. Sumina, 3 vsum@rambler.ru
1 Department of information and control systems, Voronezh State University of Engineering Technologies, Voronezh, Russia
2 mathematic department, MESC AF «N. E. Zhukovsky and Y. A. Gagarin Air Force Academy» (Voronezh), Staryh Bolshevikov str., 54, Voronezh
3 Faculty of Applied Mathematics, Informatics and Mechanics,Voronezh State University, 1 Universitetskaya pl., Voronezh, Russia_
Summary. Analytical study of the processes of heat conduction is one of the main topics of modern engineering research in engineering, energy, nuclear industry, process chemical, construction, textile, food, geological and other industries. Suffice to say that almost all processes in one degree or another are related to change in the temperature condition and the transfer of warmth. It should also be noted that engineering studies of the kinetics of a range of physical and chemical processes are similar to the problems of stationary and nonstationary heat transfer. These include the processes of diffusions, sedimentation, viscous flow, slowing down the neutrons, flow of fluids through a porous medium, electric fluctuations, adsorption, drying, burning, etc. There are various methods for solving the classical boundary value problems of nonstationary heat conduction and problems of the generalized type: the method of separation of variables (Fourier method) method; the continuation method; the works solutions; the Duhamel's method; the integral transformations method; the operating method; the method of green's functions (stationary and non-stationary thermal conductivity); the reflection method (method source). In this paper, based on the consistent application of the Laplace transform on the dimensionless time 0 and finite sine integral transformation in the spatial coordinates X and Y solves the problem of unsteady temperature distribution on the mechanism of heat conduction in a parallelepiped with boundary conditions of first kind. As a result we have the analytical solution of the temperature distribution in the parallelepiped to a conductive mode free convection, when one of the side faces of the parallelepiped is maintained at a constant temperature, and the others with the another same constant temperature._
Keywords: Analytical solution, finite integral transform, heat conduction, boundary conditions of first kind.
Для цитирования Битюков В. К., Хвостов А. А., Сумина А. В. Нестационарное температурное поле в параллелепипеде в режиме теплопроводности при граничных условиях первого рода // Вестник ВГУИТ. 2016. № 2. С. 65-68. (М:10.20914/2310-1202-2016-2-65-68
For citation
Bityukov V. K., Khvostov A. A., Sumina A. V. Nonstationary thermal field in the parallelepiped in the mode of heat conduction under boundary conditions of first kind. Vestnik VSUET [Proceedings of VSUET]. 2016. No. 2. pp. 65-68 (in Russ.). doi:10.20914/2310-1202-2016-2-65-68
Введение
При анализе кондуктивно-ламинарной свободной конвекции в замкнутых объёмах уравнения Обербека-Буссинеска могут быть представлены в несопряжённом виде [1]. Для плоских и осесимметричных геометрий в [2, 3] имеется множество решённых задач теплопроводности, но следует отметить, что анализ пространственных постановок представлен недостаточно. В [4] приведены методы получения аналитических решений, однако метод конечных интегральных синус-преобразований является наиболее эффективным [5].
В связи с этим, на примере области в форме параллелепипеда при граничных условиях первого рода демонстрируется алгоритм получения аналитического решения с помощью преобразования Лапласа по временной переменной и конечного интегрального синус-преобразования по геометрическим координатам.
Рассматривается задача нестационарного распределения температуры по механизму теплопроводности в параллелепипеде при граничных условиях первого рода.
д/ 1
^д 2t
дг pCp
д 2t
д 2t ^
дх2 дУ2 дzг
t (0, y, z, в) = t{, t (x, y, z,0) = 0; t(hl, y, z, в) = t(x ,0, z, в) = t (x, h2, z, в) = = t(x, y,0, в) = t (x, y, h3, в) = 0. Введём безразмерные координаты:
(1)
(2) (3)
II ;y и Z =-; a hu
h2 ; b ii h h ;0 = атт h12
Л
Рср
T = -L ti
тогда система (1)-(3) примет следующий вид:
dT д 2T ' '
дв
- + A
2 д2t + b2 д2t
(4)
дх 2 ду 2 дг 2
Т (Х ,У, г ,0)= 1;т (х, У, г ,0) = 0; (5)
т (1, у , г, в) = т (х ,0, г , в) = т (х ,1, г , в)= = т (х , У ,0, в)=т (х , У ,1, в)=0. Применение преобразования Лапласа по безразмерному времени 9 переводит (4)-(6) в систему для изображения
(6)
TL (X ,Y, Z, 5 )=L-1 [T (X ,Y, Z ,0) ]
d 2T
дХ2
- + A'
d T
дГ 2
+ B
d 2T
3Z
L = sT ■
2 Л L -
TL (0, Y, Z, s) = -;
s
TI (1, Y, Z, s) = TI (X ,0, Z, s) = TI (X ,1, Z, s) = = TI (X, Y ,0, s) = TI (X ,Y ,1, s) = 0.
(7)
(8) (9)
Применение конечного интегрального синус-преобразования по безразмерной коор-
Y •
динате
Fr [TL (X, Y, Z )] = Fr (X, Z ) = j T (X ,Y, Z )sin(/Y )dY,
0
где ^ - корни уравнения sin л = 0 ; с учётом
того, что
Fr
d % dY2
= "Л2Ф Y; Fy
Fv
d 2Tl dZ 2
d T
dX2
d2 Ф Y dZ2
d2 Ф Y dX2
1
Fy [l] =--(cos / —1);
Fy [s— ]= —— (cos/ — l); V
переводит исходную задачу в изображение:
d 2Ф„
- + B2
d 2 Ф¥
-(A2/2 + s )фг = 0; (10)
dX 2 dZ2
ФY (0, Z, s)=--— (cos / — l); (11)
V
ФY (l, Z, s )= ФY (X ,0, s)= ФY (X ,1, s) = 0. (12) Повторное применение синус-преобразования по безразмерной координате Z преобразует систему (10)—(12) в систему обыкновенных дифференциальных уравнений следующего вида:
d 2ф" — (B2y2 + AV + s )Ф XY = 0;
dX2
ФXY (0, s ) = -
(cos л — 1) (cos у -1)
VYs
Фyz (1, s) = 0.
Решение (13)—(15) таково:
2 ( cos / —1)( cosf —1)
(13)
(14)
(15)
ФYZ ( X, s ) = --
Л7
shy¡B У + A2 л2 + s (х -1) (16)
s * sh (2^1 B Y + A2/2 + s
Применим обратное преобразование Лапласа. Для этого воспользуемся теоремой Ващенко-Захарченко [2]:
L-1
p(X, y, s) w(л,y, s) _
p(X, л, Y,0) ((Л, Y,0)
.¿(wj x ехр(5й0).
((Л,Y, sn )
В нашем случае:
L"1
2 (cos W-1)( cos^-l)
Ф( Хф,У, s)
s) J WY
'sh^B'y2 + A2^ (x -1) sh (blB2r2 + A2W)
в sb,(X-1)
+2£(Д'+4ВУ+4 AW) cos A
xexp
-A - в2у2 - a2w
4
где fík - корень уравнения sin pk = 0. Переходя последовательно к оригиналу в (17), получим окончательное решение (рисунок 1):
(eos Мт - l)(cos Yn -1) ^
T ( X, Y, Z ,в) = ~
2 2
(17)
Ish^B2/2 + A2W (x -1)
sh (2^/B2y2 + A2^2)
A (X -1)
Pk sin
(18)
+2£(Pkk + 4B2 y2 + 4 A2 W) cos Pk
<exp
-A - b 2Г2 - A V2 |*é>
< sin (/JmY) x sin (y„Z)
(a) (b) (c)
Рисунок 1. Распределение температур при фиксированных значениях h1 = 0,1; h2 = 0,5; h3 = 20; в = 1 и: a) X = 0,05; b) Y = 0,25; c) Z = 0,25
Figure 1. Temperature distribution at fixed values h1 = 0,1; h2 = 0,5; h3 = 20; в = 1 и: a) X = 0,05; b) Y = 0,25; c) Z = 0,25
Заключение
В работе продемонстрировано применение конечного-интегрального преобразования для получения решения поставленной задачи новой структуры.
ЛИТЕРАТУРА
1. Latif M. Heat convection. New York: Springer, 2009. 552 р.
2. Цветков Ф. Ф., Григорьев Е. А. Тепломассообмен. М.: Изд-во МЭИ, 2011. 550 с.
3. Полянин А. Д., Зайцев В. Ф., Журов А. И. Методы решения нелинейных уравнений математической физики и механики. М.: Физ-матлит, 2011. 256 с.
4. Courant R., Hilbert D. Methods of mathematical physics. V. 2. Partial differential equatins. Singalore: Wiley - VCH, 1989. 896 р.
5. Полянин А. Д., Зайцев В. Ф. Справочник по нелинейным уравнениям математической физики: Точные решения. М.: Физмат-лит, 2002. 432 с.
6. Duffy D. G. Transform methods for solving partial differential equations, second edition. Chapman and Hall/CRC, 2004
REFERENCES
1 Latif M. Heat convection. New York, Springer, 2009. 552р.
2 Tsvetkov F. F., Grigoriev E. A. Tep-lomassoobmen [Heat and mass transfer] Moscow, MEI, 2011. 550 p. (in Russian).
x
x
Для связи с редакцией: vestnikvgta@mail.ru
67
3 Polyanin A. D., Zaitsev V. F., Zhurov A. I. Metody resheniya nelineinykh uravnenii matematich-eskoi fiziki i mekhaniki [Methods of solving nonlinear equations of mathematical physics and mechanics]. Moscow, Fizmatlit, 2011. 256 p. (in Russian).
4 Courant R., Hilbert D. Methods of mathematical physics. V. 2. Partial differential equatins. Singalore, Wiley - VCH, 1989. 896 р.
СВЕДЕНИЯ ОБ АВТОРАХ
Виталий К. Битюков д.т.н., профессор, президент, Воронеж. гос. ун-т. инж. техн., пр-т Революции, 19, г. Воронеж, Россия, president@vsuet.ru Анатолий А. Хвостов д.т.н., профессор, кафедра математики, ВУНЦ ВВС «ВВА имени профессора Н. Е. Жуковского и Ю. А. Гагарина» (г. Воронеж), ул. Старых Большевиков, 54 «А», г. Воронеж, Россия, khvtol1974@yandex.ru Александра В. Сумина студент, факультет прикладной математики, информатики и механики, Воронеж. гос. ун-т., Университетская пл. 1, г. Воронеж, vsum@rambler.ru
КРИТЕРИЙ АВТОРСТВА
Анатолий А. Хвостов обзор литературных источников по исследуемой проблеме, провёл эксперимент, выполнил расчёты
Виталий К. Битюков консультация в ходе исследования Александра В. Сумина написала рукопись, корректировала её до подачи в редакцию и несёт ответственность за плагиат
КОНФЛИКТ ИНТЕРЕСОВ
Авторы заявляют об отсутствии конфликта интересов. ПОСТУПИЛА 24.03.2016 ПРИНЯТА В ПЕЧАТЬ 21.04.2016
5 Polyanin A. D., Zaitsev V. F. Spravoch-nik po nelineinym uravneniyam matematicheskoi fiziki [Handbook of nonlinear equations of mathematical physics: Exact solutions]. Moscow, Fizmatlit, 2002. 432 p. (in Russian).
6 Duffy D. G. Transform methods for solving partial differential equations, second edition. Chapman and Hall/CRC, 2004
INFORMATION ABOUT AUTHORS
Vitaly K. Bityukov D. t. sc., professor, president, Voronezh state university of engineering technologies, Revolution Av., 19 Voronezh, Russia, president@vsuet.ru Anatoly A. Khvostov D. t. sc., professor, mathematic department, MESC AF «N. E. Zhukovsky and Y. A. Gagarin Air Force Academy» (Voronezh), Staryh Bolshevikov street, 54 «A», Voronezh, Russia, khvtol1974@yandex.ru Alexandra V. Sumina student, Faculty of Applied Mathematics, Informatics and Mechanics, Voronezh State University, 1 Universitetskaya, Voronezh, Russia, vsum@rambler.ru
CONTRIBUTION
Anatoly A. Khvostov review of the literature on an investigated problem, conducted an experiment, performed computations
Vitaly K. Bityukov consultation during the study Alexandra V. Sumina wrote the manuscript, correct it before filing in editing and is responsible for plagiarism
CONFLICT OF INTEREST
The authors declare no conflict of interest.
RECEIVED 3.24.2016 ACCEPTED 4.21.2016
