книги / Потенциальная теория неустановившихся сверхзвуковых течений
..pdfдля вычисления давления, задаются выражениями
|
|
1 |
|
|
|
|
|
|
: ^ |
= 2 |
^ 1|,<»-М| = 2 т - - ^ т = |
(9.7.10а) |
|
|
|
|
X<c“>= 2 - t . |
|
(9.7.10b) |
|
Если |
произвести |
наложение влияний |
обеих |
кромок |
||
для 1 < |
т < 2 |
(см. § 7.5 и |
рис. 7.5), то выражения (9.7.10) |
|||
оказываются |
справедливы |
во всем интервале 0 < г < 2 . |
||||
Граничная задача для потенциала |
получающаяся |
в результате подстановки выражения (9.7.9) в граничное условие (9.7.7), может быть решена как методом Эвварда (см. § 7.3), так и путем наложения конических течений. Выкладки, ведущие к цели, несколько длинны, но окон чательный результат, т. е. интеграл от потенциала по размаху и его первая производная по времени (см. Лагерстром [^^°]), имеют вид
= 4 |
- |
(т^ - о |
1 |
+ |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
+ i ( |
1 + 4 т ) (Т»- 4)*+ (I - т) arch ( 4 |
т)| |
(9.7.11а) |
||||||
И |
|
|
|
|
|
|
|
|
|
*"’= 7 |
{ ( T |
- о |
T + т |
У |
( |
т |
^)} |
||
|
|
|
|
|
|
|
|
(9.7. Ilb ) |
|
Величина |
|
обращающаяся в нуль при т = |
2, должна |
||||||
быть сложена с х‘“* при 2 < |
т < 4 . |
|
|
|
|
||||
Переходя к асимптотическому разложению для случая |
|||||||||
больших X, |
но сохраняя |
предположение, |
что зависимость |
||||||
скоса потока |
v от времени |
задана |
соотношением |
(9.7.8), |
|||||
рассмотрим |
следующую |
граничную задачу для |
преобразо |
||||||
вания Лапласа: |
|
|
|
|
|
|
|||
|
|
|
1 И | ,_ о = - а - ‘ |
( I I K l ) , |
|
(9.7.12) |
|||
|
|
|
r U - o = |
0 |
(| | | > 1). |
|
|
(9.7.13) |
|
|
|
|
1|>й + |
'И а -а > * --= 0 . |
|
|
(9.7.14) |
Очевидно, что решение этой задачи может быть получено из результатов § 9.3 простой заменой параметра X на —is
(посколь1{у |
заменено |
на е” ). Подставляя |
в |
выражение |
|||
(9.3.3) V — S"^ |
и используя |
разложение (9.3 |
.7), |
получим |
|||
а,^ = |
^ |
^ |
(п — нечетное), |
(9.7.15) |
|||
откуда (так как при четных /г коэффициенты |
обращаются |
||||||
в пуль) из (9.3.5) для ф* следует выражение |
|
|
|
||||
H)»(cosv) = |
s - i 2 |
|
|
|
|
|
|
|
«=0 |
|
|
|
|
|
|
х |
( |
^2.4.1 ( - « ) s e 2 „ , i (^v, - J S ' - * ) |
|
(9.7.16) |
Вычисляя интеграл по размаху крыла от выражения (9.7.16), C помощью разложения (9.3.7) получим
Х» = |
4 |
- я 5 - 1 2 [ s ‘r |
“ ’ ( - T S ' “) ] * f , „ . , ( - i s ) . (9.7.17) |
||
|
|
п=0 |
|
|
|
Влияние изменения знака у параметра |
<7= — |
рассмот |
|||
рено |
в |
монографии |
Маклеклена |
2.18 |
и 13.41, сле |
дуя результатам которого выражения (9.7.17) и (9.3.6,10) примут вид (см. Майлс 1“®^1)
Х* = 4 - ж ' ‘ 3 |
(9.7.18) |
|
^ 2 n + l( - « )= |
- |
(9.7.19а) |
|
Gekjn+i ^0. — 4^**^ |
|
= |
( 2 |
4 ? f f ' ( | 5 = ) { ( 2 r + i ) + |
(9.7.19b)
Коэф4)ициенты Л м о ж н о разложить в степеппой ряд
(^no аргументу ~ |
, радиус сходимости |
которого (см. Мак* |
|
леклен |
2.17) |
будет не менее ^/g, а |
опыт показывает, |
что ряд сходится в радиусе, близком к 4. Искомые разложе ния, полученные с помощью соотношений (2.151), (2.21),
(3.10) и |
(3.33) |
цитированной |
монографии, |
имеют вид |
|||
|
|
|
|
|
|
|
(9.7.20) |
.+,5-.L , (Д. 0+ (| - - • - ) ( ! = ) + ( ' Ц - ¥ -•) ( т |
OVo [ (| .) ‘ |) • |
||||||
|
|
|
|
|
|
|
(9.7.22) |
^2ri+l ( |
^®) = |
|
|
|
|
|
|
|
|
1 + (п»+-)-‘ ( 4 0 |
- ю |
[ ( | г У |
] |
|
|
= ( 2 r t + i r |
, + з ( „ . + „ ) - . ( | , ) + о [ ( i . y j j |
( о > 1 ) , |
|||||
|
|
|
|
|
|
|
(9.7.23) |
|
|
2 |
= -^s®, |
L = y + ln(^-^s^ |
(9 .7 .2 4 а ,b) |
||
и у — постоянная |
Эйлера. |
|
|
|
|
||
Предполагается, что каждая |
из функций |
имеет |
по крайней мере один полюс, так что интеграл %* имеет их бесчисленное множество. Но приведенные выше разложе ния (9.7.20—23) служат лишь для определения с приемлемой точностью наименьшего из этих полюсов, а именно
G e k ; ( |
о |
. |
.=0, |
|
Ч |
|
l3=rSo, во |
|
|
5„= ^ 0 ,6 8 2 + |
11,114= |
1,306е'2>120^ |
(9.7.25) |
где черточкой сверху обозначена комплексная сопряженная величина. Выделяя из выражения для %*(5) сингулярные части (особенности), связанные с этими полюсами, и, пре небрегая особенностями более высокого порядка, мы можем разложить выражение (9.7.18) в ряд по возрастающим сте пеням параметра s. Этот ряд имеет вид
X*(s) = 4-«»“‘ |
3 |
( v ) ”] + |
|
Wi=O |
|
|
4 |
|
+ ^-1 |
V 2 |
J ^^ |
+( ^ - т ) ( 1 0 ' +(^'-Т^ +^)(т^> +О(^1п’4 |
* |
|
|
(9.7.26) |
|
где символом а _1 обозначен вычет, связанный с полюсом |
^), |
|
который равен |
|
|
а_1 = 0,760ё-»о-72з_ |
(9.7.27) |
Обратное преобразование выражения (9.7.26) дает асимпто
тические |
разложения |
|
|
|
||
X (T) ~ 4 |
л { |
1 + |
2Re t0 .,s; v » 4 + 4 |
т'^ - |
|
|
|
- |
I |
[|п (4т)— g - ] х-- + |
0(т-«1п «т)} |
(9.7.28а) |
|
И |
|
|
|
|
|
|
Xt(T) |
4 я |
{ 2 Re |
- 4 т ' “+ |
|
|
|
|
+ 3 [ In (4т) - |
.g. ] т-“ + О (Т-’ 1п> т)} |
(9.7.28Ь) |
|||
Хотя члены порядка |
при т > 1 , |
входящие в выражение |
||||
(9.7.26), |
при г > о не влияют на асимптотические разложе |
ния (9.7.28), следует отметить, что они могут влиять на эти разложения, если функция v* содержит члены типа 5 '“,
и т. д. (соответствующие членам тип? t, t- и т. д. в выра жении для функции о). Отметим также, что хотя, в принципе,
1) При вычислении а.] учитывались член
МЫ никак не обосновали здесь сохранение экспонен циальных членов и в то же время пренебрежение членами порядка t " “ ln*t, но расчеты показывают, что экспоненциаль ные члены доминируют в асимптотических приближениях на большей части диапазона, где эти приближения приме няются.
Результаты численных расчетов по формулам (9.7.10, 11,28) приведены на рис. 9.6. Мы видим, что приближение
Рис. 9.6. Интеграл от потенциала по размаху для единично го скачка скоса потока в момент времени т = 0.
(9.7.28) в пределах точности графического изображения
совпадает с более точным результатом |
,+ х "’ |
во всем |
||
диапазоне 2 < t < 4 . |
Следовательно, |
|
необходимо |
вычис |
лять лишь функцию |
при малых |
значениях х |
(т < 2). |
|
В работе Майлса |
приближение Кирхгофа, соответ |
|||
ствующее преобразованию Лапласа |
от выражения (9.5.3), |
сравнивалось с вышеприведенными результатами, и была показана его непригодность за исключением очень малых
значений т, где более точный результат |
проще прибли |
жения Кирхгофа. |
|
§9.8. Крылья квазималого удлинения и метод исследования низкочастотных колебаний
при околозвуковых скоростях
Результат теории крыла малого удлинения (§ 9.2) может рассматриваться как первый член разложения решения граничной задачи о крыле в ряд по степеням удлинения (см, примечание на странице 156),
Адамс и Сирс (*], (Ч получили следующий член этого
разложения для |
случая |
стационарного обтекания ^). |
||
В этом разделе мы дадим |
обобщение |
результатов |
Адамса |
|
и Сирса на случай нестационарного |
обтекания, |
предпо |
||
лагая, что ArM=O(I). Еслиже ЛМ > 1 |
и удлинение крыла |
|||
не очень мало, то |
задача |
будет уже |
относиться |
к слу |
чаям В[ или В'1 из таблицы |
1. |
|
|
|
Наши рассуждения здесь будут весьма сходны с рассуж |
||||
дениями § 9,4, но, |
поскольку мы теперь сохраняем |
произ |
водную по д: в дифференциальном уравнении в частных производных, вместо (|, Ti) следует воспользоваться коорди натами {у, г). Будем исходить из дифференциального урав нения (3.3.4), предполагая гармоническую зависимость от времени (3.10.2).
Тогда =)
Tvt, + |
Ф .. = Р Ч х |
+ |
2»/гМ“Ф , - |
( 9 .8 . 1 ) |
|
Преобразование |
Лапласа |
от |
этого |
уравнения |
имеет вид |
|
ф5у4-Ф?г — |
= |
|
(9.8.2) |
|
где в отличие от выраокения (9.1.10) |
параметр Я определен |
||||
соотношением |
|
|
|
|
|
|
X2 = p V = 2 tW s -/ fe W .. |
(9.8.3) |
Уравнение (9.8.2) формально идентично уравнению (9.1.9), если в последнем Я предполагать мнимым, а переменные (|, г|) заменить на (у, г). Тогда, по аналогии C (9.4.1), решение уравнения (9.8.2), обращающееся в ф* (s, у)
при 2 = |
о - f и в нуль на бесконечности, молено представить |
|||||
в виде |
|
|
|
|
|
|
|
|
Ф* (S, у , |
2) = ^ |
5 ^ '" ^ 1 |
|
|
1) Дальнейшее |
развитие |
этого направленн |
в работах Лен- |
|||
дала |
[**•) |
и Эйхельбремнсра (®®]. |
|
уравнения |
||
-) |
За |
исключением специально оговоренных случаев, |
||||
этого |
раздела могут быть |
записаны или через <р |
и о, |
или через |
комплексные амплитуды ф п v, поскольку множитель, свидетель ствующий о гармоническом характере зависимости от времени, можно либо сохранить, либо опустить.
где теперь уже
и интегрирование по "п проводится в интервале (— оо, + оо). Последнее обстоятельство обусловлено тем, что потенциал Фо (s. "П) зависит, в общем случае, от фо (^, л) при всех зна чениях (Пределы интегрирования в выражении (9.8.4) могут быть заменены на ± бтах, однако пределы ± со не менее удобны.) Если выражение (9.8.4) продифференциро
вать по 2, а обратное преобразование Лапласа |
от |
получен |
ного результата приравнять — v при Z = О, |
(у) |
< Ь{х), |
то получится интегральное уравнение, аналогичное урав нению (9.4.5).
Однако, как и в уравнении (9.4.7), целесообразно сначала выделить сингулярную часть ядра, и тогда, опуская
нулевой |
индекс |
при |
значении |
потенциала |
ф на |
крыле, |
получим |
|
|
|
|
|
|
+ |
J g*i{s, |
y — r\)(p*{s, 7\)dr] = v, |
\y\< b{x), |
|||
|
|
|
|
|
|
(9.8.6) |
a L — оператор, |
определенный |
соотношением |
(9.2.15) |
и соответствующим образом видоизмененный применительно к рассматриваемой задаче:
|
1ф = |
% (X, Ч) 4т1 |
|
(9.8.8) |
|
. (У— П) |
|
||
|
|
|
|
|
Интегральное уравнение (9.8.6) можно решать методом |
||||
итерации подобно (9.4.24) ^). Тогда первая |
итерация |
будет |
||
1) Чтобы получить полное и систематическое |
разложение реше |
|||
ния по степеням |
полуразмаха |
6, можно, вообще |
говоря, воспользо |
|
ваться методом, |
описанным |
соотношениями (4 .4 .28 — 31), |
однако |
|
тогда нужно будет включить наряду с членами |
вида |
и члены |
||
Iog Ь. |
|
|
|
|
где <p‘“’ — решение, полученное по теории крыла малого удлинения ILcp" = v\. Первая итерация (9.8.9) дает решение с точностью до членов порядка о\х-Ь^), если только функция g* аппроксимирована первым членом разло жения в ряд выражения (9.8.7), т. е. имеет вид (ср. 9.4.25)
g',(y)= — |
+ |
( т »'1!'| )] + 0 [x y in (» .|< /|)]. |
|
|
(9.8.10) |
Подставляя это выражение в формулу (9.8.9) и ввделяя член C InX, после перемены порядка операций интегриро
вания пот] и преобразования Лапласа получим ^)
-Ах)
Ь(Л) |
|
|
X \ |
1]) |
(9.8.11) |
- Ш |
|
|
Пользуясь результатами § 9.2, можно показать, что инте грал от по хорде стремится к нулю, при х - ^ 0 , как Ь\х), а следовательно, по условию (9.1.1), как х®. Обратное преобразование первого члена (9.8.11), таким образом, будет иметь вид
Нх)
-Ь(*)
Ь(х)
= ^ I [ - Т + 1 п (у | у - 1 1 | )]ф "’Ч ^ Ч )‘'Л- (9.8-12)
-Н х )_______
Заметим, что при получении решения в ряде случаев бывает целесообразно сохранить функцию в форме (9.8.11) и выполнить операции преобразования Лапласа после проведения необходимых интегрирований по хорде крыла с помощ,ью преобразования Лапласа.
где символом Л обозначен оператор
Рассмотрим обратное преобразование последнего члена выражения (9.8.11). Обозначим
Ь(.х) |
|
|
X" (■«) = \ |
(^1 л) |
(9.8.14) |
-« *) |
|
|
Принимая, что P"i = 0 ( l) и пользуясь соотношением (9.8.3), разложим InX в ряд
1пХ = 1п(М + у 1 п [ |
] + |
|
+ 4 l n [ ^ ± i ^ M |
± l T L ] |
(9.8.15) |
Теперь, применяя теорему Фальтунга и ссылаясь на монографию Кемпбелла и Фостера [^®] (893, 894,2), можно показать, что
i?-4?^My + ln(Ps)]^X<“)} = |
|
|
= AiS-^ { |
], 5д а о )| _ |
(9 8 .16а) |
X |
|
|
= Л [(1пР)х<">-| 5 In(X-I)Xhhdldl] |
(9.8.16b) |
|
Jg-I |хЧ п |
j;x i» )} = |
|
= |
( 1 - е х р [ - « Ш ( М ± 1 ) - ‘ ( * - | ) 1 } - ~ (9.8.17) |
откуда
-Н х) |
^ |
^ |
|
X |
|
X Ф<»'(X,Ч) йлIr S (*- 1)х<« (У <И+
О
+ J [1 - |
V-6) COS (Л мр'" {X - 1))1 |
I . |
|
|
(9.8.18) |
Если в выражение (9.8.18) подставить разложение для функции ф<®) из (9.2.17) и выполнить интегрирование по ij, получим
оС О ..»+ 1 л { [| п (I f » ) - у |
] ft*(х )< ’ (X)- |
|
о |
|
|
JC |
|
|
+ J l l - e-i'‘M23-8(,:_6)cos (ШР*2 (^x - 1))] |
^ |
|
m=l |
|
|
X COS |
arccos-^^J I |
(9.8.19) |
Решение интегрального уравнения (9.8.9) с правой частью заданной выражением (9.8.18) или (9.8.19), получается теперь непосредственно по аналогии с решением, данным в § 9.2, в виде
ф(‘)(д?, у) |
(х) 2 |
(л) Sin [ п arccos |
(9.8.20) |
n=5sl
где
п
= |
J *)(;:, ^COS G)sin (п0) SinGлЮ. (9.8.21) |
|
о |