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

книги / Течения газа около тупых тел. Метод расчета и анализ течений А. Н. Любимов, В. В. Русанов

.pdf
Скачиваний:
22
Добавлен:
12.11.2023
Размер:
14.64 Mб
Скачать

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

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

При проведении экспериментальных исследований с большими сверхзвуковыми скоростями трудно обеспечить моделирование течений, близких к реальным. По­ этому экспериментальным путем в таких случаях нелегко получить даже распреде­ ление давления вдоль поверхности тел и форму ударных, волн. Например, при Моо = 20 -ь 25, так же как и при Моо = 2 -ь 6, не представляет больших принци­ пиальных затруднений процесс измерения давления на поверхности тел и фотогра­ фирование ударных волн при проведении экспериментов в аэродинамических трубах. Однако пересчет полученных экспериментальных данных на «натуру» вызывает серьезные затруднения.

В последние годы некоторое развитие получили оптические методы визуализа­ ции потока [186, 187]. Они позволяют получить общую картину поля плотности, что очень важно для выяснения всей картины течения. Детальное исследование поля плотности оптическими методами требует сложной аппаратуры и трудоемкой мето­ дики обработки результатов измерений. Заметим к тому же, что такие исследования затруднены тем, что надежные данные о деталях поля плотности оптическими мето­ дами могут быть получены при значительных градиентах плотности.

Исследования осесимметричных течений около тупых тел проводятся и путем отстрела моделей (см., например, [188]). Однако в этих экспериментальных данных, как правило, содержится очень мало информации о поле течения и точность ее мала.

2. Для исследования осесимметричных течений газа около тупых тел широко используются аналитические методы. Большое развитие получили аналитические методы, основанные на различных предположениях: о малой толщине ударного слоя, о постоянстве плотности в ударном слое и т. д. Подробный перечень и анализ этих методов содержится в работе [19]. Они иногда весьма полезны для проведения раз­ личных оценок. Например, аналитическая теория Буземана, точная для предель­ ного случая течения М» = <х>, к = 1, когда ударная волна совпадает с поверхно­ стью тела, дает полезные оценки давления и при конечных, очень больших, значе­ ниях числа Моо.

Значительное развитие получили аналитические методы, которые дают прибли­ женное решение основных уравнений газовой динамики. Они обычно основаны на представлении решения в виде рядов. В результате основная система уравнений га­ зовой динамики сводится к более простой системе обыкновенных или алгебраиче­ ских уравнений. В некоторой области течения с помощью методов, основанных на представлении искомых функций в виде рядов, можно найти значения функций, близкие к точным (см., например, [34, 38]). Основной недостаток этих методов со­ стоит в том, что область сходимости полученных рядов, как правило, неизвестна. Поэтому пока не получено каким-либо другим методом точное решение задачи, неизвестна и область течения, в которой могут быть использованы полученныеаналитические решения.

Аналитические методы применяются для изучения многих отдельных вопросов механики течений около тупых тел: для определения вида линий тока, формы удар­ ных волн и звуковых линий, коэффициентов отражения возмущений от ударных волн и т. д. (см., например, работы [189—191]). Аналитические методы в таких случаях позволяют делать обобщения на целый класс задач. В этом заключается их основное преимущество перед экспериментальными и численными методами. С помощью последних обобщения можно делать только после проведения большой серии экспе­ риментов или расчетов. Отметим, что на эффективное аналитическое решение отдель-

1 Подробнее о связи экспериментальных и численных исследований сказано в § 1.

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

3. Основные сведения об осесимметричных смешанных течениях газа около ту­ пых тел к настоящему времени получены с помощью численных методов. В § 2 наз­ ваны работы, посвященные численному исследованию обтекания тупых тел идеаль­ ным газом. Отметим еще раз основные, на наш взгляд, работы, в которых не только предложены численные методы расчета обтекания тупых тел, но и получены кон­ кретные и ценные сведения о течениях [23, 46, 48, 52, 57, 74—77]. В то же время опуб­ ликован целый ряд работ, в которых или почти не приводятся результаты расчетов, или полученные результаты имеют невысокую точность. Поэтому, к сожалению, иногда существуют значительные отклонения при оценке характеристик одного и того же случая обтекания с помощью численных результатов, полученных различ­ ными авторами. Приведем один пример.

В табл. 13.1 даны расстояния от критической точки до ударной волны при обтекании воздухом = 1.4) сферы и эллипсоида. Они определены с помощью раз­ личных методов. Отметим, что расстояние от критической точки до ударной волны находится более точно, чем все остальные параметры течения, с помощью любого из этих численных методов. Тем не менее, как видно из табл. 13.1, при определении расстояния отхода ударной волны от тела с помощью различных численных методов совпадают только первые два десятичных знака.

 

 

 

Т а б л и ц а 13.1

Источник

Сфера, Мод = 4

Сфера, Мдо = 6

Эллипсоид,

Ъ/а= 0.5, Моо = 6

[58]

0.175

0.149

0.159

[56]

0.175

0.147

[21]

0.17402028

0.1577

[23]

0.174

0.148

Настоящая [ работа

0.1753

0.1488

0.1605

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

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

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

§14. Осесимметричные и плоские течения около носков тел

Из физических соображений и анализа общей структуры смешанного течения (см. § 8) следует, что течение около носовой части тела определяется только тем участ­ ком его образующей, который примыкает к дозвуковой и, возможно, небольшой части сверхзвуковой области. Предполагая непрерывную зависимость течения от формы тела, можно ожидать, что для тел, образующие которых близки на определяющем участке, структура течения будет аналогична. Если образующая тела имеет конеч-

1Г.З

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

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

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

Если 2, г — координаты в плоскости течения, то уравнение конического сечения

имеет вид

г* =.2р*1+да»,

 

(14.1)

 

 

 

где р — радиус кривизны в вершине и д — параметр,

определяющий тип кривой.

При д =

О получаем параболу, при д < 0 — эллипс с полуосями — р/д и р / ]/"—

при д =

1&2е > 0 — гиперболоид с углом между асимптотами 2е. В частности, при

д = — 1 получаем сферу радиуса р.

 

 

Предположим, что уравнение тела вблизи вершины может быть представлено

в виде 2 = а2г2 + я4г4 +

О (г6). Тогда, положив р =

д = —

мы полу­

чим коническое сечение,

аппроксимирующее образующую вблизи точки ъ = г = О

с точностью О (гв) = О (23), в то время как соприкасающаяся сфера

аппроксими­

рует ее, вообще говоря, с точностью О (г4).

 

 

В этом и следующем параграфах рассматриваются осесимметричные и плоские течения около тел с образующей вида (14.1), а также около тел с образующей вида г = 2а (а = 0.25 и а = 0.125).

Заметим, что для равновесных течений газа изменение параметра р, т. е. радиуса кривизны, эквивалентно изменению масштаба и не влияет на само течение. По-

1.0

2.9

3?0

Фиг. 14.1

этому задача

обтекания для тел*

г

с образующей

вида (14.1)

имеет

 

фактически

 

двухпараметрическое

 

семейство

решений,

зависящее

 

только от Моо и

 

Параметр р мо­

 

жет

быть

произволен,

хотя, ко­

03

нечно,

для

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

нужно приводить

к одному значе­

 

нию р.

 

 

 

 

 

 

 

 

 

 

На фиг. 14.1 приведены гра­

 

фики

изменения

 

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

 

кривизны образующих сферы, па­

 

раболоида

и гиперболоидов

(е =

 

= 10°,

е =

25°

и

е =

35°)

при

о.2

р =

1. По

оси абсцисс

отложено

произведение длины дуги на

зна­

 

чение кривизны в вершине тел, а

 

по оси ординат — отношение кри­

 

визны

образующей

к

кривизне

 

в вершине. На фиг.

14.2 приведе­

 

ны графики

изменения

кривизны

 

образующих

тел

с

уравнениями

 

г =

г0*25 и г —г0125. По оси абсцисс

 

отложена

длина

дуги

обра­

 

зующей от начала

системы

коор­

 

динат.

 

 

 

 

 

 

 

 

 

 

Изменение углов наклона ка­

 

сательных к поверхностям тел по­

 

казано на фиг. 14.3. По

оси

орди­

 

нат отложено значение угла на­

0

клона касательной к образующей,

а по

оси абсцисс — длина

дуги

 

образующей от начала координат.

 

Анализ

поля

течения прове­

 

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

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

1. Изменение газодинамических функций в зависимости от значений ^ и со. Вначале рассмотрим изменение осевой и и радиальной Vкомпонент вектора скорости, давления р и плотности р в поле течения около сферы, а затем и около других тел.

ний

График и (|) при со ^ 45° близок к линейной зависимости. Для больших значе­

со (со

45°) и (!•) является существенно нелинейной функцией. При малых

Моо

и (I) монотонно убывает по !•, а при средних (Моо^5) и больших Моо имеет ми­

нимум. Чем больше Моо, тем ближе этот минимум сдвигается

к значению | = 0.

 

На фиг. 14.4 приведены графики функции и (5) при Моо =

20 и со = 0, 45, 90°

для совершенного газа (сплошные линии), а также с учетом равновесных физико­ химических превращений в газе. Условия невозмущенного потока в этих случаях приблизительно соответствуют высоте Я = 20 км (пунктир) и Я = 50 км (штрих-

пунктир). Учет физико-химических превращений в газе приводит к существенному изменению величины осевой компоненты вектора скорости, особенно на поверхности

сферы.

Характер зависимости и от 5 при этом качественно остается таким же, как

и для

совершенного газа.

малых

М«>

(Мто ^

2)

монотонно

Радиальная компонента V (|) при со = 45° и

убывает, а при Моо= 4, 6 имеет внутренний минимум. Для М» >

6

V (!•) уже мо­

 

нотонно

 

возрастающая

функция,

 

причем

 

д2у/д%2 >

0.

При

со =

90°

 

V (|)

также монотонно возрастает по

 

но

д2г>/д^2<с0.

 

Для

малых

М»

 

при

со

=

90°

функция

V (^)

имеет

 

внутренний максимум. На фиг. 14.5

 

приведены

 

графики

функции

V (^)

 

при со =

45,

90°

для

Моо =

4.

Из

 

этого графика,

в

частности,

видно,

 

что

величина

V— наибольшая

при

 

значениях со, близких к со

=

45°. Од­

 

нако изменения Vпо !*при со =45° не

 

такие большие, как при со

=

90°.

 

 

 

Давление вдоль критической ли­

 

нии тока (со = 0 ) монотонно убывает

 

по ^ для всех чисел М». Для малых

 

Моо

при

со

=

45° р (|) имеет

 

макси­

 

мум, а для Моо>3

монотонно

воз­

 

растает

 

по

к

причем

зависимость

 

р (^) близка

линейной. При

со =

 

=90° и малых Моо (Моо<12) д2р)д^<^

 

<С 0, а при больших

Моо

д2р/д%2^>*0.

 

Д Л Я

СреДНИХ

ЗНаЧеНИЙ

Моо (Моо =

 

= 4

 

6) зависимость

р (^)

близка

 

к прямой.

 

 

 

 

 

 

 

 

 

 

 

 

На фиг. 14,6 приведены графики

 

функции р (^)

для

Моо = 2 0 . Пунк­

 

тиром и штрих-пунктиром нане­

 

сены значения давления

с

учетом

0.2

О.Н

0.6

0.8

I

Фиг. 14.6

Фиг. 14.7

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

 

Из фиг.

14.6 видно,

что

учет фи­

 

зико-химических процессов более суще­

 

ственно влияет на величину

давления

 

около поверхности ударной волны. От­

 

метим, что в дозвуковой

области

тече­

 

ния изменение р от 5 близко

к

линей­

 

ному

или

 

квадратичному

 

закону.

 

В сверхзвуковой области течения функ­

 

ция

р (|)

имеет

более сложный закон

 

изменения.

 

 

 

 

 

 

 

 

 

 

Плотность вдоль критической ли­

 

нии тока

(со =

0) монотонно возрастает

 

по

достигая

в

критической

точке

 

максимального

значения,

 

причем

 

2р/д|2< 0 .

Для малых М» (М» < 2 )

 

при со =

45°

р (5)

имеет

внутренний

 

максимум.

При

со = 90°

для

малых

 

Моо д2р/д|2с

 

 

а Для

 

 

 

больше

5 = 1 близки при разных со,

0.

они, как

известно,

все

 

Дл

а для Моо = <х>

равны

при любом значении со. В этих же случаях на поверхности сферы плотность

изме­

няется весьма сильно. На фиг.

14.7 представлены графики р (5) для Мо©= ©о. От­

метим также, что в этом случае плотность на критической линии тока при всех 5 <С 1

превосходит р = 6.

 

 

 

 

 

 

 

 

 

 

 

 

Функции и (со) при всех значениях 5 и Моо являются возрастающими. Графики

этих функций имеют перегибы,

а график и (со) для 5 =

0 пересекает при больших

значениях со графики и (со) для

ё >■ 0. Следовательно, при малых со

значения осе­

вой компоненты при 5 = 0 меньше, чем при 5 =

1, а при больших со — наоборот.

На фиг. 14.8 приведены графики функций и (со) для различных 5 при Моо =

 

20.

Графики г>(со) для всех 5и Моо имеют внутренние максимумы, причем чем больше 5,

тем дальше вправо по со смещены максимумы г; (со). На фиг. 14.9 приведены графики

функций V (со) для различных 5 при Моо = 20.

 

 

 

 

 

 

 

 

 

 

 

Давление и плотность уменьшаются с увеличением со при всех 5 и Моо. Графики

р (со) и р (со) имеют перегибы. Для различных значений

5 они пересекаются. Пере­

сечение графиков происходит при тем меньших значениях со, чем больше число Моо. Таким образом, для малых 5 давление и плотность на поверхности сферы больше, чем на ударной волне, а при больших со, наоборот, меньше.

На фиг. 14.10 приведены графики р (со), а на фиг. 14.11 — графики р (со) при Моо = 20. Сплошной линией нанесены значения давления и плотности для совер­ шенного газа, а пунктирной и штрих-пунктирной — значения плотности для иде­ ального газа с учетом физико-химических превращений при условиях невозмущен­ ного потока, примерно соответствующих высотам Я = 20 км и Я = 50 км.

Учет физико-хими*:еских превращений весьма сильно влияет на значение плот­ ности, особенно при малых со. Это влияние гораздо существеннее сказывается на плотности газа, чем на давлении (см. фиг. 14.6) и компонентах вектора скорости (см. фиг. 14.4.). Кроме того, при учете физико-химических реакций качественно* изменяется характер зависимости р от со. У совершенного газа при Моо < оо плот­ ность монотонно уменьшается при увеличении со, причем величина плотности не­ значительно отличается от 6. У идеального газа с учетом физико-химических прев­ ращений значение плотности не монотонно уменьшается при увеличении значения со, а величина плотности весьма существенно отличается от 6 (фиг. 14.11).

Характер изменения основных газодинамических функций в поле течения около

Фиг. 14.11

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

На фиг. 14.12 и 14.13 представлены графики V (!) и г; (со) для течения около кругового цилиндра при Моо = 4 и Мсо= 20 соответственно. Характер других функ­ ций от ^ и со для кругового цилиндра качественно меньше отличается от соответ­ ствующих функций для сферы.

Рассмотрим изменение и1 V, р, р около параболоидов и параболических ци­ линдров.

Компонента скорости и (!) при всех со и М« — возрастающая функция, и толь­ ко при Моо <; 2 и больших со м(!) имеет внутренний максимум. В отличие от пара­ болоидов, у параболических цилиндров при больших со функция н(!) может иметь минимумы и перегибы.

На фиг. 14.14 представлены графики н(!) для параболического цилиндра, р = 0.25, Моо = 20. Они, в частности, показывают сложный и различный характер зависимости и (!) при со = 45° и со = 90°.

Радиальная компонента скорости течения около параболоидов при малых Моо является убывающей функцией координаты ! для всех со, а при Моо >►4 — возра­ стающей, причем даг;/д!2 с 0. Для параболического цилиндра эта функция имеет качественно несколько другой вид.

На фиг. 14.15 показаны функции V (!) для параболического цилиндра, р = 0.25 при Моо = 6 (пунктир) и параболоида (сплошные линии).

Функция и (со) в поле течения около параболоида является неубывающей. Гра­ фики и (со) имеют перегибы. При различных ! они, как правило, не пересекаются. Отличие графиков функции и (со) для параболического цилиндра состоит в том, что при различных ! они пересекаются, подобно графикам и (со)*для сферы*(см. фиг. 14.8).

О 20

40

ВО 80 ш

Фиг. 14.13

о.2

ом

о.б

о .з

е,

 

Фиг. 14.14

 

 

 

Функции г; (со) в поле течения около параболоида для всех ^ имеют максимумы,

причем при Моо<12 на

поверхности

параболоидов (^ = 0) величины V больше,

чем на поверхности ударной волны, а при Моо > 6, наоборот, меньше. При Мм = 4 (и близких к ним значениях) графики функций V (со) для различных 5 пересекаются.

На фиг. 14.16 приведены графики функций г; (со) для параболоида, Моо = 2, р = 0.25. Для параболического цилиндра функции V (со) имеют качественно такой же характер, как и для параболоида. Однако величина V (при тех же числах М» и в со­ ответствующих точках), как правило, меньше из-за того, что процесс расширения га­ за около параболического цилиндра менее интенсивен, чем около параболоида.

Характер изменения давления и плотности около параболоидов и параболиче­ ских цилиндров в зависимости от значения со подобен характеру изменения функций р (со) и р (со) около сферы (см. фиг. 14.10 и 14.11).

Рассмотрим изменение компонент век­

 

 

 

 

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

0.8

 

 

 

ле течения

около гиперболоидов и гипер­

 

 

 

болических

цилиндров (е = 10,

25, 35°).

 

 

 

 

Функция и (!) есть возрастающая функ­

 

 

0.5

 

ция !. При со =

0 она близка к линейной,

0.6

 

 

 

 

 

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

 

 

Го

 

! не отличается от характера этой зависи­

 

 

^

мости около сферы и параболического ци­

 

 

 

 

линдра (см. фиг. 14.4 и 14.14). При больших

ОМ

 

 

 

значениях со отличия становятся более су­

 

 

 

 

щественными (качественно и количествен­

 

 

 

 

но), причем

графики функции

и (!) для

0.2

 

 

 

гиперболоида и гиперболического цилинд­

 

 

 

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

 

 

 

 

Это хорошо видно из фиг. 14.17, где пред­

 

 

 

 

ставлены графики функций и (!)

для ги­

20

40

60

80 о>

перболоида и гиперболического цилиндра,

р = 0.25, е =

25°, Мое = 4.

 

 

Фиг. 14.16

 

 

Давление

и плотность при малых со

 

 

 

 

обычно являются убывающими функциями

 

 

 

 

а при

средних (со = 30 -г- 40°) и больших (со = 50 ч- 70°) — возрастающими.

При со =

90° (и близких со) функции р (!) и р (!) могут иметь внутренние макси­

мумы. Графики функций р (!) для различных значений со могут пересекаться.

На фиг. 14.18 приведены графики функции р (!) для различных со в поле тече­ ния около гиперболоида, р == 0.25, е = 35°, М» = 4. Хорошо видно, что с увеличе­ нием значения со плотность уменьшается и на поверхности гиперболоида и на по­ верхности ударной волны, а в середине потока она сначала уменьшается, а затем

 

начинает возрастать. Этот интересный эф­

 

фект

объясняется тем,

 

что

при

е =

 

35°,

 

со =

90°

уже

начался

заметный

про­

 

цесс перестройки потока от течения около

 

затупления к асимптотическому

течению.

 

Асимптотические значения

плотности

на

 

фиг. 14.18 нанесены пунктиром.

 

 

 

 

Зависимость осевой компоненты векто­

 

ра скорости от со для гиперболического ци­

 

линдра при различных Моо,

е

и

!

каче­

 

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

на

фиг.

 

14.8 графикам функции

и (со)

для сферы.

 

То же имеет место и

для

гиперболоида,

 

однако графики функции и (со) для различ­

 

ных ! у гиперболоида,

 

как

правило,

не

 

пересекаются, в отличие от графиков и (со)

 

для сферы и гиперболического

цилиндра.

 

Изменения

функции г; (со)

для гипер­

 

болоида и

параболоида

(см. фиг. 14.16)

 

подобны. Следует заметить,

что при со зна­

 

чительно больших, чем

со =

90°,

зависи­

 

мости г; от со для гиперболоида и парабо­

 

лоида будут существенно различны вслед­

 

ствие

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

 

этих

двух

течений.

Графики

функции

 

г; (со) в поле течения около гиперболоида и

Фиг. 14.17

гиперболического цилиндра,

построенные

Соседние файлы в папке книги