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

Научный журнал Российского НИИ проблем мелиорации, № 2(10), 2013 г., [130-156]

Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования «Саратовский государственный технический университет им. Гагарина Ю. А.» (ФГБОУ ВПО «СГТУ им. Гагарина Ю. А.»)

К вопросу об автомодельности течений в продольно-однородных турбулентных потоках

УДК 532.542

Высоцкий Лев Ильич - доктор технических наук, профессор, профессор, заслуженный деятель науки и техники РСФСР

Контактный телефон: 8-845-228-89-38

E-mail: vysotli@yandex.ru

Аннотации

В статье показано, что вопреки общепринятому мнению продольно-однородные турбулентные потоки наилучшим образом удовлетворяют принципу кусочной автомодельности. Выявлено, что в пределах каждой из пяти условных зон, на которые делятся указанные потоки, распределение осредненных скоростей является автомодельным. В то же время оказалось, что границы условных зон определяются значениями глобального числа Рейнольдса и относительной шероховатости. Указанные утверждения проверены с привлечением результатов обширных численных экспериментов. Диапазон чисел Рейнольдса составлял от 4Ч10³ до10¹¹. Относительная шероховатость при этом менялась от 0 до 0,1.

Ключевые слова: турбулентность, продольно-однородные потоки, осредненная скорость, автомодельность, относительная шероховатость.

Vysotskiy Lev Ilich - Doctor of Technical Sciences, Professor, Federal State Budget Educational Establishment of Higher Professional Education “Saratov State Technical University of Gagarin” (FSBEE HPE “SSTU Gagarin”), Professor, Honored Scientist of the RSFSR.

Contact telephone number: 8-845-228-89-38.

E-mail: vysotli@yandex.ru

ON THE ISSUE OF SELF-SIMILAR STREAMS IN THE LONGITUDINALLY UNIFORM TURBULENT FLOWS

The paper shows that longitudinally uniform turbulence flows give the best fit to the principle of piecewise self-similarity contrary to popular belief. It is revealed that the flows can be divided into five conventional zones. The averaged velocity distribution is self-similar within each of the zone. At the same time, it was found that the boundaries of the conditional zones are defined by values of global Reynolds number and by relative roughness. The specified statements are proved by a lot of numerical simulations. Range of Reynolds numbers was from 4Ч10³ to 10¹¹. Herewith the relative roughness changed from 0 to 0.1.

Keywords: turbulence, longitudinally uniform flows, averaged speed, self-similarity, relative roughness.

Введение

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

Тем не менее, экспериментальные данные вскрывали все новые недостатки предлагаемых зависимостей, что вынуждало исследователей вводить в них поправки. Как бы то не было, но достаточно прочно утвердилась в качестве как бы оптимальной трехслойная модель, в которой толщины условных зон турбулентного потока приняты фиксированными и в единицах внутреннего масштаба составляют 5/, 30/ и 70/. Причем в области значений расстояния от твердой стенки 70/ распределение осредненных скоростей считается автомодельным, то есть не зависящим от числа Рейнольдса. Опять-таки и это положение не выдержало экспериментальной проверки.

В данной статье автор вновь возвращается к изложенной проблеме, пытается на базе уточненного представления о процессах, происходящих вблизи твердой стенки на очень малых расстояниях от нее, предложить принципиально иное строение модели продольно-однородного турбулентного потока и иной взгляд на проблему автомодельности, во многом противоречащий современным представлениям. Автор надеется, что новая модель позволяет получить более надежные расчетные зависимости (для распределения осредненных скоростей, температуры и др.).

Исторические аспекты развития проблемы анализа распределения осредненных скоростей в турбулентных потоках излагались неоднократно. Существуют прекрасные обзоры в таких широко известных капитальных трудах, как «Механика жидкости и газа» Л. Г. Лойцянского [1], «Теория ламинарного пограничного слоя» Г. Шлихтинга [2], «Турбулентность, ее механизм и теория» И. О. Хинце [3] и многих других. С течением времени в подобного рода обзорах, которые чаще всего являются введением к изложению оригинального материала, появляются новые сведения, более современные трактовки, новые точки зрения, дополнительные критические замечания, новые данные, полученные с помощью ультрасовременных измерительных средств. К их числу можно отнести оценки Г. И. Баренблатта [4], Д. Остерлунда [5], М. Загаролы [6], И. Незу [7], В. Роди [8] и др. Здесь автор не преследует цель расширения оценки имеющихся в этой области достижений, однако предпослание основной части некоторого краткого введения считает необходимым.

Известно, что заметное поступательное движение в изучении турбулентности началось после обоснования и введения ряда основополагающих понятий. К ним относятся:

- разделение действительной скорости на осредненные и пульсационные составляющие;

- введение понятий о турбулентных (кажущихся) касательных и нормальных напряжениях [9];

- введение удачной (по крайней мере, для продольно-однородных турбулентных течений) связи касательных турбулентных напряжений - с градиентом осредненной скорости через коэффициент турбулентной вязкости [10];

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

- получение О. Рейнольдсом осредненных по времени уравнений Навье-Стокса, описывающих турбулентные течения;

- обоснование «энергоснабжения» турбулентных течений за счет отбора энергии у осредненного течения.

Однако набор перечисленных и разработка других понятий и представлений долгое время не позволяли решить даже простейшие практические задачи, связанные с турбулентным течением жидкостей и газов. Пожалуй, основной из них являлась задача о распределении осредненных скоростей по сечению труб и каналов.

Первоначально были предложены эмпирические формулы для распределения осредненных скоростей - знаменитая парабола Х. Базена [10] для широких каналов и формула Дарси [11] для расчета распределения осредненных скоростей в круглых трубах.

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

Тем не менее, по мере совершенствования измерительных средств и лавинообразного накопления экспериментальных данных появилась возможность сделать некоторые существенные выводы, разработать феноменологические модели строения турбулентного потока и предложить на их основе инженерные методы расчета турбулентных потоков в наиболее простых случаях, так называемых продольно-однородных течениях (течениях в круглых трубах и широких каналах). Эти случаи до сих пор продолжают оставаться важными объектами для изучения как с прикладной точки зрения, так и с научной, выступая в качестве своеобразного оселка для оттачивания новейших идей.

Большое значение имело установление того факта, что влияние твердых поверхностей, ограничивающих турбулентный поток, простирается лишь на сравнительно тонкий пристенный слой.

Такую гипотезу в 1930 г. высказал Т. Карман [12]. Он показал, что вне вязкого подслоя влиянием вязкости на распределение осредненных скоростей можно пренебречь. Это предположение дало ему возможность описать закон распределения скоростей простой формулой:

(1)

где ; ;

- постоянная интегрирования;

- динамическая скорость;

и - плотность и кинематическая вязкость жидкости;

- касательное напряжение на стенке;

- безразмерная скорость;

= - число Рейнольдса;

- расстояние от стенки;

- константа Кармана.

По Т. Карману, получалось, что закон (1) является универсальным для всего потока, кроме вязкого подслоя. Это означает, что входящие в него параметры не зависят от числа = (или = - в других случаях), где - средняя скорость; - диаметр трубы; - характерная длина для других потоков ( = для продольно однородного потока со свободной поверхностью, например).

Исходя из других положений, в 1932 г. Л. Прандтль [13] получил тот же закон. В «Механике сплошных сред» приведен вывод этого же закона на основе анализа размерностей, выполненного Л. Д. Ландау [14].

Сказанное позволило условно разделить весь турбулентный поток на основную его часть, свободную от влияния вязкости и называемую турбулентным ядром, и пограничный слой, где влияние вязкости велико.

Х. Людвиг и В. Тильман [15] подтвердили существование «закона стенки», который является универсальным, то есть независящим от изменения числа Рейнольдса и других параметров турбулентности. Универсальными масштабами при этом являются и

И. О. Хинце и Теннекес обнаруживали в своих опытах отклонение от универсального значения. Г. И. Баренблатт усомнился в полной автомодельности турбулентных течений и развил оригинальную теорию неполной автомодельности, что, в известном смысле, эквивалентно утверждению о влиянии вязкости жидкости при больших числах Рейнольдса. Последнее объясняется тем, что «турбулентный поток представляет собой совокупность громадного множества вихрей, пронизывающих движущуюся жидкость. С ростом числа Рейнольдса количество вихрей возрастает. Как бы велико не было число Рейнольдса, вязкость остается существенной вблизи «ядер» вихрей, и таким образом, ее динамическое воздействие на поток не исчезает». К этому можно добавить, что сам процесс каскадного дробления вихрей не является идеально организованным процессом. На границе соприкасающихся больших вихреобразований возникает множество мелких (вплоть до колмогоровских) вихорьков, что приводит к некоторой диссипации энергии турбулентности (сравнительно малой) и в толще турбулентного потока.

Далее при рассмотрении проблем теплопроводности и турбулентности Л. Прандтлем и Т. Карманом были введены понятия о «длине смешения» и даны формулы для ее вычисления. Эта идея оказалась настолько заразительной, что, несмотря на жесточайшую (уничтожающую) критику этого понятия, оно продолжает существовать, а сама формула для ее расчета снабжается все новыми уточнениями.

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

Так как уравнение (1) не удовлетворяет граничным условиям ни на твердой границе (при = 0), ни на оси потока (при = ), пришлось выходить из этого положения остроумным и, казалось бы, логически обоснованным путем, а именно, введением предположения о разделении турбулентного потока в круглой трубе на две зоны: турбулентное ядро и ламинарный пограничный слой.

В тонком ламинарном пограничном слое (его толщина была оценена в 11,6 /, где соотношение / играет роль внутреннего масштаба длины, а - динамическая скорость - масштаба скорости) существенно проявляется влияние вязкости. Обычно в гидроаэродинамике принято записывать получаемые соотношения в безразмерных координатах с нормировкой с помощью указанных масштабов длин и скоростей:

.(2)

Л. Прандтль и многие другие полагали (и полагают), что в ламинарном пограничном слое распределение осредненных скоростей является линейным:

= или = .(3)

Важное значение имело установление универсального логарифмического закона распределения осредненных скоростей, совпадающего с достаточно надежными экспериментальными данными, полученными И. Никурадзе [16] в опытах с круглыми трубами в диапазоне чисел Рейнольдса ( где - средняя скорость) от критического до 3,24•10⁶. Этот факт подтверждается графиком (рисунок 1), из которого следует, что опытные точки при всех числах удовлетворительно размещаются в окрестности прямой, соответствующей зависимости:

(4)

или

.(5)

Рисунок 1 - Универсальный логарифмический закон распределения скоростей в гладкой трубе: кривая (1) соответствует уравнению = , т. е. ламинарному течению; кривая (2) - переходу от ламинарной формы течения к турбулентной; кривая (3) - турбулентному течению при любых числах Рейнольдса; кривая (4) - турбулентному течению при < 10⁶, кривая (5) - уравнению = 11,5^1/10

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

При изменении числа Рейнольдса, естественно, изменяются как профиль эпюры осредненных скоростей, так и размеры областей, занятых пограничным слоем и турбулентным ядром. Если достаточно велико, то обычно считают взаимодействие между этими регионами пренебрежимо малым, но тогда в них должны существовать независимые решения
[, ]. По Л. Прандтлю около стенки профиль осредненных скоростей определяется лишь :

,(6)

а в турбулентном ядре лишь отношением , то есть:

.(7)

Х. Базен [10] в своих опытах установил фактически независимость закона распределения скоростей по сечению трубы от причин, обусловливающих касательное напряжение, выразив этот постулат формулой:

,(8)

где - скорость на оси трубы;

- скорость в текущей точке;

- динамическая скорость;

- постоянный коэффициент;

- расстояние от стенки;

- радиус трубы.

По-видимому, этот вывод Х. Базена послужил толчком к получению «универсальных» формул для распределения осредненных скоростей.

Это функциональное соотношение было указано также Стантоном в 1911 году [10] и известно как «закон дефицита скорости», а закон (6) известен под названием «закона стенки».

Таким образом, получается, что, если в обоих областях и осредненные скорости не зависят от числа Рейнольдса, то в них существует полное подобие, в противном случае существует неполное подобие.

С. В. Милликен в 1938 году предположил, что при достаточно больших числах Рейнольдса могут существовать перекрывающиеся области, где внутренний и внешний законы могут иметь место одновременно. Указанная область ограничивалась значениями 1 < <. При использовании внутреннего масштабирования переменных, получим:

(9)

где - аддитивная константа.

В терминах других масштабов тех же переменных окажется, что:

,(10)

где - также аддитивная константа.

Все три константы (, и ) определяются экспериментально. Часто используются их значения = 0,41; = 5,2; = 6,5. Обычно существование логарифмического закона распределения осредненных скоростей предполагается при значениях [17] в диапазоне:

.(11)

Позже Т. Карман [12] предложил трехслойную модель строения турбулентного потока. В ней примыкающий к гладкой стенке слой с линейным распределением осредненных скоростей имел толщину, равную , а скорость распределялась по закону:

.(12)

Буферная часть пограничного слоя ограничивалась расстояниями от стенки, равными -. В пределах нее осредненные скорости задавались формулой:

.(13)

При (турбулентное ядро) распределение скоростей подчиняется закону:

(14)

Из закона Дарси следует, что коэффициент Дарси определяется соотношением:

= 8 ,(15)

где - средняя скорость.

Для гладких труб Л. Прандтль получил (с учетом экспериментальных данных И. Никурадзе [16]) выражение:

(16)

Оно известно как прандтлев универсальный закон трения для гладких труб.

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

Однако нараставший объем экспериментальных данных заставил усомниться в непререкаемости мнения Л. Прандтля как в области распространения логарифмического закона, так и в обоснованности его теории «длины пути смешения».

Функциональное соотношение (7) называют законом дефекта скорост». Если не зависит от числа Рейнольдса, то в этой области существует полное подобие.

Между этими двумя областями могут существовать промежуточные или пересекающиеся области, где расстояние от стенки является большим по сравнению с внутренним масштабом /, но малым - по сравнению с другим масштабом длины .

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

Первым, кто исследовал эту проблему, был Г. И. Баренблатт [4]. Выполненный им анализ в предположении о неполной автомодельности течения по локальному числу Рейнольдса при отсутствии автомодельности по глобальному числу Рейнольдса , привел его к выражению для распределения осредненных скоростей в виде степенного закона:

(17)

где - показатель степени, зависящий от глобального числа Рейнольдса;

- характерный поперечный размер потока (диаметр трубы или глубина потока , например).

В качестве подтверждения справедливости этого результата Г. И. Баренблатт привлекает известные опытные данные И. Никурадзе, который получил хорошее соответствие в широком диапазоне чисел Рейнольдса опытных данных со степенным законом (рисунок 2) , где является функцией числа и изменяется в диапазоне от 6 до 10 при < 3,2·10⁶.

Известно, что степенные формулы описывают распределение осредненных скоростей ничуть не хуже, чем логарифмические [10]. Эти обстоятельства дали возможность Г. И. Баренблатту заключить, что теория неполной автомодельности является в такой же степени теоретическим обоснованием степенных формул, в какой для логарифмических - теория полной автомодельности.

В более поздних работах Г. И. Барентблатт и др. пришли к выводу, что при очень больших числах Рейнольдса степенной закон переходит в логарифмический.

Рисунок 2 - Графики зависимости по И. Никурадзе

При рассмотрении течения вблизи стенки Г. И. Баренблатт [4] вообще исключает область, в пределах которой вязкие и турбулентные касательные напряжения соизмеримы.

Еще раз напомним, что как логарифмический, так и степенные законы распределения осредненных скоростей не удовлетворяют граничным условиям ни на стенке, ни у оси потока и не могут описывать характер распределения скоростей вблизи стенки.

Г. И. Баренблатт [4] и другие, развивая теорию неполного динамического подобия, предложили ввести в формулу для распределения осредненных скоростей множитель и экспоненту , которые в виде асимптотического разложения по малой величине е представляются в виде:

;(18)

,(19)

где ;

и - эмпирические константы.

В одной из работ Г. И. Баренблатт используется концепция неполного подобия для зон турбулентных течений, включающих вязкий подслой и внешнюю часть потока [17]. С использованием анализа размерностей градиент осредненных скоростей был представлен в следующем виде:

.(20)

Это соотношение предполагается существующим при больших значениях в случае больших чисел Рейнольдса. При этом, если > ? и > ?, то (,) > 1/. То есть при очень больших значениях и пределом введенной безразмерной функции является 1/, что соответствует гипотезе Кармана о полном подобии.

Из предположения о неполном подобии следует, что конечный предел функции не существует, а она может быть представлена в виде [4]:

.(21)

В конечном счете зависимость для представляется в виде закона с масштабом, зависящим от числа [10]:

,(22)

где , и - универсальные константы.

С использованием данных Никурадзе для труб было определено:

.

При этих значениях формула (22) принимает вид:

(23)

или

.(24)

По Г. И. Баренблатту [17] показательный закон имеет место для подавляющей части профиля скорости:

40 < < 0,85 или (25)

Используя данные И. Никурадзе, Г. И. Баренблатт и В. М. Простокишин [18] нашли, что величины , и равны 0,577; 2,50 и 1,5 соответственно при пренебрежении всеми остальными. В этом случае показательный закон записывается следующим образом:

.(26)

На рисунке 3 М. В. Загарола и другие представили сопоставление расчетов по формуле (26) с их экспериментальными данными. По мнению указанных авторов, согласование этих результатов достаточно плохое, особенно при больших числах . Ими же делается предположение, что это может быть обусловлено недостаточно большим диапазоном чисел в опытах И. Никурадзе ( = 3,2•10⁶).

Из зависимости (26) легко усмотреть влияние чисел на величину , что как бы подтверждает выводы Г. И. Баренблатта.

Рисунок 3 - Сравнение расчетных и опытных значений осредненных скоростей

Однако, использование более полных, а главное, более точных опытов позволили М. В. Загароле [19] не согласиться с этим выводом при признании существования перекрывающихся зон. М. В. Загарола и другие полагают, что результаты, выполненные в Принстонском университете на установке «Большая труба» в широком диапазоне чисел Рейнольдса (на порядок перекрывающем диапазон опытов И. Никурадзе) от 31•10³ до 35•10⁶ экспериментов, тщательная подготовка к ним и обработка опытных данных, свидетельствуют о том, что при существенно больших числах Рейнольдса профили осредненных скоростей больше соответствуют логарифмическому закону их распределения, чем показательному. Это наводит авторов на мысль, что теория полного подобия предпочтительнее теории неполного подобия. По мнению авторов, результаты их опытов опровергают теории, с недавних пор развиваемые Г. И. Беренблаттом и другими Известно, что уже в 1982 г. Г. И. Баренблатт опубликовал обобщение своих воззрений [4]..

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

Пристенный слой имеет толщину , где - параметр, зависящий в общем случае от числа Рейнольдса и относительной шероховатости :

,

где = или (для круглой трубы или плоского потока).

В случае гладких стенок , а при шероховатых стенках .

Пристенный слой при гладких стенках делится на три зоны:

- примыкающий к гладкой стенке слой, толщиной , в пределах которого коэффициент турбулентной вязкости отрицателен, то есть < 0;

- среднюю часть пристенного слоя, ограниченного расстоянием между и от стенки, в пределах которого _т увеличивается от нулевого значения по зависимости , обращенной выпуклостью вниз;

- внешнюю часть пристенного слоя, ограниченную расстоянием между и от стенки, в пределах которого продолжает увеличиваться по зависимости , обращенной выпуклостью вверх.

Значения и также определяются параметрами и , то есть = (;) и = (;).

Принято, что соотношение : : = 1,5 : 30 : 70.

Принятие нового строения трехслойной модели пограничного слоя допустимо рассматривать в качестве развития представлений Т. Кармана. В частности, в предлагаемой модели слой жидкости толщиной с линейным распределением осредненных скоростей заменен слоем жидкости с переменной толщиной , в котором существует отрицательная турбулентная вязкость. Следующий слой, ограниченный по Т. Карману расстоянием от стенки в , заменен слоем с границей, отстоящей от стенки на . Граница всего пограничного слоя, часто определяемая его толщиной в , заменена на , при .

В отличие от модели Т. Кармана и других в данной модели толщины слоев пограничного слоя не являются фиксированными с помощью постоянных коэффициентов 5, 30, 70, а являются переменными, зависящими в общем случае от числа Рейнольдса и относительной шероховатости, и определяемыми переменными параметрами , , , которые подлежат отысканию.

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

,(27)

где = - число Кармана.

Можно записать, что:

(28)

и получить отсюда соотношение:

.(29)

Что касается турбулентного ядра, то оно делится на две части: нижнюю с логарифмическим и верхнюю с параболическим распределением осредненных скоростей.

Это означает, что для нижней части коэффициент турбулентной вязкости принят в интервале изменяющимся по линейному закону , а для верхней - постоянной величиной, то есть при ? • принято, что:

.(30)

Описание предлагаемой модели приведено выше. Для удобства выпишем полученные универсальные формулы для распределения осредненных скоростей. Они имеют вид:

- для - задается в табличной форме;

- для гладких стенок:

при :

(31)

при ? ?

;(32)

- для шероховатых стенок:

при ? ? :

2,5+,(33)

;(33а)

где =;

- высота эквивалентной шероховатости.

Для предлагаемых формул, определяющих значения осредненных скоростей, был выполнен численный эксперимент с целью выявления их универсальности по участкам: пристенный слой 0 ? ? , внутренняя часть турбулентного ядра ? ? и внешняя часть турбулентного ядра ? ? . С этой целью были вычислены в соответствии с изложенной методикой параметры для продольно-однородных плоских турбулентных потоков в плоских каналах при широком диапазоне чисел Рейнольдса (от 3000 до 10⁶) при относительных шероховатостях / = 10-²-10-⁶.

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

Затем была произведена выборка значений осредненных скоростей сначала по зонам пристенного слоя 0 ? ? ; ? ? ; ? ? = , которые для каждой из них были нанесены на график в координатах соответственно (рисунки 4 а, б, в):

.

а)

б)

в)

Рисунок 4 - Совмещенные эпюры распределения осредненных скоростей при диапазоне изменения чисел Re_H = 300010⁶ и относительных шероховатостей /H в интервале значений
/Н = 10-²10-⁶ а - для зоны 0 z _(-); б - для зоны z_с z z_в; в - для зоны z_в z _ст

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

Перейдем к внутренней зоне турбулентного ядра. Определим разность скоростей в точках с ? (при > ) и = , то есть его внешней границе. Получим:

Тогда получим:

(34)

Нанесем опытные точки на график в координатах и в пределах от до (рисунок 5).

Поскольку все точки ложатся на одну кривую, то отсюда следует вывод о существовании масштаба скорости и линейного масштаба /.

Повторим тот же прием для случая шероховатых стенок, когда ? :

.

Рисунок 5 - Совмещенный график зависимости и для всех значений Re_H и Д/Н

В этом случае дефицит скорости также равен:

.(35)

Результат совпадает с предыдущим, что логично, так как при = пристенный слой исчезает полностью.

Рассмотрим, наконец, внешнюю часть турбулентного ядра. Запишем выражения для и при · ? ? :

;

Тогда дефицит скорости составит:

(36)

Построим график в координатах и (рисунок 6).



Рисунок 6 - Совмещенные эпюры распределения осредненных скоростей во внешней зоне турбулентного ядра для данных в диапазоне изменения Re_H = 300010⁶ и относительной шероховатости /Н = 10-²10-⁶

Зависимости кусочно универсальны как во внутренних, так и во внешних переменных в рассмотренных зонах турбулентного потока.

Полученный график наглядно демонстрирует наличие масштаба для осредненной скорости в виде и линейного масштаба .

Итак, как же ответить на вопрос, влияет ли число Рейнольдса на масштабирование скорости и длины в проанализированном введении трех основных зон турбулентного продольно-однородного потока?

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

Это означает, во-первых, что значения и влияют на расстояния от стенки, в пределах которых осредненные скорости подчиняются той или иной универсальной закономерности. Такое подобие турбулентных потоков логично назвать кусочно автомодельным.

Во-вторых, при сравнительно малых числах Рейнольдса, когда = , второй зоны не существует вообще. Если = 1, то весь продольно-однородный поток представляется лишь пристенным слоем с = .

По мере увеличения числа параметр уменьшается.

Дальнейшее увеличение числа Рейнольдса приводит к уменьшению толщины пристенного слоя пока она остается большей высоты выступов шероховатости > . При условии равенства указанных значений зона сопротивлений становится квадратичной , а пристенный слой разрушается.

Выводы

1. Анализ результатов численного эксперимента, выполненного в широком диапазоне значений числа Рейнольдса и относительной эквивалентной шероховатости показал, что в пределах различных зон пятислойной модели продольно-однородного турбулентного потока имеет место автомодельность в распределении осредненных скоростей. Этот вывод распространяется как на продольно-однородные потоки так и на погранслойные течения при нулевом градиенте давления.

2. Толщины указанных зон вопреки общепринятому мнению, напротив, оказываются зависимыми как от глобального числа Рейнольдса, так и от относительной шероховатости.

3. Представленную модель продольно-однородного потока предлагается назвать кусочно автомодельной.

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

1. Лойцянский, Л. Г. Механика жидкости и газа / Л. Г. Лойцянский. - М.: Наука, 1970. - 904 с.

2. Шлихтинг, Г. Теория пограничного слоя / Г. Шлихтинг. - М.: Наука, 1974. - 712 с.

3. Хинце, И. О. Турбулентность, ее механизм и теория / И. О. Хинце. - М.: Физматгиз, 1963. - 612 с.

4. Баренблатт, Г. И. Подобие, автомодельность, промежуточная асимптотика / Г. И. Баренблатт. - Л.: Гидрометеоиздат, 1982. - С. 256-612.

5. Osterlund, J [Электронный ресурс]. - Stockholm, 1999. -Режим доступа: http:www.mech.kth.se/~jens/zpg/.

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

7. Nezu, I. Open channel flow measurements with a laser Doppler anemometer / I. Nezu, W. Rodi // Journal Hydr. Engrg. ASCE. - 1986. - № 112(5). - P. 335-355.

8. Роди, В. Модели турбулентности окружающей среды / В. Роди // Методы расчета турбулентных течений. - М.: Мир, 1984. - С. 227-322.

9. Lorentz, H. A. Uber die Enstehung turbulenter Flussigkeitsbewegungen und uber den Einflyss dieser Bewegungen bei der Stromung duren Rohren / H. A. Lorentz // Abh. Theor. Phys. - 1907. - Lgz. Bd 1. - P. 43-71.

10. Bazin, Х. Experiences nouvelles sur la distribution des des vitesses dans les tuyaux / Х. Bazin // Memoires d'Acad. des Sciences de l'Inst. de France. - 1902. - № 6.

11. Darcy H. Recherches expйrimentales relatives aux movement de l'eau dans les tuyaux des conduits / H. Darcy. - Paris, 1858.

12. Karman, Т. Y. Mechanical similarity and turbulence (in German) / Т. Y. Karman // Nachrichten von der Gesellschaften der Wissenschaften zu Giittingen. Mathematisch Physikalische Кiasse. - 1930. - Р. 56-76.

13. Prandtl, L. Nenere Ergebnisse der Turbulenzforschung / L. Prandtl // V. D. I. - 1933. - № 5. - P. 107-110.

14. Ландау, Л. Д. Механика сплошных сред / Л. Д. Ландау, Е. М. Лифшиц. - М.: Гостехтеориздат, 1953. - 736 с.

15. Tollmien, W. Ein allgemeines Kriterium der Instabilitat laminarer Geschwindigkeitsverteilung / W. Tollmien // Nachr. Ges. Wiss. Gottingen. Math. Phys. Klasse. Fachgruppe I. 1. - 1935. - P. 79-114.

16. Nikuradse, J. Strdmungsgeseize in rauhen Rohren / J. Nikuradse // Forsch. Geb. Ing.-Wes. Heft 361. - Berlin, 1933.

17. Barenblatt, G. I. Self-Similar Intermediate Structures in Turbulent Boundary Layers at Large Reynolds Numbers / G. I. Barenblatt, A. J. Chorin, V. M. Prostokishin // Journal Fluid Mech. - 2000. - № 410. - P. 263-283.

18. Barenblatt, G. I. Scailing laws in fully daveloped schear flows. Part 2. Processing of experimental data / G. I. Barenblatt, V. M. Prostokishin // Journal Fluid Mech. - 1993. - № 248. - P. 521-529.

19. Zagarola, M. V. A new mean velocity scaling for turbulent boundaru layers / M. V. Zagarola, A. J. Smits // Prouedings of FEDSM'98 DC. A. SME. FEDSM. -1998. - 4950.

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

...