Научный журнал
Фундаментальные исследования
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,074

ОБОБЩЕНИЕ ФОРМУЛЫ КОЛБРУКА – УАЙТА НА ТЕЧЕНИЯ ЖИДКОСТИ В ТРУБЕ С ПРОИЗВОЛЬНОЙ ПЕСОЧНОЙ ШЕРОХОВАТОСТЬЮ СТЕНКИ

Кондратьев А.С. 1 Ньа Т.Л. 1 Швыдько П.П. 1
1 Московский политехнический университет
Проведенный анализ инженерных методов расчета коэффициентов гидравлического сопротивления в круглой трубе и профиля скорости показал, что все рассмотренные методы расчета течения при переходном режиме направлены на сопряжение расчетных зависимостей для гидравлически гладкой и предельно шероховатой поверхности. В получаемых расчетных зависимостях не учитывалось наблюдаемое в опытах минимальное значение коэффициента гидравлического сопротивления, величина которого зависела от числа Рейнольдса и относительной шероховатости поверхности. В данной работе на основе анализа известных опытных данных Никурадзе и формулы Колбрука – Уайта предложена новая эмпирическая зависимость, которая достаточно хорошо согласуется с опытными данными, в том числе в области минимальных значений коэффициента гидравлического сопротивления. Получена простая зависимость, аппроксимирующая результаты численных расчетов, позволяющая определить профиль скорости при произвольной степени шероховатости стенки трубы.
турбулентный режим течения
гладкая стенка
песочная шероховатость
гидравлическое сопротивление
профиль скорости
1. Альтшуль А.Д. Гидравлические сопротивления. – М.: Недра, 1982. – 224 с.
2. Байков В.Н. Универсальное распределение скоростей в водных потоках при различных режимах гидравлического сопротивления // Вест. МГСУ. – 2009. – № 1. – С. 19–184.
3. Брянская Ю.В. Уточнение кинематических характеристик турбулентного течения // Инженерно-строительный жур. – 2013. – № 6. – С. 31–38.
4. Брянская Ю.В. Перемежаемость течения при переходном режиме гидравлического // Вест. МГСУ. – 2013. – № 1. – С. 177–184.
5. Кондратьев А.С., Ньа Л.Т., Швыдько П.П. Инженерный метод расчета коэффициента гидравлического сопротивления и профиля скорости при произвольной песочной шероховатости стенки трубы // Гидравлика. Сетевой журн. – 2016. – № 12. – С. 41–49.
6. Кондратьев А.С., Овсянников В.М., Олофинский Е.П. и др. Транспортирование водоугольных суспензий: гидродинамика и температурный режим. – М.: Недра, 1988. – 213 с.
7. Лобанов И.Е. Теория гидравлического сопротивления в прямых круглых трубах с шероховатыми стенками // Отраслевые аспекты технических наук Науч. – практ. журн. – 2012. – № 4(16). – С. 4–13.
8. Храмцов И.В., Бульбович Р.В., Павлоградский В.В. Влияине шероховатости на профиль скорости в трубах квадратного сечения // Приволжский науч. вест. – 2013. – № 6 (22). – С. 5–12.
9. Шлихтинг Г. Теория пограничного слоя. – М.: Физматгиз, 1974. – 712 с.
10. Янышев Д.С. Применение функции Ламберта в теории турбулентного трения // Электр. журн. Труды МАИ. – 2015. – Вып. 50. – 11 c; URL: www. mai.ru/science/trudy.

Необходимость повышения точности расчета коэффициентов гидравлического сопротивления при течении жидкостей в трубах с различной степенью песочной шероховатости внутренней поверхности стенки трубы требует дальнейшего обобщения опытных данных и, на их основе, совершенствования инженерных методов расчета. Основой для анализа являются опытные данные Никурадзе по движению ньютоновской жидкости в гладких и шероховатых трубах, приведенные, например, в работе [9]. На основании анализа этих опытных данных, с привлечением гипотез Прандтля, Кармана и других исследователей о физической природе турбулентного движения жидкости, получены достаточно простые аналитические выражения, позволяющие определить коэффициенты гидравлического сопротивления и профиля скорости в двух предельных случаях: гидравлически гладкой стенки (kond01.wmf) и режима течения с полным проявлением шероховатости kond02.wmf. Здесь λ – коэффициент гидравлического сопротивления; υτ – динамическая скорость; ν – кинематическая вязкость жидкости; ks – величина песочной (эквивалентной) шероховатости стенки; R – радиус трубы; kond03.wmf – число Рейнольдса, определенное по средней скорости Um и диаметру (2R).

В значительно меньшей степени исследованы режимы течения жидкости при переходном режиме

kond04.wmf.

Проанализируем работы, выполненные в последние годы. В работах [2, 3], используя, формулы Прандтля и Кармана для коэффициентов гидравлического сопротивления гидравлически гладкой и предельно шероховатой поверхности

kond05.wmf (1)

kond06.wmf (2)

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

kond07.wmf, (3)

kond08.wmf (31)

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

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

В работе [7] расчетная схема включает прилегающий к стенке вязкий подслой, переходящий в турбулентную зону смешения с постоянной длиной смешения. В области выше элементов шероховатости реализуется турбулентный режим течения с длиной пути смешения, зависящий от расстояния от поверхности. Заметим, что такая схематизация течения жидкости вдоль шероховатой поверхности ранее использовалась в работе [6] для расчета течения не только в ньютоновской жидкости, но и для неньютоновских жидкостей. Автор сравнивал результаты расчетов не с опытными данными Никурадзе, а с расчетом по формуле Колбрука – Уайта, которая, как известно [10], применима не к песочной, а к технической шероховатости. Кроме того, сравнение расчетных данных проводилось при больших числах Рейнольдса от области минимального значения коэффициента гидравлического сопротивления до максимального, соответствующего режиму с предельным проявлением шероховатости. Например, выпал из рассмотрения диапазон чисел Рейнольдса kond09.wmf для частиц с максимальной песочной шероховатостью R/ks = 15, а для частиц с минимальной шероховатостью R/ks = 507 не рассмотрен диапазон чисел Рейнольдса kond10.wmf.

В работе [4] автор предлагает рассчитывать коэффициент гидравлического сопротивления при переходном λt, гладком λ и квадратичном λs режимах сопротивления по формуле, выражающей принцип суперпозиции:

kond11.wmf, (4)

где γ – коэффициент перемежаемости, равный отношению времени существования турбулентного течения у стенки к общему времени наблюдения. Именно такая зависимость ранее использовалась в [6], где величина γ физически толковалась как доля поверхности, течение над которой соответствует режиму течения над гидравлически гладкой поверхностью, а kond12.wmf течению над поверхностью с предельным проявлением шероховатости. Затем, не имея достаточных оснований, выражение (2) записывается в виде [4]

kond13.wmf. (41)

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

В большинстве из рассмотренных работ [2–4, 7] использовались логарифмические профили скорости для гидравлически гладкой и предельно шероховатой стенки трубы и соответствующие им зависимости коэффициентов гидравлических сопротивлений, а затем они обобщались на случай переходной режим. В работе [8], в рамках традиционной схемы течения над шероховатой поверхностью [9], выражение для скорости с произвольной степенью песочной шероховатости поверхности записывается в виде

kond14.wmf, (5)

где величина B зависит от режима обтекания поверхности.

При переходном режиме используется известная зависимость

kond15.wmf;

при 0,025 < ks/ δ< 5 и 3 < υτks/ν < 70. (6)

Также как в работе [7] закон сопротивления при переходном режиме задается формулой Колбрука – Уайта

kond16.wmf, (7)

что, как отмечалось выше, физически некорректно [9]. Поэтому расчеты течения в круглой трубе приведут к тем же ограниченным результатам, которые были получены в работе [7].

Следует отметить, что при использовании выражения (7) и подобной, наиболее простой в части процедуры вычислений, формуле Альтшуля [1]

kond17.wmf (8)

величина λ плавно уменьшается с ростом числа Рейнольдса, то есть полностью игнорируется наличие минимума величины λ, зависящее от обеих величин ks/(2R) и Re. Игнорирование этого опытного факта представляется весьма существенным, в связи с чем все конкретные рекомендации по расчету профиля скорости в этой зоне представляются недостаточно обоснованными.

В работе [10] автор показывает, что формулы Прандтля и Колбрука – Уайта, соответственно, могут быть выражены в виде функции Ламберта, зависящей только от величин ks/(2R) и Re. Использование функции Ламберта в чисто вычислительном аспекте не имеет никакого преимущества перед использованием классических выражений в области их применимости.

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

kond18a.wmf (9)

где f – поправочная функция, зависящая от безразмерного комплекса (ksυτ/ν). Основанием для такого предположения является то, что, согласно опытным данным [9], величина B зависит от безразмерной переменной величины lg(ksυτ/ν).

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

kond19.wmf (10)

На рис. 1 показано сравнение опытных и расчетных величин λ по формулам (9), (10). В диапазоне чисел Рейнольдса kond20.wmf, расчетные значения λ хорошо согласуются с опытными данными в трубах с песочной шероховатостью стенки трубы, причем это относится и к положению минимального значения величины λ при величине относительной шероховатости в пределах kond21.wmf. Расчетные значения величины λ для технической шероховатости R/ks = 1300 также удовлетворительно согласуются с опытными данными.

Число Рейнольдса, при котором шероховатую поверхность надо считать гидравлически гладкой, определим из условия равенства величин λ, рассчитываемых по формулам (1) и (9) [5]. В таблице приведены значения чисел Рейнольдса при различной степени песочной шероховатости стенки трубы.

kondr1.tif

Рис. 1. Закон сопротивления для гладких и шероховатых труб. Кривая 1 соответствует ламинарному режиму течения. Кривая 2 соответствует формуле Блазиуса. Кривая 3 соответствует формуле Прандтля. Кривые 4–10 расчет по формуле (9)

Граница перехода течения с гидравлически гладкого на переходной режим

R/ks

15

30,6

60

126

252

507

1300

Re

4, 58×103

104

2,32×104

5,15×104

1,12×105

9,5×104

kondr2.tif

Рис. 2. Зависимость величины В от lg(υτks/ν).. Кривая 1 соответствует режиму без проявления шероховатости. Кривая 2 соответствует режиму с полным проявлением шероховатости. Кривая 3 – расчет по формуле (6)

kondr3.tif

Рис. 3. Расчетные значения относительной скорости U/υτ, определенные по формуле (5), и их аппроксимация по формуле (11)

Используя расчетные значения λ можно рассчитать величину В по формуле (6). На рис. 2 представлены результаты расчетов величины B по (8) при различной степени песочной шероховатости.

Профиль относительной расчетной скорости U/υτ в трубе с произвольной песочной шероховатостью стенки трубы определяется по формуле (5), в которой величина B рассчитывалась по (6), а величина λ по (9), (10). В качестве абсциссы использовался предложенный в [1] безразмерный параметр kond22.wmf

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

Результаты численных расчетов с погрешностью 10 % аппроксимируются зависимостью

kond23.wmf (11)

Отметим, что ранее нами была рассмотрена эта же задача [5], при решении которой вначале была предложена эмпирическая зависимость для расчета величины В, а затем рассчитывался коэффициент гидравлического сопротивления и профиль относительной скорости. В вычислительном плане предложенная аппроксимация (10) значительно упрощает процедуру расчета.


Библиографическая ссылка

Кондратьев А.С., Ньа Т.Л., Швыдько П.П. ОБОБЩЕНИЕ ФОРМУЛЫ КОЛБРУКА – УАЙТА НА ТЕЧЕНИЯ ЖИДКОСТИ В ТРУБЕ С ПРОИЗВОЛЬНОЙ ПЕСОЧНОЙ ШЕРОХОВАТОСТЬЮ СТЕНКИ // Фундаментальные исследования. – 2017. – № 1. – С. 74-78;
URL: http://www.fundamental-research.ru/ru/article/view?id=41318 (дата обращения: 13.11.2019).

Предлагаем вашему вниманию журналы, издающиеся в издательстве «Академия Естествознания»
(Высокий импакт-фактор РИНЦ, тематика журналов охватывает все научные направления)

«Фундаментальные исследования» список ВАК ИФ РИНЦ = 1.074