Курс механики сплошных сред
..pdfа сложную функцию 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. В тон |
|||||
ких |
же |
трубах |
медленные |
течения, |
как |
правило, ламинарные. |
||
Введение в течение подкрашенной |
струйки |
позволяет |
выявить |
фундаментальное различие между турбулентными и ламинарными режимами. В ламинарном режиме струя долгое время остается
прямолинейной и диффундирует очень слабо; этим подтверждается тот факт, что теоретическая схема, постулированная вначале, хорошо описывает траектории частиц. В тур булентном режиме диффузия проис ходит очень быстро, и жидкость вско ре окрашивается полностью (рис. 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, положительным или отрицательным, схемы тече ний и переходных режимов могут сильно различаться.
Из сказанного прежде всего вытекает, что решение рассматри ваемой задачи не единственное. Рассмотренное выше ламинарное течение, являющееся к тому же самым простым, существует всегда для каждой данной точки из дважды заштрихованной области, но наблюдаемое на практике течение оказывается не таким, так как
* Точный |
фиаический смысл числа Рейнольдса и понятие подобия течений |
будут даны в |
рамках механики жидкостей. |