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

Размещено на http://www.allbest.ru/

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

Л.И. Выcоцкий

Саратовский государственный технический университет им. Гагарина Ю. А., Саратов, Российская Федерация

В статье излагаются результаты сопоставительного анализа формул для определения значения коэффициента Дарси в продольно-однородных турбулентных потоках (круглых трубах). С этой целью были использованы старые (традиционные) их варианты и новейшие, полученные на установках с использованием «принстонской» и «орегонской» труб, а также с применением современных высокоточных измерительных методик. Сопоставление с опытными данными, полученными на уникальных установках, позволило сформулировать соответствующие рекомендации по применению формул тех или иных авторов.

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

The article presents the results of comparative analysis of formulas to determine the value of Darcy coefficient in longitudinally uniform turbulent flows (round pipes). For this purpose the old (traditional) formulas and the latest ones obtained by the Princeton “Superpipe” and Oregon device, as well as modern precision measuring techniques were used. Comparison with the experimental data obtained by the unique installations has allowed laying down guidance for the use of formulas created by different authors.

Ключевые слова: smooth pipes, Darcy coefficient, distribution of averaged velocities, wind tunnel, liquefied gases.

В предыдущих публикациях автора в соавторстве с И. С. Высоцким [1, 2] предложен новый метод расчета распределения осредненных скоростей в продольно-однородных турбулентных потоках (ПОТП), в частности круглых трубах. Вычислительная процедура предусматривает введение в качестве исходных данных значения числа Рейнольдса и относительной эквивалентной шероховатости , где - высота выступов эквивалентной шероховатости, - радиус трубы.

Существенную роль в вычислительном процессе играет значение коэффициента Дарси . Известно, что в общем случае при равномерном течении в круглых трубах коэффициент Дарси зависит от числа и .

В случае гладких труб, то есть при = 0, как широко известно, коэффициент Дарси определяется лишь значением числа Рейнольдса. Естественно, пользователь, желающий произвести расчет распределения осредненных скоростей в ПОТП, нуждается в рекомендациях относительно использования формул для определения коэффициента Дарси, соответствующего исходным данным (при гладких трубах это значение числа Рейнольдса ). На протяжении более двухсот лет многие исследователи пытались установить связь между и при условии = 0. Подобные исследования продолжаются до сих пор. В последнее время были проведены дополнительные уникальные опыты по определению зависимости коэффициента Дарси от числа Рейнольдса с использованием «принстонской» и «орегонской» труб. Результаты этих опытов опубликованы в трудах R. J. Donnelly, G. G. Ihas, B. J. McKoen, A. J. Smits, C. J. Swanson, M. V. Zagarola [3-6].

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

- формула Блазиуса:

турбулентный поток коэффициент дарси

; (1)

- формула И. Никурадзе [7], которую предлагается использовать при больших числах:

; (2)

- формула Кольбрука [8], которая при условии = 0 имеет вид:

; (3)

- формула, приведенная в работах B. J. McKoen, A. J. Smits, M. V. Za-garola [4]:

(4)

- уточненная формула для больших чисел Рейнольдса [5]:

. (5)

Формула Блазиуса была получена им в 1908 г. в предположении, что осредненные скорости распределяются в круглой трубе по степенному закону с показателем степени = 1/7. Имеются различные рекомендации, ограничивающие область применимости формулы Блазиуса для гладких круглых труб по диапазону допустимых значений числа Рейнольдса. Наиболее распространенные рекомендации сводятся к следующим указаниям:

а) 2300 10;

б) 4000 10;

в) 4000 40;

г) 4000 100000.

Заметим, что указания типа а), б), в) плохо обусловлены. На самом деле, при гладких трубах, то есть при = 0, получается, что эти условия справа переходят в требование . Другими словами, выходит, что формула Блазиуса справедлива во всем диапазоне значений числа Рейнольдса, в то время как сам Блазиус ограничил применимость своей формулы значениями 100000.

Формула Кольбрука - Уайта [8] была предложена в 1939 г. и относится к универсальным формулам, то есть к формулам, которые можно использовать как при гладких, так и при шероховатых стенках и при доквадратичной зоне сопротивления. Однако эта формула, как и ей подобные (А. Д. Альтшуля например), имеет асимптотический характер и в силу этого должна иметь условие, ограничивающее ее применение по числу Рейнольдса. Как известно, таким условием является значение критического числа Рейнольдса, соответствующего наступлению квадратичной зоны сопротивления, например в виде:

= 500.

В предыдущей публикации автора (2011) [1] это условие, полученное при обработке результатов вычислительного эксперимента, имеет схожую структуру:

= 400.

Оно было несколько уточнено и представлено в виде:

(6)

Все критерии удовлетворяют естественному условию: при = 0 .

Результаты сопоставительного анализа

Наиболее известными исследованиями зависимости коэффициента Дарси, выполненными в последнее время, являются опыты, осуществленные группами исследователей в Орегоне [5] и Принстоне [3, 4] с использованием аэродинамических труб. Разница в их размерах значительна. Если в Принстоне масса аэродинамической трубы (superpipe) составляла примерно 25 тонн, то в Орегоне масса трубы составляла всего около 1 унции ( 30 граммов). В Принстоне в качестве рабочей среды использовался сжатый воздух. В Орегоне при сравнительно небольших значениях числа Рейнольдса использовались гелий, кислород, азот и др., а при больших - жидкий гелий [5]. Несмотря на грандиозное различие в размерах установок, результаты исследований имеют достаточно малые расхождения. Полагая, что при получении цитируемых результатов использовалась наиболее совершенная техника проведения опытов и обработки полученных данных, проведем сопоставительный анализ применения приведенных ранее формул (1)-(6) для расчета значений коэффициента Дарси. Сопоставление выполним в табличной форме с указанием процентов расхождения расчетных и опытных данных (таблицы 1, 2).

Анализ табличных данных позволяет сделать ряд заключений:

- сделанное в 1908 году Блазиусом предположение о степенном законе распределения осредненных скоростей при = 1/7 обеспечивает расхождение с орегонскими опытными данными в пределах 2 % при диапазоне изменения значений числа от 3000 до 140000. По принстонским данным та же формула дает тот же разброс до = 145000 (в таблице 2 условные границы выделены жирными линиями);

- формула И. Никурадзе, как и предполагал ее автор, дает результаты с разбросом не более 3 % в диапазоне значений чисел Рейнольдса от 30000 до 587000 (при сопоставлении с орегонскими данными). В случае принстонских данных - с тем же разбросом в диапазоне чисел Рейнольдса от 74000 до 1000000;

- формула Кольбрука - Уайта с той же точностью соответствует орегонским данным при значениях чисел Рейнольдса более 3000, а принстонским данным - при значениях чисел Рейнольдса от 31000 до 30000000;

- формула (4) при том же разбросе соответствует диапазону значений числа Рейнольдса от 3000 до 587000 с орегонскими данными и диапазону от 31000 до 35000000 - с принстонскими данными, причем при
31000 30000000 расхождение с опытными данными не превышает 1 %;

- формула (5) рекомендована для применения при числах Рейнольдса , больших 74000, в этом случае она приводит к погрешности, меньшей 1 %; неплохие результаты (с погрешностью в пределах 3 %) она дает и при значениях , превышающих 6000 (по орегонским данным).

Размещено на http://www.allbest.ru/

Таблица 1 - Сопоставление расчетных и опытных данных применительно к условиям опытов в Принстоне

Таблица 2 - Сопоставление расчетных и опытных данных применительно к условиям опытов в Орегоне [5]

Размещено на http://www.allbest.ru/

В результате выполненного сопоставительного анализа по установлению соответствия рассмотренных формул новейшим опытным данным, полученным с использованием современных высокоточных измерительных методик, выявлено, что формулы (4) и (5), предложенные авторами новейших опытов [3-6], имеют наивысшую точность. При числах Рейнольдса, больших 145000, некоторое преимущество имеет формула (5).

Одновременно можно выразить восхищение результатами, полученными авторами формул (1), (2) и (3), скрупулезностью и ответственностью их подходов к решению проблемы.

При реализации предложенного автором в 2005 г. метода расчета распределения осредненных скоростей в круглых трубах [2] рекомендуется использование зависимости коэффициента Дарси от числа Рейнольдса, выраженной формулой Кольбрука - Уайта, так как изложенный выше анализ подтверждает правильность ее выбора для этой цели. Она лишь немногим уступает формуле (4), но применима не только для гладких труб, но и для квадратичной зоны сопротивления. Особо подчеркнем, что при расчете распределения скоростей по предлагаемому методу можно использовать и результаты непосредственного определения коэффициента Дарси опытным путем.

Список использованных источников

1 Высоцкий, Л. И. Продольно-однородные осредненные турбулентные потоки / Л. И. Высоцкий, И. С. Высоцкий. - Саратов: СГТУ, 2011. - 560 с.

2 Высоцкий, Л. И. Построение сквозной для всех зон сопротивления формулы для распределения осредненных скоростей в продольно-однородных турбулентных потоках / Л. И. Высоцкий // Совершенствование методов гидравлических расчетов водопропускных и очистных сооружений: межвуз. науч. сб. / СГТУ. - Саратов: СГТУ, 2005. - С. 7-65.

3 Nikuradze, J. Gezetzmassgkeiten dez turbulenten Stromung in glatten Rohren / J. Nikuradse // VDI. Forschungsheft. - 1932. - № 356 // Рус. пер. в сб.: Проблемы турбулентности / под ред. М. А. Великанова, Н. Т. Швейковского. - М. - Л.: ОНТИМ, 1936. - С. 75-150.

4 Colebrook, C. F. Turbulence flow in pipes with particular Reference to the transition region between the smooth and rough pipe lines / C. F. Colebrook // Journal of Institute of Civil Engineering. - 1939. - P. 133-156.

5 Friction factor for smooth pipe flow / B. J. McKoen, C. J. Swanson, M. V. Zagarola, R. J. Donnelly, A. J. Smits // J. Fluid Mech. - 2004. - V. 511. - P. 41-44.

6 McKoen, B. J. A new friction factor relationship for fully developed pipe flow / B. J. McKoen, M. V. Zagarola, A. J. Smits // J. Fluid Mech. - 2005. - V. 538. - P. 429-443.

7 Pipe flow measurement over a wide range of Reynolds numbers using liquid helium and various gases / C. J. Swanson, B. Julian, G. G. Ihas, R. J. Donnelly // J. Fluid Mech. - 2002. - V. 461. - P. 51-60.

8 Zagarola, M. V. Mean-flow sca1ing of turbulent pipe flow / M. V. Zagaro1a // Journal of F1uid Mechanics. - 1998. - Vol. 373. - Р. 33-79.

Размещено на Allbest.ru