CC BY
72
Видин Ю.В., Казаков Р.В. Расчет распределения температуры...
УДК 536.24
Ю.В. Видин, Р.В. Казаков
Расчет распределения температуры в ребре постоянного поперечного сечения при радиационном отводе тепла с его поверхности
Yu.V. Vidin, R.V. Kazakov
Calculation of temperature distribution in rib with constant cross-section and radiative heat transfer on surface
Аннотация. Разработан приближенный аналитический метод расчета распределения температуры вдоль стержня постоянного поперечного сечения при радиационном теплообмене на его поверхности. Предложенный метод основан на получении нижней и верхней оценок искомого поля температуры. При этом используется интегральное линеаризующее преобразование, благодаря которому удалось существенно уменьшить влияние нелинейного члена в исходном дифференциальном уравнении. Полученные аналитические выражения позволяют сравнительно просто оценить максимальное и минимальное значение температуры по длине ребра. При умеренных значениях радиационного числа Старка интервал между верхней и нижней границами оказывается небольшим. При повышенных числах Старка расхождение между граничными значениями несколько возрастает. Однако средняя арифметическая величина между этими граничными температурами получается не слишком отличающейся от действительной.
Abstract. The approximate analytical calculation method of temperature distribution along the rib with constant cross-section and radiative heat transfer on surface has been developed. The method is based on calculation of the lowest and highest parameters of the temperature field. By using the integral linearizing transformation the authors considerably decrease impact of the nonlinear term in the basic differential equation. The found analytical formulas have allowed comparatively easy calculate the maximum and minimum temperature indexes along the rib. The interval between the minimum and maximum index is minor with the minor index of the Stark number. With the Stark number increasing the interval between the minimum and maximum increases. But the arithmetical mean of these boundary temperatures is not too different from the real one.
Ключевые слова: температурное поле, радиационный теплообмен, граничные значения, аналитический метод, ребристая поверхность
Key words: temperature field, radiative heat transfer, boundary conditions, analytical method, ribbed surface
1. Введение
В монографии (Керн, Краус, 1977) представлен всесторонний математический анализ процессов переноса тепла через различные стенки, оснащенные ребрами разной конфигурации. При этом рассмотрены главным образом задачи, относящиеся к классу линейных. Однако на практике имеют место случаи, когда тепловой поток, отдаваемый ребристой поверхностью, является нелинейной функцией его температуры. Так, например, возможен процесс радиационного охлаждения в среде с температурой близкой к нулю (в условиях вакуума).
Такая задача, представляющая определенный инженерный интерес, может быть теоретически сформулирована следующим образом
dx2 X Р ’ (1)
Т= Т0 прих = 0, (2)
dT п ; — = 0 при х = 1, dx (3)
где Т = T(x) - искомое распределение температуры вдоль оси ребра, К; Т0 - абсолютная температура основания ребра, К; X - коэффициент теплопроводности материала, Bm/мК; с - видимый коэффициент теплообмена излучением, Bm/^K4; P, f - периметр и площадь поперечного стержня соответственно, м, м2; l - длина ребра, м.
130
Вестник МГТУ, том 18, № 1, 2015 г.
стр. 130-133
2. Результаты исследований
Целесообразно систему (1-3) привести к безразмерному виду.
Если ввести безразмерные комплексы Т = —; X = — ; Sk- SEl— . где Sk - число Старка, то
Т
/
У
система (1-3) преобразуется к виду
= Sk&
d2$ dX2
А = 1 приХ = О,
d& „ v л --= 0 приХ = 1.
dX
(4)
(5)
(6)
Так как дифференциальное уравнение (4) является существенно нелинейным, получить его строгое аналитическое решение весьма затруднительно. Поэтому целесообразно применить к задаче (4-5) наиболее эффективный приближенный аналитический подход. Для ослабления нелинейности в уравнении (4) может быть применен метод интегрального линейного преобразования, предложенный в работе (Видин, 1992). Для этого нужно ввести новую зависимую переменную, связанную с искомой температурой fl (X) соотношением
(7)
1 Л 3
Тогда система (4-6) запишется следующим образом
-* = о,
dX2 \dX )
U = 0 приХ= О,
dU
--= 0 приХ= 1.
dX
dU
(8)
(9)
(10)
Если в первом приближении функцией F(X) = 4!4 ----- , входящей в уравнение (8)
dX
пренебречь, то нетрудно получить решение в форме Отсюда с учетом (7) следует, что
и = SkX| у-1
1
1 + ЗЖ| 1-f
(11)
(12)
Так как комплекс F(X) = 403f---1 является на отрезке (X = 0-1) положительной величиной,
\dX )
т.е. представляет собой некоторый условный дополнительный тепловой источник в стержне, то, очевидно, что выражение (12) дает заниженное значение искомого распределения температуры fl (X) при X > 0. Наименьшая величина безразмерной температуры ребра будет на его вершине (при X = 1). Эта температура согласно (12) равна
1
(13)
3
l + -Sk 2
Если теперь в уравнение (8) подставить вместо комплекса /' (X) = 4ХI-I приближенное
dX
значение F(X) = 4Q3imJ^
оно примет вид
d2U
dX2 , 3
1 + — Sk
— I -Sk = 0. dX
(14)
2
2
4
131
Видин Ю.В., Казаков Р.В. Расчет распределения температуры...
Обозначим
W =
dU dX '
причем
(15)
(16)
(17)
dX 2 + 3Sk
Интегрируя дифференциальное уравнение (17) с учетом условия (16) можно получить следующее аналитическое решение
W = 0 при X = 1.
Для новой переменной W дифференциальное уравнение (14) запишется
dW
8
-W -Sk = 0.
W = .
[(2 + 3Sk)Sk 6ХР 4^| q ч i
2 + 3 Sk(l }
8 ехр- 4^ 2Sk q ч ! ^
2 + 3 Sk
(18)
Формул (18) также можно представить в более удобном виде через гиперболический тангенс
(Бронштейн, Семендяев, 1965)
w_J(2 + 3Sk)Skth
(19)
Затем на основе выражения (15) легко установить зависимость для функции U (Бронштейн, Семендяев, 1965)
ТТ 2 + 3 Sk,
U =--------In
8
(20)
| 2 + 3Sk
Используя это решение, удается окончательно установить нижнюю границу для искомого температурного поля рассматриваемого тела
1
(21)
chi-8^(l-X)
_ \2 + 3 Sk _
Аналогичным образом находится верхняя граница для функции д = д (X). С этой целью уравнение (8) записывается в виде
d2U л(dU
---+ + 4 -
dX2 {dX
-Sk = 0.
(22)
т.е. фиктивный тепловой источник в стержне F(X) берется по максимально возможной величине. Интегрируя систему (22), (9-10) подобно тому, как это было изложено ранее, находится новое соотношение для переменной U
U = —In 4
ch2y[Sk (1 - X)
(23)
ch2yfsk
Отсюда следует, что формула для расчета верхних граничных значений искомого распределения температуры будет иметь вид
ch2*fsk chl^Ski 1-X)
(24)
Проведение практических расчетов по зависимостям (21) и (24) не представляет затруднений, так как гиперболические функции подробно затабулированы (Сегал, Семендяев, 1962).
В таблице представлены результаты вычислений безразмерной температуры в сечении X = 0,5 и X = 1 для чисел Старка Sk = 0,1; 0,5; 1,0.
2
132
Вестник МГТУ, том 18, № 1, 2015 г.
стр. 130-133
Таблица
Y Sk = 0,1 Sk = 0,5 Sk = 1,0
д KJ наим. д наиб. д KJ наим. д наиб. д KJ наим. д наиб.
0,5 0,963 0,968 0,881 0,892 0,813 0,843
1,0 0,957 0,957 0,848 0,858 0,767 0,795
3. Заключение
На основе результатов, представленных в таблице, можно сделать ряд выводов:
1. Различие между д наим. и д наиб. при небольших числах Sk почти отсутствует.
2. С ростом величины Sk различие между в Н£ШМ и в наиб возрастает. Однако по абсолютной величине невязка Д9 = &наиб - 0-наиш оказывается с инженерной точки зрения небольшой.
3.
Очевидно, что среднее арифметическое значение температуры 9
,9 +9
^Наиб. ' ^Hi
2
будет обладать
высокой точностью при любом значении числа Sk.
Литература
Бронштейн И.Н., Семендяев К.А. Справочник по математике. М., Наука, 1965. 609 с.
Видин Ю.В. Инженерные методы теплопроводности. Красноярск, Изд-во Красноярского государственного университета, 1992. 96 с.
Керн Д., Краус А. Развитые поверхности теплообмена. М., Энергия, 1977. 461 с.
Сегал Б.И., Семендяев К.А. Пятизначные математические таблицы. М., Государственное изд-во физико-математической литературы, 1962. 450 с.
References
Bronshtein I.N., Semendyaev K.A. Spravochnik po matematike [Mathematical guide]. M., Nauka, 1965. 609 p. Vidin Yu.V. Injenernie metodi teploprovodnosti [Engineering approach in thermal conductivity]. Krasnoyarsk, Izd-vo Krasnoyarskogo gosudarstvennogo universiteta, 1992. 96 p.
Kern D., Kraus A. Razvitie poverhnosi teploobmena [Developed heat exchange surfaces]. M., Energia, 1977. 461 p.
Segal B.I., Semendyaev K.A. Pyatiznachnie matematicheskie tablici [Five-place mathematical tables]. M., Gosudarstvennoe izd-vo fiziko-matematicheskoy literaturi, 1962. 450 p.
Информация об авторах
Казаков Роман Владимирович - Политехнический институт Сибирского федерального университета, кафедра теплотехники и гидрогазодинамики, канд. техн. наук, доцент, e-mail: [email protected]
Kazakov R.V. - Polytechnic Institute, Siberian Federal University, Department of Heat Engineering and Hydraulic Gas Dynamics, Cand. of Tech. Sci., Associate Professor, e-mail: [email protected]
Видин Юрий Владимирович - Политехнический институт Сибирского федерального университета, кафедра теплотехники и гидрогазодинамики, канд. техн. наук, профессор, e-mail: [email protected]
Vidin Yu.V. - Polytechnic Institute, Siberian Federal University, Department of Heat Engineering and Hydraulic Gas Dynamics, Cand. of Tech. Sci., Professor, e-mail: [email protected]
133
