Коэффициент теплопередачи вода алюминий

Расчет коэффициентов теплоотдачи

Интенсивность теплоотдачи зависит от динамического вида течения, определяющего структуру пограничного слоя у поверхности теплообмена, который в свою очередь зависит от скорости потока. Увеличение скорости потока ведет к уменьшению пограничного слоя, повышает турбулентность и приводит к увеличению интенсивности теплоотдачи.

Теплоотдача так же зависит от характеристик теплоносителя. Высокая теплопроводность уменьшает термическое сопротивление пограничного слоя и увеличивает теплоотдачу.

Снижение вязкости жидкости уменьшает пограничный слой, что так же благоприятно влияет на теплообмен между поверхностью и потоком теплоносителя.

Уменьшение пограничного слоя происходит так же в случае повышения кинематической вязкости или увеличения плотности рабочей среды, что так же повышает теплоотдачу.

Так же интенсивность теплоотдачи зависит от теплоемкости жидкости. При повышении теплоемкости повышается и теплоотдача, поскольку жидкость с большей теплоемкостью способна переносить большее количество теплоты.

Дополнительными факторами, влияющими на теплоотдачу, являются форма поверхности теплоотдачи, химические реакции и фазовые переходы в теплоносителе.

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

Расчет коэффициента теплоотдачи плоской стенки

Вычислить коэффициент теплоотдачи плоской поверхности можно с помощью уравнения подобия:

Nu_l = 0,66×Re_l 0,5 ×Pr 0,33 ; при ламинарном пограничном слое

Nu_l = 0,037×Re_l 0,8 ×Pr 0,43 ; при турбулентном пограничном слое

Re_l — число Рейнольдса, Pr — число Прандтля.

Исходные данные:

L — размер поверхности в направлении потока, миллиметрах;

w — скорость потока, метрах в секунду;

μ — динамическая вязкость теплоносителя, в паскаль×секунда;

ρ — плотность теплоносителя, в килограммах / метр 3 ;

λ — коэффициент теплопроводности теплоносителя, в ваттах / метр×°C×сек;

C_p — удельная теплоемкость теплоносителя, в джоулях / килограмм×°C.

КОЭФФИЦИЕНТ ТЕПЛООТДАЧИ ПЛОСКОЙ СТЕНКИ

Размер поверхности L, мм

Скорость потока, w, м/c

Динамическая вязкость, μ, Па*с

Плотность теплоносителя, ρ, кг/м 3

Теплопроводность, λ, Вт/(м* 0 C×сек)

Удельная теплоемкость, С_p, Дж/(кг* 0 C)

Источник

Коэффициент теплоотдачи поверхность — воздух

В статье рассмотрен расчет мощности теплового потока от горизонтальных и вертикальных плоских поверхностей тела, помещенного в «безразмерное» воздушное пространство при принудительной и естественной конвекции с учетом радиационной составляющей теплоотдачи.

Зная коэффициент теплоотдачи на поверхности (α), разделяющей твердое тело и окружающее это тело воздушное пространство, очень просто определить мощность теплового потока (Q) по известной разности температур (Δt).

Q=α*A*Δt, Вт – мощность теплового потока от или к поверхности тела.

• α=α_к+α_р, Вт/(м 2 *К) – суммарный коэффициент теплоотдачи на границе воздух – поверхность тела
  • α_к=?, Вт/(м 2 *К) – коэффициент конвективной теплоотдачи
  • α_р=ε*5,67*10 -8 *((t_п+273,15) 4 — (t_в+273,15) 4 )/(t_п-t_в)), Вт/(м 2 *К) – коэффициент радиационной теплоотдачи (теплоотдачи излучением), ε – степень черноты поверхности
• А, м 2 – площадь поверхности
• Δt=|t_п-t_в|, К – разность температур поверхности и воздушной среды
  • t_п, °C – температура поверхности
  • t_в, °C – температура воздуха

Основная сложность расчета заключается в определении коэффициента конвективной теплоотдачи (α_к)! Автоматизировать в первую очередь решение этой трудоемкой задачи поможет Excel.

Нестабильность процесса естественной конвекции у поверхностей различной формы и расположения в пространстве породила большое разнообразие эмпирических формул для вычисления коэффициента конвективной теплоотдачи (α_к). Неизбежные погрешности экспериментальных данных привели к тому, что результаты вычислений для одних и тех же поверхностей и условий по формулам разных авторов отличаются друг от друга на 20% и более.

После тщательного детального ознакомления с материалами современных западных изданий по теплообмену (список литературы – в конце статьи) были выбраны формулы, рекомендованные к применению большинством авторов, для использования в представленной далее программе в Excel.

Схемы теплообмена:

На представленных ниже рисунках показаны 8 вариантов схем, для которых программа может выполнить вычисления.

Розовый цвет пластин свидетельствует о том, что они горячее окружающего воздуха. Голубой цвет – пластины холоднее воздуха.

На схемах 1а и 1б воздух принудительно движется (вентилятор, ветер) вдоль поверхности пластины независимо от её ориентации в пространстве. На всех остальных схемах окружающий воздух находится в спокойном состоянии (помещение, полный штиль), а положение пластин сориентировано в пространстве.

Расчет в Excel:

Формулы алгоритма программы:

t₀=(t_в+t_п)/2

l₀=L – для схем 1а и 1б

l₀=(B*L)/(2*(B+L)) – для схем 2а, 2б, 3а, 3б, 4а, 4б

Для определения теплофизических параметров воздуха при определяющей температуре (t₀) в диапазоне -70°C … +1200°C использованы формулы из предыдущей статьи на сайте.

Re=w*l₀/ν

Gr=g*β*|t_п— t_в|*l₀ 3 /ν 2

Ra=Gr*Pr

α_к=Nu*λ/l₀

α_р=ε*0,00000005670367*((t_п+273,15) 4 — (t_в+273,15) 4 )/(t_п-t_в)) – при t_в *) α_р=0 – при t_в>t_п

α=α_к+α_р

q=α*(t_п-t_в)

Q=q*B*L

*) Нагрев поверхностей Солнцем или иными источниками теплового излучения программой игнорируется.

Вычисление теплофизических параметров воздуха и числа Нуссельта, как видно из вышеприведенных формул, являются ключевыми и самыми трудоемкими при определении конвективного коэффициента теплоотдачи.

Тестирование программы проводилось на примерах из книг, представленных в конце статьи. Отклонения результатов в основном не выходили за пределы ±5%.

Замечание:

В отечественной теплотехнической литературе для решения рассмотренных задач широко используются формулы второй половины прошлого века М.А. Михеева и В.П. Исаченко, которые в современной западной литературе не упоминаются. Беглый сравнительный анализ результатов расчетов по формулам разных авторов дал противоречивые и неоднозначные ответы. Если при принудительной конвекции результаты фактически идентичны, то при естественной конвекции отличаются порой на 30% и более, но иногда почти совпадают…

Литература:

1. John H. Lienhard IV, John H. Lienhard V, A Heat Transfer Textbook (Fifth Edition), 2019.
2. Frank Kreith, Raj M. Manglik, Mark S. Bohn, Principles of heat transfer (Seventh Edition), 2011.
3. Adrian Bejan, Convection Heat Transfer (Fourth Edition), 2013.
4. Michel Favre-Marinet, Sedat Tardu, Convective Heat Transfer, 2009.
5. Harlan H. Bengtson, Convection Heat Transfer Coefficient Estimation, 2010.
6. Rajendra Karwa, Heat and Mass Transfer, 2017.
7. Stuart W. Churchill, Humbert H. S. Chu, Correlating equations for laminar and turbulent free convection from a vertical plate, International Journal of Heat and Mass Transfer, Volume 18, Issue 11, November 1975.
8. http://people.csail.mit.edu/jaffer/SimRoof/Convection/
9. И. И. Кирвель, М. М. Бражников, Е. Н. Зацепин ЭНЕРГОСБЕРЕЖЕНИЕ В ПРОЦЕССАХ ТЕПЛООБМЕНА, 2007.

Прошу уважающих труд автора скачать файл с программой после подписки на анонсы статей!

P. S. (01.11.2020)

Дополнение по естественной конвекции у вертикальной поверхности:

Если построить графики по вышеприведенным формулам Черчилля и Чу для числа Нуссельта при естественной конвекции у вертикальной изотермической поверхности (схемы 2а и 2б), то можно увидеть, что при Ra=10 9 кривые не совпадают!

По этому поводу авторы формул Черчилль и Чу дают примерно следующее пояснение: «уравнение, основанное на исследованиях Черчилля и Усаги Nu=(0,825+0,387*Ra 1/6 /(1+(0,492/Pr) 9/16 ) 8/27 ) 2 дает хорошие результаты для средней теплопередачи при свободной конвекции у изотермической вертикальной пластины во всем диапазоне значений Ra и Pr от 0 до ∞, даже если оно не работает для обозначения дискретного перехода от ламинарного к турбулентному потоку». Линхарды в [1] отмечают, что рассматриваемое уравнение чуть менее точно для ламинарных условий при Ra 9 и рекомендуют в этом диапазоне использовать первое уравнение тех же авторов Nu=0,68+0,67*Ra ¼ /(1+(0,492/Pr) 9/16 ) 4/9 . Хотя, судя по графикам, в диапазоне Ra 7 для воздуха обе функции чрезвычайно близки друг к другу.

Еще один нюанс, который встретился только у Линхардов в [1]: «свойства флюида следует оценивать при t₀=(t_в+t_п)/2 за одним исключением, если флюид – газ, то коэффициент объемного расширения β следует определять при t₀=t_в». Но сами авторы зависимостей Черчилль и Чу о таком условии ничего не пишут. По этому поводу в их статье [7], говорится, что «для больших температурных перепадов, когда физические свойства существенно различаются, Ид рекомендует оценивать физические свойства как средние значения температуры поверхности и объема, а Уайли дает более подробные теоретические указания для режима ламинарного пограничного слоя».

Максимальная относительная ошибка для Nu=(0,825+0,387*Ra 1/6 /(1+(0,492/Pr) 9/16 ) 8/27 ) 2 , если β=1/t_в вместо β=2/( t_в+t_п), составляет в процентах:

ε=(((t_в+t_п)/(2*t_в)) 1/3 -1)*100%, или

ε=((|(t_п— t_в)|/(2*t_в)+1) 1/3 -1)*100%

Как видно из графика при температуре среды — воздуха t_в=20°C=293,15K и при перепаде температур поверхности и воздуха Δt=|t_п— t_в| 90 °C расхождение результатов быстро нарастает.

Правы Линхарды или множество других авторов, рассчитывающих все свойства флюидов при одном значении определяющей температуры t₀=(t_в+t_п)/2? Однозначного ответа у меня нет.

(По материалам Обри Джаффера [8].)

Эмпирические уравнения для суммарного коэффициента теплоотдачи:

В инженерных расчетах для быстрого приближенного определения суммарного коэффициента теплоотдачи, учитывающего и конвекцию, и излучение на границе поверхность тела – среда, можно использовать более простые зависимости, приведенные в [9].

При расчете тепловых потерь через наружные поверхности тел, которые находятся в спокойном воздухе закрытых помещений, можно применить нижеприведенные формулы. Результаты вычислений по этим формулам достаточно близки к результатам более точных расчетов.

α=9,74+0,07*(t_п-t_в), Вт/(м2*°C) при t_п On-line калькуляторы для расчетов коэффициентов конвективной теплоотдачи от плоских, цилиндрических и сферических поверхностей:

Инструменты представлены Группой исследований теплопередачи (HTRG). Группа была создана в 2014 году преподавателями Лаборатории теплотехники и жидкостей факультета машиностроения инженерной школы Сан-Карлоса (EESC) Университета Сан-Паулу (USP) для проведения передовых, качественных фундаментальных и прикладных исследований по вопросам теплопередачи для многофазных и однофазных систем.

Точность результатов вычислений не проверял.

Источник