книги из ГПНТБ / Теория ветвления и нелинейные задачи на собственные значения
..pdf160 |
ДЖ . ДЖ . СТОКЕР |
|
||
уравнение неразрывности |
имеет вид |
|
|
|
|
их + ѵу = |
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как координаты