книги / Потенциальная теория неустановившихся сверхзвуковых течений
..pdfПри |
этом |
гиперболическ |
|
преобразуется |
||
к |
виду |
|
|
|
|
|
|
|
|
|
Д'ф = 2ф £,1-ф „. |
(3.7.5) |
|
§ |
3.8. |
Общее гиперболическое |
преобразование координат |
|||
|
Рассмотрим вкратце преобразование координат следую |
|||||
щего |
вида: |
|
|
|
|
|
x' = x'{q„ |
q^), |
у = IJiq^ |
^7г. <7з),5 |
z = z(<7i, q., |
||
|
|
|
|
|
|
(3.8.1) |
при котором линейный гиперболический элемент изме няется согласно соотношению
{ds'Y = |
{ d x ' f - { d i j f - { d z f = |
|
= |
{K dq^ ^ -{iudq^ Y -{h^ dq.;)\ |
(3.8.2) |
где (/ii, /ig. ^з) есть набор действительных положительных коэффициентов, определяющих метрику (см. Курант и Гиль берт [^’ ], т. 2, стр. 433). Мы можем найти преобразованный вид гиперболического лапласиана (AV)i при соблюдении условия (3.8.2), путем сравнения с преобразованием Ляме (см. Уиттекер иВатсон I®"’ ], стр. 401) обычного лапласиа на. -В преобразованном виде гиперболический лапласиан примет вид
=CssrJ
Преобразования, определяемые соотношениями (3.8.1—3), представляют интерес главным образом в связи с реше нием гиперболического уравнения Гельмгольца (3.9,2) мето дом разделения переменных. Классическое (эллиптическое) уравнение Гельмгольца допускает разделение переменных лишь в одиннадцати системах координат (см. Робертсон [2« ] и Эйзеихарт [^^l). Однако Хейс в 1957 году указал, что благодаря анизотропности оператора Д' (в отличие от изо тропности классического лапласиана) в гиперболическом случае каждая из этих одиннадцати систем координат
в какой-либо заданной части пространства может иметь один, два или три аналога. (При этом, правда, может ока заться необходимым переход из одной системы координат в другую при пересечении характеристической поверхно сти, такой как конус Маха, проходящий через начало коор динат.)
Исследование, проведенное в направлении, указанном Хейсом, показало, что гиперболическое уравнение Гельм гольца должно допускать разделение переменных в двад цати четырех системах координат ^), при этом каждая группа родственных систем, необходимая для того, чтобы покрыть все пространство, считается за одну систему координат.
Говоря о том, что гиперболическое волновое уравнение разделяется лишь в двадцати четырех системах координат, следует пояснить, что речь идет только о преобразовании геометрических пространственных переменных {х, у, г). В противном случае, вследствие инвариантности уравнения (3.4.4) к взаимным заменам t' с у, или z, можно было бы получить ряд новых систем, включающих i'. Более того, если (3.4.4) записать в виде
ф/Ч' + ф1Л/+ фгг = ф*'д:'. |
(3.8.4) |
ТО К координатам Г , у, z применимы классические преобра зования волнового уравнения в евклидовом пространстве.
§ 3.9. Преобразование Фурье
Введем преобразование Фурье по аналогии с § 2.6:
ф |
= |
J |
(3.9.1а) |
Ф = |
-2д |
\ e^^^'Фdx, |
(3.9. Ib) |
1) Одиннадцать нз них были найдены автором |
(Майлс Г “ ]), |
однако его утверждение о том, что гиперболическое уравнение Гельм гольца допускает разделение переменных лишь в этих одиннадцати системах координат, было ошибочным.
где, как и раньше, х может быть комплексным, а путь инте грирования в плоскости X выбран соответствующим обра зом. Д алее мы везде будем пользоваться соответствующими заглавными буквами для обозначения преобразования Фурье.
Преобразуя гиперболическое волновое уравнение (3.4,4), получим уравнение
Д'ф + х2ф = 0, |
(3.9.2) |
которое мы назовем гиперболическим уравнением |
Гельм- |
гольца, по аналогии с классическим уравнением Гельмголь ца. Уравнение (3.9.2) аналогично уравнению (2.6.2) с той лищь разницей, что в последнем вместо х стоит параметр А.= хМ, вследствие выбора в качестве характерного вре мени величины IjU, а не 1/а^. Характеристиками уравнения
(3.9.2) |
являются конусы Маха (3.4.13). |
|
|||||
Применяя преобразование Фурье к выражению для рас* |
|||||||
пределения |
подъемной |
силы |
и к |
граничным |
условиям |
||
(3.4.5—8), |
получим |
|
|
|
|
|
|
|
|
L = 4р-1 (Ф:,, + /М-^кФ),=о*. |
(3.9.3) |
||||
|
|
OA^=O = V, |
' |
(x',y)^ S \ |
(3.9.4) |
||
|
(Ф.^ + /М-^нФ)г=о = |
О, |
(х', у) е W', |
(3.9.5) |
|||
|
|
Ф Iz=O = |
О, |
|
( х ' y)^R ' |
(3.9.6) |
|
|
|
§ 3.10. Гармоническое движение |
|
||||
Как и в случае дозвукового движения (§2.7), рассмотрим |
|||||||
колеблющееся крыло, для |
которого |
|
|
||||
|
|
о = Ке{оехр(//гМ0}. |
(3.10.1) |
||||
|
|
Ф = Ке{фехр (1‘Ш О }. |
(3.10.2) |
||||
|
|
|
|
|
|
|
(3.10.3) |
Здесь |
V и ф — комплексные (не зависящие от t) амплитуды, |
k — приведенная частота и ш — истинная угловая частота.
Для того чтобы связать о и (р с рассмотренными в пре дыдущем разделе преобразованиями Фурье, выразим, пользуясь соотношением (3.4.IOb), t через /' и х' Тогда
Ш/ = Ш р - '( М х '- 0 - |
(3.10.4) |
Для того чтобы получить выражение экспоненциальной зависимости от времени — ехр (/я/'), приравняем коэффи циент при t' величине я
(3.10.5)
Затем введем видоизмененные амплитуды с помощью соотиощений
V = Uexp (/яМх'') = оехр(^^^^^^ |
(3.10.6) |
|
Ф = ^ ехр ( — Ы Ш ') = Ф ехр |
(3 |
.10.7) |
(вспоминая, что, согласно (3.4.IOb), x — ftx'). |
Теперь |
гра |
ничная задача определения Ф по заданным значениям V формулируется соотношениями (3.9.2—6). Зная решение этой задачи, комплексное значение амплитуды давления получим (см. 3.4.5 и 3.9.3) в виде
T = L ехр— ikfAH |
|
(3.10.8а) |
= 4P■^[Ф.. |
‘ |
(3.10.8b) |
§ 3.11. Преобразование Лапласа
Если начало системы координат {х, у, г) или {х\ у, г) расположено в наиболее передней точке плоскости кры ла S , то, очевидно, при х < .0 возмущения отсутствуют вследствие того, что скорость движения превышает скорость звука. В соответствии с этим мы можем поставить условие
Ф =*0 при X (или х ' Х 0. |
(3.11.1) |
Это условие естественно наводит на мысль о применении преобразования Лапласа по координате х\ которое можно применить как к исходному гиперболическому волновому уравнению (3.4.4), так и к гиперболическому уравнению Гельмгольца (3.9.2). Более целесообразно сделать это по отношению к последнему, поэтому применим преобразо вание Лапласа
Ф * = J |
g-SA.-' ф |
= |
(3.11.2а) |
= ^ |
d x ' ^ |
|
( 3 .1 1 .2 b ) |
Соответствующие обратные преобразования имеют
c+ioo
(3.11.3а)
ф = е т 5 |
J |
(3.11.3b) |
Преобразуя выражения (3.9.2 и 3), получим^)
Фуу + Ф*г - |
= 0, |
(3.11.4) |
L* = 4р-1 (S + tM -'x) Ф* |.=о*. |
(3.11.5) |
|
X = (s2 + |
x2)2. |
(3.11.6) |
|
Выбор разрезов для Xв плоскостях s и и будет определяться условиями на бесконечности (§ 1.2), согласно которым, в общем случае, действительная часть X должна быть поло жительна (см., например, (5.2.6)). В связи с граничными условиями на бесконечности следует отметить, что преобра зованное дифференциальное уравнение (3.11.4) является эллиптическим.
^) Более |
подробное рассмотрение, допускающее наличие |
разры |
|
вов.в |
см. |
в работах франкеля .[<>4] и Хислета и Ломэкса |
[®*]. |
Поскольку |
различные граничные условия для ср и Ф |
(в областях 5, |
W и R) являются функциями, зависящими |
от X , для произвольной формы крыла в плане 5 не пред ставляется возможным записать в явном виде граничные условия для функции Ф*”. В отдельных же частных случаях, таких, как, скажем, прямоугольное крыло (глава 7), это может быть сделано. Даже в тех случах, когда не удается выразить в явном виде смешанные граничные условия для Ф*, использование преобразования Лапласа может ока заться целесообразным в качестве промежуточного шага в анализе конкретных задач.
ГЛАВА 4
СВЕДЕНИЕ К ЗАДАЧЕ О СТАЦИОНАРНОМ ТЕЧЕНИИ
§ 4.1 .. Введение
Существует ряд способов, с помощью которых гранич ная задача, связанная с нестационарным обтеканием крыла заданной формы в плане сверхзвуковым потоком, может быть сведена к одной или ряду эквивалентных стационар ных задач. Для некоторых крыльев, таких как прямоуголь ное (§ 7 .3 и последующие), такой путь дает существенные преимущества, но во многих случаях общая задача о ста ционарном обтекании ничуть не проще соответствующей нестационарной задачи, и поэтому существующая возмож ность свести задачу к стационарной не имеет особого прак тического значения.
§ 4.2. Метод Магнарадзе
Как показал Магнарадзе (см. работу Галина [70]), если ¥ (дс', у, г) есть решение уравнения установившегося тече
ния |
|
|
|
|
Л ''Г = |
0, |
(4.2.1) |
то решение гиперболического уравнения Гельмгольца |
|||
|
А 'Ф +х2ф =0 |
(4.2.2) |
|
можно представить |
в виде |
|
|
Ф (^'. У, 2) = '^{х\ |
у, г) — |
|
|
- |
у, Z) I |
у,[к |
- V f \d i (4.2.3а) |
1I)(.V', », г) = j Л (и (л :'» - |
» . Z) d i. (4.2.3b) |
О |
|
где в выражении (4.2.ЗЬ) функция Ф предполагается непре рывно стремящейся к нулю с приближением к передней границе зоны влияния крыла (эта граница лежит в области х' > 0 ) . Магиарадзе получил выражения (4.2.3), сравнивая решения, представленные в виде рядов Фурье, но они могут быть получены и с помощью преобразования Лапласа ^), которое, одновременно, дает возможность связать между собой граничные условия для 'Г с граничными условиями для Ф (см. Майлс 1^"®]).
Поскольку для крыльев простой формы в плане общее решение может быть получено в явном виде (глава 6), основной проблемой, для которой преобразование Магнарадзе дает существенные результаты, является задача о прямоугольном крыле (ср. § 7.4). Эта же идея может быть обобщена путем введения недекартовой системы коорди нат, как это было сделано Стюартсоном [-'®] для гиперболаческих координат (гиперболический аналог сферической полярной системы координат).
§ 4.3. Метод Гарднера
Иной метод был предложен Гарднером [’П, который свел задачу о нестационарном обтекании к двум стационарным задачам. Этот метод, в принципе, применим к крыльям со сверхзвуковыми кромками ‘) произвольной формы, а также к крыльям C дозвуковыми не взаимодействующими кром ками, составленными из отрезков прямой. Однако практи ческое применение этого метода ограничивается случаем
прямоугольного крыла |
(Гудмен [®П, |
[®®], Ломэкс, |
Фуллер |
|||
|
1) Преобразования Лапласа от уравнений |
(4 .2.1,2), |
которым |
|||
удовлетворяют функции |
и Ф**, отличаются |
лишь заменой 5 на |
||||
|
|
Таким образом, выражение (4.2.3) |
следует из обш,его |
|||
соотйошепия, связывающего |
и |
|
(см. Маг |
|||
нус |
II |
Оберхеттиигер [i6°], |
стр. 123). |
|
|
|
|
*) |
См. замечание о |
сверхзвуковых |
дозвуковых |
|
|
на |
стр. 66. |
|
|
|
|
и Сладер В дальнейшем мы ограничимся формулиро ванием метода Гарднера для прямоугольного крыла, распо
ложенного в области X > 0 |
и у > 0 . |
|
|
||
Гарднер |
вводит |
два |
вспомогательных |
потенциала |
|
^ (I. у, 2; X , |
t') и X |
{х , I, t'\ у), связанных |
соотношениями |
||
|
¥,|.=o^ = X (i/> 0 ) |
|
(4.3.1) |
||
Ф (х', у, |
Z, t') = |
'F (О, у, z; |
t'). |
(4.3.2) |
Функция X должна быть решением граничной задачи, формулируемой следующим образом:
|
= 0 |
(g < 0 ). |
(4.3.3) |
X l ^ = O = - ( J t ' > 0 ) , |
(4.3.4) |
||
Xl^=O = O |
( х '< 0 ) . |
(4.3.5) |
В эту задачу у входит лишь как параметр.
Отметим также, что граничные условия для функции X не являются смешанными, поэтому, в принципе, ее опреде ление не вызывает затруднений.
После того как функция |
X |
найдена, граничная задача |
||
для Y формулируется следующим |
образом: |
|
||
|
|
|
= |
(4.3.6) |
’FUo = X(x',E, |
Г; |
г/) |
(у>0). |
(4.3.7) |
IFUo = O |
( у < 0 ) |
(4.3.8) |
||
4^ = 0 |
( K - X ' ) . |
(4.3.9) |
Здесь х' и t' выступают лишь в качестве параметров. Таким образом, определение 'F свелось к решению экви
валентной задачи стационарного обтекания (потоком, направленным вдоль оси I) прямоугольного крыла с перед ней кромкой —х'. Гарднер дает явное решение задачи для W в функции V^).
Трудность обобщения метода Гарднера на крылья не прямоугольной формы в плане заключается в том, что гра ничные условия для функции X (4.3.4) и (4.3.5) не удается
1) Если зависимость от времени имеет вид ехр (Ш ’), решение Гарднера сводится к решению, полученному в § 4.2.
выразить в явном виде, если уравнение боковой кромки содержит X Это означает, что на плоскости ^ = O будет существовать область, где величина X останется неизвест ной. Естественно, в тех случаях, когда граничная задача, связанная с крылом не прямоугольной формы, может быть представлена в виде, эквивалентном случаю прямоуголь ного крыла (см. главу 8), описанный метод применим и к случаям крыльев другой формы в плане. Так, Гарднер высказал предположение, что метод может быть обобщен на случай любого крыла, к которому, в стационарной за даче, применим метод Эварда (см. § 7.3) и при условии, что все дозвуковые кромки образованы ломаными прямыми линиями. Однако применение метода к крыльям с числом углов в плане более четырех приводит к чрезмерно громозд ким вычислениям, и даже для этой формы метод не сулит очевидных преимуществ перед другими.
Гарднером p -j был предложен аналогичный метод и для конических тел. Mo область применимости его весьма ограггичена и практического развития и прилолсений он не получил (см. также Жермен и Бадёр [’®J).j
§4.4. Метод разложения по частоте
Втех случаях, когда характер движения крыла изме няется медленно (например, гармонические колебания низ кой частоты), целесообразно пытаться рещать задачу методом преобразования Фурье, представляя рещение в виде разло- л<ения в ряд по параметру частоты х. Рассмотрим вначале крыло, задние кромки которого полностью сверхзвуковые,
апотому в граничные условия параметр х явно не входит. Тогда подлежащая решению граничная задача формули руется системой соотношений (3.9.2, 4 и 6), а именно:
А'ф-Ьх2ф = 0, |
(4.4.1) |
|
Ф, !.= O = - V , |
(x ',t/ )6 S ' |
(4.4.2) |
ф|,=о = 0 |
(х', t/)6/?'. |
(4.4.3) |
Граничное условие (3.9.5) здесь молшо не учитывать, поскольку течение в следе не влияет на обтекание самого крыла, вследствие предположения о том, что задние кром к и — сверхзвуковые.