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

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

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

Таблица I. Ограничения на параметры для различных случаев упрощения задачи обтекания тонкого крыла!)

Характеристика случаев

Дифференциальное уравиеинс потенциала скорости

Л 1

Нелинейное, околозвуко­ [ (М = - 1) - | - ( ^ ) Ф ,] Ф^ ^ + 2 ( ^ ) Ф:с1-Ч > ^ -'Р :^ = < '

 

вое течение

Az Д озвуковое и сверхзву­ ковое течение

-As Нестацнонарное около­ звуковое течение

A i Крыло малого удлинения или квазистационарное околозвуковое течение

B i Высокие частоты колеба инй, сжимаемая среда

B i Высокие частоты колеба­ ний, крыло малого удли­ нения

( М * - ' ) Ф х * + 2 ( ^ ) , - Ь ^ Ф „ - Ф „ « - Ф « = 0

о* '’’(I Ф х* Фу1/ Фг2 ®

1:Г ^ и -Ч > уу-Ч > г2= 0

B i

Высокие частоты

колеба­

 

 

ний, крыло квазнмалого

 

 

удлинения

 

 

B i

Теория поршня

 

- ^ Ф , , - Ф „ = 0

Bz

Высокие частоты колеба*

Ф хх+ Ф «у+Ф 2г= 0

 

инй, несжимаемая

среда

 

 

B s

Высокие частоты

колеба­

Ф уу+Ф 22=0

 

ний, крыло очень малого

 

 

удлинения

Давление на крыле

"Pz2= 0-f

~ (^ Q (f^ )

Соотиошенне подобия для ( ^

(УК.

б“/я

36Vox,

 

( M - I )

(Y -I-I)V s

(Y-t-l)^/36Vs’ (Y-|-1)®/S6V3

 

 

 

 

 

---------^

F [1 М2-

1

I V s^,

feM2

 

IМ2- I I

 

IM2- 11V2 L

 

 

 

6/*'(X. Л)

6 kF (k)

A 26F(X, Ш )

m l F ikM l)

co n st /гбМ I

ф, <

Л2б Г (М

con st

 

 

 

Ограничения на параметры

 

 

 

 

 

 

Удлинение Я-

 

Число М аха

Приведенная частота

Относительная

 

толщина илБ

(1 »

6)

 

M

 

к

 

амплитуда

 

 

 

 

 

 

 

6

 

Ti-I=O

(6V 3)

 

I M - I 1=0(б^/з)

k = 0 (б*/з)

 

6 <SI

 

X -I= O d М 2-1 |V2)

 

I M - I 1 » б^/з,М «б-1

 

-A=OCI 1 -М 2 |)

6 «

I,

6 «

M -I

X -I= O (/«V2)

 

I M - I

I « Л

й = 0 (1 )

; А » б “/з

6 «

I .

6 «

fe^/2

 

_

2

 

 

 

U » (МЯ.)-2

 

 

 

 

X « 6 - V 3 , Х<'<1 М 2 -1 I

2

I М 2 -1

К<Х -2

 

 

6 « I

 

X -I= O (ЛМ)

 

М -1|= 0(Л ), М «(Лб)-1

 

■ ft«(M 6) - i

6«(ftM-i, 6«ft-l)

1 «

I

 

M=O (ft-iX-l)

 

ft «

6 -1

 

 

 

 

X -I= O

(кЬЛ)

 

M »

-1

 

ft» M-I

 

 

 

 

 

 

 

 

1 ^

ft ^ 6-1

2

 

 

 

 

 

 

 

 

 

 

 

 

 

^

6 « 4

 

X -I « Ш

 

M » fe-i

 

ft » M-t

 

 

 

 

X -I= O (1)

 

M « !fe -l

 

ft «

M -I

 

 

 

 

Х «

1

 

M <<(/гХ)-1

 

ft «

(M X)-I

 

 

 

 

I) Координаты H потенциал скорости в этой таблице н с б е з р а з м е р н ы , однако теми ж е символами в последующих главах обозначены беа)рйзмернЫе величины (см . главы 2 и 3 ).

получим

фГ =o± = ± - у Г ‘

(^. '/)6-5,

(1.2.35b)

где через S обозначена проекция крыла на плоскость 2 = 0. Теперь следует установить граничные условия в осталь­ ной части плоскости z = 0. Вообще говоря, необходимо

Рис. 1.1. Проекция

поверхности крыла

(S)

на плоскость z =

0,

вихревая пелена

(1¥^,

ос­

тальная

часть плоскости (Р).

 

 

делать различие между вихревой пеленой, ограниченной задней кромкой области 5 , и двумя линиями, сходящими C концов крыла, параллельными направлению полета и про­

стирающимися

вниз по потоку

(область W на рис. 1.1),

C одной стороны, и остальной

частью плоскости {R на

рис. 1. 1) — C

другой.

 

Всимметричной задаче мы непосредственно обратимся

кфизическому условию непрерывности вертикальной составляющей скорости (антисимметричной относительно

плоскости Z = O), из которого следует,

Ф1*^г=0 = 0 .

{X, у) ^ R + W

(1.2.36)

Характерная простота (в принципе) симметричной задачи

является прямым следствием

того,

что па всей

плоскости

2 = 0 известны значения ср.,

т. е.

граничные

условия не

являются смешанными (в отличие от антисимметричной задачи, см. ниже).

Симметричные граничные условия в нестационарном потоке реализуются на практике в случае, когда в нежест­ ком крыле возникают симметричные, изменяющиеся во времени деформации поверхностей, или когда мы рас­ сматриваем задачу о coпpotивлeиии, связанном с ускоре­ нием в направлении движения. Метод решения будет анало­ гичен изложенному в § 6.2 этой книги, однако первая из названных задач, видимо, не представляет практиче­ ского интереса, тогда как во второйобычно достаточно рассматривать течение в каждый момент времени так же, как если бы мгновенная скорость полета была постоянной (т. е. предполагая, что изменение скорости полета при прохож­ дении пути протяженностью в одну хорду мало в сравнении C мгновенной скоростью полета; ср. с § 1.5). Мы подробно рассмотрим лишь движения, поперечные к направлению полета и идентичные на верхней и на нижней поверхностях крыла. Таким движениям будут соответствовать антисим­

метричные граничные условия.

 

В антисимметричной задаче

автоматически полу­

чается непрерывной при z = 0, и необходимо найти другие граничные условия. Мы можем поставить условие, что в области R непрерывны тангенциальные составляющие скорости, откуда получим

ф^“Мг=о = 0,

{х, y)^ R .

(1.2.37)

C другой стороны, касательные скорости могут иметь разрыв при переходе через вихревую пелену, и мы обязаны потребовать, чтобы давления здесь оставались непрерыв­ ными. Воспользовавшись соотношением (1.1.13), получим

(1.2.38)

Интегрирование этого выражения в системе координат, жестко связанной с поверхностью S, даст нам

i = //), { x , y ) O V , (1.2.39)

где 2/ — разрыв потенциала ф, движущийся вниз по потоку со скоростью полета U (и потому стационарный в системе координат, неподвижной в пространстве). Подчеркнем, что если условие (1.2.38) распространяется наряду с открытой областью W также и на заднюю кромку (разделяющую

области S и

U/), то давление на задней кромке остается

непрерывным,

и условие К утта— Жуковского — выпол­

ненным.

 

Остается установить граничные условия на бескоиечFIOCTH. В задаче с начальными данными требуется лишь чтобы возмущения оставались ограниченными на больших расстояниях от их источника (условие !юнечности возму­ щений). Мы будем рассматривать целый ряд задач о квазйустановившихся движениях тела (с гармоническим законом изменения параметров по времени), и тогда необходимо (ср. Бекер и Копсон [ “], стр. 25) добавить требование, чтобы в системе координат, по отношению к которой воздух поко­ ится, возмущения на бес1сонечности представлялись в виде уходящей волны [условие радиации). Иными словами, можно считать, что частота имеет небольшую отрицательную мни­ мую часть и наложено лишь условие конечности возмуще­ ний (Бекер и Копсон I®], стр. 154, ср. с §2 .6, ниже в этой книге). Если система координат движется со сверхзву­ ковой скоростью, условие радиации следует заменить ус­ ловием отсутствия распространения возмущений вперед по потоку.

В загслючение заметим, что если в окрестности тела линеа­ ризация допустима, то задача о несущих свойсхвах тела может быть отделена от симметричной задачи, даже если последняя приводит к нелинейным уравнениям в некоторых других областях поля потока. В соответствии с этим усло­ вия, рассматриваемые в этом разделе как достаточные для того, чтобы провести линеаризацию в окрестности тела, можно считать также достаточными и во всем поле потока,

если нашей конечной целью является получение равномерной аппроксимации для поперечных сил, действующих на тело.

§ 1.3. крылья малого удлинения

Граничную задачу для крыла можно существенно упро­ стить, если удлинение крыла достаточно мало, чтобы считать

член C г|)т)11 главным в

правой части

уравнения (1.2.8). Эта

возможность

соответствует случаю

(1.2.20), если ^ = O (I),

или случаю

(1.2.23),

если А » 1 .

 

Сначала рассмотрим ограничения, связанные со случаем At (см. (1.2.20)). В дополнение к ограничениям, налагаемым

условиями

(1.2.32), получим

 

е = о б « 1 ,

о З б с 1 , I M ' ^ - l l o ' » ^ ! ,

(1.3.1)

 

Vvv + Vzz

(1.3.2)

Уравнение для давления (1.1.13) остается тем же. Ограничениям (1.3.1) можно удовлетворить различными

путями. Так, если |M — 1 1и ^ достаточно малы, то требует­ ся лишь, чтобы величина а была малой по сравнению

е й 3 (а не C единицей), а течение является мазистацштар' ным, трансзвуковым. C другой стороны, если а < 1, M и могут иметь средние значения и тем не менее продолжать удовлетворять условиям (1.3.1). Следует также напомнить, что для того, чтобы задача оставалась плоской (см. ниже), условие а > б должно выполняться даже если а < 1.

Сведение уравнения для потенциала к двумерному урав­ нению Лапласа, в котором координата, направленная по потоку (л:), и время (i) входят лишь как параметры, позво­ ляет применить такой мощный аппарат, каким является конформное отображение, и существенно расширить диапа­ зон форм крыла в плане, которые могут быть исследованы. Граничные условия остаются теми же, что и в § 1.2, только геометрическая форма границ сводится к одной или несколь­ ким прямым линиям в плоскости 2 = 0.

Аппроксимация, называемая нами крылом малого удли­ нения, вначале была получена Джонсом 1^°®] из рассмотре­ ния изменения количества движения воздуха, сообщаемого ему поперечным сечением крыла, взятым вдоль его размаха, подобно тому, как это было сделано Мунком [^^’1 в его рабо­ те об обтекании корпуса дирижабля. Тот же метод приме-

ним также и к задачам о нестационарных течениях (Мдйлс [108]^ [18G]) 5олее строгое изложение теории обтекания как крыльев, так и не плоских тонких тел в установившемся

сверхзвуковом потоке было дано Уордом

а

затем

обобщено

на нестационарные

течения

Майлсом

 

[100J

[107]^

 

 

 

 

Если условие того, что тело

плоское,

а > б,

заменить

условием

пространственности, а = О(б),

граничные

усло­

вия следует задать при некотором среднем положении поверхности тела, поскольку решение уравнения (1.3.2) в этом случае будет иметь особенность вдоль оси тела. Более того, при этом обычно необходимо учесть члены с фу и ф| в выражении для давления (в плоской задаче в этом не было нужды благодаря условию а > б, хотя Спрейтер 1-®®1, не обратив внимания на это ограничение, высказал противоположное мнение); ср. § 12.1, в частности уравне­ ние (12.17) и последующие.

Если А > 1, как это часто бывает в практически интерес­ ных задачах о нестационарных течениях, то рассматривае­ мая аппроксимация для крыла малого удлинения должна основываться на уравнении для потенциала в форме (1.2.8Ь),

а не (1.2,8а), Простейшим

является

случай

для кото­

рого, согласно условиям

(1.2,23),

получаем

ограничения:

= k a b < 1,

I ,

Ш а < 1 ,

(1.3.3)

при этом уравнение для потенциала сохраняет форму (1.3.2). Еще более упростится выражение для давления BTOM отнощении, что членом ф,. можно пренебречь по сравне­ нию C ф^, но это не приводит к заметным упрощениям задачи.

Ограничение А Mo < 1 является весьма сильным, поэтому представляется интересным рассмотреть такие варианты случаев или B j, где членом фц в уравнении (1.2.8Ь) можно пренебречь [полное уравнение (1.2.8Ь) в случае B j ничуть

не проще уравнения (1.1.12)]. Рассмотрим два

возможных

случая,^) а

именно:

 

 

 

 

в;':

Cs = O(I),

C ,=

!;

(1.3.4)

Cfi= 1,

с, =

0

(1 ).

(1.3.5)

1) Третий случай рассмотрен в § 1.4.

Отсюда в дополнение к ограничениям (1.2.32) получим соот­ ветствующие условия:

В [:

е = / габ < 1,

а < 1 ,

Шсг = 0 ( 1 ); (1.3.6)

В [ :

б

1.

ш:

 

 

 

В обоих случаях уравнение для потенциала сводится к виду

Ф„„ + Ф.з = -^Ф/|-

(1.3.8)

Это двумерное волновое уравнение, в котором направлен­ ная по потоку координата (jc) является лишь параметром.

Выражение для давления

аналогично выражению (1.1.13),

хотя,

как и в

случае

в нем можно пренебречь членом

C

благодаря

условию

А > 1.

Из предыдущего очевидно, что если оба условия — <т<1 и М(У < 1 — удовлетворяются, то уравнение (1.3.8) справед­ ливо при любых А, хотя при A = O (I) оно и сводится к урав­ нению (1.3.2). При больших А, естественно, остаются усло­ вия А6 < 1 и АМб < 1, однако если величина Аб равномерно ограничена (ср. примечание на стр. 19), то эти условия обычно не являются слишком стесняющими.

В заключение заметим, что теория крыла малого удли­ нения может быть использована в качестве отправной точки при построении решений, имеющих вид разложения по степеням а, как это было сделано в посвященной стационар­ ным течениям работе Адамса и Сирса [“Ч, которая легко обобщается на нестационарный случай (§9 .8).

Дальнейшее развитие методов теории крыла малого удлинения будет рассмотрено в главе 9.

§ 1.4. Теория поршня

Если предположить, что произведение Ш а в случае В\ достаточно велико, то в дифференциальном уравнении для потенциала член с может быть опущен, и тогда мы полу­

чим новый частный случай В", когда

е =

А М » 1 .

А М о » 1 ,

(1.4.1)

 

аоФ** = ф».

 

(1.4.2)

 

 

 

Это приближение можно назвать (линеаризованной) теорией поршня, поскольку возмущенное давление в точке на поверхности тела, где скорость, направленная вдоль внешней нормали, равна v, определяется соотношением (ср. C уравнением (9.1.13))

P - Po = - Qo^oU,

(1.4.3)

т. е. в соответствии с классическим результатом для давле­ ния на поршне в одномерном движении, или, что эквива­ лентно, C соотношением между давлением и скоростью

плоской звуковой волне.

Первое применение теории поршня к задачам о крыле, видимо, принадлежит Хислету и Ломэксу [^*-], а возмож­ ность ее использования в задачах флаттера широко исследо­

вана в работах Молло —Христеисеиа и Эшли

а также

Эшли и Зартаряна

р]. В

1953 г. Лайтхилл

развил

нелинейную форму

этой

аппроксимации применительно

к гиперзвуковым задачам. Этот последний метод будет рас­

смотрен в § 13.4,

 

 

§

1.5. Ускоренный

полет

В предыдущих

разделах речь

шла преимущественно

о поперечных по отношению к равномерному установивше­ муся полету возмущениях.

Ко второй категории задач о нестационарных течениях, которую мы рассмотрим здесь лишь очень бегло, относятся задачи о движении с ускорением, направленным вдоль тра­

ектории

полета.

 

 

 

 

 

Характерным

безразмерным

параметром

ускоренного

полета

является

величина

 

 

 

 

 

 

 

и Л

 

 

 

 

M

 

 

 

(1.5.1)

где I — длина тела (в

направлении полета). Очевидно, что

M в большинстве

практически

возможных случаев

мало,

а следовательно, как показано в работе Майлса

влия­

ние ускорения может

быть существенным

лишь

когда

к трансзвуковой параметр

[1 = \2М\ ' ^ ( M - I )

(1.5.2)

Переходя к изложению второй и третьей глав, отметим, что безразмерное время в ускоренном полете удобнее опре­ делять как чем как Ut/l, где U является величиной, зависящей от времени. Кроме того, зависимость U от време­ ни неудобна при переходе в уравнениях движения к системе координат, движущейся вместе с телом. Так, в частности, потенциал скоростей ускоренного полета удовлетворяет уравнению (3.2.2) и не удовлетворяет (3.3.3).

ПРЕОБРАЗОВАНИЕ УРАВНЕНИЙ ДОЗВУКОВОГО

ОБТЕКАНИЯ КРЫЛА

§2.1. Введение

Вэтой главе мы рассмотрим различные преобразования линеаризованного уравнения для потенциала дозвукового

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

§ 2.2. Введение безразмерных параметров

Чтобы сделать безразмерными наши обозначения, отне­ сем все длины, скорости и давления к следующим характер­

ным величинам: I — длина, till — время и у “ ско­

ростной напор.

Величину I можно выбрать произвольно, но в большин­ стве случаев в качестве I выбирается величина хорды или размаха крыла. U, как и в главе 1, — скорость полета.

Характерным временем в этой системе является время, потребное крылу для прохождения пути, равного характер­ ной длине, но C не меньшим основанием в. качестве характер­ ного времени можно принять величину //По* т. е. время рас­

пространения

малого возмущения

на характерную дли­

ну. Основным

возражением

против

последнего

варианта

является неопределенность,

возникающая при

переходе

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