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

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

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

Рис. 2.14. Сравнение распределения давления на те­ лах вращения, вычисленного по методу касательных конусов и методу характеристик

С качон

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

§ 2.7. МЕТОДЫ ИСПОЛЬЗОВАНИЯ СООТНОШЕНИЙ НА СКАЧКЕ УПЛОТНЕНИЯ И В ВОЛНЕ ВОЗМУЩЕНИЙ. МЕТОД «СКАЧОК — РАЗРЕЖЕНИЕ»

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

Такие зависимости для сверхзвукового потока даны известными решениями Донова [83] и Буземана, рассматривавших течения рас­ ширения и сжатия с двух сторон профиля в предположении, что по­

токи не оказывают влияния друг на друга.

поверхности

профиля

Если обозначить через v угол наклона

к потоку, выражения для давления представятся в виде:1

Рраеш =

CjV +

+ c 8v3 + • ■•.

(2-29)

Раж —

+

Q v2 + (Сз — D) Vs

+ ...

(2.30)

Равенство (2.29) "Получено разложением коэффициента.давления в ряд из условия соотношения в волне возмущений, а выражение (2.30) — из условий на скачке.

Коэффициенты разложения определяются в виде:

 

c i

УМ* -

1

 

 

 

 

 

 

r _

k W + (М! -

2)s

 

3 =

 

2(М3 — I)2

 

к

/М-1 ЛЛ8 , 2k- — 7/г — 5 ЛД6,

 

—-— М8 +

------- я------- М6 +

(М2 - 1)Л - 1

 

 

 

1(2.31)

 

+ — "з+ П W - 2М3 +

D

(\ * т г 1М“

 

М4+ ( к - 3)М 2 + 2

(М2— 1)' г—- 1

 

 

 

1 Появление коэффициента D, отличающего коэффициент давления за скач­ ком от коэффициента давления в волне возмущений, объясняется отражением возмущений от поверхности скачка уплотнения,

61

Если применить условие гиперзвукового потока, то

 

 

 

 

 

 

 

 

 

 

 

k

 

/ М 2-

1«М ;

 

(М2 — 2)2 - М4;

(М2 — 1)*- 1 ~М 7

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

C l~ M

 

'

 

 

 

 

 

 

 

 

 

 

k +

\

 

 

 

 

 

k + \

АД ,

2k27k — 5 АД ,

, . ГАД

 

(2.32)

 

 

— 1—

M + -----------------м - 1+ О [М~3],

 

^

(£ +

1)2

 

 

5£3-

16А -

13

АД . ,

- ...

31

( С 3— £>) ^

-—

 

М + ---------—---------М-1 +

О [М '31

 

16

 

 

 

 

12

 

 

 

 

 

Для &= 1,4 формулы (2.29)

и (2.30)

примут вид:

 

Р р а ш = \

V +

1,2v- +

0,40Mv3 - f 0[M _1] Vs

(2.29а)

Реме =

^

 

+ 1,2 v3 +

0.36MV® +

О [М-Ч V3.

(2.30а)

Полученные

выражения

могут

быть

представлены в виде

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

К M»v.

 

 

Q

 

имеем:

 

 

 

 

 

 

 

Для любого k = —

 

 

 

 

 

 

 

 

 

 

Cv

 

 

 

 

 

 

 

 

 

I

pjm . =

 

+ - Ш

-

+

 

/Г4-0[А Г 2] ,

(2.296)

 

 

= Т

 

+^ Т ^ + ^ Т б ^ ^ +

0[/С2]’

(2,30б)

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

Если предположить, что первое соотношение учитывает оба фак­ тора, т. е,-волны расширения и скачки уплотнения, то совершаемая при этом ошибка будет составлять не более 10% от третьего члена равенства (2.306) в его интервале сходимости при А=1,4. Фактиче­

62

ски обобщение первого соотношения на волны сжатия означает пред­ положение об изэнтропичиости процесса сжатия, что оправдано при М >3.

Сравнение расчетов для коэффициента давления по первому со­ отношению (2.296) во всей области течения с обычной теорией для скачков и волн возмущений дает хорошее совпадение (рис. 2.16).

Рис. 2.16. Сравнение расчетов распределения давления по двум первым членам равенства 2.296

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

С помощью данного приближения были получены аналитические формулы для определения аэродинамических коэффициентов профилей при гиперзвуковых скоростях (см. § 4.2). Однако при больших значениях числа М даже при малых углах v совпаде­

ние приближенных значений р с точными ухудшается; причина состоит в том, что радиус сходимости рядов для коэффициента

давления пропорционален величине 1

Если участок контура профиля вблизи передней кромки прямоли­ неен, распределение давления можно рассчитать точно, используя

63

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

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

'Скачан волны возмущений. В ’уплотнения этом случае взаимодейст­

вие сводится к тому, что волна возмущений ис­

 

 

кривляет скачок, но ска­

Возмущения

Отраженные

чок не влияет

на течение

от тела

в волне возмущений.

 

возмущения

 

 

 

В случае, если интен­

 

 

сивность головного скач­

 

 

ка достаточно

велика

и

 

 

поверхность профиля

по­

 

 

падает в область влияния

 

 

скачка, необходимо учи­

Рис. 2.17. Взаимодействие возмущений с телом

тывать

изменения

давле­

ния, возникающие вслед­

 

 

ствие

отражения

от

по­

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

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

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

На рис. 2.18 дана форма скачка и распределение давления по 10% чечевицеобразному профилю при М = оо , определенное этими двумя методами. Из рисунка видно, что результаты расчетов хорошо совпадают.

64

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

1 дв_

дх » -/ М 2 — 1 ду

В качестве примера при­ ведены результаты расчетов распределения давления по методу «скачок — разреже­ ние» и сравнение их с мето­ дом характеристик для осе­ симметричных тел (рис. 2.19).

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

Рис. 2.18. Сравнение расчетов по методу характеристик и методу „скачок — разре­ жение* для 1 0 % чечевицеобразного про­

филя

метр подобия К^> 1 (например, для тела вращения К =

Следует иметь в виду, что этот метод дает удовлетворитель­

ную точность, пока скачок остается

присоединенным

и k > 1,2 .

На рис. 2.20 показана область применимости

метода

„скачок —

разрежение* при &= 1,4. Поясним этот график подробнее.

Мерой интенсивности отраженных волн является отношение

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

угла наклона вдоль линий возмущения, от-

ходящих

от

поверхности,

д'*

 

д'>

этого

-з— к скорости изменения -з--

угла вдоль

 

 

{/£]

 

OC<i

потока,

отраженных волн. Здесь v — угол поворота

сх— характеристика

1-го семейства,

ct — характеристика

2-го се­

мейства

(рис. 2.17).

При

чистом

течении

Прандтля — Мейера

— ■ — 0. Исследованием потока за косыми скачками уплотнения

5 М. Л. Гофман

65

 

 

 

 

 

 

дч

 

установлено, что при всех числах

М > 1

отношение

дс1

значи-

 

 

 

 

 

 

дч

 

тельно меньше

единицы,

если

 

 

 

дс2

макси-

угол у не приближается к

 

2,0

 

 

 

 

 

 

'

i.s

 

 

 

 

 

 

/L

1,0

 

 

 

 

 

 

 

0,5

 

 

 

 

/

 

 

оL

 

 

 

 

 

 

0.5

1,0

1.5

2.0

 

 

 

о

 

 

 

 

Конус-цилиндр

 

 

 

 

 

 

 

скачок-разреж ение

 

 

 

 

 

 

скачок-разрежение

 

 

Р_

 

 

с поправкой

 

 

 

 

метод характеристик

 

 

А»

 

 

 

 

 

 

 

ОL ___________ —------------------ -

 

 

 

О

0,5

1.0

1,5

2,0

 

 

 

 

Цилиндр с ожибалом

 

 

 

Рис. 2.19. Сравнение расчета распределения дав­ ления по методу характеристик и методу „ска­ чок— разрежение" для осесимметричных^ тел при М = 3

мально возможному углу отклонения потока. Показано, что даже

в предельном случае при М - * о о ,

k =

1,4 и у < 40° это

отноше­

ние никогда не бывает больше

чем 0,06. Однако при

/г->1

оно

быстро растет по величине и метод

скачков — разрежений

ста­

новится неточным.

Метод «скачок — разрежение» был применен для случая обтека­ ния крыла конечного размаха с присоединенными скачками уплотне­ ния [84]. В качестве примера было рассмотрено обтекание тонкого плоского треугольного крыла со сверхзвуковыми передними кромка­ ми при малых положительных углах атаки.

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

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

66

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

Рис. 2.20. Область применимости метода .скачок—раз­ режение”

Обзор методов и техники эксперимента при гиперзвуковых скоро­ стях представлен, например, в работах [85, 86]. Описание отдельных гиперзвуковых установок и методов экспериментирования на них да­ но в [87—92]. Результаты экспериментов по определению аэродина­ мических характеристик несущих поверхностей и тел вращения при больших числах М изложены в делом ряде работ (например, [93—98]).

5*

Г Л А В А III

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

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

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

§ 3.1. РАСПРЕДЕЛЕНИЕ ДАВЛЕНИЯ НА НЕСУЩИХ ПОВЕРХНОСТЯХ

Результаты измерения давления на пластинках и профилях [94] указывают на значительное отклонение экспериментальных данных от результатов, полученных теоретически без учета вязкости. Кроме обычно имеющих место в сверхзвуковом потоке расхождений опыт­ ных и теоретических данных в области задней кромки, при гиперзву­ ковых скоростях вследствие взаимодействия пограничного слоя со скачками уплотнения заметно значительное возрастание давления в области передней кромки и в области за угловыми точками поверх­ ности крыла (рис. 3.1). Подобное расхождение имеет место вообще сразу за максимальной толщиной профиля даже при отсутствии уг­ ловых точек, как, например, у линзообразного профиля ( рис. 3.2).

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

68

принять, что поверхность изменяет свою форму на величину, рав­ ную толщине вытеснения пограничного слоя.

На рис. 3.3 приведены теоретическое распределение давления на плоской пластинке при наличии пограничного слоя и результаты экс­ перимента при числе М= 6,9. График показывает, что получено хо­ рошее согласование результатов, хотя здесь пренебрегалось взаимо­

действием

скачка

уплотне­

ния с пограничным слоем.

Значительное

увеличение

давления в области у пе­

редней кромки

происходит

главным

образом

из-за

очень быстрого возрастания

толщины пограничного слоя

при

гиперзвуковых

скоро­

стях.

Это

иллюстрируется

графиком

влияния

числа М

на изменение толщины вы­

теснения

ламинарного по­

граничного слоя

на плоской

пластинке (рис. 3.4).

 

Переход от малых сверх­

звуковых скоростей к гипер­

звуковым

скоростям

ведет

к возрастанию толщины по­

граничного слоя для данно­

го числа Рейнольдса.

Утол­ Рис. 3.1. Типичное распределение давле­

щение настолько значитель­

ния по прямому крылу при М -

6,9

но, что оно существенно ме­

на носовой части тела (рис.

3.5).

няет распределение давления

На этом графике дано опытное распределение относительного ста­ тического давления в зависимости от расстояния от передней кром­ ки на пластинке и клине. Видно, что вблизи передней кромки отно­ сительное давление увеличивается в 4—7 раз.

Для определения давления на пластинке вблизи передней кромки можно использовать следующее приближенное выражение [99, 100]:

- £ - = 0,555

+

0,489 -j- . . . ,

(3.1)

Р.

У R e,

 

 

 

 

R e t

V „ x

 

 

 

 

 

где V» — ск о р о сть

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

п о то ка,

 

V.. — ки н ем ати ч еск ая в я зк о с ть ,

 

 

х расстоян и е вдоль пластины .

69

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