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

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

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

р

Г

 

 

 

 

 

Р

^

конуса. Сравнивая эти графики,.

2=100

2=1

 

 

 

 

 

2=5

можно сделать выводы

о влия­

 

 

 

 

 

 

 

 

 

нии параметров

Моо, рк,

а на

 

 

 

 

 

 

 

г

 

характер

зависимости

Рг

от 0.

71

 

 

 

 

г=юо

 

 

Например, уменьшение

разно­

 

 

 

 

 

 

 

 

сти рк — а (при а ^

рк) приво­

 

 

 

 

х-

 

 

 

 

дит к

появлению экстремумов

 

 

[

X. * _ /г

 

Ч - ----------4 .4 0 -\о .7 0

функции

Рг (0)

при

меньших

 

-1.725

 

‘-х- г

— — - Ч \

значениях 0 и сближению

друг

 

 

 

 

 

т

У

\

 

с другом

дальше от

полуплос­

 

 

 

 

 

 

 

 

кости 0 =

я, чем в случае боль­

76 Н

 

 

 

\

 

—1

шей разности рк — а.

 

отме­

------ ------ -------—/~ Ч\---- — -А\

Вполне понятно,

что

 

 

 

 

 

 

 

 

 

ченный характер функции/1! (0)

 

 

 

 

 

 

 

 

 

и связанные с ним явления за­

 

т о

 

 

 

/

 

 

- 0.69

трудняют расчет течений около

 

 

 

 

 

4.30

тупых конусов,

начиная

с не­

 

 

 

 

 

 

 

7

 

которого значения 2.

 

 

 

 

 

 

 

 

1 /

 

А

 

На

фиг.

18.6

приведены

 

 

 

 

 

 

Л

 

графики функции Р (0) для ту­

69

 

 

 

 

г~

 

 

 

пого конуса, Моо = 7, рк =

30°,

 

 

 

 

 

 

5 у

 

а — 10°. В этом случае при ма­

 

 

 

 

 

 

 

лых ъ расстояние до ударной

 

 

 

 

 

 

 

\ V

 

 

 

 

 

 

 

 

___ I 4.20

 

волны возрастает при увеличе-

 

1475

 

 

 

 

 

0.68

НИИ 0.

В полуплоскости 0 = 0

 

 

 

 

 

60

 

120

 

оно минимально,

а в полупло­

 

 

 

 

 

Фпг.

18.6

 

скости

0 =

я — максимально.

 

 

 

 

 

 

 

 

 

С увеличением значения ъ ха­

 

 

 

 

 

 

 

 

 

рактер функции Р (0)

становит­

ся другим: большие расстояния до ударной волны имеются с наветренной стороны конуса, а меньшие — с подветренной. При этом максимум функции Р (0) располага­ ется не в плоскости симметрии течения. На фиг. 18.6 расстояние до ударной волны максимально при 0 ^ 60°, а при 0 = я имеется местный минимум. Расстояние до ударной волны около острого конуса для а << Рк определено в работе [15]. Там же бы­ ло обнаружено, что это расстояние у острого конуса может быть максимально не в плоскости симметрии течения. Таким образом, при обтекании тупого конуса максимальное расстояние до волны может быть и не в плоскости симметрии течения.

На фиг. 18.6 также приведены графики функции Рг (0) для тупого (пунктир) и острого (крестики) конусов. Как видно, в этом случае график Рг (0) при больших 2 имеет вид, подобный виду графика Р (0), а значения наклонов ударной волны

К ОСИ 2 у

острого и тупого конусов при

ъ = 100 отличаются очень мало.

На

фиг. 18.7 приведены графики

функции Р (0) для случая обтекания ту­ пого конуса, Мьо = 6, рк= 25°, а = 10°. На оси ординат показано начало отсче­ та Р для каждой кривой. Масштаб всех графиков функции Р (0) одинаков. Здесь представлен интересный эффект — по­ явление на поверхности ударной вол­ ны, образованной тупым конусом, впа­ дины и перемещение ее при увеличе­ нии значения г к плоскости симметрии

течения (к полуплоскости 0 = я).

Фиг. 18.7

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

2. Изолинии в плоскостях г = сопз!. Рассмотрим изолинии газодинамических функций в плоскостях 2 = сопз!, которые для конуса имеют более сложный вид по сравнению с изолиниями около параболоидов, рассмотренными в § 17.

На фиг.

18.8 приведены

изобары р = сопз! в плоскости ъ = 1.0, Моо = 2, рк =

= 15°, а =

10°. Это сечение

расположено недалеко от плоскости, проходящей

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

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

поверхности слабого разрыва не было и все поле течения рассчитывалось по единому алгоритму.

В области влияния конуса, как видно из фиг. 18.8, давление по нормали к его поверхности вблизи конуса изменяется мало, а в области влияния сферы — сущест­ венно. С наветренной стороны конуса с увеличением значения ъ излом изобар проис­ ходит под все более острым углом, затем они принимают вид характерных «языкооб­

разных» кривых

и, наконец,

разрываются в полуплоскости 0 = 0. Про дальней­

шем увеличении значения ъ с

наветренной стороны около тупого конуса постепенно

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

картина изобар, соответствующая течению около острого конуса.

На фиг. 18.9 приведены изобары для того же случая течения, что и на фиг. 18.8, но в плоскости 2 = 5 . Пунктиром нанесены изобары для острого конуса, рассчитан­ ные при тех же значениях Моо, Рк» оь. Для сравнения поверхности обоих конусов совмещены (в острый конус «вложен» тупой). С наветренной стороны конусов изоба­ ры имеют примерно один и тот же вид, хотя количественно еще отличаются друг от друга. С подветренной стороны изобары отличаются весьма существенно. Об­ ратим внимание также на следы ударных волн тупого и острого конусов. Они отли­ чаются значительно и с наветренной и с подветренной сторон. Поэтому нельзя считать, что при 2 = 5 10 течение около тупого конуса не отличается от течения около соответствующего острого конуса [138].

При дальнейшем увеличении значения ъ различия в картинах изобар для тупо­ го и острого конусов будут все меньше и меньше. На фиг. 18.10 приведены изобары в плоскости 2 = 50 для того же случая течения, что и на фиг. 18.8 и 18.9. Ударные вол­ ны и изобары у обоих конусов совпадают с наветренной стороны (при 0 <; х/2 я). Однако с подветренной стороны имеются еще количественные отличия, хотя качествен­ но изобары для тупого и острого конусов имеют один и тот же характер.

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

Из сравнения фиг. 18.8 и 17.15, на которых представлены изобары около тупого конуса и параболоида, следует, что в этих двух случаях характер изобар качественно не отличается при малых 2 (при 5^ 0 .3 ). При больших ъ отличия весьма существенны.

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

На фиг. 18.11 приведены изохоры в плоскости 2 = 5, М» = 2, рк = 15°, а = = 10°. Видно, что изохоры «подходят» к поверхности тупого конуса при 0<; у4я почти по касательной (направление изохор при 5 = 0 блиэко к направлению касательной к поверхности конуса в этой же точке плоскости 2 = 5 ) . Таким образом, в этой об­ ласти течения градиенты энтропии по нормали к поверхности уже значительны. На фиг. 18.11 пунктиром нанесены изохоры для соответствующего острого конуса при тех'же значениях М®, рк, Видно, что изохоры около поверхности острого конуса имеют совсем другой характер. В частности, около поверхности острого конуса при 7вл с 0 с У2я плотность почти не изменяется по направлению нормали. При 0 > У 2я изохоры имеют противоположное направление «подхода» к поверхности остро­ го конуса по сравнению с тупым.

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

поверхности слабого разрыва не было и все поле течения рассчитывалось по единому алгоритму.

В области влияния конуса, как видно из фиг. 18.8, давление по нормали к его поверхности вблизи конуса изменяется мало, а в области влияния сферы — сущест­ венно. С наветренной стороны конуса с увеличением значения г излом изобар проис­ ходит под все более острым углом, затем они принимают вид характерных «языкооб­

разных» кривых

и, наконец,

разрываются в полуплоскости 0 = 0. Про дальней­

шем увеличении значения 2 с

наветренной стороны около тупого конуса постепенно

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

картина изобар, соответствующая течению около острого конуса.

На фиг. 18.9 приведены изобары для того же случая течения, что и на фиг. 18.8, но в плоскости 2 = 5. Пунктиром нанесены изобары для острого конуса, рассчитан­ ные при тех же значениях Моо, Рк> Для сравнения поверхности обоих конусов совмещены (в острый конус «вложен» тупой). С наветренной стороны конусов изоба­

ры имеют

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

и тот же вид, хотя количественно еще отличаются друг

от друга.

С подветренной

стороны изобары отличаются весьма существенно. Об­

ратим внимание также на следы ударных волн тупого и острого конусов. Они отли­ чаются значительно и с наветренной и с подветренной сторон. Поэтому нельзя считать, что при 2 = 5 ч- 10 течение около тупого конуса не отличается от течения около соответствующего острого конуса [138].

При дальнейшем увеличении значения ъ различия в картинах изобар для тупо­ го и острого конусов будут все меньше и меньше. На фиг. 18.10 приведены изобары в плоскости 2 = 50 для того же случая течения, что и на фиг. 18.8 и 18.9. Ударные вол­ ны и изобары у обоих конусов совпадают с наветренной стороны (при 0 < х/2я). Од­ нако с подветренной стороны имеются еще количественные отличия, хотя качествен­ но ивобары для тупого и острого конусов имеют один и тот же характер.

Фиг. 18.12 Фиг. 18.13

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

Из сравнения фиг. 18.8 и 17.15, на которых представлены иэобары около тупого конуса и параболоида, следует, что в этих двух случаях характер изобар качественно не отличается при малых 2 (при? ^ 0 .3 ). При больших 2 отличия весьма существенны.

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

На фиг. 18.11 приведены изохоры в плоскости 2 = 5, Моо = 2, рк = 15°, а = = 10°. Видно, что изохоры «подходят» к поверхности тупого конуса при 0 < г/гл почти по касательной (направление изохор при | = 0 близко к направлению касательной к поверхности конуса в этой же точке плоскости 2 = 5 ) . Таким образом, в этой об­ ласти течения градиенты энтропии по нормали к поверхности уже значительны. На фиг. 18.11 пунктиром нанесены изохоры для соответствующего острого конуса при тех'же значениях М », рк» Видно, что изохоры около поверхности острого конуса имеют совсем другой характер. В частности, около поверхности острого конуса при 1/вл с 0 С 1/гл: плотность почти не изменяется по направлению нормали. При 0 > 7 2л;изохоры имеют противоположное направление «подхода» к поверхности остро­ го конуса по сравнению с тупым.

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

На фиг. 18.12 приведены изо­ линии числа М в плоскости 2 = 1 ,

Моо = 2, рк =

15°, ос =

10°. Здесь,

так же как

и

на картине изобар

при малых

2,

хорошо

видна об­

ласть влияния конуса на картину течения. При увеличении ъ вид изолиний М изменяется. На фиг. 18.13 приведены изолинии числа М для того же случая обте­ кания, что и на фиг. 18.12, но в плоскости 2 = 5. Пунктиром на­ несены эти же изолинии около острого конуса. Хорошо видно резкое изменение характера изо­ линий вблизи поверхности кону­ са. Отметим, что максимум М достигается внутри области тече­ ния в полуплоскости 0 = я.

Изолинии окружной компо­ ненты вектора скорости т около тупого конуса, |3К= 15° при М» = = 2, а = 10° качественно не от­ личаются от ш-изолиний около па­ раболоида. На фиг. 18.14 приве­ дены и;-изолинии около тупого и

острого

конусов

в

плоскости

2 = 5 ,

Моо = 2 ,

Рк=15°,

а =

= 10°. Они несколько

отличаются

друг от друга,

особенно

около

ударной

волны

с подветренной

стороны течения. В этом же при­ мере, но при 2 = 30 изолинии около острого и тупого конусов совпадают полностью (с точностью до 2—3 единиц третьего знака ко­ ординат), хотя с подветренной сто­

роны ударные волны

отличаются

еще существенно.

 

 

Покажем, как влияет на поле

течения изменение

числа

Моо.

Для этого

рассмотрим изолинии

различных функций около того же

самого тупого конуса и при том же угле атаки, но при Мто =

7. В плоскости 2 = 1

и близких к ней плоскостях качественная картина изолиний

практически остается

той же самой, что и при Мсо = 2.

7, рк =

15°, а =

10°.

На фиг. 18.15 приведены изобары в плоскости ъ = 5, М«> =

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

чения продолжают оставаться весьма существенными, в то время как с наветренной стороны изобары полностью совпадают. Это иллюстрирует фиг. 18.16, где показаны изобары в плоскости 2 = 30, Мю = 7, [Зк= 15°, а = 10°. Поле давления с подвет­ ренной стороны при Моо = 7 качественно отличается от поля давления около парабо­

лоида (фиг.

17.16)

и от поля давления около этого же тупого конуса при М» = 2

(фиг. 18.9).

18.17

приведены изохоры в плоскости 2 = 5 ,

= 7, рк = 15°, а =

На фиг.

= 10°. Пунктиром нанесены изохоры для острого конуса. Отметим различие в ха­ рактере поля плотности около острого и тупого конусов с подветренной стороны, которое сохраняется вплоть до

2 = 5 .

На фиг. 18.18 для этого же случая приведены изолинии числа М в плоскости 2 = 5 около тупого и острого конусов. Здесь опять следует обратить внимание на совершенно раз­ личный характер линий по­ стоянных значений числа М у тупого и острого конусов. Эти отличия весьма существенны (качественно и количественно) во всем поле течения и с подвет­ ренной и с наветренной сторон.

Ф и г . 18.16

Фиг. 18.18

Рассмотрим изолинии в поле течения около более «толстого» конуса рк = 30е при Моо = 7 и а = 10°. В этом случае асимптотические значения устанавливаются при значительно меньших я, чем в предыдущем. Все отличия поля течения около ту­ пого конуса от соответствующего острого локализуются в более тонком, но более мощ­ ном энтропийном слое (в энтропийном слое имеются еще большие, чем в предыдущем случае, градиенты функций по нормали к поверхности конуса).

На фиг. 18.19 приведены изобары в плоскости ъ = 10 для тупого и острого (пунк­ тирные линии) конусов. При 0 < 2/зя эти изобары совпадают (с точностью до единицы третьего знака координат). Как видно из фиг. 18.19, поля давления около тупого и ост­ рого конусов отличаются лишь количественно в небольшой области с подветренной стороны.

На фиг. 18.20 для этого же случая течения приведены изохоры в плоскости ъ = = 10 для тупого и острого (пунктирные линии) конусов. Эти изохоры с наветренной стороны течения в области, прилегающей к ударной волне, не отличаются. В области энтропийного слоя и с подветренной стороны течения они отличаются значительно.

Ж

в

На' фиг.

1&21

приведены

изолинии ю

 

 

плоскости

2 =

10, Моо = 7,

Рк =

30°,

 

 

а = 10°. Сплошной кривой показаны изо­

 

 

линии для тупого- а

пунктирной — для

 

 

острого

конусов. При

6 < : 7 г

Л Изолинии

50

 

в этих ,двух течениях совпадают. С

подвет-

 

ренной* стороны’течения имейтОя су1цествен-

 

 

ные качественные иколичественйые отличия,

 

 

которые х!орошо виднй н а 'фйг. 18:21: Отме­

 

 

тим, что в плоскости ъ =

30 в этом же случае

2.5

 

обтекания окружные кЬ^пЪненты

скорости

 

у ‘острого и тупого конуЬЬв практически пол­

 

 

ностью совпадают.

приведены

изоэнтропы

 

 

 

На

фиг.

18.-22

 

 

в плоскости ъ =

 

50, Мое = 7 , рк =

30°. Во

о

 

всем поле потока, кроме'гнеболыпой области

 

с подветренной

стороны

течения, характер

 

 

изоэнтроп для

тупого конуса соответствует

 

 

течению около острого конуса

(см.,

напри­

 

 

мер, [15]). Отличия имеются только

с под­

г‘5

 

ветренной стороны течения и в энтропийном

 

слое. Однако эйропийный слой

в плоскости

 

 

2 =

50 очень тонкий и на рисунке практиче­

 

 

ски совпадает

с

 

поверхностью

тела. Есте­

^

 

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

на поверх-

 

ности тупого конуса

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

 

 

энтропии на поверхности острого

конуса. На

 

 

фиг. 18.22 заметно, что расстояние от поверх­

 

 

ности тупого конуса до ударной волны в по­

ц

 

луплоскости 0 =

 

л меньше, чем

в полупло-

 

скости 0 = 0. Такой вид ударной

волны из­

 

 

вестен и для острых конусов с большим

уг­

 

 

лом полураствора

[15].

 

 

 

 

М

 

 

 

На фиг.

18.23

приведены изолинии

 

Фиг. 18.19

в плоскости г =

30

для

течений

около

ту­

 

 

пого и острого (пунктирные линии)

конусов, М« = 6, (Зк = 25°, а = 10°. В этом при­

мере также видны существенные различия в течениях

около тупого и острого кону­

сов

с подветренной

стороны.

 

 

 

 

 

11 для случая обтекания ту­

 

На фиг. 18.24 приведены изобары в плоскости 2 =

пого конуса,

Моо

=

6,

= 25°, а =

17°30'

С подветренной стороны изобары об­

разуют

седловую

точку,

которая при увеличении 2 перемещается к ударной волне.

Обратим внимание на то, что эта особенность в поведении изобар проявляется в третьем знаке давления. В наветренной области течения при 2 = 1 1 поле давления име­ ет качественно такой же характер, как и у острого конуса.

Выше были приведены примеры течений около тупых конусов при а < рк. Рас­ смотрим примеры течений около тупых конусов при а > рк.

При углах атаки, близких к углу полураствора конуса, течение с подветренной стороны приобретает еще более сложный характер. На фиг. 18.25 приведены изоба­ ры в плоскости 2 = 21, Моо = 10, (5К= 10°, а = 10°. Как и на фиг. 18.24, изобары образуют здесь седловую точку. Она появилась у поверхности конуса при 2 ^ 5 . С увеличением 2 седловая точка перемещается по направлению к ударной волне.

Иа фиг. 18.26 приведены изобары в плоскости 2 = 1.5 для случая течения около тупого конуса, Моо = 6, рк = 13°, а = 15°. По характеру изобар хорошо просле­ живается граница области влияния сферы на поле течения. С подветренной сто­ роны конуса еще не образовалась седловая точка. Она образуется при больших

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