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

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

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

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

Добавлен: 29.10.2024

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

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

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

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

/

Рис. 25. Спектральное

рас­

/

J

S

7

пределение амплитуд при п = 10

 

 

 

 

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

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

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

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

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

2

85

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

 

 

 

 

йо =

41)уіуа +

/о.

 

+

/х,

 

 

 

 

г>! =

с^ѵ.2ѵ3

+

c[z\ v .2

 

 

 

 

ѵ„

=

сО)п3к4 -f- с|2)у1к3 +

4% 0f lt

 

 

4

 

...............................(5.7)

 

 

 

 

 

=

^Чч-Лч-З + Ч'-Л+l + <fЧ-Л-2.

 

 

4, =

C»VA +2 + Ф - Л - а — K Vn>

 

Vj

4+i =

4

v

»

- i - \

, A tr

 

 

Здесь

 

— скорость

г-й

моды,

с.

— параметры взаимодействия;

 

 

на достаточно высокие моды, начиная с тг-й, действуют силы внут­ реннего трения, описываемые коэффициентами Хл; /0 и /х — внеш­ ние силы.

- В системе (5. 7) выписаны лишь такие нелинейные члены, ко­ торых достаточно для выполнения в отсутствие внешних сил и вязкости двух законов сохранения, выражаемых квадратичными интегралами:

*=4 2*

(5-8)

где кі — волновое число (Ь . —і/к— масштаб г-й моды). Инте­ гралы Е и Й являются аналогами интегралов энергии и квадрата вихря для двумерного движения.

Для того, чтобы d E ld t = 0 и dQ,ldt—0, необходимо выполнение следующих условий:

ЩсѴ

Щ Ср

+

С(2)і + С(3)2 =

0,

 

 

+

+1

с?+1 +

Щ+2с(%

= О,

2

= 0,

1, . . . .

Отсюда легко выразить

все коэффициенты,

например, через

 

 

-с.‘- 1

— frf+i.

с(3)__

п

 

k

— 4-1

сО) = с, cf) =

Щ - А?+1

 

Щ-2 U - Щ

і3-,

 

 

 

 

^-2

 

 

 

і = 1,

2, . ..

 

г = 2, 3, . . .

 

 

 

 

 

 

 

 

 

 

 

И система (5. 7) примет вид

 

 

— ci-i

А?_і - А?+1

+

.

vt = с іѵш ѵі+2

/c;+i

+

.

A?_2 — AJ_t

 

.

. ,

+ 4 - 2

 

U 4 - l 4 - 2 —

1,

U

 

г = 0, 1,

. . . , гг +

 

(5-9)

(5.10)

86


где, по определению,

можно считать,

что

 

 

 

 

 

 

 

 

 

 

c_j

=

с_

о

с„

=

ся+1 =

О,

 

 

 

 

 

 

 

 

 

 

 

 

гі

 

 

 

 

 

 

 

 

 

 

 

 

/,• — 0,

 

 

 

=

2,

3,

4, . .

 

 

 

 

 

 

 

 

 

 

X. =

0, Cf

>

 

= 0 ,

1, .

. re — 1,

 

 

 

 

 

 

" ,> 0 ,

 

 

 

 

0,

ft<+1> Ä r

 

 

и

ѵ-_2

возрастает

Из (5. 10) легко видно, что при возбуждении

 

 

 

скорость

Vf.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Найдем стационарное решение системы (5. 10) при условии

самоподобия

 

 

п с»+і

 

 

^

 

 

>

 

 

1

,

 

^.+і

 

 

 

 

 

(5.11)

= 2 ,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

ДляС:гV.J.

3,

. .с,

 

* ? -і -

1

 

имеем

 

 

 

 

 

 

ѵі

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

0.

 

 

 

1 А І-Й + , у,--Л+і + с,._. *5-, - X? У ,

 

Или, используя (5. 11),

 

 

 

1 +

g2

 

у . ~

Ly»,

 

 

 

 

 

Отсюда

 

 

 

 

1)

 

p =

g3

 

V»,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

q

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

2)

 

Р =

 

3_1.

и

 

 

 

— = WE

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Ws

 

Рассмотрим потоки энергии

 

WE

 

квадрата вихря

 

 

 

 

W E =

f ^ 0 +

 

fivv

 

 

 

 

dt

 

 

E’

 

 

 

 

 

 

^ в =

^

 

о

 

+

 

*!/ Л ,

 

 

 

 

 

 

 

Из (5. 7),

 

 

 

 

 

 

 

5 F = ^ e -

 

 

 

(5. 10),

(5.WИ )=

имеем

 

 

 

№8] ,

 

 

 

 

 

 

 

 

 

 

 

e

 

 

coV? vI \ ^ —

 

 

 

 

 

 

 

 

 

 

We =

 

— с0%

 

 

 

q*.

 

 

 

 

 

 

 

 

 

 

 

 

3у3 [gjj.3 _

1]

 

 

 

 

 

Таким образом, ѵв случае 1)

 

^ = 0 ,

 

а для

ѵ{

получим

 

 

 

 

 

 

 

 

 

\

=

а2к.

7»,

 

 

а‘

 

 

 

 

g3cp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

.(g2—1) *0

 

 

 

 

 

 

в случае

2)

 

 

и

 

 

 

 

 

 

 

 

 

 

 

 

 

PF,e= =0

 

 

 

 

 

 

 

 

 

 

(g2 - 1 ) со

 

 

 

 

 

 

 

 

V2. Ъ2к~\

 

 

 

 

 

 

 

 

Q3h

 

 

 

 

 

 

 

 

Мы нашли,

что для системы (5.

7),

сохраняющей два

квадратич­

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

87


энергии по спектру. Второй режим связан с передачей квадрата вихря по спектру. Ему соответствует известный «закон минус 3» для энергетического спектра в плоском турбулентном потоке.

Пусть теперь, кроме сил /„ и /х, на моды / — 1, / и / + 1 с до­ статочно большим у действуют силы f /\+1. Тогда для системы (5. 7), (5. 9) нетрудно получить следующее стационарное решение:

И ^ С О ,

JEs =

0 ,

и^а > 0,

РЕв =

0.

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

Лилли (1972)).


Глава VI

СИСТЕМЫ С КОМПЛЕКСНЫМИ ХАРАКТЕРИСТИКАМИ

И ЗАДАЧА О ВЗАИМОДЕЙСТВИИ КОЛЕБАНИЙ

§ 1. Комплексное расширение систем гидродинамического тина

Для описания процессов нелинейного взаимодействия колебаний, в частности воли, интересно обобщение понятия СГТ на случай

координат, представляемых комплексными числами.

 

 

 

Пусть

S 2m

 

 

 

 

 

 

 

 

 

 

система

чет­

— некоторая квадратично-нелинейная

ного порядка, имеющая

квадратичный интеграл.

Мы будем

го­

ворить, что система

S 2m

допускает комплексное

представление,

если при

разбиении

координат

состояния иа

пары

{хъ

г/Д,

(х2,

г/2), . . .,

(хт, ут)

(в'некоторой подходящей системе координат) и

введении

комплексного вектора

состояния

т,

 

 

 

 

 

zk =

xk +

ilJk

>

к = і ,

2, . . .,

 

 

 

 

уравнения движения могут быть записаны в форме

 

 

 

 

 

 

 

dz'

=

г

$ **'**,

 

 

 

(6. 1)

 

 

 

 

dt

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

причем существует квадратичный интеграл движения, представ­ ляемый некоторой положительно определенной эрмитовой формой

E = (Z<gikzk),

(6- 2)

ë k i = S i k

В силу существования квадратичного интеграла коэффициенты взаимодействия (в общем случае комплексные числа) удовлетво­ ряют известным соотношениям

r w * + r /rt< + r fciy = o.

Совершенно очевидно, что далеко не всякая система четного порядка допускает представление в комплексной форме. Это следует из простого подсчета числа параметров (коэффициентов взаимодействия) системы. Так, например, ге=6, 7тг=3, число вещественных параметров системы, допускающей комплексное представление, равно 16 (удвоенное число параметров общего триплета без требования регулярности), в то время как наиболь­

89


шее число независимых параметров СГТ шестого порядка (с уче­ том регулярности системы) оказывается равным 64, т. е. в 4 раза больше.

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

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

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

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

Простейшая нетривиальная система S 2 всегда может быть за­ писана в форме

_____ьг Т

Гdt --

(6. 3)

Эта система имеет интеграл

Е — И-*? + А-

Переходим теперь к ее комплексному расширению, полагая

zi = хі + г'?/і,

Z2 = X2 + ІУі-

Уравнения движения получаем, пользуясь общей формулой

(6. 1)

p - ^ — —kz2Z1,

* =

(6.4)

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

90