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

гидромеханика нефти

.pdf
Скачиваний:
30
Добавлен:
02.01.2021
Размер:
20.67 Mб
Скачать

СПБГУАП| Институт 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− 2L1 , получаем

 

 

 

 

 

 

 

 

 

 

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

 

 

12

 

 

формулы ДарсиВейсбаха и Вейсбаха

 

 

 

 

 

 

 

 

 

 

 

 

 

h

=

λ l w2

,

hm =

ζ w2 .

 

 

 

(16.77)

 

 

τ

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d 2g

 

 

 

 

 

2g

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Так как уравнение неразрывности не содержит вязкостных характерис- тик жидкости, то оно имеет одинаковый вид как для линейно-вязкой, так и для любой неньютоновской жидкости. Уравнение Бернулли, представ- ляющее собой закон сохранения механической энергии, также, очевидно, сохраняет свой вид, однако коэффициенты Кориолиса α и величины по- терь h12 будут иными, чем в случае линейно-вязкой жидкости. Действи- тельно, значение величины α определяется законом распределения скорос- тей по сечению трубы, а величина потерь h12 зависит от вязкостных ха- рактеристик среды. Формулы ДарсиВейсбаха и Вейсбаха получены из общих соображений теории размерностей. Поэтому их структура сохраняет- ся и для неньютоновских жидкостей, однако зависимости коэффициента гидравлического сопротивления λ и коэффициента местных сопротивле- ний ζ от критериев подобия будут иметь свой вид для каждого типа неньютоновской жидкости. Из изложенного следует, что все схемы расче- та трубопроводов, опирающиеся на соотношения (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