книги из ГПНТБ / Романков, П. Г. Гидромеханические процессы химической технологии
.pdfAwx, Аwy, Аwz. Их значения и характеризуют турбулентность по тока.
Различают изотропную и анизотропную турбулентность. При изотропной турбулентности пульсации одинаково вероятны во всех направлениях, т. е. для Аwx, Awy и Аw2 имеется одинаковое число положительных и отрицательных значений. Другими словами, сред
ние квадратичные составляющие пульсационной скорости (У (Ада)2 имеют одинаковые значения в направлении осей координат: Awx — = Awy = Awz. Кроме того, при изотропной турбулентности сред нее значение произведения двух различных составляющих пульса ционной скорости равно нулю. При анизотропной турбулентности пульсационные скорости различны во всех направлениях и не рав новероятны.
Турбулентность можно характеризовать [3] так называемым турбулентным расстоянием (путь смешения) /т, представляющим
Рис. 3-5. Схема ламинарного режима |
Рис. 3-6. Схема турбулентного ре |
потока. |
жима потока. |
собой максимальное расстояние между точками, движущимися в турбулентном потоке с одинаковой скоростью. Турбулентное рас стояние постоянно изменяется во всем объеме жидкости и служит масштабом турбулентности. Максимальное значение оно будет иметь в центре потока, где скорость будет наибольшей.
Турбулентное расстояние /т является гидродинамическим ана логом длины свободного пробега молекулы /м в кинетической тео рии газов. Поэтому такое свойство жидкости, как вязкость (внут реннее трение), будет зависеть в турбулентном потоке от турбу лентного расстояния.С увеличением интенсивности турбулентности турбулентное расстояние Іг оказывается значительно превышаю щим длину свободного пробега /м, что приводит к существенному увеличению касательного напряжения т в потоке. Таким образом, в одномерном турбулентном потоке касательное напряжение по аналогии с законом вязкостного трения Ньютона будет выражать ся зависимостью
т = ± (Им + Ит) |
(3-30) |
где цм — коэффициент внутреннего (молекулярного) трения; цт —! коэффициент турбулентного трения (турбулентная вязкость); dwldx — градиент скорости.
57
Турбулентная вязкость изменяется пропорционально градиенту скорости от нуля у стенок трубы до относительно больших вели чин в центре потока. При развитой турбулентности можно прене бречь значением рм.
Характеристика режима движения зависит от средней скорости ©ср потока, его плотности р и вязкости р, а также от диаметра трубы d. Эти величины входят в безразмерный комплекс — кри терий Рейнольдса Re = wdp/ц — см. стр. 35.
Переход от ламинарного режима движения к турбулентному происходит при определенном, так называемом критическом значе нии критерия Рейнольдса (ReKp).
Ламинарный режим при движении потока в прямой трубе на
блюдается при значениях Re ^ |
2300 |
(рис. 3-5). В трубах с глад |
|||||||
|
|
|
|
кими |
стенками |
неустойчи |
|||
|
|
|
|
вый ламинарный режим дви |
|||||
|
|
|
|
жения |
может |
существовать |
|||
|
|
|
|
и при значениях Re, превы |
|||||
|
|
|
|
шающих ReKp, |
причем даже |
||||
|
|
|
|
небольшие |
возмущения |
или |
|||
|
|
|
|
начальная |
турбулентность |
||||
|
|
|
|
во входящем потоке вызы |
|||||
|
|
|
|
вают переход к турбулент |
|||||
|
|
|
|
ному режиму (рис. 3-6). Раз |
|||||
|
|
|
|
витый турбулентный режим |
|||||
|
|
|
|
наступает |
при |
Re > |
ІО4. |
||
Рис. 3-7. |
Зависимость ReKp от отношения |
Критическое значение крите |
|||||||
djD в |
змеевиках |
(d — диаметр трубы, |
рия Рейнольдса |
характерно |
|||||
D — диаметр витка змеевика). |
|
для каждой группы процес |
|||||||
прямых |
трубах |
ReKp = 2320, |
|
сов. Если для движения в |
|||||
то для движения |
в |
змеевиках |
ReKp = f(d/D), как показано на рис. 3-7. Следует отметить, что, на пример, для процесса осаждения ReKp = 0,2, для перемешивания
ReKp = 50 и т. п.
РАСПРЕДЕЛЕНИЕ СКОРОСТЕЙ ПО СЕЧЕНИЮ ПОТОКА И ЗАКОНЫ СОПРОТИВЛЕНИЯ
Скорость движения жидкости в трубах различного сечения в процессах химической технологии значительно меньше скорости звука, поэтому жидкость можно считать несжимаемой.
Рассмотрим установившееся движение вязкой несжимаемой
жидкости |
в длинной |
цилиндрической трубе. В условиях л а м и |
н а р н о г о |
р е ж и м а |
жидкость поступает в трубу с равномерной |
скоростью, причем на стенках образуется пограничный слой пере менной толщины (см. стр. 110). Вследствие этого на протяжении так называемого начального участка трубы профиль скорости ме няется от сечения к сечению. На некотором расстоянии от входа толщина пограничного слоя становится равной радиусу трубы и в дальнейшем профиль скорости уже не меняется.
58
Длина начального участка LHa4 приближенно определяется тео ретическим расчетом по уравнению [4]:
^нач — О»03d Re |
(3-31) |
где d — диаметр трубы; Re — значение критерия Рейнольдса. Таким образом, при Re = 1000 основной участок располагается
на расстоянии 30 диаметров от входа в трубу.
Система уравнений движения жидкости на основном участке, состоящая из уравнения неразрывности и уравнения Навье— Стокса
Уравнение (3-35) можно переписать в виде
Рис. 3-8. К выводу уравне ния движения жидкости в цилиндрической трубе.
причем для установившегося |
движения др/дх = dpjdx = |
const. |
|
Здесь под р следует понимать динамическое давление. |
|
||
Из уравнений (3-33) и (3-34) следует, что |
|
||
w = |
f(y, |
г, т) |
(3-36) |
Р = |
Г |
( * , т ) |
|
Решение задачи облегчается при использовании цилиндрических координат (рис. 3-8)
у = г cos Ѳ
(3-38)
2 = Г sin Ѳ
ИЛИ
r = V y 2 + z 2
(3-39)
Ѳ = arc tg —
У
Используя эти соотношения, преобразуем уравнение Навье — Стокса:
dp ^ |
I d2w |
1 dw |
1 |
d2w \ |
(3-40) |
|
dx |
^ \ дг2 |
г ’ дг ' г2 |
<ЗѲ- ) |
|||
|
Так как движение симметрично относительно оси х, то dzw!dQ2 = О и уравнение (3-40) упростится:
1 |
dp |
d 2w |
, 1 |
dw |
(3-41) |
|
]Г |
dx |
dr2 |
+ 7 |
dr |
||
|
69
Проинтегрируем уравнение (3-41) |
в |
граничных |
условиях |
|||||||||
dw/dr = |
0при г — О, |
а затем |
при w = |
0 и г = |
R |
и получим |
рас |
|||||
пределение скоростей в сечении трубы: |
|
|
|
|
|
|
|
|||||
|
|
|
|
4ц |
|
|
|
|
|
|
(3-42) |
|
|
|
|
|
|
|
|
|
|
|
|
||
Профиль скоростей (3-42) |
в поперечном сечении круглой ци |
|||||||||||
линдрической трубы представляет параболоид вращения |
(рис. 3-9). |
|||||||||||
Максимальная |
скорость достигается при |
г = 0 |
на |
оси трубы: |
||||||||
|
|
|
|
|
|
|
_ Я2 (_ dp |
(3-43) |
||||
|
|
|
|
|
|
® ш к с - |
4(1 ( |
|
d x |
|
||
|
|
|
|
|
Величина |
градиента |
давления |
|||||
|
|
|
|
|
зависит |
от объемного |
расхода |
ѴСек |
||||
|
|
|
|
|
жидкости, протекающей через трубу: |
|||||||
|
|
|
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
|
Усек = 2я J" wrdr = |
|
|
|
|
|
||
|
|
|
|
|
|
о |
|
|
|
|
|
|
Рис. 3-9. Профиль |
скоростей |
|
dp |
я |
|
|
|
|
|
|||
в поперечном сечении |
круглой |
|
dx |
2ц |
|
|
|
|
|
|||
цилиндрической трубы. |
|
|
|
лЯі |
( |
dp |
|
|||||
|
|
|
|
|
|
|
|
(3-44) |
||||
|
|
|
|
|
|
|
|
8ц |
V |
dx |
||
|
|
|
|
|
|
|
|
|
||||
Средняя скорость движения жидкости при ламинарном режиме |
||||||||||||
определяется зависимостью [с учетом уравнения |
(3-43)]: |
|
|
|||||||||
|
|
... |
|
V сек |
dp ■Я2 |
|
|
|
|
|
|
|
|
|
|
d x |
|
|
|
|
|
|
(3-45) |
||
|
|
иср ■ |
лЯ2 |
8ц |
2 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|||
Определим силу, действующую со стороны жидкости на участок |
||||||||||||
трубы круглого |
поперечного |
сечения |
длиною |
I |
и |
радиусом R |
||||||
(рис. 3-10). Из уравнения количества |
|
|
|
|
|
|
|
|||||
движения эта сила будет равна произ |
|
|
|
|
|
|
|
|||||
ведению |
разности |
давлений |
( р і — р 2) |
|
|
|
|
|
|
|
||
на площадь сечения. |
|
цилиндр |
|
|
|
|
|
|
|
|||
Выделенный |
в |
жидкости |
|
|
|
|
|
|
|
|||
радиусом г на этом же участке трубы |
|
|
|
|
|
|
|
|||||
находится в равновесии под действием |
|
|
|
|
|
|
|
|||||
силы сопротивления и касательных на |
|
|
|
|
|
|
|
|||||
пряжений т, приложенных к его боко |
|
|
|
|
|
|
|
|||||
вой поверхности: |
|
|
|
|
|
|
|
|
|
|
||
|
т • 2ягі = (р2 — рі) ял2 |
|
|
|
|
|
|
|
|
|||
Отсюда |
|
|
|
|
|
Рис. 3-10. |
К |
выводу закона |
||||
|
|
Р1— Р2 Г |
|
|||||||||
|
|
(3-46) |
сопротивления |
при движении |
||||||||
|
|
|
I |
’ 2 |
||||||||
|
|
|
|
потока в круглой трубе. |
Из уравнения (3-46) следует, что касательное напряжение про порционально расстоянию от оси трубы и достигает максимального
60
значения у стенки |
Ар R_ |
|
||
|
(3-47) |
|||
|
/ |
2 |
||
Коэффициент трения f может быть определен как отношение |
||||
силы сопротивления ApnR2 к скоростному напору рwlp/2 |
и к пло |
|||
щади участка боковой поверхности трубы: |
|
|||
f- |
АpnR2 |
RAp |
(3-48) |
|
Pwlp 2nRl |
21 Pwlp |
|||
|
|
|||
Заменим R через диаметр трубы d: |
|
|
||
|
d |
Дp |
|
|
|
' - 7 - T ' |
Pwcp12 |
(3'49) |
|
Отсюда потеря давления на преодоление трения при движении |
||||
потока в круглой цилиндрической трубе |
|
|||
|
/ |
Pwtp |
(3-50) |
|
|
Ар = 4/- |
|
||
Уравнение (3-50) известно как |
уравнение Фаннинга |
(/ — коэф |
фициент Фаннинга), однако в отечественной литературе оно более
распространено как урав- |
л |
|
||||||
нение Дарси — Вейсбаха |
|
|
||||||
с заменой 4f на Я: |
|
|
|
|
||||
|
|
|
РЩ |
|
|
|
|
|
|
Ар = |
Я ■ |
|
ср |
(3-51) |
|
|
|
|
|
|
|
|
||||
где |
Я — коэффициент |
со |
|
|
||||
противления. |
в |
уравне |
|
|
||||
Подставив |
|
|
||||||
ние |
(3-51) |
значение |
|
|
||||
Ар/1 — wCp-8p/R2, |
соглас |
|
|
|||||
но уравнению (3-45) |
по |
|
|
|||||
лучим для Я выражение: |
|
|
||||||
Я = |
2d |
8|д®ср |
|
32ц |
Рис. 3-11. Зависимость коэффициента |
сопро |
||
|
Р®ср |
|
|
Р^ср |
|
|||
|
|
|
|
|
(3-52) |
тивления для гладких круглых труб |
при ла |
|
|
|
|
R = |
|
минарном режиме движения. |
|
||
Обозначив |
|
Д/2 |
и |
Сплошная кривая построена по уравнению (3-53), |
||||
wCpdp/p = 'R.e, |
получим за |
точки нанесены по опытным данным Гагена. |
||||||
|
|
кон сопротивления при ламинарном движении в круглой цилинд рической трубе с гладкой внутренней поверхностью:
Я = |
64 |
(3-53) |
|
Re |
|
Зависимость Я = /(Re) приведена |
на рис. 3-11. |
Ламинарный режим движения был впервые изучен более 100 лет тому назад Гагеном и Пуазейлем [2, 5].
61
Уравнение для определения потери напора на трение при ла минарном режиме, полученное преобразованием зависимостей (3-51) и (3-52), носит название уравнения Гагена — Пуазейля (см.
стр. 30):
ДрТр • |
32[ли)ср/ |
(3-54) |
|
Т* |
|||
|
|
Уравнение (3-54) справедливо при малых значениях критерия Рейнольдса (Re < ReKP) и особенно важно при исследовании дви жения жидкости в трубах малого дйаметра, а также в капиллярах и порах.
Знание распределения скорости и перепада давлений в потоке
при т у р б у л е н т н о м д в и ж е н и и |
для решения многих инже |
||||||||
|
|
нерных |
задач |
представляет |
|||||
|
|
еще большую важность. Одна |
|||||||
г |
— |
ко точного решения уравнения |
|||||||
Навье — Стокса для этого слу |
|||||||||
|
|||||||||
|
|
чая не существует. Имеющие |
|||||||
|
|
ся |
рекомендации основаны |
на |
|||||
|
|
экспериментальных |
данных |
и |
|||||
|
|
подкреплены |
теоретическими |
||||||
|
|
рассуждениями. |
|
режи |
|||||
|
|
|
При |
турбулентном |
|||||
|
|
ме (Re > |
ReIip) движения жид |
||||||
Рис. 3-12. Изменение усредненной по |
кости в |
трубе |
также |
следует |
|||||
времени скорости при турбулентном |
учитывать длину |
начального |
|||||||
движении потока в круглой трубе: |
участка. По данным Никурад- |
||||||||
1—ламинарный подслой; 2—переходная зона; |
|||||||||
3—область развитой турбулентности (турбу |
зе |
[6], І„ ач = (254-40) d*. |
|
||||||
|
лентное ядро). |
|
Для |
получения |
уравнения |
||||
|
|
распределения |
скорости в круг |
лой трубе при развитом турбулентном режиме можно разделить область движения на турбулентное ядро и ламинарный подслой вблизи стенки (рис. 3-12). В ламинарном подслое скорость жид кости мала, пульсации скорости практически отсутствуют, но вследствие прилипания жидкости к обтекаемым стенкам имеют место очень большие поперечные градиенты скорости, которые вы зывают значительные напряжения силы трения (в полном соот
ветствии с законом Ньютона т = ± ц - 4 ^ • В турбулентном ядре
вследствие большой извилистости и сложности траекторий частиц жидкости уравнения движения заменяют зависимостями между осредненными величинами и ищут их решение, используя пара метры, описывающие мгновенное состояние движения потока (в частности, осредненные уравнения количества движения приме няются для получения так называемых уравнений Рейнольдса, устанавливающих связь между турбулентными напряжениями в йотоке). Таким образом, в случае турбулентного режима движения
* По данным Кирстена [7], Z.H»ч = (60 -t- 100)rf.
62
закон распределения скорости может быть получен только на осно вании анализа экспериментальных данных. Менаду ламинарным подслоем и турбулентным ядром находится переходная зона, для которой одинаково важны и молекулярная вязкость и турбулент ность.
В ламинарном подслое распределение скоростей можно считать линейным:
где г — расстояние |
от оси трубы (в направлении, |
перпендикуляр |
|||
ном стенке); |
бл — толщина ламинарного подслоя |
(порядка 1 мм). |
|||
В турбулентном |
ядре |
рас |
|
||
пределение |
осредненных |
ско |
|
||
ростей в |
пределах |
изменения |
|
||
значений |
критерия |
Рейнольд |
|
||
са от ІО4 до ІО5 хорошо опи |
|
||||
сывается |
степенной зависи |
|
|||
мостью: |
|
|
|
|
|
W |
|
Г_\4п |
(3-56) |
|
|
®макс |
|
R } |
|
|
где п зависит от величины критерия Re и в данных пре делах может быть принято равным 7 (по эксперименталь ным данным).
Таким образом, прибли женно для турбулентного тече ния *
®макс 5
Рис. 3-13. Сравнение распределения скоростей при ламинарном (1) и тур булентном (2) режимах движения в круглой трубе.
Качественное сравнение распределения скоростей при ламинар ном и турбулентном режимах показано на рис. 3-13.
В гидромеханических процессах при решении задач, связанных с определением режима транспортирования жидкостей или газов в трубопроводах обычно пользуются графической зависимостью между отношением аУср/^макс и значением критерия Рейнольдса
(рис. 3-14).
Для определения коэффициента сопротивления Я при турбу лентном режиме движения в пределах изменения значений кри терия Re от 4 -ІО3 до ІО5 для гидравлически гладких труб можно пользоваться формулой Блазиуса:
Я = 0,316 |
(3-58) |
Re1-'4 |
|
* С помощью уравнений пограничного слоя может быть рассчитано, что
ZfiVp = 0,817Шмакс [8].
63
Более точная зависимость (для больших значений Re) между коэффициентом сопротивления Я и режимом движения может быть получена при использовании логарифмического закона распределе ния скоростей. При выводе логарифмического (универсального) профиля Re^- оо, так как пренебрегают молекулярной вязкостью ц по сравнению с турбулентной (хт (см. стр. 57).
Рис. 3-14. Зависимость между отношением а'ср/іИмакс и критерием Рейнольдса;
, - “’ср/и,макс= Г(а’м а к с ^ ); 2- й'ср/“’макс“ Г К р ^ )
Для значений критерия Рейнольдса Re > ІО5 коэффициент со противления можно рассчитать по формуле:
у = - = 2 1 g (R e К Л ) — 0.8 |
(3-59) |
Сравнение применимости формул (3-58) и (3-59) приведено на рис. 3-15.
Исследованиями Никурадзе [6], Шиллера [9] и других ученых
установлено, что |
коэффициент |
сопротивления |
Я |
в значительной |
||||
А |
|
|
степени зависит также |
и от |
||||
|
|
шероховатости труб: |
|
|||||
|
|
|
|
|||||
|
|
|
Я = f (Re, е) |
|
(3-60) |
|||
|
/ 2 |
|
где е — эффективная высота |
|||||
|
|
выступов на внутренней по |
||||||
|
|
|
верхности |
|
трубы. |
Обычно |
||
|
Г |
|
для характеристики шеро |
|||||
|
|
ховатости |
используют |
так |
||||
103 |
ш5 |
ю7 |
называемую |
относительную |
||||
шероховатость ejd |
или |
d/e |
||||||
|
Re |
|
||||||
Рис. 3-15. Значения |
коэффициента Я для |
(d — диаметр трубы). |
|
|||||
Если высота выступов е |
||||||||
гладких круглых труб при турбулентном |
||||||||
режиме: |
|
в трубе |
меньше |
толщины |
||||
/ —по формуле (3-58); 2—по формуле (3-59). |
|
ламинарного подслоя б, то |
||||||
|
|
|
шероховатость стенок не ока |
зывает влияния на величину коэффициента сопротивления Я при турбулентном режиме движения потока. Такие трубы носят на
звание |
гидравлически гладких. При |
большой высоте выступов |
(е > б) |
турбулентность потока увеличивается и сопротивление воз |
|
растает |
(при этом профиль скоростей |
изменяется), Наконец, при |
64
некотором предельном значении шероховатости сопротивление ста новится постоянным и такие трубы называют вполне шерохова тыми.
Так как толщина ламинарного подслоя 6 = /(Re), то труба может быть вполне шероховатой при одном расходе жидкости и гидравлически гладкой при другом.
Рис. 3-16. Значения относительной шероховатости для промышленных труб:
трубы\ 1, 3— клепаные стальные; 2, 4—бетонные; 5—чугун
ные; |
6— деревянные; 7— оцинкованные железные; 8— чугун |
ные, |
покрытые битумом; 9—из мягкой сварочной стали; |
|
/0-— цельнотянутые. |
При ламинарном режиме движения влияние шероховатости сте нок трубы на сопротивление очень незначительно и им обычно пренебрегают.
Относительная шероховатость промышленных труб приведена на рис. 3-16 [10].
Зависимость коэффициента сопротивления К от режима тече ния (Re) и относительной шероховатости e/d (или d/e) графически представлена на рис. 3-17 и 3-18. Для некруглых труб d — dg.
3 Зак. 840 |
65 |
л
4,0
2,0
500
700
1000
1500
2500
4 000
Рис. 3-18. Зависимость коэффициента трения Я в технических трубах от критерия Re при различной шероховатости стенок.