книги из ГПНТБ / Бекнев В.С. Газовая динамика газотурбинных и комбинированных установок учеб. пособие
.pdfВ результате первое уравнение примет вид
1 dp _ d-vx
цdx dz-
Проннтегрнруем полученное уравнение 2 раза по z, тогда
г2 dp |
, |
г |
г |
'°х ~~ 2ц dx |
' |
l Z i |
- |
Определим произвольные постоянные интегрирования С\ и С 2 .
Из граничных условий скорость ѵх= |
0 на стенке |
при z = |
±/г, |
|||||||
значит |
= 0 и С, |
= • — С л е д о в а т е л ь н о , |
скорость |
|||||||
|
|
|
|
_ |
(г» - /.») |
dp |
|
|
|
|
|
|
|
А |
' — |
2ц |
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 |
2ц |
V dx |
)J |
|
~ |
12ц |
dx |
|
|
—Ii |
|
|
|
|
|
|
|
Отсюда интегрированием |
по |
х |
в |
предположении, что при |
|||||
X = 0 давление р |
= |
р 0 , получаем |
искомую потерю давления, т. е. |
||||||
|
Ро — р = Ар = |
- ^ - |
и с р |
A.v. |
|
Следовательно, при установившемся ламинарном течении вяз кой несжимаемой жидкости потеря давления пропорциональна средней скорости течения, длине пройденного участка канала, вязкости жидкости, а также зависит от размеров канала. Этот результат хорошо согласуется с данными опытов при ламинарном течении.
180
Подобное течение в случае движения электропроводящей жидкости в поперечном магнитном поле называется течением Гартмана. Рассмотрим плоское стационарное ламинарное течение вязкой несжимаемой электропроводящей жидкости между двумя изолированными пластинами в поперечном внешнем магнитном поле. Будем полагать (рис. 87), что скорость потока ѵ направ лена по оси X, приложенное постоянное и однородное магнитное поле В — по оси z, а электропроводность а постоянна. Начало' координат находится посредине между пластинами, которые будем считать изоляторами. Примем, что 2/і много меньше 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). Определим подъем ную силу, действующую на колодку.
В реальных подшипниках скольжения толщина 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.