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

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

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

потока растет, а на участке отвода теплоты — падает, что непо­ средственно следует из уравнения (139).

Полное давление р* при подводе теплоты уменьшается (в этом состоит эффект теплового сопротивления), а при отводе теплоты — растет. Это может быть показано на основе рассмотрения полного

импульса Ф потока,

который

постоянен

в

любом сечении трубы

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

не действуют никакие

силы.

 

 

 

 

 

 

 

Рассмотрим два

соседних

сечения,

для

которых уравнение

импульсов

через газодинамические функции

имеет вид

 

 

Фг

=

Ф2

 

 

 

или

 

 

 

 

 

 

 

Так как

поперечное сечение sx

= s2,

то

 

Характер изменения функции / (к) (см. рис. 42)

таков, что

для к <С 1 величина f

(к)

возрастает с увеличением к,

а для к > 1

с ростом

к уменьшается

f (к). В

дозвуковом

потоке

подводом

теплоты)

при

кх <[ к2

отношение

рг/рі < 1,

а полное

давление

вдоль трубы

снижается.

 

 

 

 

 

В сверхзвуковом потоке (с отводом теплоты) при kt

>• Х2 отно­

шение рУрі >

1 и полное давление вдоль трубы растет.

 

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

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

Обозначим через Ѳ = T*DUJTDX степень подогрева газа и опре­ делим предельно допустимое число квх. к р на входе в трубу, вызы­ вающее критическую скорость на выходе из нее для фиксирован­ ного значения 0. Запишем отношение коэффициентов к для вход­ ного и выходного сечений:

^вых

ѵвых

акР-

вх

квх

ѵвх

акр-

вых

170

Используя уравнение неразрывности, состояния и газодина­ мическую функцию т (X) — Т/Т* для k = const, представим отно­ шение скоростей для цилиндрической трубы в виде

 

 

 

 

 

 

^вых

g

т (Ä-вых)

Рвх

 

 

 

 

 

 

 

 

 

 

 

 

 

t>Bx

 

 

т а в

х )

 

F вых

 

 

 

 

 

 

 

Используя

газодинамическую

функцию

 

г (X) — ps/Ф

и

по­

стоянство

полного

импульса

Ф вдоль трубы,

получим

 

 

 

 

 

 

 

 

 

 

Рвх

_

r(Xßx)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Рвых

f

(^вых)

 

 

 

 

 

 

 

Отношение

aKp

в х / а к р і

В Ь | Х

=

1/]/Ѳ.

Учитывая

конкретную

за­

висимость

функций г (X) и т (X) от

величины

X,

найдем

оконча­

тельно,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

^вх

 

 

1 +

с

т/Ѳ.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

1 +

^ в х

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Откуда при Хвых =

1 полу­

 

 

 

 

 

 

 

 

 

 

 

 

чаем

квадратное

уравне­

 

 

 

 

 

 

 

 

 

 

 

 

ние

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

"вх. нр "

2 ] / 0 Х в х . к р + 1 = О,

 

 

 

 

 

 

 

 

 

 

 

 

решение

 

которого

 

 

 

 

 

 

 

 

2

 

 

 

 

 

 

^ х . к Р

=

> / Ѳ ± / 9 ^ Т

 

 

 

 

2

4

(

 

 

6

8

g _ К,,

 

 

 

 

 

 

 

показано

на рис.

82.

 

 

 

 

 

 

 

 

 

 

 

 

 

Ts,

 

 

Рис. 82. Зависимость предельно допустимой

Рассматриваемый

слу­

приведенной

скорости

на

входе в трубу от

чай движения газа не пред­

 

 

 

 

степени

подогрева:

 

 

 

ставляет

собой

энергоизо-

/

— с в е р х з в у к о в о й поток;

2

д о з в у к о в о й

поток

лироваиный процесс.

По­

 

 

 

 

 

 

 

 

 

 

 

 

казатель

 

политропы

и

этого процесса не

постоянен, что

видно

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

уравнений

(164)

и (165).

Если

эти

уравнения

разделить

одно

на другое,

то

получим,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

d In p

 

: Ш 2 .

 

 

 

 

 

 

 

 

 

 

 

 

 

 

П = d In p

 

 

 

 

 

 

 

 

Проанализируем

качественно

ход процесса в Ts-диаграмме

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

до сверхзвукового

(линии

Релея).

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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

числах

М) значение

п > 0,

т. е.

процесс

близок к

изобарному

(рис. 83). Ускорение потока вызывает увеличение показателя

политропы

п,

который становится равным единице при

M =

= l/^/fe,

что

соответствует изотермическому процессу

вблизи

этого режима. В момент достижения потоком критической ско­ рости показатель политропы п = k, а процесс происходит адиаба­ тически. При дальнейшем ускорении потока необходимо заменить

171

подвод теплоты на отвод, что соответствует протеканию процесса с уменьшением энтропии s потока. Ускорение потока до больших чисел M приближает процесс к изохорическому с п —> оо.

На этом же рисунке обозначен

процесс торможения потока

при тепловом воздействии.

 

Разгон

Торможение

T=const,n=i

Отвод теплоты

 

Подвод теплоты

 

Отвод теппотьі

 

 

Подвод теплоты

Рис. 83. Линии Редея для течения в трубе с теплообменом

§ 22. ОДНОМЕРНОЕ ДВИЖЕНИЕ ЭЛЕКТРОПРОВОДЯЩЕГО ГАЗА ПРИ ВОЗДЕЙСТВИИ ЭЛЕКТРОМАГНИТНЫМ ПОЛЕМ

 

Рассмотрим

изолированное магнитогидродинампческое

воздей­

ствие на поток

электропроводящего

газа, понимая под этим

ста­

ционарное квазиодномерное движение электропроводящего

газа

 

 

 

 

в канале постоянного сече­

 

 

 

 

ния

при отсутствии

тре­

 

 

 

 

ния,

теплообмена

с внеш­

 

 

 

 

ней средой, изменения рас­

 

 

 

 

хода газа и подвода или

 

 

 

 

отвода механической энер­

 

 

 

 

гии. Газ при своем движе­

 

 

 

 

нии

пересекает

 

внешнее

 

 

 

 

электромагнитное

поле Е

 

 

 

 

и В (рис. 84), а

через

газ

 

 

 

 

течет

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

ток

Рис.

84.

Схема течения проводящего газа

плотностью /.

 

 

 

_

_

в скрещенных полях Е и В

Будем для

определен-

_

 

ности

полагать

 

векторы

Е, В и v взаимно перпендикулярными

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

в соот­

ветствии с рис. 84, так что

 

 

 

Е=Ё)орп

B = Bk,

~ѵ = ѵі.

(167)

172

Магнитогндродинамическое воздействие — воздействие не эле­ ментарное. В электропроводящей среде, имеющей конечное зна­ чение электропроводности ст, оно слагается из двух неразделимых факторов: джоулева нагрева и действия пондеромоторной силы. Джоулев нагрев является необратимым воздействием, а действие

пондеромоторной силы fm = -^-(jxB)

обратимо и зависит от на­

правления векторов магнитного поля В и электрического тока /.

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

] = У/орті

а направление действия пондеромоторной силы [,п совпадает с направлением скорости ѵ потока (табл. 3). Если изменить на­ правление электрического тока / на обратное, то направление действия пондеромоторной силы /,„ будет противоположно ско­ рости v, что соответствует каналу МГД-генератора.

3. Схематическое изображение электромагнитного воздействия

на поток

 

 

в канале

постоянного

сечения

 

 

Х а р а к т е р

Д о з в у к о в о й

поток

M < 1

С в е р х з в у к о в о й

поток

M > 1

воздействия

 

 

 

 

 

 

 

dv < 0

MТДнасос

-

МГД-генератор

 

0

3*

'/

1

J

 

Торможение

1

 

 

 

 

 

 

1

V

1

 

fm ~~ *

"

Ѵ

 

 

 

 

г—

~1

 

 

 

МГД-генерашор

 

МГД-насос

 

Ускорение

1

 

1

 

!

1,

 

dv>

0

/ —

 

V

 

 

 

 

 

 

 

 

 

 

1

 

1

 

!

1

 

Для вывода уравнений воздействия удобнее обратиться к ис­ ходной системе (27), (137), (138), отбросив члены, не учитывающие электромагнитное воздействие:

° - Ѣ — - т & + ^

( 1 6 8 )

pvs = G; р = - ^ г -

 

173

Согласно рис. 84,

Im = 7 ( 1 ' Х ^ ) = 7 ßl = îmJ

и направлена по движению газа. С учетом этого мощность

Величина Q,n представляет собой джоулеву теплоту, выде-

1 у*

ленную в единице массы газа и равную — -^-.

В результате уравнения (168) и (169) примут вид

ѵ

4

-

=

-

(

1

7

°

)

dx

 

 

p dx

'

р 1

 

 

 

ѵ

'

dT ,

d

( v2

\

1 . 0

I

1 j -

 

 

, w n

Исключим из уравнения (171) энергии механические величины, для чего вычтем из него уравнение (170), умноженное на ѵ. В ре­ зультате получим

°PV dx ~~ +

p dx ~т~ p a " '

Если уравнения расхода и состояния записать в дифферен­ циальной форме (учтя, что dm = 0, ds = 0, где s — площадь сечения потока), то окончательно будем иметь следующую систему уравнений:

" 4 L

=

L

4 L

+

— / В ;

 

(172)

dx

 

p

dx

1

p 1 '

K

'

- ^ =

- L

- ^

+1

A

4

(173)

d.v

pCp

dx

 

opcp a

 

'

» - £ + » т = °;

 

( ' « )

Выразим изменение величин v, p, T, p под влиянием электро­ магнитного воздействия. Для этого используем соотношения:

cp — cv = R, -^- = k,

Р

Для определения изменения ѵ необходимо подставить в урав­ нение (175) значение dp/dx из уравнения (174) и dTldx из уравне-

174

ния (173). Затем значение dpldx, найденное из уравнения (175), следует подставить в уравнение (172), в результате получим

( M » - 1 ) J *

+

 

( 1 7 6 )

v

; V dx

va2 pa 1 о2 p

v

'

Чтобы найти изменение величины р, необходимо использовать

уравнение (172), подставив в него

значение dv/dx из

уравне­

ния

(176):

 

 

 

 

( M 3

_ ] ) ^ ^ ( ^ _ l l z ^ _ . ß _

( 1 7 7 )

Для определения

изменения р

надо использовать

уравне­

ние

(174), подставив

в него значение

dv/dx:

 

( M ' - D - g — ( 1 7 8 )

Чтобы найти изменение Т, нужно использовать уравнение (173), подставив в него значение dpldx из уравнения (177), тогда

ѵ ' dx Сррѵ a рср1 4 '

Следует отметить, что аналогичное соотношение можно полу­

чить для

производной dM/dx [соотношения

(176)—(179)

можно

получить

и из выражений (147)—(150), если

величины

dQm,BK,

и dLm, мех трактовать согласно § 17].

Рассмотрим качественные результаты, следующие из системы уравнений (176)—(179). Первый член справа представляет собой необратимое воздействие вследствие джоулева тепловыделения. Характер этого действия аналогичен действию подводимой извне

теплоты

dQrlltBli,

разобранному

в § 21. Например, в дозвуковом

потоке,

где

М <

1, джоулево

тепловыделение ускоряет

поток

(dv >> 0), а

в сверхзвуковом,

где M >• 1 тормозит (dv

<С 0).

Принцип обращения воздействия к данному воздействию не при­ меним в силу его необратимости.

Второй член справа описывает обратимое воздействие от пондеромоторной силы (jxB)x. Знаки перед вторыми членами соответ­ ствуют рассмотренному на рис. 86 движению электропроводящего газа в МГД-насосе, для перехода к режиму МГД-генератора необ­ ходимо их изменить на обратные, что соответствует изменению на обратное направление электрического тока /'. На рис. 87 дана таблица возможных воздействий пондеромоторной силы. Следует обратить внимание, что разгон и торможение потока не всегда соответствуют действию пондеромоторной силы fm вдоль или соответственно против потока.

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

чительно

джоулево

тепловыделение, т. е. при очень больших

значениях

чисел Rem

= av\iml.

175

Для рассмотрения количественных зависимостей необходимо добавить к системе соотношений (176)—(179) уравнение (94) индукции.

В случае стационарного одномерного движения вдоль оси д уравнение (94) принимает вид

d (vB)

d*-B

,, о л ч

Уравнение (180) показывает, что внутри движущейся жидкости величина электромагнитного поля не равна значению приложен­ ных извне полей Е и В, а целиком определяется уравнениями Максвелла.

Рассмотрим в дальнейшем более простой, но технически важ­ ный случай движения потока с малыми числами Rem . Это озна­ чает, что внутри движущейся жидкости можно пренебречь наво­ димыми магнитными полями вследствие текущих в ней токов При этом уравнения Максвелла игнорируются, т. е. игнори­ руется связь между полями Е п В и электрическим током /. Счи­ тается, что электрическое Е и магнитное В поля в жидкости тожде­ ственно равны внешним приложенным полям Е и В, которые вдоль оси X можно задавать произвольным образом с помощью технических устройств. Величина электрического тока в жидкости определяется только законом Ома [формула (8)J.

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

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

_ Tj

Эта величина может меняться вдоль оси х, так что следует рассматривать местный электрический к. п. д. г|. Для условий, изображенных на рис. 84, было получено, что

~fvv = jBv.

176

Согласно соотношениям (167) получаем, что Ej — Ej и, сле­ довательно,

Е

В МГД-насосе Ej затрачивается на джоулево тепловыделе­ ние /2/ст и преодоление работы fvv, т. е.

В МГД-канале баланс рассмотренных величин таков:

 

 

 

 

 

}ѵѵ

=

 

 

 

 

 

 

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

МГД-насоса электрический

к. п. д.

 

 

 

 

 

 

 

 

Л =

=

 

 

 

 

 

 

 

 

 

 

 

 

 

fvv

 

 

 

 

 

 

всегда

больше

единицы, а для МГД-генератора

величина

 

 

 

 

 

 

 

 

4

=

 

 

 

 

 

 

заключена

между

единицей

(когда

величина

fvv =

Ej

и

нет

джоулева тепловыделения) и нулем (когда вся

работа

fvv

=

/ 2

затрачивается на джоулево тепловыделение).

 

 

 

 

Исходя из векторной записи закона Ома, согласно обозначе­

ниям

рис.

86,

имеем

 

 

 

 

 

 

 

 

 

 

 

j = jy =

оЕ

аѵВ =

оѵВ

1).

 

 

(181)

При т) >

1 ток

/ >

0

и

расположение векторов на рис. 84

действительно

соответствует

каналу МГД-насоса; при

0 < г )

•< 1

ток /

< 0 , как

это

и должно быть

для

канала

МГД-генератора.

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

Учитывая, что а2 =

получим выражение

 

(M» - 1 ) - g -

= -

<|І (v -

Wl) (v - wu),

(182)

в котором комплексы шх

и w.2

имеют

размерность

скорости:

k1

 

Е

 

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

12

в. С. Бекнев

177

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

временно с

равенством M = 1 выполнялось

одно из

условий

и = Wx или

V = w2. Использование этого и остальных

уравне­

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

изменения

полей Е

и В вдоль оси X, а также построить диаграмму (рис. 85),

наглядно

иллюстрирующую возможное сочетание параметров при электро­

магнитном

взаимодействии.

 

 

 

 

 

 

 

 

 

Отметим, что на диаграмме вдоль линии бб приращение числа M

равно нулю

(d.M. =

0). Направления

возможных переходов через

 

 

 

 

 

M = 1,

в точках

Пх и Я 2 ,

разре­

 

 

'%

 

шаемые уравнением

(182), обозна­

3,5

1

д

чены

стрелками.

 

 

 

 

 

 

 

По

 

заданным

исходным

пара­

 

 

У

 

метрам

на входе

в трубу

опреде­

 

 

 

 

ö

ляется

 

пара

значений M

и

vlwx-

 

 

 

 

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

того,

в

какую

 

 

D

Б

 

область

диаграммы

они

попадут,

 

 

 

 

 

возможно

различное изменение

 

 

 

 

 

параметров

потока при его движе­

снии вдоль оси трубы X.

 

 

 

 

)

 

 

 

 

Так,

переход

через

скорость

 

г *

 

 

 

звука

(М = 1)

с разгоном

потока

 

 

 

 

 

А,

 

 

 

 

 

 

 

возможен

только

из

областей

Аг

 

 

 

 

 

 

 

и Сх- Переход через скорость

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

звука

с

непрерывным

торможе­

Рис. 85.

Схематическое изображе­

 

нием потока возможен только из

ние характера

 

изменения

парамет­

областей

Аг

и

С2

-

Движение

ров в канале

постоянного

сечения

вдоль

оси X трубы

соответствует

при электромагнитном

воздействии

на

диаграмме

движению

по

на­

 

 

 

 

 

 

 

 

правлению указанных стрелок.

 

По оси ординат рис. 85 отложено значение

 

 

 

 

 

 

 

 

 

 

V

k

 

ѵВ

_

k

l

 

 

 

 

 

 

 

 

 

 

 

Wx

ft 1

E

 

Г Tj

"

 

 

 

 

 

Для

МГД-насоса

r\ > 1

(при

значении

k = 1,4),

что соот­

ветствует

области значений

ниже линии дд; для

МГД-генератора

0 - O l <

1, что соответствует области значений

выше

линии

дд.

Таким

образом,

непрерывный

переход

через

скорость

звука

в точке Пх требует изменения

знака

воздействия.

Непрерывный

переход через скорость звука в точке П.г

изменения

знака

воздей­

ствия не требует (в этом проявляется

неэлементарность

МГД-

воздействия).

 

 

 

 

 

 

 

 

 

 

 

 

 

Возможно также движение потока в пределах данной длины

трубы и без перехода

через скорость звука. Возможные

варианты

весьма разнообразны и нами не рассматриваются. Уместно только отметить, что при интенсивном джоулевом тепловыделении можно, разгоняя поток, уменьшить число M вдоль оси х. Это объясняется

178

•///////

тем, что в формуле M = vlYkRT знаменатель растет более интенсивно, чем числитель.

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

чается из уравнений (177)

и (178):

 

 

_ d\nP

_

А (fe I ) (л — 1 ) M 2

— 1

 

dlnp

 

(k 1) (ii — 1) — 1 "

§ 23. ТЕЧЕНИЯ

ПУАЗЕЙЛЯ

И ГАРТМАНА

 

Рассмотрим плоское стационарное ламинарное течение вязкой

несжимаемой

жидкости вдоль оси х между

двумя бесконечными

параллельными пластинами (рис. 86). Подобное течение назы­ вается плоским течением Пуазейля. Анализ течения проведем, используя уравнения движения (65) и неразрывности (26); при этом массовыми силами пренебрегаем. Полагаем, что составля­ ющие скорости ѵу и vz

вдоль осей у и z отсутст­ вуют. Тогда из уравнения

неразрывности следует, что

 

дѵхІдх =

0 и скорость

ѵх

 

вдоль оси x

остается

по­

 

стоянной. Далее, согласно

 

рис. 86, видим, что силы

 

вязкости

обуславливают

Рис. 86. Схематическое изображение течения

изменение

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

Гартмана

ѵх вдоль

оси

2, т. е.

ѵх

 

является функцией от z : ѵх — vx(z). Вдоль оси у величина ѵх изменений не претерпевает вследствие симметрии течения.

Проектируем уравнение движения на координатные оси. В силу сделанных выше предположений в них остаются только следующие члены:

0 =

1

др

д°-ѵх

р

дх

 

 

 

 

0 =

1

dp

 

Р

ду

 

 

Из этих уравнений заключаем, что давление р постоянно вдоль осей у и z и является только функцией х, т. е. р = р (х). Отметим, что градиент дрідх = const, так как мы предположили, что ѵх не зависит от х.

12*

179

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