CC BY
39
УДК 536.253
Посохин В.Н. - доктор технических наук, профессор
E-mail: posohin@kgasu. ru
Зиганшин А.М. - кандидат технических наук, доцент
E-mail: amziganshin@kgasu. ru
Горохова А.Ю. - студент
E-mail: [email protected]
Казанский государственный архитектурно-строительный университет
Адрес организации: 420043, Россия, г. Казань, ул. Зеленая, д. 1
О естественной конвекции над горизонтальными теплоисточниками
Аннотация
Численным методом исследуется конвекция над горизонтальным протяженным теплоисточником заделанным заподлицо с окружающей поверхностью. Исследование проводится с использованием различных вариантов замыкания системы уравнений турбулентного движения, а также при помощи ламинарной модели. Находятся значения критических высот, разделяющих зоны ламинарного и турбулентного режимов, в зависимости от значения числа Рэлея (Ra) определяемого по характеристикам источника. Определяются осевые параметры - скорость, избыточная температура.
Ключевые слова: естественная конвекция, плоский теплоисточник, режимы течения, критическое значение числа Рэлея, модели турбулентности.
Схемы конвективных струйных течений вблизи вертикальной и горизонтальной пластин показаны на рис. 1. Вначале струя ламинарна, затем следует зона переходного режима течения и далее зона развитого турбулентного течения.
Рис. 1. Схемы конвективных течений: а) пристенная струя у вертикальной пластины; б) свободная струя над горизонтальной пластиной
Закономерности естественной конвекции у теплоисточников (интенсивность теплоотдачи, распределение температур и скорости, протяженность характерных зон течений) определяются значением числа Рэлея:
,Ж3
Ra = A-
T3
(1)
где АТ = ТП - Тш - разность температур поверхности теплоисточника и окружающего воздуха, К; или АТ = Тх - Тш - разность температур на оси струи в сечении х и окружающего воздуха, К; I - характерный размер, м, в качестве которого можно принять размер горизонтального теплоисточника Ь или координату х. Т - определяющая
температура, К; обычно принимают Т = (ТП + Тш)/2; А - размерный коэффициент, включающий в себя физические константы и характеристики воздуха. При атмосферном давлении Р = 101325 Па и Тш = 293 К, А « 2-1015 К2/м3.
Учитывая сказанное выше, мы можем сформировать число Рэлея в некотором произвольном сечении струи в виде:
Иах = Яа0 + АИа,, (2)
где
(3)
Ra0 = 2-10 --
DRax = 2-1015 X • (4)
Число Ra0 есть отношение сил гравитации и сил вязкости в области прилегающей непосредственно к горизонтальному теплоисточнику; добавка ARax характеризует изменение числа Рэлея по высоте струи. Для вертикальной пластины, очевидно, Ra0 = 0. Запишем еще выражение (1) в виде:
Ra x
где
Ra* = ^ =1 + AT*X , (5)
— T - T x
_ x ¥ _ x
x , x = — •
ТП - T¥ b
П
Сведения о влиянии числа Ra0, приводимые в известных руководствах по теплообмену, касаются, в основном, интенсивности теплоотдачи на поверхности горизонтально расположенного теплоисточника [1-6].
Данные о влиянии Ra0 на распределение скорости и температуры, протяженность зон ламинарного и переходного режимов, пороговые значения числа Rax, разделяющих эти режимы для течений над горизонтальными теплоисточниками практически отсутствуют.
Далее будем рассматривать струйное течение над горизонтальным протяженным теплоисточником - пластиной шириной 2b, заделанной заподлицо в ограничивающую поверхность.
Система уравнений двумерного ламинарного движения, дополненная уравнением газового состояния, лучистого теплообмена и граничными условиями, решалась численно с помощью пакета программ Fluent. На части границы моделирующей теплоисточник, в качестве граничного условия принималось равномерное распределение теплового потока Q0 (Вт/м), оставшаяся часть нижней горизонтальной границы -адиабатическая стенка, а остальные границы - проницаемы для потоков воздуха и тепла. В ходе численного эксперимента характерный размер теплоисточника оставался неизменным - 2b = 0,01 м, изменялся тепловой поток Q0 и значение числа Rao, которое находилось в пределах 8-102 ^ 104. Такие значения Ra0 согласно рекомендациям [1-6] соответствуют ламинарному режиму теплоотдачи. Ниже в табл. приведены условия численного эксперимента для четырех характерных случаев. Здесь введены обозначения: Qk - конвективная мощность теплоисточника, Вт/м;
им = ? —^ Q':T ~ ' q1 3 - скорость течения, м/с, принятая за масштаб;
\ Cp ' Р ¥ '
g - ускорение свободного падения, м/с2;
cp - удельная теплоемкость воздуха Дж/кг- К;
рш - плотность воздуха окружающего струю, кг/м3).
Таблица
Qo, Вт/м Qk, ВТ/М tст, К Мм, м/с Ra0
100 37,4 577,1 0,1 9,8-103
20 11,1 386,8 0,067 6,8-103
10 5,9 346,4 0,054 4,6-103
1 0,56 300,5 0,025 8,0-102
В процессе вычислений находились распределения продольной скорости их, избыточной температуры АТх и соответствующие значения Яах.
На рис. 2 показаны результаты вычислений относительных осевой скорости
их = их1им , осевой избыточной температуры АТх = АТх/ ДТП и числа Рэлея Яах.
i 6) irt I n
1 N ll
л ■4 \ 1 R a,
i i -----4,6^10' (J 4
•ч i' - - - SUM ■ - LR] i * j
\
-
iV
jj II t №
J ■^ Sv; ■ ч -. ._ /
90
КО
Ы1
50
V/
¿11
10
0 ^ 4 и.
0,1
fi2
AT, 10' 10 Ю1
Ra
Рис. 2. Изменение относительных характеристик по высоте струи: а - осевая скорость, б - осевая избыточная температура, в - текущее число Рэлея
Точками показаны кривые построенные по формулам И. А. Шепелева [8]: - осевая скорость:
§ • а
u3 =-
Cp Р® T¥
■ f \Х
(6)
где
fA-b\-
1 + Pr„„
Pr„„
i 42■c
erf
4Пс 43 • b \ 42 • cx
43 • ь
- осевая избыточная температура:
DTx -
T • £f,ix
2 2 Cp ^ Р»
(7)
где
f2 \ bx\-
8 ■ fi
erf
4 + Pr„ ■ b
42 ■ cx
Pr„ - турбулентное число Прандтля (согласно рекомендациям [7] принято равным 0,85); с - экспериментальная константа, рекомендуемое значение которой - 0,082.
Во всех исследованных случаях осевая скорость плавно возрастает с высотой. При некотором значении x = хкр возникают возмущения, которые далее нарастают, течение становится неустойчивым. Можно предположить, что такие возмущения соответствуют началу переходного режима, который затем переходит в развитый турбулентный. На
этом же рисунке показано как изменяется текущее значение числа Рэлея Ra x.
Здесь неустойчивость возникает на большей высоте. Объяснение этому факту состоит видимо в том, что Rax зависит только от избыточной температуры - формула (5). Из рис. 2 видно, что кинематическая неустойчивость течения наступает заметно раньше, чем начинаются колебания температуры, поэтому и колебания Rax начинаются на большей высоте.
На рис. 3 показаны связи между критической высотой хкр и значением числа Rao, определенные по показателям неустойчивости характеристик ux и Rax (DTx ) .
9Ü Sli
то
60
40
за
-щ—
8s-
о
■о - ПО 11, Ф ■ по ШмйТ,}
Кя.-ю'
Рис. 3. Изменение Хкр
В обоих случаях зависимость эта близка к линейной:
- по показателю кинематической неустойчивости:
хкр =-0,0062 • Ra0 + 99,
- по показателю температурной неустойчивости:
Хр =-0,0073 • Ra0 +115.
Приведем еще результаты расчета осевой скорости в струе, полученные при решении системы уравнений турбулентного движения. Для замыкания системы уравнений использовались три модификации k-s модели турбулентности - «стандартная» (SKE), «реализуемая» (RKE) и «ренормализованных групп» (RNGKE), и две модификации модели Рейнольдсовых напряжений устанавливающих связь между корреляциями давления и напряжений - линейной (RSM LPS) и квадратичной (RSM QPS). Результаты вычислений при Rao = 7,8-103 приведены на рис. 4, из которого видно, что использование RSM моделей приводит к существенному несоответствию с аппробированной формулой (6).
10
10' 10 10' Ra.
Ъ !
U i Í 4 1
Ra, iL Щ i i
1
ь N i i
[J ri l i
Ч1Т
-— RNUKE RKE FlSM LR RäM QP3 И Г1- ■V i 1 i
— j: i i
Г г J 1 i ¡ 1 J i i
j[i Iii (ft 1 t ■ V
¿f .
х
90 90 70 60 50 40 30 20 10
О
0,5 IJ0
lh5 2,0
Рис. 4. Сравнение результатов численного расчета с использованием различных моделей замыкания системы уравнений турбулентного движения
Расчеты, проведенные при разных значениях числа Рэлея (4-102 < Rao < 1-1011) показали, что скорость в струе возрастает с увеличением Ra0, однако зависимость эта не
так существенна. В практических расчетах можно воспользоваться формулой (6), в которой параметр Ra0 отсутствует и которая довольно хорошо соответствует результатам численных расчетов по k-s моделям. Заметим, что формула (6) хорошо соответствует результатам с использованием ламинарной модели при x < 10, далее расхождение становится все заметнее. При использовании турбулентных моделей все наоборот -формула (6) довольно сильно расходится с численным решением при x < 30, после -результаты аналитического и численного решений практически совпадают.
Список библиографических ссылок
1. Батурин В.В., Эльтерман В.М. Аэрация промышленных зданий. 2-е изд., испр. и доп. - М.: Госстройиздат, 1963. - 317 с.
2. Кутателадзе С. С. Теплопередача и гидродинамическое сопротивление: Справочное пособие. - М.: Энергоатомиздат, 1990. - 367 с.
3. Попов И. А. Гидродинамика и теплообмен внешних и внутренних свободноконвективных вертикальных течений с интенсификацией. Интенсификация теплообмена: монография / Под общ. ред. Ю.Ф. Гортышова. -Казань: Центр инновационных технологий, 2007. - 326 с.
4. Martynenko O.G., Khramtsov P.P. Free-Convective Heat Transfer // Free. Heat Transf. With Many Photogr. Flows Heat Exch. Berlin/Heidelberg: Springer-Verlag, 2005. - 516 p.
5. Corcione M. Heat transfer correlations for free convection from upward-facing horizontal rectangular surfaces // WSEAS Trans. Heat Mass Transf., 2007, Vol. 2, № 3. - P. 48-60.
6. Уонг Х. Основные формулы и данные по теплообмену для инженеров / Справочник. - М.: Атомиздат, 1979. - 216 с.
7. ANSYS FLUENT 6.3 Documentation / 12.3.8 Convective Heat and Mass Transfer Modeling. URL: https://www.sharcnet.ca/Software/Fluent6/html/ug/node476.htm (дата обращения: 27.02.2016).
8. Романов С.В., Посохин В.Н., Зиганшин А.М. О влиянии способа подвода тепла на результаты численного решения задачи о свободной конвекции над объёмным теплоисточником // Известия высших учебных заведений. Проблемы энергетики, 2013, № 3-4. - С. 48-53.
Posokhin V.N. - doctor of technical sciences, professor E-mail: posohin@kgasu. ru
Ziganshin A.M. - candidate of technical sciences, associate professor E-mail: amziganshin@kgasu. ru Gorohova A.Y. - student E-mail: [email protected]
Kazan State University of Architecture and Engineering
The organization address: 420043, Russia, Kazan, Zelenaya st., 1
About natural convection over the horizontal heat sources
Resume
This work presents results of numerical research of natural convection above the horizontal plate flush with the surrounding surface. Arising jet can be divided into three sections - initial where the flow has a laminar flow regime, continue to increase the height of the flow begins to lose stability - a zone of transition flow, and with a certain height - developed turbulence jet. The concept of the Rayleigh number for an arbitrary section Rax, as the sum of Ra0 - shows the ratio of gravity forces and viscous forces in the area adjacent to the horizontal heat source, and ARax additive characterizes the change in the number of Rayleigh on the jet height.
Conducted a series of numerical studies using laminar model i.e. solve a system of equations of fluid motion without the use of RANS approach. Solving the tasks with changing of heat flux of the source, so that Ra0 was within 8-102^104, which corresponds to a laminar
mode of heat transfer. At the same time, Rax values depending on the height changed significantly. According to the numerical results plotted the graphs of changes of relative axial
velocity and excessive temperature and relative Rayleigh number Ra* = Rax/Ra0 . There was
some critical height xcr - the beginning of the transient zone. Also determined the dependence
*KP = f (Rao).
Further results of numerical modeling using RANS approach are shown. The graphs of the relative change in axial velocity using different turbulence models are plotted. It is shown that, k-s model reproduce the kinematics of the flow better than the other on the main zone of the jet, and on the initial section - laminar model.
Keywords: natural convection, flat heat source, flow regimes, the critical value of the Rayleigh number, turbulence models.
Reference list
1. Baturin V.V., Elterman V.M. Aeration of industrial buildings. 2nd ed.. and ext. - M.: Gosstroiizdat, 1963. - 317 p.
2. Kutateladze S.S. Heat transfer and flow resistance: A Reference Guide. - M.: Energoatomisdat, 1990. - 367 p.
3. Popov I.A. Hydrodynamics and heat transfer of external and internal natural convection vertical flows with intensification. Enhancement of heat transfer: a monograph / Ed. Y.F. Gortyshov. - Kazan: Center for Innovative Technologies, 2007. - 326 p.
4. Martynenko O.G., Khramtsov P.P. Free-Convective Heat Transfer // Free. Heat Transf. With Many Photogr. Flows Heat Exch. Berlin / Heidelberg: Springer-Verlag, 2005. - 516 p.
5. Corcione M. Heat transfer correlations for free convection from upward-facing horizontal rectangular surfaces // WSEAS Trans. Heat Mass Transf., 2007, Vol. 2, № 3. - P. 48-60.
6. Wang H. Basic formulas and data on heat transfer for engineers / Directory. - M.: Atomizdat, 1979. - 216 p.
7. ANSYS FLUENT 6.3 Documentation / 12.3.8 Convective Heat and Mass Transfer Modeling. URL: https://www.sharcnet.ca/Software/Fluent6/html/ug/node476.htm (reference date: 02.27.2016).
8. Romanov S.V., Posokhin V.N., Ziganshin A.M. The effect of the method of heat input to the results of numerical solution of the problem of free convection on volumetric heat sources // Izvestija vysshih uchebnyh zavedenij. Problemy energetiki, 2013, № 3-4. - P. 48-53.
