Движение газа в каналах простой формы. Учебное пособие
.pdf
Ядро постоянных скоростей
Пограничный |
слой |
1 |
2 |
Рис. 2.1. Скорость в прямой трубе:
1 – начальный участок; 2 – основной участок
с законом Ньютона возникает касательное напряжение трения, которое равно произведению динамической вязкости на градиент скорости по нормали к поверхности:
τ=µ |
∂c |
(н/м2). |
(2.1) |
∂n |
На элементарном участке канала длиной dl возникает сила трения dRw = τ dl U (н), где U – смоченный периметр канала, равный πd (м) для круглой трубы. Сила сопротивления движения преодолевается за счет уменьшения механической энергии газа.
Так как за пределами пограничного слоя скорости остаются равномерными, сила вязкости отсутствует. Поток разделяется на вязкий пограничный слой и невязкое ядро потока. В поперечном сечении, поскольку поток не отклоняется от прямолинейного направления, поле статических давлений остается равномерным. Поэтому в пограничном слое полное давление меньше полного
давления в ядре потока: pпс* = p +0,5c2 p +0,5C2 = pя* . Здесь и далее при сопоставлении параметров потока скорость в погранич-
номслоебудемобозначатьстрочнымибуквами,авядре–пропис- ными.
Преодоление сопротивления трения происходит за счет снижения кинетической энергии в пограничном слое. Так как по мере продвижения вглубь канала потеря механической энергии нарас-
l
тает (Nw =(c Rw )l =U ∫cτdl (Вт)), масса газа с пониженной кине-
0
11
тической энергией увеличивается, т. е. по мере продвижения вглубь канала пограничный слой занимает все большую часть поперечного сечения. Толщина пограничного слоя – это условная величина.Допустимоприниматьееравнойрасстояниюотповерхности стенки до точки в потоке, в которой значение скорости меньше скорости основной части (ядра) потока на 1 %. Более строгим считается определение с учетом интегральной величины – толщины вытеснения [3]:
δ*
∫(C −c)dn
δ* = |
0 |
|
. |
(2.2) |
|
C |
|||
|
|
|
|
Помереудаленияотначалаканалапроисходитувеличениетолщины пограничного слоя. Так как пограничный слой в середине канала смыкается, постепенно невязкое ядро потока перестает существовать. Под начальным участком подразумевают участок трубы, в котором равномерный профиль скоростей на входе уступает место профилю развитого течения, где невязкое ядро отсутствует.
Второе название начального участка – разгонный участок – связано с тем, что в середине канала максимальная скорость возрастает из-за увеличения толщины вытеснения. При ламинарном течении на основном участке скорость в середине трубы вдвое больше равномерной скорости на входе в канал.
На основном участке силы вязкости проявляются по всему сечению. Деление потока на ядро и пограничный слой здесь неправомерно. В турбомашинах длина каналов обычно меньше длины разгонного участка. В этих каналах, как правило, можно выделить ядро потока и пограничные слои на стенках. Поэтому простейшей схемой каналов ТК следует считать именно начальный участок трубы. Аналогично движению газа на начальном участке трубы:
–вканалахтурбомашиннеизбежнапотерямеханическойэнергии газа из-за трения о стенки проточной части;
–пограничный слой оказывает влияние на поле скоростей невязкого ядра, «загромождая» проходные сечения каналов тем
12
сильнее, чем дальше от входа в канал находится рассматриваемое поперечное сечение;
– в особо узких каналах проточной части, характерных, например,дляцентробежныхкомпрессороввысокогодавления,течение может принять развитый характер.
Ламинарный и турбулентный характер течения. При разных ско-
ростях движение газа в прямой трубе свидетельствует о существовании двух разных режимов течения:
–при увеличении скорости от нуля до предельного значения газовые частицы движутся строго поступательно, отдельные слои газа – с разной скоростью, но не перемешиваются между собой. Это ламинарный (слоистый) режим течения;
–при дальнейшем повышении скорости движение газовых частиц становится нестационарным и хаотическим. Значение продольной скорости меняется во времени относительно среднего значения. Кроме того, появляются поперечные составляющие скорости,вызывающиевзаимныеперемещениягазовыхчастициз одного слоя в другой. Это турбулентный режим течения.
Хаотический турбулентный характер течения удобно рассматривать в виде наложения на вектор средней продольной скорости турбулентных пульсаций (рис. 2.2).
Впродольномнаправлениидажевневязкомтурбулентномядре частицы соседних слоев постоянно смещаются друг относительно друга, что вызывает потерю трения. Измерения показывают, что в проточнойчаститурбомашинэтипотерималы.Наличиепродольных пульсаций скорости – нехарактерное отличие турбулентного течения от ламинарного.
Существенное различие вносит наличие поперечных турбулентных пульсаций. При этом в пограничном слое или в развитом потоке без невязкого ядра ламинарные и турбулентные режимы существенно отличаются.
wrn wrs wrs
w
wrn
Рис. 2.2. Схема турбулентного потока как результат сложения осредненнойскоростиспродольной и поперечной пульсацией
13
При турбулентном течении градиенты скорости ∂c/∂n у стенок значительно больше, чем при ламинарном потоке. Это связано с тем, что частицы быстро движущихся слоев, попадая в более медленные слои, отдают его частицам часть своей кинетической энергии. И наоборот, частицы из более близких к стенкам канала слоев, попадая в быстро движущиеся слои, приобретают часть кинетической энергии его частиц.
Ламинарный или турбулентный характер течения определяется соотношением между силами инерции и вязкости, действующими на газовые частицы. Это соотношение устанавливает критерий подобия газовых потоков, носящий имя английского ученого Рейнольдса. Для анализа развитого течения в каналах критерий Рейнольдса имеет вид
Re= |
cdг |
= |
cdг |
ρ , |
(2.3) |
ν |
|
||||
|
µ |
|
|
||
где с – средняя расходная скорость в сечении канала; dг = 4f/U – гидравлический диаметр (f – площадь сечения; U – смоченный периметрканала);ν–кинематическаявязкость;μ–динамическая вязкость; ρ – плотность среды.
Так как при турбулентном движении имеет место беспорядочное перемешивание частиц, касательное напряжение намного
выше, чем при ламинарном движении: τ=µ dydc +ρl2 dydc dydc , где
l – длина пути перемешивания. В результате наличия высоких касательных напряжений и турбулентного перемешивания слоев газанаразгонномучасткетурбулентныйпограничныйслойстановится тоньше. Различаются и профили скорости на основном участке течения при ламинарном и турбулентном режимах (рис. 2.3). При ламинарном режиме стабилизированный профиль скорости устанавливается по параболическому закону, при турбулентном режиме – приближенно по логарифмическому или степенному закону.
Расчет потери напора в основном участке канала. На основном участке прямолинейных каналов профили скорости по длине не
14
w
r0
y
l ≥ lнач
1
2
Рис. 2.3. Профили скорости на основном участке трубы: 1 – ламинарный режим; 2 – турбулентный режим
меняются. Сопротивление трения преодолевается за счет снижения статического давления*.
Это делает экспериментальное определение закона сопротивленияоченьпростым.Посоображениямразмерностикасательное напряжение пропорционально динамическому давлению:
|
c2 |
|
|
τ=λ1ρ |
ср |
τ , |
(2.4) |
2 |
где τ – касательное напряжение на поверхности стенки канала;
cср = 1 ∫cdf – средняя расходная скорость (индекс «ср» обычно f f
опускают); λ1 – безразмерный коэффициент пропорциональности.
* Это справедливо при отсутствии илипренебрежимо малом влиянии сжимаемости. При движении газа по очень длинной трубе (магистральный газопровод) давление и плотность заметно уменьшаются. При этом скорость газа увеличивается вплоть до скорости звука в гипотетическом случае.
15
Так как на основном участке каналов эпюры скоростей по длине не меняются, то везде наблюдается одно и то же касательное напряжение и один и тот же безразмерный коэффициент трения. При этом на участке канала длиной l сила сопротивления равна произведению касательного напряжения на смоченную поверхность, которая равна произведению длины участка на периметр
c2
поперечного сечения канала: Rтр = τUl =λ1ρ 2ср Ul .
Сила трения преодолевается силой давления – произведением потери давления ∆p на поперечное сечение канала F: Rтр = R дав,
c2
т. е. Rтр = Rдавλ1ρ 2ср Ul = ∆pF . Отсюда на участке канала длиной l потеря давления
∆p =λ ρ |
cср2 |
|
l |
. |
(2.5) |
2 |
|
F U |
|||
1 |
|
|
|
Отношение поперечного сечения к его периметру называется гидравлическим радиусом: rг = F/U. Так как для канала круглого сечения гидравлический диаметр равен его диаметру, удобнее пользоваться учетверенной величиной гидравлического радиуса, котораяназываетсягидравлическимдиаметром:dг = 4F/U.Сучетом этого потеря давления и коэффициент потерь основного участка канала имеют вид
∆p =λρ |
cср2 |
|
l |
|
; |
|
(2.6) |
|||||
2 |
|
|
|
|
||||||||
|
|
|
|
dг |
|
|
|
|||||
ζтр = |
∆p |
|
|
=λ |
l |
, |
(2.7) |
|||||
|
c2 |
|
|
d |
г |
|||||||
|
|
|
|
|
|
|
|
|
||||
|
ρ |
ср |
|
|
|
|
|
|
|
|
||
|
2 |
|
|
|
|
|
|
|
|
|
|
|
где λ = 4λ1 – эмпирический коэффициент трения.
Итак, коэффициент потерь трубы пропорционален безразмерному коэффициенту трения, длине трубы и обратно пропорционален гидравлическому диаметру. Потеря давления пропорциональна коэффициенту потерь и динамическому давлению ρc2/2.
16
Зависимость коэффициента трения от числа Рейнольдса и относи-
тельной шероховатости. В соответствии с теорией подобия безразмерный коэффициент трения зависит от формы канала F и кри-
териев подобия: λ= f(F, k, M, Re).
Форму прямого канала определяет форма его поперечного сечения. Опыты показывают, что входящее в уравнения (2.6) соотношение длины и гидравлического диаметра канала, непосредственно отражая влияние формы на потери, не влияет на коэффициент трения λ. Он практически одинаков у каналов с круглым, прямоугольным и другим поперечным сечением. На коэффициент трения при определенных условиях оказывает влияние
относительная шероховатость поверхностей |
|
|
ш = |
kш |
– отноше- |
|
k |
||||||
|
|
|||||
|
|
|
|
dг |
||
ниевысотынеровностейкгидравлическомудиаметру.Дляпрямой трубы относительная шероховатость – это геометрический параметр, по которому можно определить форму трубы.
При рассмотрении каналов ограниченной длины критерии сжимаемости k, M не влияют на газодинамические параметры по определению. В них потеря давления не столь велика, чтобы заметно влиять на снижение плотности газа по длине канала.
Определяющим критерием подобия является критерий дина-
мического подобия – число Рейнольдса |
|
||||
Re= |
cdг |
=ρ |
cdг |
. |
(2.8) |
ν |
|
||||
|
µ |
|
|
||
Критерий Рейнольдса характеризует соотношение между силами инерции и вязкости в потоке. Первые пропорциональны произведению характерного линейного размера на характерную скорость движения, вторые – кинематической вязкости. При развитом течении в прямом канале характерный линейный размер – гид- равлическийдиаметрканала,характернаяскорость–средняярас- ходная. Для ТК вместо кинематической вязкости удобнее пользо-
ваться динамической вязкостью µ= ρν .
17
Экспериментальное определение коэффициента трения не вызывает трудностей. Для этого достаточно измерить уровень падения давления на участке трубы с развитым течением длиной l с известным гидравлическим диаметром dг. Полученная плотность газа (жидкости) и измеренный массовый расход позволяют определить коэффициент потерь по формулам (2.6) и (2.7) и соответствующий данному эксперименту критерий Рейнольдса по формуле (2.8).
На рис. 2.4 представлена экспериментальная зависимость коэффициента трения от числа Рейнольдса и от относительной шероховатости по Никурадзе.
Следует обратить внимание на то, что обычно характерные для течения в турбомашинах значения критерия Рейнольдса большие, они соответствуют турбулентному режиму. Влияние критерия Рейнольдса наблюдается при очень значительном его изменении. Поэтому для наглядности графиков по оси абсцисс используется логарифмическая шкала.
λ
0,100
0,090
0,080
0,070
0,060
0,050
0,045 Ламинарноетечение
0,040
0,035
0,030
0,025
0,020
0,018
0,016 5 103
|
|
|
Турбулентное |
|
|
||
|
|
|
течение |
|
|
|
|
= |
|
|
λ |
|
0,3164 |
|
|
λ |
|
|
= |
|
|
|
|
|
|
|
|
|
|
||
0,3164 |
|
|
Re |
|
|
||
Re |
|
|
|
|
0,25 |
|
|
|
|
|
|
|
|
|
|
0,25 |
|
|
|
|
|
|
|
2 |
3 |
4 5 |
104 |
|
2 |
3 |
4 5 |
|
r |
|
=15 |
|
|
|
|
Kш |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
30,5 |
|
|
|
|
60 |
|
|
|
|
||
|
|
126 |
|
|
|
|
|
|
252 |
|
|
|
|
|
|
507 |
|
|
|
|
105 |
2 |
3 |
4 5 |
106 |
||
|
|
|
|
|
|
Re |
Рис. 2.4. Зависимость коэффициента трения λ для круглой трубы от числа Рейнольдса и относительной шероховатости по Никурадзе [3]
18
Экспериментальные данные показывают, что при ламинарном режиметечения(Re 2300)коэффициенттрениянезависитотшероховатости труб и связан только со значением числа Рейнольдса:
λ = 64/Re. |
(2.9) |
Отсутствие влияния шероховатости на сопротивление связано с тем, что неровности поверхностей закрыты слоем практически неподвижного газа (жидкости). Сопротивление определяется трением движущегося газа о неподвижный пристеночный слой.
Так как поперечные турбулентные пульсации способствуют возрастанию градиента скорости у стенок канала, в зоне перехода от ламинарного режима к турбулентному режиму сопротивление увеличивается (см. рис. 2.3). При этом касательное напряжение
τ=µ' ∂∂nc ρ возрастает. Легко показать, однако, что безразмерный
коэффициент трения λ обратно пропорционален числу Рейнольдса. Поэтому при дальнейшем увеличении Re во время движения газа в канале с гидравлически гладкой поверхностью коэффициент трения уменьшается.
Так как для поперечных турбулентных пульсаций поверхность непроницаема, даже при турбулентном течении в непосредственной близости от поверхности они не наблюдаются. У поверхности существует ламинарный подслой. Поверхность гидравлически гладкая при условии, что этот подслой покрывает неровности поверхности. Сопротивление определяется скольжением турбулентного потока по ламинарному подслою, и шероховатость не оказывает влияния на коэффициент трения. Для турбулентного течения в канале с гидравлически гладкими стенками коэффициент трения обычно оценивают по формуле Блазиуса
λ = 0,31641 . |
(2.10) |
Re 4 |
|
Формула (2.10) точна в ограниченном диапазоне критерия Re, но существуют более универсальные соотношения [3].
С ростом числа Рейнольдса толщина ламинарного подслоя уменьшается. Любая гидравлически гладкая поверхность стано-
19
вится шероховатой, когда вершины неровностей поверхности начинают выступать за пределы ламинарного подслоя. При этом вершины неровностей обтекаются турбулентным потоком. При обтекании неровностей наблюдается отрыв потока и образование вихрей, кинетическая энергия которых безвозвратно преобразуется в тепло. Сопротивление определяется потерей механической энергии в результате вихреобразования и перестает зависеть от числа Рейнольдса. Область течения, в которой коэффициент сопротивления не зависит от числа Рейнольдса, называется обла-
стью автомодельности по числу Рейнольдса.
Коэффициент трения λ шероховатых каналов в области автомодельности зависит от относительной шероховатости:
kш = kdш , (2.11)
г
где kш – средняя высота неровностей поверхности.
Для расчета коэффициента трения канала с шероховатой поверхностью рекомендуют использовать формулу, аппроксимирующую экспериментальные данные Никурадзе [3]:
|
|
|
|
|
|
|
|
2 |
|
|
|
|
1 |
|
|
|
|
||
λ = |
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
. |
(2.12) |
|
|
1 |
|
|
|
|||||
|
|
|
|
|
|
||||
|
2,0 lg |
|
|
|
|
+ 1,74 |
|
||
|
|
|
|
|
|||||
|
|
|
2kш |
|
|
|
|||
|
|
|
|
|
|||||
Чем больше относительная шероховатость, тем при меньшем числе Рейнольдса гидравлически гладкая поверхность становится шероховатой.Процесспереходаотгидравлическигладкойкшероховатой поверхности соответствует значению Рейнольдса
R |
= |
ckш |
ρ = 100. |
(2.13) |
|
||||
0 |
µ |
|
|
|
У критерия Рейнольдса в форме (2.13) за характерный линейный размер принимается средняя высота неровностей. Граница между гидравлически гладким и шероховатым режимами обтека-
20
