Файл: Нелинейные системы гидродинамического типа..pdf

ВУЗ: Не указан

Категория: Не указан

Дисциплина: Не указана

Добавлен: 29.10.2024

Просмотров: 86

Скачиваний: 0

ВНИМАНИЕ! Если данный файл нарушает Ваши авторские права, то обязательно сообщите нам.

нѳй» координаты при этом определяется инвариантным способом

z2 =

 

Легко][видеть, что функция

(6. 21)

Н = Im (zjz2z3) = рС°5рС°)рС0>sin cp= const,

dt

 

т. e. являетсяdlинтегралом движения, так как

 

(ZjZ^g) = — I z,z312 + I z3Zl I* — I zxz21

вещественно. Пользуясь функцией H , можно записать уравнения движения в гамильтоновской форме. Для этого следует положить

dz

 

2і — Р і + г'?і> .

z2 — ?2 + Фг>

 

z3 =

Рз~\~ Чз’

1

. /

I

\

 

dH . .dH

\ I

- д т =

1д£1=(РіЯз — ЯгРз) +

г (?г?з +

Р-гРз) = — 2г2з-

Другие уравнения проверяются аналогично.

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

1. Н = 0. Это означает, что ср=0, т. е. можно выбрать систему координат, в которой в некоторый момент времени все координаты zx, z2, z3 вещественны. Поскольку параметры системы также веще­ ственны, то мы приходим к классической задаче движения гиро­ скопа, решение которой можно найти в любом курсе механики (см., например, книгу Ландау и Лифшица (1958)). При фиксиро­

ванном

Е

и уменьшении

S

период возрастает.

Случай 5 = 0

от­

вечает

апериодическому движению по сепаратрисе.

 

 

 

 

 

 

ff=ft=

 

 

 

 

 

 

 

 

 

 

 

 

2. 5 = 0Z,j = z 3, 0. В этом случае [ zx | = | z31, выполнив соответству­

ющее фазовое преобразование, можно добиться того,

чтобы

при i= 0

 

а

благодаря инвариантности

уравнений

движе­

 

 

 

 

 

>zx

 

 

 

 

 

 

 

 

 

ния при замене zx—>z3, z3— это равенство сохраняется для лю­

бого момента времени.

 

 

 

 

 

 

 

 

 

 

 

Полагая

z3=2w ,

z1= z 3=\/2z, приходим к уже исследованному

уравнению (6. 11)

 

 

dz)dt

 

z2, wz.

 

 

 

 

 

 

 

 

 

 

 

 

 

dwjdt =

 

 

 

 

 

 

 

 

В общем

случае, когда

 

— —2

 

 

 

проводится

ана­

5=^=0,

исследование

логично. Существует максимальное значение

Н ,

зависящее от

Е

и

5, при котором все координаты меняются\

по

 

чисто гармониче­

скому закону и обмена энергией при этом не происходит;

частота

при этом максимальна. При уменьшении

Н\

период увеличива­

 

 

ется и превращается в период колебания обычного гироскопа при //->0.

Заметим, что энергии «крайних компонент» меняются синфазно и противофазно по отношению к энергии средней ком­ поненты z2.

96


Детальный анализ соответствующих колебаний имеется в статье Цельмана (1971).

Связь рассмотренной выше системы с задачей о колебаниях связанных осцилляторов становится очевидной, если рассмотреть

величины

 

и =

 

 

 

 

и = гZje1'“'*,

(6.22)

где

и,шх,V, (Dg

произвольны,

3еіш>(,

w — г2еіш‘1,

 

а ш2= ш 1 + 0)з (все частоты положитель­

ные) и комплексные амплитуды удовлетворяют уравнению (6. 19),

то

w

подчиняются

уравнению гармонических

колебаний

 

 

 

 

с дополнительными членами, описывающими нелинейные воз­

действия (в квадратичном

приближении)

 

j f =

ia>і» +

Р > “0,

(6.23)

^ = io )3y +

;•(£?«;),

^ = m9w + q (ну).

Если при этом pr^>0, то возможен «распадный режим», когда колебание с наибольшей частотой ш2 в течение длительного времени (по отношению к периоду Т 2=2п/щ) отдает свою энергию колеба­ ниям с меньшими частотами.

Основная задача при исследовании конкретных колебательных процессов (включая волны) заключается в расчете параметров р, г, описывающих нелинейность системы.

Конкретные примеры рассмотренных процессов обмена энер­ гии при нелинейном резонансе ( ш2= -[- ш3), представляющих интерес для геофизики, можно найти в сборнике переводных статей «Нелинейная теория распространения волн» (см., напри­ мер, Филлипс (1967)), а также в статьях Бреховских (1966, 1968). Заметим, что уравнение триплетного взаимодействия, по-види- мому, раньше всего было использовано в задачах нелинейной оптики (Бломберген (1966)) (см. также обзор Кадомцева и Карп-

мана (1971)).

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

7 Нелинейные системы

97


то ее всегда можно преобразовать к виду (6. 9) и, следовательно, записать в комплексной форме в виде системы второго порядка с комплексными координатами. Наличие такой группы обычно можно установить, исходя из условий физической задачи. Так, например, оно имеется в задаче о движении жидкости в кольце­ вом канале и в других задачах с аксиальной симметрией.

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

В качестве конкретного примера применения вышеизложенной методики рассмотрим, следуя Лонге-Хиггипсу и Гиллу (1967), резонансное взаимодействие планетарных воли (воли Россби— Блиновой).

§ 2. Резонансное взаимодействие планетарных волн

«Планетарными волнами» или волнами Россби—Блиновой (см. работы Россби (1939), Блиновой (1943)) называются гиро­ скопические волны, вызывающие глобальные движения, проис­ ходящие в слое жидкости или газа, покрывающей вращающийся земной шар.

Пусть (х, у) — координаты в касательной плоскости сферы, причем X отсчитывается па восток, а у — на север. Основным динамическим уравнением, описывающим планетарные волны, является уравнение сохранения потенциального вихря

 

 

^ [ ( Д - « 2)ф] +

Р Й = 0 ,

(6.24)

где

a r = f2/gh,

ß—dfldy,

g

и / обозначают соответственно ускорение

силы тяжести

и параметр Кориолиса,

h

— средняя

глубина

жидкости, ф — функция тока, а через

D/Dt

обозначено дифферен­

цирование в системе,

движущейся

вместе с жидкостью.

 

Уравнение (6. 24) является приближенным, а именно, в нем оставлены только наиболее важные нелинейные члены, которые определяют горизонтальную адвекцию вихря. Полный вывод динамического уравнения, выражающего сохранение потенциаль­ ного вихря, содержится, например, в работе Обухова (1962) и монографии Кочина, Кибеля, Розе (1963).

При рассмотрении решений уравнения (6. 24) удобно перейти к приближению ß-плоскости, в котором а и ß принимаются по­ стоянными. Тогда, выбирая единицы измерения длины и времени

соответствующим образом, можно положить в (6.24)

a = ß = l и

переписать уравнение в виде

 

l , [ ( 4 - l ) t J + g = ° .

(6.25)

98



Линеаризированное уравнение (6. 25) при этом имеет вид

* Г ( Д -1 ) ф ] +

** = 0.

(6. 26)

Это уравнение имеет гармоническое решение

(6. 27)

ф = а ехр (г (кх +

Іу °t)}

при условии, что о удовлетворяет дисперсионному соотношению:

а (/с2 —f- Z2 —J—1) — А: = 0.

(6.28)

При этом компонента фазовой скорости в направлении с востока на запад определяется соотношением

а __

1

(6.29)

кк - - |- / “ - f - 1

Так как выражение в правой части (6. 29) всегда отрицательно, то волны всегда распространяются на запад. Эти волны и называ­ ются волнами Россби—Блиновой.

Рассмотрим три свободные планетарные волны, определяе­ мые в первом приближении выражением

Ф., = ѳхр (і (Л:„ж + г„і/ — ant)}, (п = 1, 2, 3).

(6.30)

Для каждой волны удовлетворяется дисперсионное соотношение

(6. 28), т. е.

°п (^1 + + 1) + /сп = 0. (г а = 1) 2, 3). (6.31)

Ддя резонансного взаимодействия этих волн необходимо, чтобы,

кроме того, выполнялись бы еще условия

 

 

/Сі -f-

к2

-f-

к3 =

0,

 

 

l\

Z2 ~Ь

=

 

 

 

+

 

(6.32)

 

°i 4* а2+

аз — 0.

Отметим, что условий, накладываемых на

девять величин аи,

кп, Іп,

всего шесть. Поэтому имеется большой набор волн, для кото­

рых выполняются соотношения

(6. 31), (6.

32). При заданном

волновом векторе, например к2=(/с2, Z2), подробный анализ вол­ новых векторов Iq, к3, при которых имеется резонансное взаимо­ действие, содержится в цитированной выше работе Лонге-Хиг- гинса и Гилла.

Рассмотрим теперь динамику взаимодействия этих волн. Рас­ крывая производную D/Dt в явном виде, можно записать урав­ нение (6. 25)

<в- аз)

Уравнение (6. 33) имеет ряд интегралов движения. Это, во-первых

! [т(ё) + т © = const (6.34) dxdlJ

7* 99