Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Методическое пособие 780

.pdf
Скачиваний:
14
Добавлен:
30.04.2022
Размер:
10.22 Mб
Скачать

 

p

p

r2

 

v

1

2

C.

(3.8)

 

4l

 

 

 

 

Постоянную интегрирования C определим из условия на стенке v 0 при r r0 . Подставляя в выражение для v , получим формулу Пуазейля (1840)

v

p1 p2

r02 r2 .

(3.9)

4 l

 

 

 

Согласно формуле Пуазейля эпюра скоростей в поперечном сечении трубы имеет форму параболы (рис. 3.4). Максимальная скорость наблюдается при r 0 , здесь

v

 

p

p

r2

(3.10)

1

2

0 .

max

 

4l

 

 

 

 

 

 

 

Расход в трубе можно определить интегрированием по сечению трубы элементарных расходов, которые равны произведению скорости (3.10) на площадь элементарного кольца 2 rdr

 

 

r0

 

 

r4

p

p

 

Q v2 r

0

 

1

2

.

 

(3.11)

 

8

 

 

 

 

0

 

 

 

 

l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Средняя скорость в трубе

 

 

 

 

 

 

 

 

 

 

v

 

Q

 

p1

p2 r02

 

vmax

.

(3.12)

r

 

8 l

 

 

ср

 

2

 

 

2

 

 

 

0

 

 

 

 

 

 

 

 

 

 

 

Используя (3.12) легко определить величину гидравличе-

ского коэффициента трения тр

 

по формуле Дарси (3.2). Дей-

ствительно, принимая, во внимание, что

 

 

 

p1 p2 hl ,

 

 

 

 

(3.13)

 

 

 

 

,

 

 

 

 

(3.14)

 

 

 

 

 

 

 

g

 

 

 

 

 

 

 

 

 

 

r

D

,

 

 

 

 

 

(3.15)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0

2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

101

получаем

 

 

64

 

64

.

(3.16)

 

 

тр

 

vD

 

Re

 

 

 

 

 

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

3.3. Некоторые частные случаи ламинарного течения жидкости

3.3.1. Ламинарное течение жидкости в круглой трубе

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

Рис. 3.5

Уравнение Бернулли для выбранных сечений будет иметь

вид

P

 

P

h

 

 

1

2

,

(3.17)

 

 

 

 

 

тр

 

 

 

 

 

 

где hтр – потеря напора на трение. Тогда

102

h

P P

 

Pтр

 

 

1 2

 

.

(3.18)

 

 

тр

 

 

 

 

 

 

 

 

 

В потоке жидкости выделим контрольный объём в форме цилиндра, как показано на рисунке, и составим уравнение равномерного движения его. Уравнение будет иметь следующий вид

 

 

( P

P ) r2 2 rl 0 .

 

 

 

1

2

 

 

 

Здесь ( P

P ) r 2

– суммарная составляющая сил давления;

1

2

 

 

 

 

 

2 rl – силы сопротивления.

 

 

 

Отсюда

 

 

 

 

 

 

 

 

 

Pтр r

.

(3.19)

 

 

 

 

 

 

 

 

2l

 

Выразим касательное напряжение по закону Ньютона, т.е.

 

dv

.

(3.20)

 

 

dr

 

Знак минуса обусловлен тем, что направление отсчёта радиуса противоположное ранее принятому направлению отсчёта нормали. Получаем

 

Pтр r

 

dv

,

(3.21)

 

2l

 

dr

 

 

 

 

 

 

 

 

 

dv

Pтр

 

rdr,

(3.22)

2 l

 

 

 

 

 

 

 

 

v

 

Pтр

 

r2

 

C.

(3.23)

2 l

2

 

 

 

 

 

 

 

Постоянную интегрирования определим из условий, заданных на стенке трубы, где при r r0 ; v 0.

C

Pтр

r2 .

(3.24)

4 l

 

0

 

103

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

 

 

 

 

v

 

Pтр

 

r02 r2 ,

 

 

 

(3.25)

 

 

 

 

 

4 l

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v

 

 

 

 

Pтр

r2 .

 

 

 

(3.26)

 

 

 

 

 

 

 

4 l

 

 

 

 

 

 

 

max

 

 

 

0

 

 

 

 

Следовательно, расход жидкости составит

 

 

 

Q

 

dQ

 

vds

 

 

Pтр

( r

2 r 2

)2 rdr

Pтр

r 4

. (3.27)

 

 

 

 

 

 

 

 

4 l

0

 

8l

0

 

 

 

 

 

 

 

 

 

 

( Q )

 

( S )

( r )

 

 

 

 

 

 

 

 

 

 

Среднее по сечению значение скорости будет

v

 

Q

 

 

Pтр

r2

,

r

 

8l

cp

 

2

 

0

 

 

 

0

 

 

 

 

т.е. vcp 12 vmax .

Получим закон сопротивления, т.е. зависимость расхода и геометрии трубы. Из (3.28) имеем

Pтр 8 l Q ,

r04

h

 

Pтр

 

8 lQ

 

128vlQ

.

 

r4

gd 4

тр

 

 

 

 

 

 

 

0

 

 

 

(3.28)

hтр от

(3.29)

(3.30)

Формула (3.30) получена Пуазейлем экспериментально в 1840 г. Получим так называемый закон Пуазейля-Гагена, который используется для расчёта трубопроводов с ламинарным режимом течения. Приведём уравнение (3.30) к виду

h

 

l

 

vcp2

 

 

 

 

.

(3.31)

 

 

тр

 

d 2g

 

 

 

 

104

С этой целью в формуле (3.30) заменим расход через про-

изведение 4d 2 vcp , после соответствующих преобразований бу-

дем иметь

 

 

 

 

 

 

32vlvcp

 

64v

 

L

 

vcp2

 

64 l vcp2

 

 

l

 

vcp2

 

 

 

h

 

 

 

 

 

 

 

 

 

 

 

 

 

, (3.32)

 

 

 

gd 2

 

 

 

 

 

 

 

 

 

 

 

 

 

 

тр

 

 

 

v d d 2g Re d 2g

 

л d 2g

 

 

 

 

 

 

 

 

 

cp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

64

.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

л

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Re

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Здесь следует обратить внимание, что hтр при ламинарном

течении пропорциональна скорости в первой степени, а квадрат скорости в формуле (3.32) получен искусственно.

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

 

v3ds

 

1

 

r0

 

 

r2

3

rdr

 

( S )

 

 

v3ds 16 1

 

 

 

 

.

3

3

2

2

 

vcp s

 

vcp s

( s )

0

 

 

r0

 

r0

После замены переменной

1 r2 z ,

r02

получим

(3.33)

(3.34)

8 10 z3dz 2 .

(3.35)

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

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

105

 

v2ds

 

4

.

(3.36)

v2

s

3

( S )

cp

 

 

 

 

 

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

3.3.2. Начальный участок ламинарного течения

Формирование параболического профиля скоростей на начальном участке ламинарного течения в круглой трубе под действием сил вязкости происходит как показано на рис. 3.6.

Рис. 3.6. Формирование параболического профиля скоростей на начальном участке ламинарного течения в круглой трубе

Расстояние от начала трубы, на котором происходит стабилизация параболического профиля скоростей, называется

начальным участком течения ( lнач. ). За пределом этого участка

мы имеем стабилизированное ламинарное течение. Для определения длинны начального участка используется формула Шиллера

 

lнач

0, 029 Re .

(3.37)

 

 

 

d

 

Тогда максимально возможная длина начального участка

составит

 

lнач.max d 0, 029 Reкр 66, 5d .

(3.38)

106

Т.к. значение градиента скорости у стенки на начальном участке больше, чем на участке стабилизированного ламинарного течения, то сопротивление начального участка будет больше и составит

h k

 

x vcp2

(3.39)

 

 

 

 

,

 

 

тр

л d 2g

 

x

где k 1 и является функцией Re d , с увеличением которой он

уменьшается и при x lнач становится равным 1,09, т.е. сопро-

тивление всего начального участка на 9 % больше, чем сопротивление такого же участка трубы, взятого в области стабилизированного ламинарного течения.

3.3.3. Ламинарное течение жидкости в зазорах

Рассмотрим ламинарное течение жидкости в зазоре, образованном двумя параллельными твёрдыми стенками, расстояние между которыми равно a . Выделим контрольный объём жидкости в зазоре в форме прямоугольного параллелепипеда, как показано на рис. 3.7, и запишем условие его равномерного движения вдоль оси x

Рис. 3.7. Ламинарное течение потока в зазоре

107

Т.к. движение

2 yP

 

dv

2l ;

 

тр

 

dy

1

 

 

 

P P P .

 

тр

1 2

 

 

Т.к. при y a2

получаем

является равномерным вдоль оси x, то

dv Pтр ydy;

l

v 0, то c Pтр a2

2 l 4

 

P

a2

 

 

v

тр

 

 

y2

.

 

 

 

2 l

4

 

 

v Pтр y2 C .Здесь

2 l

, откуда окончательно

(3.40)

А теперь получим значение расхода,

Q

dQ vds

P

a2

 

 

 

P a3

 

тр

 

 

y2

 

2dy

тр

. (3.41)

 

 

 

Q

s

y 2 l

4

 

 

 

12l

 

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

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

P

 

Q 12l

 

12lvcp

.

(3.42)

a3

a2

тр

 

 

 

 

Рассмотрим случай, когда одна из стенок движется, а давление вдоль зазора постоянно. При этом возникает так называемое фрикционное безнапорное движение. Из уравнения равновесия элемента жидкости следует d 0 Тогда const . Из уравнения Ньютона получаем

 

dv

const ,

(3.43)

dy

 

 

 

v

С

y С .

(3.44)

 

 

 

 

1

 

 

 

 

 

108

Рис. 3.8. Случай с движущейся стенкой

Для определения постоянных интегрирования C и C1 ис-

пользуем следующие граничные условия: при

y

a

;

v 0 и

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

при

 

 

y

a

 

v V . Отсюда получаем

С

V

;

 

С

V

;

 

 

 

 

 

 

 

 

 

 

 

 

2

 

 

 

 

a

 

1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

v

1

 

y

 

 

 

 

 

 

 

 

 

 

 

 

 

 

V .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2

 

2

 

 

 

 

 

 

 

 

 

 

 

Расход жидкости при этом составит

 

Q

V

a .

(3.45)

 

2

 

 

Т.е. имеет место течения Куэтта. Если перемещение стенки происходит при наличии перепада давления, то закон распределения скорости как алгебраическая сумма ранее полученных выражений, а именно

v

Pтр

a2

y

2

 

 

1

 

y

 

 

 

 

 

 

 

 

 

.

(3.46)

 

 

 

 

 

 

2 l

4

 

 

 

 

2

 

a

 

Соответственно выражение для расхода будет иметь вид

Q

Pтр a3

V

a .

(3.47)

12 l

 

2

 

 

 

109

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

3.3.4. Течение с теплообменом

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

ются в уточнении.

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

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

110