Конечномерные модели динамики вихревых течений идеальной жидкости в квадратной области Текст научной статьи по специальности «Математика»

Изв. вузов «ПНД», т. 17, № 6, 2009

УДК 532.5.031

КОНЕЧНОМЕРНЫЕ МОДЕЛИ ДИНАМИКИ ВИХРЕВЫХ ТЕЧЕНИИ ИДЕАЛЬНОЙ ЖИДКОСТИ В КВАДРАТНОЙ ОБЛАСТИ

Н.В. Петровская

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

Ключевые слова: Идеальная несжимаемая жидкость, двумерные вихревые движения, уравнение типа Рауса, метод Галеркина.

Введение

В работах [1,2] предложен новый, нетрадиционный подход к исследованию динамики вихревых движений идеальной несжимаемой жидкости в односвязной области. Кратко изложим схему этого подхода. Известно, что уравнение Эйлера в лагранжевой форме обладает бесконечным набором интегралов циркуляции скорости. Процедурой исключения этих интегралов лагранжевы уравнения гидродинамики приводятся к системе интегро-дифференциальных уравнений, которая вместо интегралов циркуляции имеет бесконечную серию интегралов - якобианов det(dx/da) для каждой жидкой частицы (здесь a, x - «метки» жидкой частицы, a = (ai,a2) - начальное положение жидкой частицы при t = 0, x(a,t) = (xi,x2) - ее положение в момент времени t). Для нашего случая будет интересна рассмотренная в данных работах задача о двумерном течении в односвязной области D, когда уравнение несжимаемости

det(dx/da) = 1

локально может быть разрешено с помощью одновременной замены независимых переменных ai , a2 и функций xi , x2

ai = ai + Фа2, a2 = a2 - Фах, xi = ai - Фа2, x2 = a2 + Фах, ф|4=0 = 0. (1)

159

Здесь а = (ai, a2), a £ D, Ф(а,Ь) - новая неизвестная функция. В силу этой замены уравнение несжимаемости выполняется автоматически для любой гладкой функции Ф(а, t). При t = 0 замена (1) является тождественной. При t = 0 для ее взаимной однозначности требуется, чтобы якобиан

J(a,t) = det ^ д(а а2)^ = 1 + Фа1а1 Фа2а2 — Фаха2 (2)

был отличен от нуля. Для функции Ф(а, t) получено замкнутое интегро-дифферен-циальное уравнение типа Рауса

2Ф* = - J G(x,y)^o(b)\J (|3,t)|d|3 (3)

D

с начальным условием

Ф\*=о = 0. (4)

Здесь b = (b1 ,b2) и y = (y1,y2) связаны с в = (въ |32) (в £ D) соотношениями, аналогичными (1),

b1 = в1 + Фв2 > b2 = в2 - Фв1> yi = в1 - Фв2 > y2 = в2 + Фв1. (5)

Якобиан J(в, t) определяется аналогично (2). В двумерном случае вихрь имеет только одну ненулевую компоненту шз = ш с начальным распределением шо(Ь). G(x, y) - функция Грина краевой задачи

-ДФ = ш, = 0 (6)

для функции тока ф(х) поля скорости жидкости v = (v1, V2).

Автором работ [1,2] высказано предположение, что уравнение (3) может оказаться эффективным средством численного анализа течений идеальной жидкости. Конечно, требование невырожденности замены (1) является очень сильным ограничением, а именно: решения уравнения (3) имеют смысл лишь в достаточно малой окрестности точки t = 0.

В настоящей работе для исследования приближенных решений задачи (3), (4) в квадратной области D используется метод Галеркина в сочетании с методом малого параметра. Полученные конечномерные модели наследуют ряд свойств исходной задачи и, в частности, сохраняют поле вихря, если его начальное распределение является одной из собственных функций задачи

-ДФ = Хф, MdD = 0 (7)

где Х - соответствующее собственное значение. При помощи численных экспериментов изучается эволюция малых возмущений таких решений. Результаты вычислений сравниваются с аналогичными, полученными непосредственным применением метода Галеркина к уравнению Эйлера.

160

1. Конечномерные аппроксимации Галеркина

Рассмотрим двумерные движения идеальной несжимаемой жидкости в квадратной области D: 0 ^ x1 ^ 1, 0 ^ x2 ^ 1. Функция Грина G(x, у) краевой задачи (6) определяется соотношением

G(x,y)

4 i

п42ЕЕ Vij (x)Vij (у) — i=1 j=1

i2 + j2 ’

(8)

где

Vij (x) = sin (nix1) sin (njx2)

собственные функции задачи (7). Подставляя (8) в (3) и переходя к интегрированию по b, приходим к уравнению

Ф*

ЕЕ

i=1 j=1

2

2

л

f+xj2f Vij (у)ю° (b)db-

D

Для его исследования применим метод Галеркина в сочетании с методом малого параметра. Пусть P - некоторое конечное множество пар индексов (m,n). Приближенное решение разыскивается в виде

Ф(а, t) = Е cmn (t)V mn (а). (9)

(m,n)£P

Используя стандартную процедуру метода Галеркина, получим систему обыкновенных дифференциальных уравнений для определения cmn(t)

Cmn

_8_

п2

Е

(i,j)eP

Kn Ai

i2 + j 21 ^i^i3mn’

1

(m, n) <E P,

(10)

где

A

ijmn

Vij(x)Vmn(a) da, Kj = Vij(y)^o(b) db

(11)

DD

Уравнения (10) существенно нелинейны, так как аргументы x и у базисных функций Vij в (11), согласно (1) и (5) зависят от Ф. Для получения более простых конечномерных моделей воспользуемся тем, что уравнение (3) нужно решать при нулевом начальном условии (4). Начальное поле вихря rno(b) возьмем в виде

®o(b)

wpqVpq(b), w°Pq = 4/ mo(b)Vpq(b) db,

(pq)eP D

тогда

Kij wpq Vij (y)Vpq(b)db-

(p,q)eP D

Решение системы (10) должно удовлетворять нулевым начальным условиям

(12)

cmn(0') — °-

(13)

161

Так как при t = 0 замена (1) является тождественной (x = а = a, у = в = b), то

cmn(0) —

wm

2n2(m2 + n2) ’

Поскольку, не теряя общности, можно считать, что (pq)eP(wpq)2 = 0(1), то

cmn(t) = 0(e) на временах порядка е (е ^ 1). Поэтому для вычисления правых частей системы (10) как функций коэффициентов Галеркина cmn(t) можно использовать метод малого параметра, аппроксимируя правые части дифференциальных уравнений отрезками ряда Тейлора. Если в простейшем случае ограничиться линейными относительно неизвестных cmn функциями, то получим следующую систему линейных дифференциальных уравнений с постоянными коэффициентами:

wm

cmn —

2n2(m2 + n2)

+

E

(k,i)eP

Ckl

E

(p,q)eP

w

pqEpqklmni (m,n) £ P-, (14)

где

Epqklmn — 0 8

rpqklmn rmnklpq rpqklmn

m2 + n2

p2 + q2

rijklmn — kjsimkslnj ilskmisjnl y

Slmn = +1 при n = \l— m\, Slmn = -1 при n = l+m, Slmn = 0 при n = \l ± m\.

Отметим, что решение задачи (14), (13), так же как и исходной задачи (3), (4), имеет смысл только при малых t, пока замена (1) остается невырожденной. Чтобы преодолеть это ограничение и исследовать поведение решений на достаточно больших временных промежутках, можно воспользоваться способом, предложенным

В.И. Юдовичем [1]. При заданном начальном распределении завихренности юо(Ь) приближенное решение уравнения (3) ищется на промежутке t £ [0, h]

(h - фиксированная константа) посредством решения задачи Коши (14), (13). После этого определяется поле вихря при t = h и делается новая замена типа (1), где за новый начальный момент времени принимается значение t = h. При этом на каждом шаге разыскивается решение системы (14) с нулевыми начальными данными cmn(0) = 0, но с новым начальным полем вихря. Конечно, константа h должна быть достаточно малой, чтобы исключить обращение в нуль якобиана (2).

Получим формулы для определения поля вихря rn(x, t) при t > 0. Будем разыскивать w(x,t) в виде

rn(x,t) = wij(t)Фij(x). (15)

(i,j)eP

Для определения коэффициентов Wij (t) используется выражение для функции тока

ty(x,t) = J G(x,y(b,t))w0(b)db.

D

Используя (8), (11)-(12) при интегрировании в формуле для ^x, t) и учитывая соотношение ю = — Дф, находим

wij (t) = wij + n2/2 £ w°qY Ckl(t)Tijklpq. (16)

(p,q)eP (k,l)eP

162

Отметим, что при этом закон сохранения вихря

w(x(a, t),t) = wq(a)

(17)

выполняется с точностью до нелинейных относительно переменных Cki слагаемых, отброшенных при выводе уравнений (14). Чтобы определить, насколько сильно это повлияет на свойства приближенных решений, вычисляемых на больших временных промежутках, проведем сравнение приближенных решений задачи, рассчитанных двумя различными способами: на основе линейных уравнений (14), (16) и нелинейных моделей (см. ниже), которые выводятся непосредственным применением метода Галеркина к уравнениям Эйлера движения идеальной несжимаемой жидкости (см., например, [3]). Напомним, что в двумерном случае уравнения Эйлера для функции тока ty(x,t) приводятся к виду

дДф дф дДф дф дДф

dt дх1 дх2 дх2 дх1 ’

Разыскивая функцию тока и вихрь в форме

(18)

Ф(x,t)

пЧ^+i 2) Фу(х), w(x,t)= Vi3 (t)Vj (x)>

(i,j)eP ( j ) (ij)eP

подставляя эти выражения в (18) и требуя ортогональности невязки функциям Vmn(x), выводим уравнения для определения неизвестных vmn(t)

Vmn = - 4 ^2 .2 + .2^2 VklTijklmn, (m, n) Е P. (19)

(i,j)eP г + j k,i

Уравнения нелинейных моделей (19) нужно решать при начальных условиях

vmn(0) = wQ,n.

2. Простейшая модель

Исследование аппроксимаций Галеркина (14), (16) естественно начать с самого простого случая. Простейшая нетривиальная модель получается, если положить

ф = Сифи + С12Ф12 + С21Ф21, «0 = ^°1фц +

Уравнения (14), (16) в этом случае можно привести к виду

c = — D0w0 + L(w0) • c, (20)

w(t) = w0 + [c, w0 ], (21)

где c = (3л2/2)(сп,С12,С21), w = (^11,^12,^21), w0 = (Щ, w^, w°1), а Dq и

163

L(w0) - матрицы

L(w0)

Do = diag(3/8, 3/20, 3/20),

0 — (3/10)w01 (3/10)w102

О? о Ж юо 0 (3/80)w101

V — (3/40)w02 —(3/80)w101 0

Отметим, что след матрицы Sp L(w0) и ее определитель det L(w0) равны нулю для любого начального поля вихря, и система (20) инвариантна относительно вращений

(c, w0) м (Rc, Rw0), где

R

/10 0 \

0 cos у — sin у \ 0 sin y cos y J

Y - произвольный угол. Правые части уравнений (20) - однородные функции коэффициентов w0, и система инвариантна относительно преобразований

(t, w0) м (kt,kw0), k = 0.

Решение задачи Коши с нулевыми начальными данными c(0) = 0 для системы (20) определяется явными формулами. В частности, значениям w0 = к ■ е^, (ei - координатные орты) отвечают решения, линейно зависящие от времени

c = —к ■ D0e; ■ t.

(22)

Однако решения системы (20) имеют смысл только локально, на достаточно малом промежутке времени [0, h], величина которого ограничена требованием невырожденности замены (1). Условие положительности якобиана (2) приводит к оценке h < 8/(3п2). Таким образом, определено отображение

(0, w0) м (c(h), w(h)). (23)

Для решений вида (22) вихрь сохраняется: w(h) = w0. Более того, ввиду инвариантности простейшей модели (20) относительно вращений, сохраняются все распределения вихря вида

ю0 = w0^ sin (nai) sin (na2),

ю0 = w02 sin (na1) sin (2na2) + w01 sin (2na1) sin (na2).

(24)

Далее можно следить за итерациями отображения (23), причем достаточно рассмотреть начальные распределения вихря

ю0 = sin (na1) sin (na2) + e sin (na1) sin (2na2), e> 0. (25)

Для любого начального поля вихря (25) и любого достаточно малого h итерации сходятся к одному из полей вида (24); при этом (w0! монотонно убывает, а

164

(w02)2 + (w21)2 монотонно растет. Рис. 1 иллюстрирует этот процесс для е = 0.02: a - начальное распределение, в - финальная фаза движения. Процесс всегда развивается одинаково: сначала единственный вихрь начинает прецессировать вдоль границы области, затем в одном из углов квадрата появляется второй вихрь (с противоположным знаком завихренности), и прецессирует уже пара вихрей. По мере выравнивания интенсивности вихрей их вращение замедляется, и устанавливается одно из стационарных распределений вида (24).

Что касается уравнений (19), то в рассматриваемом простейшем случае они приводятся к виду

V = [v, D0v], v = (vil, V12) V2l) (26)

и по сути совпадают с уравнениями Эйлера движения свободного твердого тела (причем осесимметричного). Уравнения (26) имеют два независимых первых интеграла: это квадрат нормы вихря ||v||2 = (v, v) и кинетическая энергия

H=1/2(v, D0v). Имеется два семейства равновесий: v= (w01,0, 0) и v= (0, w02, w01); им отвечают поля завихренности (24). Решение задачи Коши с начальными данными (25) определяется формулами

9 9

vii(t) = 1, v 12(f) = еcos —1, V2i(t) = еsin— t,

то есть все движения являются либо стационарными, либо периодическими по времени. Эта ситуация качественно отличается от предыдущей, так как для любого решения уравнений (26) величина v11(t) постоянна. Причина такого различия в том, что в формуле (21) сохранены только линейные относительно Cj слагаемые.

165

Рис. 2. Зависимости ш12 (сплошные линии), полученные посредством итерирования отображения (23) при различных значениях шага h. Начальные данные (25) при е = 0.02. Точки соответствуют значениям функции v12(t) = е cos(9/40) • t

Количественная близость решений на конечном интервале времени может быть обеспечена выбором достаточно малого шага h, что иллюстрируется кривыми рис. 2. Отметим, что обе модели, линейная и нелинейная, определяют одну и ту же частоту колебаний. Амплитуда колебаний в линейном случае растет с увеличением t тем сильнее, чем больше h. Это является следствием пренебрежения нелинейными членами в (16) и приводит к парадоксальному результату: для начального поля ю0,

положительного всюду в области D, с течением времени возникают области с отрицательными значениями завихренности. Однако на конечном промежутке времени разность \wij(t) — Vij(t)| с уменьшением h стремится к нулю. Это иллюстрируется

таблицей 1, где приведены зависимости max |w12(t) I, значения квадрата нормы

0^1^200

вихря ||w||2 = (w, w) и энергии H(w) = 1/2(w,Dow) от h. Для сравнения в последнем столбце даны значения этих же величин при t = 0 (в начальный момент времени эти величины одинаковы для обеих моделей, и на решениях системы (26) они постоянны, так как являются ее первыми интегралами).

Таблица 1

Зависимость характеристик приближенного решения от значений параметра h. Линейная модель (23), начальные данные (25) при е = 0.02, t Е [0, 200]

h 0.02 0.01 0.005 0.0025 t = 0

max W12 (t)\ 0^4^200 0.0319 0.0252 0.0224 0.0212 0.0200

||w||2 (t = 200) 1.000538 1.000452 1.000423 1.000411 1.000400

2H (w) (t = 200) 0.499957 0.500034 0.500060 0.500071 0.500080

3. Компьютерный эксперимент

Уравнения (14), (16) определяют семейство конечномерных моделей, зависящих от выбора множества P. Рассмотрим серию линейных моделей (14), (16), в которых P определяется условием m + n ^ N при N = 4 ... 11 (простейшая модель (20), (21) получается при N = 3). Их решения приходится изучать численно. Для контроля точности получаемых решений удобно использовать нелинейную модель (19) при том же фиксированном P, что и линейная. В обоих случаях для получения решения задачи Коши используется одна и та же модификация метода Рунге-Кутты.

При заданном N можно по предыдущей схеме определить векторы c, w, w0 и записать уравнения (14), (16) в векторной форме

C = —D0w0 + L(w0) • c, (27)

x 0 + 0 О (28)

166

где Do - диагональная матрица с числовыми элементами; L(w0) и M(w0) - матрицы с элементами, однородно и линейно зависящими от w0j, причем M(w0) - кососимметричная матрица для любого вектора w0. Далее определяется отображение вида (23) и изучаются его итерации. Все такие отображения сохраняют начальные распределения вихря

®0(я) = wij Wij (a). (29)

Напомним, что для конечномерных моделей (19), так же как для уравнений Эйлера (18) в квадратной области D, это - стационарные решения.

Цель компьютерных экспериментов - исследование эволюции во времени вихревых движений, в начальный момент времени близких к распределению (29). Начальные поля вихря в численных экспериментах выбирались в форме w0(a) = q>ij(a) + eq>pq(a), где параметр e порядка 10_1 — 10_3. Ниже приведены результаты исследования малых возмущений стационарных движений с полями завихренности 911(a), 912(a) и 921(a).

Найденные приближенные решения сравнивались с решениями задачи Коши с теми же начальными данными для соответствующей нелинейной модели (19). Заметим, что для всех таких моделей квадрат нормы вихря ||v||2 = ^i+j^N и кине-

тическая энергия H(v) = 1/2 ^i+j^N vij/(i2 + j2) являются первыми интегралами, и их удобно использовать для контроля точности приближенных решений.

С увеличением N наблюдается улучшение соответствия между приближенными решениями, полученными разными методами, что иллюстрируется рис. 3. Величина

max \w12(t) — v 12(t)| ~ 2 • 10 4,

0^t^400

относительной погрешности остается в пределах 1% после 80000 итераций отображения. В таблице 2 приведены значения ||w||2 = (w, w) и энергии H(w) = 1/2(w, D0w) при t = 200 в зависимости от N, последний столбец -эти же величины при t = 0 (начальные данные (25) при e = 0.02). Отклонения этих характеристик от начальных значений на фиксированном временном промежутке за-

Таблица 2

Рис. 3. Зависимость ш12 (сплошная линия), полученная посредством итерирования отображения (23) при h = 0.005, начальные данные (25) при е = 0.02. Точки соответствуют значениям функции v12(t) нелинейной модели (19). Для обеих моделей N = 6

Зависимость характеристик приближенного решения от значений параметра N. Линейные модели (27) и (28), начальные данные (25) при e = 0.02

t 200 0

N 3 5 7 9 11

||w||2 1.000423 1.000423 1.000422 1.000422 1.000421 1.000400

2H (w) 0.500060 0.500083 0.500071 0.500074 0.500077 0.500080

167

висят от величины N и h. По-видимому, для ||w||2 более существенной является зависимость от h, а для энергии - от N.

Слабо возмущенное распределение фц (а) эволюционирует одним из двух возможных способов в зависимости от выбора возмущения фрд(а); эти сценарии представлены на рис. 4, 5. Характер движения зависит от свойств четности суммы индексов p + q. Дело в в том, что линейные модели (14), (16) обладают следующим свойством: если Wj (t) = 0 при i + j нечетных, то они равны нулю и при всех t > 0 (аналогично, для нелинейных моделей (19) подпространство, определяемое условиями Vij = 0 при i + j нечетных, является инвариантным).

Рис. 4. Эволюция начального поля вихря фц(а) + £фРЧ (a) при нечетных p + q. Расчет для N = 6, начальное распределение ю0(а) = фц(а) + 0.1 • ф14(а), h = 0.005. Изолинии завихренности после 77000 итераций отображения в последовательные моменты времени t: а - 386, б - 390, в - 394

Рис. 5. Эволюция начального поля вихря фц(а) + £фРЧ(а) при четных p + q. Расчет для N = 8, начальное распределение ю0(а) = фц(а) + 0.1 • ф1з(а), h = 0.005. Изолинии завихренности после 15000 итераций отображения в последовательные моменты времени t: а - 77.17, б - 82.5

168

Если p + q нечетно, то после первой же итерации отображения коэффициенты W12 (t) и W2i(t) становятся отличными от нуля. Поэтому движение имеет тот же характер, что и для трехмерной модели. Вихрь прецессирует вдоль границы области, и такое движение может продолжаться долго - несколько десятков полных оборотов вихря вдоль границы. На рис. 4 приведены изолинии вихря в последовательные моменты времени. Сильная деформация вихря обусловлена сравнительно большой величиной начального возмущения.

Если же p + q четно, то в процессе движения остаются равными нулю все коэффициенты Wij (t) с нечетными i + j и, в частности, w12(t) и w21(t). Движение имеет качественно иной характер: наблюдаются деформации растяжения-сжатия вихря вдоль диагоналей квадрата (рис. 5). Следует еще раз подчеркнуть, что речь идет о свойствах решений на конечном, хотя и большом промежутке времени.

Для двухвихревых режимов движения ф12(а) и ф21(а) наиболее интересные эффекты наблюдаются при малых возмущениях вида e<fpq(a) при p и q нечетных (то есть, при возмущениях с ненулевой средней завихренностью). Наиболее типичный сценарий - почти равномерное вращение пары вихрей вокруг общего центра. Например, на рис. 6 представлены изолинии вихря в последовательные моменты времени, направление вращения пары вихрей - против часовой стрелки.

На рис. 7, а приведены графики w12(t) и w21(t) - кривые A и B, соответственно, это почти гармонические колебания. Кривая C - график величины 5 W (t), W(t) = w‘21(t) + ^3<i+j^N w2j (t) (множитель 5 использован для удобства, в этом режиме движения величина W(t) мала и не превышает 0.02). Отметим, что W(t) достигает минимума при w12(t)w21 (t) = 0 и максимума при \w12(t)\ = \w21 (t)|.

Рис. 6. Изолинии завихренности в последовательные моменты времени t: а - 574.2, б - 636.9, в - 710.1. Начальное распределение юо(а) = ф!2(а) + 0.05фц(а). Расчет для N = 10, t е [0,1600]

169

Другой типичный сценарий - вращательные колебания пары вихрей относительно начального состояния: система переходит из состояния рис. 6, а в состояние рис. 6, в и обратно. На рис. 7, б приведены графики w 12(f), w21(t) и W(t) - кривые A, B и C, соответственно. В этом режиме движения величина W (t) существенно больше, чем в предыдущем случае, а w 12(f) не меняет знак.

Наблюдаются также движения смешанного типа: перемежаются колебания около одного из стационарных состояний ±ф12(а), ±ф21(а) и режим вращения, описанные выше.

Наиболее любопытными представляются движения, в которых наблюдаются колебания пары вихрей относительно одной из диагоналей квадрата. При этом величина \w12(t)\ — \w21(t)\ совершает колебания небольшой амплитуды около нуля, и движение близко к периодическому. На рис. 7, в отражены два таких периода движения: при t Е [10000, 23500] колебания происходят в окрестности состояния wi2(t) = w21(t), а при t Е [25000, 37000] - в окрестности состояния w12(t) = —w21(t). В ряде расчетов наблюдается «застревание» фазовой точки в области \w12(t)\ & \w21(t)\, то есть, возникает близкий к стационарному режим движения с парой почти симметричных относительно диагонали квадрата вихрей. Вообще, создается впечатление, что состояния с \w12(t)\ = \w21(t)\ играют особую роль в динамике системы. Это согласуется с результатами работы [4], где изучались вихревые движения слабо вязкой жидкости в квадрате с условиями периодичности на его границе. В описанных компьютерных экспериментах кусочно-постоянное начальное распределение вихря с близким к нулю средним значением эволюционирует к двухвихревому режиму движения, причем вихри стремятся расположиться по диагонали квадрата.

в

б

Рис. 7. Графики w12(t), w21(t) и W(t) -кривые A, B и C, соответственно: а - ыо(а) = ф12(a) + 0.05фц(а), N = 10; б - ы0(а) = ф12(а) + 0.02фи(a), N = 10; в - ю0(а) = ф12(а) + 0.025ф15(а), N = 8

170

Заключение

В данной работе представлен новый подход к численному исследованию двумерных вихревых движений идеальной жидкости, основанный на приближенном решении интегро-дифференциального уравнения, полученного в работах В.И. Юдо-вича 2000-2005 годов. Метод Галеркина в сочетании с методом малого параметра применен для вывода семейства линейных конечномерных моделей. Для контроля точности приближенных решений использованы нелинейные аппроксимации Галер-кина уравнений Эйлера движения идеальной жидкости в двумерном случае. Компьютерные эксперименты показали, что с ростом размерности моделей улучшается соответствие между приближенными решениями, полученными двумя различными методами. Поэтому можно предположить, что линейные модели высокой размерности могут оказаться эффективным вычислительным средством (количество слагаемых в правой части уравнений для них растет пропорционально квадрату количества уравнений, а для нелинейных моделей - пропорционально кубу). Конечно, проверка этого предположения - предмет отдельного исследования.

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

Автор благодарит М.Ю. Жукова за постановку задачи и А.Б. Моргулиса за полезные обсуждения.

Работа выполнена в рамках европейского научного объединения «Регулярная и хаотическая гидродинамика» (грант РФФИ 07-01-92213 НЦНИЛ) и аналитической ведомственной программы «Развитие научного потенциала высшей школы» (гранты 2.1.1/554 и 2.1.1/6095). Исследование поддержано грантами РФФИ 07-01-00389 и 08-01-00895, а также грантом АФГИР (CRDF) RUM1-2842-RO-06.

Библиографический список

1. В.И. Юдович. Косимметрия и консервативные системы II. Ростовский государственный университет. Ростов-на-Дону, 2000. 24 с. Деп. в ВИНИТИ 06.09.00, № 2772-В00.

2. Yudovich V.I. Topics in an ideal fluid dynamics // Journal of Mathematical Fluid Mechanics. 2005. Vol. 7, Suppl. 3. P. S299.

3. Гледзер Е.Б., Должанский Ф.В., Обухов А.М. Системы гидродинамического типа и их применение. М.: Наука, 1981.

4. Segre E., Kida S. Late states of incompressible 2D decaying vorticity fields // Fluid Dynamics Research. 1998. Vol. 23. P. 89.

Поступила в редакцию 18.09.2008 После доработки 25.06.2009

171

LOW-ORDER DYNAMICAL MODELS FOR VORTICAL FLOWS OF INVISCID FLUID IN SQUARE AREA

N.V. Petrovskaya

The Galerkin method together with the method of small parameter is applied for study of Routh-like equations describing the dynamics of two-dimensional inviscid incompressible fluid flows. A set of simple models for some vortical flows of such fluid in rectangular area has been derived and analysed.

Keywords: Inviscid incompressible fluid, 2D vortical flows, Routh-like equations, Galerkin method.

Петровская Наталья Владимировна - родилась в 1949 году в Ростовена-Дону, окончила Ростовский государственный университет (1971). С 1971 года работает в Ростовском государственном университете (с 2006 года - Южный федеральный университет). Защитила диссертацию на соискание ученой степени кандидата физико-математических наук (1987) по специальности механика жидкости, газа и плазмы в Московском государственном университете. Область научных интересов: компьютерное моделирование течений жидкости, асимптотическое и численное исследование бифуркаций в гидродинамических моделях. Автор 46 научных публикаций.

344090 Ростов-на-Дону, ул. Мильчакова, 8а

Южный федеральный университет, факультет математики, механики и компьютерных наук E-mail: [email protected]

172