ВУЗ: Не указан
Категория: Не указан
Дисциплина: Не указана
Добавлен: 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