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

книги из ГПНТБ / Бекнев В.С. Газовая динамика газотурбинных и комбинированных установок учеб. пособие

.pdf
Скачиваний:
48
Добавлен:
25.10.2023
Размер:
14.39 Mб
Скачать

В результате первое уравнение примет вид

1 dp _ d-vx

цdx dz-

Проннтегрнруем полученное уравнение 2 раза по z, тогда

г2 dp

,

г

г

х ~~ 2ц dx

'

l Z i

-

Определим произвольные постоянные интегрирования С\ и С 2 .

Из граничных условий скорость ѵх=

0 на стенке

при z =

±/г,

значит

= 0 и С,

= • — С л е д о в а т е л ь н о ,

скорость

 

 

 

 

_

(г» - /.»)

dp

 

 

 

 

 

 

 

А

'

d.v "

 

 

 

 

Эпюра скоростей в сечении канала представляет собой пара­

болу

с вершиной при

z =

0.

Поскольку скорость ѵх^> 0,

a z2

— / г

< 0 ,

получаем,

что

градиент давления

вдоль

оси

х

будет

-£-«>•

 

 

 

 

 

 

 

 

 

Это уменьшение давления

в направлении

движения жидкости

обусловлено действием сил вязкости, единственными рассматри­ ваемыми нами силами.

Уменьшение давления Ар на длине Ах возможно найти, интег­ рируя полученную формулу для скорости ѵх. Однако пользо­ ваться полученным результатом затруднительно, так как не­

известна истинная скорость ѵх. Поэтому

введем

в рассмотрение

среднюю скорость

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

j рох

dz

 

 

 

 

 

 

У с р = - Ѵ -

 

 

( 1 8 3 )

С учетом формулы для ѵх

получим,

что

 

 

1_

f

z*-—h* h-

(dpdp

\ w ,.

 

(2hf_dp

 

 

+ h

 

 

 

 

 

 

 

СР

2/1

.J

V dx

)J

 

~

12ц

dx

 

 

—Ii

 

 

 

 

 

 

 

Отсюда интегрированием

по

х

в

предположении, что при

X = 0 давление р

=

р 0 , получаем

искомую потерю давления, т. е.

 

Ро — р = Ар =

- ^ -

и с р

A.v.

 

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

180

Подобное течение в случае движения электропроводящей жидкости в поперечном магнитном поле называется течением Гартмана. Рассмотрим плоское стационарное ламинарное течение вязкой несжимаемой электропроводящей жидкости между двумя изолированными пластинами в поперечном внешнем магнитном поле. Будем полагать (рис. 87), что скорость потока ѵ направ­ лена по оси X, приложенное постоянное и однородное магнитное поле В — по оси z, а электропроводность а постоянна. Начало' координат находится посредине между пластинами, которые будем считать изоляторами. Примем, что 2/і много меньше и /, так

Рис. 87. Направление индуцированного магнит­ ного поля от токов, теку­ щих в плоскости z = О вдоль оси у:

J

— изолятор; 2

э л е к т р о д ;

3

— с о п р о т и в л е н и е н а г р у з к и

 

или источник

энергии

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

Пусть боковые пластины являются электродами. Если они соединены с внешней нагрузкой, то жидкость движется под дей­ ствием приложенного перепада давления (преодолевая пондеромоторную силу), что соответствует каналу МГД-гнератора. Если электроды соединены с источником электрического тока, то жид­ кость движется под действием развивающейся в ней пондеромотор­ ной силы, преодолевая внешний перепад давления, что соответ­ ствует каналу МГД-насоса. Если боковые пластины выполнять изоляторами, то ток в жидкости, движущейся под действием

перепада

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

В силу

сделанных предположений подобно течению Пуазейля

скорость v имеет ^-компоненту, которая зависит только от ко­ ординаты z:

• Ъ[ѵх{г), 0, 0].

Магнитное поле внутри канала будет слагаться из наложен­

ного поля В — const и из

индуцированного поля Б и н д

от токов,

текущих по жидкости. На

рис. 88 схематично показаны

индуци-

181

рованные магнитные силовые линии для случая, когда электри­ ческий ток в плоскости z = 0 течет перпендикулярно плоскости zx. В рассматриваемом случае это будет именно так, ибо при сделан­ ных предположениях нет причин для отклонения тока от указан­ ного направления: эффект Холла считаем отсутствующим, сплош­ ные электроды выравнивают электрический потенциал вдоль оси X. Следовательно, можно считать магнитное поле внутри канала состоящим из двух компонент: постоянного наложенного

поля Вг = const и индуцированного

поля Вх, зависящего

только

от координаты z, т. е. В [Вх (z), 0,

Вг].

 

Запишем основные уравнения с

учетом того, что в

нашем

случае уравнение (26) неразрывности имеет вид div ѵ = 0, а массо-

 

 

 

 

 

вая сила Jm =

±-Qx~В).

 

 

 

 

Электрические токи

 

 

Р

 

 

 

 

 

 

 

 

 

 

Тогда уравнение движе­

 

 

 

 

 

 

 

 

 

ния примет вид

 

 

 

 

 

 

 

 

 

 

(о ѵ ) о =

1•gradp

+

 

 

 

 

 

 

 

 

 

+

| ( / х 5 ) +

ѵЛ5.(184)

 

 

 

 

 

 

 

 

 

 

Уравнения Максвел­

 

 

 

 

 

 

 

 

 

ла

используем в виде

Рис

Распределение

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

осевых

уравнений

(1)—(4),

а

обобщенный

закон Ома

 

скоростей по сечению канала

 

 

 

 

 

 

 

 

 

 

 

 

в виде выражения (8).

 

Продолжим рассмотрение возникающих упрощений. Из урав­

нения (4)

 

 

дВх

 

дВг

 

 

 

 

 

 

 

 

 

 

 

=

0,

 

 

 

 

 

 

 

 

дх

 

dz

 

 

 

 

 

 

следует подтверждение, что Вг

постоянно по координате z, так

как

поле Вх

от координаты х не зависит.

 

 

 

 

Из уравнения (1) следует, что

электрический ток имеет со­

ставляющую

только

по оси у:

 

 

 

 

 

 

 

 

 

 

.

_

1

дВх (г)

 

 

 

 

 

 

 

l

y ~

fima

dz

'

 

 

 

 

т, е. является только функцией z, а именно / [0, \у (г), 0] .

Перейдем к электрическому полю и найдем вектор Е из урав­ нения (8):

 

 

Ё =

±-]-(дхВ).

 

Векторное произведение (ѵхВ)

имеет только

^-составляющую,

равную —vz (z) Вг.

Следовательно, и поле

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

//-составляющую:

_

1

дВЛг)

 

 

р

vx(z)Bz.

 

и

~

ѵт

dz

 

 

182

От координат х и у входящие в это выражение величины не зависят, поэтому и производная [дЕуІдх = 0. Однако и от коорди­ наты z поле Еу тоже не зависит. Это видно из уравнения (2), имеющего в данном случае вид

dz

' дх

Поскольку дЕу/дх = 0, то Ец от координаты z не зависит. Значит поле Еу const, что также понятно из физических усло­ вий задачи: нет токов Холла и имеются сплошные электроды. Итак, электрическое поле имеет следующую структуру:

 

 

Ё

[0, Еу,

0] .

 

 

Из

уравнения (3), в

котором

для

рассматриваемого

случая,

div Е =

0, следует, что

в

жидкости

нет избыточных

зарядов.-

Таким образом, для практического использования остаются уравнения (1), (8), (184), которые проектируются на координат­

ные оси следующим

образом:

 

 

 

 

^ -

Т ^

+

І

^

+ ѵ ^ ;

(185)

 

° = — Н ? - :

 

( 1 8 6 )

 

0 =

- - |

-

- / А

;

(187>

Mm a

u z

іу = аЕу~аахВг.

(189)

Из уравнения (186) следует, что как и для течения Пуазейля, давление не зависит от координаты у, однако по высоте канала вдоль оси z теперь существует градиент давления

др

 

 

dz

~

 

—ІуВх,

 

 

 

 

который

является

искомой

величиной

наряду

со скоростью

ѵхг

магнитным полем Вх и электрическим током \у . Магнитное поле

Вг

и характеристики

среды р, er, та,

ѵ следует

задать.

Градиент

давления дрідх,

согласно уравнению (185), постоянен в силу

наших

допущений.

 

 

 

 

 

 

 

Для

решения

дифференциальных уравнений (185), (187)—

(189) в

качестве

граничных условий

следует

принять: ѵх = 0

на стенках при z = ±h; Вх

=

0

на оси х при

z = 0.

 

Последнее условие очевидно

из

рассмотрения рис.

87.

 

183-

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

'Ох

где ѵср

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

по уравнению (183);

 

 

z =

 

 

рѵ ср

 

Вх

5,Re„

 

 

 

 

 

 

 

 

 

 

 

 

 

Re,„ =

iimaavcph;

RE

- E:i

 

 

 

 

 

 

 

 

 

 

vcpBz

 

 

 

 

 

 

 

 

rt 2

a

D 2 , 2

 

lu

 

 

 

 

 

На

= —

BJi

 

oBzvcp

 

 

 

 

 

 

Решение

уравнения

(185)

 

 

 

 

будет

иметь

вид

 

 

 

 

 

 

 

 

На fch На —ch

(Наг)]

(190)

 

 

 

 

 

 

 

На ch На — sh На

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Из

уравнения

(189)

имеем

 

 

 

 

решение для

плотности

тока

 

 

 

 

 

 

 

] Y

= RE~vx.

(191)

 

 

 

 

 

 

Наведенное магнитное поле

 

 

 

 

Вх

из

уравнения

(188)

 

~1,0

 

 

 

 

 

 

sh (На г) — г sh На

(192)

 

 

 

 

 

 

 

Рис. 89.

Распределение

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

 

 

На ch На — sh На

 

 

 

 

 

 

 

 

 

плотностей тока по сечению

канала

 

Градиент

давления поперек

 

 

 

 

 

 

канала

переменен

и

из

уравнения

(187)

 

 

 

 

 

• _ â £ . =

_ J ^ L

H a

2

в

 

 

(193)

 

 

dz

 

Re

 

1 , J

х

 

 

 

Графически полученные результаты представлены на рис. 89— 91. Графики иллюстрируют влияние числа Гартмана на поведение поля осевых скоростей, поля токов и наведенного магнитного поля.

На графиках распределение параметров показано для трех

характерных случаев течения

жидкости:

1) критерий электрического

поля ^ £ = 1, что означает разрыв

внешней электрической цепи (режим холостого хода). Весь гене­

рируемый электрический

ток замыкается внутри канала;

2) R E — 0. внешняя

электрическая цепь замкнута накоротко

и генерируемый в канале электрический ток беспрепятственно течет по ней в режиме короткого замыкания;

184

3) RE

=

2, течение тока по каналу обеспечивается

питанием

от внешнего источника (режимы МГД-насоса).

 

Масштаб шкалы токов на рис. 90 соответствует течению с пара­

метром

RE

= 1, для течений с параметрами RE = 0

и RE = 2

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

Из рис. 88 видно, что с ростом числа Гартмаиа профиль ско­ ростей в средней своей части становится более плоским, а распре­ деление скоростей у стенок — круче. Объяснение этому факту будет дано ниже на примере

Рис.

90.

Распределение безразмерной

Рис. 91. Распределение осевых скоро-

напряжеиности

индуцированного

маг-

стей при

течении Пуазейля

нитного поля

по сечению канала

 

 

 

Разберем действие пондеромоторной силы/Ѵ на элемент жидко­

сти

в

канале. Учитывая,

что

электрический

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

«/-составляющую, а магнитное поле х- и 2-составляющие, получим,

что сила /у = / X В = jyBzi

jyBxk, так

что

f\'x — lißz, fvy — 0.

 

fVz

= -juBx.

(194)

Рассмотрим, каковы знаки составляющих магнитного поля

иэлектрического тока.

Вэтом режиме течения эпюра распределения индуцированного магнитного поля Вх (рис. 91) такова, что в верхней части канала

при z > 0 поле Вх < 0 (в нижней его части при z <

0 поле Вх

>

0,

а для двух других режимов течения наведенное магнитное поле

Bx

существенно больше). Составляющая Вг всегда

больше

нуля.

Распределение электрического тока в рассматриваемом режиме холостого хода (RE = 1) таково, что (см. рис. 90) у стенок канала электрический ток течет вдоль оси у, т. е. \ у >> 0, а внутри ка-

185-

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

Действительно, из закона Ома в нашем случае электрический

ток

 

 

jy = оЕу оѵхВг.

Проинтегрируем

это

выражение по высоте канала, учтем

что в направлении

оси

у полный

ток

 

 

 

 

 

{ jy dy

= О,

 

 

—h

 

и введем в рассмотрение среднюю скорость из уравнения (183). Тогда

 

 

 

Еу -

Вгѵср

 

 

 

:и выражение

для

электрического

тока примет

вид

 

 

 

 

 

Іу = оВгср

— ѵх).

 

 

 

Итак, при

ѵср

>

ѵх (у стенок)

электрический

ток

} у >• О,

т. е. ток течет вдоль

оси у; при ѵср

•< ѵх (в ядре потока)

\ у < 0 ,

т. е. ток течет вдоль оси — у.

 

 

 

 

 

Возвращаясь

к составляющим

 

выражения

(194)

пондеромо-

торной силы /V, с учетом полученных результатов можем сделать

•следующее заключение:

у

стенок,

где

ѵср

>

ѵх,

составляющая

пондеромоторной силы

fVx

> 0 и

ускоряет

движение жидкости,

в ядре потока — наоборот сила fVx

< 0

и тормозит это движение.

Действие составляющей

fVx

показано

на рис.

89.

 

Рассмотренное изменение эпюры

скорости

движущейся элек­

тропроводящей жидкости под действием приложенного магнит­ ного поля и составляет главное содержание эффекта Гартмана. С ростом градиента скорости у стенки увеличивается сила трения: приближенно можно считать, что отношение коэффициентов тре­

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

равно ѴзНа.

§ 24. ЭЛЕМЕНТЫ

ГИДРОДИНАМИЧЕСКОЙ

 

ТЕОРИИ СМАЗКИ

 

 

Уравнения движения вязкой жидкости позволяют рассчитать

подшипник скольжения.

 

Ограничимся

упрощенной моделью явления

применительно

к течению масла в щели между неподвижной колодкой и враща­ ющейся опорой упорного подшипника скольжения. Будем счи­ тать, что движение несжимаемой вязкой жидкости между пло­ скостью, двигающейся со скоростью U0, и неподвижной криво-

186

Рис. 92. Схематическое изображение принципа действия гидродинамиче­
(195) ского подшипника

линейной стенкой установившееся (рис. 92). Определим подъем­ ную силу, действующую на колодку.

В реальных подшипниках скольжения толщина h масляного слоя весьма мала по сравнению с размером / колодки и незна­ чительно изменяется вдоль длины т. е. угол а мал. В связи с этим можем считать, что отношение hit стремится к нулю и

составляющая

скорости

ѵц

0;

будем

также

полагать,

что

в

направлении оси

z

движения масла

нет и

скорость

ѵг

=

0.

 

 

 

 

 

 

 

 

В силу этих предположений уравнение неразрывности для

несжимаемой

жидкости div v

=

0

примет вид

 

 

дх 0.

Следовательно, скорость ѵх может зависеть лишь от коор­ динаты у : ѵх = ѵх (у). Если пренебречь также массовыми силами, то уравнения (65) дви­ жения запишутся в виде

0 =

]_

 

д2ѵх

•-

4 г # - v

 

 

 

Р

дх

ду2

 

0

=

др_

др

 

 

 

 

 

ду

=

0.

 

 

 

 

 

Следовательно, как и для течения Пуазейля, получено уравнение, связывающее изме­ нение давления масла в масля­ ном слое с изменением скоро­ сти в этом слое:

1

dp

d2vx

(x

dx

dy2

Ось Вала

0К \ \ \ \ \ \ \ Ч ѵ

Вал

и0

где р = р (х).

Отметим, что градиент давления dp/dx целиком определяется

вязкостью масла

и. и производной dv2Jdy2,

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

только силу трения. Силы инерции

в

полученное

уравнение

не входят, хотя

скорость с д в и ж е н и я

масла не мала;

силы инер­

ции не рассматривались в силу сделанного предположения о ма­ лости расстояния h между колодкой и подшипником.

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

Обозначим скорость ѵх через U и проинтегрируем 2 раза уравнение (195) по у, тогда получим

187"

Обозначая

 

 

— b (.г), найдем

из

граничных

условий С±

и С2 . Граничные

 

условия

имеют вид:

1)

при у =

0 скорость

U = (70 ; 2) при

у

=

h (х)

скорость

(7 =

0.

Следовательно,

 

С2

— (70,

С! —

^

 

.

 

Подставим эти

значения

в формулу для (7, тогда

 

 

U

=

^f(y-h)--^-(y-h).

 

Через каждое сечение с высотой h проходит одно и то же количество масла

q=\udy

= -

12

' 2

ö

 

 

 

Отсюда находим величину

 

 

Ъ{х)

12Q

,

6(7,

h3

1

2

 

т. е. градиент давления можно представить в виде

J_

_dp_ _

6U0 _

12Q

j.i

dx

h-

h3

Поскольку h уменьшается с ростом x, перепад давления вдоль колодки будет не постоянным. Интегрирование от 0 до х дает:

X

 

p (x) = ро + 6u.<70 J - g - -

12д.(? J -£jL.

(196)

об

Характер

распределения давления р вдоль колодки

показан

на рис. 92.

Очевидно, что давление р 0 перед колодкой

на входе

в щель и за колодкой на выходе из щели должно быть одинаковым. Из условия, что р (х) = ро при x = / находим зависимость рас­ хода масла от величины зазора, которая является заданной величиной:

п — J^L. f -ÉL

V ~ 2 J Л2

Таким образом, вычислив расход Q, можно подсчитать рас­ пределение давления р (х) вдоль колодки по формуле (196). Подъ­ емная сила, действующая на колодку, определится выражением

Р = \[р{х)~

р 0 ] dx.

188

Можно также найти значение /г0 =• /г (л'„), при котором dp/dx = 0, т. е. определить сечение, в котором давление масла достигает наибольшего значения. Приравнивая dp/dx к нулю,

получим, что

^Я - =

0> т -

е -

ho — 2Q/t70 .

 

Сравним

это

выражение

с формулой для

расхода

Q, тогда

 

 

 

J

dxlh-

 

 

 

 

 

j

dx/h3

 

 

 

 

 

о

 

 

 

 

Отметим,

что

подъемная

сила

возникает

в первую

очередь

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

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