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

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

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

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

Добавлен: 29.10.2024

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

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

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

Интеграл энергии. Первые два члена подынтегрального выраже­ ния представляют собой плотность кинетической энергии, а тре­ тий — плотность потенциальной энергии (в безразмерной форме). Кроме того, имеются интеграл уравнения (6. 33), выражающий закон сохранения потенциального вихря

 

 

 

 

 

Лdxdij (Дф — ф) = const,

 

 

 

(6. 35)

и интеграл,, выражающий закон сохранения массы

 

 

 

 

 

 

 

J

dxdyè

=

 

const.

 

 

 

 

(6. 36)

Представим решение для слабо

взаимодействующего триплета

волн в

виде

Ф

а, cos Gj +

а.2

соя 02 + а3 cos Ѳ3,

 

 

(6. 37)

где Ѳн обозначают фазы волны —

Qn= h Hx -\-lny—

 

+ е„ + = 1 ,

2, 3),

 

t.

 

 

 

а амплитуды

ап (t)

представляют

собой

медленно меняющиеся

функции

одного только

времени

 

 

Подстановка

(6. 37) в (6. 33)

приводит к уравнению

(у.| +

1)

а2

cos Ѳ2 +

(-/•! + 1)

а3

cos 03 =

 

 

(-/.2 + 1)0! cos Oj +

 

 

 

 

 

=

С ха2а3

[спя (02 + 03) — cos (0„ — 03)] +

С 2а3ах

[cos (03 + Ѳх) —

 

— cos (Ѳ3 — 0j)] +

С3аха2

[cos (0j + Ѳ2) — cos (0X— 02)J,

(6. 38)

 

 

 

 

где у-„=\кп\2=Щ-\-11, а через C x, C2, C 3 обозначены коэффициенты взаимодействия, так, например,

С г = \ (4 4) (z tk2 X М

и т. д.,

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

(6. 32),

поэтому

z • [k2 X k3] =

z ■ [k3 X kj] — z • [kx X

k2] = 2b,

(6.39)

Сг =

Ь № - * * ) ,

C 2 =

b ( 4 - - 4 ) ,

C3 =

b ( 4 - 4 J .

Если, далее,

фазы en связаны соотношением

 

 

 

 

 

ei +

е2 + £з = О,

 

 

 

то три члена в левой части

(6. 38) уравниваются тремя членами

в правой части при условии, что

 

 

а2а3,

 

 

 

 

(1 +

у.2)

йг = Ь (4 —

у-І)

 

 

 

 

(1 +

4) â2 — b (;4

— у?) а3а ,,

 

(6. 40)

 

 

(1

+ ' 4 ) ä3 = b (4 — 4) аій2-

 

 

 

 

 

 

 

100


Остальные члены в правой части в общем случае не соответствуют свободным волнам и поэтому не будут непосредственно участ­ вовать в резонансном переносе энергии. Уравнения (6. 40) можно получить также следующим путем. Умножая, например, (6. 38) на соз0х и усредняя по быстро меняющимся фазам, мы и прихо­ дим при условии (6. 32) для медленно меняющейся величины

к первому из уравнений системы (6. 40) и т. д. Система (6. 40) представляет собой гидродинамический триплет, рассмотренный в предыдущих главах. Система уравнений (6. 40) имеет два ин­

теграла движенияа:\

+

(1 +

-4) of + (1 + х|) а* =

const,

(6. 41)

(1 +

-/.?)

(1 +

*!)* а? +

(1 +

4)2 а* + (1 + *2)2

аі =

const,

(6. 42)

 

 

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

Если обозначить через Дн изменение величины а2 в течение времени от момента £„ до момента t^>tQ, то с помощью интегралов (6. 41), (6. 42) получаем соотношения

(1 +

4) \ +

(1 +

4)

+

(1 +

4) Аз =

о,

(6.43)

(1 +

х?)2 Ді +

(1 +

х|)2 Д2 + (1 +

х2)2 Д3

= 0,

из которых следует:

4 ± 4 д 2:

1 +

3-2 ад3-

 

(6. 44)

 

1 ± А а

 

— *.?

4

 

 

 

 

 

Пусть, например, хj <Ѵ.|<У з.

Тогда

знак

Д2 всегда

противополо­

жен знакам Дх и Д3, т. е. поток энергии в сторону двух крайних волновых чисел хх и х3 всегда имеет один и тот же знак, противо­ положный знаку потока энергии в сторону промежуточного

волнового

числа х2.

Однако что касается знака потока энергии,

например,

в сторону волнового числа х, то он зависит не только

от знака величины аха2а3, но также от знака величины

b

(4 —

*§)■

 

 

 

С {

 

 

 

Отметим, что в случае, когда два из волновых чисел xls х2, х3

равны между собой,

соответствующий коэффициент

 

обраща­

ется в нуль.

 

 

 

 

 

Уравнения гидродинамического триплета подробно обсужда­

лись в предыдущих главах, и поэтому мы не будем останавли­ ваться на их решении.

Отметим, что аналогичным образом можно рассматривать резонансное взаимодействие волн Россби—Блиновой и на сфере. В этом случае разложение следует проводить по сферическим гармоникам. Таким же образом можно рассматривать и задачу, представляющую большой интерес для метеорологов, о резонанс­ ном взаимодействии волн как между собой, так и об их взаимо­ действии с зональным потоком (Галин (1974)).

101


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

В заключение отметим также, что аналогичным образом можно рассматривать и другие волновые задачи, представляющие боль­ шой геофизический интерес, например, задачи взаимодействия внутренних гравитационных волн (Филлипс (1967)), поверхност­ ных волн (Бреховских (1966)).

Глава VII

СТРУКТУРНЫЕ СВОЙСТВА КВАДРАТИЧНО-НЕЛИНЕЙНЫХ СИСТЕМ

§ 1. Аффинные инварианты и критерий существования квадратичного интеграла для систем второго порядка

В гл! II мы уже исследовали некоторые свойства простейших квадратично-нелинейных систем при линейной замене переменных Естественно ввести общие понятия эквивалентности квадратичнонелинейных систем. Две такие системы называются эквивалент­ ными, если путем некоторой замены переменных, осуществляе­ мой с помощью обратимого линейного преобразования (детерми­ нант соответствующей матрицы отличен от нуля), можно добиться численного равенства компонентов взаимодействия Г* .fc, опреде­

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

Регулярность и существование квадратичного интеграла движения являются примером структурных свойств, с помощью которых в гл. II был выделен специальный класс СГТ. Число

независимых линейных интегралов системы также является струк­

турным свойством. Мы будем называть это число

индексом вы­

рождения

системы. Если линейных интегралов движения вообще

ие существует, то система называется невырожденной.

симме­

Примером однократно вырожденной системы является

тричный гироскоп, для которого / 1= / .2= / и уравнения

Эйлера

имеют вид

 

 

7 = — (7з — у)

(7- й

Поскольку в этом случае существует интеграл момента, то (о3= =const, и уравнения движения становятся фактически линейными.

103


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

~ = х „ А І х к + 1

/, /, к < п - 1,

(7. 2)

где матрица /1гА= Г ‘Л.. В силу постоянства хн А), определяет линей­ ную часть прогностического оператора. Встречающиеся в при­ ложениях смешанные системы, уравнения движения которых содержат линейную и квадратичную части, формально можно рассматривать как чисто квадратичные нелинейные системы в про­ странстве состояний, размерность которого на единицу больше порядка данной системы, но допускающие линейный интеграл.

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

 

 

 

 

 

 

 

 

=

% х Ѵ

S = lU аікх'хк,

(7.3)

и мы хотим отыскать квадратичную форму

 

 

 

 

являю­

щуюся интегралом движения, т.

е.

такую,

для которой

d S/dt=

О,

то из уравнений движения для неизвестных коэффициентов

а.к

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

 

 

 

 

 

 

 

 

 

 

 

 

aJm

 

akmV[nj =

 

 

 

 

 

 

 

 

 

 

 

 

Іѵ =

I'7,--+

 

 

0-

 

 

(7. 4)

М

 

п

 

«„„г;, +

 

 

 

п

 

 

таких уравнении

 

 

-----

 

 

 

 

 

 

 

Число

 

 

 

 

а число неизвестных

 

=

 

 

 

N

 

М

 

небольших

 

значений

 

 

получаем

 

 

 

^ . В частности, для

 

 

 

следующие

значения

 

и

 

:

 

 

 

 

 

 

 

 

 

 

 

 

 

п

 

 

2

 

3

4

 

5

 

 

 

 

 

 

 

 

 

N

 

 

4

10

20

35

 

 

 

 

 

При

 

2

М

 

 

3

 

6

10

15

 

 

 

 

 

 

число уравнений

всегда

больше числа неизвестных»

и условие совместности системы (7. 4) накладывает существенные

ограничения

на

коэффициенты

взаимодействия ІЧ,..

Система

(7. 3) допускает

положительно

определенный интеграл

только

при условии,

что система (7. 4)

оказывается совместной и среди

104


возможных решений удается отыскать такие, для которых форма S является положительно определенной. Следует заметить, что такой прямой способ исследования приводит к крайне громозд­ ким вычислениям уже при п = 3.

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

инвариантным способом.

системы можно построить

аффинор:

По уравнениям движения

Г (х), представляемый в любой системе координат матрицей

ГНх) =

^ = ]% * ’".

(7.5>

 

ОХІ

 

С помощью афинора I1(х) можно записать уравнения возмущен­ ного движения вблизи заданной траектории x' (t). Полагая

.г* =£'(£) -ҢЛ где x'{t) удовлетворяет исходным уравнениям движе­ ния и пренебрегая квадратичным по у членом, получаем

У =

Г (х)у .

-

(7.6>

Отсюда следует, в частности,

что собственные значения аффинора

Г (х) являются некоторыми скалярными функциями вектора х и имеют размерность обратного времени.

Для квадратично-нелинейных систем характеристический аф­ финор зависит от векторного аргумента линейно. Условия сим­

метрии тензора ГД

по пижним индексам можно записать в следу­

ющей форме:

Г (х)у = Г(у)х.

(7.7)

Линейное аффинорное поле Г (х) в фазовом пространстве, удов­ летворяющее условию (7. 7), взаимно однозначно связано с урав­ нениями движепия системы. Вместо характеристических чисел

матрицы ||Г*.|| (аффинора

Г)

можно рассматривать величины

—а

=

Sp Г2 (х), Хз = Sp Г 3 (х), . . . (7.8)

(х) = Sp Г (х), X =

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

определяющиепсимметричные тензоры\+П соответствующего порядка.

Среди этих форм функционально

независимыми

 

могут быть

не более чем

форм, так

как '

^н+2> • •

• можно

выразить

через формы Xj

Х2,

. . . .,

Хп,

 

0.

т

 

 

Для регулярных

систем Sp Г (х) = — о =

 

при исполь­

Дифференцирование однородной формы порядка

 

зовании квадратично-нелинейных уравнений движения

приводит

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

105