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

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

.pdf
Скачиваний:
6
Добавлен:
23.10.2023
Размер:
7.63 Mб
Скачать

ка сумме скоростей среды и и скорости звука а в данной точке. Чем больше давление в рассматриваемой области, тем больше и + а. Давление среды в точке 1 больше, чем в точке 2. Поэтому «i + a i> ti2+ Ö2 - Если рассматриваемая волна сжатия является зву­ ковой и амплитуда ее очень мала, то разница скоростей движения возмущения от точек 1 и 2 также невелика и форма волны по мере'

 

С

С

f -

 

"1—

 

у'

/ £*

 

 

2 1 .

~Р'

 

 

 

в о л н а с я с а т а я

 

г,

 

Н е в о з м о ж н а я д з о р м а

 

д в и ж

е н и

я

У с т о й ч и в а ?

с р о р м а

д в и ж е н и я ( у д а р н а я в о л н а )

Рис. 5-3.

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

круче,

а участки спада давления — все положе

(на рис. 5-3 кри­

вые А

II В). Деформируясь, волна сжатия как

бы захлестывает

сама себя. Это явление аналогично тому, которое наблюдается-при набегании морских волн па пологий берег, когда основание волн тормозится трением о грунт, а гребни движутся с большей ско­ ростью. Для сильных звуковых волн выбегание гребней вперед при­ водит к изменению тембра звука: появляются гармоники более высоких частот. Чем больше амплитуда волны сжатия, тем больше разница величин и + а для различных точек волны и тем быстрее происходит такое изменение ее формы. Выбегание гребней вперед и отставание впадин не может продолжаться бесконечно, иначе форма волны должна была бы принять вид, изображенный кри­ вой с на рис. 5-3. Пунктирная кривая, соответствующая этому слу­ чаю, определяет неоднозначную зависимость давления в волне, т. е. в одной и той же точке пространства должно одновременно наблюдаться несколько значений давления. Это физически невоз­ можно. Поэтому изменение формы волны сжатия заканчивается тем, что участок кривой, соответствующий подъему давления, пре­ вращается в вертикальный отрезок, т. е. образуется резкая перед­ няя' граница зоны сжатия или фронт ударной волны, распростра­ няющаяся со скоростью, большей скорости звука в невозмущенном газе. При дальнейшем распространении волны ее структура с рез-

50

RUM фронтом остается устойчивой. Если даже по какой-то причине фронт ударной волны разрушится, например повышение давления от р ] до р2 произойдет не одним, а двумя скачками, то второй ска­ чок, движущийся в среде, сжатой и приведенной в движение пер­ вым скачком, догонит его и снова образует единый фронт.

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

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

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

§5-2. Основные уравнения для фронта ударной волны

'Рассмотрим картину движения фронта ударной волны (рис. 5-4). Пусть невозмущенная среда характеризуется давлением р0, массо­ вой плотностью Ро, температурой Т0 и скоростью Uо=0. Параметры газа «а фронте ударной .волны обозначим соответственно рф, рф,

иПф, а скорость фронта по отношению к иевозмущенной среде D. За время Ат фронт волны перемещается на расстояние ПДт, а ча­ стицы газа, находившиеся в начальный момент времени то во фрош те волны, — на расстояние Е7фДт. Среда, занимавшая при т0 слой

невозмущенного

газа толщиной

ДДт,

займет

слой толщи­

ной (D — 6ф) Дт.

Давление среды

в' этом

слое увеличится от ро

до Рф, плотность

от Ро до рф, температура

от Т0 до

7 ф, скорость

от 0 до £/ф. Будем записывать соотношения для единицы площади фрощ*а волны. Масса среды в слое, толщина которого изменилась от ПДт до (D — Ц,) Дт, осталась неизменной:

p0DA~ = рф (D — Uф) Д т.

4*

51

Таким образом, имеем уравнение сохранения массы или нераз­ рывности:

РоД = p,i (Д — //ф).

(5-1)

или

 

Р ф — P u

(5-Г)

Uф = D

Масса вовлеченного в движение газа составляет р0/Мт, а при­ обретенное этой массой количество движения вследствие увеличе­ ния скорости от нуля до 0ф равно р0£)Дт£Уф. Это изменение коли­ чества движения произошло под действием разности давлений рф— р о за время Дт. Поэтому уравне­

Яние изменения количества дви­

 

 

 

 

S.

жения

может быть записано так:

 

 

 

 

 

 

 

 

 

 

 

 

 

Т„

 

РоОі/фДL== ( Рф

Ра) Ат,

 

 

 

 

 

или

 

 

 

 

 

 

 

 

 

 

 

PoDL/ф^Рф — р().

(5-2)

 

 

 

 

 

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

кинети­

 

 

 

 

 

ческой и внутренней энергии слоя

 

 

 

 

 

газа толщиной /Мт. До возмуще­

 

 

 

 

 

ния

газ в этом слое

находился

б покое и обладал только внутренней энергией Ео:

 

 

 

 

Д0 =

Д Д0РоD^- = 6'

°г>

Ро РоДдэ

 

 

 

'Г'

/^0

 

 

' ' р

. «I

 

 

где

внутренняя энергия единицы массы.

 

 

 

 

 

После возмущения энерпня рассматриваемого слоя Дф будет

равна сумме внутренней энергии Е и кинетической Дкш,:

 

 

 

Д = ^ДфроDAт =

V

 

• — РоДАт;

 

 

 

 

 

 

 

 

 

Р ф

 

 

 

 

Дкпн= РоДДх^

 

 

V

 

 

 

 

 

 

 

 

 

 

Полная энергия рассматриваемого слоя газа составит

 

 

/ 7

___ / 7

I

С 4 ____

 

рф

. Д ф

 

 

 

 

 

“ Г

^ кііи

Сѵ

Р ф

Ь ДГ РоД^-

 

 

 

 

 

 

 

 

 

 

52

Изменение энергии слоя произошло в результате работы внеш­ ней /?ф Силы на расстоянии 1 7 . Поэтому будем иметь уравне­ ние энергии в виде

или

Еф

Д, = РфОф^1-,

 

 

Рф , Щ _

Zu Р »О = рфі'ф.

' й ,

9

P o

 

 

Подставляя в последнем уравнении к = — и производя неслож­ ные преобразования, получаем

1 ( г

- — )р»0 = Р^и ф

- ^ ° ио--

'(5-3)

К' — 1 Р ф

Р о /

-

 

Уравнения неразрывности (5-1), количества движения (5-2) и энергии (5-3) содержат четыре величины, характеризующие фронт ударной волны: рф, рф, І7ф и D. Задание любой из этих величин однозначно определит остальные с помощью этих трех уравнений. Таким образом, между газодинамическими параметрами на фронте существует однозначная связь, выражающаяся в виде алгебраи­ ческих соотношений.

§ 5-3. Основные соотношения между параметрами на фронте ударной волны

Произведем некоторые преобразования выведенных уравнений. Определяя £7ф из уравнения (5-1) и подставляя его значение

б уравнение (5-2), получим

 

 

 

1 Ръ~Рп

 

(5-4)

£>2 = Дг

 

 

Po _L__L:

 

 

Р о

рФ

 

 

и, = у 0

 

1

(5-5)

 

 

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

Выразим отношение

через рф и ро. Для этого запишем

уравнение (5-3) следующим образом:

гО (^ _ fir) p"s = Ui ( Pt ~ Т p"D ü t )

(5'30

Подставляя в уравнение

(5-3') значение р0О£/ф из

(5-2), будем

и меть

 

 

 

 

 

 

1

(Р*

Р±

P Q D

= и ф PФ + Ро

 

(5-3")

к -

1 VРф

P o

J

2

 

 

53

Введя в

(5-3") отношение

и*

из (5-1)

н сделав

преобразова-

ния, найдем

 

 

 

 

D

 

 

 

 

+

1 _ Р

ф

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ф

 

 

( К

 

рѵ _

к

 

 

 

Р

_

( к + 1 ) Р ф +

 

1)

 

 

/С—

1

 

 

 

 

Р ^ "

( К — 1 ) р ф 4 - \ К + 1 ) Ра ~ Р *

к + 1

 

 

 

 

 

 

 

 

 

 

 

рф н

Рм

' К

1

 

 

 

Уравнение

(5-6),

связывающее

 

р|(1

в

процессе

сжатия

в ударной волне, носит название уравнения

Погонно,

или

уравне­

ния ударной адиабаты.

 

 

 

 

 

 

 

 

 

 

 

 

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

волны D. Подставим

рфр„

из

(5-6)

в уравнение

(5-4). После

не­

сложных преобразований

 

 

 

 

 

 

 

 

 

 

 

 

 

 

о ’ - ' У р . Ч

 

 

* > + « £

 

 

 

(5-7)

и далее

 

 

 

 

 

 

 

 

 

°п = ,\ / / К +

 

 

 

 

 

 

 

 

 

 

 

 

 

1 .

 

 

 

 

2

 

 

(5-8)

 

 

 

 

2 р ,

 

 

 

 

 

 

 

 

 

 

легко определяются

І7ф

и

 

Гф:

 

 

 

 

 

 

 

 

 

 

, ,

 

Р ф — Po

 

 

 

 

 

(5-9)

 

 

 

 

 

 

Po D

'■

 

 

Р ф

 

 

 

 

 

 

 

Р ф Р о

'

 

Р ф

 

IC -}- 1

'

 

 

 

 

 

 

Ф

"

 

К

1

Р о

 

 

(5-10)

 

 

" Р ф Р о

"""

Р о

К +

1

Р ф

Ь

 

 

 

 

 

 

 

 

 

 

 

 

К ~ \

Р о

 

рф, D . и

С1ф

Уравнения

(5-6),

(5-8) и (5-9) дают

зависимость

1

от Рф в явной форіме іи поэтому удобны для расчетов.

Для ударных волн, распространяющихся в воздухе при стандартных условиях, имеем: ро=Ю321 кГ/м2, Т 0 = 288“ К, ро=0,125 кГ • сек2/м4, к=1,4, ао= 342 м/сек. Подставляя эти значе,- ния в выражения (5-6) —(5-10), получаем расчетные формулы (давления в кГ/см2) \

Рф

6/^ф

р Р(\

(5-60

Рч

Р ф

+

бР о 5

 

£> =

3421/1 +0,822 (рф - р п);

(5-80

t f * = 8 n 0 ^ ^ ;

(5-90

Гф “

°

Рф ■

(5-100

 

54

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

состояния газа, а также предполагалось постоянство к = — .

Поэтому полученные результаты должны быть справедливы при всех условиях, когда можно считать справедливым характеристи­ ческое уравнение и постоянство к. Это означает, что полученные соотношения должны достаточно хорошо выражать связь между параметрами во фронте ударной волны независимо от вида, мощ­ ности и формы источника ударной волны и формы поверхности фронта. Правильность этих соотношений подтвердилась при экспе­ риментальном изучении ударных волн при различных скачках уплотнения, при взрывах обычных взрывчатых веществ и- при взрыве ядерных зарядов для давлений во фронте до нескольких десятков кГ/см2. Для больших давлений во фронте необходимо учитывать изменение к вследствие заметного изменения тепло­ емкости, а также появления диссоциации газа.

§ 5-4. Дополнительные соотношения между параметрами на фронте ударной волны

Кроме полученных выше соотношений между параметрами на фронте ударной волны, удобных для определения всех параметров по Рф, можно получить еще ряд соотношений,- используемых при

решении различных задач,

в частности

при изучении течения

; в ударных трубах.

 

 

 

 

 

Найдем отношение Рф'рй. Для этого запишем уравнение (5-7)

в вице

 

 

 

 

 

 

 

D2 =

к + 1 „ Рп

 

Po

(5-70

Подставим

Po

2к

р„ Po

)

Р

к ---

ао2

и решим уравнение относительно

----Ф .

Тогда

Po

 

 

 

 

Ро

 

Рф _ к — 1 / 2к D2

 

(5-11)

 

~РІ~ к + \ W — 1 ' ^

 

 

 

 

Получим другое выражение для І7ф. Запишем выражение

(5-2)

следующим образом:

 

 

 

 

 

 

 

D U b = P * —

"

 

(5-20

Из уравнения (5-7)

 

Po

 

 

 

 

 

 

 

 

Рф— Ро

2

(D2-

at).

(5-7")

 

 

Po

к + 1

55

Уравнения (5-2') и (5-7") определяют выражение для £/ф:

 

 

 

 

2

 

 

 

2

 

 

 

(5-12)

 

 

U,

 

 

D - —

 

 

 

 

 

Ф

к - И І

 

 

D

 

 

 

 

Найдем другое выражение для — .

Запишем уравнение Гюго-

нио (5-6) в виде

 

 

 

 

 

 

рФ

 

 

 

 

 

Р ф

 

+

1

 

 

 

/

а

к

 

Рп

, «к — 1

к — 1

\ 1

 

:-И

к

 

Г

К + \ Рф

1

к + 1

к;Р±— 1+ К ~+ Х1 I _

 

 

 

 

 

 

Ро

к +

1

к — \ рп

 

 

 

 

 

 

к

- 1

1 +

 

 

 

4k

 

 

 

(5-6">

к

1

 

 

-f- 1

/?ф

 

 

Подставив в уравнение

(5-6")

 

 

 

 

Рф

(5-11), получим

значение —

 

 

 

к — 1

 

 

 

 

РкI

 

 

 

 

Р о

 

 

 

 

 

оіі

 

 

(5-13)

 

 

 

К +

1 1+

к — I

D2

 

 

 

 

 

 

 

 

Представим выражение D через

Uф, оф и о„. Для этого запишем

уравнение (5-7) так:

 

 

 

 

 

 

 

 

 

 

 

 

D2

 

Рф-Р»

,

 

k

Рф + ра

 

(5-7'">

 

 

 

2р„

'

 

ГА,

 

2

 

 

Из уравнений (5-2)

и (5-3")

 

 

 

 

 

 

 

 

 

Рф— Ро = DU,

 

 

 

 

 

 

 

 

 

 

 

 

Ф>

 

 

 

 

 

 

 

 

Рф+Ро

1

 

 

D

(Рф

Р^

 

 

 

2ру

 

к -

1

 

и Ф \ Рф

Р о

 

 

Подставляя полученные выражения в уравнение (5-7'"), будем

иметь

 

 

U,

 

 

 

 

 

 

 

 

 

 

 

 

 

1

 

 

2

2

 

 

 

 

 

 

+

I

 

fl-ф

CLQ

 

 

(5-14>

 

D = —

— 1

и *

 

 

 

 

 

2

 

 

 

 

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

связь между параметрами

С7ф, оф и а0. Для его получения

под­

ставим значение D (5-14)

в (5-12). После преобразований

 

и ф4 -]----- (о* + #о)

к {к — I)2 (Оф — Оо)2 = 0.

(5-15).

56

Отношение давлений

можно выразить через отношение плот-

p.

Р»

 

уравнением

Гюгонио. В ре-

ностей — , связанных

между собой

Ро

 

,г г,

'

Рф

получим

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

(5-6) относительно •—

FU1,

Ро

или

Ро

( к + 1)Рф — ( « — 1) Ро

(5-16)

( « + 1 ) Р о — ( « — 1 ) Р ф

рф

к + 1

Р ф

 

к — \

Р о

(5-16')

 

 

 

 

к — 1

 

Т,

на основании формул (5-10)

и (5-160

Отношение температур

Л)

 

'

 

может быть записано следующим образом:

 

 

к + 1 _£о_

 

2*

к

1 Рііі

 

То

к +

1— — .

(5-1б")'

 

__________________ _ Р ф _

 

 

К------- 1

р „

 

§ 5-5. Ударная адиабата

!

Выше было получено соотношение между плотностью и давле­ нием на фронте ударной волны (5-6), носящее название ударной адиабаты, или адиабаты Гюгонио. При рассмотрении обратимого адиабатного процесса с идеальным газом получают уравнение: адиабаты в виде

_ Р _

( 5 - т

Р о

 

Это уравнение называют просто уравнением адиабаты, или уравнением адиабаты Пуассона.

При сравнении выражений (5-6) и (5-17) можно заметить, что они дают различную связь между плотностью и давлением. Выяс­ ним особенности этого различия.

При сжатии по адиабате Пуассона, как следует из уравне-

ния (5-17), отношение — неограниченно растет при неограничен-

о

Р о

 

ном увеличении

— . Согласно уравнению (5-6) имеем

 

Р о

 

 

 

 

П т — =

(5-18)

 

 

р»

* _ 1

 

 

Рс

 

5 7

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

к + 1

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

может превзойти отношение ----- г

К— I

 

 

ударной волны. Следовательно, при ударном сжатии

 

Рф

К + 1

(5-19)

1 < — <

к — 1

Р о

 

Из курса технической термодинамики известно, что при пере­ ходе газа от состояния, определенного давлением ро и плот­ ностью ро, к состоянию с р и р изменение энтропии определяется выражением (1-9):

As = с ( In —- -f ln

-■

= cv \n P_

P o

(5-20)

V A .

P

Po

p

 

Обратимый процесс сжатия на основании (5-20) и (5-17) идет при As = 0, т. е. при постоянной энтропии.

Выясним, как изменяется энтропия при ударном сжатии. Ока­ зывается, что при ударном сжатии всегда

Р_

р

\ Л

(5-21)

Ро

 

)

Р о

 

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

деляя — из уравнения (5-6) и подставляя найденное значение

Ро

в (5-21), будем иметь

 

 

 

 

 

 

 

 

 

 

 

 

 

к — 1

 

>

 

 

 

(5-22)

 

 

 

 

к

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

к — 1

 

 

 

 

 

 

где х = -

 

 

 

 

 

 

 

 

 

 

 

Ро

 

 

 

 

 

 

 

 

 

 

Так как для ударного сжатия любой интенсивности

 

 

 

1 < * <

к

 

1

 

 

 

 

 

 

л — Г

 

 

 

 

то неравенство (5-22)

равносильно неравенству

 

 

 

 

 

1

 

 

 

 

к + 1

< 0 ,

 

 

>’= 1 - к — 1X

X

\X — к — 1

 

или

 

^ "I-

 

 

 

 

 

 

 

 

 

 

1

I

I ^

Д* I

к

X

к-И

г\

/сг оо\

 

У = 1 — -— гX

+ -— Г

X

<0.

(5-23)

 

 

1

К

 

1

 

 

 

 

 

Когда

х=1, г/= 0.

Поэтому для

неравенства

(5-23)

при усло­

вии (5-21)

необходимым и достаточным является

<

0.

58

Произведем дифференцирование по

ли

 

dy _

к + 1

1 -f кх

- 1

( « - 1)л'‘

(5-24)

dx

к — 1

к

 

 

 

 

При ,ѵ = 0 функция

dy

= 0. Определим знак производной

d2y

 

 

И х 2 :

 

d2y

 

= к (к -|- 1)

х к"2 (1 — х).

(5-25)

І х 2

Так какх>1, то

-^— -всегда меньше нуля. Это доказывает спра­

ведливость неравенств (5-25) и (5-21).

Вследствие выполнения неравенства (5-19) значение квадрат­ ной скобки в уравнении (5-20) всегда больше единицы и поэтому при ударном сжатии Д5>0, т. е. при ударном сжатии энтропия всегда возрастает.

Увеличение энтропии в результате ударного сжатия объяс­ няется необратимым характером изменения состояния газа. В ре­ зультате такого процесса часть кинетической энергии газа необра­ тимо переходит в тепло. При отсутствии теплообмена с внешней средой внутренняя энергия, а значит и температура, возрастает. Графики изменения плотности и температуры при изоэнтропичеокоім т ударном сжатии для газа к= 1,4 приведены на рис. 5-5.

При формальном распространении ударной адиабаты (5-6) на

случай разрежения ( — < 1 и — <1) на основании уравнения (5-20)

• Р о Ро 1

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

Рассмотрим уравнение ударной адиабаты для слабой ударной

волны. В уравнении

(5-6) положим

рп = р; рф = Р dp\ р0=р; рф =

= Р + dp:

(« +

1) +

dp) +

-

1) р

 

р + dp

(5-60

р

(к — 1) +

dp) +

+

1) р

 

или

 

 

 

 

 

 

 

dp _

 

dp______

 

(5-6")

 

р —

, к — 1 ,

 

 

 

 

 

 

 

'

кр 4----- dP

 

 

 

59

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