Научная статья на тему 'Трансзвуковое течение в окрестности точки излома профиля со свободной линией тока'

Трансзвуковое течение в окрестности точки излома профиля со свободной линией тока Текст научной статьи по специальности «Математика»

CC BY
84
24
i Надоели баннеры? Вы всегда можете отключить рекламу.
Журнал
Ученые записки ЦАГИ
ВАК
Область наук

Аннотация научной статьи по математике, автор научной работы — Диесперов В. И.

В работе получено асимптотическое разложение решения системы уравнений Эйлера, описывающее течения в окрестности точки излома профиля со свободной линией тока, когда в ней достигается скорость звука. В качестве главного члена этого разложения берется автомодельное решение уравнения Кармана с показателем авто-модельности п = 6/5.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Текст научной работы на тему «Трансзвуковое течение в окрестности точки излома профиля со свободной линией тока»

УЧЕНЫЕ ЗАПИСКИ Ц А Г И Том XIII 1982

М 2

УДК 533.6.011.35:629.7.025.73

ТРАНСЗВУКОВОЕ ТЕЧЕНИЕ В ОКРЕСТНОСТИ ТОЧКИ ИЗЛОМА ПРОФИЛЯ СО СВОБОДНОЙ ЛИНИЕЙ ТОКА

В. Н. Даесперов

В работе получено асимптотическое разложение решения системы уравнений Эйлера, описывающее течения в окрестности точки излома профиля со свободной линией тока, когда в ней достигается скорость звука. В качестве главного члена этого разложения берется автомодельное решение уравнения Кармана с показателем автомодельности « = 6/5.

Рассматривается точка излома на поверхности тела, образованная пересечением двух гладких кривых АО и ОВ (рис. 1). Предполагается, что скорость набегающего вдоль АО потока идеального газа дозвуковая и в точке О достигается скорость звука. Известно [1], что в некоторой окрестности звуковой точки все характерные особенности течения описываются уравнением Кармана. Оно допускает класс автомодельных решений с показателем автомодельности п [1]. В случае прямолинейной образующей АО решение уравнения Кармана, удовлетворяющее на ней условию непротекания, можно продолжить различными способами, приемлемыми с физической точки зрения, в область, лежащую за точкой О. Требование о том, чтобы оно переходило в окрестности положительной оси х в решение, описывающее течение Прандтля — Майера, приводит к спектру значений п = == (6/ — 1)/(6/ — 2) (/ = 1, 2, 3, . . .) [2, 3] и значению п — 2. Главное значение спектра п = 5/4 впервые было найдено в работе [2]. Если п — 2, то возмущения из сверхзвуковой части потока достигают угловой точки О вдоль предельной (^-характеристики [4]. В настоящей работе ставится требование, чтобы решение уравнения Кармана в области за угловой точкой описывало в ее окрестности течение со свободной линией тока ОВ (см. рис. 1), на которой скорость постоянна и* равна звуковой. Это требование приводит к

Рис. 1

спектру решений при п — 6//(6/ — 1)(/=1, 2, 3, . . .). Главное решение реализуется при п = 6/5. С использованием его строится решение системы уравнений Эйлера для задачи обтекания точки излома профиля со сходящей с нее свободной линией тока.

Набор „идеальных" решений интересен также с точки зрения теории вязких течений, так как общая проблема единственности решения уравнений Навье — Стокса для сжимаемого газа до сих пор не решена. Описанные выше „идеальные" решения (кроме случая п — 2) имеют благоприятный градиент давления, стремящийся к бесконечности при приближении к угловой точке. Выбор необходимого из ряда возможных невязких решений с бесконечными локальными градиентами давления проводится с помощью разработанной в последнее время асимптотической теории отрывных вязких течений [5 — 7]. Отрыв несжимаемой жидкости от угловой точки контура тела в рамках этой теории был исследован в [8].

1. Рассмотрим сначала обтекание точки излома профиля, когда образующая АО прямолинейна. Введем декартову систему координат (х, у) с началом в точке О, ось х которой совпадает с АО, и систему координат Мизеса (х, Ф). Далее через цх и <7у обозначим проекции вектора скорости на оси х и у соответственно, р-—плотность. В качестве характерных значений всех параметров потока берутся их критические значения. Термодинамические переменные связаны уравнением состояния идеального газа. Набегающий поток считается безвихревым. Ниже все параметры течения и уравнения, их связывающие, предполагаются безразмерными. Они берутся в виде

*>Мх--±-гп — Г).

дх дф ‘ дф ’

д$у дду да

си

здесь у-—отношение удельных теплоемкостей.

Обозначим через *ох и компоненты возмущенных скоростей ъу — ду. Окрестность угловой точки О, в которой

1 ух I ^ Ь \ 4°у \ ^ 1» обозначим через 6. Требуется найти такие

решения системы (1) в области в при ^>0, чтобы:

а) при ф -* 0, х < 0 удовлетворялось условие непротекания

Уу = 0;

б) при Ф 0, д: > 0 скорость потока была равна звуковой

<1+гд* + ®* = 1.

В области О систему уравнений Эйлера (1) и краевые условия „а“ и „6“ можно упростить и в первом приближении представить в виде

_ П V) V Л?Х_ , _ о х_ _ о /оч

и + и дг дх + ^ дх ^ —О, (2)

+ о при Ф 0, х < О, (3)

их 0 при -Ь 0, х > 0. (4)

9

Рис. 2

Система уравнений Кармана (2) обладает классом автомодельных решений

= vy = Ф3"—3So(Е). е = .■■■ ■ •

Т '() Y

Если через 9 обозначить потенциал, скоростей, то ср == *рЗл—2 Для определения функций /о(£), &)(£), фоО) получим обыкновенное дифференциальное уравнение

(ф'о - п2'?) + 5я (л — 1) £Фо — 3 (я — 1) (Зя - 2) ф0 = 0. (5)

Если теперь перейти к фазовым переменным

= ,==|п1М.

то уравнение второго порядка для функции /0 сведется к уравнению первого порядка

d4T — 6F— 5nW + 6^2 -f 7 Fff + ¥2

a?/7 ~ (ns-fjy ’ ^

Уравнение (6) имеет четыре особые точки А (0, 0), В{0, 1), С[п2, —/г(/г + 1)] и D[n2y —6п(п—1)] в конечной части плоскости (F, W) и три особые точки, лежащие в бесконечности—-В, <5, Q. Поведение интегральных кривых уравнения (6) и расположение особых точек на плоскости (F, ЧР1), отображенной сначала на полусферу, а затем на круг, показано на рис. 2 при п = 6/5.

2—„Ученые записки" № 2 17

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

/2 \ I 3 13/2 :

(4г+<г]=ч^/;'+'г1 +•••• <7>

Точка В является седлом. Достижение точки В по кривой с отрицательным наклоном соответствует переходу к течению Прандт-ля — Майера. Точки С и О лежат на предельной линии Е = п?. Нижняя из них всегда будет седлом, верхняя — узлом. Они соответствуют предельным характеристикам, выходящим из точки О. Особая точка 6 является седлом. На физической плоскости она соответствует оси у, на которой достигается скорость звука. Ускорение потока при этом бесконечно. Особая точка Е является аналитическим узлом. В ее окрестности интегральные кривые бесконечно дифференцируемые. На физической плоскости она соответствует оси у. Точка С2 является узлом, одной из осей которого является предельная линия. При приближении к точке <3 ускорение потока стремится к бесконечности. При ф -+■ 0, х < 0, т. е. при 5 ->—оо, асимптотическое поведение решения уравнения (5), удовлетворяющего граничному условию (3), имеет вид

/

Ф0 = А0 (- 03~2/Л-(ая-^УСЯ~1) Д2(_ 5)3—4/я 4__ (8)

Решение (8) назовем симметрическим и обозначим через Кх' Оно соответствует условию СА —0. При х > 0, т. е. при

5 --»• + оо, асимптотическое поведение решения (5), удовлетворяющего (4), будет следующим:

Ф„ В„ + 4- ~3) +. . .. (9)

Решению (9) в (7) соответствует значение Сл = со. Назовем его антисимметрическим и обозначим через К2. При п — 4/3 кривая становится прямой 4? =—3/2/7 и попадает в точку О. При всех п£( 1, 4/3) кривая /С1(Е< 0) проходит через точку Е. При п > 4/3 она попадает либо в точку (2, либо в точки С или О. На основании теоремы Брио и Буке в работе [9] было доказано, что существует единственная интегральная кривая /С1(/г<0) при я = 2, проходящая аналитически через эти особенности и точку Е. При 2, пройдя точки О и Еу кривая К\ достигает точки С. Далее, можно с разрывом второй производной попасть в точку А по кривой К2(Р>0), но скорость V будет при этом положительной.

Таким образом, задача (2) —(4) свелась к задаче отыскания таких значений я£(1, 4/3), при которых симметрическая кривая К1(Е<0)у выходящая из точки <4, достигала бы сначала точки Еу в ней переходила бы аналитическим образом в кривую /С2(Г< 0) и по ней обратно возвращалась в точку А. При таком движении вертикальный компонент скорости всюду отрицателен.

Разрешимость такой задачи удобнее всего исследовать в плоскости годографа. Уравнения (2) можно обратить и в результате для функции тока ф == (1 -Ь т)13ф получить уравнение Трикоми [1].

Требуется найти такие его решения Ф>0 в области « = (1+^)1/3Х X ъх < 0, у = ъу<0, чтобы удовлетворялись условия

Ф = 0 при =а 0, (10)

Ф = 0 при м = 0. (11)

Решение этой задачи будем искать в виде

~ ^^У* ^)* ^ 1Р*» 1 — 6(п-1) '

Функция Z^JІ С) удовлетворяет высшему гипергеометрическому уравнению. Одно из его линейно-независимых решений, удовлетворяющее граничному условию (10), имеет вид

ь = о? V (- 4- в»)'"'* Г (-/ + -у + ±, ±; с»), о;жо. (13)

Аналитическим продолжением решения (13) в окрестность

и = 0, ^ < 0 будет

+ —У ч---------2“» —^ т*’ ^)] ’

_ Г (3/2) Г ( 1/3) _ Г (3/2) г (1/3) .

£*11 ---- “ ” 1 /Л\ тГ> • т с, ,ЛЧ * £>12

Г (— / + 1/2) Г О' —|— 2/3) ’ Г(_у + 5/6)Г(/+ 1},

(14)

здесь через Г (а) обозначена гамма-функция, а через Г (а, р, х)~ высшая гипергеометрическая функция. Отсюда следует, что решение (14) будет удовлетворять граничному условию (11), если£12 = = 0. Вспоминая, что я£( 1, 4/3), получим

>=./ + -§-, « = /-0,1,2,....

Таким образом, получился спектр значений п, при которых решения уравнения Трикоми удовлетворяют граничным условиям (10), (11). Обозначим их через ф,. В силу линейности задачи ее

решением будет также любая линейная комбинация фш — X! Ф/ •

1—0

Главный член в сумме Фщ соответствует значению /==0, п = 6/5. В области О течение всюду дозвуковое, поэтому якобиан отображения в области годографа и < 0, г><0 всюду отличен от нуля.

Обратимся теперь к решению с п = 6/5 и найдем связь между постоянными Л0 и Д) в разложениях решения (8) — (9). В плоскости годографа решение при /1 = 6/5 имеет вид

а

Г „.2

(*)

Используя теперь (8) — (9), сразу получим

Форма свободной линии тока дается уравнением

у = 4 рт в0 х™ х>0, Р = (1 н- т)-1;

Интересно отметить, что уравнение свободной линии тока (16) получается таким же, как и в случае обтекания точки излома несжимаемой жидкостью [8, 10].

2. Рассмотрим теперь решение уравнения Трикоми, представленное в виде

¥ =Ф» + ь©*’ + §,- + ФП; (17)

здесь — малый параметр, ^ дается формулой (13) и при ^ = 0, й<0 обращается в нуль. Функция Ф/а) определяется с помощью второго линейно-независимого решения гипергеометрического уравнения для 2Г(/, ч):

ур{—;, -У+4-, 4

. г~~і і

’ " 1

При 'У = 0 функция ф{а)^гО. Решение (17) взаимно однозначно отображается на физическую плоскость и дает решение системы уравнений Кармана. Структура функций $5) позволяет утверждать, что при достаточно малом существует некоторая область Е(и, ъу еу), принадлежащая сегменту и<0, ^<0, в которой I ф0| |Ф/1 * Решение системы уравнеийй Кармана (2), соответст-

вующее решению (17), будем искать в виде

Г./»

рі

Фр,(5)|- (18)

Члены в разложении (18), объединенные в первую сумму, соответствуют решению (17) с Они являются малыми добав-

ками к основному решению Ф8/5Ф0(!)- Так как отображение плоскости годографа на физическую плоскость взаимно однозначно, то существует область О^. с О, в которой члены в разложении (18),

связанные с ф/ .являются малыми добавками к основному решению. Функция ФР/(5) является решением однородного дифференциального уравнения

(-1-м;) - [4-(2р, - -у-) 5 + Фо] ®;,+р,(р, - о ф„,=о.

(19)

В отличие от Фр. функции Ф^ и ФР1 удовлетворяют неоднородным уравнениям (19) и являются их частными решениями. При

*!> -> 0, *<0, т. е. при £ -*■—ос, асимптотическое поведение компонентов скоростей %)х И будет иметь вид

% = АЦ (- г) • (- 5)«^ 4 -

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Р) = - "Г Р/ (- р*)6 " +

+ 4- (Ро -1) л^> (-рЛ)'і‘<р"2Ч +..

©у. р. = ву

•§■ Р) [р> - 4) * (^’"2,4> +

<20)

Положим теперь в решении (17) постоянную £)/а) = 0. На физической плоскости это соответствует тому, что <р будет определять симметрическое решение. Найдем теперь по решению (17) поведение ъх в окрестности образующей АО, которая совпадает при £)/а) = 0 с отрицательной частью оси х. Используя формулы перехода от плоскости годографа к физической, в окрестности АО получим

(“ «)■

4*>)

1/3

(—л;)1'3

4Ї/3 3/—7/6

ЗУ + 1/2 ф<*>)Ж/2

(— х)*-1/2 + •

Сравнение с (20) позволяет определить связь между показателем р; в (18) и значением параметра у в (17). Второе линейнонезависимое антисимметрическое решение уравнения (19), связанное с постоянной А^ в (20), соответствует в плоскости годографа В результате получаем

Рі

(2У + 1).

(21)

Аналитическое продолжение решения (17) в окрестность « = 0, ^<0 позволяет исследовать поведение линейно-независимых решений уравнения (19) в зависимости от значения параметра рОно имеет вид

4 = Ф0 + V» \{ВП йТ + В22 £>5“»] (- СГ Р(- /• + 4-, -/ + -§- , 4-; ■с) +

+И*,г12 + о}‘,)в,,]/г(-у, -/+4-, с)}>

где Вп и Вп даются формулами {И), а

„ Г (I/2) Г (1/3) „ . .. Г (1/2) Г (— 1/2) /от

#21 — 7 рг 7 рГ , #22 “ 7 рр • ^

г(-у+-т)г(/+т) г<-Лг('+т)

При ф->0, я>0, т. е. при |-^4-оо, асимптотическое поведение решения Фр (Е) уравнения (19) аналогично (20)

фр. (?) = в{р. ет ^+4“} (/>г 1} + — ' (23>

Решение уравнения (19) при 5^- — оо и А^ = 0 назовем симметрическим слева, а при 5 -► + оо и 5^=0 — симметрическим справа. Если же £ -»• — оо и — О, то решение назовем антисим-метрическим слева, а при £-*-Н-оо и Вр] — 0—антисимметрическим

справа. Будем говорить, что решение уравнения (19) принадлежит классу (5 -► 5), если симметрическое решение слева переходит в симметрическое решение справа. Точно также определяются классы (5 А), (А 5) и (А ->• Л). Решение принадлежит классу (5

-*5 4-Л), если симметрическое решение слева переходит в линейную комбинацию симметрического и антисимметрического решений справа. Аналогично понимаются классы (5 4-Л->-Л),

(Л-5 +Л).

Обратимся теперь к формуле (22) и рассмотрим случай, когда &\$) В\ъ 4* 0\а) Вгх — 0. На физической плоскости это соответствует Фр. £ (5 4- А -*■ А). Последнее равенство будет выполняться, если

5,2 = 0, £>(/” = 0, &Р ф О, либо 5„ = 0, =0, й^фО. Эти урав-

нения и соотношение (21) вместе дают

/>,. = -§-(3/+ 4), —§-(2/4-1), (/=0,1,2,...), Ф,уС(5-А); (24)

1

Р^Щг’ = (/ = 0,1,2,...), ФР£(А~А). (25)

Далее положим в (22) равным нулю коэффициент

Ви + В,а ОТ' = 0. (26)

На физической плоскости это соответствует случаю + ^ 5)- Равенство (26) будет выполняться, еслиЯц—0,

= 0, О^фО, фО, либо В„ = 0, ор =0, В{рфО. Отсюда

и из (21) получаем

/»,= -!-(/+1), />,= -^М(/ = 0, 1, 2, ...),

^.= 4-(2/ + 1), Л=-§-(2/-1)(/=0,1,2, ...), Ф,уе(А-5).

(27)

Функции Фр.(1), принадлежащие спектру (24), являются собственными решениями задачи (2) — (4). Связь между постоянными Ар} и В(р] задается формулой

Таким образом, в этом пункте статьи изучено поведение фундаментальных решений уравнения (19) в зависимости от значений параметра рПолученные соотношения будут использованы ниже при рассмотрении краевой задачи „а“ и „б“ для системы уравнений Эйлера (1).

3. Ее решение в области С/ будем искать в виде разложений

= Г’5/о (Е) + £' “7* (5), гг, = .р &,({) + £ Ч,Р*~' (У- (28)

Л = 1 А—1 •

Они так же, как и решение (18) краевой задачи (2) — (4), содержат две группы членов. Первая группа, возникающая благодаря нелинейности задачи, состоит из членов с показателями Ръ — **8/5 +2/5 Л (£—1,2,3...); вторая —из собственных решений задачи (2) — (4), каждое из которых дает произвольную постоянную в разложениях (28).

Подставляя (28) в систему (1), получим уравнения для определения функции /А и gk. Они имеют вид

-4 !/«/*] + (/’*-!)£*—г5-Ж-[(л-т)л—

Правые части Рки Рк2 зависят от функций /0, gь и предыдущих приближений. Для удобства исследования введем функцию Ф и в ее разложении возьмем следующие члены:

ф „ ф8,5 ф0 (5) + фВДб Ф1 (Ё) + ъ ф + ^14/5 Хз ф +- (29)

Функция потенциала (18) для системы уравнений Кармана начинает отличаться от разложения (29) со второго приближения. Для определения Ф^!) получаются уравнение и краевые условия

Симметрическое решение слева однородного уравнения (30) принадлежит классу (5^5 + А). Если не положить его равным нулю, то невозможно удовлетворить второму граничному условию, которое должно выполняться на свободной линии тока ОВ. Поэтому в качестве решения уравнения (30) берется его частное решение. Оно имеет следующее асимптотическое поведение:

tfi + О [(— £)-5/3], £ °°»

^ + 0(?-5/2), £-+оо;

здесь аи Ьх — постоянные, зависящие от А0.

Функции второго приближения х2у Л» ё2 находятся из уравнений

(ж «2 -ф;) й - (-1-5 +ф;) й+ж * -

— Р21

^ ® f 1 . X422 Хс~аГ“Г

lim ф7;5 g2 = 0, Urn ft [2^Хг + Wfbfx -f ё1\ = I

'ji^-0, дг<0 <L-vO, .*>0 j

Симметрическое слева решение однородного уравнения (31) принадлежит классу (S S). Асимптотическое поведение Ха (£)♦ удовлетворяющее первому граничному условию, дается формулой

/л —•

Чтобы удовлетворить второму граничному условию, необходимо ПОЛОЖИТЬ Вг):$ — — В^(40).

Для отыскания х3, /з» £з будем иметь

(-1-52 -®;)й - (-§-6 +ф;)й+¥х.■- рз.-[4

^32 _ р ^ _ <*Хз „ __ 14 б , £*Хв , п . ■

^ *32» /з— ^ 1 £з — “5— Хз 5 5 £ Ц|32>

Иш ф9 5 g3 = 0, 11т ф8''5 [2р/3 + р2 (/2 + 2/0/г) + 2^0 ^1 = 0.

^->0, *<0 ф-*0, .*>0

Симметрическое слева решение Хз принадлежит классу (5->Д) и является собственным решением задачи (1), „а“ и „б“. Выпишем его асимптотическое поведение

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

у=>--f Р1'2 в, л| р5» [в$5 + 4- So] ^ + ■ ■ ■ ■

4. Рассмотрим случай обтекания точки излома, когда образующая АО задается уравнением у = f{x). Считается, что в области G функция f{x) раскладывается в ряд Тейлора

00

+ ^2<0.

ft=0

Чтобы удовлетворить условию непротекания на АО, необходимо в разложение (29) ввести члены yf>l^pl{\) с показателями Рг = (6/5)/5 (^ = 1, 2, ...). Функции 1(.р1 с такими показателями принадлежат спектру (25), a £ (Л Л). Это означает, что симметрическое слева решение нужно положить равным нулю.

В заключение автор благодарит С). С. Рыжова за обсуждение статьи.

ЛИТЕРАТУРА

1. Гу дер лей К. Г. Теория околозвуковых течений. М., Изд. иностр. лит.. I960.

2. Vaglio-Laurin R. Transonic rotational flow over a convex corner. „J. Fluid Mech.*, vol. 9, N 1, 1960.

3. Есин А. И., Чернов И. А. К вопросу о локальном течении вблизи выпуклого угла. Аэродинамика. Межвузовский научный сборник, вып. 6(9). Изд. Саратовского ун-та, 1978.

4. Диесперов В. Н. Трансзвуковое течение разрежения в окрестности выпуклого угла. ПМ.М, т. 45, вып. 4, 1981.

5. Нейлаид В. Я., Сычев В. В. Асимптотические решения уравнений Навье— Стокса в областях с большими локальными возмущениями. „Изв. АН СССР, МЖГ“, 1966, № 4.

6. Н е й л а н д В. Я. К теории отрыва ламинарного пограничного слоя в сверхзвуковом потоке. „Изв. АН СССР, МЖГ“, 1969, № 4.

7. Сычев В. В. О ламинарном отрыве. „Изв. АН СССР, МЖГ*,

1972, № 3.

8. Рубан А. И. О ламинарном отрыве от точки излома твердой поверхности. „Ученые записки ЦАГИ“, т. V, № 2, 1974.

9. Лифшиц Ю. Б., Рыжов О. С. Об асимнтотическом типе плоскопараллельного течения в окрестности центра сопла Лаваля.

ДАН СССР, т. 154, № 2, 1964.

10. Akkerberg R. С. Boundary-layer separation at a free stream line. Part L Two-dimensional flow. „J. Fluid Mech,1*, vol. 46, pt. 2., 1970.

Рукопись поступила 26jfX 1980 г.

i Надоели баннеры? Вы всегда можете отключить рекламу.