книги / Потенциальная теория неустановившихся сверхзвуковых течений
..pdfК случаю несжимаемого потока, когда |
полагается беско |
нечной величиной. |
|
Обозначим безразмерные декартовы |
координаты в непо |
движной системе (рис. 2.1) через X , у, |
2, где X измеряется |
вниз по потоку от самой передней точки крыла в начальный момент времени, у положителен в сторону правого крыла,
©
(T)
Рис. 2 .1. Система координат, неподвижная в про странстве.
а положительным направлением оси z принято направле ние вверх. Соответствующее безразмерное время в этой системе обозначим через Т. Теперь, если безразмерны" потенциал скоростей ф таков, что вектор возмущенной скорости q задается соотношением
д = Ш ц>, |
( 2.2. 1) |
то волновое уравнение (1.1.10) приобретает вид
Дф = M V т . |
(2-2.2) |
где оператор А ( = V^) записан в координатах X, у, z. Иначе уравнение (2.2.2) можно переписать в виде
Граничные условия иа крыле получим, приводя к безраз мерному виду выражение (1.2.35Ь)
ф2 Iz=O = ■ дТ = — V, |
(2.2.4) |
где Vопределена как заданная величина скоса потока, отне сенная к и . (Следует помнить, что v не является компо нентой скорости потока.)
Поскольку нас интересует лишь давление на поверх ности крыла, введем безразмерную величину перепада дав лений (которая будет характеризовать местную (удельную)
подъемную силу) |
|
I-=O-- ' ’ I-» + |
(2.2.5) |
Подставляя р из (1.1.11) и пользуясь антисимметрией |
|
течения, получим^) |
|
/=4срт Iz=O- |
(2.2.6) |
§2.3. Преобразование Галилея
Всвязи C трудностями, возникающими при формулиро вании граничных условий на движущихся телах, представ ляется удобным ввести систему координат, движущуюся вместе с поверхностью 5 . Проще всего сделать это с помощью преобразования Галилея
(2.3.1а)
(2.3.1b)
для которого якобиан
‘{Vr) |
(2.3.2) |
|
I) Так как при безразмерном рассмотрении характерная длина никуда явно не входит, то смысл обозначения / всегда будет ясен из контекста.
Координата Jt измеряется теперь от самой передней точки по верхности S (рис. 2.2). (Преобразование записано в матрич ной форме, для того чтобы подчеркнуть свойства симмет рии и упростить сравнение с преобразованием Лоренца, рассмотренным в § 2.4.) Обычно в практических задачах величина v задана в функции х и t, а потому и окончатель ное решение для ф выразится через те же переменные.
Вводя преобразование (2.3.1) в уравнение (2.2.2), полу чим (ср. C (1.1.12)):
= |
+ |
(2.3.3) |
или, в более распространенной форме,
(1 - M"=) Ф*.., + ф„„ + ф „ - 2 M 4 i - М^Фп = 0. (2.3.4)
Соответствующее выражение для перепада давлений будет иметь вид (ср. C (1.1.13)):
Граничные условия в координатах {х,у, г, i) запишутся сле дующим образом:
+г‘“>, (X. y ) ( S . (1 3 .6 )
(л:, y )iW
(х. y ) iR ,
где границы областей S, W и R не зависят больше от Главным достоинством преобразования Галилея является
то, что оно связывает координаты с поверхностью S и оставляет время инвариантным. Основной недостаток — в том, что преобразование не оставляет волновое уравнение инвариантным.
§ 2.4. Преобразование Лоренца
Для того чтобы оставить инвариантным волновое урав нение и одновременно связать систему координат с поверх ностью S, обратимся к преобразованию Лоренца (Страттон [280]^ стр. 77; Курант и Гильберт I"*’ !, т. 2, стр. 429).
' 6 9 -
(2.4.1а)
(2.4. Ib)
(2.4.1с)
(2.4.2)
Уравнения (2.3.3—8) в преобразованном виде запишутся следующим образом:
Аф =M^cprr. |
(2.4.3) |
|
Фх'х' + ф„„+ф=2- М®фгг= О, |
(2.4.4) |
|
/ = 4р-1(Фх' + ф/')г=о+, |
(2.4.5) |
|
ф х М = — |
{х \ y ) ^ S \ |
(2.4.6) |
(фхЧ-фг)2=о = 0, |
(х’ , у )^W ', |
(2.4.7) |
Ф М = 0. {х \ у ) ^ R' |
(2.4.8) |
Подчеркнем, что при переходе от (S, U^, R) к (5 ', W , R') происходит изменение масштаба, как это показано на рис. 2.3. Поэтому удлинение преобразованного крыла (часто именуемое эффективным удлинением) связано с удлинением исходного крыла следующим соотношением ^):
V = рХ. |
(2.4.9) |
Преобразование Лоренца было введено в дозвуковую теорию крыла Кюсснером 1^^®]. Оно сводится к комбинации
(to
Р и с.-2 .3 . Система координат, связанная с крылом, после преобразования Лоренца.
преобразований Галилея и Прандтля — Глауэрта (послед нее — просто изменение масштаба) в установившемся тече нии.- Сравнивая преобразование Лоренца с преобразова нием Галилея, видим, чтр свойство инвариантности времени здесь сменилось инвариантностью волнового уравнения,
1) _3десь мы используем стандартное определение величины удли
нения Я как квадрата |
|
максимального размаха крыла, |
поделенного |
|
на его площадь. Более |
общим |
является соотношение |
о '~ ^ а , где |
|
определение величины |
а см. в § |
1.2. |
|
а это, в свою очередь, привело к тому, что t' уже не явл51ется истинным временем.
Между координатами (2.3.7) и (2.4.1) существует очевид ная связь
Р ) =
f ) =P"(-м м ) (Я
Что касается симметрии матрицы преобразования (2.4.Га)< заметим, что выбор в качестве характерной скорости вели
чины «о Дзл бы нам вместо матрицы C .:) симметричную
(\ |
M |
(ср, |
C (3.4.1)). |
матрицу |
^ ^ |
||
|
§ |
2.5. |
Несжимаемая жидкость |
Классический предельный случай течения несжимаемой жидкости получается, если положить ао стремящимся к бес конечности в линеаризованных соотношениях предыдущих
разделов. |
При этом число |
M становится равным нулю, |
а потому |
преобразования |
Галилея и Лоренца совпадают |
и получающаяся граничная задача отличается от сформу лированной в § 2.3 лишь тем, что уравнение (2.3.3.) сво дится к уравнению Лапласа:
Дф= Фхх+ |
Ф„о + Ф.г = О- |
(2.5.1) |
Следует отметить, что в |
предельном случае |
M -^O не |
обязательно предполагается течение несжимаемой жидкости. Возможной альтернативой является случай (У = 0, как- в акустике, но тогда уравнение (2.3.3), хотя и с очевидной неопределенностью (поскольку t отнесено к UU), сводится к обычному волновому, уравнению. Эти два предельных случая, AQ- > и I/—>0, хотя обычно и отличаются сущест венным образом физически, могут иногда совмещаться, как, например, в случае быстрого ускорения .тела в почти несжи маемой жидкости (см. Майлс Лонгхорн
§ 2.6. Преобразование Фурье
При рассмотрении случаев произвольной зависимости от времени целесообразно ввести преобразование Фурье по координате t' Из многочисленных возможных форм инте грала Фурье мы выберем ту, которая принята в таблице
Кемпбелла и Фостера |
1*“], |
а именно: |
|
Ф |
= ^ |
g -ix rjp c ((' |
(2 .6 .1 а ) |
со.
( 2 .6 . Ib )
-OO
Вобщем случае х может быть комплексным, и путь инте грирования в плоскости X должен выбираться в соответ ствии C определенными физическими условиями.
Преобразуя выражения (2.4.3—8) и обозначая преобразо вания Фурье от различных функций соответствующими заглавными буквами, получим
АФ-ЬА,2ф = 0, |
A- = Mx, |
(2.6.2) |
L = 4p-’^(Ox' + txO)j=o+, |
(2.6.3) |
|
Ozlz=O=-- V , |
|
(2.6.4) |
(0.v + ix0),==0 = 0, |
y ) ^ W \ |
(2.6.5) |
О Iz=O = 0, |
(A:'. y)^R'. |
(2.6.6) |
В некоторых случаях, особенно когда нужно включить начальные условия, более естественным, чем преобразова ние Фурье, является преобразование Лапласа. Хотя между этими двумя преобразованиями нет существенной разницы (поскольку X рассматривается как комплексное), преобразо вание Лапласа позволяет обойти такой парадокс^), как
См. серию |
заметок |
Ван де Воорепа, Лайтона, Джонса, |
а также Денглера |
и др. в J . |
Аегоп. Sci 19, 209 — 14 (1952). |
возникающий при исследовании затухающего движения, когда невозможность введения начальных условий, соот ветствующих конечному движению, подразумевает наличие бесконечного источника энергии при i = — со. Тем не менее нами избрано преобразование Фурье благодаря его глубо кой связи C важнейшим частным случаем нестационарного движения — гармоническими колебаниями.
§ 2.7. Гармоническое движение
Мы рассмотрим здесь случай, когда движение колеблю» щегося крыла задается уравнением
O |
= Re {ое*"'},(2.7.1) |
где |
|
k - — |
(2.7.2) |
V — комплексная амплитуда (не зависящая от /), k — при веденная частота и ы — истинная (а не безразмерная) угло вая частота. Потенциал записывается в форме
(P = Re |
(2.7.3) |
изадача приводится к определению ф в функции от о. Для того чтобы свести эту задачу к задаче, сформули
рованной в предыдущем параграфе для преобразования Фурье, установим, исходя из (2.4.IOb), что
kt = ifep-i { f - M V ) = k^-Ч’ - |
(2.7.4) |
Отсюда, чтобы получить зависимость от времени ехр (tx/') в виде, представленном в уравнении (2.6. Ib), положим
X = P-^A:, |
(2.7.5) |
и введем видоизмененную амплитуду колебаний
T/ |
'ikWx |
(2.7.6) |
к = O ехр — - Р* • |
так, что величины Ф и V будут связаны решением задачи, сформулированной уравнениями (2.6.2 — 6). Таким образом,
зная 1», мы определим V из (2.7.6), решим задачу (2.6.2 — 6)
И найдем ф из соотношения
Ф = Ф е х р ^ ^ |
(2.7.7) |
Подчеркнем, что видоизмененные амплитуды V и (D не являются преобразованиями Фурье функций о и ф в том смысле, в каком об этом говорилось в предыдущем разделе. Поэтому ф и Ф, определенные соотношениями (2.7.3) и (2.7.7), не удовлетворяют равенством (2.6.1).
Вбольшинстве случаев нас интересует не потенциал ф,
аперепад давлений /. Если мы запишем
/ = Re{7e*'‘'} |
(2.7.8) |
и сравним это выражение с (2.6.3) и (2.7.7), то получим
|
|
(2.7.9а) |
|
|
-(Ф .' + 1р-'А:Ф),=о+= |
(2.7 |
.9Ь) |
= 4 ехр |
(Ф.., 4- ф-УЩ г= |
(2.7 |
.9с) |
Приведенные выше соотношения позволяют свести решение задачи о колеблющемся крыле в сжимаемой жид кости (М # 0) к решению аналогичной задачи в несжимаемой жидкости, если членами порядка (ЛМ/Р)“ можно пренебречь, и получить таким образом обобщение хорошо известного для стационарных задач правила Прандтля—Глауэрта на нестационарные течения (Майлс 1^®’^]).
Г Л А В А 3
ПРЕОБРАЗОВАНИЕ УРАВНЕНИЙ СВЕРХЗВУКОВОГО
ОБТЕКАНИЯ КРЫЛА
§ 3.1. Введение
Обратимся теперь к преобразованию линеаризованных уравнений, описывающих обтекание тонких крыльев сверх* звуковым потоком. Построение этой главы, в основном, аналогично построению предыдущей, однако в ней будут рассмотрены некоторые дополнительные виды преобразо ваний. Благодаря ограниченности области влияния любой точки в поле све])хзвуков9го течения, круг возможных задач и методов их рещения, вообще говоря, расширяется. В част ности, сверхзвуковая задача существенно упрощается, если компонента скорости потока, нормальная к задней кромке крыла в плане, повсюду остается сверхзвуковой, что позво ляет, при определении потенциала на крыле, не рассматри вать течения в его следе.
§ 3.2. Введение безразмерных параметров
За характерные величины здесь мы выберем (ср. с § 2.2) длину /, время //Но и давление QOC/“/2.
Отметим, что существует ряд практически важных задач, например, определение нагрузки при входе в порыв ветра крыла бесконечного размаха, или треугольного крыла, где нельзя установить характерную длину, и которые поэтому являются инвариантными к преобразованиям масштаба. В таких случаях естественно искать однородное решение типа конических течений Буземана [®Ч, [^"]^).
1) См. также работы Хейса [в®], Лагерстрома |
Жермена [тв|, |
Гольдштейна и Уорда [в®] и Уорда fs®®). |
|