Цепочки Гюгонио-Маслова и особые вихревые решения системы уравнений мелкой воды тема автореферата и диссертации по математике, 01.01.03 ВАК РФ

На правах рукописи

Семенов Евгений Сергеевич

Цепочки Гюгонио-Маслова и особые вихревые решения системы уравнений мелкой воды

Специальность 01.01.03 - Маюматическая физика

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико -математических наук

Москва - 2005

Работа выполнена в лаборатории математических методов механики Института проблем механики Российской Академии Наук

Научный руководитель:

Официальные оппоненты:

Ведущая организация:

доктор физико-математических паук, профессор Доброхотов Сергей Юрьевич доктор физико-математических наук, профессор Данилов Владимир Григорьевич доктор физико-математических наук, профессор Радкевич Евгений Владимирович Санкт-Петербургское отделение математического института им. В.А. Стеклова Российской Академии Наук

Зашита диссертации состоится _ &.2005 г. в ^ _

часов в ауд АЛ.___на заседании диссертационного совета К 212.133.01

в Московском государственном институте электроники и математики по адресу: Москва 109028, Б. Трехсвятительский нер., 3/12.

С диссертацией можно ознакомиться в библиотеке Московского государственного института электроники и математики.

Автореферат разослан « ^ ^_» сХ^Ь-^я«^_2005 г.

Ученый секретарь диссертационного совета а^Ч Е. Р. Хакимуллин

ЗЪо8>

Общая характеристика работы

Актуальность темы. В 1980 году В.П. Маслов сформулировал гипотезу1, согласно которой широкий класс квазилинейных систем, включая гидродинамические уравнения, имеет лишь несколько типов решений с особенностями, обладающих следующими свойствами' во-первых, структура особенности сохраняется в течение некоторого интервала времени (свойство структурной самоподобности); во-вторых, структура особенности не меняется при малых возмущениях (свойство структурной устойчивости) К указанным типам относятся ударные волны, "бесконечно узкие" солитоны и особые вихревые решения с особенностью "типа корня квадратного из квадратичной формы." Важными общими характеристиками таких решений являются бесконечные системы обыкновенных дифференциальных уравнений - цепочки Гюгонио-Маслова. В случае ударных волн первое уравнение цепочки совпадает с хорошо известным условием Гюгонио Систематическое изучение цепочек Гюгонио-Маслова для особых вихревых решений было начато сравнительно недавно в работах В.В, Булатова, Ю В. Владимирова, В.Г. Данилова и С.Ю. Доброхотова2. В работах3 были обнаружены весьма любопытные свойства этих цепочек. В указанной работе1 В П. Маслов также выдвинул гипотезу, согласно которой решения последнего типа могут описывать крупномасштабные вихри в атмосфере и океане. В связи с этим интересно изучить в деталях структуру и динамику особых вихревых решений и соответствующих цепочек для некоторых уравнений газо- и гидродинамики.

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

Общая методика исследования основана на использовании асимптотических методов и метода квазиклалсического приближения

'Маслов В П Три алгебры, отвечающие негладким решениям систем квазилинейных гиперболических уравнений// Усп маг наук. - 1980 - 35, № 2 - С 252-253

2БулатопВ В , Владимиров Ю В,ДанилопВ Г , Доброхотов С Ю Пример вычисления "I лаза" тайфуна иа основе гипотезы В П Маслова//Докл РАН. - 1994. - 338, № 1 - С 102-105

3Доброхотов С Ю Редукция к уравнению Хилла цепочки Гюгонио-Маслова для траекторий уединенных вихрей уравнений мелкой поды// Теор и мат физ - 1997. - 112, Л* 1. - С 47-66

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

Научная новизна определяется следующими основными результатами теоретического исследования'

• описана локальная структура особых лихревых решений и выведены цепочки Гюгонио-Маслова для уравнений мелкой воды на /?-плоскости;

• описаны траектории оборванной цепочки Гюгонио-Маслова в предположении малости параметра /3, характеризующим изменение силы Кориолиса с широтой;

• в разумном классе функций установлена единственность структуры особых вихревых решений для уравнений мелкой воды на /?-плоскости;

• найдено многообразие параметров цепочки, в окрестности которого могут существовать разумные с точки зрения приложений медленные траектории в случае, когда ¡3 не является параметром возмущения;

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

Научная и практическая ценность. Работа носит теоретический характер Результаты работы также могу']' иметь практическое применение в проблеме описания и предсказания траекторий крупномасштабных вихрей в атмосфере и океане. С математической точки зрения, работа интересна тем, что проясняет ряд вопросов о динамике особенностей квазилинейных уравнений.

Личное участие автора Результаты диссертации, касающиеся доказательства единственности вида особых вихревых решений для уравнений мелкой воды на /3-плоскости были получены совместно с научным руководителем С.Ю. Доброхотовым и К.В. Панкрашкиным,

описание траектории оборванной цепочки Гюгонио-Маслова, а также описание локальной структуры вихревых особых решений уравнений мелкой воды, учитывающих энергетический обмен, были получены совместно с научным руководителем и профессором Б.Тироцци. Вклад ! автора заключается в проведении части доказательств и выкладок, а также в численном интегрировании оборванных цепочек Гюгонио-Маслова. Вывод уравнения эйконала для фазы особого решения, поправок к условиям Коши-Римана и вывод цепочек Гюгонио-Маслова для уравнений мелкой воды, учитывающих энергетический обмен, был проведен автором самостоятельно.

Положения, выносимые на защиту:

• Детальное описание локальной структуры особых вихревых решений, вывод цепочки Гюгонио-Маслова для уравнений мелкой воды.

• Справедливость уравнения эйконала для фазы в любом порядке тейлоровского разложения.

• Единственность структуры особых вихревых решений для уравнений мелкой воды на /3-плоскости

• Описание траекторий оборванной цепочки Гюгонио-Маслова. в предположении малости параметра /3, характеризующего изменение силы Кориолиса с широтой.

• Медленные траектории оборванной цепочки Гюгонио-Маслова в случае, когда р не является параметром возмущения

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

Апробация работы Основные результаты работы были представлены на семинаре профессора С. Альбеверио в Институте прикладной математики Боннского Университета, на международном семинаре "Дни Диффракции" в 2002, 2003, 2004 годах (Санкт-Петербург), на международном семинаре "Asymptotic Analysis and the Physics of Atmosphere and Ocean" в 2004 году (Италия, Рим), на

международном семинаре "Mathematical aspects of tsunami" в 2005 году (Италия, Рим).

Публикации. Основное содержание работы отражено в 6 публикациях, список которых приведен в конце автореферата.

Структура и объем диссертации. Диссертация состоит из введения и шести глав. Материал диссертации изложен на 107 страницах машинописного текста. Список литературы содержит 55 наименований.

Краткое содержание работы

В введепии проводится обоснование актуальности темы диссертации, содержится постановка целей исследования, формулируются основные результаты работы, а также определяется система уравнений мелкой воды на ^-плоскости, которая используется для моделировании различных физических процессов, включая распространение крупномасштабных вихрей в атмосфере и океане4:

+ ") = 0, + (и, У)и — а/Ги + У?7 = 0. (1)

где х — 1(х1, х-}) 6 1&2, а неизвестными являются двумерный вектор и(ж, — и ((.27, иг (.'Л, ¿)) и функция 77(2г, ~ гео потенциал атмосферы

или высота поверхности океана, Т — ^ ^ 0 ) ' V = '(з^! ¿")>

-- и + /За?2 есть удвоенная частота Кориолиса на ^-плоскости. Наложение построено следующим образом В главе 1 приводится постановка задачи и формулируются результаты, касающиеся структуры особого решения и вида цепочек Гюгонио-Маслова. Ищутся решения системы уравнений (1) вида:

u(ar,i) =u(x,t)+ v,(x,t), u{x,i) - U(x,t)F(S(x, t)),

(2)

Здесь 1 е [0,Т]; /), Я и 5 гладкие скалярные функции, и = й(г>, и>) и II = V2) гладкие двумерные вектор-функции. Предполагается, что функции удовлетворяют одним или нескольким из следующих условий

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

(i) для любого t S(x,t) > 0, и равенство S(x,t) = 0 выполняется в любой момент времени t б [О, Т] в единственной точке х — X(t) = 1(Х1(г), X2(t)). Множество Г — {л; — X{t),t € [О, Г]} назовем траекторией особого решения (2) уравнений (1) на интервале времени [О ,Т].

Функции S, U и R находятся в общем положении, в том смысле, что

(и) Матрица H(t) =

вторых производных

г

= Hess S

невырождена и, как следует из условия (i), положительна на Г.

(iii) Собственные значения матрицы HessS|r различны (показано, что если существует решение (2), такое, что S(x,t) удовлетворяет сформулированному предположению при t = Iq, то это предположение выполняется также для любого t > to).

(iv) разложение U и R по степеням (х — X(t)) начинается с наименее возможных степеней (Это предположение иногда можно опустить, но оно является естественным с точки зреьия возможных приложений и, кроме того, позволяет упростить большинство рассуждений.)

(v) Функция F удовлетворяет следующим условиям, (v-a) Функция F(t) непрерывна при т > О, F(0) — 0; (v-b) F(t) является гладкой при г > 0, linv->+o F'(t) — 00 •

Рассматриваются только те решения системы (1), для которых г) является строго положительной (это следует из ее физического смысла). Пусть gij(t) — i\j\dl+]д{х,t)/дх\дх2\т для любой гладкой функции g(x, t). Положим также uig = и рд = ту|г. Справедлива

Теорема 1. Пусть система (1) vмест решение (2), удовлетворяющее условиям (i)-(v). Тогда имеют место следующие утверждения. 1. Без ограничения общности можно считать, что в (2) F — у/т. 2а. Траектория X(t) вморожена в ноле скоростей и (а т,акже в и):

X(t) = u(X(t),t) = u(X(t), t) Ш V(t) = г(Ш V2(t)) (3)

2b. С точностью до функций, обращающихся в нолъ вместе со всеми своими производными на тра.ект,ории. Г, функция S является решением уравнения эйконала (или Гамильтона-Якоби)

St + (и, VS) = 0. (4)

2с. Комплексные скорости u{x,t) — v(x, t) + iw(x,t) (и Ui(x,i) + гщ(х, t)) на траектории X(t) удовлетворяют условиям Коши—Римана

dv dw de} dv . dw de} , . , ,

= = î(0. = r=P(0 (5)

и

P20\r = №|r + pVi/2, pnlr - W2. (6)

2d. Вторая поправка к условиям Коши-Римана имеет вид

(2Н - tr(H)I) ( Vw + + ~ U tr(H) ( Vw - V02 ~Wn).

\ w20 + W02 — ^cpw J V W02 - W2Q ~ )

(7)

2e. "Потенциальный вихрь" П сохраняется вдоль траектории Г:

m wo - 2р П г = ——— = с = const.

Ро

2f. Функции р uu = Uy/S определяются формулами

где гладкие функции S, tp и вектор-функция Р имеют вид

S = ев&(х - X{t),U(t)Bir(t)(x - X(t))) + 0(\х - X(t)\3), P = VS = Po{t)Xl{t)BU*{t)(x - X{t)) + 0(1® - X(t)]2,

TT/m ( COS в sm в \

11(6') = I ^^ cosö / MamPu4a повоРота na угол

m = + f'p(t)dt = 0O - § /¿rot3u(X(t),t)dt, <p = A + 0Qx-X(t)j), A = const,

B fl+Ъ 0 \ , . .

tS = I p j jl, b, во и А - вещественные константы,

характеризующие начальную структуру вихревого решения, |6| < 1

Интересен тот факт, что возникают условия Копш-Римана (5), а также "поправки" (7) к ним. Эти условия, как правило, не выполняются для произвольных траекторий поля скоростей, и, следовательно, они описывают влияние существования вихря па "гладкий фон" и(х, t).

Цепочка Гюгонио-Маслова для особых решений вида (2) -бесконечная цепочка уравнений, связывающая коэффициенты в разложении функций фона u(x,t) и р(х, t) в ряд Тейлора по ж в окрестности траектории X(t) с учетом соотношений (5) и (6). Для описания траекторий особенностей решений системы (1) цепочка замыкается в духе работы Ф. Прасада и Р. Равиндран5: слагаемые, содержащие коэффициенты полиномов третьего порядка в разложении р и и приравниваются к нулю. Это приводит к системе 17 обыкновенных дифференциальных уравнений, связывающих коэффициенты тейлора в разложении р и и до второго порядка включительно и включающих уравнения X = V, u>q = ßV'i- Цепочка Гюгонио-Маслова приводится ниже (см. формулы (11)-(19)) в новых комплексных переменных.

В главе 2 доказывается пункт 1. Теоремы 1. Схема доказательства выглядит следующим образом. Прежде всего, удобно перейти к подвижной системе координат (x',t), где х' = х — X(t), и положим и' = и — V. Тогда система (1) записывается в виде (мы опускаем штрихи на новых переменных х)

rit+ (V,Tju) = 0, U( + (u,V)u +Vr? +V'-wT(u-f-y) = 0. (8)

После подстановки решения вида (2) в уравнения (8), получаем уравнение вида

F'{T)A + F(T)B + F(T)F'(T)C + F\r)D + E = 0.

с гладкими векторными коэффициентами А, В, С, D и Е. Здесь г = S(x,t). Анализ этого уравнения приводит к следующему утверждению Лемма (о модельных уравнениях). Пусть выполнены условия (ii)- (iii), тогда справедливы следующие утверждения А. Существуют гладкие функции а, ß, 7, S, зависящие от (х, t) и такие, что F удовлетворяет уравнению Риккати

а{х, t)F'{r) + ß(x, t)F2(r) + 7(®, t)F(r) + 5{х, t) = 0,

где производные по х функции а не равны нулю при х = 0.

В. Функция F удовлетворяет одному из следующих трех "модельных" уравнений:

'Ravindran R., Prasad Р A new theory of shock dynamics. Part I (П)// Appl Math Lett - 1490 -3, Л» 3. - С 77-79

- квадратному уравнению: a(x,t)F2(r) + b{x,t)F{r) + c(x,t) = О,

линейному дифференциальному уравнению: a(x,t)F'(r) + b(x, t)F(r) + с(х, t) = О,

- или кубическому уравнению: а(х, t)F^(r)+b(x, t)F2(T)+c(x, t)F(r) + d(x,t) — 0. Здесь a, b, с, d также гладкие функции, некоторые производные которых по х-ие равны нулю при х — О.

Анализ модельных уравнений приводит к следующему утверждению. Лемма. Без ограничения общности можно считать, что существует всего четыре возмоо/сные функции F в (2): (Fl) F = тк, (F2) F = г log г, (F3) F = т1'3 + ат2'3+п, (F4) F = г2/3 + ат^+п, где 0 < к < 1, о" = ±1 или 0, п > 0 целое число. Квадратное уравнение имеет лишь решение т,ипа (F1) с к = 1/2; линейное дифференциальное уравнение имеет решения типа (F1) и (F2); кубическое уравнение имеет решения (F1) с к = 1/2,'(F3) и (F4).

Далее решения вида (2) с функциями F из (F1)-(F4) подставляются в уравнения (8). При к ф 1/2 получившиеся системы уравнений оказываются несовместны. Таким образом справедливо следующее утверждение.

Лемма. Решений с особенностью вида (Fl)-(F4) при к ф 1/2 не существует

Из перечисленных лемм следует утверждение пункта 1. Теоремы 1.

В главе 3 изучается негладкая часть решения и доказываются пункты 2а. -2f. Теоремы 1. После подстановки решения вида (2) с F = у/т в (8) уравнения для (и, р) и (й, р) распадаются следующим образом:

| +V.(fm + pu)= О, ■

§ +(u,V)u + Vp+V~wT(u + V) + (u,V)u = 0;

I + 0, f .

Щ +(и, V)u + (и, V)u + Vp - cjTU = 0. v '

Показывается, что хотя система (10) линейная относительно и и р, нельзя напрямую применить методы разложений по гладкости или лучевых разложений6 для задачи Коши p\t=a — R(x,0)y/S(x, 0),u\t=o = U(x, 0) -\fS(x, 0). По существу, используются идеи, восходящие к

6Вабич В М , Булдырев В С., Малоткоя И А Прогтрапстпегато-временной лучевой метод -Ленинград: ЛГУ, 1985

этим методам, однако имеются некоторые существенные отличия, существенные с точки зрения лучевого метода. Прежде всего, использование разложения по гладкости не гарантирует структуры (2) решения, даже если эта структура имеет место при Ь = 0. Кроме того возникает ситуация, известной п теории гиперболических уравнений как "негладкая смена кратности характеристик": собственные частоты V = 54 + А± = + (и,У)5| ± р|У5|, соответствующие

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

рассматривается вся система (9)—(10) и изучитаются ее свойства. С первого взгляда может показаться, что система (9)—(10) недоопределена' имеются шесть уравнений относительно семи неизвестных функций v, ги, р и Vь и2, Я, 5. Показывается, что заданные свойства функции 5 и гладкость всех названных функций фактически приводит к возникновению уравнения эйконала для фазы 5 и система замыкается.

Процедура аналитических вычислений состоит в разложении всех функций в ряды Тейлора, подстановки этих рядов в систему (10) и приравнивании коэффициентов при соответствующих степенях х. Указанная процедура приводит к некоторым равенствам и уравнениям для коэффициентов в рядах Тейлора и т.д. При этом

прослеживается прямая аналогия с методом ВКБ: для фазы выводится уравнение эйконала (4), а условия Коши-Римана (5) получаются из уравнения, аналогичному уравнению переноса. Показывается также, что условия Коши-Римана можно получить из условия сохранения потенциального вихря. Также, выводятся поправки к условиям Коши-Римана следующего порядка (7), которые интересны тем, что в уравнения (7), дающие условия на гладкий фон входят параметры особой части решения (2).

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

г = XI - ХхЦ) + г(ж2 - Х3(<)), г = х1 - Х^) - г(х2 - Х3($).

После введения новых комплексных переменных X' = Х\ + гХг, V' = У\ + ¿^2, У, и вещественных переменных /х и А (эти переменные

выражаются через коэффициенты в разложении функций р и и в ряд Тейлора по г, I в окрестности траектории Х{£) до второго порядка

включительно) цепочка Ггогонио-Маслова принимает вид

X = V, (И)

V + iuQV+ W-2Z) = Q, (12)

Y = i(p-Uo)Y-.^^ + fl, (13) ¿ = 1(зр-Шо)г + Щ£ + /2, (14)

W = -фЖ, (15)

¿г-г(р + Ш0)С/ + /з, (16)

Л = 2ßn3(f ImY - ^ImZ - f/mW), (17)

w0 = ßlmV = (18)

gf + 6V = + 2|c\Re((Z - f) W + §ЯУ))- (19)

Здесь слагаемые fi, h и /з содержат коэффициенты Тейлора 3-го порядка. Обрыв цепочки (замыкание) означает, что /ь /г и /з приравниваются к нулю.

В главе 5 изучается оборванная цепочка (11)—(19) в трех случаях: ß = О, ß - мало, и случай, когда ß нельзя считать параметром возмущения. В работах С.Ю. Доброхотова7 было показано, что оборванная цепочка (11) -(19) точно интегрируется и сводится к уравнению Хилла в случае ß = О

4* + 0(ф)* = О-

Это сведение позволяет выявить ряд свойств траекторий Г связав их со свойствами уравнения Хилла (такими, как свойства устойчивости, рассмотрением зон и лакун в спектре и т.д.). Физически разумное предположение о достаточной гладкости траекторий Г и вполне элементарное исследование уравнений для скорости и функций У, Z. W приводят к требованию малости малости одной из следующих частот: vi = (p~cjо) ) мала, тогда Z ~ О, W а 0 (случай I); ¡/2 ~ (uio - Зр) ) мала, тогда Y ~ О, W cz 0 (случай II); Уз ~ —р) мала, тогда Y ~ О, Z ~ 0 (случай III). Это требование должно быть выполнено как в случае ß — 0, так и в случае малых ß.

7 Доброхотов С. Ю. Редукция к уравнению Хилла цепочки Гюгопио-Маслова для траекторий уединен ных вихрей уравнений мелкой воды// Теор. и мат физ. - 1997. - 112, Л* 1. - С- 47-66

В случае /3 — 0 у системы (11) (19) легко получаются траектории в виде окружности, если положить в случаях I, II, III последовательно (2 = о, IV = 0), (У = 0, IV = 0), (г = 0, У = 0). Если считать параметр /3 малым, то слагаемые с /3 в (11)—(19) играют роль адиабатического возмущения, которое, в частности, "деформирует" указанные траектории. Одним из эффектов влияния этих слагаемых состоит в том, что' движение по "деформированным" окружностям может менять направление и появляются траектории Г в виде зигзагов. Такие "деформированные" траектории могут быть описаны с помощью гамилътоновых систем с одиой степенью свободы. Для всех случаев I, II и III они имеют одинаковую структуру

¿0 = /^г, Ф - уо1<=о = ^о(О), ф\{=о = Ф(0),

Н = сое фМ + N,

но с разными функциями /, М и N. Здесь гамильтоновыми переменными являются частота шо и угол ф направления скорости V Например, для случая, когда мала частота Зр — и/о,

N = 7М3/2 г7 У^У1 - 6) .

4 к у/(уг + 3)(21/7 - I)3

f - ^У^О^У - I)3 _ WQ 7Г1у/Г] + 3 ' U>2' где Mj и ui2 константы интегрирования Ее фазовые портреты на плоскости (шо, Ф) определяются интегралом системы Н.Н— const, и для отдельных областей почти совпадают с хорошо известными фазовыми портретами траекторий физического маятника (см. Рис. 1).

В каждом из 3-х случаев, траектория X(t) зависит от шести параметров М„ с, -фу,(так как шо = О + flX^)- Эти шесть "существенных" параметров могут быть восстановлены при помотци полученных в диссертационной работе формул через гнесть имеюгцих физический смысл величин' начальные положение (Xj, X®) и скорость

Рис. 1: Пример фазового портрета гамильтоповой системы (а) и траектории центра вихря на /^-плоскости в случае когда мала частота 3р — шо (Ь). Траектория, отмеченная жирной линией, соответствует траектории гамильтоповой системы с одной стспспыо свободы, траектория, отмеченная тонкой линией, - решению 17 уравнений оборванной цепочки Гюгопио-Маслопа до усреднения.

(V,0, У§) центра вихря, потенциальный вихрь с и третья компонента 1Ыз и в центре вихря.

При выводе указанных в конце предыдущего пункта гамильтоновых систем существенным является предположение об относительной малости слагаемых, пропорциональных /3 ("вынуждающих сил") в правых частях уравнений для У и 2. Если эти слагаемые не малы, то ясно, что их нельзя рассматривать как адиабатическое возмущение, и более того, они могут играть существенную роль в описании траекторий Г Имея ввиду явные выражения для правых частей в системе (11)—(19), предполагается, что для интересующих нас решений в первом приближении всеми производными по t можно пренебречь, что приводит к уравнению для определения критических точек системы. Элементарный анализ показывает, что эти критические точки образуют два двумерных (пересекающихся) гладких "критических" многообразия, характеризуемых условиями (а) — 0 и (Ь) р = 0. Таким многообразиям отвечают траектории с равномерным движением, параллельным оси х\ (т.е. движениям вдоль параллели), причем случай (Ь) не дает подходящих для описания мезомасштабных вихрей траекторий. В первом приближении движение определяется линеаризованными уравнениями в некоторой окрестности

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

Л? + 22а;02Л^ + (105 + ^)^ + (148+^)^ + (64-81<5-^)а;08 = О,

аЪ

где безразмерный параметр о =

Медленные движений в окрестности критического многообразия возможны тогда, когда существуют "маленькие" корни Л этого полинома Одно из собственных значений обращается в ноль тогда и только тогда, когда слагаемое (64 — 815 — в последнем уравнении обращается

в ноль. Это даст 6 = ~ 0,465. Например, в случае, когда

¡3 плоскость касается земли в точке, расположенной на широте 20° (а; и 0,18 1/час) такая ситуация реализуется, если принять значение геопотенциала на траектории Г соответствующее высоте слоя атмосферы ~ 10км (77|г ~ 860000км2/час2). Иными словами, в окрестности широты 20° возможно наличие медленных траекторий, если предположить, что воображаемая высота слоя атмосферы равна 10 км. Заметим, что в задачах о мезомасштабных вихрях в атмосфере этой величине часто приписываются значения 7 Ч-12 км. 1 В главе 6 описывается цепочка Гюгонио-Маглова для более; сложной

I, системы уравнений, описывающих динамику двумерных вихрей в

* атмосфере с учетом энергетического обмена между атмосферой и

океаном8: После обезразмеривашш эта система выглядит следующим образом.

д

— + и • Уи - шТи + Чг) + (?? + 6)Чв - 0. (20)

^ + V ■ (И + <5го13и - Г)к{е - в) = 0. (21)

~ + и • V© + к{е - 0) = 0. (22)

В отличие от уравнений мелкой воды, добавляется еще одно уравнение для скалярной функции 0(аг, £) (плотность или температура), юЬ^и = д\12]дх\ - дщ/дхг, <5, к - физические константы, характеризующие приток массы из пограничного слоя атмосферы (слоя Экмана) и нагрев вихря, 0 соответствует значению температуры поверхности океана.

8Должанский Ф В., Кримов Б А , Мании Д Ю Устойчивость и вихревые структуры квазидвумерных сдвиговых течений// Усп. физ нлук - 1990 - 160, № 7. - С. 1-47

Показывается, что для особых решений типа корня квадратного из квадратичной формы выполнены условия вмороженности и уравнение Гамильтона-Якоби вида (4) для фазы S. Условия Котпи-Римана (5) тоже сохраняются, если потребовать дополнительные условия на гладкую составляющую 6(x,t) функции 0(х, ¿), а именно V9(x, i)|r - О Численный анализ соответствующей цепочки Гюгонио-Маслова показал, что у оборванной цепочки существует семейство траекторий близких по параметрам к описанному в предыдущей главе режиму, но при этом центр особенности может двигаться в меридиональном направлении.

Работы, опубликованные по теме диссертации

1 Доброхотов С, Ю , Панкрашкин К. В., Семенов Е. С. О гипотезе Маслова о структуре слабых точеч иых особенностей уравнений мелкой воды// Докл. РАН. - 2001. - 379, № 2 - С. 173-176

2. Dobrokhotov S Yu., Pankrashkiri К. V., Semenov Е S. Proof of Maslov's conjecture about the structure of weak point singular solutions of the shallow water equations// Russ. J. Math. Phys. - 2001. - 8, № I. - C. 25-54

3. Семенов E. С. Об условиях Гюгонио-Маслова для вихревых особых решений системы уравнений мелкой воды// Мат. заметки. - 2002. - 71, № 6. - С. 902-913

4. Доброхотов С Ю., Семенов Е. С., Тироцци Б. Цепочки Гюгонио-Маслова для сингулярных вихревых решений квазилинейных гиперболических систем и траектории тайфунов// Современная математика. Фундаментальные направления. - 2003. - 2. - С 5-44

5. Dobrokhotov S. Yu., Pankrashkin К. V., Semenov Е. S. On Maslov's conjecture about square root type singular solutions of the shallow water equations, In: Nonlineai Algebraic Analysis and Applications, edited by A Delcroix ct al. - Cambridge Sci. Publ., 2004

6. Доброхотов С. К)., Семенов Е. С,, Тироцци Б. Вычисление интегралов цепочки Гюгонио-Маслова для особых вихревых решений уравнений мелкой воды// Теоретичесая и математическая физика. -2004. - 139, № 1. - С. 500-512

i

Цепочки Гюгонио-Маслова и особые вихревые решения системы уравнений мелкой воды

Семенов Евгений Сергеевич Подписано к печати 20.07.05. Заказ №6-2005. Тираж 70 экз.

Отпечатано на ризографе, ИПМсх РАН 119526 Москва, проспект Вернадского, д. 101, к. 1.

»Î3824

РНБ Русский фонд

2006-4 9308

«i i

Введение

1 Цепочки Гюгонио-Маслова и особенности типа квадратного корня

1.1 Цепочки Гюгонио-Маслова как необходимые условия существования особых вихревых решений.

1.2 Свойства сингулярной вихревой составляющей решения

1.3 Замыкание цепочки.

2 Единственность особенности типа квадратного корня

11 2.1 Вспомогательные утверждения.

2.2 Уравнение для функции F с особенностью в морсовых координатах

2.3 Модельные уравнения.

2.4 Возможные решения модельных уравнений.

2.5 Исходная система уравнений и особенности перечисленных типов.

3 Анализ негладкой компоненты решения и возникновение условий Коши-Римана на траектории

3.1 Негладкая (вихревая) компонента решения.

3.2 Особенность, "вмороженная" в поле скоростей

3.3 Уравнение эйконала для функции S

3.4 Уравнение переноса для амплитуды и возникновение условий Коши-Римана.

3.5 Вычисление поправки U^.

3.6 Порядок остатка для негладкой составляющей решения

3.7 Закон сохранения для потенциального вихря и условия Коши-Римана. $ 3.8 Поправка к условиям Коши-Римана.

4 Анализ гладкой составляющей решения

4.1 Вывод цепочки Гюгонио-Маслова для гладкой составляющей, завершение доказательства теоремы

4.2 Цепочки Гюгонио-Маслова в новых комплексных переменных.

5 Интегрируемость оборванной цепочки: редукция к уравнению хилла и одномерные гамильтоновы системы

5.1 Новые зависимые переменные и интегралы оборванной цепочки.

5.2 Уравнение Ермакова и редукция к уравнению Хилла

5.3 О влиянии устойчивости уравнения Хилла и наличия постоянной силы Кориолиса на траектории.

5.4 Критические режимы.

5.5 Критические режимы в случае /3 = 0.

5.6 Влияние /б—эффекта и гамильтоновы системы в критических режимах.

5.7 Медленные траектории оборванной цепочки в случае, когда /3 не является параметром возмущения.

6 Цепочки Гюгонио-Маслова для системы уравнений мелкой воды с учетом энергетического обмена

В 1980 году В.П. Маслов [28, 29] сформулировал гипотезу, согласно которой широкий класс квазилинейных гиперболических систем, включая гидродинамические уравнения, допускает лишь несколько типов решений с особенностями, обладающих следующими свойствами. Во-первых, структура особенности сохраняется в течение некоторого интервала времени (свойство структурной самоподобности); во-вторых, структура особенности не меняется при малых возмущениях (свойство структурной устойчивости). Такой выбор свойств соответствует наличию в уравнениях нелинейности: в случае линейных гиперболических систем структура любой особенности в начальных данных сохраняется для решения (по крайней мере в течение малого интервала времени). К указанным типам особых решений принадлежат ударные волны, "бесконечно узкие" солитоны и вихревые особые решения "типа корня квадратного из квадратичной формы." Такие решения могут быть описаны формулой, аналогичной "нелинейным" (уиземовским) решениям и искаженным волнам Римана (см. [36]) w = f(x,t) + g(x,t )F(S(x,t)), (1) где w - векторная (или скалярная) функция, х € Kn, F{t) - скалярная функция, гладкая вне т = 0 и имеющая особенность при г = 0, а фаза S(x,t), векторный (или скалярный) фон f(x,t) и амплитуда д(х, t) - гладкие функции. Особенность может соответствовать, например, разрыву первого рода (тогда мы имеем дело с ударными волнами), а может принадлежать классу С1. Очевидно, особенности w(x,t) определяются нулями X(t) функции S(x,t). Например, для ударных волн в одномерном (n = 1) случае имеем F = ©(г), где ©(г) функция Хевисайда (0 = 0 при г < 0 и 1 при т > 0), S = х — X (t). Для другого типа особенности по-прежнему S = х — X(t), но F = Sol (г), где Sol (г) = 0 при т = 0 и 1 при т = 0. В этом случае функция w описывает бесконечно узкий солитон на фоне и(х, t). Решения такого вида возникают, например, как предельные решения уравнения Кортевега—де Фриза с дисперсией е2 при е —» 0 (см. [14, 30]).

Рассмотрим еще один пример в двумерном случае, который будет служить основным примером в этой работе. Пусть F = гг, (0 < г < 1), и для каждого t неотрицательная функция S обращается в нуль в единственной точке *(xi = Xi (t), х2 = X2(i)). Тогда в общем случае функция S с точностью до слагаемых более высокого порядка представляет собой положительную квадратичную форму по переменным х с различными собственными значениями. Нулем функции S является точка х = X(t) = *(Хх(4),Х2(г)), двигающаяся вдоль траектории Г = {х = '(.X^), X2(f))} (задающей особенность). В данном случае мы имеем слабую "точечную" особенность: сама функция w непрерывна и даже обращается в нуль в особой точке, но некоторые производные функции w не являются непрерывными. Конечно, можно рассматривать и другие функции F(t) с особенностями, обладающие похожими свойствами, например г log г, (const + log г)-1, а также различные их линейные комбинации.

Свойства "самоподобности и устойчивости" структуры означают следующее. Пусть решения имеют при некотором <0 ВИД (1) с заданной функцией F(t). Тогда, во-первых, решение сохраняет вид (1) с той же функцией F(t) и при t > to, по крайней мере при достаточно малых t — to; во-вторых, малые изменения начальных данных S(x,to), /(ж, t0), g(x,t0) и коэффициентов исходного уравнения не влияет на структуру особенности функции w, определяемой функцией F. Гипотеза Маслова заключалась в том, что для многих квазилинейных гиперболических уравнений, имеющих физический смысл, почти все возможные особенности с указанными свойствами имеют описанный в приведенных выше примерах вид. Более того, в последнем примере возможна особенность только с F = у/т (см. [29, 22, 14]). И хотя соответствующее уравнение может иметь частные решения (например, радиально-симметричные, т.е. такие, для которых S ~ xl+xQ, отличные от описанных выше, эти особые решения, по-видимому, пропадают при малых возмущениях. Доказательство этой гипотезы совсем не тривиально; в настоящей работе оно приведено для слабой точечной особенности системы уравнений мелкой воды. (По поводу другого типа особенностей см. [14].) Приведем лишь некоторые алгебраические соображения, обосновывающие существование особенностей упомянутых типов.

Поскольку исходная система нелинейна, необходимо определить произведения и степени компонентов функции (1). Функции вида f(x,t) + g(x,t)F(S(x,t)), где F = ©(г), F = Sol(r) и F = у/т, обладают весьма важным свойством: каждая из них порождает алгебру функций с особенностью над кольцом гладких функций, имеющую только две образующих: 1 и F. Это свойство справедливо в силу того, что F2 = F в первых двух случаях и F2(S) = G{x,t) есть гладкая функция в третьем случае. Для того чтобы изучить свойства решения (1), его следует подставить в исходную систему, сгруппировать слагаемые при различных типах особенностей, а затем и регулярные (гладкие) слагаемые и приравнять их к нулю. При этом исходная система распадается на две системы с функциями S, /, д (которые, вообще говоря, подлежат определению). Если рассматривать решения (1), основываясь на алгебрах с более чем двумя генераторами (например, если F = г1/3, образующие суть 1,т1/3,т2/3), то исходная система будет распадаться на три системы, являющимися переопределенными, и т. д. Это лишь эвристическое объяснение: в случае бесконечного числа образующих не очевидно, как можно корректно обосновать группировку слагаемых, соответствующих различным особенностям, и дальнейшее разложение исходной системы.

Цепочки Гюгонио—Маслова и их замыкание

Другое очень важное наблюдение, сделанное в [29], заключалось в том, что, независимо от физической природы рассматриваемого явления, решения с особенностью (1) и их описания очень близки с математической точки зрения. Одной из таких общих характеристик являются новые бесконечные системы обыкновенных дифференциальных уравнений (цепочки), определяющие динамику решений. Хотя решения (1) имеют особенности, они определяются гладкими скалярными и векторными функциями S(x,t), f(x,t), g(x,t). Следовательно, с этими функциями, а также и с самими решениями (1), можно связать как множество коэффициентов их рядов Тейлора в (особых) точках х = X(t), так и сами траектории х = X(t). Указанные цепочки возникают как необходимые условия существования решения (1). В случае ударных волн первое уравнение в такой цепочке совпадает с хорошо известным условием Гюгонио, а последующие уравнения могут рассматриваться как поправки к первому уравнению. Несмотря на то что цепочки для особенностей других типов могут иметь иной вид и, вероятно, имеют другой физический смысл, в [29] было установлено, что их математическая сущность оказывается одинаковой. Такие цепочки называются цепочками Гюгонио—Маслова.

На цепочки Гюгонио—Маслова в нелинейных задачах можно также посмотреть с другой точки зрения. Решения вида (1) являются обобщенными решениями, и, следовательно, методы исследования таких решений связаны с построением алгебр распределений и их приложениями в теории нелинейных уравнений (см. [39, 21], а также [14], где эти уравнения рассматриваются достаточно подробно). Понятно, что структура обобщенных решений может зависеть от выбора соответствующей обобщенной функции. При этом оказывается, что одних лишь математических соображений недостаточно для того, чтобы указанный выбор был единственным. Обычно введение того или иного распределения в нелинейных задачах основывается на предшествующей регуляризации обобщенной функции, зависящей от некоторого малого или большого параметра. Вообще говоря, конечный результат зависит также и от используемой регуляризации, что "нехорошо" с точки зрения возможных физических приложений. С другой стороны, в физических задачах гиперболические системы не существуют сами по себе; как правило, они возникают в качестве некоторого приближения уравнений с малой вязкостью или дисперсией. Таким образом, гиперболические системы следует рассматривать как пределы таких "прообразов", их решения должны быть пределами решений "прообразов", и решения "предельной" гиперболической системы (так же как и соответствующий малый параметр) должны быть согласованы с исходными уравнениями. Грубо говоря, это всегда (и автоматически) выполняется для регулярных решений гиперболических систем. Однако выяснение того факта, выполняются ли указанные условия для решений с особенностями, представляет собой серьезную проблему. Часто оказывается весьма сложно или даже невозможно изучить описанное соответствие. Более того, во многих реальных физических ситуациях соответствующий «прообраз» с вязкостью или дисперсией принимается не всеми специалистами в качестве подходящей модели. Тем не менее, существуют некоторые характеристики обобщенных решений, инвариантные относительно выбора способа регуляризации обобщенных функций. Цепочки Гюгонио—Маслова как раз и являются такого рода характеристиками, поскольку они возникают как необходимые условия.

Цепочки для ударных волн и решений уравнений газодинамики были изучены много лет назад (см. [28, 30, 42], а также [52, 13]). Систематическое изучение цепочек Гюгонио-Маслова для особых вихревых решений было начато сравнительно недавно в работах [22, 10].

Цепочки Гюгонио—Маслова незамкнуты, поскольку первые N уравнений содержат более N неизвестных. Поэтому положение особенности, вообще говоря, нельзя однозначно определить из таких систем. Замыкание такой системы возможно при помощи использования глобальных свойств решений или некоторых дополнительных предположений [28, 30] (см. также [47]), например, предположения о малости негладкой амплитуды. Незамкнутые цепочки уравнений возникают в различных задачах статистической физики и механики. Хорошо известны цепочки Боголюбова—Борна— Грина—Кирквуда—Ивона, цепочки моментов, возникающие в статистической гидромеханике (см., например, [11, 25]). Проблема замыкания этих цепочек является одной из наиболее интересных и трудных в математической физике. Напомним, что после замыкания ББГКИ-цепочек получается кинетическое уравнение Больцмана.

Для цепочек типа Гюгонио, возникающих при описании ударных волн для простейшего нелинейного уравнения (уравнения Хопфа для простой волны) в работе [47] был предложен метод замыкания, основанный на приравнивании лишних компонент решения (с большими номерами) к нулю. Похожая процедура замыкания использовалась в статистической механике и гидродинамике (см., например, [25, 11]). Цель настоящей работы состоит в использовании цепочек типа Гюгонио, которые обрываются в духе работы [47], для описания траекторий вихревых особенностей решений системы (2) и в применении этого подхода к исследованию задачи о траекториях мезомасштабных вихрей в атмосфере и океане.

Ясно, что системы уравнений, полученные обрывом цепочек Гюгонио-Маслова, не могут адекватно описывать распространение особенности на больших интервалах времени, так как обрыв цепочки неизбежно влечет за собой локализацию задачи: эволюция особенности определяется только тем, что происходило вблизи данного момента времени. Однако если рассматриваются не слишком большие интервалы времени и если оборванная цепочка обладает некоторым свойством устойчивости, то использование такой оборванной цепочки представляется вполне оправданным. С другой стороны, например, в задаче о траектории глаза тайфуна невозможно достаточно достоверно измерить поле скорости и и геопотенциала 7? в начальный момент времени для того, чтобы получить корректную задачу Коши для соответствующей системы уравнений в частных производных. Однако траектория глаза тайфуна может быть при этом измерена достаточно точно, например при помощи спутников, и последующая траектория может быть предсказана на основании известной предыстории посредством решения задачи экстраполяции. С этой целью можно в свою очередь пользоваться формулами для траекторий вихря, полученными интегрированием (оборванных) цепочек. Таким образом, локализация оказывается вполне обоснованной с точки зрения физической постановки задачи.

Уравнение мелкой воды и мезомасштабные вихри

В статье [29] содержится также идея, согласно которой решения (1) могут применяться для описания некоторых природных явлений: волн цунами в океане (случай «узких ссшитонов») и мезомасштабных (или «крупных») вихрей в атмосфере и океане (тайфунов, ураганов, рингов и т.д., что соответствует случаю особенности типа квадратного корня). Таким образом, оборванные "вихревые" цепочки Гюгонио-Маслова могут с определенной точностью описывать динамику мезомасштабных вихрей.

Система уравнений мелкой воды также хорошо известна как простейший пример двумерного бездисперсионного приближения с нулевой вязкостью в моделировании различных эволюционных физических процессов, включая распространение мезомасштабных вихрей в атмосфере и океане [20, 34, 18]. Такая система с переменной силой Кориолиса в так называемом /^-плоскостном приближении имеет вид д дух + (V,77u) = 0, -^ + (u,V)u-u/Tu + V77 = 0. (2) где х = *(х1,х2) Е R2, а неизвестными являются двумерный вектор u(x, t) = f(ui(x, t), и2(ж, t)) и функция r](x,t) - геопотенциал атмосферы (или уровень свободной поверхности в теории волн на воде), Т = -1 0 ) ' V = А)' и = й + @х2 есть удвоенная частота

Кориолиса на ^-плоскости, наконец и, /3 суть параметры (физические константы), Р предполагается достаточно малой. Обозначим также через (,) скалярное произведение; индекс t слева вверху будет обозначать транспонирование матриц и векторов.

Система (2) обладает рядом важных свойств, таких как наличие закона сохранения, возможность гамильтонова представления и т. д. Одним из таких свойств является существование лагранжева инварианта — так называемого потенциального вихря или инварианта Россби п U2X1 ~ Ц1Д2 + Ц /дЧ V

Свойство инвариантности означает, что П остается неизменным вдоль траекторий поля вектора скорости и: П4+ (и, \7)П = 0.

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

Приведем теперь соображения Маслова относительно связи слабых особенностей уравнений мелкой воды (2) с моделями движения мезомаспггабных вихрей. Если мы предположим, что решения системы (2) описывают крупномасштабные явления в атмосфере или океане (например распространение тайфунов, ураганов, рингов) и мезомасштабные вихри отвечают решениям, имеющим в некоторой точке слабую точечную особенность и обладающим свойствами структурной устойчивости и самоподобности, то в силу единственности структуры такой особенности, траектория центра вихря («глаза» тайфуна) должна быть близка к траектории центра слабой точечной особенности типа квадратного корня. Отметим еще, что траектория центра вихря и геопотенциал rj, соответствующие природному мезомасштабному вихрю, должны обладать следующими свойствами, разумными с физической точки зрения. Траектории не имеют слишком больших «дрожаний и петель», а г] есть положительная функция, достаточно медленно изменяющаяся вдоль траектории. Эта дополнительная информация играет важную роль при изучении решений (1) с точки зрения их пригодности к описанию динамики мезомасштабных вихрей.

Интегрируемость оборванных вихревых цепочек Гюгонио—Маслова для уравнений мелкой воды

В настоящей работе показано, что после ряда преобразований оборванная цепочка для вихревых решений системы (2) сводится к сложной системе, содержащей 17 обыкновенных дифференциальных уравнений. Интересным и довольно неожиданным оказывается тот факт, что интегрирование этих 17 уравнений в случае постоянных частот Кориолиса (/3 = 0) приводит в точности к линейному уравнению Хилла, хорошо известному в теории колебаний, небесной механике, теории солитонов и т.д.

Наличие /^-эффекта приводит к появлению дополнительных слагаемых порядка /J в соответствующей системе 17-ти обыкновенных дифференциальных уравнений. Если считать /3 малым параметром, то наблюдается медленная («адиабатическая») эволюция постоянных интегрирования в невозмущенной оборванной цепочке. Отметим еще одно интересное и неожиданное явление: при наличии упомянутых выше разумных предположений о траектории центра вихря и геопотенциале усредненные уравнения оказываются гамильтоновыми системами с одной степенью свободы, интегрируемыми в квадратурах и аналогичными уравнению физического маятника. В случае, когда параметр нельзя считать параметром возмущения, у указанной системы 17-ти обыкновенных дифференциальных уравнений удается найти семейство двумерных многообразий параметров цепочки, в окрестности которых могут существовать разумные с точки зрения приложений медленные траектории.

В конечном счете мы имеем дело с семейством решений, зависящих от некоторых параметров. Запишем эти параметры как Г = (71,72 ■ • .)■ Предположим, что мы построили семейство таких решений и траектория движения описывается некоторой функцией X(t, Г) = (Xt (t, Г), X2(t, Г)). Тогда, исходя из предположения о том, что функции (1) описывают мезомасштабные вихри, и зная траекторию центра вихря (глаза тайфуна) при t € [О, Т], мы можем попытаться предсказать последующее движение центра следующим образом. Пусть Г^ = (xi = Х^ (t), х$ = -^гЧ*))» t £ £2])) — траектория центра настоящего природного вихря. Выберем параметры Г = (71,72.) из того условия, что траектории Г и Г7' должны быть близки на интервале [О, Г0], например, в смысле минимизации среднего значения квадратного корня из разности между теоретической и наблюдаемой траекториями вихря. Тогда, если нам известны параметры 71,72,., мы можем однозначно определить траекторию Г при t > t2- Таким образом, приходим к классической задаче оптимизации.

В связи с соображениями, высказанными в этом пункте, еще раз отметим роль единственности структуры типа квадратного корня. Если бы утверждение относительно единственности отсутствовало, другой выбор функции F в (1) (например, F = т1/3 или г log г и т.д.) мог бы привести к совершенно другим формулам для возможных траекторий, скоростей и пр. Таким образом, можно сказать, что формулы для траекторий мезомасштабных вихрей, скоростей, геопотенциалов и т.д., которые получены в этой работе, являются необходимым условием (хотя, конечно, они получаются с помощью довольно грубого приближения).

Краткое содержание и структура работы

Изложение построено следующим образом. В главе 1 приводится постановка задачи и формулируются результаты, касающиеся структуры особого решения и вида цепочек Гюгонио-Маслова. Ищутся решения системы уравнений (2) вида: v(x,t) = p(x,t) + p(x,t), p{x,t) = R(x,t)F(S{x,t)), u (®,t) = u{x,t) + u{x,t), u{x,t) = U{x,t)F(S{x,t)),

4)

Здесь t e [0,T]; p, R и S гладкие скалярные функции, и — *{v,w) и U = l(Ui, С/2) гладкие двумерные вектор-функции. Предполагается, что функции удовлетворяют одним или нескольким из следующих условий i) для любого t S(x,t) > 0, и равенство S(x,t) = 0 выполняется в любой момент времени t G [0,Т] в единственной точке х = X(t) = {t),X2(t)). Множество Г = {х = X(t),t Е [0,Т]} назовем траекторией особого решения (4) уравнений (2) на интервале времени [О, Г].

Функции S, U и R находятся в общем положении, в том смысле, что

И) Матрица Н(«) = Hess S вторых производных невырождена и, как следует из условия (i), положительна на Г. iii) Собственные значения матрицы Hess S|r различны (показано, что если существует решение (4), такое, что S(x,t) удовлетворяет сформулированному предположению при t = to, то это предположение выполняется также для любого t >tQ). iv) разложение U и R по степеням (х — X(t)) начинается с наименее возможных степеней. (Это предположение иногда можно опустить, но оно является естественным с точки зрения возможных приложений и, кроме того, позволяет упростить большинство рассуждений.) v) Функция F удовлетворяет следующим условиям: (v-a) Функция F(r) непрерывна при т > О, .F(O) = 0; (v-b) F(t) является гладкой при т > 0, limr>+0 F'(t) = 00.

Рассматриваются только те решения системы (2), для которых 77 является строго положительной (это следует из ее физического смысла). Пусть gij(t) = i\j\dt+^g(x, t)/dx\dx32\r для любой гладкой функции д(х, t). Положим также щ = о;|г ир0 = 77|г- Справедлива

Теорема 1. Пусть система (2) имеет решение (4), удовлетворяющее условиям (i)-(v). Тогда имеют место следующие утверждения.

1. Без ограничения общности можно считать, что в (4) F = у/т.

2а. Траектория X(t) вморожена в поле скоростей и (а также в и):

X(t) = u(X{t),t) = u(X{t),t) Wv(t) = \Vx{t),V2{t)) (5)

2b. С точностью до функций, обращающихся в ноль вместе со всеми своими производными на траектории Г, функция S является решением уравнения эйконала (или Гамильтона-Якоби)

St + {u,VS) = 0. (6)

2с. Комплексные скорости u(x,t) = v(x,t) + iw(x,t) (и Ui (х, t) + iu.2(x, t)) на траектории X (t) удовлетворяют условиям Коши—Римана dv dw def dv . dw *./ и

Р2о\т = Po21 г + £Vi/2, Pair = PV2/2. (8)

2d. Вторая поправка к условиям Коши-Римана имеет вид

2Н - tr(H)I) ( V20 + + Щ01 ~Я) = tr(H) ( ~ ~Wl1 V v V UJ20 + W02 ~ TqCPW J ' \W02-W20-Vn J

9)

2e. "Потенциальный вихрь" П сохраняется вдоль траектории Г; и0 - 2р

П г =-— с = const.

Ро

2f. Функции р uu — Uy/S определяются формулами где гладкие функции S, (р и вектор-функция Р имеют вид s = ed&(x - X(t), n(t)BW(t)(x - X(t))) + 0(1® - X(t)I3), P = VS = p0(t)n(t)Bn*(t)(x - X(t)) + 0(jx - X(t)I2, j-j//)\ ( cos0 sine \

Щ") ~ I cos0 J матРиЧа поворота на угол

6(t) = 60 + Jq p(t)dt = 0O~I f*rot3u(X(t), t)dt, <p = A + 0{\x-X{t)\), A = const, rt ( 1 + b О Л , . jd = ( ^ 1 — 6 /' ' vq и A - вещественные константы, характеризующие начальную структуру вихревого решения, |Ь| < 1.

Интересен тот факт, что возникают условия Коши-Римана (7), а также "поправки" (9) к ним. Эти условия, как правило, не выполняются для произвольных траекторий поля скоростей, и, следовательно, они описывают влияние существования вихря на "гладкий фон" и(х, t).

Цепочка Гюгонио-Маслова для особых решений вида (4) -бесконечная цепочка уравнений, связывающая коэффициенты в разложении функций фона u(x,t) и р(х, t) в ряд Тейлора по i в окрестности траектории X(t) с учетом соотношений (7) и (8). Для описания траекторий особенностей решений системы (2) цепочка замыкается в духе работы Ф.Прасада и Р. Равиндран1: слагаемые, содержащие коэффициенты полиномов третьего порядка в разложении р и и приравниваются к нулю. Это приводит к системе 17 обыкновенных дифференциальных уравнений, связывающих коэффициенты тейлора в разложении р и и до второго порядка включительно и включающих уравнения X = V, <jjQ = fiV2. Цепочка Гюгонио-Маслова приводится ниже (см. формулы (13)—(21)) в новых комплексных переменных.

В главе 2 доказывается пункт 1. Теоремы 1. Схема доказательства выглядит следующим образом. Прежде всего, удобно перейти к подвижной системе координат (x',t), где х' = х — X(t), и положим и' = и — V. Тогда система (2) записывается в виде (мы опускаем штрихи на новых переменных х)

7?t + (V, т)и) = 0, ut + (и, V)u + V77 + V - шТ(и + V") = 0. (10)

После подстановки решения вида (4) в уравнения (10), получаем уравнение вида

F'{t)A + F(T)B + F(T)F'(T)C + F2(T)D + E = 0. с гладкими векторными коэффициентами А, В, С, D и Е. Здесь г = S(x,t). Анализ этого уравнения приводит к следующему утверждению Лемма (о модельных уравнениях). Пусть выполнены условия (ii)-(iii), тогда справедливы следующие утверждения А. Существуют гладкие функции а, Р, 7, S, зависящие от (х, t) и такие, что F удовлетворяет уравнению Риккати а{х, t)F'{r) + Р(х, t)F2(r) + 7 (ж, t)F(r) + 6(х, t) = 0, где производные по х функции а не равны нулю при х — 0.

1 Ravindran R., Prasad P. A new theory of shock dynamics. Part I (П)// Appl. Math.

Lett. - 1990. - 3, № 3. - С 77-79

В. Функция F удовлетворяет одному из следующих трех "модельных" уравнений:

- квадратному уравнению: а(х, t)F2(r) + b(x, t)F(r) + с(х, t) = О, линейному дифференциальному уравнению: a(x,t)F'(r) + b(x,t)F(r) + c(x,t) = О,

-1или кубическому уравнению: a(x,t)F3(T) + b(x,t)F2(r)+c(x,t)F(T) + d(x,t) = 0. Здесь а, Ь, с, d также гладкие функции, некоторые производные которых по х-не равны нулю при х = О.

Анализ модельных уравнений приводит к следующему утверждению.

Лемма. Без ограничения общности можно считать, что существует всего четыре возможные функции F в (4): (Fl) F = тк, (F2) F = г log г, (F3) F = т1/3 + от2/3+п, (F4) F = г2/3 + <тт4/3+п, где 0 < к < 1, а = ±1 или 0, п > 0 целое число. Квадратное уравнение имеет лишь решение типа (F1) с к = 1/2; линейное дифференциальное уравнение имеет решения типа (F1) и (F2); кубическое уравнение имеет решения (F1) с «=1/2, (F3) и (F4).

Далее решения вида (4) с функциями F из (F1)-(F4) подставляются в уравнения (10). При к ф 1/2 получившиеся системы уравнений оказываются несовместны. Таким образом справедливо следующее утверждение.

Лемма. Решений с особенностью вида (F1)-(F4) при к, ф 1/2 не существует.

Из перечисленных лемм следует утверждение пункта 1. Теоремы 1.

В главе 3 изучается негладкая часть решения и доказываются пункты 2a.-2f. Теоремы 1. После подстановки решения вида (4) с F — л/г в (10) уравнения для (и, р) и (и, р) распадаются следующим образом:

Щ +V • (ри + рг!) = 0, g + = 0, ( . |f +{u,V)u+(u,W)u + Vp-uTu = 0. К '

Показывается, что хотя система (12) линейная относительно и и р, нельзя напрямую применить методы разложений по гладкости или лучевых разложений2 для задачи Коши p|t=0 = R(x,0)y/S(x, 0),tt|t=o = U(x, 0) y/S(x, 0). По существу, используются идеи, восходящие к этим методам, однако имеются некоторые существенные отличия,

2Бабич В. М., Булдырев В. С., Молотков И. А. Пространственно-временной лучевой метод. - Ленинград: ЛГУ, 1985 существенные с точки зрения лучевого метода. Прежде всего, использование разложения по гладкости не гарантирует структуры (4) решения, даже если эта структура имеет место при t = 0. Кроме того возникает ситуация, известной в теории гиперболических уравнений как "негладкая смена кратности характеристик": собственные частоты Ло = St + (u,V)5, А± = |St + («,V)S| ± p\VS\, соответствующие гидродинамической и акустической моде одновременно обращаются в ноль на траектории центра особенности X(t). Таким образом, рассматривается вся система (11)—(12) и изучитаются ее свойства. С первого взгляда может показаться, что система (11)-(12) недоопределена: имеются шесть уравнений относительно семи неизвестных функций v,w,pnUi, U2, R, S. Показывается, что заданные свойства функции S и гладкость всех названных функций фактически приводит к возникновению уравнения эйконала для фазы S и система замыкается.

Процедура аналитических вычислений состоит в разложении всех функций в ряды Тейлора, подстановки этих рядов в систему (12) и приравнивании коэффициентов при соответствующих степенях х. Указанная процедура приводит к некоторым равенствам и уравнениям для коэффициентов в рядах Тейлора и т.д. При этом прослеживается прямая аналогия с методом ВКБ: для фазы выводится уравнение эйконала (6), а условия Коши-Римана (7) получаются из уравнения, аналогичному уравнению переноса. Показывается также, что условия Коши-Римана можно получить из условия сохранения потенциального вихря. Также, выводятся поправки к условиям Коши-Римана следующего порядка (9), которые интересны тем, что в уравнения (9), дающие условия на гладкий фон входят параметры особой части решения (4).

В главе 4 выводится цепочка Гюгонио-Маслова. Показано, что цепочку удобнее всего получить, если вместо вещественных переменных xi, Х2 в исходных уравнениях мелкой воды, а также в описании функций р, и использовать комплексные переменные z = xx- Xx{t) + г(х2 - Х2(f)), z = x1-Xl(t) - i(x2 - X2(t)).

После введения новых комплексных переменных X' = Хх + iX2, V' = V\ +iV2 ,Y,Z,U,W и вещественных переменных ц и А (эти переменные выражаются через коэффициенты в разложении функций р и и в ряд Тейлора по г, г в окрестности траектории X(t) до второго порядка включительно) цепочка Гюгонио-Маслова принимает вид

X — V, (13)

V + iu0V - + W — 2Z) = 0, (14)

У = i(p - u0)Y - my^ + fu (15)

Z = i(3p-u0)Z+^ + f2, (16)

W = -ipW,; (17)

U = i(p + u0)U + f3, (18)

Л = 2^3(f ImY - Z^ImZ - flmW), (19) u>q = filmV = J3V2, (20)

Mx+- +lZY))- (21)

Здесь слагаемые /i, /2 и /3 содержат коэффициенты Тейлора 3-го порядка. Обрыв цепочки (замыкание) означает, что /1, /2 и /3 приравниваются к нулю.

В главе 5 изучается оборванная цепочка (13)—(21) в трех случаях: = 0, уб - мало, и случай, когда р нельзя считать параметром возмущения. В работах С.Ю. Доброхотова3 было показано, что оборванная цепочка (13)-(21) точно интегрируется и сводится к уравнению Хилла в случае /9 = 0

Это сведение позволяет выявить ряд свойств траекторий Г связав их со свойствами уравнения Хилла (такими, как свойства устойчивости, рассмотрением зон и лакун в спектре и т.д.). Физически разумное предположение о достаточной гладкости траекторий Г и вполне элементарное исследование уравнений для скорости и функций У, Z, W приводят к требованию малости малости одной из следующих частот: vi = (р — u}0) ) мала, тогда Z ~ 0, W ~ 0 (случай I); i/2 — (ujo — 3р) ) мала, тогда Y ~ 0, W ~ 0 (случай И); vz = -p) мала, тогда У ~ 0, Z ~ 0 (случай Ш). Это требование должно быть выполнено как в случае (3 — 0, так и в случае малых /5.

В случае /? = 0 у системы (13)—(21) легко получаются траектории в виде окружности, если положить в случаях I, II, Ш последовательно (Z = 0, W = 0), (У = 0, W = 0), {Z = 0, У = 0). Если

3 Доброхотов С. Ю. Редукция к уравнению Хилла цепочки Гюгонио-Маслова для траекторий уединен ных вихрей уравнений мелкой воды// Теор. и мат. физ. - 1997. -112, № 1. - С. 47-66 считать параметр /3 малым, то слагаемые с (3 в (13)—(21) играют роль адиабатического возмущения, которое, в частности, "деформирует" указанные траектории. Одним из эффектов влияния этих слагаемых состоит в том, что движение по "деформированным" окружностям может менять направление и появляются траектории Г в виде зигзагов. Такие "деформированные" траектории могут быть описаны с помощью гамилътоновых систем с одной степенью свободы. Для всех случаев I, П и Ш они имеют одинаковую структуру дН • дН Uo = f-щ, Ф = woU = wb(0), = т,

Н = cos фМ + N, но с разными функциями /, М и N. Здесь гамильтоновыми переменными являются частота и>0 и угол ф направления скорости V. Например, для случая, когда мала частота Зр — и0,

N = П7 у9(5у7 - 6)

4 Л yj[у7 + 3)(2у7 — I)3

4^/0^0(277 — I)3 Шо 1 ~ 7t7v^T3 ,T?~W2' где Mi и Шг - константы интегрирования. Ее фазовые портреты на плоскости (w0, ф) определяются интегралом системы Н: Н = const, и для отдельных областей почти совпадают с хорошо известными фазовыми портретами траекторий физического маятника (см. Рис. 1).

В каждом из 3-х случаев, траектория X(t) зависит от шести параметров Mi: с, ф0, X®, (так как ш0 = О, 4- (ЗХ®)- Эти шесть "существенных" параметров могут быть восстановлены при помощи полученных в диссертационной работе формул через шесть имеющих физический смысл величин: начальные положение (Х^Лз) и скорость (V?, У2°) центра вихря, потенциальный вихрь с и третья компонента rot3 v в центре вихря.

При выводе указанных в конце предыдущего пункта гамильтоновых систем существенным является предположение об относительной малости слагаемых, пропорциональных /3 ("вынуждающих сил") в правых частях уравнений для Y и Z. Если эти слагаемые не малы, то ясно, что их нельзя рассматривать как адиабатическое а)

2.5 2 .25 2

1 .75 1. 5

1.25 1

О .75

Ь) х2,кт

700

-125 -100 -75

Рис. 1: Пример фазового портрета гамильтоновой системы (а) и траектории центра вихря на /3-плоскости в случае когда мала частота 3р — шо (Ь). Траектория, отмеченная жирной линией, соответствует траектории гамильтоновой системы с одной степенью свободы, траектория, отмеченная тонкой линией, - решению 17 уравнений оборванной цепочки Гюгонио-Маслова до усреднения. возмущение, и более того, они могут играть существенную роль в описании траекторий Г. Имея ввиду явные выражения для правых частей в системе (13)—(21), предполагается, что для интересующих нас решений в первом приближении всеми производными по t можно пренебречь, что приводит к уравнению для определения критических точек системы. Элементарный анализ показывает, что эти критические точки образуют два двумерных (пересекающихся) гладких "критических" многообразия, характеризуемых условиями (а) W = 0 и (b) р = 0. Таким многообразиям отвечают траектории с равномерным движением, параллельным оси х\ (т.е. движениям вдоль параллели), причем случай (Ь) не дает подходящих для описания мезомасштабных вихрей траекторий. В первом приближении движение определяется линеаризованными уравнениями в некоторой окрестности этих критических многообразий. Показано, что собственные значения соответствующей матрицей Якоби находятся из уравнения

А? + 22^ + (105 + jН4Аг4 + (1484- ™)«0вА? + (64 - 815 - ^>08 = 0, яЧ где безразмерный параметр S = .

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

Одно из собственных значений обращается в ноль тогда и только тогда, когда слагаемое (64 — 815 — в последнем уравнении обращается в ноль. Это дает 6 = ~ 0,465. Например, в случае, когда

5 плоскость касается земли в точке, расположенной на широте 20° (и ~ 0.18 1/час) такая ситуация реализуется, если принять значение геопотенциала на траектории Г соответствующее высоте слоя атмосферы ~ 10км (т]\г ~ 860000км2/час2). Иными словами, в окрестности широты 20° возможно наличие медленных траекторий, если предположить, что воображаемая высота слоя атмосферы равна 10 км. Заметим, что в задачах о мезомасштабных вихрях в атмосфере этой величине часто приписываются значения 7-5-12 км.

В главе 6 описывается цепочка Гюгонио-Маслова для более сложной системы уравнений, описывающих динамику двумерных вихрей в атмосфере с учетом энергетического обмена между атмосферой и океаном4: После обезразмеривания эта система выглядит следующим образом. о u • Vu - wTu + Vfj + (rj + 6)VQ = 0. (22)

С/V + V • (T)VL) + £rot3u - 77* (0 - в) = 0. (23) dQ u- V©+ *(©-©) = 0. (24)

В отличие от уравнений мелкой воды, добавляется еще одно уравнение для скалярной функции Q(x,t) (плотность или температура), rot3U = дщ/дх1 — дщ/дх2, S, к - физические константы, характеризующие приток массы из пограничного слоя атмосферы (слоя Экмана) и нагрев вихря, © - соответствует значению температуры поверхности океана.

Показывается, что для особых решений типа корня квадратного из квадратичной формы выполнены условия вмороженности и уравнение Гамильтона-Якоби вида (6) для фазы S. Условия Коши-Римана (7) тоже сохраняются, если потребовать дополнительные условия на гладкую составляющую 6(x,t) функции Q(x,t), а именно V0(x,£)|г = 0. Численный анализ соответствующей цепочки Гюгонио-Маслова показал, что у оборванной цепочки существует семейство траекторий близких по параметрам к описанному в предыдущей главе режиму, но при этом центр особенности может двигаться в меридиональном направлении.

4Должанский Ф. В., Крымов В. А., Манин Д. Ю. Устойчивость и вихревые структуры квазидвумерных сдвиговых течений// Усп. физ. наук. - 1990. - 160, №

7. - С. 1-47

Доклады и публикации

Основные результаты работы были представлены на семинаре профессора Альбеверио в Институте прикладной математики Боннского Университета, на Международном семинаре "Дни Диффракции" в 2002, 2003, 2004 годах, на международном семинаре "Asymptotic Analysis and the Physics of Atmosphere and Ocean" в 2004 году, на международном семинаре "Mathematical aspects of tsunami" в 2005 году, а также отражены в 6 публикациях [1, 2, 3, 4, 5, 6].

Благодарности

Результаты этой работы были получены в рамках проектов РФФИ 9900415 и 01-02-00850. Результаты главы 5, б были получены совместно с профессором Б.Тироцци из Университета "La Sapienza" (Рим). Результаты главы 2 были получены совместно с К.В.Панкрашкиным. Автор диссертации выражает им свою благодарность.

В ходе работы над диссертацией автор имел чрезвычайно полезные дискуссии с В.Г.Даниловым, В. Ф. Должанским, В. П. Масловым, А.И.Шафаревичем, А.Сперанса, которым автор также выражает свою признательность.

Особую благодарность автор хочет выразить своему научному руководителю и наставнику Сергею Юрьевичу Доброхотову за неоценимую помощь и поддержку в научной работе.

1. Доброхотов С. Ю., Панкрашкин К. В., Семенов Е. С. О гипотезе Маслова о структуре слабых точен ных особенностей уравнений мелкой воды// Докл. РАН. - 2001. - 379, № 2. - С. 173-176

2. Dobrokhotov S. Yu., Pankrashkin К. V., Semenov Е. S. Proof of Maslov's conjecture about the structure of weak point singular solutions of the shallow water equations// Russ. J. Math. Phys. 2001. - 8, № 1. - С 25-54

3. Семенов E. С. Об условиях Гюгонио-Маслова для вихревых особых решений системы уравнений мелкой воды// Мат. заметки. 2002. - 71, № 6. - С. 902-913

4. Доброхотов С. Ю, Семенов Е. С., Тироцци В., Цепочки Гюгонио-Маслова для сингулярных вихревых решений квазилинейных гиперболических систем и траектории тайфунов// Современная математика. Фундаментальные направления. 2003. - 2. С. 5-44

5. Доброхотов С. Ю., Семенов Е. С., Тироцци В., Вычисление интегралов цепочки Гюгонио-Маслова для особых вихревых решений уравнений мелкой воды// Теоретичесая и математическая физика. 2004. - 139, №1.- С. 500-512

6. Арнольд В. И., Козлов В. В., Нейштадт А. И. Математические аспекты классической и небесной механики// Итоги науки и техн., сер. Совр. пробл. мат., Фундам. напр. М.: ВИНИТИ, 1985. - 3

7. Бабич В. М. Фундаментальное решение гиперболических уравнений с переменными коэффициентами// Мат. сборник. 1960. - 52 (94), JV® 2. -С. 709-738

8. Бабич В. М., Булдырев В. С., Молотков И. А. Пространственно-временной лучевой метод. Ленинград: ЛГУ, 1985

9. Булатов В. В., Владимиров Ю. В., Данилов В. Г., Доброхотов С. Ю. Пример вычисления "глаза" тайфуна на основе гипотезы В. П. Маслова// Докл. РАН. 1994. - 338, Л» 1. - С. 102-105

10. Вишик М. И., Фурсиков А. В. Математические проблемы статистической механики. М: Наука, 1980

11. Гордин В. А. Математические задачи гидродинамического прогноза погоды: аналитические аспекты. JL: Гидрометеоиздат, 1987

12. Гринфельд М. А. Лучевой метод вычисления интенсивности волновых фронтов в нелинейном упругом материале// Прикл. мат. и мех. 1978. -42, JV* 5. - С. 883-898

13. Данилов В. Г., Маслов В. П., Шелкович В. М. Алгебры особенностей обобщенных решений строго гиперболических систем квазилинейных уравнений первого порядка// Теор. и мат. физ. 1998. - 114, JV® 1. - С. 3-55

14. Доброхотов С. Ю. Цепочки Гюгонио-Маслова для траекторий точечных вихревых особенностей урав нений мелкой воды и уравнение Хилла// Докл. РАН. 1997. - 354, № 5. - С. 600-603

15. Доброхотов С. Ю. Редукция к уравнению Хилла цепочки Гюгонио-Маслова для траекторий уединен ных вихрей уравнений мелкой воды// Теор. и мат. физ. 1997. - 112, X* 1. - С. 47-66

16. Доброхотов С. Ю., Тироцци Б. О свойстве гамильтоновости укороченных цепочек Гюгонио-Маслова для траекторий мезомасштабных вихрей// Докл. РАН. 2002. - 384, № 6. - С. 741-746

17. Должанский Ф. В., Крымов В. А., Манин Д. Ю. Устойчивость и вихревые структуры квазидвумерных сдвиговых течений// Усп. физ. наук. 1990. - 160, № 7. - С. 1-47

18. Дубровин Б. А., Новиков С. П. Гамильтонов формализм одномерных систем гидродинамического типа и метод усреднения Боголюбова-Уизема// Докл. АН СССР. 1983. - 270, № 4. - С. 781-785

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

20. Егоров Ю. В. К теории обобщенных функций// Усп. мат. наук. 1990. -45, № 5. - С. 3-40

21. Зайцев В. Ф., Полянин А. Д. Справочник по обыкновенным дифференциальным уравнениям. М.: Физматлит, 1995

22. Захаров В. Е. , Кузнецов Е. А.// Успехи физ.наук. 1997. -167, № 11. - С. 405-432

23. Danilov, V.G., Omel'yanov, and Rozenknop, D. Calculation of the hurricane eye motion based on singularity propogation theory// Electr. Jour, of Differential Equations, 2002. 16. - С 1-17

24. Dobrokhotov S. Yu. Hugoniot-Maslov chains for solitary vortices of the shallow water equations, I, II// Russ. J. Math. Phys. 1999. - 6, JV® 2. - С 137-173, № 3. - С 282-313

25. Hongyan Zhu, Roger K. Smith, and Wolfgang Ulrich. A Minimal Three-Dimensional Tropical Cyclone Model// Journal Of The Atmospheric Sciences.- 2001. 58. - C. 1924-1944

26. Lions J.-L., Temam, R., and Wang, S. New formulation of the primitive equations of atmosphere and applications// Nonliniarity. 1992. -5, - С 237-288

27. Milnor J. Morse theory// Ann. Math. Stud. 1993. - № 51

28. Ooyama К. V. Conceptual evolution of the theory and mod-eling of the tropical cyclone// J. Meteor. Soc. Japan. 1982. - 60. - C. 369-379

29. Ravindran R., Prasad P. A new theory of shock dynamics. Part I (II)// Appl. Math. Lett. 1990. - 3, № 3. - С 77-79

30. Reznik G. M., Grimshaw R. Ageostrophic dynamics of an intence localized vortex on a beta]-plane// J. Fluid Mech. 2001. - 443. - С 351-376

31. Rogers C., Schief W. K. Multi-component Ermakov systems: structure and linearization// J. Math. Anal. Appl. 1996. - 198, № 1.- 0 194-220

32. Rogers C., Elliptic Warm-Core Theory: The Pulsrodon// Physics Letters A.- 1989. 6, 7. - C 138

33. Shapiro L. J. Potential vorticity asymmetries and tropical cyclone evolution in a moist three-layer model// J. Atm. Sc. 1999. - 57, JV» 21. - C 3645-3662

34. Shugaev F. VM Shtemenko L. S. Propagation and Reflection of Shock Waves.- Singapore: World Scientific, 1998

35. Smith R. K., Ed. The Physics and Parameterization of Moist Atmospheric Convection. Kluwer, 1997

36. Trev F. Propagation of singularities and semi-global existence theorems for (pseudo-) differential operators of principal type// Ann. Math. Stud. 1978.- № 108. C. 569-609

37. Wang Y. An Explicit Simulation of Tropical Cyclones with a Triply Nested Movable Mesh Primitive Equation Model: TCM3. Part I: Model Description and Control Experiment// Monthly Weather Review. 2000. - 129. - C. 13701394