Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

35782

.pdf
Скачиваний:
13
Добавлен:
21.03.2016
Размер:
7.48 Mб
Скачать

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

Аналитические методы не позволяют описать отрывные течения, поэтому, формулы строятся на обработке экспериментальных данных. График α = f (φ) представлен на рис.8.6. Для количественных расчетов можно рекомендовать формулы [6,7]:

Re f ,d = 5 ÷103;

 

0,38

 

 

Prf

 

0,25

 

 

 

 

 

 

 

Nu f ,d = 0,5 Re f ,d Prf

 

 

 

 

 

;

 

Pr

 

 

 

 

w

 

Re f ,d =103 ÷2 105;

 

 

 

 

 

 

0,6

0,38

 

 

Prf

 

0,25

 

 

 

 

 

 

Nu f ,d = 0,25Re f ,d

Prf

 

 

 

.

Pr

 

 

 

 

 

w

 

где индексы f и w говорят о том, что соответствующие физические параметры, входящие в числа Nu, Pr, Re, рассчитываются при температурах жидкости (f) и стенки (w).

8.4. Полуэмпирическая теория турбулентности

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

t = t + t ,

(8.25)

 

 

 

где t - средняя температура; t- пульсационное изменение температуры. Подставим в дифференциальное уравнение энергии значение

температуры (8.25) и используем зависимость

acρ = λ ,

тогда уравнение

энергии после преобразований примет вид

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

t

 

 

t

 

 

t

 

 

 

 

 

v

 

(8.26)

ρcvx

 

+ ρcvy

 

=

 

 

λ

 

 

 

 

+

 

(ρcvyt

) ρvx vy

 

x .

 

 

 

 

 

 

 

 

 

x

 

 

 

 

 

 

 

 

 

 

 

y

 

 

 

 

y

 

 

 

y y

 

y

 

 

 

 

 

 

Проделаем те же преобразования с уравнением Навье-Стокса, т.е. примем ряд допущений. Первое допущение: представим произведение ρcvyt в

форме, аналогичной первому члену в правой части уравнения (8.26):

 

 

 

 

 

 

 

 

 

 

 

 

t

(8.27)

ρcvyt ′ = Aq

 

 

,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

y

 

где Aq - коэффициент турбулентной теплопроводности. Тогда уравнение энергии (8.26) примет вид

81

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

t

 

t

 

t

 

ρcvx

+ ρcvy

=

(λ + Aq )

+ (μ + Aσ )

vx

.

(8.28)

 

 

 

 

 

 

 

 

 

x

 

y

 

y

y

 

 

 

 

 

 

 

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Допущение №2 сохраняет силу и выражается равенством:

σT =τW .

Далее идёт ряд допущений которые позволяют в итоге получить значения скоростей при турбулентном движении и гидравлического сопротивления. Коэффициент теплоотдачи рассчитывают по формуле (8.19)

St = Nu /(Re Pr) = 0.5R,

которая, как уже указывалось, удовлетворительно работает для случая, когда Pr 1.

Глава 9. Теплообмен в разрежённых газах

9.1. Описание физических процессов

а) Критерий Кнудсена. При обычном рассмотрении переноса теплоты в газах структура газа считается сплошной и поэтому не требуется привлечения представлений о молекулярном строении газа. Поток и явления переноса теплоты при таких условиях непрерывности среды могут быть адекватно выражены через критерии Рейнольдса, Нуссельта и Прандтля. Однако при малых абсолютных давлениях газ частично теряет характерные свойства непрерывности, и появляются явления, которые могут быть объяснены только в том случае, когда будут приняты во внимание представления о молекулярном строении газа.

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

Для характеристики разреженности газа вводят критерий Кнудсена

Kn =

Λ

,

(9.1)

 

L

 

 

где Λ - средняя длина свободного пробега молекул, L - характерный размер тела.

Если Kn>>1, то газ считается разреженным, при Kn<<1 говорят о сплошной среде (более детальная классификация режимов будет приведена в дальнейшем).

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

82

Таблица 9.1

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

б) Коэффициент скольжения. Введем коэффициент

f =

mvτ п mvτ o

,

(9.2)

 

 

mvτ п

 

Рис. 9.1. Скачки скорости и температуры

который характеризует полноту обмена тангенциальным количеством движения. Здесь m - масса молекулы, vτ о и vτ п - тангенциальные

составляющие скорости до и после отражения (рис.9.2). Если соударение упругое, то f = 0, при абсолютном неупругом взаимодействии f=1 (при неупругом соударении молекула поглощается стенкой, т.е. временно задерживается, и в дальнейшем может иметь любую тангенциальную составляющую).

Для газового потока (большого числа молекул) коэффициент f принимает значение

f =

v u

,

(9.3)

v

 

 

 

83

Рис. 9.2. Здесь vτ о и vτ п тангенциальные составляющие скорости потока до и после взаимодействия со стенкой.

В большинстве случаев молекулы, исходящие от стенки, могут разлетаться в разные стороны беспорядочно и u – близко к 0, т.е. f≈1. В табл.9.2 приведены значения f для некоторых систем. В большинстве случаев в первом приближении можно считать f=1.

Таблица 9.2

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

 

dv

y

 

 

(9.4)

v |x=0

vw = ξ

 

 

= v,

 

 

 

 

dy

 

 

 

 

 

x =0

 

 

где ξ =

2 f

Λ - коэффициент скольжения; v |x=0

-

скорость газа у

f

 

 

 

 

 

 

 

поверхности

стенки; vw - скорость поверхности;

dv

y

 

- градиент

 

 

 

 

 

 

 

 

 

 

 

 

 

скорости у поверхности.

 

dy x =0

 

 

 

 

 

 

в) Коэффициент аккомодации. Эффект температурного скачка был впервые обнаружен в 1898г. Смолуховским. В дальнейшем Кнудсен предложил ввести коэффициент аккомодации, который характеризует долю истинного энергообмена по сравнению с максимально возможным:

84

a =

Eпад Еотр

,

(9.5)

 

 

Епад Еw

 

где Eпад , Eотр - энергия падающей и отраженной молекулы, Ew

- энергия

молекулы при температуре стенки.

Для одноатомного идеального газа точно, а для других газов с хорошим приближением равенство (9.5) можно записать в виде

a =

tпад tотр

.

(9.6)

 

 

tпад tw

 

Значения а для некоторых систем приведены в табл.9.3.

Таблица 9.3

Из молекулярно-кинетической теории газа Максвелл показал, что температурный скачок на поверхности описывается выражением

T = 2 a

2k

 

Λ

 

T

 

 

(9.7)

 

 

 

,

 

 

 

a

k +1 Pr y

 

 

 

x =0

 

 

где k = Cp - постоянная адиабаты.

Cv

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

9.2. Режимы течения газа

При исследовании теплообмена необходимо знать характер протекающих процессов и доминирующих явлений. Это можно сделать с помощью классификации режимов течения газа, которая проводится по численным значениям следующих критериев: критерия Маха M, Re, Kn. Укажем взаимосвязь между ними:

Kn =1,26 k

M

.

(9.8)

 

 

Re

 

85

Для одноатомного и двухатомного газов формула (9.8) имеет соответственно вид

Kn =1,66

M

;

Kn =1,48

M

.

Re

 

 

 

 

Re

В отдельных случаях в критерий Кнудсена в качестве определяющего размера вводят толщину пограничного слоя δ

Knδ = δΛ = ΛL δL = Kn δL .

При обтекании пластины ламинарным потоком

δ

=

4,64

.

(9.9)

L

 

 

Re

 

Классификация режимов, предложенная Тзяном, приведена в табл.9.4.

Таблица 9.4

Другая классификация режимов предложена Девиеном: Kn <10-3 - непрерывная среда;

103 < Kn < 0,25 - режим скольжения;

0,25 < Kn <10 - переходный режим от течения со скольжением к свободно-молекулярному;

Kn >10 - свободно-молекулярный режим.

При Kn >10 межмолекулярными столкновениями можно пренебречь. Область аэродинамики, которая занимается этой проблемой, называется супераэродинамикой, или молекулярной аэродинамикой.

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

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

86

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

Запишем выражение для теплового потока q в условиях непрерывного течения (индекс "0"):

 

T

 

(9.10)

q = −λ

 

 

=α0 (Tf Tw0 ).

 

 

x x=0

 

 

Если газовый поток разрежен, а критерий Re=idem (соблюдается гидродинамическое подобие), то возникает температурный скачок, который вызывает дополнительное тепловое сопротивление. Запишем выражение, аналогичное (9.10):

 

T

 

(9.11)

q = −λ

 

 

=α(Tf Tw ).

 

 

x x =0

 

 

где α - учитывает как тепловое сопротивление пограничного слоя, так и температурный скачок.

При одинаковых тепловых потоках в пограничном слое у поверхности градиенты температуры одинаковы, т.е. одинаковы углы наклона касательной к температурной кривой T=T(x) у стенки (рис.9.3); тогда для температур справедливо соотношение

Tw = Tw0 T.

(9.12)

Рис. 9.3. Температурное поле в пограничном слое при отсутствии ( Т=0) и наличии скачка температуры на поверхности

Приравняем правые части (9.10), (9.11) с учетом (9.12)

α0 (Tf Tw0 ) =α(Tf Tw0 ) +α T.

Разделим обе части на α0 (Tf Tw0 ) :

1 =

 

α

+

 

α

 

 

 

T

.

α

 

α

 

T

 

 

 

0

 

0

 

f

T

 

 

 

 

 

 

 

 

w0

87

Согласно выражению (9.7)

T = β

Λ

T

;

β′ =

2 a

 

2k

 

.

 

 

 

 

 

 

 

a

k +

1

 

Pr

x x =0

 

 

 

Тогда

1 = α + α β′Λ.

α0 λ Pr

Здесь использовано соотношение (9.10)

 

T

=

α

0

(T

 

T

).

 

 

 

 

 

 

 

 

 

x x =0

 

λ

 

f

w0

 

Окончательно в критериальном виде

Nu =

 

 

 

1

 

 

.

(9.13)

 

1

 

+ β

Kn

 

 

 

Nu

0

Pr

 

 

 

 

 

 

 

 

 

 

 

 

 

Формула (9.13) учитывает температурный скачок на одной поверхности, что достигнуто за счет рассмотрения условий на границе раздела. Эта схема была предложена Кавенау и группой американских ученых в конце 50-х годов. Сравнение с экспериментальными данными показывает, что формула (9.13) в целом правильно отражает влияние температурного скачка, однако при больших M наблюдается сильное расхождение (т.к. скачок скорости не учитывается) с экспериментом.

9.4. Эффективные значения коэффициентов переноса (метод Р.С. Прасолова)

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

Рассмотрим перенос теплоты в плоской газовой прослойке, ограниченной поверхностями с температурами T1 и T2 (рис.9.4).

Рис. 9.4. Ход температуры в плоской газовой прослойке

88

Плотность теплового потока в прослойке можно выразить через истинную теплопроводность газа, увеличив прослойку на толщину ϕ1 + ϕ2 :

 

 

q = λ0

 

 

T1 T2

 

.

 

 

(9.14)

 

 

δ +ϕ +ϕ

2

 

 

 

1

 

 

 

 

 

 

 

 

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

q, запишем

 

 

 

T1 T2

 

 

 

 

 

 

 

 

 

 

 

q = λэф

.

 

 

 

 

 

(9.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

δ

 

 

 

 

 

 

Сравнив выражения (9.14) и (9.15), найдем

 

 

 

λэф = λ0

δ

 

 

 

.

 

(9.16)

 

 

δ +ϕ +ϕ

2

 

 

1

 

 

 

 

 

Для нахождения ϕ1 и ϕ2 воспользуемся равенством (9.7), где

производную dT

 

x=0 для линейного

 

распределения

температуры можно

 

 

dy

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

записать в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dT

 

 

 

=

 

T .

 

 

 

 

 

 

dy

 

x=0

 

 

 

ϕ

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

ϕ = 2 a a k2+k1 PrΛ .

Используя последнее соотношение для нахождения ϕ1 и ϕ2 , получим выражение для λэф :

λэф

 

2k

 

 

 

2 a

 

2 a

 

 

Kn 1

(9.17)

 

 

 

 

 

 

1

 

 

 

2

 

 

 

λ

= 1+ k +

1

+

a

 

Pr

.

 

a

2

 

 

0

 

 

 

 

 

1

 

 

 

 

0

 

 

Здесь λ0 и Pr0 берутся из таблиц при нормальном давлении (рис.9.5).

Рис. 9.5. Зависимость λэф от критерия Кнудсена

λ0

89

Если рассматривается теплообмен одиночной пластины, когда имеет место только один температурный скачок, то

λэф

=

1

+

2k

 

 

2 a

 

Kn 1.

(9.18)

λ

k +1

a

 

 

 

 

 

Pr

 

 

0

 

 

 

 

 

 

 

 

0

 

 

Аналогичное выражение можно получить для эффективного значения динамической вязкости с учетом скачка скорости. С одной

стороны, касательное напряжение

τ = μ0

v

 

 

, с другой стороны,

δ +ε

1

+ε

2

 

 

 

 

 

τ = μ δv .

Приравняв правые части последних выражений, с учетом (9.4) получим

μэф

 

 

2 f

1

 

2 f

2

1

(9.19)

 

 

 

 

+

 

 

 

 

 

 

 

 

 

 

μ

= 1+

f1

 

f2

 

Kn .

 

 

 

 

 

 

 

 

При теплообмене одиночной пластины

μэф

=

1

+ 2

f

Kn 1.

(9.20)

 

 

μ

 

 

f

 

 

 

 

 

 

9.5. Расчёт коэффициента теплоотдачи по Р.С. Прасолову

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

Коэффициент восстановления при обтекании пластины при ламинарном режиме течения определяется по формуле r = Pr . Найдем значение Pr с учетом скачка скорости и температуры v и T:

 

 

 

 

 

 

 

 

1+

2 a

 

2k

 

 

Kn

 

 

 

 

 

 

ν

 

νcρ

 

μc

a

 

k +1

Pr

 

 

μ

c

 

Pr =

=

=

=

 

 

 

 

 

 

=

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

0

 

 

a

 

λ

λ

 

1+

2

f

Kn

 

 

λ0

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1+

2

+ a

 

2k

 

 

 

Kn

 

 

 

a

k +

1

Pr

Pr0 .

(9.21)

 

 

 

 

 

 

 

 

 

 

 

 

1

+

2

f

Kn

 

 

 

 

 

 

 

 

 

 

 

 

 

 

f

 

 

 

 

 

 

 

 

 

Сравним это с литературными данными. Если Kn>>1, то

 

 

 

 

 

 

 

 

 

 

 

1

 

 

 

 

 

 

r =

 

 

(2 a)2kf

 

 

2

 

 

 

 

 

 

Pr =

 

 

.

 

 

 

 

 

a(k

+1) Pr (2

f )

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

 

 

Молекулярно-кинетическая теория для данного режима дает

соотношение

r =

2k

 

 

, которое при

a=0,9

и

f=1 находится в хорошем

k +

1

 

 

 

 

 

 

 

 

 

 

 

соответствии с последним выражением.

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

90

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]