Потери напора при турбулентном течении жидкости
Рис. 4.3. Пульсация скорости в турбулентном потоке
При турбулентном режиме движения жидкости в трубах эпюра распределения скоростей имеет вид, показанный на рис. 4.4.
Рис. 4.4. Характер линий тока
в турбулентном потоке
В тонком пристенном слое толщиной δ жидкость течет в ламинарном режиме, а остальные слои текут в турбулентном режиме, и называются турбулентным ядром (рис. 4.5). Таким образом, строго говоря, турбулентного движения в чистом виде не существует. Оно сопровождается ламинарным движением у стенок, хотя слой δ с ламинарным режимом весьма мал по сравнению с турбулентным ядром. Кроме того, турбулентное течение всегда является неустановившимся, так как значения скоростей и давлений, а также траектории частиц, изменяются по времени. Однако его можно рассматривать как установившееся течение при условии, что осреднённые по времени значения скоростей и давлений, а также полный расход потока не изменяются со временем. Такое течение встречается на практике достаточно часто.
Рис. 4.5. Модель турбулентного режима движения жидкости
Основной расчетной формулой для потерь напора при турбулентном течении жидкости в круглых трубах является эмпирическая формула, называемая формулой Вейсбаха-Дарси и имеющая следующий вид:
Различие заключается лишь в значениях коэффициента гидравлического трения λ. Этот коэффициент зависит от числа Рейнольдса Re и от безразмерного геометрического фактора - относительной шероховатости Δ/d (или Δ/r0, где r0 - радиус трубы).
Первые систематические опыты для выявления влияния различных параметров на величину λ были проведены Никурадзе[3] под руководством Людвига Прандтля в 20-х г.г. XX века в Германии.
Эти опыты проводились в латунных трубах, гладких, что достигалось шлифовкой и с искусственной однородной шероховатостью, которая создавалась наклеиванием зерен песка определённого размера на внутреннюю поверхность труб. В трубах с полученной таким образом определённой шероховатостью при разных расходах измерялась потеря напора и вычислялся коэффициент λ, значения которого наносились на график в функции числа Рейнольдса. Результаты опытов Никурадзе представлены графически на рис. 4.6. На этом графике по горизонтальной оси отложены величины lgRe, а по вертикальной оси — lg(l00 λ). Кривые построены по данным опытов с трубами относительной шероховатости от ε=∆/d= 0,001 (самая нижняя кривая) до ε=0,033 (самая верхняя кривая).
Анализируя представленный график, можно сделать следующие выводы: cуществуют четыре различные области.
Рис. 4.6. График Никурадзе
1. Область ламинарного режима (I). В области ламинарного режима (т.е. при Re < 2300, чему соответствует lg Re < 3,36) опытные точки, независимо от шероховатости стенок, уложились на одну прямую линию I. Следовательно, здесь λ зависит только от числа Рейнольдса и не зависит от шероховатости, т.е. λ =f (Re).
Остальные участки кривых (II, III, IV) относятся к турбулентному движению.
В области перехода от ламинарного движения к турбулентному Re = 2000-4000 (3,3< lgRe< 3,6) наблюдается большой разброс опытных точек и кривая между I и II на рис. 4.6 проведена условно.
2. Область гидравлически гладких труб (II). В этой области опытные точки для труб с различной шероховатостью располагаются в некотором диапазоне чисел Re на одной прямой II, отрываясь от неё в сторону возрастания коэффициента λ тем раньше, чем больше шероховатость стенок. Таким образом, при некоторых условиях шероховатость не оказывает влияния на потери напора также и при турбулентном движении, т.е. и здесь λ =f (Re).
3. Область смешанного трения (III). Здесь каждая кривая относится к определённому значению относительной шероховатости и величина также меняется с изменением числа Рейнольдса, т.е. коэффициент гидравлического сопротивления зависит как от числа Re, так и от ε (λ =f(Re,ε)).
4. Область «вполне шероховатых труб» (IV), При увеличении числа Re кривые области III переходят в линии, параллельные оси lg Re, т,е. коэффициент λ в этой области не зависит от числа Re и определяется только относительной шероховатостью. Полуэмпиричекая теория турбулентности позволяет предложить выражение для коэффициента λ, исходя из распределения скорости в живых сечениях потока.
Предложенная полуэмпирическая теория не отражает особенностей сопротивления в области смешанного трения.
Опыты Никурадзе проводились в трубах с одной искусственной шероховатостью. Трубы же, применяемые на практике, имеют шероховатость неоднородную и неравномерную. Поэтому долгое время оставалось неясным, насколько правильны будут выводы, полученные Никурадзе на трубах с искусственной шероховатостью, в применении к обычным промышленным трубам с естественной шероховатостью и каковы численные значения шероховатости для подобных труб. Выяснению этих вопросов был посвящен ряд проведенных экспериментальных исследований (работы Колбрука, И.А. Исаева, Г.А. Мурина, Ф.А. Шевелева).
Наибольший интерес представляют опыты Г.А. Мурина по исследованию гидравлических сопротивлений в обычных промышленных стальных трубах, законченные в 1948 г. Результаты этих опытов представлены на графике, изображенном рис. 4.7, показывающем изменение коэффициента λ в зависимости от числа Рейнольдса для стальных труб.
Рис. 4.7. Диаграмма Г.А. Мурина
Подтвердив основные закономерности, установленные Никурадзе, эти опыты показали, что для труб с естественной шероховатостью коэффициент λ в переходной области имеет всегда большие значения, чем в случае вполне шероховатых труб (а не меньше, как у Никурадзе), Поэтому кривые на диаграмме Мурина не имеют впадины, характерной для кривых Никурадзе.
Рис. 4.6.
График Никурадзе
Итак, потери напора, определяемые по формуле Вейсбаха-Дарси, можно определить, зная коэффициент гидравлического сопротивления, который определяется в зависимости от числа Рейнольдса Re и от эквивалентной относительной шероховатости. Для удобства сводные данные по определению λ представлены в таблице 4.1. Пользоваться приведенными в табл. 4.1 формулами для определения коэффициента λ не всегда удобно. Для облегчения расчетов можно воспользоваться номограммой Колбрука-Уайта (рис.4.8), при помощи которой по известным Re и Δэ/ d весьма просто определяется λ.
Таблица 4.1
Таблица для определения коэффициента гидравлического трения
Рис. 4.8. Номограмма Колбрука-Уайта для определения коэффициента гидравлического трения