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

книги из ГПНТБ / Романков, П. Г. Гидромеханические процессы химической технологии

.pdf
Скачиваний:
29
Добавлен:
22.10.2023
Размер:
10.97 Mб
Скачать

или с учетом симплекса геометрического подобия

 

 

 

 

 

 

Еи = СТ1

 

 

 

 

(3-142)

где Ti — 1/d;

С — коэффициент, определяемый

экспериментально.

Из уравнения Дарси — Вейсбаха

(3-51)

 

 

 

 

 

 

 

 

 

 

Ы

w2

 

 

 

 

получим

 

 

 

 

= ТГР

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3-143)

т, е. значение

коэффициента С Х/2.

Опытным путем

при

Re =

= 250 000 -f- 300 000

в условиях

средней относительной шерохо­

ватости труб найдено зна-

Ей

 

 

 

 

 

 

 

чение

коэффициента

тре­

 

 

 

 

 

 

 

ния X = 0,0286

[6].

 

А

 

 

 

 

 

 

 

Таким

образом,

для

 

* con st

 

 

 

 

у

La

 

 

 

 

второй автомодельной об-

\

 

 

 

 

 

 

 

ласти

 

 

 

 

0,046 - -V -

 

 

 

 

 

 

 

Еи = 0,0143Гг

(3-144)

 

\

 

 

1

 

Eu-const

или

 

 

 

 

 

0.СИ4 ___Х-!

 

1

 

 

 

 

 

 

 

 

 

F

 

 

 

 

 

 

 

 

ßl

 

 

 

До = 0 ,0 1 4 3 -^ --4 -

(3-145)

 

 

 

 

1

 

_______ , „

______і_____________I1

 

 

 

 

 

 

 

 

 

2300

>250000

Re

т. е. в этом случае ско­

Рис. 3-31.

Зависимость

Eu =

/(Re) для дви

рость

потока

 

изменяется

 

жения однофазных потоков.

 

строго по

квадратичному

 

 

 

 

 

 

 

 

закону.

 

 

 

 

 

 

 

 

 

 

 

 

В

гидромеханических

 

 

 

 

 

 

 

 

процессах

закон сопро­

 

 

 

 

 

 

 

 

тивления в форме уравне­

 

 

 

 

 

 

 

 

ний (3-139),

(3-144)

или

 

 

 

 

 

 

 

 

в общем виде (3-134) ши­

 

 

 

 

 

 

 

 

роко используется для оп­

 

 

 

 

 

 

 

 

ределения Ар или w. При

 

 

 

 

 

 

 

 

этом

необходимо

знать

 

 

 

 

 

 

 

 

зависимость X от Re.

 

 

 

 

 

 

 

 

 

На рис. 3-31 и

3-32

Рис. 3-32.

Зависимость X = f'(Re)

для

движе­

приведены

 

характерные

 

 

ния однофазных потоков.

 

для

движения

однофаз­

Eu — f(Re)

и X = f'(Re).

При

малых

ных

потоков

зависимости

значениях скорости на движение потока сильно влияют силы вяз­

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

зависит

от

величины критерия Re

в первой степени (участок АВ

на рис.

3-31

и 3-32). Этот закон ха­

рактерен для ламинарного режима движения потоков в трубах и каналах, а также для движения мелкодисперсных твердых мате­ риалов в потоке при осаждении (см. гл. 4).

Увеличение влияния инерционных сил и уменьшение влияния молекулярной вязкости характерны для перехода от ламинарного

37

режима движения к турбулентному и коэффициент трения X зави­ сит от Re в степени меньше 1, причем чем меньше степень у Re, тем больше турбулентность потока (участок ВС на рис. 3-32 и CD на рис. 3-31). При уменьшении степени у Re до 0 режим движения перестает зависеть от действия сил вязкости (Re = const) и насту­ пает режим развитой турбулентности (автомодельный), при кото­ ром коэффициент X не зависит от Re (CD на рис. 3-32 и EF на рис. 3-31). Такой режим характерен для движения потоков с боль­ шой скоростью в шероховатых трубах при перемешивании (см.

гл. 4).

Закон сопротивления в однофазном потоке в соответствии с об­ общенным уравнением (3-134) можно выразить как зависимость Ар от скорости потока следующим образом:

1)ламинарный режим (Ар : : w)\

2)турбулентный режим (Ар : : wl>75) ;

3) инерционный режим (режим развитой турбулентности,

Ар :: w2) .

Согласно рис. 3-31, между ламинарным и турбулентным режи­ мами может существовать переходный режим движения потока.

ТЕЧЕНИЕ НЕНЬЮТОНОВСКИХ ЖИДКОСТЕЙ

Ньютоновские жидкости характеризуются линейной зависи­ мостью между касательным напряжением т и градиентом скорости

dw

dw/dn ^согласно закону Ньютона, т = ± ц - ^ - | при постоянных

температуре и давлении в условиях независимости их характери­ стик от времени. К ньютоновским жидкостям относятся все газы и чистые жидкости (а также их смеси) с низкой молекулярной массой. Для неньютоновских жидкостей эта зависимость нелиней­ на. По характеру зависимости

 

 

(3-146)

неньютоновские

жидкости разделяют на несколько классов.

К п е р в о м у к л а с с у

относятся вязкие (или стационарные)

неньютоновские

жидкости,

для которых функция (3-146) не зави­

сит от времени (рис. 3-33). По виду кривых течения в этом классе различают следующие группы.

Бингамовские пластичные жидкости (кривая 2 на рис. 3-33),

которые начинают течь только после приложения напряжения то (то — начальное напряжение сдвига или предел текучести).

При напряжениях, меньших то, бингамовские жидкости ведут себя как твердые вещества. Очевидно началу течения соответствует разрушение их прочной структуры. При напряжениях, больших то, поведение бингамовских жидкостей не отличается от поведения ньютоновских жидкостей, т. е. зависимость т от dw/dn для них также линейна. Структура бингамовских жидкостей при обратном уменьшении напряжения (т < т0) восстанавливается. Примером

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

Уравнение кривой течения в этом случае имеет вид

,

dw

(3-147)

т = т0 +

Ц п л -^ -

(при т > т0)

Здесь |іпл — постоянная (аналогичная вязкости обычной жид­ кости), называемая пластической вязкостью.

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

Сравним уравнение (3-146) с уравне­ нием, описывающим вязкое течение по аналогии с ньютоновской жидкостью:

^ — И-каж

dw

(3-148)

dn

 

Тогда

Ц каж -------

f (dw/dn)

(3-149)

 

 

 

 

dw/dn

 

 

 

 

 

 

 

 

 

Таким образом, если соединить ка­

 

 

 

 

кую-либо точку на кривой течения с на­

 

 

 

 

чалом координат, то тангенс угла накло­

 

 

 

 

на проведенной прямой (пунктир на рис.

 

 

 

 

3-33) будет представлять собой кажу­

 

' dw

 

 

щуюся вязкость.

 

 

V = f

для

неньюто-

 

 

dn

Для бингамовских жидкостей кажу­

 

 

 

новских жидкостей:

щаяся вязкость

 

 

( — ньютоновская

жидкость:

 

dw

 

2—бингамовская

жидкость;

 

 

3 — псевдопластичная

жидкость;

_ ^пл ЧгГ + т°

(3-150)

4дилатантная

 

жидкость;

Цкаж —

dw/dn

5—кажущаяся вязкость.

 

 

 

 

 

Из уравнения (3-150) следует, что кажущаяся вязкость умень­ шается с ростом градиента скорости сдвига dw/dn.

Псевдопластичные жидкости (кривая 3 на рис. 3-33) соста­ вляют наиболее представительную группу в классе вязких ненью­ тоновских жидкостей. К ним относятся растворы полимеров, цел­ люлозы и суспензии с асимметричной структурой частиц, вытяну­ тые молекулы которых перепутываются и при малых напряжениях тормозят течение.

Для этих жидкостей касательное напряжение, по Оствальду, определяется уравнением, которое 'часто называют степенным за­ коном реологии:

где k a m — постоянные.

Коэффициент k зависит от консистенции жидкости и увеличи­ вается с ростом вязкости. В табл. 3-1 приведены значения k и m

89

 

 

 

 

Т а б л и ц а

3- 1

К о э ф ф и ц и е н т ы

к и т у р а в н е н и я ( 3 - 1 5 1 )

д л я

р а з л и ч н ы х

 

с у с п е н з и й

п р и к о м н а т н о й

т е м п е р а т у р е

п о д а н н ы м

 

 

 

 

М е т ц н е р а

[ 2 7 ]

 

Суспензия

 

Концентрация,

k ,

т,

 

 

масс.

%

(кгс-сш/ м2) безразмерный

Глина в воде ...............................................

 

23,3

 

0,566

0,229

Карбоксиметилцеллюлоза в воде . . .

0,67

 

0,0308

0,716

 

 

 

1,5

 

0,318

0,554

Известь в

воде

 

3,0

 

0,945

0,566

 

33

 

0,732

0,171

Бумажная

масса в воде

.......................

4

 

2,04

0 5 /5

Цемент в воде ...........................................

 

54,3

 

0,256

0,153

для ряда жидкостей. Показатель степени т < 1 для псевдопластичных жидкостей (г. е. угол наклона кривой течения меньше 45°), однако часто приближается к этой величине (когда свойства жидкости приближаются к свойствам ньютоновских жидкостей). Кажущаяся вязкость псевдопластичных жидкостей уменьшается с увеличением градиента скорости dw/dn:

(при т < 1)

Кривая течения постепенно переходит в прямую с предельным зна­ чением ц/і = оо при бесконечно большом градиенте dw/dn.

В логарифмических координатах зависимость т = f (dw/dn) для псевдопластических жидкостей является линейной.

Дилатантные жидкости (кривая 4 на рис. 3-33) содержат жид­ кую фазу в таком количестве, чтобы она могла заполнить пустоты между частицами твердой фазы, находясь в состоянии покоя (или при очень медленном течении), и при этом обнаруживают свойства, близкие к ньютоновским жидкостям. При увеличении градиента скорости частицы твердой фазы начинают быстрее перемещаться относительно друг друга и объем суспензии начинает увеличи­ ваться. При этом жидкости уже недостаточно для заполнения уве­ личившейся порозности и кажущаяся вязкость увеличивается. При­ мером таких жидкостей являются суспензии крахмала, силиката калия, различные клеи (с большим отношением Т : Ж).

Касательное напряжение также определяется по формуле (3-151), но при значениях m > 1. Для дилатантных жидкостей ка­ жущаяся вязкость увеличивается с возрастанием градиента dw/dn, и для ее расчета может быть использовано уравнение (3-152) при значениях показателя степени пг > 1.

Ко в т о р о м у к л а с с у относятся иеньютоновские жидкости, характеристики которых зависят от времени. Для этих жидкостей кажущаяся вязкость определяется не только градиентом скорости сдвига dw/dn, но и его продолжительностью, что усложняет анализ

90

явления, так как для определения вязкости надо знать историю обработки жидкости.

Обычно для проведения классификации в жидкости проводят увеличение касательных напряжений, а затем их постепенное уменьшение. Если обе кривые течения совпадут, то можно утвер­ ждать, что время не влияет на вязкость испытываемой системы. Если кажущаяся вязкость изменяется со временем, то на графике x — f(dw/dn) будет получена петля, т. е. две отдельные кривые

(рис. 3-34).

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

ческой и тиксотропной жидкостей:

а —реопектическая жидкость, б —тиксотропная жидкость,' штриховые линии — кажущаяся вязкость.

Кажущаяся вязкость реопектических жидкостей увеличивается со временем (см. рис. 3-34, а).

Некоторые исследователи [27] предлагают использовать пло­ щадь петли на графике х = f(dw/dn) в качестве показателя реопектичности. К реопектическим жидкостям относятся суспензии бентонитовых глин и некоторые коллоидные растворы. После стоя­ ния реопектические жидкости приходят в первоначальное со­ стояние.

Тиксотропные жидкости отличаются тем, что в противополож­ ность реопектическим жидкостям их кажущаяся вязкость умень­ шается со временем (рис. 3-34, б).

Примером тиксотропных жидкостей являются многие краси­ тели. При продолжительном воздействии касательных напряжений структура этих жидкостей разрушается и увеличивается текучесть. Однако после снятия напряжений и стояния первоначальная струк­ тура постепенно восстанавливается и вязкость возрастает. Это свойство у некоторых красителей особенно заметно при окрашива­ нии ими вертикальных поверхностей — тиксотропность задержи­ вает стекание краски, так как ее вязкость со временем увеличи­ вается.

91

К т р е т ь е м у к л а с с у относятся вязкоупругие, или максвел­ ловские, жидкости. Их кажущаяся вязкость уменьшается под воз­ действием напряжений, однако после снятия последних вязкоупру­ гие жидкости частично восстанавливают свою форму. Примером вязкоупругих жидкостей являются некоторые смолы и пасты.

Для неньютоновских жидкостей характерно одновременное на­ личие свойств, типичных для твердого тела и для жидкости.

Предложено большое число эмпирических уравнений и «моде­ лей» для описания неньютоновских гидравлических систем. Устано­ вившееся состояние зависимости между т и dw/dn при постоянных температуре и давлении удовлетворительно описывают пять мо­ делей:

1)Бингама — уравнение (3-147);

2)Оствальда — уравнение (3-151);

3)Эйринга —■уравнение, полученное на основе кинетической теории жидкостей:

т = А sinh 1

dm \

I

(3-153)

du )

В .

 

 

где А и В — константы. Уравнение расхода по Эйрингу [11, 29]:

2nRB R2 cosh е —

(е sinh 8 — cosh е + 1)

(3-154)

Профиль скоростей:

 

 

да = BRe cosh е — cosh т

(3-155)

где

 

 

е _

*£*

(3-156)

21А

.4) Эллиса — обобщенное уравнение для скорости сдвига [28]

dw

— т + Ст5 (3-157) dn Но

на основе которого выведены уравнения расхода

т, _

nR3

R A p

,

nR3C ( R A p \ s

(3-158)

сек

4цо '

21

+

s +

3 \ 21

)

 

и профиль скоростей

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(3-159)

где С и s — константы.

 

 

 

 

(3-157)

превращается в урав­

В предельных случаях уравнение

нение (3-151) или уравнение Ньютона.

обобщенную зависи­

5) Райнера — Филиппова — представляет

мость в пределах изменения вязкости от 0 до оо [29].

 

Течение неньютоновских жидкостей является предметом изуче­ ния науки о деформациях и течении — реологии *. В последние

* От греческого rheos — поток.

92

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

Наиболее подробно исследованы реологические модели для вяз­ ких (стационарных неныотоновских) жидкостей, характеристики которых не зависят от времени [28].

Рассмотрим движение этих жидкостей при ламинарном и тур­ булентном режимах.

Ламинарный режим. Изотермический установившийся поток не­ ньютоновской жидкости в цилиндрической трубе находится под действием тех же сил, что и рассмотренный ранее (стр. 58) поток ньютоновской жидкости.

В общем виде можно записать, что градиент скорости сдвига пропорционален касательному напряжению:

4 ^ = / ( Д

(3 -1 6 0 )

Для потока, движущегося в цилиндрической трубе

dw

Hr)

(3 -1 6 1 )

HF

 

 

Если выделить в движущемся потоке жидкости цилиндр длиной I и радиусом г, то этот цилиндр будет находиться в равновесии под действием силы сопротивления, определяемой разностью давлений на обоих концах выделенного элементарного объема, и касатель­ ного напряжения, приложенного к его поверхности:

 

Аряг2 — 2nrh = 0

 

Откуда

Максимальное напряжение будет у стенки трубы:

Тст==~ Іг "’ где

^ — радиус трубы. Таким

образом:

 

т = тст-^-

(3-162)

Подставим зависимость (3-162) в уравнение (3-161)

dw

(3 -1 6 3 )

dr

Проинтегрировав уравнение (3-163), получим распределение ско рости в трубе:

оR

W

Г

или

R

(3 - 1 6 4 )

г

ѳз

Объемный расход жидкости, проходящей через элементарный объем:

dVсек — w dr

(3-165)

Интегрируем уравнение (3-165) по частям:

! dVQgK1”

оIw • 2r dr я

R

 

r2 dw о

(3-166)

Так как при г = R w — 0, то уравнение (3-166) упрощается. Под­ ставив w из уравнения (3-164), получим:

R

 

 

Тсек “ я j" r2f

j dr

 

О

 

 

Так как г — R —— (из уравнения

3-162), то уравнение

расхода

Тех

 

 

в общем виде:

 

 

я R3

х2f (т) dx

(3-167)

 

Подставляя в уравнение (3-167) ской моделью различных групп

 

 

du)

й

т

77F

f(x) в соответствии с реологиче­ неньютоновских жидкостей, полу­ чим расчетные уравнения для ла­ минарного движения в цилиндри­ ческой трубе.

1. Для бингамовских пластич­ ных жидкостей (см. уравнение

3-147)

_ dw

т - То — Цпл

Откуда

Рис.

3-35. Распределение скоростей,

!(х ) — —----

при

т0< т < т ст (3-168)

касательных напряжений и гради­

ента скорости в потоке бингамов­

 

Цпл

 

 

 

ской

жидкости при ламинарном

причем

f (т)

= 0

при 0 <

т < то-

 

режиме.

На

рис.

3-35

приведена

схема

 

 

движения бингамовской жидкости в трубе, причем следует отметить, что хотя напряжение т падает до 0 в центре трубы, однако существует зона (вблизи осевой ли­ нии), в которой отсутствует скорость сдвига. В этой зоне жидкость движется как твердый поршень.

При подстановке

функции (3-168) в уравнение (3-167) полу­

чим:

 

Vсек

(3-169)

94

После интегрирования:

 

 

я R3

X

------Тп

 

 

(3-170)

 

 

 

4

**'£ф

 

 

ИЛИ

 

[-f

 

 

 

 

 

 

 

 

 

 

 

4 Ц п л

L

 

То3 + -т с х

 

3

(3-171)

Подставив тст —

, получим

уравнение,

которое

носит назва­

ние уравнения Бэкингема, выражающее связь Ѵсек бингамовских жидкостей (при

их ламинарном

течении по цилинд­

ДО ‘V

рической трубе)

и перепадом давле­

 

ния Др:

 

о,6

между расходом

1

'•4

 

я А р

 

2/Тп

 

 

 

*

 

 

>

 

 

 

 

 

 

АpR

 

 

 

0,2

 

V

N

 

 

 

+

1

( 2/т° У]

 

(3-172)

о

 

4

з \

 

 

 

 

 

/ /

 

где Цпл —

3

I ДрР j

J

 

0,2

 

 

/

 

 

 

 

 

У

 

 

 

пластическая

вязкость;

 

 

 

А

'

 

то — предел текучести.

 

 

 

 

 

 

 

 

 

 

 

0,6

 

 

 

 

В тех случаях, когда область, за­

 

 

 

 

 

 

 

 

 

 

 

нимаемая ядром потока (поршнем),

 

 

 

 

 

 

мала,

уравнение

 

(3-172)

 

 

можно

 

1,0

10

 

20

30

упростить. Если rn/R = 0,5,

то мож­

 

 

 

 

 

 

 

 

ф/ uJcp

но отбросить последний член пра­

 

 

 

 

 

 

 

 

 

 

 

вой части уравнения и ошибка со­

Рис. 3-36.

Профиль скорости для

ставит 5,9%; для r„/R — 0,4 ошибка

жидкостей,

поведение

которых

составит только

1,8%. Приближен­

характеризуется

степенным рео­

логическим. законом при ламинар­

ное уравнение расхода будет иметь

 

 

ном режиме:

 

вид:

 

 

 

 

 

 

 

/ —псевдопластичная жидкость (т = ]/3);

..

 

яR3 ( hpR

4 То

 

 

2—ньютоновская

жидкость

(m=l);

 

 

(3-173)

 

3—дилатантная жидкость ( т = 3).

 

 

 

 

 

 

 

 

 

 

4Цпл \

2/

 

 

 

 

 

 

 

 

 

 

Профиль скоростей (рис. 3-36) можно получить при подста­ новке в уравнение (3-164) функции

/_дpR_

М-пл \ 2/ R То

Тогда

 

R

Д рг

 

 

 

ч

т )0d r

 

 

 

Мтіл 2

~1Г

 

 

 

После интегрирования получим

 

 

 

 

W ■

А Р (R2

П - f i R -

■г)

(3-174)

4 ІЦп

 

 

 

 

при гя < г < R. При граничных условиях 0 <

г -< гя

 

 

Ар

гя)2

 

(3-176)

 

2Lun (R -

 

 

 

 

 

1

 

 

95

Скорость сдвига у стенки трубы:

dw

ДpR

To

(3-176)

dr

2/Цпл

И

 

2. Для псевдопластичных и дилатантных жидкостей в соответ­ ствии со степенным законом

%= k dwdu

получим

Нт) =

Подставив функцию (3-170) в уравнение (3-167), запишем:

 

 

тст

2 Ч----

 

 

 

nR3

Г

т т

dx

 

 

J

k_і_

 

 

 

СТ

о

m

 

 

После интегрирования

 

 

 

 

З т + 1

 

 

 

 

 

 

 

пR3

 

 

т

/ Тст\ г

У сек —

3т +

 

Зот +

1 \ k j

 

 

kк mxzТст

(3-177)

(3-178)

(3-179)

Подставив значение тст, после перегруппировки получим расчетное уравнение:

г -

т

лп3( R äp\ k

(3-180)

сек

Зт + 1

R \ 2Ik )

 

Профиль скоростей будет характеризоваться зависимостью:

R 1 R I

(3-181)

После интегрирования получим:

 

W ■

( т+І

m+i \

т

(3-182)

+

 

или

т+ 1

 

 

 

W =

wv

(3-183)

Уравнение (3-182) можно преобразовать с помощью средней ско­

рости:

Г т+1

“' = ^ p ( ^ Lf 1L) [ l - ( - ^ - ) m j

(3-184)

Скорость сдвига у стенки для псевдопластичных и дилатантных жидкостей рассчитывается по уравнению:

dw

 

Зт -f-

(3-185)

dr

nR3

т

 

96

Соседние файлы в папке книги из ГПНТБ