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

книги из ГПНТБ / Гофман М.Л. Аэродинамика гиперзвуковых скоростей и супераэродинамика

.pdf
Скачиваний:
42
Добавлен:
30.10.2023
Размер:
17 Mб
Скачать

центрации, т. е. содержание свободных электронов в одном кубиче­ ском сантиметре.

Заметное увеличение ионизации атмосферы наступает на высоте 90—100 км, т. е. сразу же за верхней границей мезосферы. Значи­ тельная концентрация ионов имеется и на больших высотах. Так, да­ же на высоте 800 км обнаружено 160 тыс. ионов в 1 см3 и нейтраль­ ных частиц до одного миллиона.

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

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

Раньше полагали, что атмосфера простирается примерно до вы­ соты 1000 км. Эта условная граница определялась исходя из того, что на этой высоте плотность воздуха становилась равной плотности рассеянных газов в космическом пространстве (плотности «межпла­ нетного» газа). Однако из наблюдений над торможением советских, и американских искусственных спутников 1 установлено, что земная атмосфера простирается примерно до 3000 км [40].

За последнее время изменились наши представления о межпла­ нетном газе. Если несколько лет тому назад считали, что плотность межпланетного газа ничтожна и концентрация частиц в нем не пре­ вышает нескольких единиц в 1 см3, то в последнее время возникла другая точка зрения, по которой межпланетный газ является значи­ тельно более плотной средой с концентрацией до 1000 частиц в 1 см3.

Исследования верхней атмосферы с помощью космических ракет и спутников показали [41, 42], что за пределами воздушной атмосфе­ ры землю окружает новая, радиационная атмосфера, состоящая из протонов и электронов. Радиационная атмосфера состоит из двух зон земной «короны»: внутренней, простирающейся от 1000 до 6000 км

(предполагают, что она заполнена протонами, которые

возникают

в верхних слоях атмосферы под влиянием космических

лучей) и

внешней, заполненной электронами. На высоте около 60 тыс. км внешняя зона короны переходит в межпланетный газ.

Состав атмосферы на больших высотах значительно отличается от состава атмосферы вблизи уровня моря. До высоты порядка 80 км

1 Существование первого советского спутника три месяца вместо предпола­ гавшегося срока в 1 год объясняется тем, что на высоте порядка 600 км плот­ ность воздуха оказалась в 20 раз больше расчетной.

30

молекулярный вес остается примерно постоянным. На высотах 90— 180 км основной причиной изменения молекулярного веса является диссоциация кислорода. На высотах более 180 км условиями, опре­ деляющими изменение молекулярного веса смеси атмосферных га­ зов, являются диффузионное разделение газов и диссоциация азота.

За последнее время методы исследования физических свойств ат­ мосферы больших высот значительно расширились благодаря на­ блюдениям с помощью ракет и спутников [43]. Так, анализ материа­ лов наблюдений над спутниками показывает, что плотность на высо­ те 270 км в 10 млрд, раз меньше, чем у земли, а при увеличении вы-

Рис. 1.14. Изменение плотности воздуха с высотой

соты еще на 100 км она уменьшается еще в 10—12 раз. Вообще ат­ мосферу на определенной высоте от поверхности земли можно харак­ теризовать двумя параметрами: ее плотностью и высотой так назы­ ваемой однородной атмосферы. Высота однородной атмосферы (иногда называемая геопотенциальной) прямо пропорциональна температуре воздуха и обратно пропорциональна его молекуляр­ ному весу.

Так, у поверхности земли, где средняя масса молекул составляет 4,8-10~23г, а 7’=273°К, высота однородной атмосферы Нодн =7,9 км. От этого параметра зависит скорость уменьшения величины плотно­ сти с высотой. Чем меньше величина однородной атмосферы, т. е. чем больше температура воздуха и чем меньше его молекулярный вес, тем быстрее происходит уменьшение плотности атмосферы с высотой.

На рис. 1.14 представлено изменение плотности с высотой на ос­ нове согласования экспериментальных данных (полученных при по­ мощи советских и американских искусственных спутников) с теоре­ тическими оценками изменения плотности. Кривая осреднена, так как имеется разнос экспериментальных данных, обусловленный ши-

31

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

Более трудной задачей является определение теплового режима атмосферы на больших высотах (выше 100 км). Прежде всего возни­ кает вопрос о возможности применения к этой области понятия тем­ пературы.

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

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

Кинетическая температура связана со средней кинетической энергией движения молекул известным соотношением

M v J

2

где М — масса молекулы,

v M— скорость молекул,

Тм— кинетическая температура молекулы.

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

Так, вследствие чрезвычайно малой плотности среды и большой длины свободного пробега молекул температура тела в основном оп­ ределяется лучистой энергией, которую оно поглощает от Солнца, Земли и нижних слоев атмосферы. Например, расположенные на по­ верхности спутника кремниевые батареи в зависимости от того, на­ ходятся ли они в земной тени или на солнце, меняют свою темпера­ туру всего лишь в пределах от 16 до 30°.

Касаясь методов определения температуры на больших высотах, следует отметить, что, несмотря на интенсивное развитие за послед­ нее время так называемых прямых 1 методов исследования атмосфе­

1 Прямыми они называются потому, что в этом случае определяется темпе­ ратура той точки пространства, где находится измерительная аппаратура.

32

ры с помощью ракет и спутников, до сих пор значительная часть ин­ формации относительно строения верхней атмосферы получена с по­ мощью косвенных методов [41].

К наиболее распространенным косвенным методам можно отне­ сти акустический метод, метод наблюдения метеоров, оптические ме­

тоды, радиометоды, наблю­

 

 

дения

атмосферных

прили­

 

 

вов

и др.

 

 

 

 

Прямые методы исследо­

 

 

вания

атмосферы

больших

 

 

высот чрезвычайно разнооб­

 

 

разны. К ним относятся: из­

 

 

мерение температуры на ра­

 

 

кетах с использованием тер­

 

 

мометра

сопротивления и

 

 

манометрических

измере­

 

 

ний;

 

измерение

скорости

 

 

звука и др., а также опреде­

 

 

ление температуры по наб­

 

 

людениям спутников, опре­

 

 

деление

плотности

атмосфе­

 

 

ры по наблюдаемому тормо­

 

 

жению

 

спутников

и другие

 

 

методы.

 

 

 

 

 

На рис. 1.15 дан график

 

 

изменения реальной кинети­

 

 

ческой температуры с высо-

р ис

i . 1 5 Осредненный график измене-

той.

На

графике

заметен

ния

реальной кинетической температуры

рост

температуры на

боль-

 

с высотой

ших

высотах. Так,

 

напри­

 

 

мер,. в ионосфере температура растет, достигая величины порядка 2000 °К 1 и более. Основную роль в нагревании верхней атмосферы на высотах более 100 км играют два источника: ультрафиолетовое излучение солнца и ролнечные корпускулярные потоки. Для раз­ ных областей атмосферы существенными являются разные источ­ ники. Так, в пределах ионосферы важную роль играет фотоиониза­ ция. За пределами ионосферы фотоионизация утрачивает свое зна­ чение и причиной нагревания становятся корпускулярные потоки и земное корпускулярное излучение.

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

1 Считается, что в ионосфере температура возрастает на 4 °К на 1 к м высо­

ты [44].

3 М. Л. Гофман

33

 

 

 

 

 

 

Таблица 2

 

Осредненная модель атмосферы больших высот

 

Н .

Р .

fj,

т° к

g -

М ,

к м

к г \ с м 2

к г - с е к - / м 4

м \ с е к 2

г / м о л ь

1

2

3

4

5

6

0

1,033

1,25- 10-1

288

6,81

28,97

1

9,16- 10-1

1,13- ю - 1

281,5

9,81

28,97

2

8,11 - 10-1

1,03- ю - 1

275

9,81

28,97

3

7,15.

10-1

9,2710~2

268,5

9,81

28,97

4

6,28- 10-1

8,35-10“ 2

262

9,81

28,97

5

5,51 10-1

7,51 - 10~2

255,5

9,81

28,97

6

4,81- Н Г 1

6,73-10-2

249

9,81

28,97

7

4,18- Ю '1

6,01 -10—2

242,5

9,81

28,97

8

3,63

-10-1

5,36-10“ 2

236

9,81

28,97

9

3,13-10—1

4,76.10-2

229,5

9,8Д

28,97

10

2,69

-10-1

4,21 - 10-2 *

223

9,81

28,97

11

2,31

-10-1

3,71 10-2

216,5

9,81

28,97

12

1,97- 10-1

3,1710~2

216,5

9,81

28,97 .

13

1,68-10-1

2,71 - 10-2

216,5

9,81

28,97

14

1,44

-10-1

2,3 L-10 ~2

216,5

9,81

28,97

15

1,23

-10-1

1,97.10~2

216,5

9,81

28,97

16

1,05

-10-1

1,69-10-2

216,5

9,81

28,97

17

8,95

-10-2

1,4410“ 2

216,5

9,81

28,97

18

7,64

-10~2

1,23-10-2

216,5

9,81

28,97

19

6,53

-10~2

1,05.10- 2

216,5

9,81

28,97

20

5,57- Ю-2

8,97-10"3

216,5

9,81

28,97

22

4,07

-10-2

6,54- 10-3

216,5

9,81

28,97

24

2,96

-10-2

4,77-10-3

216,5

9,81

28,97

26

2,16-10~2

3,4810~3

216,5

9,81

28,97

28

1,58-10"2

2,54- И)"3

216,5

9,81

28,97

30

1,15-10~2

1,85- 10-3

216,5

9,81

28,97

32

8,75

-10~3

1,39-10_3

224

9,80

28,97

34

6,65

-10~3

1,04-10_3

232

9,79

28,97

36

5,05

-10-3

7,76- Ю-4

240

9,78

28,97

38

3,84

-10~3

■5,80-10-4

247

9,77

28,97

40

2,92

-10-3

4,35- Ю-4

255

9,76

28,97

42

2,22 -ю -3

3,26-10-4

262

9,75

28,97

44

1,69-10-3

2,44-10-4

270

9,74

28,97

46

1,28-10“ 3

1,83-10-4

278

9,73

28,97

48

9,75-10-4

1,36- ю - 4

285

9,72

28,97

50

7,44

-10-4

1,02-10“ 4

293

9,71

28,97

52

5,64

-10-4

7,65-10~5

287

9,70

28,97

54

' 4,29

-10-4

5,72-10“ 5

281

9,69

28,97

56

3,26

-10“ 4

4,29-10“ 5

274

9,68

28,97

58

2,48

-10~4

3,20-10-5

268

9,67

28,97

34

 

 

 

 

 

 

 

Продолжение

1

2

3

 

4

5

6

60

1,8910~4

2,40-10-5

262

9,66

28,97

62

1,43-10~4

1,80-10-5

256

9,65

28,97

64

1,09-10-4

1,35-10-5

250

9,64

28,96

66

8,28-10-5

1,07-10-5

243

9.63

28,96

68

6,30-10~5

7,54-10

6

237

9,62

28,95

70

4,81-10-5

5,65-

10-8

231

9,61

28,95

72

3,6610~5

4,23-

10

8

225

9,60

28,94

74

2,7910“ s

3,16-

10-8

219

9,59

28,94

76

2Д 2-10—5

2,36-

ГО'"6

212

9,58

28,94

78

1,61 - 10“ 5

1,77-

10-е

206

9,57

28,93

80

1,23-10-5

1,33-

ю -8

200

9,57

28,93

82

9,3410- 6

9,94-

10

7

200

9,56

28,93

84

7,1010 6

7,44. 10-7

200

9,55

28,93

86

5,40-10-6

5,55-

10-7

200

9,54

28,92

88

4,10-10-6

4,15-

10-7

200

9,54

28,91

90

3,12-10~6

3,12-10-7

200

9,53

28,9

92

2,38- 10-6

2,33-10-7

200

9,53

28,64

94

1,81- 10-6

1,74- 10-7

200

9,52

28,38

96

1,37-10-6

1,31 - 10-7

203

9,52

28,12

98

1,04-10-®

9,78-

10-8

209

9,51

27,86

100

7,95- 10-7

7,28-10-8

215

9,51

27,6

110

2,29-10-7

1,92-10-8

250

9,48

26,4

120

1,10- ю - 7

6,71 • 10—9

285

9,45

25,0

130

4,57ГО"8

2,83-10-9

319

9,42

24,6

140

2,40-10-8

1,44- 10-9

354

9,39

24,4

150

1,32-10-8

6,53-10

10

389

9,36

24,2

160

7,95-10-9

3,98- Ю -10

424

9,33

24,0

170

4,79-10—9

2,29- Ю -10

459

9,30

23,8

180

3,16-10-9

1,44.10—10

493

9,27

23,6

190

2,0910~9

9,1210~п

528

9,24

23,0

200

1,51 - 10~9

6,42- 10-и

563

9,22

22,4

250

3,63-10—10

1,26

-10-11

737

9,08

20,1

300

1,1510 -ш

3,63-10"12

911

8,94

18,6

350

5,01 • 10-11

1,17

-10~12

1088 .

8,81

17,7

400

2,40- 10-п

4,90-10-13

1260

8,68

17,05

450

1,20-10-11

2,29-10-13

1434

8,55

16,6

500

6,91-10-12

1,20

-10 -13

1609

8,43

16,25

550

3,98-10-12

6,30-10—14

1783

8,31

15,95

600

2,51 • 10~12

3,63

-10 -14

1958

8,19

15,7

650

1,58-10 -12

1,99

-10 -14

2132

8,07

 

700

1,05-10-12

1,05

-10“ 14

2307

7,96

 

750

7,25-10-13

5,75

-10~!5

2431

7,88

 

800

5,01 ■10-13

3,16

-10“‘15

2655

7,74

 

3*

35

Количественные расчеты распределения плотности, давления и температуры атмосферы на больших высотах проводились многими учеными. Модели верхней атмосферы, предложенные разными иссле­ дователями [45—51], не совпадают вследствие различия исходных предположений. Поэтому целесообразно пользоваться «сглаженны­ ми» данными по имеющимся моделям.

В таблице 2 представлены осредненные нами физические свойст­ ва атмосферы больших высот, составленные по экспериментальным и расчетным данным различных моделей верхней атмосферы.

ГЛАВА II

НЕКОТОРОЕ МЕТОДЫ РЕШЕНИЯ ЗАДАЧИ ОБТЕКАНИЯ ПРИ ГИПЕРЗВУКОВЫХ СКОРОСТЯХ

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

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

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

Ниже (§2.1 и 2.2) будет сделан обзор некоторых наиболее рас- ’• пространрнных методов расчета обтекания тел при гиперзвуковых скоростях; подробное рассмотрение отдельных методов представле­

но в § 2.3—2.7.

§ 2.1. МЕТОДЫ РАСЧЕТА ОБТЕКАНИЯ ТОНКИХ ТЕЛ

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

37

нии потоков, обтекающих тонкие тела [55], т. е. тела, которые вносят в поток лишь небольшие возмущения. Теория малых возмущений при гиперзвуковых скоростях является характерной теорией такого класса.

Слово «малых» применимо потому, что скорости возмущений ма­ лы по сравнению со скоростью невозмущенного потока и давления малы по сравнению с дина­

 

Дкмерет

Г и п е р и у н

мическим давлением

невоз­

 

 

 

мущенного потока. Однако,

 

 

 

как уже указывалось в § 1.4,

 

 

 

скорости возмущений не ма­

 

 

 

лы по сравнению со ско­

 

 

 

ростью

звука

невозмущен­

 

 

 

ного потока, а возмущения

 

 

 

давлений не малы по срав­

 

 

 

нению со статическим дав­

 

 

 

лением

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

по­

 

 

 

тока.

Поэтому

возмущения

 

 

 

вообще нельзя считать ма­

 

 

 

лыми в том смысле, в каком

 

 

 

они применяются в лине­

 

 

 

аризированной

теории,

и

 

 

 

теория

малых

возмущений

 

 

 

существенно нелинейна.

 

 

 

 

Таким образом, линеа­

 

 

 

ризированное

.

уравнение

 

 

 

для

потенциала

скорости

Рис. 2.1. Коэффициенты су и сх плоской

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

гиперзву­

ковых скоростей,

тем самым

пластинки,

вычисленные по теории Акке­

становится. неприменимым

рета

и по нелинейной

теории

для таких скоростей

реше­

 

 

 

ние

Аккерета

для профиля

и метод источников и стоков для тел вращения. Действительно, можно показать, что для гиперзвукового потока коэффициент давления изменяется нелинейно с изменением толщины и угла атаки по сравнению с известными линейными законами для тон­ ких тел в плоском сверхзвуковом потоке. Для иллюстрации этого на рис. 2.1 дано сравнение коэффициентов су и сх плоской пластинки, вычисленных по теории Аккерета и по нелинейной теории; начиная с Mst;5 заметна нелинейность кривых коэффи­ циентов су и сх .

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

38

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