Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
Теория электронагрева (курс лекций).doc
Скачиваний:
147
Добавлен:
05.11.2018
Размер:
7.82 Mб
Скачать

Теплоотдача при вынужденном продольном омывании плоской поверхности

Такая задача может возникнуть в случае принудительного охлаждения (с помощью воздуха, нагнетаемого вентилятором) какой-либо плоской поверхности.

Считаем, что жидкость несжимаема. За пределами гидродинамического и температурного пограничных слоев скорость и температура постоянны и равны 0 и Т0.

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

Интегральные уравнения пограничного слоя

Для теплового пограничного слоя уравнение энергии имеет вид:

х  у  а

Поскольку qy = -  Т/у и   Т2/у2 = - qy/у  уравнение энергии имеет вид

ср(х у) - qy/у (1)

Проинтегрируем это уравнение от у = 0 до у = . Заметим, что за пределами пограничного слоя производные, входящие в уравнение (1), равны 0. Поэтому увеличение верхнего предела от k до  не дает изменение интеграла.

Интегрирование правой части уравнения (1) дает

- dy = qc, (плотность теплового потока на стенке)

здесь учтено, что (qy)у= = -  (Т/у) у= = 0.

В результате интегрирования получим

х(T0 - T)dy = - – (2)

  • интегральное уравнение теплового потока для теплового пограничного слоя

Здесь k – толщина теплового пограничного слоя;

T0 – температура невозмущенного потока.

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

k = k(х)

Уравнение движения в проекциях на ось Ох для рассмотрим здесь течения в приближении пограничного слоя:

Wx + Wу   (А)

Учитывая, что касательная сила трения S =  xy) и   , представим уравнение (А) в

динамический коэф.вязк. кинематич.коэф.вязкости

виде (Wx + Wу) = Sy (3)

Видим, что уравнение энергии (1) и уравнение движения (3) аналогичны. Поэтому, выполняя аналогичные преобразования над уравнением (3) получим интегральное уравнение импульсов для гидроэлектродинамического пограничного слоя

(4)

Здесь Sс – касательное напряжение трения при у = 0, то есть на поверхности стенки; - толщина гидродинамического пограничного слоя.

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

Теплоотдача при ламинарном пограничном слое

Для расчета теплоотдачи при ламинарном слое используем интегральное уравнение теплового потока для теплового пограничного слоя:

Чтобы рассчитать теплоотдачу, необходимо знать распределение скорости в слое. В случае ламинарного пограничного слоя распределение скорости в нем близко к параболе (кубическая парабола):

х = а + bу + су2 + dу3 (a)

Уравнение распределения скорости должно удовлетворять граничным условиям.

При у = 0 х = 0 (условие «прилипания»), полагаем также, что (2х/у2) = 0.

При у = (на внешней границе слоя) 0 = х и (х/у) = 0.

Уравнение (а) будет удовлетворять этим требованиям, если

а = 0; с = 0; b = 3/20/; d = - 1/20/3.

Распределение скорости при этом примет вид:

(б)

При распределении скорости согласно (б) из интегрального уравнения импульсов

-

можно получить, что толщина гидродинамического пограничного слоя определяется выражением:

или в безразмерной форме

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

Пограничные условия для температур:

При у = 0  (здесь и ниже принято  = Т – Тс) и (у)у=0 = const и  (2у2) = 0, если учесть, что в жидкости, непосредственно прилегающей к плоской стенке, теплота передается по у только теплопроводностью.

При у = k (на внешней границе теплового слоя)  = 0 = const и (у)у=k = 0.

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

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

0  1,5 (у/k) – 0,5 (у/k)3 (А)

После подстановки в интегральное уравнение теплового потока и упрощений получили:

,

где Reх = 0х/; Рr = .

Тогда – то есть отношение k/ не зависит от х, а зависит только от числа Рr.

Для капельных жидкостей Рr  1,  k.

Для газов Рr = 0,6 – 1 (для воздуха Рr  0,7). При этом k.

Определим коэффициент теплоотдачи из уравнения   у=0 зная, что (это следует из производной), опуская знак - получим:

(а)

Коэффициент теплоотдачи обратно пропорционален толщине пограничного слоя (чем меньше толщина k, тем больше  лучше теплоотдачи).

Уравнение (а) можно привести к безразмерному виду (умножив левую и правую части на х/ и подставляя значение ):

(б)

Здесь Nux = x/ = / x/ = NuеX; Rex = 0x/ = 0/ x/ = ReеX; Pr = /a; – длина пластины вдоль потока.

Формула (б) справедлива только при условии, что температура поверхности пластины постоянна, физические параметры жидкости, не зависят от температуры и в начале пластины нет необогреваемого участка.

Коэффициент теплоотдачи капельной жидкости зависит от рода жидкости, ее температуры, направления теплового потока. Особенно существенное влияние оказывает зависимость вязкости от температуры. При охлаждении жидкости, вязкость увеличивается, течение замедляется. Поэтому для жидкостей вводится поправка (Prж/ Prс)0,25, которая учитывает влияние на теплообмен изменения вязкости. Для (Prж/ Prс)0,25  1.

Индекс «с» - означает, что Prс вычислен при температуре стенки; «ж» - при температуре жидкости вдали от стенки.

Переход ламинарного течения в турбулентное

1 – ламинарный пограничный слой; 2 – переходная область; 3 – турбулентный пограничный слой; 4 – вязкий (ламинарный) подслой.

Переход ламинарного течения в турбулентное происходит на некотором участке. Течение на этом участке имеет нестабильный характер и называется переходным.

Законы теплообмена при ламинарном и турбулентном режимах различны, поэтому определение их границ имеет большое значение.

О режиме течения судят по критическим значениям числа Re:

Reкр1 = 0хкр1/ и Reкр2 = 0хкр2/,

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

Reкр = 104 - 4106

На переход влияют такие характеристики внешнего источника, как степень турбулентности, масштаб турбулентности. При ускорении потока переход затягивается, при замедлении - наступает при меньших значениях х (или Reх).

Кроме параметров внешнего потока на переход из ламинарного течения в турбулентное влияют параметры, связанные с омываемым теплом. Значения Reкр1 и Reкр2 зависят:

  1. от интенсивности теплообмена;

  2. от волнистости, шероховатости поверхности;

  3. обтекаемости передней кромки пластины;

  4. вибрации тела.

На рис. представлена зависимость критических чисел Рейнольдса от степени турбулентности набегающего потока

,

где 2x, 2y, 2z – средние во времени квадраты трех составляющих пульсации скорости; 0 – скорость внешнего потока.

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

Увеличение Tu приводит к уменьшению Reкр.

На практике сечение перехода можно определить по изменению осредненной во времени х (у). При турбулентном течении х резко увеличивается вблизи стенки, на удалении от нее х (у) становится более выровненной. Выравнивание объясняется турбулентным переносом кол-ва движения.

Теплоотдача при турбулентном пограничном слое

Перенос теплоты и количества движения поперек турбулентного пограничного слоя можно описать уравнениями:

q = - (  т     срq,

S = (  т    s,

где S – касательная сила трения; т, т – коэффициенты турбулентного переноса теплоты и кол-во движения; q = т/ср ,s = т/ - кинематические коэффициенты турбулентного переноса теплоты и кол-ва движения.

Эти уравнения можно переписать в виде:

q = - (1 + ) (1)

S = (1 + s/) (2)

Здесь Prт = s/q; Pr = /а = //с

Величину Prт – называется турбулентным числом Прандтля. С учетом (1) и (2) дифференциальное уравнение энергии и движения для турбулентного пограничного слоя примут вид

 а 1 + ), (3)

  1 + s/). (4)

Если Pr = 1 (а = ) и Prт = 1, то уравнения (3) и (4) идентичны. В этом случае при идентичных граничных условиях поля температуры  и скорости х будут подобны.

Для того, чтобы проинтегрировать уравнение (3) и (4), необходимо иметь сведения о коэффициентах турбулентного переноса теплоты s и q.

Опыты показывают сложность движения в турбулентном слое:

А

Б

II

п I

Рис. Турбулентный пограничный слой

А – внешняя область; Б – пристенная область (I – вязкий подслой (ламинарный); II – промежуточный слой).

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

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

Аналогично вязкому подслою непосредственно у стенки наблюдается тепловой подслой.

Он характеризуется преобладанием теплоты теплопроводностью над турбулентным переносом.

При Pr = 1 толщина вязкого подслоя п и теплового kп совпадают.

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

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

На основе эмпирических данных и вычислений был получен безразмерный комплекс – число Стантона

St

При Pr = 1 и Prт = 1 – получаем аналогию переноса теплоты и кол-ва движения. На основе опытных данных получена формула

Nuжх = 0,0296Reжх0,8Prж0,43(Prж/Prс)0,25

За определяющую принята температура жидкости вдали от тела Т0 (за исключением Prс, выбираемого по температуре Тс – стенки). Определяющим размером является координата х, отсчитываемая от начала участка теплообмена.

Рис. а – ламинарное течение; б – смешанное (переходное); в - турбулентное течение

Посмотрим как изменяется коэффициент теплоотдачи вдоль пластины.

Если вся пластина занята турбулентным слоем (в случае высокой стенки турбулентности набегающего потока, плохой обтекаемости пластины и т.п.), то коэффициент  изменяется в соответствии с кривой 1. при наличии на передней части пластины ламинарного пограничного слоя коэффициент изменяется по более сложному закону. В этом случае теплоотдачу необходимо рассчитывать отдельно для участков с различными режимами течения.