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

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

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

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

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

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

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

In the paper new analytical dependences describing the distribution of averaged velocities in longitudinal homogeneous flow validated for a wide range of the reynolds number and the relative equivalent roughness. The methodology for calculation the dimensionless averaged velocities accompanies by explanations. To implement the input there are two key parameters, the reynolds number and the relative equivalent roughness. In order to receive dimensional values the additional input of characteristic length ( - radius for round pipes, - depth for flat flow) and dynamic velocity , is required.

Keywords: negative turbulent viscosity, averaged velocity, maximum velocity, distribution law.

Одной из основных и, казалось бы, простейших задач гидравлики является получение закона распределения осредненных скоростей в турбулентных продольно однородных потоках. Однако до сего времени точного решения этой задачи не существует, что связано со сложностями, вызванными турбулентностью потоков жидкости и газа. Тем не менее, продолжается разработка все новых моделей продольно однородных потоков, и предлагаются все новые формулы для расчета распределения осредненных скоростей по нормали к стенке в них. Авторами предлагается новая модель продольно однородного турбулентного потока, основанная на открытии (диплом № 376), в которой в составе пристенного слоя потока жидкости предполагается наличие очень тонкого подслоя с отрицательной турбулентной вязкостью [1]. Следует упомянуть, что предлагаемая модель является пятислойной. Весь поток делится на пристенный слой и турбулентное ядро. Пристенный слой делится на три зоны: слой с отрицательной турбулентной вязкостью , толщиной ; последующие два слоя отличаются характером изменения в них коэффициента . Турбулентное ядро делится на две зоны. Примыкающая к пристенному слою зона с линейным изменением :

(1)

Где - постоянная кармана;

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

- расстояние по нормали от стенки.

В этой зоне имеет место логарифмический закон распределения скоростей.

Во второй зоне турбулентного ядра коэффициент турбулентной вязкости принят постоянным:

(2)

Где - расстояние от стенки до границы раздела двух зон в турбулентном ядре;

= const;

- характерная длина ( - радиус для круглых труб, - глубина для плоских потоков).

В этой зоне распределение осредненных скоростей является параболическим.

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

(3)

(4)

Реализация данных зависимостей привела к установлению формулы для распределения осредненных скоростей [2].

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

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

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

- в первой зоне ;

- во второй зоне ;

- третьей зоне

Соотношение расстояний от стенки принято равным 1,5 : 30 : 70.

Данные для расчета осредненных скоростей в пределах пристенного слоя приведены в таблице 1.

Таблица 1 - значение безразмерных осредненных скоростей в пристенном слое

2 в турбулентном ядре формулы для распределения осредненных скоростей имеют вид:

- при :

; (5)

- при :

; (6)

- при

; (7)

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

- при :

; (8)

- при :

; (9)

- при :

; (10)

Где ;

= 0,4 - константа кармана;

(или ) - характерные длины для потоков в трубе или со свободной поверхностью;

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

(или ) - расстояние от стенки по нормали к ней.

Из найденных соотношений легко получить формулы для определения относительной средней скорости :

- при :

; (11)

- при :

; (12)

- при = 1:

. (13)

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

- при

(14)

А если принято =0,4, то:

(15)

- при

(16)

В случае плоских потоков те же формулы имеют вид:

, (17)

А формула для максимальной скорости имеет тот же вид, что и в случае круглой трубы:

(18)

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

- длина, входящая в размерность кинематической турбулентной вязкости , представляет собой просто расстояние от стенки (как это первоначально и было предложено л. Д. Ландау [3]);

- у гладкой стенки ламинарный пограничный слой отсутствует;

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

- универсальная константа кармана принимается таковой и равной = 0,4;

- влияние шероховатости учитывается эквивалентной шероховатостью ;

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

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

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

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

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

- турбулентное ядро условно делится в соответствии со схематизацией эпюры распределения турбулентной вязкости или на две части. В примыкающей к пристенному слою части, где принято, что (или ), эпюра распределения осредненных скоростей естественным образом становится логарифмической, а в остальной части, где принято const, - параболической. Эти части эпюры смыкаются с обеспечением непрерывности в распределении и . Первая часть простирается от стенки на расстояние (или ), причем при больших =0,1735;

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

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

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

- предлагаемые формулы получены интегрированием исходного уравнения: скорость итерация ньютон уравнение

или .

В пределах пристенного слоя, то есть при изменении , интегрирование велось численно, а при - аналитически;

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

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

Задача выявления указанной зависимости решается с последующим использованием вычислительного эксперимента с опорой на следующие наиболее надежные теоретические и экспериментальные факты:

- средняя скорость связана с коэффициентом дарси соотношениями:

при = ; (19)

при = ; (20)

- дефицит средней скорости составляет:

при = ; (21)

при = , (22)

Далее имеем:

, (23)

Где ,

Или

 при = ; (24)

при = , . (25)

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

и . (26)

В подготовленных для решения этих задач нелинейных уравнениях методом ньютона исходные выражения имеют вид:

(27)

(28)

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

 (29)

(30)

; (31)

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

Итерационная формула ньютона имеет вид:

; (32)

, (33)

Где

и - нулевые приближения и соответственно.

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

. (34)

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

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

Используем далее в рассмотренном случае уравнение для . Для построения итерационного процесса по ньютону можно использовать одну из производных от , то есть или . Примем .

Выражение для имеет вид, несколько отличный от выражения (27), так как , а в выражении для пропадает .

Тогда итерационные значения определяются по формуле:

. (35)

При определении и расчеты по приведенным формулам не вызывают затруднений.

На основании вышеизложенного можно сделать следующие выводы:

- впервые предложена универсальная формула для расчета распределения осредненных скоростей в продольно-однородных потоках;

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

- формулы прошли широкую апробацию.

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

1 высоцкий, л. И. Явление возникновения течения с отрицательной турбулентной вязкостью в продольно-однородном турбулентном потоке жидкости / л. И. Высоцкий, и. С. Высоцкий // научные открытия 2009: сборник кратких описаний, научных открытий, научных гипотез. - м.: раен, 2010. - с. 25-27.

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

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

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

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