
гидромеханика нефти
.pdfСПБГУАП| Институт 4 группа 4736
ЛАМИНАРНЫЕ ТЕЧЕНИЯ НЕНЬЮТОНОВСКИХ ЖИДКОСТЕЙ |
325 |
Из формулы (16.49) видно, что для определения констант η |
и τ 0 не- |
|||||||||||||||||||||||||||||||||||||||||
обходимо измерить две пары значений M, Ω . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
III. Степенная жидкость. Для степенной жидкости в соответствии с фор- |
||||||||||||||||||||||||||||||||||||||||||
мулой (16.17) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
f(τ ) = |
τ |
|
1 |
, s(τ a) = 0 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
γ = |
|
n |
|
|
|
|
|
|
|
|
|
|
(16.50) |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
k |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Подставив выражение (16.50) в формулы (16.23) и (16.24), с учетом |
||||||||||||||||||||||||||||||||||||||||||
равенства (16.20), получим, соответственно, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
1+ n |
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
1 |
|
|
|
|
n+ 1 |
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
v(r) = |
n |
|
|
τ a n |
|
τ n |
|
|
|
|
n |
|
|
+ 1 |
|
∆ |
p n |
r n |
|
|||||||||||||||||||||||
|
|
= |
|
|
|
|
n |
|
||||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 − |
|
|
|
|
|
|
|||
|
|
a |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
a |
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
n + |
1 |
|
k |
|
|
|
τ a |
|
|
|
|
|
|
|
n + 1 |
|
|
|
2kl |
|
|
|
a |
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
τ a |
|
1 |
|
|
|
|
|
|
|
|
|
3n+ 1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
||||||
|
|
Q = |
|
|
n |
3 |
n |
= |
|
|
|
n |
|
|
|
∆ p n |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
π a |
|
|
|
|
|
|
|
|
π a |
|
|
|
|
|
|
|
, |
|
|
|
(16.51) |
|||||||||||||
|
|
|
|
|
+ 1 |
|
|
|
3n + 1 |
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
|
|
|
|
|
3n |
|
k |
|
|
|
|
|
|
|
|
|
|
2kl |
|
|
|
|
|
|
то есть получим формулы для распределения скоростей и расхода при те- чении по круглой трубе.
Для ротационного вискозиметра в соответствии с формулами (16.30), (16.34), (16.35) имеем
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|||
|
v n τ i |
|
τ |
|
n |
n |
|
M |
|
n |
|
|||||||||||||||
|
n |
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
1 − |
|
|
|
|
|
= |
|
|
|
|
|
|
|
− |
||
= |
|
|
|
|
|
|
|
τ i |
|
|
|
2 |
|
|
2 |
1 |
||||||||||
|
r 2 |
k |
|
|
|
|
|
|
|
|
|
|
2π kRi |
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
|
|
|
1 |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
n τ i |
n |
|
|
τ e |
|
n |
n |
M |
|
n |
|
||||||||||||
Ω = |
|
|
|
|
|
|
1 − |
|
|
|
|
|
= |
|
|
|
|
|
|
1 − |
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
τ i |
|
|
|
|
|
|
2 |
||||||||||||||
|
|
|
2 |
k |
|
|
|
|
|
|
|
|
2 |
|
|
2π kRi |
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
||
|
Ri |
n |
|
|
||||
|
|
|
|
, |
||||
|
||||||||
|
r |
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
||
|
Ri |
|
n |
|||||
|
|
|
|
|
. |
|||
|
|
|
||||||
|
|
|
|
|
||||
|
Re |
|
|
|||||
|
|
|
|
|
|
|
|
IV. Ряд Рейнера. Предполагается, что для жидкости с начальным на- пряжением сдвига τ 0 функция f (τ − τ 0) может быть разложена в степен- ной ряд. Так как эта функция нечетная, то ряд может содержать только не-
четные степени τ |
− |
τ 0 . Следовательно, |
|
|
||||||
|
|
|
|
|
0 , |
|
τ ≤ |
τ 0 , |
|
|
γ |
= |
f(τ − |
τ |
0 ) = |
|
|
|
|
|
|
n |
|
bk (τ |
− τ ) 2k+ 1, |
τ ≥ τ , s( τ a) = 0 , |
(16.52) |
|||||
|
|
|
|
|
∑ |
|
0 |
0 |
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
k= |
0 |
|
|
|
|
где bk, τ 0 |
– |
реологические параметры жидкости. Подставив выраже- |
ние (16.52) в формулу (16.23), после рассуждений, совершенно аналогич-
Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736
326 ГЛАВА XVI
ных приведенным при рассмотрении жидкости Бингама–Шведова, получа- ем
v(r) |
|
|
a |
n |
|
b |
|
|
(τ a − |
τ 0) 2 |
k |
+ 2 = |
|
|
a |
n |
|
|
b |
|
|
τ a2k+ |
|
|
|
|
r |
2k+ 2 |
|
|
|
|
|||||||||
= |
|
|
|
|
∑ |
k |
|
|
|
|
|
|
∑ |
|
|
k |
|
|
1 1 |
− |
|
|
0 |
|
|
= |
const, |
0 ≤ r ≤ r0, |
|||||||||||||
2τ a |
k + |
1 |
|
|
2 |
|
k + |
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
k= |
0 |
|
|
|
|
|
|
|
k= 0 |
|
1 |
|
|
|
|
|
|
|
a |
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
v(r) = |
a |
|
n |
|
bk |
|
|
[(τ a − τ 0)2k+ 2 − ( τ − τ |
|
0) 2k+ 2] = |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
∑ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
2τ a |
k + |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||
|
|
k= |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
a |
n |
b |
|
|
|
|
|
|
r |
|
2k+ 2 |
|
|
r |
|
2k+ |
2 |
|
|
|
|
r |
|
2k+ 2 |
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
2k+ |
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
k |
|
|
|
|
0 |
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
r |
r |
a |
||||||
|
|
= |
|
|
|
∑k + |
|
|
τ a |
− |
|
|
|
|
|
|
− |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
2 |
|
|
1 |
|
|
|
|
|
|
|
|
|
|
1 − |
|
|
|
|
|
|
|
, |
0 ≤ |
|
≤ . |
||||||||||||||
|
|
|
1 |
|
|
|
a |
|
|
|
a |
|
|
|
|
|
|
|
r0 |
|
|
|
|
|
|||||||||||||||||
|
|
|
|
|
|
k= |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Из полученных формул видно, что картина распределения скоростей по радиусу качественно аналогична представленной на рис. 16.10, то есть и в этом случае существует ядро течения радиуса r0 .
Подставив выражение (16.52) в формулу (16.24), получим
|
π a |
3 |
n |
bk (τ a − |
|
0)2k+ 2 |
|
(τ a |
− |
τ 0) |
2 |
|
2τ 0 (τ a − τ 0) |
|
τ |
2 |
|
|
Q = |
∑ |
τ |
+ |
+ |
|
|||||||||||||
|
|
|
0 |
|
. (16.53) |
|||||||||||||
τ a3 |
|
|
4 |
|
|
2k + |
|
|||||||||||
|
|
|
|
2k + |
|
|
2k + 3 |
2 |
||||||||||
|
|
|
k= 0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
При течении в кольцевом зазоре так же, как в случае жидкости Бинга-
ма–Шведова, |
при τ i < |
τ 0 |
течения не происходит. При M > |
M0 = |
2π Ri2τ 0 |
||||||||||||||||
сдвиговое течение происходит в зазоре |
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Ri < |
r < r0 = |
|
M |
, |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2πτ |
0 |
|
|
а при r0 |
< r < |
Re |
|
жидкость вращается с постоянной угловой скоростью, то |
|||||||||||||||||
есть как твердое тело. При M ≥ |
2π Re2τ 0 |
сдвиговое течение охватывает всю |
|||||||||||||||||||
область |
Ri |
< |
r < |
|
Re . В соответствии с формулами (16.34) и (16.49) при |
||||||||||||||||
τ i > τ > |
τ e |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
v |
|
|
1 |
|
n |
|
τ i (τ − τ 0) 2k+ 1 |
|
|
|
||||||
|
|
|
|
|
|
= |
|
|
|
∑bk ∫ |
|
|
|
|
dτ , |
Ri < r < r0, |
|
|
|||
|
|
|
|
|
r |
2 |
|
|
τ |
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
k= 0 |
|
τ |
|
|
|
|
|
|
|
|
|
v |
|
|
|
1 |
|
n |
|
τ i (τ − τ 0 ) 2k + 1 |
|
|
|
|
|
|||||||
|
|
= |
Ω |
= |
|
|
∑bk ∫ |
|
|
|
dτ = const, r0 < r < |
Re . |
(16.54) |
||||||||
|
r |
2 |
|
τ |
|||||||||||||||||
|
|
|
|
|
|
|
|
k = |
0 |
τ 0 |
|
|
|
|
|
(Ri < |
r < Re ) |
||||
При τ e |
> |
τ 0 |
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
распределение скоростей во всем зазоре |
будет определяться формулой (16.54), где интеграл необходимо брать в пре- делах τ 0,τ i , а угловая скорость вращения внешнего цилиндра Ω будет
Контакты | https://new.guap.ru/i03/contacts
СПБГУАП| Институт 4 группа 4736
ЛАМИНАРНЫЕ ТЕЧЕНИЯ НЕНЬЮТОНОВСКИХ ЖИДКОСТЕЙ |
327 |
равна
|
1 |
n |
τ i |
(τ |
− τ 0) 2k+ 1 |
|
|
|
Ω = |
|
∑bk ∫ |
|
|
|
dτ . |
(16.55) |
|
2 |
|
τ |
|
|||||
|
|
k= 0 |
τ e |
|
|
|
|
|
Реологическое уравнение (16.52) содержит |
n + |
2 реологических па- |
раметра: τ 0 , b0 , b1, ..., bn . Из формул (16.53) и (16.55) видно, что для их определения необходимо сделать n + 2 измерения пар значений ∆ p, Q или M, Ω .
§7. Коэффициент гидравлического сопротивления
Рассмотрим жидкость с реологическим уравнением |
|
γ = f(τ ,α 1,α 2 , ...,α n ) , |
(16.56) |
где, как и раньше, α 1,α 2 , ...,α n – реологические параметры. По аналогии с соображениями, приведенными при выводе формулы Дарси–Вейсбаха
(5.30), можно утверждать, |
что перепад давления ∆ p |
на длине l в трубе |
диаметром d представляется соотношением вида |
|
|
∆ p = |
ϕ (l, d, ρ , w,α 1,α 2 , ...,α n ) . |
(16.57) |
Приняв величины d, ρ , w в качестве параметров с независимыми размерностями, после рассуждений, тождественных использованным при получении формулы (5.30), из формулы (16.57) получаем
∆ |
p = λ |
l ρ w2 |
|
|||||
|
|
|
, |
|
||||
d |
2 |
|
||||||
|
|
|
|
|
|
|||
где |
|
|
|
|
|
|
|
|
λ = λ (Π 1, Π 2 , ..., Π n ) , |
(16.58) |
|||||||
причем величины |
|
|
α |
|
|
|
|
|
Π |
i = |
|
i |
|
|
|||
d |
β ρ |
γ wδ |
|
|||||
|
|
|
представляют собой критерии подобия. Из формул (16.56) и (16.58) следу- ет, что число критериев подобия равно числу реологических параметров жидкости.
Рассмотрим в качестве примера жидкость Бингама–Шведова. В этом
случае формула (16.57) принимает вид
∆ p = ϕ (l, d, ρ ,η ,τ 0, w) ,
а выражение (16.58) – |
|
λ = λ (Π 1, Π 2) , |
(16.59) |
Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736
328 |
|
|
|
|
|
ГЛАВА XVI |
где |
ρ wd |
|
|
τ 0 |
|
|
Π 1 = |
|
Π 2 = |
. |
(16.60) |
||
η |
, |
ρ w2 |
||||
|
|
|
|
|
Для получения аналитического вида зависимости (16.59) рассмотрим формулу (16.47) . С помощью равенства (16.20) ее можно представить в виде
2 |
|
π a4 |
∆ p |
|
|
4 2l τ 0 |
|
|
1 |
2l |
|
τ 0 |
4 |
|
||||||||||||||||||
Q = π a w = |
|
|
|
1 − |
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
(16.61) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
8η |
l |
|
|
3 a ∆ p |
|
|
|
|
|
|
a |
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
3 |
|
∆ p |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d2 |
∆ p |
|
4 4l |
|
|
τ 0 |
|
|
|
|
1 |
4l τ 0 |
|
4 |
|
|
|
||||||||||||||
w = |
|
|
1 − |
|
|
|
|
|
|
|
|
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
(16.62) |
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
32η l |
|
3 d ∆ p |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
3 |
d ∆ p |
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Очевидно, что для получения соотношения вида (16.37) необходимо |
||||||||||||||||||||||||||||||||
соотношение (16.62) разрешить относительно ∆ |
p . Положим |
|
||||||||||||||||||||||||||||||
|
|
|
|
∆ p = |
4lτ 0 |
z |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(16.63) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
d |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
и подставим это выражение в формулу (16.62). Тогда после элементарных преобразований получим
|
|
|
|
|
|
|
|
z4 − |
4 |
α z3 + |
|
1 |
= |
|
0 , |
|
|
|
|
|
(16.64) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
τ 0d . |
|
|
|
|
|
||
|
|
|
|
|
|
|
α |
= |
1 + |
|
6 |
|
|
|
|
A = |
|
|
|
|
(16.65) |
||||||||||
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
A |
|
|
|
|
|
|
|
|
η w |
|
|
|
|
|
|
|||
Используя стандартную методику решения уравнений четвертой сте- |
|||||||||||||||||||||||||||||||
пени, получаем выражение для корней уравнения (16.64) |
|
||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
= |
c |
|
± |
|
|
1 − |
|
|
|
|
3b |
|
|
|
|
|
(16.66) |
||||||
|
|
|
|
|
z1,2 |
|
|
1 |
|
|
c(c − |
α ) |
, |
|
|
|
|
||||||||||||||
|
|
|
|
|
3 |
|
|
|
|
|
|
||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
z3,4 = |
− |
|
|
c − |
2α |
|
± |
|
|
1 − |
|
|
|
|
|
|
3b |
|
|
(16.67) |
|||||||||||
|
|
|
|
|
|
1 |
|
|
(c − α )( c − |
|
, |
||||||||||||||||||||
|
|
|
3 |
|
|
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2α ) |
|
|||||||||||||
где |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||
|
b = 3 α 2 + |
α 4 − 1 + 3 α 2 − |
α 4 − 1 = |
|
|||||||||||||||||||||||||||
|
|
2 |
3 |
|
|
|
|
|
|
− 4 + |
3 |
|
|
|
|
|
|
|
|
|
|
2 |
|
(16.68) |
|||||||
= α |
|
|
+ |
1 − α |
1 |
− |
1 − α |
= α 3 β , |
|
||||||||||||||||||||||
|
3 |
1 |
|
|
− 4 |
|
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c = |
α |
+ |
|
|
3 b + |
α |
2 . |
|
|
|
|
|
(16.69) |
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736
ЛАМИНАРНЫЕ ТЕЧЕНИЯ НЕНЬЮТОНОВСКИХ ЖИДКОСТЕЙ |
|
|
329 |
||||||||||||||||||||||||
Рассмотрим подкоренное выражение в формуле (16.67). В соответст- |
|||||||||||||||||||||||||||
вии с равенством (16.69) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
3b = 2(c − α ) 2 − 2α 2 , |
|
|
|
|
|
|
|
|
|||||||||||||
то после элементарных преобразований имеем |
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||
|
|
|
1 − |
|
|
|
|
3b |
|
|
|
= − |
c + α |
. |
|
|
|
|
|
||||||||
|
|
|
(c − α )( c − 2α ) |
|
|
|
|
|
|
||||||||||||||||||
|
|
|
c − α |
|
|
|
|
|
|||||||||||||||||||
Так как в соответствии с формулой (16.65) |
|
α |
|
> 1, то из равенств (16.68) |
|||||||||||||||||||||||
и (16.69) следует, что b и с величины вещественные, причем b > |
0, c > α . |
||||||||||||||||||||||||||
Таким образом, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c + |
|
α |
|
|
|
|
|
|
|
|
|
|
|
1 − |
|
|
|
|
3b |
|
|
= − |
|
|
< 0 , |
|
|
|
|
||||||||||
|
|
(c − α )( c − 2α ) |
|
|
|
|
|
|
|||||||||||||||||||
|
|
c − α |
|
|
|
|
|
||||||||||||||||||||
и корни z3,4 – комплексные. |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
Перейдем к рассмотрению корней |
|
z1,2 . Непосредственной проверкой, |
|||||||||||||||||||||||||
используя равенство (16.68), можно убедиться, что |
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
b3 − 3b − 2α 2 = 0 , |
|
|
|
|
|
|
|
|
|
|||||||||||
и из формулы (16.69) имеем |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
Ò = α |
+ |
|
b3 |
|
|
|
|
|
|
|
|
|
|
|
(16.70) |
|||
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставив это выражение в формулу (16.66), получаем |
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z1,2 |
= |
|
|
|
1 − |
|
|
|
|
|
|
|
|
|
|
|
(16.71) |
|||||
|
|
|
|
|
3 |
1 ± |
2b |
. |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
|
|
|
|
|||||
Из формул (16.68) и (16.69) следует, что при α |
= |
1 |
b= |
2, c= |
3 и |
||||||||||||||||||||||
db |
= |
2α |
3 |
α |
2 |
+ |
α |
4 |
− 1 |
− |
3 |
α |
2 |
− |
α |
4 |
− |
|
|
||||||||
dα |
|
α |
|
|
|
|
|
|
|
|
|
1 , |
|
||||||||||||||
|
3 |
4 − 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
откуда ddbα > 0 |
при α |
> |
|
1. Таким образом, функции b(α ) |
|
и c(α ) |
монотон- |
||||||||||||||||||||
но возрастают с ростом α и |
|
|
6 |
|
< |
1 . |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
c |
2b |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
Итак, корни z1,2 |
– вещественные. Для их дальнейшего анализа пере- |
пишем, используя формулы (16.63), (16.65), (16.68) и (16.70), соотноше-
ние (16.71) в виде
∆ pd |
|
1 |
|
6η w |
|
|
|
β |
3 |
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
= |
τ 0 |
+ |
|
|
1 |
+ |
|
|
|
1 |
± 1 |
− |
|
|
. |
(16.72) |
4l |
|
3 |
|
d |
|
|
2 |
|
|
|
c |
|
|
|
|||
|
|
|
|
|
|
|
|
|
2b |
|
Контакты | https://new.guap.ru/i03/contacts

СПБГУАП| Институт 4 группа 4736
330 |
|
|
|
|
ГЛАВА XVI |
Переходя в равенстве (16.72) к пределу при τ 0 → |
0 , получим |
||||
|
∆ pd |
= 4 |
η w |
(1 ± 1) . |
|
|
4l |
|
|
||
|
|
d |
|
Из выражения (16.61) видно, что этот предельный переход должен привести к формуле Пуазейля. Следовательно, в формулах (16.71) и (16.72) необходимо выбрать знак «+» и окончательно
|
|
|
|
c |
|
|
|
|
|
|
6 |
|
|
|
|
|
|
z = |
|
|
|
|
|
1 |
− |
|
|
, |
|
|
|
|
|
3 |
1 + |
2b |
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
c |
|
|
|
|
||
или, с учетом равенства (16.63), |
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
4 l |
|
|
|
|
|
|
6 |
|
|
|
||
|
|
∆ p = |
|
|
τ |
|
|
+ |
1 |
− |
|
|
|
(16.73) |
|
|
|
3 d |
0c 1 |
2b |
. |
||||||||||
|
|
|
|
|
|
|
|
c |
|
|
|||||
Как следует из формул (16.65), (16.68) и (16.69), |
b = b(A), c = c( A) . |
||||||||||||||
Сравнивая выражение (16.73) с формулой Дарси–Вейсбаха, получаем |
|||||||||||||||
|
|
|
|
λ |
= |
|
8 |
Bϕ (A) |
, |
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
||
где B = |
τ 0 |
– безразмерный параметр, |
|
|
|
|
|
||||||||
ρ w2 |
|
|
|
|
|
||||||||||
|
|
ϕ (A) = |
|
|
|
|
6 |
|
|
|
|
||||
|
|
|
|
|
1 − |
|
|
|
|
|
|||||
|
|
c 1 + |
|
c 2b |
. |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
Таким образом, коэффициент гидравлического сопротивления при течении жидкости Бингама–Шведова есть функция двух независимых критериев подобия А и В, причем В совпадает с Π 2 в формуле (16.60),
а A = Π 1Π 2 . |
|
|
|
||
Численные значения функции ϕ (A) приведены в таблице. Можно по- |
|||||
казать, что при |
1 |
= |
η w ≥ 0,1 функция ϕ (A) может быть аппроксимиро- |
||
|
|||||
|
A |
τ 0d |
|
|
|
вана с погрешностью менее 2% выражением |
|||||
|
|
|
|
6 |
|
|
|
|
ϕ (A) = 4 1 + |
|
. |
|
|
|
|
||
|
|
|
|
A |
Контакты | https://new.guap.ru/i03/contacts
СПБГУАП| Институт 4 группа 4736
ЛАМИНАРНЫЕ ТЕЧЕНИЯ НЕНЬЮТОНОВСКИХ ЖИДКОСТЕЙ |
|
331 |
|||||||||
|
|
|
|
|
|
|
|
|
|
|
Таблица |
1/А |
ϕ (A) |
|
1/А |
ϕ (A) |
|
1/А |
ϕ (A) |
|
|
1/А |
ϕ (A) |
|
|
|
|||||||||
0,0000 |
3,00 |
|
0,0060 |
3,53 |
|
0,0250 |
4,25 |
|
|
0,0700 |
5,52 |
0,0005 |
3,14 |
|
0,0080 |
3,63 |
|
0,0300 |
4,40 |
|
|
0,0800 |
5,78 |
0,0010 |
3,20 |
|
0,0100 |
3,71 |
|
0,0350 |
4,55 |
|
|
0,1000 |
6,29 |
0,0020 |
3,29 |
|
0,0120 |
3,79 |
|
0,0400 |
4.70 |
|
|
0,1500 |
7,54 |
0,0030 |
3,36 |
|
0,0140 |
3,87 |
|
0,0450 |
4,84 |
|
|
0,2000 |
8,76 |
0,0040 |
3,42 |
|
0,0160 |
3,94 |
|
0,0500 |
4,98 |
|
|
0,2500 |
9,97 |
0,0050 |
3,48 |
|
0,0200 |
4,08 |
|
0,0600 |
5,25 |
|
|
0,3000 |
11,18 |
В качестве следующего примера рассмотрим степенную жидкость. Для такой жидкости формула (16.54) принимает вид
∆ p = ϕ (l, d, ρ , k, n, w) .
Приняв в качестве параметров с независимыми размерностями вели- чины d, ρ , w , используя Π -теорему и учитывая, что в соответствии с фор-
мулой (16.17) [k] = MTn− 2L− 1 , получаем |
|
|
|
|
|
|
|||||||
|
|
|
|
dn |
ρ w2− n |
|
l 2 |
||||||
∆ p = |
|
|
|
|
|
|
|
|
|
|
|
|
ρ w , |
f n, |
|
k |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
d |
|
|||||
откуда |
|
|
|
|
|
|
|
|
|
|
|
|
|
λ |
= |
|
|
|
dnρ w2− n |
||||||||
|
|
|
|
|
|
|
|
|
|
|
|||
2f n, |
|
k |
|
|
|
|
|
. |
|||||
|
|
|
|
|
|
|
|
|
|
|
|||
Безразмерными критериями подобия являются величины |
|||||||||||||
n, |
|
dnρ w2− n |
= |
|
Re′, |
||||||||
|
|
|
|
|
|
|
k
где Re′ – аналог числа Рейнольдса для линейно-вязкой жидкости. Для вы- яснения вида зависимости λ (n, Re′) рассмотрим выражение (16.51) или
|
|
|
|
|
|
n+ 1 |
|
∆ |
1 |
|||
w = |
n |
|
d |
n |
|
|
p |
n |
||||
|
|
|
|
|
|
|
|
|
||||
3n + |
|
|
|
|
||||||||
|
1 |
2 |
|
2kl |
Разрешив соотношение (16.74) относительно
∆ p = |
2klw |
n 3n + |
1 n |
2 n+ |
||||
|
|
|
|
|
|
|
||
|
|
|
||||||
|
|
|
n |
|
|
d |
. (16.74)
∆ p , получим
1
.
Контакты | https://new.guap.ru/i03/contacts
СПБГУАП| Институт 4 группа 4736
332 ГЛАВА XVI
Сравнивая это выражение с формулой Дарси–Вейсбаха, получаем
λ = |
|
n |
+ 3 |
3n + |
1 n |
k |
|||
2 |
|
|
|
|
|
|
. |
||
|
|
|
|
|
|||||
|
|
|
|
|
n |
|
|
ρ w2− nd |
§8. Дополнительные замечания к расчету течения неньютоновских жидкостей по трубам
Основными соотношениями, описывающими установившееся движе- ние вязкой жидкости по трубам являются:
уравнение неразрывности
|
|
|
|
Q = |
|
wS = |
const , |
|
|
|
|
|
|
|
|
(16.75) |
||||||
уравнение Бернулли |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
z + |
p1 |
+ |
α |
|
w12 |
|
= |
z + |
p2 |
+ |
α |
|
|
w22 |
|
+ |
h |
, |
(16.76) |
|||
|
1 2g |
|
2 2g |
|||||||||||||||||||
1 |
ρ g |
|
|
|
2 |
ρ g |
|
|
1− 2 |
|
|
|||||||||||
формулы Дарси–Вейсбаха и Вейсбаха |
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
h |
= |
λ l w2 |
, |
hm = |
ζ w2 . |
|
|
|
(16.77) |
||||||||||||
|
|
τ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
d 2g |
|
|
|
|
|
2g |
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Так как уравнение неразрывности не содержит вязкостных характерис- тик жидкости, то оно имеет одинаковый вид как для линейно-вязкой, так и для любой неньютоновской жидкости. Уравнение Бернулли, представ- ляющее собой закон сохранения механической энергии, также, очевидно, сохраняет свой вид, однако коэффициенты Кориолиса α и величины по- терь h1− 2 будут иными, чем в случае линейно-вязкой жидкости. Действи- тельно, значение величины α определяется законом распределения скорос- тей по сечению трубы, а величина потерь h1− 2 зависит от вязкостных ха- рактеристик среды. Формулы Дарси–Вейсбаха и Вейсбаха получены из общих соображений теории размерностей. Поэтому их структура сохраняет- ся и для неньютоновских жидкостей, однако зависимости коэффициента гидравлического сопротивления λ и коэффициента местных сопротивле- ний ζ от критериев подобия будут иметь свой вид для каждого типа неньютоновской жидкости. Из изложенного следует, что все схемы расче- та трубопроводов, опирающиеся на соотношения (16.75) – (16.77), с уче- том приведенных выше оговорок могут быть использованы для расчета течения неньютоновских вязких жидкостей.
Контакты | https://new.guap.ru/i03/contacts
СПБГУАП| Институт 4 группа 4736
Глава XVII
ДВУХФАЗНОЕ ТЕЧЕНИЕ В ТРУБАХ
С течением двухфазных (многофазных) сред в трубах приходится сталкиваться почти во всех отраслях нефтегазового производства. При бу- рении скважин – это течение аэрированных промывочных и тампонажных жидкостей и вынос шлама. При эксплуатации нефтяных и газовых место- рождений – эрлифтная добыча нефти, течение газоконденсатных, водо- нефтяных и газоводонефтяных смесей в стволе скважины. В сборных сетях месторождений также могут иметь место многофазные течения. Этот спи- сок может быть существенно продолжен.
Под фазой подразумевается отдельная часть однородной системы, ог- раниченная поверхностью раздела. Так, смесь нефти и воды представляет собой двухфазную систему жидкость-жидкость. Смесь газа и конденсата, или газа и нефти – двухфазная система газ–жидкость. Смесь воды, нефти и газа – трехфазная система.
Фаза может состоять из одного вещества, например, воды. Такая фаза называется однокомпонентной. Если фаза состоит из нескольких химиче- ских веществ, например, смеси углеводородных газов, то она называется многокомпонентной.
Истинные растворы (соли в воде, смеси газов и т.д.) представляют со- бой однофазные многокомпонентные системы.
При описании движения многофазных сред обычно вводятся следую- щие предположения.
1.Размеры включений или неоднородностей в смеси (отдельных час- тей неоднородной системы) много больше расстояний между молекулами, длин свободного пробега молекул и т.д. Иначе говоря, размеры включений таковы, что к каждой отдельной части неоднородной системы приложимы методы механики сплошной среды.
2.Размеры указанных включений много меньше расстояний, на кото- рых макроскопические параметры смеси или фаз меняются существенным образом, то есть эти размеры много меньше характерных размеров рас- сматриваемой системы.
Указанные ограничения позволяют использовать для описания дви- жения многофазных сред модель многоскоростного континуума. Много-
Контакты | https://new.guap.ru/i03/contacts
СПБГУАП| Институт 4 группа 4736
334 |
ГЛАВА XVII |
скоростной континуум представляет собой совокупность N континуумов, каждый из которых относится к своей составляющей (фазе или компонен- те) и заполняет один и тот же объем, занятый смесью. Таким образом, в каждой точке многоскоростного континуума имеется N плотностей, N скоростей и т.д.
§1. Уравнения законов сохранения
Исходя из общих принципов, использованных в главе II для получе- ния уравнений законов сохранения для однофазной среды, и модели мно- госкоростного континуума, можно выписать уравнения законов сохране- ния массы, изменения количества движения и энергии для каждой состав- ляющей смеси.
В интегральном представлении эти законы имеют вид: закон сохранения массы
|
|
|
|
∂ (α i ρ i ) |
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
||||
|
|
|
∫ |
|
dV + ∫α i ρ ivindS = ∫∑Jji dV, i = |
1,2, ..., N ; |
(17.1) |
||||||||||||||||
|
|
|
∂ t |
||||||||||||||||||||
|
|
|
V |
|
|
|
|
|
|
S |
|
|
V |
j = 1 |
|
|
|
|
|
|
|||
закон изменения количества движения |
|
|
|
|
|
|
|
|
|
||||||||||||||
|
∂ (α ρ v ) |
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
∫ |
|
|
i i |
i |
dV + |
∫ |
α ρ v v dS= |
α ρ |
|
FdV+ |
∫ |
p dS+ |
∫ |
∑ P dV, |
|
||||||||
|
|
|
|
|
|
|
|||||||||||||||||
|
|
∂ t |
|
|
|
|
|
i i i in |
∫ |
|
i i i |
ni |
|
ji |
(17.2) |
||||||||
V |
|
|
|
|
|
|
|
S |
|
V |
|
|
|
|
|
S |
|
V j= |
1 |
||||
i = |
1, 2, ... ,N; |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
закон сохранения энергии |
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
∂ (α i ρ i Ei ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
∫ |
|
|
|
|
|
dV + ∫α i ρ i EivindS = ∫α i ρ i FividV + |
∫ pnivi dS + |
|
|||||||||||||
|
|
|
|
∂ t |
|
|
|
||||||||||||||||
|
|
V |
|
|
|
|
|
|
|
|
S |
|
|
|
|
V |
|
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ |
|
∫∑Eji dV − |
∫qi(n) dS, |
i = |
1,2, ..., N, |
|
|
(17.3) |
||||||||
|
|
|
|
|
|
|
|
|
|
V |
j = 1 |
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Ei |
= |
ui + |
|
vi2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
В уравнениях (17.1) – (17.3) i |
– номер фазы (компоненты), α i |
≥ 0 – |
доля объема смеси, занимаемая фазой в данной точке, остальные обозна- чения имеют тот же смысл, что в гл. II. При этом, очевидно,
N |
|
∑α i = 1. |
(17.4) |
i= 1 |
|
Контакты | https://new.guap.ru/i03/contacts