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

книги / Потенциальная теория неустановившихся сверхзвуковых течений

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

где интеграл понимается

в

 

ния Коши.

 

 

 

Интегральное уравнение (9.2.15) решается с помощью

тригонометрической

подстановки

 

I =

CosQ,

I ' = COsQ'

(9.2.16)

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

 

■^o(cosQ)=

rt“^anSin(/zQ),

(9.2.17)

и использования следующего интегрального соотношения

COS (nO') d

sin(«0)

(9.2.18)

(cosO' —COS 0)

SinO

 

Этим путем можно показать, что коэффициенты а,^в разло* жении (9.2.17) имеют вид (9.2.6).

§9.3. Решение уравнения Гельмгольца

вэллиптических координатах

Задача, формулируемая соотношениями (9.1.5,6,9), а именно двумерное уравнение Гельмгольца с заданными граничными условиями на обеих сторонах разреза т) = О ± ,

|||<1, эквивалентна задаче дифракции на ленте. Эта последняя была впервые рассмотрена Релеем который воспользовался методом {релеевское рассеяние), нашедшим впоследствии свое применение в аэродинамических задачах (Уорд 1^®^]) под названием теории обтекания тонких удли­ ненных тел. (Метод Уорда, в свою очередь, был использован Майлсом [*^®1 в задачах дифракции). Зигер еще в 1908 г.[®“®1 дал решение этой задачи в эллиптических цилиндрических координатах, однако численные результаты были получены лишь тридцать лет спустя (Морзе и Рубинштейн I®®**]). Метод сведения задачи дифракции к интегральному уравне­ нию использовался как для получения приближенного (Майлс [^'®1), так и точного (Майлс Р®°]) ее решений. В этом и следующем разделах будет рассмотрено применение последних двух методов. В эллиптических цилиндриче­ ских координатах (9.2.8) уравнение (9.1.9) примет вид

(Магнус и Оберхеттингер [^^], стр. 157)

“Фиц + ^vv+ (ch^р.— COS-v) Ijj = О,

(9.3.1)

граничное же условие на крыле (р = 0) будет задано соотношением (9.2.9). Разделяя переменные, найдем, что решением уравнения (9.3.1), удовлетворяющим условию (9,2.9) и ведущим себя как расходящаяся волна, будет

 

 

(9.3.2)

^ v(cosv)

(V , <7)sinvrfv,

(9.3.3)

9 = 4-?^^

(9.3.4)

и функции Nej?’ и se,^ являются рещениями

уравнения

Матье, аналогами которых в полярных координатах будут функции п ’ (^ г) и sin (/10). Здесь используются обозначения, принятые в монографии Маклеклэна [“ ’ ] для функций' Матье; в частности, нормализация этих функций проведена

как в

§ 2.21 этой монографии.

 

В

практике нас интересуют значения потенциала ф на

крыле (р = 0), где можно записать его в виде

 

 

фо (cos v) =

(к) SB^ (v, q).

(9.3.5)

Здесь

 

 

 

 

F п \ -

Ne<?»(0. q)

(9.3.6)

 

 

N e '« '(0 , «7)*

 

 

 

При вычислении интегралов (9.3.3), a также интегралов от фо, удобно представить фупкцшо se,j в виде разложения в ряд Фурье (Маклеклэн, 2.17)

SCrt (V, <?) = 2

Sin (mv).

(9.3.7)

m=l

 

 

в вычислениях функций

полезно

воспользоваться

соотношениями 1) (см. там же,

13. 40 (И ),

(14)1:

NeS..

Ч) = C,,,.. S ( - ) '

(я) X

 

 

г=0

 

 

X [

Я'" (I Яе») -

J,,. ( i Хе-») H f ( i

Xei*) ]

 

 

 

(9.3.8)

NeS.. (9. Я) = С.„., S ( - ) ' ■SS"" (9) X

 

 

г=0

 

 

 

X [ A d X e - . * ) ^ ^ ^ ) * ^ ! * ) -

 

 

- y , . , ( i X e - i * ) w t “ ( l x e i * ) ] ,

(9.3.9)

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

го, связывающих функции Бесселя (Ватсон

3.63)*)

получим

 

f .„..(>•) = ( i ( - г Bflii" (я) {(2г + 1 ) + -^^^ X

г=0

 

X [^.(ix) я?>( i х) +

(1 х) я«> (| х)-

( т 0

] } ) ~ ‘ 2 ( - ) 'B S - :i" ( 4 ) .

 

T=O

 

(9.3.10)

1)Разложение ())ункций Ne^?i42. Данное в работе Мербта и Ландала [^^“), содержит небольшую ошибку.

2)После номера моногра(}жн Ватсона в скобках указан номер соответствующего раздела.

( ^ ) {2 i-)4r+i)KT(4) X

r=0

 

X[•'<■?O ( T O +

C^ O-

- 2 / „ , ( | я , ) я " ; ( 4 x ) ] } - ' 2 ( - ) Ч г + i j B S r t o ) -

 

 

r=0

 

 

 

(9.3.11)

Для численных расчетов полученные выражения можно

использовать в

том

виде, в каком они здесь даны,

или

функции Б ‘г"* и

можно разложить в ряды по степеням q

(причем ряды

для

будут содержать логарифмические

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

§ 9 .7 .

Сходимость рядов является, в общем, удовлетворительной для значений Я, меньших 8. Можно ожидать, что функции

F-n (^) будут иметь полюса в комплексной

плоскости Я (ср.

C (9.7.25)1, которые приведут к сильной неустойчивости

результатов при малых значениях Я, как

это имело место

у Морзе и Рубинштейна

 

При больших значениях Я сходимость ряда (9.3.5) можно улучшить, выделив из него предельное значение фо[приЯ-?-оо, которое равно

Iim фо =

(cosv) =

6л Se,I (v, 9),

(9.3.12)

K =

\ o(cosv) se,j (v, 9) dv.

(9.3,13)

 

0

 

 

Тогда выражение (9.3.5) можно переписать в следующей форме:

Фо (cos v) = (1Я)'1 V(cos v) + '2 [а,^Fn (Я) - 6,J se,^ (v, 9).

в качестве примера численного расчета, иллюстрирую­ щего полученные выше результаты, мы. найдем здесь иите-

Р и с . 9 . 3 . а) Д е й с т в и т е л ь н а я ч а с т ь и н т е г р а л а п о р а з м а х у X , н о р м и ­ р о в а н н а я п о т о ч н о м у з н а ч е н и ю п р и Я = 0 ; б ) м н и м а я ч а с т ь и н т е ­ г р а л а п о р а з м а х у X i н о р м и р о в а н н а я п о т о ч н о м у з н а ч е н и ю п р и Я -+.сзо

грал от потенциала по размаху крыла в важном частном случае, когда величина v постоянна по размаху. (Не умень­ шая общности, можно принять, что P = 1, а затем конем*

ный результат умножить на действительное значение v, причем последнее может включать также й зависимость от координаты х.) Значительно более обширные результаты численных расчетов приводятся в работах Мербта и Лан-

дала

Милна [--°] и Мазельского

 

Полагая в выражении (9.3.3) у = 1

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

разложение

(9.3.7), получим

 

 

(q).

(9.3.15)

Теперь подставим а,^т (9.3.15) в разложение (9.3.5) и про­

интегрируем

последнее по размаху

крыла,

что

даст нам

X = ^

5 il’o sin Vrfv =

я

3

(7)]®

(^)-

—1

о

 

п=1

 

 

 

 

 

 

 

(9.3.16)

Действительная и мнимая части х> нормированные так, чтобы прн X = O и OO они соответственно имели единичные значения, показаны на рис. 9.3, а, б, где они сравниваются C результатами различных приближенных расчетов, рас­ сматриваемых в следующих разделах.

§ 9.4. Решение уравнения Гельмгольца методом сведения к интегральному уравнению

В качестве еще одного метода решения граничной зада­ чи, сформулированной в предыдущем разделе, рассмотрим метод сведения к интегральному уравнению (§ 9.215).

Решение уравнения (9.1.9), ведущее себя на бесконеч­ ности как расходящаяся волна и обращающееся в фо(1) при 11 = 0 + » имеет вид

IMl. Ч) = ( Щ Ь (Г) d f, (9.4.1)

где

1

R =

(9.4,2)

To, что выражение (9.4.1) удовлетворяет уравнению (9.1.9), следует из равенства

 

 

 

(9.4.3)

и из того факта, что функция

("KR) является фундамен­

тальным решением уравнения Гельмгольца. А то,

что оно

обращается в ф(^, 0 + )=ф о(|),

следует из того,

что

(XR)

 

(Х/?->0).

(9.4.4)

2iR

' nR^

 

 

Поскольку множитель в подъинтегральмом выражении не обращается в нуль в окрестности Ii = O-I- только если стремится к I (т. е. i? —^0), то выражение (9.4.1), при нахож­ дении интересующего нас предела, эквивалентно выраже­ нию (9.2.11). Как и в (9.2.11), пределы интегрирования по

( — 1, -1-1), вытекают из условия (9.1.6).

 

 

 

 

Дифференцируя выражение (9.4.1)

 

по

т|

и

полагая

T j- >0-1-, после подстановки в-условие

(9.1.5)

получим

соотношение

 

 

 

 

 

1

 

 

 

 

 

= V(I)

( П

К

1).

(9.4.5)

-1

 

 

 

 

 

которое является интегральным уравнением для фо с ядром (9.4.6)

Для удобства расчетов из выражений (9.4.5) и (9.4.6) удобно выделить сингулярные части (как это было сделано

в (9.2.13— 15)). Тогда

будем иметь

 

 

1

 

 

 

 

^Ч>«+( S i f t - l ') 4 > ,f t ') < f t ' =

<'(l)

(| | | < 1),

(9.4.7)

где символ L определяется из равенства (9.2.15), а функ

задана уравнением

 

 

 

 

g iis ;

2i\i\

 

^

(9.4.8а)

+ я | * ’

 

 

 

IiD V d V

(9.4.8b)

Иначе уравнение (9.4.5) можно получить из выражения

оо‘ 1

■“

(9.4,9)

которое, очевидно, является решением уравнения (9.1.9), обращающимся, в силу интегральной теоремы Фурье, в 1|'о при TJ = О H- 1). Условие конечности решения при т) = со требует, *1Тобы действительная часть радикала была неот­ рицательна. Если же действительная часть обращается

Рис. 9.4. а) Комплексная плоскость v при 1т(Я) < О, показаны

разрезы ]^v®—Х,®(см. (9.4.10)); б) путь интегрирования в выражении (9.4.11) при 1т(Я) = о— .

В нуль, то из условия радиации следует, что аргумент равен я/2. Следовательно, разрезы, выходящие из v = ± Я, должны быть линиями, на которых действительная часть радикала обращается в нуль и при переходе через которые она меняет знак. Налагая эти условия в предположении, что X имеет небольшую отрицательную мнимую часть (кото­ рую затем мы устремим к нулю), найдем, что разрезы должны быть отрезками равнобочной гиперболы:

Re (V) Im (V) = Re (X) Im (Я.),

(9.4.10)

как это показано на рис. 9.4, а.

 

^) Вид решения (9.4.9) предполагает, что ii> 0 ;

если л<0»

то знак перед решением и под радпкалрм р (9.4-Э) следует изменить.

Дифференцируя вырал«ние (9.4.9) по ц и подставляя результат в граничное условие (9.1.5), получим интеграль­ ное уравнение (9.4,5), где ядро теперь, уже представлено в интегральной форме

 

 

=

\

 

 

(In i(M < 0 ) .

(9,4.11)

Теперь,

 

если при ^ > О контур

интегрирования

проходит

вокруг

разреза

в

верхней

полуплоскости,

а при

| < О

вокруг

разреза

в

нижней

полуплоскости,

1 т (Л )—^►О— ,

как показано на

рис. 9 .4 ,6 ,

то

выражение (9.4.11)

путем

введения

соответствующего

интегрального

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

функции

Ханкеля

(Ватсон

1®“^]

6.11 (7)1

преобразуется

к виду

 

 

 

 

 

 

 

 

 

 

 

6 vs;

 

2i I ^ I

 

(Im (X) =

O - ) ,

(9.4.12)

 

 

 

 

 

 

 

 

что согласуется с (9.4.6).

Решение интегрального уравнения (9.4.5) введением тригонометрической замены переменных (9.2.16) можно

представить в виде ряда Фурье

типа (9.2.17)^),

так что

 

Фо(со5 0 ) = 2

п

" sin (л0)

 

(9.4.13)

2

\ ё (cos0 — COS0') sin (л0') sin 0' dQ' = v (cos 0).

71=1

0

 

 

 

 

 

 

 

 

 

 

 

(9.4.14)

Умножая

обе части равенства

(9.4.14)

на

(2/лт) sin (/л0)

sin 0 и интегрируя от О до я,

получим

 

 

 

 

f!c,„„a„ = m - W

( m = l, 2,

3,

. . . ) ,

(9.4.15)

Если решение представлено в виде разложения в ряд по функ­ циям Мэтье, то из иптегр.зльного уравнения непосредственно полу­ чаем решение в виде (9.3.5).

где Om определено соотношением (9.2.6), а

я я

о IiО ®■" ^

X sin(/i0')sin0'dO 'd0. (9.4.16)

Очевидно, что с„,„ = Если в выражении (9.4.16) выделить сингулярную часть

функции g или просто подставить ряд (9.4.13) в выражение

(9.4.7), то выражение (9.4.16) примет вид

 

тп =

 

 

^ ^ sin (ш0) sin

(cos 0 — COS0') X

 

 

 

0. о

 

 

 

 

 

X sin (n 0 ')sin 0 'd 0 'd 0 ,

 

(9.4.17)

где символом

обозначена дельта-функция

Кронекера

(6}Г=1

при т = п]

 

=O при т ф п ).

 

 

При

вычислении

величины

удобно в

выражение

(9.4.16)

подставить

интегральное

представление (9.4.11).

и тогда выражение примет вид

 

 

 

=

 

 

.ivcosf)

('” 0)

sin 0 d0 X

 

 

5 е»-

 

 

 

 

Sin(«90 sin 0'd

(9.4.18)

Интегрированием по частям интегралы по 0 выражаются

через функции Бесселя и тогда

формулы (9.4.18) приобре­

тают простой вид

 

 

 

т . = 2 ( - )

" f

I.(V)

(V) ^ ,

если

т — п четное,

(9.4.19а)

,^ = о, если

т — п нечетное.

(9 .4 .19Ь)

В пределе, когда А,= 0, выражение (9.4.19) сводится к про­ стому c„j„ = /n'^6rt (см. Ватсон 13.41 (7) и 13.42(1)), P соответствии с формулой (9.4.17). Сходимость интеграла

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