Файл: Шишков А.А. Газодинамика пороховых ракетных двигателей. Инженерные методы расчета.pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 09.07.2024

Просмотров: 200

Скачиваний: 1

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

 

 

 

 

 

 

Таблица 6

 

 

 

/г= 1,4

 

£= 1,25

£= 1,15

со°

М[

а°

Mi«

Р2 / Р 1

Mi«

Mi„

 

 

8

3

25,61

1,30

1,80

1,28

1,26

 

4

20,47

1,40

2,12

1,37

1,35

 

5

17,57

1,51

' 2,49

1,47

1,45

 

6

15,73

1,63

2,92

1,58

1,55

 

10

12,37

2,14

5,19

2,06

2,00

10

3

27,38

1,38

2,05

1,35

1,34

 

4

22,23

1,51

2,51

1,48

1,45

 

5

19,38

1,66

3,04

1,61

1,58

 

6

17,59

1,81

3,67

1,75

1,71

 

10

14,43

2,49

7,07

2,38

2,30

12

3

29,25

1,47

2,34

1,44

1,42

 

4

24,09

1,63

2,94

1,59

1,55

 

5

21,28

1,81

3,68

1,75

1,71

 

6

19,55

2,01

4,54

1,93

1,88

 

10

16,58

2,85

9,33

2,71

2,62

14

3

31,22

1,55

2,65

1,51

1,48

 

4

26,05

1,76

3,43

1,70

1,66

 

5

23,29

1,98

4,39

1,90

1,85

 

6

21,61

2,21

5,53

2,12

2,06

 

10

18,81

3,22

11,96

- 3,05

2,94

20

3

37,76

1,84

3,77

1,77

1,72

 

4

32,46

2,15

5,21

2,05

1,99

 

5

29,80

2,48

7,04

2,37

2,29

 

6

28,26

2,84

9,25

2,70

2,60

 

10

25,82

4,36

21,96

■4,10

3,96

30

3

52,01

2,36

6,36

2,22

2,14

 

4

45,22

2,84

9,24

2,67

2,57

 

5

42,34

3,37

13,07

3,17

3,04

 

6

40,79

3,92

17,76

3,68

3,53

 

10

38,52

6,23

45,08

5,82

5,57

27


После преобразования получаем (см., например, [88])

t g со = c t g - -а- - - - -

At?- - - -sin2- - - -а- - - - 1- - - -

- - -. - - - - -

(50)

Mj

— sin2а) +

1

 

Изменения параметров газа при переходе через косой скачок в зависимости от числа Мі набегающего потока и угла ю откло­ нения скорости в косом скачке приведены в табл. 6 для /г=1,4; 1,25 и 1,15.

Из формулы (50) видно, что при бесконечно малых углах от­ клонения потока (coä^O)

M 7 s i n 2 c c — 1;

а = a r c s i n —

,

1

At

 

т. е. косой скачок совпадает с линией Маха. При этом из тре­ угольников, показанных на рис. 6, с точностью до бесконечно ма­ лых первого порядка следует (с/щ<0)

tfco=

dibctea

dv,

,

 

 

 

—!—

— ------ L ctg arcsin—— ) = • — i- M ,

 

 

Vo

Vi

At

«1

M?

 

 

 

 

 

или

dvi

 

du>

 

 

 

(51)

 

v\

 

Ѵ м р T

 

 

 

 

Используя уравнения сохранения массы и полного импульса

при

бесконечно малом

повороте:

d (p + gy2) = 0 и

dp = Qvdv,

а также соотношение дѵ2 = рѵ2/ (RT) = к р М .2, получим

 

 

 

 

dp __

k№ d v__

АМЗДо

 

(52)

 

 

 

р

V f /

' At2 — 1

 

 

1.6.6. Расширение сверхзвукового потока

Расширение сверхзвукового потока (в противоположность сжатию) происходит изэнтропически. Наиболее простым случаем

 

расширения сверхзвукового по­

 

тока является

центрированная

 

волна разрежения Прандтля-

 

Майера,

возникающая

при

 

внешнем

обтекании

плоским

 

потоком

 

выпуклого

угла

 

(рис. 7).

 

 

 

 

 

 

Изэнтропическое изменение

 

скорости

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

пото­

 

ка при его отклонении на бес­

 

конечно малый угол определя­

 

ется уравнением (51)

(dvt> 0)

Рис. 7. Внешнее обтекание сверхзву­

dv

 

 

 

 

 

 

 

f f М2— 1

 

ковым потоком тупого угла

V

 

~

 

 

 

__

rfw

 

 

28



0

25

30

35

40

45

50

55

60

65

70

75

80

85

90

95

100

105

ПО

115

120

125

130

135

140

270

k = 1,25

 

 

 

ü)°

Л

 

r/ rKV

л(7)

0

1,000

1,000

0,5549

0

1 ° 17'

1,081

 

1,098

0,4989

23°32'

2°10'

1,114

1,146

0,4762

30°

3°13'

1,152

1,205

0,4504

34°54'

4°35'

1.194

1,277

0,4222

38°52'

6°13'

1,239

1,365

0,3926

45°24'

8°07'

1,287

1,471

0,3617

51°

10°12'

1,338

1,597

0,3299

55°50'

12°27'

1,392

1,706

0,2975

60°20'

15°02'

1,445

 

1,933

0,2672

66°24'

17°43'

1,501

 

2,158

0,2368

72°06'

20°34'

1,558

2,418

0,2078

75°42'

23°38'

1,614

2,744

0,1812

80°52'

26°47'

1,672

3,146

0,1557

85°48'

30°01'

1,731

 

3,636

0,1321

О CD СО

33°25'

1,789

 

4,257

0,1111

95°05'

36°54'

1,847

 

5,025

0,0923

101°

40°28'

1,905

 

5,992

0,0757

106°48'

44°08'

1,962

 

7,236

0,0614

112°21'

47°51'

2,018

 

8,865

0,0492

1І7°54'

51°45'

2,074 10,988

0,0388

123° 18'

55°31'

2,128

 

13,757

0,0303

126°

59°31'

2,182

 

17,495

0,0232

133°54'

63°22'

2,235

22,482

0,0175

135°10'

67°25'

2,286

29,369

0,0130

140°20'

О О СО

3,000

 

СО

0

220°27'

 

 

Таблица 7

k== 1,4

 

 

0)°

X

Г / Г к р

л(л)

 

 

 

'

0

1,000

1,000

0,528

1

1,067

1,087

0,479

2

1,107

1,147

0,450

3

1,142

1,205

0,424

4

1,172

1,262

0,402

6

1,227

1,372

0,364

8

1,277

1,498

0,330

10

1,323

1,626

0,299

12

1,367

1,772

0,271

15

1,428

2,005

0,234

18

1,486

2,291

0,201

20

1,523

2,500

0,181

23

1,576

2,858

0,154

26

1,628

3,319

0,130

28

1,660

3,647

0,116

32

1,722

4,436

0,0920

36

1,782

5,521

0,0717

40

1,838

6,919

0,0552

44

1,891

8,710

0,0419

48

1,943

11,20

0,0310

52

1,990

14,72

0,0229

54

2,014

16,90

0,0194

60

2,080

26,30

0,0115

61

2,090

28,32

0,0105

65

2,130

39,63

0,0071

130°27'

2,449

оо

0

29


Введя между линиями Маха угол ф= и + (я/2) — arcsin 1/М и приведенную скорость Х=ѵ/а1(р, получим после преобразования

<р= arcsln ] / (Х2 “ 1 )• (5 3 )

Зная зависимость ср (Л,) или А,(ср), можно с помощью газоди­ намических функций определить изменение всех параметров при изменении ф. Из уравнения неразрывности следует уравнение линии тока при внешнем обтекании тупого угла (Ä,i=l):

 

rKP= //?(X) = rsin (arcsin

q (к)=СЯМ.

($4)

или r/rKP = M./q(X).

 

 

Значения

ф, со = ф+ агсзіп(1/М) — я/2, X, г/г1{]) и л(Х)

приве­

дены в табл.

7.

 

 

Предельные значения угла поворота потока со и угла между линиями Маха ф соответствуют расширению в пустоту:

При таком режиме истечения граничная характеристика (ли­ ния Маха) совпадает с направлением линии тока отклонивше­ гося течения, так как

arcsin ( — ) = 0 .

V Mo !

Следовательно, С0т ах= фтах — (я/2) = (я/2) (/Vox — 1) ■

Формулы (53) и (54) применимы и для значений Л-і> 1, но при этом начальный угол отсчета определяется по формуле*

*1=V arcsin]/ ^ Іг 1)’

и линия отсчета располагается не перпендикулярно направлениюневозмущенного потока, а под углом фі + arcsin (1/Mi) к нему.

Приведенные примеры иллюстрируют эффективность аппа­ рата газодинамических функций как при теоретическом анали­ зе, так и при конкретном расчете газовых течений. Разложения газодинамических функций в ряд по малому параметру [(k l)/.(k+ 1)]Х2 (при Я <1) с точностью до членов порядка X2' включительно дают приближенные решения, соответствующие те­ чению несжимаемой жидкости; и в этом также заключается одно из достоинств представления параметров газа с помощью газо­ динамических функций.

30

1.7. ТУРБУЛЕНТНЫЙ ОТРЫВ ПОТОКА

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

Это объясняется тем, что движение жидкости, при котором два слоя двигались бы друг относительно друга, «скользя» один по другому, является абсолютно неустойчивым [52]. Другими словами, поверхность раздела между этими двумя слоями жид­ кости не может быть поверхностью «тангенциального разрыва», на которой скорость жидкости (направленная по касательной к поверхности) испытывала бы скачок. Дело заключается в том, что при малейшем возмущении поверхности раздела давление над гребнями уменьшается, а внутри впадин — увеличивается, что приводит к дальнейшему нарастанию возмущений и, в ко­ нечном счете, к разрушению поверхности раздела и к перемеши­ ванию газовых частиц струи и застойной зоны.

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

сеивание

механической энергии газового

потока

происходит,

в основном, лишь в области вихревого движения.

 

Таким

образом,

неустойчивость тангенциального разрыва,

смешение

потоков,

образование вихревой

зоны,

диссипация

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

Турбулентная область ограничена с какой-либо стороны ча­ стью поверхности обтекаемого газом тела. Линию пересечения границы между вихревым и потенциальным течением с поверх­ ностью тела называют л и н и е й о т р ы в а , так как именно с этой линии начинается «проникновение» ротора скорости в ос­ новной поток, а в силу закона сохранения циркуляции такое «проникновение» может произойти только путем непосредствен­ ного перемещения движущейся вблизи поверхности жидкости в глубь основного потока.

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

31