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

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

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

При использовании выражения (39) учтем, что

cos (п, у) = sin ср cos IT; COS (Л, Л') = sin ф sin •&; cos (л, z) — cos ф.

Подставляя все эти выражения в формулу (41), раскрывая •скобки и собирая члены при нормальных и тангенциальных на­ пряжениях, получим

я

р =

~

J dtp J (cry sin3 ф cos'2 ft - f аЛ. sin3

ф sin2 i3'-f-

 

 

u

б

 

-}-- а2 sin ф cos2

ф - j -

Л.у sin3 ф sin f>cos О -}- 2тt 2

sin2 ф cos ф sin ф -f-

-f- 2т, sin2 ф cos ф cos ft) dft.

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

напряжений ах< У і г

и тХУі XZi

y z в точке M (центр

сферы)

и счи­

тать их постоянными;

в результате получаем,

что

 

 

 

p =

_ o x

+ au + o,t

 

 

( 4 2 )

Таким образом,

за

величину

давления р

в

вязкой

жидкости

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

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

сательные напряжения равны нулю и тогда из

уравнения

(40)

получим, что

 

 

 

 

 

 

 

Px

= oJ,

Py = auj,

рг = аг7г.

 

 

Спроектируем

уравнение

(39) для рп

на ось х,

используя

полу­

ченное значение

для

рх:

 

 

 

 

 

—pncos(n,

х) = о"ѵ cos (я, х).

 

 

Знак минус в левой части стоит потому, что угол (п, х) тупой (рис. 31), а знак минус для cos (/?, х) справа был учтен ранее при выводе формулы (39). На другие оси проекции рп будут таковы:

—рп cos (п,

у) =

ау cos (п,

у);

 

—p;!cos(n,

z) =

a2cos(tt,

z).

 

Отсюда следует,

что

 

 

 

 

Рп

= — ст.ѵ =

—о,, = —ог =

р,

(43)

СО

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

ложено

и равно в любом направлении рп

= —рп.

от скорости

Связь

касательных

напряжений

т с производными

потока

может быть

принята как

опытный

факт для

многих (но

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

Определим касательные напряжения на гранях прямого «жид­ кого» угла, параллельных координатным плоскостям. Сначала рас-

5)

Рис. 34. Механизм появления парных касательных напряжений

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

делим в жидкости два нормальных

сечения тип

(рис. 34, а),

расположенные на малом расстоянии

b один от другого.

При движении плоскости А относительно плоскости В вправо сечение m за время At становится на место сечения т' и сечение п — на место сечения п'\ очевидно, что эти два сечения как бы проскаль­

зывают

одно относительно другого,

поворачиваясь вокруг то­

чек О и

О'.

 

Угловая скорость вращения этих

сечений со = - ^ - ( в момент,

когда сечения расположены нормально к движению жидкости, причем сечение m вращается вокруг точки О, сечение п — вокруг точки О'), а Д-у — угол, на который повернутся сечения m и п.

Относительная скорость скольжения сечения п по отношению к сечению m найдется по правилу механики. Остановим сечение m, для чего всей системе сообщаем движение, обратное движению сечения m, т. е. создадим вращение вокруг точки О в направлении против часовой стрелки. Тогда сечение m станет неподвижным, а сечение п (во всех своих точках) получит скорость, параллельную сечению т, равную ыЬ, т. е. сечение п скользит при этом вверх

61

по отношению к сечению m со скоростью ab. По формуле Ньютона касательное напряжение

ab

Угловую скорость со подсчитываем по скорости ѵм в точке М:

Подставляя это значение в формулу для т х , получим

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

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

рости

по координате,

 

 

 

дѵг

 

 

 

 

 

 

 

 

например —

 

 

 

 

 

 

 

 

Рассмотрим теперь общий случай, когда подвижны обе грани

«жидкого»

угла.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рассмотрим прямой «жидкий» угол AM В (рис. 34, б), который

спустя At секунд, переместится в положение А'М'В',

так

что его

стороны А'М'

и В'М'

образуют

малые

углы Ау1

и Ау2

с

перво­

начальными

направлениями.

 

 

 

 

 

 

 

 

 

 

Подсчитаем касательное напряжение на стороне MB,

которое

сложится в

результате влияния

двух величин:

 

 

 

 

 

1)

непостоянства

скорости

ѵх

вдоль

стороны

MA

 

(что

рас­

смотрено выше) и

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

поворота

стороны MB

(ранее

неподвижной)

с угловой

ско­

ростью AyJAt

=

dvjdx.

 

 

 

 

 

 

 

 

 

 

 

Угловая

скорость в первом случае dvj.dy

=

AyJAt.

 

 

 

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

касательное

напряжение

хІ/х на стороне

MB

 

 

 

х

-

A

V i + A v

.

- J

*>х

I

A

.

 

 

 

 

Очевидно, что касательное напряжение на стороне

MA,

 

рав­

ное хху,

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

скорости ѵи

по стороне MB

и поворотом

стороны

MA

с

угловой

скоростью

Ауг/Аі.

 

 

 

 

 

 

 

 

 

 

действуют

Следует обратить внимание на то, что напряжения

в одном направлении,

т. е. они

складываются.

 

 

 

 

 

62

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

т _ т _ „ ( д ѵ х I Д Ѵ * \ •

(дѵу

dvz

(44)

 

~2У — ~Уг Г \ д г - - Г

д у j

)

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

Рис. 35. Связь нормальных напряжений с полем скоростей

ным напряжением т, действующим на площадке, наклоненной под 45° к оси X. Выделим для этого призму с высотой h и основа­ нием ABC (рис. 35) и запишем условие ее динамического равно­ весия в проекции на направление AB.

Проекция силы инерции бѵдет равна

где

V = a2h/2 — объем

призмы;

ѵ — скорость ее

центра масс.

Проекцию плотности массовой силы на линию

AB обозначим

fmAB)

т а к ч т 0

полная

проекция

массовой силы

равна fmAäPV-

Проекцию поверхностной силы подсчитаем через отдельные ее

составляющие,

проектируя

их

на отрезок

AB — a/cos 45°. За

положительное направление

примем направление от Л к В. Тогда

 

 

xABh -f- ах a cos 45° h — <зуа cos

45°/г,

 

а остальные составляющие

сократятся.

 

 

63

Условие

равновесия

примет

вид

 

a-h

 

 

VI

V2

р — f"1

~ ("HF") \ AB

+ т а У 2 h

+ а * а " V 1 H

G ' J A " V 1 / l = °-

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

 

 

 

 

 

 

2т =

ау

— <тѵ. •

 

 

 

 

 

Рассмотрим теперь это же напряжение т

 

в

системе коорди­

нат х'у', повернутой относительно осей ху

на 45°. В координатах х'у'

обозначим его через хХ'Ѵ>,

тогда в соответствии с ранее

полученной

формулой

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

X

=

 

Хх'у'

=

д

дѵг,

 

 

 

 

дѵу,

 

 

 

 

 

 

 

 

ду'

 

 

 

дх'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Вторым этапом вывода будет переход от скоростей в системе х',

у ' к скоростям в системе х, у . Из рис. 35 следует,

что

 

X =

х'

cos 45° — у'

cos 45° =

 

(х'

 

у')

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(45)

у =

х'

cos

 

45°

+ у'

cos 45°

 

=

-4f-(x'

+

 

у').

 

Отсюда

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Х'

=

^

( х

+

у);

 

у'

 

=

 

^

( у ~

Х

)

,

 

так что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dx'

 

V2

,

,

ч

•%-

= vy.=^(v„-vx).

(46)

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

По правилу дифференцирования сложной функции получаем,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дѵх.

 

 

 

дѵх.

 

дх

 

L

&>х>

.-ËL

 

 

 

 

 

 

ду'

 

 

дх

'

дц'

 

 

 

ду

 

 

 

 

 

 

 

 

1

 

 

ду'

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а ѵ

 

__

д»„.

 

дх

,

 

 

даУ

ду

 

 

 

 

 

 

дх'

 

 

 

дх

дх'

'

 

 

 

ду

дх'

 

 

 

 

Делаем подстановки, пользуясь формулами

(45) и

(46):

дѵх'

 

 

 

1 ( дѵх _|_ àvtJ

\ у

 

\ ( дѵх

 

,

дѵу

 

ду'

 

 

 

2

V дх

'

дх

J

 

1

 

2

V ду

 

1

ду

 

à°u'

_

1

/

дѵѵ

 

дѵх\ .

 

 

1

/

дѵу

 

QOX

 

дх'

 

2

 

\

дх

 

дх

 

 

 

2

\

ду

1

 

ду

 

64

Сложив оба равенства и умножив их на \і, получим

Т = (.1

 

дѵх

 

дх

ду

 

г,

I дѵх

дѵ„

""у

Иначе

Если аналогичные выкладки сделать относительно пло­ скости хог, то получим

, = * ( • & - £ ) .

Запишем тождество

и сложим три последних равенства, тогда

За, - К + ои + аг) = 6ц

 

 

дѵ„

дѵ.

-

2 ^

1( ^^f + ~ ô ^ + dz

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

дѵ

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

<*х = Р +

2

 

 

- О -

р. divü,

 

öüj,

2

ndivu;

(47)

-P + 2 ( . i - ^ - - - g -

ду

 

 

 

vz = —Р + 2|-1

- 4 "

d i v

 

Для идеальной жидкости, лишенной сил трения, коэффициент р, = 0; тогда, естественно, получаем, что

Р = —сх = —Oy = —ог.

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

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

5 B . C . Бекпсв

6S

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

pnds на ось X. На рис. 36 показаны напряжения, дающие проек-

S

ции только на ось х.

Рис. 36. Поверхностные напряжения, действующие по оси х

На гранях, проходящих через начало координат, это будут напряжения ох, тух, xzx. Полагаем, что на других гранях напря­ жения могут быть выражены через эти составляющие исходные путем разложения их в ряд Тейлора (до второго члена) в окрест­ ности начала координат. Например, для грани A BCD имеем на­ пряжение

 

 

 

 

°Ч-ѵ

+1

- H -

dx;

 

 

 

 

 

 

 

 

-

 

 

дх

'

 

 

 

 

члены

более малых

порядков

отбрасываем.

 

 

 

Составляющая вдоль

оси

х

главного

вектора

поверхностных

сил

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

 

 

 

 

 

 

 

 

 

 

(-ах

+ ах + Щ%- dx)

dy dz +

( ух

+

xyx + ^f-dy)

dx dz +

+

 

+ T M +

dz)dxdy

 

=

( - f

+

*gL +

*e)dxdydz.

66

Предположим, что действие главного вектора поверхностных

сил можно заменить

главным

вектором

объемных

сил ^ fv dV

 

 

 

_

 

 

 

 

 

V

 

 

с объемной плотностью силы fv

(fVx,

у , /ѵг ). Составляющая

вдоль оси X главного вектора объемной силы

равна fyx

dx dy

dz.

Приравнивая силы, получаем формулу для

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

 

fyx

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

получаются

аналогично):

 

 

 

 

 

дах

I

 

1

dtzx

\

 

 

 

 

 

х-

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

дх

~~>

ду

1 .

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

да у

 

дх2у

 

 

 

 

 

 

 

у-ху

 

»»у

+ v-zydz

. I

 

 

 

4 Ö

 

 

 

 

~W

 

 

Ï

 

 

(

 

)

ІѴу — -я7дх--Г-яП-+-яГ'

 

 

 

 

 

hz

= дх

 

dTyZ

 

âaz

 

 

 

 

 

 

^

ду

 

dz

 

 

 

 

 

 

В векторном виде с учетом формул (40) получим

 

 

 

 

1ѵ = fvJ +

fvyj+

fvzk =

- | f +

+

(49)

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

^ху — ^ух]

Xvz =

T'zx'i T-yz ~

^zyi

 

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

напряжений

выполнено.

 

В качестве примера рассмотрим как выражаются

поверхност­

ные силы через объемные

для

идеальной жидкости,

лишенной

вязкости.

В этом случае все касательные напряжения равны нулю

и поэтому

из выражения (48) для плотности объемной силы получим

Из равенства

(43) имеем

 

 

b =

--%r~i-^-l--%-k

=

-&*àp.

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

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

5*

67

Пондеромоторная сила электромагнитного

поля, действующая

на движущуюся со скоростью Vk k-ю частицу

с зарядом qk,

Fk = qkË + q k ( ü j X ß ) ,

 

где E—• вектор напряженности электрического поля; В — вектор индукции магнитного поля.

Чтобы получить объемную пондеромоторную силу AFy, дей­ ствующую на объем А V с массой частиц Am, надо просуммировать силу Fh по всем частицам рассматриваемого объема

А

Согласно определению, вектор плотности объемной пондеромоторной силы электромагнитного поля

AF

 

+ 1 lim

 

/„ = lim - ^ -

= E lim - ^ -

 

-±—-)XB.

ДѴ^О Л Ѵ

АѴ->0

Л Ѵ

W - > 0

а |

/

В данном случае величина

ре =

lim

является объемной

 

 

 

АѴ->О Й Ѵ

 

_

плотностью электрического заряда, а величина

у =

lim — гт;

вектором плотности электрического тока. Таким образом,

 

}у=РеЁ

+

(]хВ).

 

(50)

При условиях, характерных для рассматриваемых магнитогидродинамических устройств, член реЕ, зависящий от электри­ ческого поля, существенно меньше члена (/ х В), зависящего от магнитного поля (отсюда и название магнитная гидрогазодинамика; если учитывать член реЕ, то рассматривается электромагнитная гидрогазодинамика).

Итак, считают, что плотность пондеромоторной силы электро­ магнитного поля

 

/„ = fx

В или

]

т

=

-L(fx

В),

(51)

а с учетом соотношения

Максвелла

/

=

rot Я

 

 

] ѵ

= (rot H) X В=

(rot В) X В = \xma

(rot H) X H.

 

 

И/на

 

 

 

 

 

 

 

Возможность отбрасывания члена р е £

 

по

сравнению с членом (/ХБ) в фор­

муле для пондеромоторной силы fv покажем путем оценки

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

этих членов.

По уравнению

Максвелла

заменим ре =

div

D, тогда

 

/і/ = e a

£ d i v £

-I

( r o t S ) X ß .

 

 

 

 

 

И"іа

 

 

 

 

68

Каждую размерную величину заменим через произведение безразмерного ее значения и выбранного масштаба. Например, проекцию вектора В на ось х запи­ шем так: ß v = В0Вх- В данном случае масштаб В0 определяет размерность Вх и ее~абсолютное значение, а безразмерная величина Вх является функцией коор­ динат x, у, z и времени t и определяет закон изменения величины Вх в простран­ стве и времени (покажите, что В — В0В при одном и том же масштабе В0 на каж­ дой координатной оси). Аналогично обозначим:

е а бо^» В

EQE

,

и

т. д.

Mma = Мо(А.

х =

1ох<

У = ІоУ, 2 =

V i

единичные векторы і, j , k, входящие в выражения для векторов, изменений не

претерпевают. При таких

подстановках

первый член силы /у примет вид

Biß div E =

e a £

dEx

dEu

 

_дЕл

 

дх

1 ду

 

dz

 

 

 

 

 

 

д(Е0Ёх)

д

(E0ÉU)

д (Е0Е2

J - и

еЁ div Е.

д(10х)

 

д(10у)

 

d(l0z)

 

 

 

 

 

Итак, при переходе к безразмерным переменным возникает масштабный ком­

плекс е 0 EQ/10, который определяет абсолютную

величину первого члена. Закон

функционального изменения первого члена по отношению к безразмерным пере­

менным останется тем же.

Во втором члене векторное произведение структурно представляет собой

сумму произведений

величин

двух

проекций,

например

 

(rot

В)х

Ви.

Подстановка

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

переменных

дает

 

 

Ита

 

 

 

в этом случае

 

 

 

 

1

(rot

В)х Ву

 

1

дВг

 

дВу

Ву

=

 

 

 

 

Ига

Игаа

ду

 

dz

 

 

 

 

 

 

 

 

 

 

 

 

 

1

"

д(В0Вг)

д (В0Ву)

 

 

 

 

(rolB)xB

 

 

 

Иои

L 5 (/„(/)

 

 

 

 

 

 

 

 

И

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

надо

сравнить

между

собой

масштабный

комплекс

 

Ч

е„ -г—

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Воспользуемся

обобщенным

законом

Ома / = аЕ -j- a (v X

В),

 

пред­

Ио'о

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

полагая ради простоты, что токи Холла не существенны, и соотношением rot

 

# = / .

Исключив величину

/, получим

уравнение

 

 

 

 

 

 

 

 

 

Е =

1

 

rotfi —

(vXB),

 

 

 

 

 

 

 

 

öUma

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

которое связывает поля Е и В и дает возможность оценить масштабы Еа

и В0.

По порядку величины поле Е равно или порядку первого члена справа

_ \

X

Вп

 

 

 

 

 

 

 

 

 

 

 

 

О^оИо

 

 

 

члена

ѵ0В0.

Это утверждение

остается верным

X —j—, или порядку второго

рассматриваемых

порядков:

 

 

 

 

 

 

 

 

и для квадратов

 

 

 

 

 

 

 

 

 

 

 

 

 

 

VQBQ

 

 

 

 

 

 

 

 

 

 

 

или

 

 

 

 

 

Hcfao'o

69

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