книги / Потенциальная теория неустановившихся сверхзвуковых течений
..pdfгде интеграл понимается |
в |
|
|
ния Коши. |
|
|
|
Интегральное уравнение (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). Сходимость интеграла