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

Курс механики сплошных сред

..pdf
Скачиваний:
19
Добавлен:
15.11.2022
Размер:
20.58 Mб
Скачать

а сложную функцию N (Я (5)) —

 

N (S) = N (Я(5)).

(52)

Очевидно, Я (5)— нечетная функция S, N (S) —четная функция S. Сравнивая решение (51) с решением, найденным в VI.2.2, заме­

чаем, что сдвиговое напряжение

уже больше

не является

ли­

нейным относительно k и что в направлениях xt

и х2 к давлению

добавляются нормальные напряжения

N (k). В некоторой точке

М

для направления е2 оси х2 вектор напряжения будет:

 

Т(М, e2) = (— p + N(k))e2 + x (k) et.

 

 

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

IX.4.2. Течение между двумя параллельными плоскостями. Пусть

(как и в IX.2.2) даны две

плоскости П и П'

с уравнениями x2 = h

и х2= — h, а

исследуемое

течение,

удовлетворяющее условию

не­

сжимаемости,

имеет вид

 

 

 

 

 

 

 

U1 = u(x2)1 U2 = U3 = 0.

 

 

(53)

В любой точке все компоненты

D{j равны

нулю,

кроме

 

 

 

__ 1

dи

k_

 

 

(54)

 

Z?i2 — D21 2

dx2

2 •

 

 

Таким образом, k является

неизвестной

функцией

переменной

х2У

и, следовательно, все %{J—функции только одной переменной х2, полученные при замене в формулах (51) k выражением k(x2).

Так как ускорение отсутствует, а внешние массовые силы пред­

полагаем равными нулю,

то

уравнения

движения,

сводящиеся к

виду (г,/, / = 0, дают

 

 

 

 

 

др_ _dTi2

др

_ С1т22

др

(55)

дх±

dx2 9

дх2

dx9 *

дхя

 

Таким образом, давление р зависит только от xt и х2 [согласно (558)]; оно является аффинной функцией хх [согласно (55x)]t ибо коэффициент при xi не зависит от х2, как это следует из (552). Функция р имеет, таким образом, вид

р = Р х,+ у(х2),

где Р—перепад давлений вдоль xt. Но тогда из (55J следует, что

т12 =

Р (х2 о) ,

 

где а—некоторая постоянная. С учетом законов поведения

имеем

S = T(k) = — Р (х2 — а)

или &= — Х{Р(х2 —а)\.

(56)

Это решение в принципе позволяет определить и (х2) в соответ­ ствии с (54), а остающаяся неизвестной функция у(х2) находится

 

из уравнения

(552):

 

 

 

 

У(^2) ^ Т22 Н“ Ро=

N {Р (.X2

и)} -f- p0t

(57)

 

где р0 — константа.

 

решения

^

 

Дальнейшая

конкретизация

определяется

заданием внешних

воздейст-

*2 вий и краевых условий.

 

 

 

 

 

Течение Куэтта. Здесь плоскость ГГ

Рис. 12

фиксирована, а плоскость

П движется

по­

 

ступательно

со скоростью V параллельно

оси Ох1- Перепад давлений

отсутствует: Р = 0.

 

 

 

Следовательно, сдвиг k

 

постоянен во всей области; м —аффинная

функция ха, определяемая

условиями прилипания1

 

 

 

 

 

u = y h t r -

 

 

 

 

(58)

Таким образом, распределение скоростей не зависит от закона поведения.

Давление р равно произвольной постоянной р0, и можно написать,

что

 

^83 =

Ро* ^12 == ^ ( 2/1)» ^11 ^

^22=

"Ь^

 

 

остальные компоненты alf равны нулю.

проявляется,

в

частности,

Неньютоновый

характер

жидкости

в «эффекте нормального давления», который

приводит к тому, что

<т88 не равно уже ни сги, ни а22.

 

 

 

 

 

 

Обе плоскости

неподвижны и задан перепад давлений. Учитывая

то,

что X (S) — нечетная

возрастающая

функция S,

видим, что функ­

ция

/г(х2), определяемая

из

формулы (562), имеет график, показан­

ный на рис. 12. Так как функция и является

первообразной для /г,

то площадь заштрихованной

части должна (алгебраически) равняться

нулю. Отсюда

необходимо следует, что а = 0. Следовательно, распре­

деление скоростей

дается формулами:

 

 

 

 

 

 

“-$ > < « )« - Г_„я (П)«л=£ мPI)di_

f“ х (s)ds.

 

 

 

 

 

 

 

 

 

 

 

(59)

 

Ненулевые

компоненты

тензора

напряжений

имеют

значения:

 

 

 

^ 3 3 5=5 ^ * 1

Ро» ^ 12 = ^21

=

 

 

 

 

 

 

 

аи =

= PXi — {p0 + N (Рхъ)\.

 

 

(60)

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

Q = 2/ш =

и (*,) dx, = [x,u (х,)]1 Й—

хй dx*.

Так как и ( +h) = и ( —h) = 0, то из формул (562) имеем

Q= П

я ( ^ , ) X, d*2 = ±

Я, (S) 5 dS =

(S) S dS. (61)

При фиксированном Л величина Q для данной жидкости

зависит

только от Р и имеем

 

 

 

 

j p { P * Q ( P ) \

= 2Ph*X(Ph).

 

(62)

Итак, изменяя Р и измеряя

расход Q (Р), можно найти

по фор­

муле (62)

функцию %(S), обратную вискозиметрической

функции

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

IX.4.3. Течение в круглой трубе. Рассмотрим вновь течение

Пуазейля.

Используя цилиндрические

координаты

х, г, 6(х1 = х2,

x2 = rcos0,

x8 = rsin0), построим установившееся

течение в виде

 

Ut = u(r), £/2 = i/8 = 0,

0 < r < P .

(63)

Для удобства используем в каждой точке естественный ортонормированный репер 91, соответствующий цилиндрическим координа­ там*. Обозначим через Схх, Схр = Срх, CXQ= CQx. .. компоненты в 91 некоторого симметричного тензора второго ранга С. Очевидно, что выражения (63) определяют локально некоторый тензор скоростей

деформаций простого сдвига; более точно, в репере 91 все компо­ ненты 3) равны нулю, кроме одного:

п 1 du 1 и 'D*r ~~2 ~37~ ~2к'

Таким образом, k является функцией переменной г. Согласно соотношениям (51) ненулевые компоненты тензора вязких напряже­

ний даются в 91 формулами

*** “ Ь г N № . Тг* = ^хг — т (* )

(6 4 >

и зависят, таким образом, от одной переменной г.

Что же касается уравнений движения, то при отсутствии уско­ рений они запишутся в таком виде:

dr ,

1

др

I I LC 1

1

др_

1_др п

(65)

dг*

г т хг 1=1 дх*

dг *

г гг ^

дг *

г Ъ6 =

 

Рассуждения и выкладки будем проводить таким же образом,

как и выше. Сначала покажем, что

 

р = _ Р х + у(г),

(66)

где Р — постоянный перепад давления, необходимый для поддержа-

* Некоторые из приведенных здесь формул даются в ГНУ.

ния стационарного течения. Уравнение (65j) тогда дает

 

 

* ( £ ) .

Имеем далее

 

 

тя, = тг, =

Я ( - £ ) ;

(68)

производная у (г) определяется из

уравнения

(65„):

£-4 жМт)}-

<69>

Скорость находится интегрированием уравнения (67) с учетом условий прилипания к трубе, которую считаем неподвижной:

“ <r>“ I > ( T ) di - T C 4 ( s >ds-

(70)

Обратим особое внимание на выражение для расхода:

Q = 2K ^ u ( r ) r d r = 2 K { [ « 4 ] oW- ^ - ^ - J - d r | .

С учетом условия прилипания и (R) = 0 и формулы (67) имеем

Q = л §* Я (-£ ) Г* dr = -g. ^ й/2Я (S) S2 dS.

(71)

Если при неподвижной трубе изменять перепад давления Р, то после дифференцирования по R получаем уравнение, которое поз­ воляет определить вискозиметрическую функцию X(S) на основа­ нии измерения расхода:

<72>

Решение для ньютоновой жидкости отсюда вытекает как част­ ный случай; в самом деле, если в выражении (71) положить

T (£ ) = p ft, M S ) = ^ - ,

то вновь приходим к решению (26).

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

Эта формула обобщает полученный выше результат (31).

Полученные выше результаты легко распространить на более широкие по сравнению с рассмотренными в разделе IX .4 классы жидкостей, например на материально простые среды, обладающие непрерывной памятью. Здесь мы рассмотрим вязкопластические жидкости и, в частности, несжимаемую жестко-вязко-пластическую среду, называемую жидкостью Бингама. Известно, что закон пове­ дения такой жидкости записывается в таком виде (VIII. 116):

ои = — рб,] + T,/f т,7 = g

+ 2|xDu (D’n > 0),

(74)

где g — предел текучести на

 

Uu

 

 

сдвиг

(по напряжениям); р,—коэффи­

циент вязкости. Кроме того,

если

Dty= 0,

то величины

т17 не оп­

ределены, но они удовлетворяют условию

(VIII.120):

 

/ 1

\

1/2 .

 

 

[т тиЪ?)

IX.5.1. Течение между двумя параллельными неподвижными плоскостями. Не производя заново расчетов и используя предыду­ щие выкладки и результаты, полученные в IX.4.2, видим, что для течения (53) все вязкие напряжения т|7 равны нулю, за исключе­ нием компонента т12, и что

Xii = x(k) = tik + g(sigak),

ЫфО,

\

/75\

т12= т (0).

ft= o,

/

 

где signfe равен либо + 1 , либо — 1 в зависимости

от того, поло­

жительно или отрицательно k.

 

 

 

Таким образом, для изучения течения жидкости Бингама до­

статочно приравнять N (k) = 0 и взять в качестве %{k)

разрывную

функцию, определяемую соотношениями (75).

 

Заметим, что обратная

функция A,(S) непрерывна:

 

( р-1 (S — g)

при

(76)

MS) = |

0

при —

t p _1(S + g)

при S < —g.

 

Если, например, плоскости неподвижны и перепад давлений Р задан, то формула (59) позволяет сразу же получить распределение скоростей. Действительно при S > g первообразная функции К (S) равна (2р)-1 (S — g)2, так что мы имеем!

u = j ( h —xt) [ p (/Ц — — g ) , если

y < x i < h ,

Ph> g

 

 

(77)

и =-|j- (ft—-fr)\

если 0 < * a< | .

Расход дается формулой (61):

 

Ph<g,

0,

если

Q =

 

и

(78)

ж ( 1 _ и ) ’ ( 2 + А ) ’ если

> « -

Рис. 13. Вискозиметрическая функция т(&) и об­ ратная функция К (S):

--------------ньютонова жидкость ( / / ) ; ----------------жидкость Бингама ( / ) ; --------------собственно жидкость ( I I I )

На рис. 13,а приведен график функции т(А), а на рис. 13,6 — функции k(S) для трех жидкостей—ньютоновой, бингамовой и не­ ньютоновой, характеристики которых—промежуточные между харак­

теристиками первых двух. На рис.

14 даны графики

первообразных

функций K(S) для всех трех жидкостей. Для каждого фиксирован­

ного Р можно построить, с одной

стороны,

профиль

скоростей,

а с другой—получить расходы. Из

рисунка

сразу

же

видно, что

наибольший расход у ньютоновой жидкости.

 

IX .4.3 также

IX.5.2. Течение в круглой трубе. Результаты

распространяются и на случай жидкости Бингама. Приведем резуль­ тат, относящийся к расходу (71). Легко находим (задача 32):

 

 

если

PR < 2g,

f - l A +

T p f c ) ’

если

(79)

P R > 2 g .

На рис. 15 приведена зависимость

расхода

Q от перепада дав-

--------------ньютонова жидкость

(//);

_________

жидкость Бингама ( / ) ; --------------

собственно жид­

кость (111).

 

 

Профиль скоростей получается сечением поверх­

ности Л (S)

плоскостями

S=P/i. Скорости даются

векторами, проведенными

ог ординат /,

I I , 111 до

кривых I,

/ / , I I I . Значения дебита

численно

равны площадям, :<аштрихованным этими векторами

Рис. 15. Дебит круглой трубы заданного сечения в зависи­ мости от перепада давления:

- ---------ньютонова

жидкость;

--------------жидкость

Бингама;

------------- собственно

жидкость

•пений Р для трубы радиусом R и для сравнения даны расходы Для ньютоновой и неньютоновой жидкостей (вискозиметрическая функция т(£) последней приведена на рис. 13, а).

IX.6 . УСТОЙЧИВОСТЬ ТЕЧЕНИЙ.

СУЩЕСТВОВАНИЕ ТУРБУЛЕНТНЫХ РЕЖИМОВ (НЬЮТОНОВЫ ЖИДКОСТИ)

В этом разделе ставится цель —сравнить найденные результаты с данными эксперимента. Это позволит определить пределы приме­ нимости излагаемой в настоящем курсе теории жидкостей. Ограни­ чимся рассмотрением ньютоновой жидкости и начнем с рассмотрения течения в круглой трубе или вообще в трубопроводе, что позволит выявить тот факт, что кроме изученных течений (называемых ла­ минарными) существуют совершенно отличные от них течения, на­ зываемые турбулентными. Затем на примере течения внутри двух концентрических вращающихся цилиндров будут приведены некото­ рые замечания—весьма общего характера—относительно устойчи­ вости течения, рассмотренного в IX ,2.4. Изучение вопросов об устойчивости течений (которое здесь не проводится) позволяет уточ­ нить границы применимости найденных решений и выявить началь­ ные механизмы перехода от ламинарного режима к турбулентному.

IX.6.1 . Турбулентные режимы течений в трубах. Изученное в IX.2.3 течение Пуазейля дает возможность легко сопоставить тео­ ретические выводы и данные экспериментов. С этой целью целесо­ образно представить полученное решение в безразмерной форме.

Пусть

диаметр трубы d = 2R; pt —р, —разность

давлений между

двумя

отмеченными точками трубы,

отстоящими

друг от друга на

Расстоянии /; и —средняя

скорость

течения. Введем два безразмер­

ных числа X и Re, определяемых по

формулам!

 

 

_ 2Pd

_ 2(pi —Pi)d

n c _ pud

(80)

 

pu2

puH

И ’

 

 

где X—коэффициент падения нагрузки; Re—число Рейнольдса. Формула (26), позволяющая оценить среднюю скорость в зави­

симости от падения нагрузки, принимает универсальный вид

64

X (81)

Re •

На рис. 16 показана схема опыта, поставленного для изучения данной проблемы. Жидкость в резервуаре (А) поддерживается на постоянном уровне, с тем чтобы сохранить постоянство давлений Pi—Рг- В воронке (С) содержится подкрашенная жидкость, питаю­ щая вводимую в течение струйку, за которой ведется наблюдение. Расход определяется путем измерения объема накопленной в сосуде

(г) жидкости. Труба изготовлена таким образом, чтобы стенки были как можно б0лее гладкими *.

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

Рис. 16. Схема эксперимента по изу­ чению течения в трубе:

Л-питающий резервуар. Заслонка R t регу­

лируется

таким

образом,

чтобы

уровень

жидкости оставался постоянным;

М —мано­

метр для

измерения

перепада

давления

рх- р г\ i -

сосуд

для

сбора

отработанной

жидкости (для измерения дебита);

С -сосуд

для введения струйки окрашенной жидкости

Рис. 17.

Изменение

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

потери

нагрузки

в

зависимости от

числа Рейнольдса

для разных

режи­

 

мов течения:

 

/ —Я=64 Не_ 1 -ламинарный режим;

2-Я =

я0,316 R e" 1/ 4 -турбулентный режим (аб­ солютно гладкие стенки); 3 -переходный режим; 4 -искаженне решения при боль­

ших числах Рейнольдса из-за шероховатости стенок

На рис. 17 приведены полученные результаты. На осях абсцисс и ординат значения Re и А, даны в логарифмических координатах. Очевидно, существуют два основных режима.

Первый режим, называемый ламинарным, отвечает малым чис­ лам Рейнольдса, для него формула (80) находится в замечательном соответствии с опытом.

Второй режим, называемый турбулентным, соответствует боль­ шим числам Рейнольдса. Приближенные значения % как функции Re даются для довольно широкого диапазона значений Re следую­ щей эмпирической формулой:

X = 0 , 3 1 6 R e - 1/4.

(82)

На основании многочисленных опытных данных

можно сделать

следующее заключение: течение всегда ламинарное

при Re < 2000

и, как правило, турбулентное при Re > 4000. Эти

цифры имеют

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

Течения в промышленных трубопроводах в основном турбулент­

ные.

При

обычных температуре и давлении,

например,

течение

воды

в трубе диаметром 10 см со скоростью

1

м/с и движение воз­

духа

в воздухопроводах диаметром 30 см при средней скорости

5 м/с

происходит

при числах

Рейнольдса порядка 100 000. В тон­

ких

же

трубах

медленные

течения,

как

правило, ламинарные.

Введение в течение подкрашенной

струйки

позволяет

выявить

фундаментальное различие между турбулентными и ламинарными режимами. В ламинарном режиме струя долгое время остается

/-ламинарный режим; 2 -переходный режим; 3 -турбулентный режим
Рис. 18. Диффузия струйки:
окрашенной

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

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

В переходном режиме при числах Рейнольдса порядка 3000 тра­ ектория окрашенной струи начинает размываться (диффузия прояв­

ляется

лишь на некотором расстоянии), ее можно еще

различить,

но она

имеет вид как

бы синусоиды, обвитой

вокруг

прямой,

по

которой

следовала бы

частица в ламинарном

режиме.

Отсюда

вы­

текает, что при возрастании числа Рейнольдса установившийся ре­ жим, рассмотренный в IX .2.3, перестает быть таковым. Несмотря на то что течение может быть исследовано в рамках механики сплошных сред, оно выходит из класса течений, рассматриваемых в теории Пуазейля. Можно предполагать, что с увеличением числа Рейнольдса количество механизмов, вызывающих незатухающие воз­ мущения, возрастает настолько, что, начиная с некоторого порого­ вого значения, скорость как функция времени в каждой точке уже не может описываться непрерывно дифференцируемой функцией, которую можно сравнить с экспериментом.

Эти замечания, очевидно, весьма схематичны. Из них, однако, видно, что если найдено какое-либо решение уравнений Навье — Стокса, то следует сразу же исследовать его устойчивость. Гидро­ динамическая устойчивость является в наши дни очень важным раз­ делом механики, ее важность для теории и практики трудно пере­

оценить.

 

 

 

между двумя

IX.6 .2. Гидродинамическая устойчивость течения

концентрическими

цилиндрами. Течение, изученное в

IX .2.4, было

предметом весьма

детальных

исследований, занявших

не один де­

сяток лет. Здесь

приведены наиболее примечательные результаты —

в виде функций

безразмерных

параметров:

 

R = , - ^ . R e . - * S

Рис.

19. Устойчивость течения

Ь

между двумя цилиндрами.

Рис. 20.

Течение

Косой

штриховкой

показана

об­

ласть

ламинарного

течения;

квад­

в виде

тороидаль­

ратной

штриховкой-течение в виде

ных

вихрей

 

тороидальных вихрей

 

— чисел Рейнольдса*, включающих заданные величины задачи. За­ метим, например, что формула (31) может быть записана в следую­ щем безразмерном виде:

Не уменьшая общности, можем считать

 

й0 и, следовательно,

Re, положительными.

 

найденное

в

IX .2.4,

Можно убедиться в том, что решение,

хорошо согласуется с экспериментом, если

значения

Re,

и Re,

соответствуют одной

из точек, находящихся

внутри заштрихованной

области на рис. 19.

Так обстоит дело на

практике, в

частности,

когда внутренний цилиндр неподвижен.

Что же касается точек из области, имеющей двойную штриховку, то они соответствуют наблюдаемому экспериментально течению с тороидальными вихрями. Траектории в этом случае уже не являются окружностями с центрами на оси Охл, они не выходят за пределы своих ячеек; их след в меридиональной плоскости, вращающейся вместе с жидкостью, образует кривые, изображенные на рис. 20.

При значениях Re, и Re,, соответствующих точке, находящейся несколько выше границы заштрихованных зон, ячейки несколько «гофрированы» и вращаются с угловой скоростью, равной средней угловой скорости обоих цилиндров. Если при неизменном Re, уве­

личивать

Re,, течение может еще

более

усложниться и перейти,

в конце

концов, в турбулентный

режим.

В зависимости от того,

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

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

* Точный

фиаический смысл числа Рейнольдса и понятие подобия течений

будут даны в

рамках механики жидкостей.