05 семестр / Книги и методические указания / Панаиотти С.С. Лекции по гидромеханике / Панаиотти С.С. Лекции по гидромеханике. Часть 3 - Учебное пособие. - Калуга КФ МГТУ им. Н.Э. Баумана, 2003
.pdf
ГЛАВА 3. ГИДРОДИНАМИКА ВЯЗКОЙ ЖИДКОСТИ
3.1. Закон вязкого трения Ньютона
При движении жидкости прямолинейными слоями, касательные напряжения, возникающие между слоями, пропорциональны градиенту скорости. Пусть жидкость движется параллельными слоями вдоль твердой стенки
(рис. 3.1).
Распределение скоростей показано на этом рисунке. В вязкой жидкости скорость на твердой стенке Ox равна нулю. Частицы жидкости прилипают к стенке – это свойство называется условием прилипания. В слое, расположенном на расстоянии y от начала координат – скорость равна V . В соседнем слое скорость V + V . Выберем на рис. 3.1 бесконечно малую частицу жидкости в виде прямоугольника abcd. За бесконечно малый промежуток времени t частица жидкости переместится в новое положение и примет форму параллелограмма a′b′c′d ′. В соответствии с законом Ньютона касательные напряжения между слоями жидкости:
τ = µ |
∂V |
, |
|
(3.1) |
|
|
|
||||
|
∂y |
|
|||
а касательная сила |
|
||||
T = Sµ |
∂V |
, |
(3.2) |
||
|
|||||
|
|
∂y |
|
||
где S — площадь слоя жидкости, — динамическая вязкость.
Рис. 3.1. К закону вязкого трения Ньютона
1
Отметим два отличия законов трения в жидкости от таковых для твердых тел:
1.Касательная сила T не зависит от нормального давления.
2.Касательное напряжение в твердом теле в соответствии с законом Гука пропорционально деформации, а в жидкости оно пропорционально скорости деформации. Например, чем больше угол закручивания твердого стержня,
тем больше напряжения. Покажем, что ∂V
∂y в уравнении — это скорость деформации.
Из параллелограмма a′b′c′d ′ очевидно, что tgΔθ ≈ Δθ = V t . y
Из этого соотношения следует, что |
∂V |
= lim |
V |
= lim |
Δθ |
= |
∂θ |
. Т.е. про- |
∂y |
|
|
|
|||||
|
y →0 y |
t →0 t ∂t |
||||||
изводная ∂V
∂y в (3.1) — это скорость скошения прямого угла или скорость угловой деформации.
3.2. Уравнения Навье-Стокса
По сравнению с идеальной жидкостью в движущейся вязкой жидкости, появляются касательные силы T , которые отсутствуют в идеальной жидко-
сти (рис. 3.2). Кроме этого имеются нормальные силы P. Условимся положительным нормальным напряжением считать напряжение, направленное в сторону противоположную положительному направлению осей координат. Это же правило применяется и для касательных напряжений.
Рассмотрим бесконечно малый объем движущейся жидкости в виде параллелепипеда (рис. 3.3а). Этот же параллелепипед изображен на плоскости на рис. 3.3б. Обозначим напряжения по граням параллелепипеда. В обозначениях нормальных и касательных напряжений, первый индекс обозначает грань, перпендикулярную некоторой оси; второй индекс обозначает проекцию на соответствующую ось. Вдоль каждой грани параллелепипеда дей-
Рис. 3.2. Касательная и нормальная поверхностные силы
2
а)
б)
Рис. 3.3. Элементарный объем движущейся жидкости
ствуют нормальные и касательные напряжения. Разложим последние по направлениям осей.
Напряженное состояние в данной точке жидкости характеризуется девятью напряжениями:
pxx |
τxy |
τxz |
|
τ yx |
pyy |
τ yz |
(3.3) |
τzx |
τzy |
pzz |
|
Т.к. напряжения рассматриваются как непрерывные функции координат, то при изменении координат на бесконечно малые величины, нормальные и касательные напряжения получат бесконечно малые приращения
∂pxx dx, ∂τxy dx и др., как показано на рис. 3.3б.
∂x ∂x
3
Мысленно остановим частицу и добавим к массовым силам силу инерции. Составим сумму моментов всех сил, действующих на эту частицу жидкости, относительно оси, проходящей через центр тяжести C параллелепипеда и параллельную оси z . Приравняем сумму нулю. Массовые силы приложены в центре тяжести частицы и моментов не дают. Моменты от сил, обусловленных нормальными напряжениями, равны нулю. Следовательно,
(τxy + ∂τxy dx)dydz dx 2 + τxy dydz dx 2 − (τ yx
∂x
−τ yx dxdz dy
2 = 0.
После деления на dxdydz
2 , получим
+∂τ yx dy)dxdz dy 2 −
∂y
τxy |
+ |
∂τxy |
dx − τ yx − |
∂τ yx |
dy = 0. |
|
2∂x |
2∂y |
|||||
|
|
|
|
Когда размеры параллелепипеда стремятся к нулю, т.е. dx, dy, dz → 0,
τxy = τ yx , τ yz = τzy , τzx = τxz . |
(3.4) |
Уравнения показывают, что касательные напряжения, действующие на взаимно перпендикулярных площадках одинаковы. Поэтому напряженное состояние в данной точке определяется только шестью напряжениями, а не девятью, как в (3.3).
Согласно второму закону Ньютона сумма проекций поверхностных и массовых сил на ось x равна произведению массы частицы на проекцию ускорения на эту же ось:
∂p |
xx |
dxdydz + |
∂τ yx |
dxdydz + |
∂τ |
zx |
dxdydz + X ρdxdydz = ρdxdydz |
dV |
x |
. |
|
|
∂y |
|
|
|
|
||||
∂x |
|
∂z |
dt |
|||||||
После сокращения на dxdydz , получим уравнение движения вязкой жид-
кости в напряжениях:
|
dV |
x |
|
|
∂p |
xx |
|
|
|
∂τ yx |
|
|
|
∂τ |
zx |
|
|
|
||||||||
ρ |
|
|
= ρX + |
|
|
+ |
|
|
|
+ |
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
dt |
|
|
∂x |
|
|
|
∂y |
|
|
|
∂z |
|
|
|
|
||||||||||
|
dVy |
|
|
∂p yy |
|
|
|
∂τ yy |
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
∂τzy |
|
|||||||||||||||||
ρ |
|
|
|
= ρY + |
|
|
|
|
|
|
+ |
|
|
|
|
+ |
|
|
|
|
|
|
(3.5) |
|||
dt |
|
|
∂x |
|
|
∂y |
|
∂z |
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
dV |
|
|
∂p |
|
|
|
|
|
∂τ yz ∂τ |
|
|
|
|
|
|||||||||||
|
z |
|
zz |
|
|
|
zz |
|
|
|||||||||||||||||
ρ |
|
|
= ρZ + |
|
|
|
+ |
|
|
|
+ |
|
|
|
|
|||||||||||
dt |
|
|
|
∂x |
|
∂y |
|
∂z |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Вторая и третья строки уравнений записаны по аналогии с первой.
Далее вводятся гипотезы о напряжениях. Под влиянием вязкости жидкости меняются нормальные напряжения и появляются касательные напряжения.
4
Нормальные напряжения. Ранее для идеальной невязкой жидкости все три нормальных напряжения были одинаковыми pxx = pyy = pzz = − p , при-
чем давление p называлось гидродинамическим давлением в идеальной
жидкости. В вязкой жидкости появляются добавочные нормальные напряже-
ния p′ |
, |
p′ |
, p′ |
, так что |
|
|
|
|
|
|
|
|
xx |
|
yy |
|
zz |
|
|
|
|
|
|
|
|
|
|
|
p |
xx |
= − p + p′ |
, p |
yy |
= − p + p′ |
, p |
zz |
= − p + p′ . |
(3.6) |
|
|
|
|
xx |
|
yy |
|
zz |
|
|||
Вводятся предположения, что добавочные напряжения пропорциональны скоростям относительного удлинения частицы. Вдоль осей координат скоро-
сти относительного |
удлинения |
равны |
∂Vx / ∂x, ∂Vy / ∂y, ∂Vz / ∂z. |
Поэто- |
||||||||||||||
|
|
∂V |
|
|
∂Vy |
|
|
|
∂V |
|
|
|
|
|
|
|
||
му p′ |
= 2µ |
x |
, p′ |
= 2µ |
|
, p′ |
= 2µ |
|
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
xx |
|
∂x |
yy |
|
∂y |
|
zz |
|
∂z |
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
pxx |
= − p + |
2µ |
|
∂V |
x |
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
∂x |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
= − p + 2µ |
|
∂Vy |
|
||||||||
|
|
|
|
|
|
pyy |
|
|
|
|
|
|
(3.7) |
|||||
|
|
|
|
|
|
|
∂y |
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
= − p + |
2µ |
∂V |
z |
|
|
|
|
||||
|
|
|
|
|
|
pzz |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
∂z |
|
|
|
|
|
Так как для несжимаемой жидкости
∂V |
x |
+ |
∂Vy |
+ |
∂V |
z |
= 0. |
(3.8) |
|
|
∂y |
|
|
||||
∂x |
|
∂z |
|
|||||
то сумма трех нормальных напряжений, действующих по взаимно перпендикулярным площадкам:
pxx + pyy + pzz = −3 p.
Вкачестве гидродинамического давления в вязкой жидкости принимается среднее арифметическое трех нормальный напряжений, действующих по взаимно перпендикулярным площадкам:
p = − |
pxx |
+ pyy |
+ pzz |
. |
(3.9) |
|
3 |
|
|||
|
|
|
|
|
Если повернуть оси координат, как показано на рис. 3.4, то гидродинамическое давление не изменится. Другими словами оно — инвариант относительного поворота осей координат.
5
Рис. 3.4. Гидродинамическое давление не зависит от ориентировки осей коорди- нат
Рис. 3.5. Скорость скошения прямого угла
Касательные напряжения. При движении жидкости только вдоль оси Ox касательные напряжения определяются уравнением (3.1). Они пропорциональны скорости скошения прямого угла. В случае сложного движения (как вдоль оси Ox, так и Oy ) касательные напряжения также считаются
пропорциональными скоростям скошения соответствующих прямых углов (рис. 3.5). Т.е. закон Ньютона (3.1) обобщается на случай сложного движении и предполагается, что
|
|
|
∂V |
x |
|
|
∂Vy |
|
||||||
τxy |
= τ yx |
= µ |
|
|
+ |
|
|
|
|
|
|
|||
∂y |
|
∂x |
|
|
||||||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
∂Vy |
|
|
∂V |
|
|
|
|||||
τ yz |
= τzy |
= µ |
|
|
|
|
+ |
|
z |
|
(3.10) |
|||
∂z |
|
∂y |
||||||||||||
|
|
|
|
|
|
|
|
|||||||
|
|
|
∂V |
z |
|
|
∂V |
|
|
|||||
τzx |
= τxz |
= µ |
|
|
+ |
|
x |
|
|
|
||||
|
|
|
|
∂z |
|
|||||||||
|
|
|
∂x |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
Нормальные напряжения по уравнению (3.7) и касательные по (3.10) подставляем в уравнение (3.5) и получаем:
|
dV |
x |
|
|
|
∂p |
|
|
∂2V |
x |
|
|
|
∂2V |
x |
|
|
|
∂2V |
x |
|
|
|
|
|
∂ |
∂V |
x |
|
|
|
∂Vy |
|
|
|
|
∂V |
z |
|
|
|
|
|||||||||||||||||||||
ρ |
|
|
= ρX − |
|
|
|
+ µ |
|
|
|
|
+ |
|
|
|
+ |
|
|
|
|
|
+ µ |
|
|
|
|
|
|
|
|
+ |
|
|
|
|
+ |
|
|
|
|
|
, |
|||||||||||||||||||||
|
|
|
|
|
|
|
∂x2 |
|
∂y2 |
|
∂z2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||||||
|
dt |
|
|
|
∂x |
|
|
|
|
|
|
|
|
|
|
|
|
∂x |
|
∂x |
|
|
|
|
|
∂y |
|
|
|
|
|
∂z |
|
|
|||||||||||||||||||||||||||||
|
dVy |
|
|
|
∂p |
|
|
∂2Vy |
|
|
|
∂2Vy |
|
|
|
∂2Vy |
|
|
|
|
|
∂ |
|
∂V |
x |
|
|
|
∂Vy |
|
|
|
|
∂V |
z |
|
|
||||||||||||||||||||||||||
ρ |
|
|
|
|
= ρY − |
|
|
|
+ µ |
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
|
|
+ µ |
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
+ |
|
|
|
|
|
, |
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
dt |
|
|
|
∂y |
|
|
∂x |
2 |
|
|
|
|
∂y |
2 |
|
|
|
|
∂z |
2 |
|
|
|
|
|
|
∂y ∂x |
|
|
∂y |
|
|
|
|
|
∂z |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
|
dV |
z |
|
|
|
∂p |
|
∂2V |
z |
|
|
|
∂2V |
z |
|
|
|
∂2V |
z |
|
|
|
∂ ∂V |
x |
|
|
|
∂Vy |
|
|
|
∂V |
z |
|
|
|
|
||||||||||||||||||||||||||
ρ |
|
= ρZ − |
|
|
|
+ µ |
|
|
|
|
+ |
|
|
|
+ |
|
|
+ µ |
|
|
|
|
|
|
|
+ |
|
|
|
|
+ |
|
|
|
|
. |
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||||
|
dt |
|
|
|
|
∂z |
|
|
∂x2 |
|
|
|
|
|
∂y2 |
|
|
|
|
|
∂z2 |
|
|
|
∂z |
∂x |
|
|
|
|
|
∂y |
|
|
|
|
∂z |
|
|
|
|
||||||||||||||||||||||
Так как для несжимаемой жидкости выполняется условие (3.8), то последние слагаемые в полученных уравнениях обращаются в нуль. Окончательно уравнения движения вязкой несжимаемой жидкости запишутся так:
dV |
|
|
|
|
1 ∂p |
|
|
∂2V |
x |
|
|
|
∂2V |
x |
|
|
|
∂2V |
x |
|
|
|||||||||||||||||||||
|
x |
|
= X − |
|
|
|
|
|
|
|
+ ν |
|
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
dt |
|
|
|
|
ρ ∂x |
|
|
∂x |
2 |
|
|
|
|
|
∂y |
2 |
|
|
|
|
|
∂z |
2 |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(3.11) |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂2Vy |
|
|
|
∂2Vy |
|
|
|
∂2Vy |
|
|
||||||||||||||||||
dVy |
|
|
|
1 ∂p |
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||
|
|
|
|
= Y − |
|
|
|
|
|
|
|
+ ν |
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
dt |
|
|
|
ρ ∂y |
|
|
∂x |
2 |
|
|
|
|
|
|
∂y |
2 |
|
|
|
|
|
|
∂z |
2 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
dV |
z |
|
|
|
|
1 ∂p |
∂2V |
z |
|
|
|
∂2V |
z |
|
|
|
∂2V |
z |
|
|
|
|
||||||||||||||||||||
|
= Z − |
|
|
|
|
|
|
|
|
+ ν |
|
|
|
|
|
+ |
|
|
|
|
|
|
+ |
|
|
|
|
|
|
, |
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
dt |
|
|
|
ρ ∂z |
|
|
∂x |
2 |
|
|
|
|
|
|
∂y |
2 |
|
|
|
|
|
|
∂z |
2 |
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
где ν = µ
ρ — кинематическая вязкость жидкости.
Эти уравнения были выведены французским ученым Навье в 1822 г. и независимо от него английским ученым Стоксом в 1845 г. и называются уравнениями Навье-Стокса. Уравнения (3.11) есть уравнения в частных производных второго порядка. В этих уравнениях известны плотность ρ , вязкость
ν и проекции внешних массовых сил X , Y , Z . Неизвестные величины: дав-
ление p и три проекции скорости Vx , Vy , Vz — всего четыре неизвестных.
Число неизвестных превышает число уравнений. Система уравнений незамкнута. К этой системе добавляется еще одно уравнение — уравнение неразрывности (3.8). Если задать граничные условия, то система уравнений (3.10) и (3.8) будет иметь единственное решение. Для неподвижной твердой поверхности граничное условие будет:
Vx = Vy = Vz = 0.
Выше это условие названо условием прилипания.
В частных случаях удается получить точные решения этих уравнений. Пренебрегая теми или иными членами в этих уравнениях, получаем прибли-
7
женные решения. В связи с широким распространением ЭВМ, возможно чис- ленное решение уравнений движения вязкой жидкости. Большое число решенных задач приведено в работе [13].
3.3. Пример точного решения уравнений Навье-Стокса
Рассмотрим поступательное плоское установившиеся течение тяжелой несжимаемой жидкости в канале с параллельными стенками в поле силы тяжести.
Найти потери напора на участке канала между сечениями 1-2 на расстоя-
нии l. p = ρghпт = p1 − p2 .
Уравнения Навье-Стокса для плоского движения:
|
∂V |
|
|
|
∂V |
|
|
|
∂V |
|
1 ∂p |
|
∂2V |
|
|
|
∂2V |
|
|
|
|
|||||||||||||||||||
|
x |
V |
+ |
|
x |
V + |
|
|
x |
= X − |
|
|
|
|
|
|
|
+ ν |
|
|
x |
|
+ |
|
|
|
x |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
||||||||||||||||
|
∂x |
x |
|
|
∂z |
z |
|
∂t |
|
ρ ∂x |
|
∂x |
|
|
|
|
∂z |
2 |
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∂2V |
|
|
∂2V |
|
|
|
|
|||||||||
|
∂V |
|
|
|
∂V |
|
|
|
∂V |
1 ∂p |
|
|
z |
|
|
|||||||||||||||||||||||||
|
z |
V |
|
+ |
|
z |
|
V + |
|
|
z |
= Z − |
|
|
|
|
|
|
|
+ ν |
|
|
|
z |
|
+ |
|
|
|
|
|
(3.12) |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||
|
∂x |
x |
|
|
∂z |
z |
|
∂t |
ρ ∂z |
|
∂x |
2 |
|
|
|
|
∂z |
2 |
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
Граничные условия: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
при z = ±h, Vx = 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
По условию задачи |
∂Vx |
|
= |
∂Vz |
= 0, X = 0, Z = g, V = 0, |
∂Vx |
= 0 |
и сис- |
||||||||||||||||||||||||||||||||
|
|
|
|
|||||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
∂t |
|
|
∂t |
|
|
|
|
|
|
|
|
|
|
|
|
z |
|
|
|
∂x |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
тема уравнений (3.12) существенно упрощается
Рис. 3.6. Плоское установившиеся движения жидкости в канале
8
|
1 ∂p |
|
|
∂2V |
|
|
|||||
0 = − |
|
|
|
|
|
+ ν |
|
x |
|
|
|
ρ ∂x |
∂z |
2 |
|
||||||||
|
|
|
|
|
|||||||
0 = g − |
1 ∂p |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
(3.13) |
|||
ρ ∂z |
|
|
|||||||||
|
|
|
|
|
|
|
|||||
При z = ±h, |
|
|
|
|
|||||||
|
|
|
|
||||||||
V = 0 |
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
x |
|
|
|
|
|
|
|
|
|
|
|
Давление p = p (x, z ), скорость |
Vx = Vx (z ), Vx ≠ Vx (x ). |
Интегрируем |
|||||||||
второе уравнение системы (3.13) по переменной z : |
|
||||||||||
p = ρgz + f (x ), |
|
(3.14) |
|||||||||
где f (x ) — функция интегрирования, зависит только от x. Если дифферен-
цировать уравнение (3.14) по z, то получим второе уравнение системы (3.13). Следовательно, несмотря на движение жидкости, давление в данном сечении x = const подчиняется гидростатическому закону: оно увеличивается пропорционально первой степени z.
Дважды интегрируя первое уравнение системы (3.13). В этом уравнении
|
∂p |
= f ′(x ) |
— функция x и V |
= V |
(z ). Частные производные можно за- |
|||||||||||||||||||||||||
|
|
|||||||||||||||||||||||||||||
|
∂x |
|
|
|
|
|
|
|
|
|
|
x |
|
|
|
x |
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
менить на обыкновенные. Тогда получим |
|
|
|
|
|
|||||||||||||||||||||||||
|
|
|
|
∂Vx |
= |
1 |
|
dp |
|
z + C , V = |
1 |
|
dp |
z |
2 + C z + C |
|
, |
|||||||||||||
|
|
|
|
|
|
|
|
|
|
2 |
||||||||||||||||||||
|
|
|
|
∂z |
µ dx |
1 |
|
|
|
x |
|
|
2µ dx |
|
1 |
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
где µ = νρ, |
а постоянные интегрирования C1 и C2 находим из граничного |
|||||||||||||||||||||||||||||
условия |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
0 = |
1 |
|
dp |
h2 + C h + C |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2µ dx |
1 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
0 = |
1 |
|
dp |
h2 − C h + C |
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
2µ dx |
1 |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
Поэтому C = 0, C |
|
= − |
1 |
|
|
dp |
(z2 − h2 ) и |
|
|
|
|
|
||||||||||||||||
|
|
2 |
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
1 |
|
|
|
|
|
2µ dx |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
V = V = |
1 |
|
dp |
(z2 − h2 ) |
|
|
(3.15) |
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
2µ dx |
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
Распределение скоростей в данном сечении |
x = const |
подчиняется пара- |
||||||||||||||||||||||||||
болическому закону. На стенках Vx = 0, максимум скорости
9
|
|
|
V |
|
= − |
1 |
|
dp |
h2 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
max |
|
|
2µ dx |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||
в середине канала. (Знак минус, т.к. |
dp dx < 0 ). |
|
|
|
|
|||||||||||
Расход жидкости |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ h |
1 dp |
h |
|
|
|
|
|
h |
|
2 |
|
dp |
|
|||
q = ∫ Vdz = |
∫ |
z2 dz − h2 ∫ |
dz = − |
|
h3 |
|||||||||||
|
|
|
|
|
||||||||||||
2µ dx |
|
|
||||||||||||||
− h |
|
|
|
|
|
|
−h |
|
3µ dx |
|||||||
|
|
|
−h |
|
|
|
|
|
|
|
|
|
||||
Средняя скорость
V = |
q |
= − |
1 |
|
dp |
h2 |
|
|
|
||||
ср |
2h |
|
3µ dx |
|||
|
|
|||||
(3.16)
(3.17)
в полтора раза меньше максимальной. Проинтегрируем (3.16) и найдем потерю напора p :
p1 |
3µVср |
x1 |
|
|
∫ dp = − |
∫ dx. |
|
||
2 |
|
|
||
p2 |
h |
x2 |
|
|
|
|
|
||
Поэтому |
|
|
|
|
p1 − p2 = |
p = |
3µVсрl |
(3.18) |
|
|
|
|||
|
h2 |
|||
|
|
|
|
|
Потеря энергии пропорциональна первым степеням вязкости , средней скорости Vср и длине l рассматриваемого участка канала.
Аналогичным способом решается задача о ламинарном движении жидкости в круглой трубе. Для потери напора получаем формулу Пуазейля
p = |
32µVсрl |
. |
(3.19) |
|
|||
|
d 2 |
|
|
3.4. Пример приближенного решения уравнений Навье-Стокса (движение шара в вязкой жидкости)
Задача о медленном движении шара с постоянной скоростью в вязкой несжимаемой жидкости решена Стоксом. Так как рассматривается медленное движение, то в уравнениях (3.11) можно пренебречь инерционными членами: dVx
dt = dVy
dt = dVz
dt = 0. Кроме того, пренебрегая внешними массо-
выми силами. Уравнения (3.11) упрощаются и удается их проинтегрировать. По формуле Стокса действующая на шар сила равна:
Px = 3πµV∞ d , |
(3.20) |
где V∞ — скорость движения,
10
