Методическое пособие 780
.pdf
|
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