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

книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения

..pdf
Скачиваний:
8
Добавлен:
24.10.2023
Размер:
9.31 Mб
Скачать

160

ДЖ . ДЖ . СТОКЕР

 

уравнение неразрывности

имеет вид

 

 

 

их + ѵу =

0,

(3.6)

а условие потенциальности требует,

чтобы исчезал вихрь

 

 

ѵх — «г/

=

0.

(3.7)

Уравнения

для верхнего

слоя получаются заменой в написанных

уравнениях

малых букв

заглавными.

На свободной

поверхности

5 (х, у) = 0

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

 

 

 

Usx + Vsv =

0.

(3.8)

Теперь сформулируем остальные условия на свободной поверхно­ сти, дне и границе раздела.

Удобно ввести функцию тока ф (х, у), которая удовлетворяет

уравнениям

(3.9)

ц = фу, о = —ф*.

Линии тока определяются уравнением

 

ф (х, у) = у = const.

(3.10)

Выберем у так, чтобы

 

ф (х, 0) = 0 на дне,

 

{ch на поверхности раздела,

(З.П)

с (h + Я) = ch (1 + г) на свободной поверхности.

Теперь предположим, что все линии тока взаимно однозначно проектируются на ось х. Тогда уравнение (3.10) можно обратить, что приводит к равенствам

У =

f

(х, у) в нижнемслое,

(3.12)

у —

F (х, у) в верхнемслое.

(3.13)

Функции / и F называются

функциями линий тока. Возьмем

х и у

в качестве новых независимых переменных и положим

 

и (х, у) = и (х, / (х, у))

— и (х, у) и т. д. в нижнем слое.

(3.14)

Для верхнего слоя воспользуемся аналогичными обозначениями. Таким образом,

их — их + иуух = их + ыѵфж= их ѵиу и т. д.

(3.15)

Уравнения (3.4) —(3.7) принимают следующий вид:

 

6иих = — рх + ѵру,

(3.16)

öuvx = — 6g — иру,

(3.17)

ux — vuy + uvy= 0,

(3.18)

vx— vvy — ішѵ = 0.

(3.19)

IX . Б И Ф У Р К А Ц И О Н Н Ы Е Я В Л Е Н И Я В ТЕО РИ И П О В Е Р Х Н О С Т Н Ы Х ВОЛН 161

Для верхнего слоя выполняются аналогичные уравнения. Условие на свободной поверхности (3.8) принимает такой вид:

и (sx vsy) -f vusy =5 usx = 0 .

(3.20)

Уравнение (3.20) тождественно выполняется, так как свободная

поверхность задается уравнением s (х, с (h + Н)) = 0. Преимущество новых переменных состоит в том, что уравнения

(3.16)—(3.19) нужно решать в фиксированной полосе: —оо ■< х < оо,

0 < у < ch\

аналогично, уравнение для верхнего слоя нужно решать

при —оо <

X < оо, ch < у < с (h + Н).

В качестве зависимых переменных вместо и, ѵ и р удобно исполь­

зовать и, f (или F) и р.

Тогда

 

Т у —

f x = 1\>x / % = ulu.

(3.21)

(Предполагается, что и ф 0.) Тогда уравнение неразрывности (3.18) удовлетворяется тождественно (и это не удивительно ввиду связи

между / и функцией токаф). Уравнения (3.16) и (3.17) принимают вид

0И*= f y P x + f x P y ,

(3.22)

б { Üf x x + Ü J x ) = — бg f y Ру.

(3.23)

Уравнение неразрывности (3.18) заменяется первым уравнением (3.21):

üfy = 1.

(3.24)

Условие потенциальности (3.19) принимает вид

т - І К Н Ш / Г ѵ І ѵ - О -

(3-25)

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

в = / — gh/c2,

(3.26)

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

1 1 — 0 1 2 8 5

162

Д Ж . Д Ж . СТОКЕР

 

 

менная растянута:

 

 

 

 

<т = ]/"гх/h,

г] = уIch,

и* = и/с,

U* =

U/ C,

V* = ѵ/с,

V* = Ѵ/с,

р* = р/6с2,

Р* =

Р/Ас2, (3.27)

f* = f/h,

F* = F/h,

r — H/h,

р =

Д/6.

Теперь опустим звездочки. Наши определяющие уравнения прини­ мают такой вид:

Wer—

/тіРст + fapr\,

 

 

 

 

 

(3.28)

б (w/ста +

wa / a) = е / л

//n

Р „ ,

 

 

 

и / л = 1 »

 

 

 

 

0 < т] < 1,

 

/іГП ~Ь е [/с/яч “Ьfr\foa 2fofr\foT)] = О,

 

 

 

 

 

U а

F ■цРа “ Ь F o P f\,

 

 

 

 

 

 

B{UFm + UaFc) = BFrllF r i- P ^

 

 

 

(3.29)

UF4= 1,

 

 

 

0 <

т] <

1 Н-г,

 

F т\т)“Ье [ F % F -f-F^Fаа

2Р0Р„Рот)] =

0.

 

 

 

 

Граничные условия имеют вид

 

 

 

 

 

 

/ (а,

0)

= 0,

 

р (а, 1) =

Р

(or, 1),

 

 

/(а ,

1)

= Р ( а ,

1),

Р (<х, 1 +

г)

=

0.

(

'

Первые два условия (3.30) означают, что жидкие частицы не отде­ ляются от дна и не проникают внутрь смежной среды через поверх­ ность раздела; третье уравнение представляет собой условие непре­ рывности давления при переходе через поверхность раздела, а четвер­ тое уравнение есть условие равенства нулю давления на свободной поверхности.

Теперь предположим, что все переменные, входящие в граничную задачу (3.28), (3.29), (3.30), допускают разложения по целым степеням е:

/= /( ц , г], е)= 2 /(h> (а,

г|) вй, F = F{o, г\, е) =

2

(а, т)) гп и т. д.

fc= 0

 

fc= 0

(3.31)

 

 

 

Законность разложений

по целым степеням по

существу доказана

в[10]. Так как при некоторых высших производных, фигурирующих

вуравнениях (3.28) и (3.29), имеется коэффициент е, эти уравнения

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

Если подставить разложения (3.31) в (3.28), (3.29) в (3.30), а затем приравнять коэффициенты при одинаковых степенях е, то мы получим

последовательность граничных задач.

== t/<°) (—оо, ті) = 1,

Если мы предположим, что

ы<°) (—оо, ц)

то решение уравнений нулевого

приближения

получается непосред­

IX . Б И Ф У Р К А Ц И О Н Н Ы Е Я В Л Е Н И Я

В ТЕО РИ И П О В Е Р Х Н О С Т Н Ы Х

ВОЛН 163

ственно и имеет вид

 

 

ц(0>=1,

£У(0> = 1,

 

/(0) = гі,

Я°> = т1,

(3.32)

р<0>= / ( 1 + pr —tj),

P<°> = / ( l + r - r i ) .

 

Это решение соответствует параллельному течению со скоростью с в обоих слоях.

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

ляет собой линейную систему

 

 

 

 

W + lfë' + p F ^ О,

 

U$> + IF$' + P2> = 0,

 

+

1 ,

 

IF“>+ P*> = \,

 

 

О < т) < 1,

 

 

1 < Л < 1 + г,

/«>+ и«>= О,

 

Л,и +

1/<1>=

0,

(3.33)

/ $ = 0 ,

 

=

0,

 

 

/<11(о,

0) = 0, / а>(о,

1), Р(1,(<4 1) =

рРа ,(°. 1),

 

 

 

 

Ра >(0, 1+ г) = 0.

 

Решения системы (3.33) даются равенствами

 

 

Ц(1)=

— Ö! (0 ),

6Г‘1>=

Л1(0 ) (ті — 1) +

аі (о),

 

/а> = Й1 (а)т],

Ра>=

— Лі (0),

 

 

(3.34)

ра>= — toi (°) + аі (о) +

Р (1) =

ІА і (о) (т) — 1) — toi (0 ) +

 

+ (ті— 1 —pr) + ^i,

 

+ Ai (а ) + ('Ч — 1 — r) + Ki,

где ki

и Ki — постоянные. В силу граничного условия для

давления

йі (0 ) и Аі (0 ) должны удовлетворять уравнениям

 

 

( 1 - / + р / К ( 0 ) - р Л ; ( а ) = О,

 

 

/ а ; ( 0 ) + ( / г - 1 ) Л ; ( 0 ) = О.

 

{ d ' J 5 )

Если обе производные а[ и А[ не обращаются в нуль тождественно (что приводит к однородному потоку), то определитель, составленный из коэффициентов системы (3.35), должен равняться нулю, откуда следует, что

1 --

1-{-г± [(1 — г)2 +4/р]1/2

(3.36)

2r (1 — р)

 

 

Это тот же результат, что и в равенствах (2.28),

выведенных из

линейной теории малой амплитуды. Таким образом, бифуркация имеет место при критических скоростях (см. рис. 8).

Чтобы найти амплитуды й4 и А і , м ы должны обратиться к прибли­ жению второго порядка:

и $ у = а[ (0 )т] — aj(0 )a;( 0 ) — I f ö ' — p S ' , 0 < т ] < 1 ,

Uо1' = А[ (о) 1) + а[ (о) — AI (о) А[ (о) IFT Р а \ 1< т ]< 1 + г .

И*

164

Д Ж . ДЖ . С Т О К Е Р

Решение этой системы включает новые функции а2 (а) и А 2 (а), кото­ рые должны удовлетворять двум следующим обыкновенным диффе-

U

СС

4

с

Рис. 8.

ренциальным уравнениям, которые выводятся из граничного условия для давления:

 

 

 

(1 — I -f- р/) й2(*г) рА2(о)

11

 

 

 

К (о)+(//•— і)л ; (о)= J 1

где

 

 

 

 

 

/ = 4

(3 -

/ +

рі) а["'(о) -f (р-

1 )аі (о) + За, (о)

(а) — ЗрЛ* (о) А' (о),

J = 4

(/ -

6г+

Зіг2) аГ(о) + 4

(/г3 - 3 г2) А"'{о) +

 

+ а\ (о) + гА[ (о) — 3А[ (о) Аі (о).

Так как определитель, составленный из коэффициентов системы (3.38), обращается в нуль, величины I я J должны быть связаны соотноше­ нием

(Іг — 1) / + p j = 0.

(3.39)

Используя (3.35) для упрощения (3.39), получаем

 

т0а["(а) = т д

(о) а[ (о) + /п2а; (а).

(3.40а)

Коэффициенты т0, ти т2, зависящие от г, р,

I я К и задаются равен­

ствами

 

 

 

 

Зт0 = г (г + I) I — (г2 + 1 + Зрг),

 

 

ті =

_ з [/2 (Г2 — 1) +

/ (1 — 2г) +

1 У(1г — 1),

(3.40Ь)

т2 =

—[2 /г (р — 1) +

г + 1 ] + 3 (1

— /) KJ(lr 1).

 

Общее решение уравнения (3.40) можно выразить через эллиптиче­ ские функции Якоби. Эти решения приводят к кноидальным волнам. Уединенную волну можно рассматривать как предельный случай кноидальных волн или ее можно найти непосредственно путем инте­ грирования уравнений (3.40) при условиях затухания. В результате

IX . Б И Ф У Р К А Ц И О Н Н Ы Е Я В Л Е Н И Я

В ТЕО РИ И П О В Е Р Х Н О С Т Н Ы Х

ВОЛН 165

получаем решение

 

 

 

Й1 (а) =

sch2 Т (т2Іто)1/2,

(3.41)

Аі (ѵ) =

- ^ г

-

(3.42)

При Н = О это решение переходит в решение для одного слоя, полученное в [9].

Для уединенных волн в двуслойной жидкости возможны два типа поведения:

(i) При большей критической скорости все линии тока проходят над соответствующими линиями для однородного потока. Отношение

амплитуды

волны на

свободной

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

по поверхности раздела

равно

1/(1 — 12г). Если б ~

Д, то это отно­

шение приближенно равно 1 +

г.

Если б =

Д и h =

Н, то это отно­

шение равно 2 .

 

 

скорости

имеются

различные воз­

(ii) При

меньшей критической

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

на

свободной

поверхности ведет себя противоположным образом.

по

В

случае

жидкости с плотностью, экспоненциально убывающей

у,

задача

разбирается в [7].

ЛИТЕРАТУРА

[1] Стокер Дж. Дж., Волны на воде, ИЛ, М., 1959.

[2]Weinstein А., Sur un problème aux limites dans une bande indéfinie, C.R.

Acad. Sei., Paris, 184 (1927), 497.

[3]Levi-Civita T., Determination rigoureuse des ondes permanentes d’ampleur finie, Math. Ann., 93 (1925), 264—314.

[4]Struik D. J., Determination rigoureuse des ondes irrotationelles périodiques

[5]

dans un canal ä profondeur finie, Math. Ann., 95 (1926), 595—634.

 

Peters A. S., Unpublished manuscript.

 

 

6

(1953),

[61

Stoker J. J.,

Unsteady waves on a running stream, C.P.A.M.,

[7]

471—481.

 

 

 

 

 

Peters A. S., Stoker J. J., Solitary waves in liquids having non-constant den­

[8]

sity, C.P.A.M., 13 (1960), 115—164.

the shallow water

theory,

Appendix

Friedrichs K-

O., On the derivation of

 

to J. J. Stoker, The formation of breakers and bores, C.P.A.M.,

1

(1948),

[9]

81—87.

The solitary wave and

periodic waves

in shallow

water,

Keller J. B.,

 

C .P.A .M ., 1

(1948), 323—339.

 

 

 

 

[10]Фридрихе К. О. и Хайерс Д. Г., Существование уединенных волн, в сб.

«Теория

поверхностных волн»,

ИЛ, М., 1959, стр.

45—183.

близ­

[П] Литтмен

У., О существовании

периодических волн при скорости,

кой к критической, в сб. «Теория поверхностных

волн»,

ИЛ, М.,

1959,

стр. 185—216.

довгих хвиль, 36.

Праць

Інст. Матем.

[12*]Лаврентьев М. А., До теоріі

АН УРСР, № 8 (1946), 13—69

 

 

 

 

\

166

ДЖ . ДЖ . СТОКЕР

[13*]Моисеев Н.

Н., Вводная статья в сб. «Теория поверхностных волн», ИЛ,

М., 1959,

стр. 5—27.

[14*]К.расовский Ю. П., О существовании апериодических течений со свобод­ ной границей, ДАН СССР, 133, № 4 (1960), 768—770.

[15*]Красовский Ю. П., К теории установившихся волн конечной амплитуды,

Журн. вычисл. матем. и матем. физики, 1, № 5 (1961), 836—855.

[ 16*]Тег-Кгікогоѵ А. М., Theorie exacte dés ondes longues stationnaires dans un liquide heterogene, J. de Mecanique, 2—3 (1963), 351—376.

[17*]Тер-Крикоров А. M. и Треногин В. А., Существование и асимптотика реше­ ний типа уединенной волны для одного класса нелинейных эллиптических уравнений, Матем. сб., 62 (104), № 3 (1963), 264—274.

[18*]Тер-Крикоров А. М. и Треногин В. А., Решения типа длинных волн для квазилинейных эллиптических уравнений в неограниченной полосе. Диф­ ференциальные ур-ния, 3, № 3 (1967), 496—508.

[19*]Тер-Крикоров А. М. и Треногин В. А., Решения типа длинных и уединен­ ных волн для квазилинейных эллиптических уравнений в неограниченной полосе, Тезисы докладов Всесоюзной межвузовской конференции по при­ менению методов функционального анализа к решению нелинейных задач, Баку, 1965.

[20*]Треногин В. А., Существование

и асимптотика

решений типа

«уединенной

волны», ДАН СССР,

156, № 5 (1964), 1033—1036.

и диференци-

[21*]Треногин В. А., Некоторые вопросы теории функциональных

альных уравнений в

банаховых

пространствах,

Докторская

диссертация,

М., 1968.

 

 

 

 

X

ТЕОРИЯ ВОЗМУЩЕНИЙ УСЛОВНО-ПЕРИОДИЧЕСКИХ РЕШЕНИЙ ДИФФЕРЕНЦИАЛЬНЫХ УРАВНЕНИЙ

Ю. Мозер

1. Введение

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

О п р е д е л е н и е . Функция / (t) называется условно-периодиче­ ской, если она может быть представлена в виде

f

(t)

= F (соit, . . ., сdnt),

(1.1)

где F (хи . . ., хп) — непрерывная функция с периодом 2л

по каж­

дому из переменных хи

.

. ., хп.

 

Такая функция имеет представление Фурье

 

 

 

/ = 2 cjeW’

(1.2)

 

 

і

 

где / = (/!, . . ., /„) есть п-мерный вектор с целочисленными

компо­

нентами и

 

 

 

(/'. ® ) = S / > v

О-3)

 

V — 1

 

Мы предполагаем, что кц, . . ., сод рационально независимы, т. е.

если (/, со) = 0,

то / = 0. Это предположение не приводит к потере

общности, так

как множество соь . . ., со„ можно всегда привести

к рационально

независимому множеству х). Мы называем

со,, і =

= 1 , . . ., п, базисными частотами.

Из (1.2) видно, что класс условно-периодических функций замкнут относительно сложения и умножения. Эти функции — почти периоди­ ческие с конечным частотным базисом. Периодическая функция является условно-периодической с единственной базисной частотой. Заметим, что множество чисел вида (/, ш) (спектр [) плотно на дей­ ствительной прямой, если п > 1 2)*. Ради простоты мы предполагаем, что все функции, встречающиеся в этой лекции, действительные аналитические.

То есть заменить базисом, состоящим из меньшего числа частот.— Прим.

перев.

соп.— Прим, перев.

2) Вследствие рациональной независимости coj..........

168

Ю. М О ЗЕР

2. Периодические решения *)

Прежде чем обсуждать вопрос об условно-периодических решениях обыкновенных дифференциальных уравнений, напомним о некоторых характерных трудностях, связанных с периодическими решениями. В качестве типичного примера рассмотрим уравнение Дуффинга

хи + ос2* + ßx3 = е/ (t).

(2 .1)

Здесь г — малый параметр и / (t) имеет период 2я по t. Уравнение (2.1)

является возмущенным по отношению

к автономному

уравнению

хи + ос2х + ßx3 = 0 .

(2 .2 )

Уравнение (2.2) имеет интеграл

 

 

x? + oc2r ä + - |x 4 = £'

(const),

(2.3)

являющийся для положительных е уравнением семейства замкнутых кривых Се в фазовой плоскости; Е — параметр семейства (рис. 1).

Рис. 1.

Для каждого Е параметры а и ß можно выбрать так, что любое решение (2.2) с начальными данными на Се имеет период 2я. Но мето­ дом Пуанкаре можно показать, что существуют только две начальные точки на Се , о т которых ответвляется 2 я-периодическое решение

возмущенной системы (2.1). Таким образом,

если х

(і, е) — решения

(2 .1) с периодом 2 я, то пределы lim х (t, е)

не образуют множества

 

 

6 —>о

 

 

 

 

2я-периодических реш ений

уравнения

(2.2)

с

начальными данными

на Се- Эти ж е

замечания

остаются в

силе,

если мы

сделаем зам ену

переменны х х =

г ^ 3г, которая приводит (2.2)

к виду

 

 

ztt + a 2z = гУЗ ( f — ßz3).

(2.4)

(Уравнение (2.4) представляет собой возмущение линейной системы.) Поэтому требование сохранения периодических решений при возму­ щении может повлечь за собой наложение некоторых дополнительных условий на параметры задачи.

Заметим, что общая неавтономная периодическая система

 

Уі = f (У, t) — F {у, at),

(2.5)

*) Обсуждаемые в этом параграфе понятия рассматриваются в стандартных учебниках по обыкновенным дифференциальным уравнениям и нелинейным' колебаниям.

X. Т Е О Р И Я ВОЗМ УЩ ЕН И Й У С ЛО В Н О -П ЕРИ О Д И Ч Е С К И Х Р Е Ш Е Н И Й J6<>

где/имеет период — по t, а F — период 2я по сat, эквивалентна авто­ номной системе

xt =

СО,

(2.6)

yt = F (у, х)

 

с начальным условием х (0) =

0. Поэтому можно рассматривать,

периодические решения неавтономных систем, изучая расширенную автономную систему (2.6). Напомним стандартную теорию возмущений для таких систем.

Рассмотрим автономную систему

 

zt = h (z, г).

(2.7)

Пусть (2.7) при е = 0 имеет 2я-периодическое решение z0 (t). Тогда для достаточно малых е > 0 возмущенная система (2.7) имеет реше­ ние г (t, е) с периодом Т (г) по t, такое, что г (t, 0) = z0 (t) и Т (0) =

= 2 я, при условии что множество 2 я-периодических решений z (t) уравнения в вариациях

Zt = hz 0 (t), 0) г

(2 .8)

одномерно, т. е. г (t) — az0 (t). Поэтому существование периодиче­ ского решения возмущенной системы зависит от условия невырожден­ ности в терминах уравнения в вариациях. Мы увидим, что подобные результаты имеют место и для условно-периодической задачи.

3. Типичные задачи для условно-периодических движений

 

Мы будем иметь дело с обобщением системы (2.6)

 

Xt =

со,

(3.1)

Уі =

■ ’

ч

F

(х, у),

 

где х и у — соответственно п- и m-мерный вектор. Функция F (х, у) имеет период 2 я по переменным xt, и соь . . ., со„ рационально неза­ висимы. Поставим вопрос: продолжают ли существовать условно­ периодические решения системы (3.1) при ее возмущениях? Ответ на этот вопрос, вообще говоря, отрицательный, потому что малые возмущения могут так изменить значения сог, что их рациональная независимость нарушится.

Рассмотрим случай п — 2. Эту задачу удобно изучать в плоскости

хи х2 (рис. 2), в которой F (хъ xz, у) двоякопериодична.

Уравнения

X\t — ©1, x2t = ©2

(3-2)

соответствуют семейству прямых в плоскости xlt х2 с иррациональным углом наклона ©2/©і (ср. рис. 2). Отождествляя линии хг = 0 с х± = = 2 я и х 2 = 0 с х 2 = 2 я, можно рассматривать xt и х 2как координаты

Соседние файлы в папке книги из ГПНТБ