WWW.DISS.SELUK.RU

БЕСПЛАТНАЯ ЭЛЕКТРОННАЯ БИБЛИОТЕКА
(Авторефераты, диссертации, методички, учебные программы, монографии)

 

Pages:   || 2 | 3 | 4 | 5 |   ...   | 8 |

«УНИВЕРСИТЕТ НАУК И ТЕХНОЛОГИЙ (ЛИЛЛЬ, ФРАНЦИЯ) RUSSIAN ACADEMY OF SCIENCES A. M. OBUKHOV INSTITUTE OF ATMOSPHERIC PHYSICS UNIVERSITE DES SCIENCES ET TECHNOLOGIES DE LILLE (FRANCE) V. P. ...»

-- [ Страница 1 ] --

РОССИЙСКАЯ АКАДЕМИЯ НАУК

ИНСТИТУТ ФИЗИКИ АТМОСФЕРЫ им. А. М. ОБУХОВА

УНИВЕРСИТЕТ НАУК И ТЕХНОЛОГИЙ

(ЛИЛЛЬ, ФРАНЦИЯ)

RUSSIAN ACADEMY OF SCIENCES

A. M. OBUKHOV INSTITUTE OF ATMOSPHERIC PHYSICS

UNIVERSITE DES SCIENCES ET TECHNOLOGIES

DE LILLE (FRANCE) V. P. Goncharov, V. I. Pavlov

HAMILTONIAN

VORTEX AND WAVE

DYNAMICS

Moscow GEOS 2008 В. П. Гончаров, В. И. Павлов

ГАМИЛЬТОНОВАЯ

ВИХРЕВАЯ И ВОЛНОВАЯ

ДИНАМИКА

Москва ГЕОС УДК 532.50 : 551.46 + 551. ББК 26. Г Гончаров В. П., Павлов В. И.

Гамильтоновая вихревая и волновая динамика. — М.: ГЕОС, 2008. — 432 с.

ISBN 978-5-89118-412- Монография отражает развитие и совершенствование моделей и методов, направленных на изучение проблем вихревой и волновой динамики. Книга содержит замкнутое и относительно элементарное изложение того, что известно в современной научной литературе как гамильтонов подход, применительно к гидродинамическим моделям вообще, и к геофизической гидродинамике, в частности. В сравнении с книгой «Проблемы гидродинамики в гамильтоновом описании» [48] тех же авторов, изданной в 1993 г., данная монография претерпела значительные изменения. Вопервых, гамильтонов подход рассматривается в более широком контексте, и прежде всего, как инструмент создания новых более эффективных моделей и асимптотических методов. Во-вторых, в книге нашли отражение последние достижения в области приложений. В частности, включены новые разделы, посвященные гамильтонову подходу к гидродинамическим моделям плазмы, жидких кристаллов, контурной динамики и анизотропным волновым системам. Особое внимание уделено изучению крупномасштабных волновых и вихревых структур, важных для понимания процессов вертикального и горизонтального перемешивания в атмосфере и океане. Теоретические результаты, изложенные в монографии, представляют общефизический интерес и могут найти применение в теории гидродинамической устойчивости, турбулентности и динамической метеорологии.

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

Издание осуществлено при финансовой поддержке Российского фонда фундаментальных исследований по проекту № 08-05-07048.

Goncharov V. P., Pavlov V. I.

Hamiltonian vortex and wave dynamics. — M.: GEOS, 2008. — 432 p.

The monograph is devoted to a development and an improvement of models and methods for studying problems of vortex and wave dynamics. The book contains an all-sufficient yet relatively elementary insight into that which is named in modern scientific literature the Hamiltonian approach as applied to hydrodynamic models in general and to geophysical fluid dynamics in particular. This monograph is a significant upgrade in comparison with the book «Problems of fluid dynamics in Hamiltonian descriptions» [48] published by the same authors in 1993. First, the Hamiltonian approach is considered in a broader context and as a tool to create new more effective models and asymptotic methods. Second, the book includes the latest advances in applications. In particular, it contains new sections devoted to the Hamiltonian approach for hydrodynamic models of plasma, fluid crystals, contour dynamics and anisotropic wave systems. Special attention is given to the large-scale wave and vortex structures, which are important for understanding processes of vertical and horizontal mixing in the atmosphere and the ocean. Theoretical results presented in the monograph are of general physical interest and can find a wide use in the theory of hydrodynamic stability, turbulence and dynamic meteorology.

The book is suitable for researchers specializing in fluid dynamics and its applications, theoretical and mathematical physics, applied mathematics as well as for post-graduate and advanced students.

c В. П. Гончаров, В. И. Павлов, c ГЕОС, Оглавление Предисловие............................................................ Часть I

ОБЩЕТЕОРЕТИЧЕСКИЕ АСПЕКТЫ

ГАМИЛЬТОНОВСКОГО ФОРМАЛИЗМА

Глава 1.

Глава 6.

§ 6.1. Преобразования симметрии в канонических системах. Теорема

ГАМИЛЬТОНОВ ПОДХОД

К СИСТЕМАМ ГИДРОДИНАМИЧЕСКОГО ТИПА

Глава 8.

Глава 9. Гидродинамические законы сохранения и интерпретация § 9.1. Субстанциональные характеристические функции и интегралы § 9.2. Топологическая и физическая интерпретации потенциалов § 10.1. Канонические импульсы как нефизические переменные и их Оглавление § 11.3. Калибровочные преобразования в гамильтоновых системах с § 11.4. Канонические преобразования, сохраняющие спиральность... Глава 12. Гидродинамические системы с некоммутативными полями Глава 13.

§ 14.2. Принцип наименьшего действия для баротропной спиновой Глава 15.

КАНОНИЧЕСКИЕ ПЕРЕМЕННЫЕ





ДЛЯ ГИДРОДИНАМИЧЕСКИХ МОДЕЛЕЙ

Глава 17.

Глава 18.

Глава 19. Канонические переменные для 2D-моделей несжимаемой § 20.4. Каноническое описание многослойных квазигеострофических Глава 23.

Глава 24.

Оглавление

ГАМИЛЬТОНОВ ПОДХОД

К МОДЕЛЯМ КОНТУРНОЙ ДИНАМИКИ

Глава 26.

Глава 27.

Глава 28.

§ 28.4. Контурная динамика для бароклинной квазигеострофической § 30.3. Гамильтонова динамика границ раздела в отсутствие вихревой

ФОРМАЛИЗМ НОРМАЛЬНЫХ МОД:

ПРИЛОЖЕНИЯ И ОБОБЩЕНИЯ

§ 31.1. Стационарные точки и метод линеаризации для гамильтоновых § 31.3. Гамильтоново описание пространственно однородной волновой § 32.1. Нормальные переменные в равновесной однородной волновой Глава 33.

Глава 35. Введение нормальных переменных в неоднородные модели Оглавление Глава 39. Метод нормальных мод в задачах излучения и рассеяния § 39.1. Типичные параметризации для гамильтониана взаимодействия

ГАМИЛЬТОНОВ ПОДХОД К ЗАДАЧАМ

ВОЛНОВОЙ И ВИХРЕВОЙ ДИНАМИКИ

Глава 41.

Глава 42.

Глава 43.

Глава 45.

ПРЕДИСЛОВИЕ

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

Основная установка книги такова: большинство интересных физических явлений так сложно, что при современном состоянии наших знаний вряд ли можно рассчитывать на построение универсальной теории, действующей во всей области параметров и на всех временах рассматриваемого процесса. Следовательно, нужно попытаться понять ведущие факторы, которые в тот или иной отрезок времени управляют процессом в той или иной области параметров. Но даже выделив существенные факторы, не менее важно сделать второй шаг — построить простую, но адекватную математическую модель процесса. При этом выбор динамических переменных, который гарантировал бы сравнительно простые уравнения движения в фазовом пространстве этих переменных и наличие эффективных аналитических методов их решения, — определяющий момент любой формулировки. На примере гидродинамики можно убедиться, что в этом отношении возможности традиционных формулировок в значительной мере исчерпаны.

На современном этапе именно эта причина стимулировала развитие и применение преимущественно численных методов и отодвинула в тень аналитические методы исследования.

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

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

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

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

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

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

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

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

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

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

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

Суть указанного подхода заключается в следующем. Любая гидродинамическая система характеризуется набором конечного числа полевых переменных, таких, как плотность жидкости, скорость, давление, температура и т.п. В традиционной формулировке эволюция идеальной системы описывается уравнениями Эйлера. Каждая полевая величина — функция времени, t, и координат, x = = (x1, x2, x3 ), т. е. имеет континуальное число степеней свободы по параметру x. По этой причине для таких систем, если мы переходим к гамильтоновой формулировке, должны быть использованы функциональные производные.

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

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

Оператор /uj означает функциональную производную по полю uj и является некоторым обобщением обычной частной производной.

Все строгие определения мы дадим ниже в основном тексте книги. Здесь же, в Предисловии, важно подчеркнуть, что возможность такой записи уравнений движения является исключительно правом так называемых консервативных моделей и означает, что динамическая система обладает гамильтоновой структурой, которая определяется гамильтонианом H, совпадающим с полной энергией жидкости, и функциональной скобкой Пуассона {ui (x), uj (x )}, выполняющей роль ядра интегрального преобразования. Вид скобки Пуассона определяется конкретным типом гидродинамической модели.

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

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

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

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

Указанные соображения дают основную идею последующей стратегии и объясняют в какой-то мере основную причину привлекательности гамильтоновского подхода: гидродинамические системы классифицируются по определенным критериям, все расчеты стандартизируются, гамильтонианы сводятся только к нескольким существенным, а следующие из них эволюционные уравнения — к нескольким базовым.

Современный этап развития гамильтоновского подхода к бесконечномерным полевым системам стартовал приблизительно лет тридцать назад. До этого, полевые версии гамильтоновского метода рассматривались только в связи с потребностями квантовой теории поля [10, 207]. Применительно к гидродинамическим моделям, первоначально, гамильтонов подход был предложен в основном как эффективный метод для нахождения нелинейных эволюционных уравнений в системе взаимодействующих волн. Было немедленно осознано, что гамильтонов подход обладает рядом преимуществ по сравнению с традиционными методами.

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

В начале семидесятых, после работ В. Е. Захарова, Д. Д. Фаддеева [93] и К. С. Гарднера [216], стало ясно, что большинство разрешаемых методом обратной задачи теории рассеяния нелинейных эволюционных уравнений являются бесконечномерными аналогами гамильтоновских уравнений классической механики. Приблизительно в то же время было понято, что многие полевые консервативные гидродинамические модели обладают скрытой гамильтоновой структурой. Построение и успешное использование канонических переменных для опиПредисловие сания взаимодействия поверхностных гравитационных волн [81, 94] дали один из наиболее впечатляющих примеров применения гамильтонова подхода в динамике жидкостей и послужили толчком для построения общей теории волновых взаимодействий в диспергирующих средах в рамках канонического гамильтоновского формализма [84, 87].

Определение канонических переменных всегда было одной из основных тем в развитии гамильтоновского метода, и долгое время оставалось искусством. Исторически, для многих гидродинамических систем канонические переменные были найдены не из первых принципов, а на интуитивном уровне (некоторые примеры приведены в работах [43–45, 85, 87, 165, 254], список которых, очевидно, может быть продолжен). Это обстоятельство, по-видимому, повлияло на то, что до сих пор встречаются утверждения (см., например, [128], с. 49), что канонические переменные могут быть только угаданы, хотя для систем гидродинамического типа этот вопрос уже давно не актуален [48]. Регулярная процедура нахождения канонических переменных для таких систем не только существует, но и, в некотором смысле, подобна выбору «системы ортогональных координат» в фазовом пространстве динамических переменных.

Значительные усилия были предприняты для разработки альтернативной версии гамильтоновского описания, основанного на принципе наименьшего действия в динамике жидкостей. Более детальное описание этого направления можно найти в работах [166, 197, 205, 226–228, 240, 244, 245, 272, 273]. В действительности, если рассматривать также публикации, в которых тема гамильтоновского подхода вообще сколь-нибудь глубоко затрагивалась, мы получим значительно более длинный список работ (см., например, работы [189, 231, 233, 248], а также представленные в них ссылки).

Формализм неканонических скобок Пуассона начал интенсивно развиваться в конце семидесятых годов (см. обзоры [139, 211, 241, 254]). На современном уровне классификация и систематическое изучение систем и скобок Пуассона гидродинамического типа проводится в рамках дифференциально-геометрического подхода, который был предложен Б. А. Дубровиным и С. П. Новиковым в [72, 73] и продолжает развиваться их последователями. Весьма качественный обзор этих работ представлен в [137].

Содержание этой книги определяется желанием авторов изложить суть проблематики во вразумительной и относительно полной форме c тем, чтобы книга могла послужить «трамплином» для последующих самостоятельных исследований читателя в области приложений гамильтоновых методов к проблемам волновой и вихревой динамики. Действительно, в настоящее время, несмотря на постоянно возрастающий интерес к гамильтоновскому подходу в гидродинамике, как хорошо известно, учебная или монографическая литература по данному предмету практически отсутствует. Необходимая информация, к сожалению, рассеяна в многочисленных статьях и обзорах. Авторы надеются, что данная книга восполнит существующий пробел.

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

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

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

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

Монография в значительной мере опирается на результаты и методы, полученные и разработанные в ходе исследований, которые проводились при поддержке грантов РФФИ (проекты №№ 06-05-64185, 03-05-64247, 01-05-64466, 00-05-6419, 96-05-64991, 94-05-16950), программы президиума РАН «Математические методы в нелинейной динамике», гранта Президента РФ НШ-4166. 2006. 5 и гранта РФФИ—Франция (проект № 07-05-92211-НЦНИЛ_а).

Пользуясь случаем, авторы выражают глубокую благодарность Г. С. Голицыну, В. М. Грянику, Ф. В. Должанскому, В. И. Кляцкину, В. М. Пономареву, Н. Н. Романовой, И. Г. Якушкину за ценные замечания и плодотворные дискуссии, которые способствовали прояснению ряда затронутых в книге вопросов.

ОБЩЕТЕОРЕТИЧЕСКИЕ АСПЕКТЫ

ГАМИЛЬТОНОВСКОГО

ФОРМАЛИЗМА

Проблема интегрирования нелинейных уравнений в частных производных, которые принято называть эволюционными, является классической и особенно актуальна в гидродинамике. Огромную роль на этом пути играет «правильный» выбор переменных, обеспечивающих сравнительно простые уравнения движения и наличие эффективных методов их решения. Как отмечал еще К. Г. Якоби [188, с. 176], «Главная трудность... состоит во введении удобных переменных, для разыскания которых нет никакого общего правила. Поэтому мы должны идти обратным путем: найдя какую-нибудь замечательную подстановку, разыскивать задачи, в которых она может быть с успехом применена».

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

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

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

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

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

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

МАТЕМАТИЧЕСКИЙ АППАРАТ

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

§ 1.1. Функционалы и вариационные производные В этом разделе мы представим некоторые основные понятия функционально анализа, а затем сформулируем правила функционального дифференцирования.

Как мы увидим (см., например: [106–108, 135, 161]), эти правила во многом аналогичны обычным правилам дифференцирования и являются в некотором смысле их обобщением.

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

где переменные (x, t) и v (x, t) — поля, характеризующие распределение плотности и гидродинамической скорости; x — радиус-вектор; t — время, а интегрирование распространяется на весь объем V, занятый жидкостью.

Прежде всего дадим общее определение функционала. Будем говорить, что задан функционал F [u] или что F [u] функционально зависит от u, если имеется правило, по которому каждой реализации u(x) из некоторой совокупности функций, зависящих от пространственной координаты x и составляющих область определения функционала, ставится в соответствие число F, называемое значением функционала на этой функции. Если величина F не зависит от детального вида функции u(x), а зависит только от ее интегрального поведения, то в обозначении F [u] можно не указывать зависимость от аргумента x. В тех случаях, когда необходимо подчеркнуть, что F не только функционал u, но также и явная функция независимой переменной x или параметра t, будут использоваться обозначения F [u; x] или F [u; t], соответственно.

Простейшим примером является линейный функционал где A(x) — заданная функция.

Для непрерывно-полевых моделей достаточно типичными являются так называемые локально-полевые функционалы. В частном случае одномерного скалярного поля u(x) такой функционал имеет вид где плотность функционала h — функция, зависящая от аргумента x, полевой переменной u и ее производных u(s) s u/xs (s = 1,..., k ). В более общем случае h может быть функцией векторного аргумента x, совокупности полевых переменных и их производных конечного порядка относительно x.

Часто возникает необходимость рассматривать более сложные функциональные объекты, примером которых может служить величина где B(x, x ) — заданная функция двух параметров.

Обозначение [u; x] подразумевает, что величина (1.1.2) не только обладает функциональной зависимостью от u, но также является функцией аргумента x.

Рис. 1.1. К определению вариаци- нулю всюду, кроме некоторой малой окрестности x точки x0 (см. рис. 1.1, который иллюонной производной По определению (см., например, [108, 135, 161]), функциональной (вариационной) производной функционала F [u] в точке x0 называется предел при условии, что этот предел существует и не зависит ни от вида u(x), ни от способа стягивания к нулю объема окрестности x, ни от закона, по которому стремится к нулю максимум модуля функции u(x).

1.2. Дифференцирование функционалов Отметим, что для функциональной производной более общепринято сокращенное обозначение F [u]/u(x0 ), в котором нормировка на элемент объема dx опускается. Во избежание недоразумений об этом полезно помнить, например, при вычислении и сопоставлении размерностей величин, в которые входят функциональные производные.

Из определения (1.1.3), в частности, следует, что вариация функционала и вариация функции связаны линейным интегральным соотношением Ограничиваясь далее функционалами, по отношению к которым операция (1.1.3) может быть заведомо выполнена, можно предложить более удобную формулу для вычисления вариационной производной, если в качестве функции u(x) взять дельта-функцию вида (x x0 ), где — малый амплитудный параметр.

Тогда, подставляя это выражение в (1.1.3), найдем Формула (1.1.5) позволяет свести операцию вычисления функциональной (вариационной) производной к обычному дифференцированию.

§ 1.2. Дифференцирования функционалов дифференцирования помощью (1.1.5) можно сформулировать все основные правила функционального дифференцирования:

Фигурирующая в последнем соотношении (1.2.3), F () — композиция функций F и, а F = F/ — производная внешней функций по аргументу.

Чтобы выполнять функциональное дифференцирование на практике, правила (1.2.1)–(1.2.3) необходимо дополнить важной формулой которую легко получить из (1.1.5), полагая F [u] = u(x).

Формулы (1.2.1)–(1.2.4) позволяют дифференцировать большинство функционалов, с которыми приходится встречаться на практике. Из (1.1.5) вытекает важное следствие, что операция функционального дифференцирования по переменной u коммутирует с любыми другими операторами, от u не зависящими. В частности, можно производить функциональное дифференцирование под знаком интеграла или производной. Это обстоятельство позволяет ввести самостоятельный символ функционального дифференцирования /u, правила обращения с которым во многом аналогичны правилам обращения с дифференциальным оператором /x x.

В качестве примера приведем правило дифференцирования сложного функционала F [f1, f2,..., fN ], в котором внутренние функции f1,..., fN, в свою очередь, функционально зависят от функции (x).

Это правило имеет вид и внешне напоминает правило обычного дифференцирования сложной функции.

В выражении (1.2.5) и далее, если не оговорено иное, по повторяющемуся дискретному индексу подразумевается суммирование, а штрихованные полевые переменные обозначают зависимость от штрихованных пространственных координат:

u = u(x ), u = u(x ).

локально-полевого функционала локально-полевого функционала (1.1.1) по полю u(x). Используя формулу (1.2.5), получим Откуда, учитывая соотношения:

найдем2, что Формулу (1.2.4) можно рассматривать как бесконечномерный (непрерывный) аналог формулы ui /uj = j, где ui — дискретные величины, а j — дельта-символ Кронекера, который определяi i ется правилом: j = 1, если i = j, и j = 0, если i = j.

Другой вывод формулы (1.2.6) может быть получен перестановкой оператора /u и интеграла с последующим вычислением производной h /u по правилу (1.2.3).

1.2. Дифференцирование функционалов Линейный оператор E Дифференцирование функционалов и частную производные.

Для примера рассмотрим сложный функционал F [u1, u2,..., uN ; t], зависящий от времени t как явно, так и неявно, т. е. через функциональную зависимость от внутренних функций ui (t, x), (i = 1, 2,..., N ).

Приведем правило дифференцирования такого функционала по времени:

Здесь слева стоит полная производная dF/dt, а справа — частная производная по явному времени t F, которая берется при фиксированных внутренних функциях ui. Частные производные по времени для внутренних функций t ui дают интегральный вклад, определяемый вариационными производными.

Для доказательства формулы (1.2.7) ограничимся случаем, когда функционал F [u; t] зависит только от одной внутренней функции u(t, x), и рассмотрим полную производную по параметру t как предел отношения Этот предел можно представить в виде суммы двух пределов:

Первый предел есть не что иное, как частная производная t F, т. е. производная по явному времени при фиксированной функции u(t, x). Второй легко вычисляется, если заметить, что, по существу, под знаком предела в числителе стоит вариация функционала F.

Так как эта вариация является откликом на вариацию функции u = tt u, она определяется функциональной производной и, согласно (1.1.4), выражается через интеграл Таким образом, после вычисления пределов из (1.2.8) для полной производной окончательно получим формулу обобщение которой на более общий случай (1.2.7) представляется очевидным.

Разложение функционалов что такого рода степенные разложения лежат в основе различных версий теории возмущений, к которым часто прибегают, например, при исследовании линейных и нелинейных проблем волновой и вихревой динамики (см. части V и VI). В частности, для простейшего однополевого функционала F [u + ], рассматривая поле u(x) как основное состояние, а поле (x) как возмущение, можно написать следующее разложение:

Если правую часть представить в операторном виде:

и заметить, что в фигурных скобках имеет место степенное разложение экспоненты, формулу (1.2.9) можно переписать в более изящном виде По аналогии с псевдодифференциальным оператором ehx, который действует на функцию (x) по правилу и сдвигает ее аргумент x на величину h, оператор exp иногда называют операdx (x) u(x) тором функционального сдвига.

ОСНОВНЫЕ ПРЕДСТАВЛЕНИЯ

О БЕСКОНЕЧНОМЕРНЫХ

ГАМИЛЬТОНОВЫХ СИСТЕМАХ

§ 2.1. Гамильтоновы системы простейшего типа Всякий раз, когда говорят о гамильтоновой структуре уравнений, описывающих динамику какой-либо непрерывной (бесконечномерной) системы, по существу, имеют в виду определенную форму записи этих уравнений.

Наиболее известна каноническая форма записи, когда система описывается четным числом уравнений вида относительно двух групп полевых переменных — обобщенных координат q 1,..., q N и обобщенных импульсов p1,..., pN, которые являются функциями времени t и пространственных координат x. Переменные, обеспечивающие такую структуру уравнений, называются каноническими, а величина H, стоящая под знаком вариационных производных и функционально зависящая от этих переменных, — гамильтонианом.

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

Иными словами, возможны и другие многочисленные версии гамильтоновых систем, которые известны как неканонические. Например, в приложениях часто приходится сталкиваться с проблемами, которые приводят к уравнениям вида где x — одна из пространственных координат в выделенном направлении, а ui — полевые переменные, которые зависят от времени t, x, но, вообще говоря, могут зависеть и от других пространственных координат.

В частности, если в качестве гамильтониана в (2.1.2) взять по очереди выражения:

28 Глава 2. Основные представления о гамильтоновых системах то из (2.1.2) получим соответственно два известных уравнения:

Если в первом уравнении, которое известно как уравнение Кортевега—де Вриза (КдВ), динамическая переменная u — функция одной пространственной переменной x, т. е. u = u (x, t), то во втором, которое известно как уравнение Кадомцева— Петвиашвили (КП), u является функцией двух пространственных переменных x, y, т. е. u = u (x, y, t).

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

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

В основе одного из способов современного определения таких систем лежит понятие функциональной скобки Пуассона (см., например: [74, 137, 139, 140]).

Скобка Пуассона для двух величин F и G на фазовом пространстве, состоящем из совокупности функционально независимых друг от друга полевых переменных ui (t, x), i = 1, 2,..., N, задается тензорным полем U ij [u; x, x ] по правилу Запись U ij [u; x, x ] отражает то обстоятельство, что ядро интеграла — так называемый тензор Пуассона U ij, вообще говоря, может являться не только функцией координат x и x, но и функционалом полей ui.

Подстановка F = ui (t, x), G = uj (t, x ) в (2.2.1) приводит к соотношению которое выражает очевидный факт: тензор Пуассона может быть задан через так называемые фундаментальные скобки Пуассона {ui (x), uj (x )} от базисных функций ui, образующих фазовое пространство.

Используя (2.2.2), можно переписать (2.2.1) в виде 2.3. Гамильтоновские системы общего вида По определению, скобки Пуассона должны обладать двумя основными свойствами. Во-первых, свойством антисимметричности Во-вторых, скобки Пуассона должны удовлетворить тождеству или свойству Якоби где сокращенная запись «ц.п.» обозначает слагаемые, полученные из первого циклической перестановкой по индексам и аргументам.

Чтобы подчеркнуть особый геометрический смысл свойства (2.2.4), иногда рассматривают так называемый тензор Якоби [225] T ijk (x, x, x ) = {ui (x), {uj (x ), uk (x )}} + ц.п. = который играет роль аналога тензора кривизны (см. § 3.2), но для фазового пространства динамической системы. Тогда гамильтоновые системы, требующие соблюдения тождества Якоби, можно охарактеризовать условием которое можно интерпретировать как отсутствие «кривизны» в фазовом пространстве динамических переменных.

§ 2.3. Гамильтоновские системы общего вида Скобки Пуассона, заданные на фазовом пространстве полей u1,..., uN, позволяют определить так называемые гамильтоновские системы как системы, которые эволюционируют по закону где H — гамильтониан системы — величина, функционально зависящая от полей ui.

Отметим, что возможность формальной записи уравнений движения в интегральном виде (2.3.1) еще не означает гамильтоновости системы, если для тензорного поля U ij не выполняется свойство антисимметричности и не соблюдается тождество Якоби Нарушение хотя бы одного из этих свойств, что возможно при некоторых аппроксимациях, — явный признак потери консервативности, т. е. разрушения интегралов движения.

Если гамильтониан системы явно не зависит от времени и, следовательно, t H = 0, закон сохранения энергии автоматически следует из данной формулировки, так как равенство выполняется в силу свойства антисимметричности (2.2.3).

С точки зрения данных выше определений канонический гамильтоновский формализм (2.1.1) соответствует скобкам Пуассона вида где j — дельта-символ Кронекера — определяется правилом: j = 1, если i = j, и j = 0, если i = j.

В другом примере (2.1.2) гамильтонова структура определяется скобками Пуассона вида Оба примера являются тривиальными в том отношении, что соответствующее им тензорное поле U ij не зависит функционально от полевых переменных, и, следовательно, свойство (2.2.4) выполняется автоматически.

К числу нетривиальных, но сравнительно простых континуальных гамильтоновых систем относятся системы, у которых тензорное поле U ij [u; x, x ] линейно зависит от полевых переменных. К этому виду, например, принадлежит уравнение Ландау—Лифшица [126] Здесь знак «» — обозначает векторное произведение. Это уравнение встречается в физике твердого тела (см.: [110, 128]) и используется для описания динамики классического спинового поля n (x, t) в ферромагнетиках в приближении сплошной среды.

Скобки Пуассона, соответствующие структуре уравнения Ландау—Лифшица (2.3.4), определяются уравнением Напомним (см. более подробно по этому поводу текст на с. 25), что полная производная d/dt от функционала [u1,..., uN ; t], который может зависеть от времени t как явно, так и неявно через функции u1,..., uN, берется следующим образом где t — частная производная функционала по явному времени берется при фиксированных значениях u1,..., uN.

2.4. Локальные скобки Пуассона и их свойства где eikj — единичный антисимметричный тензор (символ Леви—Чивита) — задается следующим правилом: eikj = 0, если среди индексов ikj есть два одинаковых, и eikj = ±1, если все индексы различны. Причем знак «+» реализуется, если упорядоченная система индексов ikj образует четную круговую перестановку чисел 1, 2, 3, и знак «», если нечетную перестановку.

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

Как было показано В. И. Арнольдом [4], на фазовом пространстве поля это уравнение можно записать в виде где функционал является кинетической энергией жидкости, а векторное поле (H /) есть не что иное, как гидродинамическая скорость v.

Нетрудно убедиться, что уравнение (2.3.6) является гамильтоновым, а фундаментальная скобка Пуассона дается выражением Здесь и далее i = /xi — операторы дифференцирования2 по пространственным координатам (i = 1, 2, 3).

§ 2.4. Локальные скобки Пуассона и их свойства Перечисленные выше скобки (2.3.2), (2.3.3), (2.3.5), (2.3.8) — типичные представители так называемых локальных скобок Пуассона.

По определению [74], локальные скобки Пуассона общего вида представляются конечной суммой Эта сумма составлена из дельта-функции (x x ) и ее производных с тензорij ными коэффициентами Bk [ui ; x], которые зависят от полей ui и их производных По умолчанию, если что-то не оговорено особо, операторы дифференцирования действуют на все, что стоит справа.

32 Глава 2. Основные представления о гамильтоновых системах в точке x. По поводу других используемых обозначений предполагается, что k = = (k1, k2, k3 ) — мультииндекс, по которому в (2.4.1) подразумевается суммироkkk вание, k = 1 1 2 2 3 3 — оператор дифференцирования по пространственным координатам i = /xi, (i = 1, 2, 3). Суммирование в (2.4.1), как уже говорилось, является конечным и идет по повторяющемуся мультииндексу k так, что суммарный индекс |k| = k1 + k2 + k3 пробегает все целые числа от 0 до K. Целое число K называется порядком скобки.

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

Последнее условие, например, выполняется для скобок (2.4.1) и локальных гамильтонианов специального вида у которых плотность гамильтониана h является аналитической функцией, зависящей от x, ui и частных производных k ui ui, где k — мультииндекс, используемый в соответствии с установленным выше правилом k = 1 1 2 2 3 3.

Как легко показать, в этом случае гамильтоновы уравнения (2.3.1) приводят к системе уравнений в частных производных с конечным числом членов. Действительно, в соответствии с правилом (1.2.5), найдем, что вариационная производная локального функционала (2.4.2) определяется выражением Подставляя (2.4.3) в (2.3.1) и принимая во внимание (2.4.1), получим:

где s, так же, как и k, — мультииндекс.

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

Так, свойство антисимметричности приводит к условию Равенство (2.4.4) необходимо понимать в смысле обобщенных функций. Собирая подобные члены, имеющие одинаковый характер сингулярности в точке x = x, и приравнивая затем полученный результат нулю, можно показать, что (2.4.4) эквивалентно конечной системе линейных соотношений на функции Bk, их частные производные по x и на частные производные Bk по полям ui.

2.4. Локальные скобки Пуассона и их свойства Дополнительные ограничения на функции Bk накладывает тождество Якоби.

Подстановка (2.4.1) в (2.2.4) приводит к условиям где s, p и t — мультииндексы.

Аналогичные рассуждения приводят нас к заключению, что равенство (2.4.5) равносильно конечной системе уравнений, квадратичных относительно функций Bk и их частных производных по x и по ui. m Несмотря на ограничения (2.4.4), (2.4.5), в выборе структурных функций Bk остается значительный произвол. При построении скобок Пуассона этот произвол преодолевается либо в рамках классификационного подхода, который опирается на теоретико-групповой анализ структуры скобок Пуассона (см., например, [137, 140] и цитируемую там литературу), либо в рамках феноменологического подхода. Последний исходит из ряда фундаментальных физических принципов, которые накладывают дополнительные ограничения на скобки Пуассона (см., например, [211]).

Некоторые из этих принципов могут быть сформулированы как требование инвариантности гамильтоновских систем относительно тех или иных преобразований. В контексте данного раздела ограничимся обсуждением лишь одного наиболее тривиального следствия, связанного со свойством однородности координатного пространства. Более подробное обсуждение проблемы инвариантности гамильтоновых систем и вытекающих отсюда ограничений на скобки Пуассона отложим на будущее (см. главы 7,8,12).

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

Аналогичным образом можно рассмотреть ограничения на выбор структурij ных функций Bk, возникающие вследствие инвариантности относительно поворотов осей координат. Как известно, это свойство, обусловленное равноправием всех направлений в пространстве, называется изотропностью, а соответствующие скобки Пуассона — изотропными. Однако следует иметь в виду, что для этого в каждом конкретном случае необходимо знать, как ведут себя при поворотах осей координат полевые переменные ui, которые, вообще говоря, могут иметь не только скалярный или векторный, но и более сложный геометрический характер [34].

ГЕОМЕТРИЧЕСКИЕ СВОЙСТВА

ФАЗОВЫХ ПРОСТРАНСТВ.

СИМПЛЕКТИЧЕСКИЕ СТРУКТУРЫ

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

основанные на концепции дифференциальных 2-форм К подобным гамильтоновым системам можно прийти естественным образом.

Рассмотрим с этой целью гамильтоновы системы (2.3.1) с невырожденной скобкой Пуассона. Невырожденность скобки означает, что интегральное уравнение не имеет относительно функционала C, который не зависит от времени явно, никаких других решений, кроме C = const. В этом случае тензор Пуассона U ik имеет обратный ik, и они связаны друг с другом соотношением Воспользовавшись (3.1.1), легко установить, что в этом случае эволюционное уравнение (2.3.1) приводится к виду где тензорное поле ik (x, x ) обладает свойствами:

3.1. Гамильтоновы системы, основанные на концепции 2-форм Эти свойства являются прямым следствием соответствующих свойств (2.2.3), (2.2.4), которыми обладает скобка Пуассона.

Не вдаваясь глубоко в теорию внешних дифференциальных форм, укажем, что в конечномерном случае с антисимметричными тензорами ik связывают так называемые внешние дифференциальные 2-формы, на языке которых гамильтонов формализм излагается столь же последовательно, как и на языке скобок Пуассона. Для более детального знакомства с этим предметом можно рекомендовать специальную математическую литературу (см., например, [5, 75, 143, 181]).

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

Напомним, что в конечномерном случае внешняя дифференциальная 2-форма, соответствующая антисимметричному тензору ik, записывается как где знак «» обозначает внешнее произведение, а ui — независимые переменные.

В бесконечномерном случае, когда ui и ik являются непрерывными функциями пространственной координаты x, эта формула должна быть переписана в соответствии с установленным выше правилом как Чтобы подчеркнуть отличие от (3.1.4), выражения, подобные (3.1.5), называют функциональными 2-формами.

Можно показать, что свойство (3.1.3) является условием замкнутости этой формы, т. е. вытекает из требования равенства нулю ее внешнего функционального дифференциала Тензорное поле Sikl определяется выражением и возникает как результат процедуры антисимметризации для тензора ik /ul под знаком интеграла (3.1.6).

Легко убедиться, что т. е. тензор Sikl представляет собой ковариантный аналог тензора Якоби T ikl, определяемого выражением (2.2.5).

По аналогии со скобкой Пуассона, функциональная 2-форма называется невырожденной, если интегральное уравнение не имеет относительно k никаких других решений, кроме k = 0.

В соответствии с общепринятой в дифференциальной геометрии терминологией [5, 75, 143, 181], замкнутая невырожденная дифференциальная 2-форма (3.1.5) называется симплектической структурой, а фазовое пространство, на котором она определена, симплектическим.

Очевидно, что для дискретных систем симплектическое пространство обязательно имеет четную размерность, поскольку в противном случае det ik = 0.

На непрерывные (полевые) системы это утверждение не переносится, так как под размерностью фазового пространства нельзя понимать только дискретную компоненту размерности, которая определяется числом полей ui, перечисляемых дискретным индексом i = 1, 2,..., N, а нужно еще иметь в виду бесконечномерность самих полевых переменных ui.

Иллюстрацией однополевого симплектического пространства может служить фазовое пространство уравнения где (x) = 1 sign(x) при x = 0 и (x) = 0 при x = 0.

Как будет показано далее в § 4.2, уравнение (3.1.8) — просто иная форма записи хорошо знакомого уравнения КдВ. Соответствующая симплектическая форма имеет вид Замкнутость формы очевидна, так как /u = 0. Остается проверить невырожденность, которая означает, что уравнение не имеет нетривиальных решений. Действительно, после дифференцирования (3.1.10) по x, получим решение (x) = 0.

3.2. Аналогия метрических и симплектических пространств Замечательным свойством симплектических пространств дискретных систем является то, что подходящей заменой динамических переменных тензор ik, а, следовательно, и U ik могут быть одновременно приведены к каноническому виду (где E — единичная матрица) в по крайней мере малой окрестности любой точки фазового пространства [143]. Переменные, в которых тензоры ik и U ik принимают специальный вид (3.2.1), как известно, называются каноническими.

Этот факт составляет содержание известной теоремы Дарбу и имеет наглядную геометрическую интерпретацию, основанную на качественной аналогии между геометрией метрических пространств и геометрией фазовых пространств гамильтоновых систем [181, 225, 273].

Чтобы провести такую аналогию, напомним вначале некоторые сведения из метрической геометрии (см., например, [157]). Хорошо известно, что метрическое пространство характеризуется симметричным невырожденным метрическим тензором gik, который полностью определяет всю геометрию этого пространства.

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

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

Вернемся к рассмотрению гамильтоновых систем. Как уже говорилось, динамика таких систем разворачивается в кососимметричных фазовых пространствах, геометрические свойства которых исчерпывающим образом характеризуются тензором ik или U ik.

Тензорный характер ik и U ik следует из ковариантности уравнений (2.3.1), описывающих движение гамильтоновых систем, при преобразованиях динамических переменных.

Чтобы убедиться в этом, рассмотрим взаимооднозначное преобразование к новым полевым переменным Легко показать, что в терминах полей ai уравнение (2.3.1) сохраняет свою гамильтонову структуру где новый тензор Пуассона U ik получен из старого U ik как результат преобразования Оказывается, что канонические переменные для описания гамильтоновых систем играют ту же роль, что и декартовы координаты в метрической геометрии.

Подобно декартовым координатам канонические переменные без изменения размерности фазового пространства могут быть введены только при соблюдении двух условий: невырожденности кососимметричной метрики (det ik = 0 или det U ik = 0) и обращении в нуль тензора T ikl. Тензор Якоби T ikl служит грубой аналогией тензора кривизны в метрическом пространстве [225].

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

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

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

Для системы динамических уравнений этот метод реализуется следующим образом.

Полагая переменные ui каноническими координатами, введем сопряженные полевые переменные pi, играющие роль канонических импульсов, и построим «гамильтониан»:

Тогда, с одной стороны, канонические уравнения (2.1.1), соответствующие такому «гамильтониану», воспроизводят исходную систему уравнений (3.3.1), а с другой — приводят к системе ассоциированных уравнений решения которой никак не влияют на решения (3.3.1).

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

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

Для реальных физических систем это означает, что гамильтониан не может содержать искусственные переменные типа pi, которые не выражались бы через исходные (физические) переменные ui. Такое свойство, тесно связанное с так называемой калибровочной инвариантностью теории, имеет глубокий физический смысл, и более подробно будет обсуждаться в § 6.3.

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

Этот способ продемонстрируем на уже знакомом уравнении КдВ для которого скобка Пуассона имеет вид Что значит ввести канонические переменные? Очевидно, один из возможных ответов на этот вопрос состоит в том, чтобы суметь выразить переменную u через две другие q и p, которые удовлетворяют следующим коммутационным соотношениям (условиям каноничности):

Таким образом, проблема сводится к поиску такой функциональной зависимости u[q, p], где q и p — канонические переменные, которая бы воспроизводила соотношение (3.3.3).

Такого типа зависимости известны как представления или преобразования Клебша, а канонические переменные q и p называют потенциалами или переменными Клебша. Исторически эти названия связаны с А. Клебшем [202], который впервые использовал подобное преобразование для гидродинамической скорости в несжимаемой жидкости 1.

В 1859 г. А. Клебш нашел, что разложение поля гидродинамической скорости позволяет переписать уравнения движения идеальной однородной несжимаемой жидкости в замкнутом виде как Нетрудно догадаться (подсказкой служит вид правой части (3.3.3)), что соответствующее представление Клебша для u[q, p] должно быть билинейным по переменным q, p, и оператору x, т. е.

где 1, 2, 1, 2 — некоторые константы.

Подстановка (3.3.5) в (3.3.3) с использованием (3.3.4) приводит к одному соотношению на четыре коэффициента 1, 2, 1, 2. Поэтому (3.3.6) ограничивает, но не устраняет полностью произвол в выборе этих констант.

Воспользовавшись оставшимся произволом в выборе констант и положив 1 = 1, 2 = 1/2, 2 = 1 = 0, можно привести (3.3.5) к более простому виду Отметим, что тот же самый результат мог быть достигнут сразу с помощью аффинного преобразования которое при условии (3.3.6) является каноническим, т. е. вводит новую координату Q и сопряженный импульс P.

Однако, кроме параметрического произвола, который удается таким образом устранить, остается еще функциональный, связанный с неоднозначностью выбора канонического базиса q и p при фиксированном «физическом» поле u.

Как увидим далее (см. более подробное обсуждение в § 6.3), такая ситуация весьма типична и отражает потенциальную возможность канонических калибровочных преобразований, оставляющих инвариантными все физические величины теории, а, следовательно, существование специфических законов сохранения. С геометрической точки зрения эти законы сохранения фиксируют поверхность, на которой сосредоточено действительное движение системы в симплектическом пространстве q, p.

с наложенным на скалярные переменные (потенциалы), µ, и дополнительным ограничением которое имеет место в силу условия несжимаемости · v = 0.

Хотя Клебш понимал, что полученная система эквивалентна некоторой вариационной задаче, о канонической сопряженности потенциалов и µ он не догадывался. Этот факт стал осознан позже (первый шаг был сделан вероятно Г. Лэмбом, 1932 г.), сразу после того, как уравнения Клебша были переписаны в виде где в качестве гамильтониана взята полная энергия несжимаемой жидкости

ВАРИАЦИОННЫЙ ПРИНЦИП

НАИМЕНЬШЕГО ДЕЙСТВИЯ.

ВЫРОЖДЕННЫЕ

ГАМИЛЬТОНОВСКИЕ СИСТЕМЫ

§ 4.1. Обобщение вариационного принципа основанные на дифференциальных 2-формах Хорошо известно (см., например, работы [9, 13, 33, 120]), что уравнения движения для весьма широкого класса физических систем следуют из вариационного принципа наименьшего действия. В самом общем виде этот принцип утверждает, что существует такой интеграл S, называемый действием, где L — лагранжиан динамической системы, что каждое действительное движение, удовлетворяющее уравнениям движения с фиксированными начальным, в момент t1, и конечным, в момент t2, состояниями, реализует экстремум этого функционала S:

Ограничиваясь рассмотрением только гамильтоновых систем, вначале выясним существование принципа наименьшего действия применительно к самому простому случаю — каноническим гамильтоновым системам (2.1.1).

В этом случае принцип (4.1.1), (4.1.2) формулируется в следующем виде:

с дополнительным условием:

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

Обобщение вариационного принципа (4.1.3), (4.1.4) на неканонические гамильтоновы системы можно получить, если в качестве лагранжиана рассмотреть величину а в качестве дополнительных условий взять равенства полагая при этом, что величины Ai u1,..., uN ; x — некоторые функционалы, зависящие явно, быть может, от координаты x, но только не от времени t.

Легко проверить, что уравнения которые были получены в рамках гамильтоновского формализма, основанного на концепции дифференциальных 2-форм, в точности следуют из обобщенного вариационного принципа, если тензор ik определить как Следует особо подчеркнуть, что правая сторона (4.1.8) автоматически удовлетворяет необходимым условиям антисимметричности (3.1.2) и замкнутости (3.1.3). Поэтому на языке внешних дифференциальных форм соотношение (4.1.8) означает, что где 1-форма определяется выражением В этом случае принято говорить, что форма является точной.

Таким образом, для существования вариационного принципа наименьшего действия у гамильтоновых систем необходимо, чтобы соответствующая замкнутая 2-форма была точной глобально, т. е. во всей области фазового пространства, на котором определены решения уравнений (4.1.7).

Является ли замкнутая форма точной глобально (локально это всегда так) зависит от топологических свойств фазового пространства. Препятствием для этого может оказаться неодносвязность. Если же область фазового пространства односвязна, т. е., если любой замкнутый путь или петля стягивается по этой области в точку, условие замкнутости (3.1.3) в данной области интегрируется соотношением (4.1.8).

Изучением топологических свойств, определяющих связь между замкнутыми и точными формами, занимается теория когомологий, краткий экскурс в которую можно найти в [133, 187].

4.1. Обобщение вариационного принципа Модель вихревой нити (более подробное обсуждение см. в § 28.5 и § 28.6) — один из примеров неканонических гамильтоновских систем, которые имеют вариационный принцип наименьшего действия.

Уравнения движения вихревой нити (см., например, [131]) формулируются в виде с гамильтонианом Здесь x(s) — пространственные координаты вихревой нити, параметрически заданной с помощью натурального параметра s; s = /s; H — гамильтониан, а — константа, характеризующая интенсивность вихревой нити.

Легко проверить, что в фазовом пространстве динамических переменных xi (s, t), (i = 1, 2, 3) уравнение (4.1.11) может быть переписано как (4.1.7) с тензорным полем ik вида Зная ik, из (4.1.8) можно найти соответствующие «потенциалы» Ai :

и построить лагранжиан (4.1.5).

В результате, для модели вихревой нити, согласно (4.1.1), можно найти интеграл действия, который описывается выражением:

Покажем, что рассматриваемая гамильтонова система (4.1.11), описывающая движение вихревой нити, вырождена. Действительно, согласно уравнению (3.1.7), которое в данном случае принимает вид:

и подстановкой (4.1.12) сводится к уравнению:

существует нетривиальное решение = f s x, где f (s) — произвольная скалярная функция натурального параметра s.

Этот пример поучителен в том отношении, что демонстрирует существование вариационного принципа не только для невырожденных гамильтоновых систем, но и для вырожденных — тех, в основе которых лежит вырожденная форма.

с вырожденными скобками Пуассона.

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

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

Как уже отмечалось в § 3.1, вырожденность скобки Пуассона означает, что интегральное уравнение обладает некоторым набором нетривиальных и независимых решений Cn, где целочисленный индекс n пробегает значения n = 1,..., K.

Независящие явно от времени функционалы Cn называются аннуляторами скобок Пуассона, или функционалами Казимира, или просто казимирами.

Из (4.2.1) следует, что казимир Cn коммутирует с любым функционалом F, заданным на фазовом пространстве переменных ui (i = 1,..., N ), т. е. их взаимная скобка Пуассона обращается в нуль В частности, в качестве F может выступать гамильтониан системы H. Поэтому казимиры обладают свойством и, следовательно, являются сохраняющимися величинами или инвариантами движения Так как в силу (4.2.1) существование сохраняющихся функционалов Cn целиком определяется структурой скобки Пуассона, следует подчеркнуть, что в отличие от обычных интегралов движения казимиры являются инвариантами движения для систем с любым гамильтонианом в данном классе скобок Пуассона.

С геометрической точки зрения совокупность K независимых равенств (4.2.2) выделяет в исходном фазовом пространстве K-параметрическую поверхность, на которой и сосредоточено действительное движение гамильтоновой системы.

Роль параметров, фиксирующих поверхность, играют константы в правой части (4.2.2). Каждому набору из K параметров соответствует своя поверхность, на которой вырожденная скобка Пуассона становится невырожденной. Геометрически это означает, что первоначальное фазовое пространство расслаивается на семейство симплектических поверхностей — фазовых подпространств, обладающих в силу теоремы Дарбу каноническим базисом.

4.2. Гамильтоновские системы с вырожденными скобками Пуассона Сказанное выше проиллюстрируем конкретными примерами.

Простейшим примером однополевой вырожденной гамильтоновой системы является уравнение КдВ со скобкой Пуассона В этом случае интегральное уравнение (4.2.1) принимает вид из которого следует, что Таким образом, в фазовом пространстве исходной системы, состоящем из одного единственного поля u, условие C = const выделяет подпространство — множество функций u(x, t), достаточно быстро убывающих при |x| так, чтобы существовал интеграл (4.2.5). На таком подпространстве скобка (4.2.4) становится уже невырожденной ввиду обратимости оператора x. Можно показать, что в этом случае уравнение (4.2.3) приводится к виду а соответствующая невырожденная форма оказывается точной, т. е. удовлетворяет соотношениям (4.1.9), (4.1.10).

Потенциал формы определяется выражением Следовательно, система (4.2.6) имеет вариационный принцип (4.1.5), (4.1.6), основанный на лагранжиане В качестве второго примера рассмотрим динамику спинового поля со скобкой Пуассона Этот пример уже рассматривался нами в § 2.3 (см. уравнение движения (2.3.4) и скобку Пуассона (2.3.5) этого раздела).

Инварианты Казимира для спиновой системы находятся из уравнения имеющего элементарное решение C = n2, которое в соответствии с условием (4.2.2) приводит к ограничению n2 = const.

По существу, это означает, что движение спиновой системы в фазовом пространстве векторного поля n происходит по поверхности сферы так, что фазовое пространство расслаивается на сферы (рис. 4.1).

Выбирая без ограничения общности радиус сферы равным единице, что соответствует C = 1, в описании можно оставить только две переменные, например, n1 и n2, которые образуют симплектический базис. Тогда третья компонента n исключается подстановкой Один из способов введения канонических переменных на симплектическом базисе n1, n2 можно осуществить с помощью преобразования:

где a и a — комплексные канонически сопряженные переменные (символ «»

обозначает комплексное сопряжение), i — мнимая единица.

Другой способ решить эту проблему состоит в том, чтобы представить единичный вектор спина n в сферических координатах:

4.2. Гамильтоновские системы с вырожденными скобками Пуассона Как показано на рис. 4.1, динамические переменные и, которые описывают динамику спинового поля, имеют простой геометрический смысл — полярного и азимутального угла, соответственно.

Переменные и так же, как и n1 и n2, образуют симплектический базис и могут быть выражены в терминах канонически сопряженных переменных a, a :

В терминах канонически сопряженных переменных a, a уравнение движения (4.2.8) принимает вид где гамильтониан H [a, a ] для канонического описания получается из исходного гамильтониана H [n] заменой переменных (4.2.9), (4.2.10).

Преобразование (4.2.10) от векторного поля n к канонически сопряженным переменным a, a является классическим аналогом преобразования Гольдштейна—Примакова [2, 59] и впервые использовалось в работе [276] для анализа нелинейных процессов в ферромагнетиках.

В качестве последнего примера рассмотрим уравнение движения несжимаемой однородной жидкости в форме Арнольда Этот пример также уже обсуждался в § 2.3 (см. уравнение движения (2.3.6) и гамильтониан (2.3.7) из этого раздела).

Скобка Пуассона соответствующая этой модели, является вырожденной, поскольку уравнение имеет нетривиальное решение Этот инвариант, как известно [250], несет топологический смысл, а именно: характеризует запутанность вихревых линий и называется спиральностью.

Существование топологических инвариантов типично для тех динамических систем, которые эволюционируют так, что топологическая классификация начальных условий переносится без изменений на решения уравнений в произвольный момент времени [185].

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

Пусть завихренность жидкости, связанная со скоростью v соотношением = v, сконцентрирована на замкнутом контуре x, который задан в параметрической виде с помощью параметра s.

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

Интегрируя равенство (4.2.13), получим контурный интеграл где в силу соленоидальности поля вихря · = 0, функция a(s) удовлетворяет условию Постоянная называется циркуляцией или интенсивностью и является индивидуальной характеристикой нити.

Рассмотрим теперь две замкнутые вихревые нити x1 и x2 с интенсивностями 1 и 2, соответственно.

Воспользовавшись представлением (4.2.14), (4.2.15), согласно принципу суперпозиции, находим, что полное поле завихренности для двух вихревых нитей может быть представлено как Подставляя (4.2.16) в (4.2.12) и интегрируя по x, найдем Интегралы в (4.2.17) есть не что иное, как циркуляции скорости v по замкнутым жидким контурам, в качестве которых служат сами вихревые нити.

В соответствии с теоремой Стокса [132] можно показать, что где m — число зацеплений или число витков, которые делает одна нить вокруг другой. Знак m зависит от взаимной ориентации этих вихревых нитей (рис. 4.2).

4.2. Гамильтоновские системы с вырожденными скобками Пуассона Рис. 4.2. Характер зацепления двух замкнутых вихревых нитей x1 и x2. Выбор знака в вариантах б и в определяется относительной ориентацией вихревых нитей Принято считать, что m является положительным, если направления обхода по замкнутым вихревым нитям x1 и x2 составляют правовинтовую систему, как на рис. 4.2, в, и отрицательным, если направления обхода по этим нитям составляют левовинтовую систему, как на рис. 4.2, б. Для практического определения знака существует простое «правило буравчика». Согласно этому правилу, если при повороте ручки буравчика (правой нарезки) по направлению обхода одного контура острие буравчика пойдет по направлению обхода другого контура, то выбирается знак плюс, если — против, то выбирается знак минус.

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

ПРЕОБРАЗОВАНИЕ

ДИНАМИЧЕСКИХ ПЕРЕМЕННЫХ

В ГАМИЛЬТОНОВЫХ СИСТЕМАХ

При изучении гамильтоновых систем бросается в глаза очень важная их черта — структурная жесткость уравнений, т. е. сохранение свойств (2.2.3), (2.2.4):

или свойств (3.1.2), (3.1.3):

при различных преобразованиях переменных типа (3.2.2):

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

В методическом отношении гамильтонов формализм представляет собой образец структурного подхода, суть которого может быть сформулирована словами Н. Бурбаки [21, с. 104]: «Структуры являются орудиями математика, каждый раз, когда он замечает, что между элементами, изучаемыми им, имеют место отношения, удовлетворяющие аксиомам структуры определенного типа, он сразу может воспользоваться всем арсеналом общих теорем, относящихся к структурам этого типа, тогда как раньше он должен был бы мучительно выковывать сам средства, необходимые для того, чтобы штурмовать рассматриваемую проблему, причем их мощность зависела бы от его личного таланта, и они были бы отягчены часто излишне стеснительными предположениями, обусловленными особенностями изучаемой проблемы».

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

Это обстоятельство следует иметь в виду всегда, когда существует альтернатива в выборе гамильтоновой формулировки.

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

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

Подобного рода примером в гидродинамике была попытка по аналогии с двумерными точечными вихрями ввести трехмерные — «вортоны» [138], которые бы служили базовыми элементами для построения конечномерных моделей, описывающих трехмерные вихревые течения в идеальной жидкости. Во всех таких случаях потерю консервативности легко объяснить, если обратить внимание на то, что свойство Якоби (2.2.4) или (3.1.3) может разрушаться, например, при конечно-разностных аппроксимациях.

Особое значение это замечание имеет в связи с широким использованием в современной теоретической и вычислительной физике дискретных моделей, возникающих из непрерывных аналогов как результат их точной или приближенной конечномерной редукции (см., по этому поводу работы [1, 14, 74, 91, 171, 172] и цитируемую там литературу).

В этом свете каноническую формулировку следует рассматривать как один из наиболее развитых в методическом отношении вариантов гамильтонова формализма, не только отвечающего требованию структурной простоты, но и имеющего в своем арсенале целый ряд стандартных подходов к решению разнообразных проблем гидродинамики, физики плазмы и т.д. [1, 84, 87–89, 92, 156, 171, 172, 273].

§ 5.1. Канонические преобразования С канонической формулировкой связан важный для приложений класс преобразований, который не изменяет структуру скобки Пуассона, а, следовательно, и структуру уравнений (2.1.1). Такие преобразования, называемые каноническими, подробно описаны в любом курсе классической механики для конечномерных систем (см., например, [5, 28, 33]). Поэтому результаты, изложенные ниже, следует считать своего рода обобщением уже известных на континуальные системы.

Как уже отмечалось в § 4.1, канонические уравнения (2.1.1) можно сформулировать в виде вариационного принципа наименьшего действия с дополнительным условием где t1 и t2 — начальный и конечный момент времени, соответственно.

Будем интересоваться лишь такими преобразованиями:

при которых новые переменные Qi и Pi опять являются каноническими, и, следовательно, имеют место уравнения:

где H — новый гамильтониан системы.

Уравнения (5.1.3), так же, как и (2.1.1), можно сформулировать в виде принципа наименьшего действия с дополнительным условием Для того чтобы оба принципа (5.1.1) и (5.1.4) были эквивалентны, необходимо и достаточно выполнение двух требований.

Во-первых, подынтегральные выражения, стоящие в скобках (5.1.1) и (5.1.4), могут отличаться лишь на полную производную по времени от некоторого функционала 1 :

где может зависеть от времени как явно, так и неявно (через функциональную зависимость по полевым переменным q i, pi, Qi, Pi ).

Доказательство этого утверждения для континуальных систем совершенно аналогично доказательству для дискретных систем (см., например, [66, 120]).

5.1. Канонические преобразования и производящие функционалы Во-вторых, будем полагать, что фиксированы не только начальные и конечные координаты, но и импульсы:

Всякое каноническое преобразование характеризуется своим функционалом [q i, pi, Qi, Pi ; t], который называется производящим, поскольку с его помощью находят искомое преобразование. Соответствующая процедура формулируется следующим образом.

Найдем выражение для полной производной d/dt в развернутом виде. Учитывая, что может зависеть от времени как явно, так и неявно (через функциональную зависимость по полевым переменным q i, pi, Qi, Pi ), согласно правилу дифференцирования (1.2.7), получим где t — частная производная производящего функционала по явному времени, которая берется при фиксированных q i, pi, Qi, Pi.

После подстановки (5.1.6) в (5.1.5) и перегруппировки подобных членов под знаком интеграла, найдем Очевидно, что для реализующихся динамических процессов производные t q i, t pi, t Qi, t Pi не равны тождественно нулю. Поэтому, приравнивая нулю коэффициенты при этих производных, приходим к формулам:

определяющим каноническое преобразование.

Рассмотрим четыре наиболее распространенных вида канонических преобразований.

Преобразования первого вида В этом случае, согласно (5.1.7), находим:

Разрешая первые две системы равенств относительно Qi и Pi, каноническое преобразование можно сформулировать в виде (5.1.2).

Преобразования второго вида После подстановки (5.1.8) в (5.1.7), находим:

Преобразования третьего вида После подстановки (5.1.10) в (5.1.7) получим:

которое, в соответствии с (5.1.7), приводит к следующим формулам:

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

§ 5.2. Канонические точечные преобразования Очень часто в приложениях возникает необходимость совершать так называемые точечные преобразования, при которых старые координаты выражаются только через новые координаты 5.3. Канонические переменные для поверхностных волн Таким каноническим преобразованиям отвечает производящий функционал третьего вида (5.1.10), для которого Следовательно, в соответствии с (5.1.11) новые импульсы Pi и новый гамильтониан H определяются как Аналогичное точечное преобразование для импульсов задается производящим функционалом и, согласно (5.1.9), приводит к преобразованию Следует иметь в виду, что в каноническом формализме различие между импульсами и координатами чисто номенклатурное, т. е. всегда с помощью производящего функционала можно произвести преобразование которое осуществляет взаимное переименование координат и импульсов.

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

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

Напомним вначале один исторически предшествующий работе [81] промежуточный результат [194], согласно которому при потенциальных движениях сжимаемой идеальной безграничной жидкости канонической координатой служит плотность, а импульсом — гидродинамический потенциал.

Рис. 5.1. Модель двухслойной без- снизу (рис. 5.1).

вихревой жидкости с контактной ющие самоиндуцированную эволюцию контактной поверхности, является использование функции Хевисайда, которая обладает свойством С помощью -функции распределение плотности и гидродинамического потенциала в такой модели можно записать в виде разложений:

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

Во-первых, такая запись в математически корректной форме учитывает существование контактной поверхности, а во-вторых, открывает возможность для прямого использования сформулированного выше правила, согласно которому для данной модели роль канонической координаты играет плотность (5.3.1), а импульса — гидродинамический потенциал (5.3.2).

Рассмотрим теперь соотношение (5.3.1) как точечное преобразование к новой обобщенной координате, роль которой играет форма поверхности. Отметим попутно, что при этом размерность задачи снижается на единицу.

Согласно формулам (5.2.1) и (5.2.2), такому преобразованию соответствует производящий функционал поэтому новый импульс определяется как 5.4. Преобразование канонических переменных при замене координат В случае свободной границы раздела, когда + = 0, в соответствии с (5.3.3) канонический импульс определяется как Полученный результат совпадает с результатом [81] с точностью до постоянного множителя, который без ограничения общности можно положить равным единице или убрать с помощью простейшего канонического неунивалентного преобразования где H — исходный гамильтониан модели, в качестве которого обычно берут полную энергию.



Pages:   || 2 | 3 | 4 | 5 |   ...   | 8 |
Похожие работы:

«ИННОВАЦИОННО-ОРИЕНТИРОВАННАЯ ПОДГОТОВКА ИНЖЕНЕРНЫХ, НАУЧНЫХ И НАУЧНО-ПЕДАГОГИЧЕСКИХ КАДРОВ С.И. ДВОРЕЦКИЙ, Е.И. МУРАТОВА, И.В. ФЁДОРОВ ИННОВАЦИОННО-ОРИЕНТИРОВАННАЯ ПОДГОТОВКА ИНЖЕНЕРНЫХ, НАУЧНЫХ И НАУЧНО-ПЕДАГОГИЧЕСКИХ КАДРОВ ИЗДАТЕЛЬСТВО ТГТУ Министерство образования и науки Российской Федерации ГОУ ВПО Тамбовский государственный технический университет С.И. ДВОРЕЦКИЙ, Е.И. МУРАТОВА, И.В. ФЁДОРОВ ИННОВАЦИОННО-ОРИЕНТИРОВАННАЯ ПОДГОТОВКА ИНЖЕНЕРНЫХ, НАУЧНЫХ И НАУЧНО-ПЕДАГОГИЧЕСКИХ КАДРОВ...»

«Министерство образования Республики Беларусь УЧРЕЖДЕНИЕ ОБРАЗОВАНИЯ ГРОДНЕНСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ИМЕНИ ЯНКИ КУПАЛЫ И.И.Веленто ПРОБЛЕМЫ МАКРОПРАВОВОГО РЕГУЛИРОВАНИЯ ОТНОШЕНИЙ СОБСТВЕННОСТИ В РЕСПУБЛИКЕ БЕЛАРУСЬ И РОССИЙСКОЙ ФЕДЕРАЦИИ Монография Гродно 2003 УДК 347.2/.3 ББК 67.623 В27 Рецензенты: канд. юрид. наук, доц. В.Н. Годунов; д-р юрид. наук, проф. М.Г. Пронина. Научный консультант д-р юрид. наук, проф. А.А.Головко. Рекомендовано Советом гуманитарного факультета ГрГУ им....»

«МЕДИЦИНСКАЯ АКАДЕМИЯ ПОСЛЕДИПЛОМНОГО ОБРАЗОВАНИЯ В. В. Афанасьев, И. Ю. Лукьянова Особенности применения цитофлавина в современной клинической практике Санкт-Петербург 2010 Содержание ББК *** УДК *** Список сокращений.......................................... 4 Афанасьев В. В., Лукьянова И. Ю. Особенности применения ци тофлавина в современной клинической практике. — СПб., 2010. — 80 с. Введение.................................»

«Особо охраняемые природные территории УДК 634.23:581.16(470) ОСОБО ОХРАНЯЕМЫЕ РАСТЕНИЯ САМАРСКОЙ ОБЛАСТИ КАК РЕЗЕРВАТНЫЙ РЕСУРС ХОЗЯЙСТВЕННО-ЦЕННЫХ ВИДОВ © 2013 С.В. Саксонов, С.А. Сенатор Институт экологии Волжского бассейна РАН, Тольятти Поступила в редакцию 17.05.2013 Проведен анализ группы раритетных видов Самарской области по хозяйственно-ценным группам. Ключевые слова: редкие растения, Самарская область, флористические ресурсы Ботаническое ресурсоведение – важное на- важная группа...»

«А.Н. КОЛЕСНИЧЕНКО Международные транспортные отношения Никакие крепости не заменят путей сообщения. Петр Столыпин из речи на III Думе О стратегическом значении транспорта Общество сохранения литературного наследия Москва 2013 УДК 338.47+351.815 ББК 65.37-81+67.932.112 К60 Колесниченко, Анатолий Николаевич. Международные транспортные отношения / А.Н. Колесниченко. – М.: О-во сохранения лит. наследия, 2013. – 216 с.: ил. ISBN 978-5-902484-64-6. Агентство CIP РГБ Развитие производительных...»

«Министерство здравоохранения Российской Федерации Тихоокеанский государственный медицинский университет В.А. Дубинкин А.А. Тушков Факторы агрессии и медицина катастроф Монография Владивосток Издательский дом Дальневосточного федерального университета 2013 1 УДК 327:614.8 ББК 66.4(0):68.69 Д79 Рецензенты: Куксов Г.М., начальник медико-санитарной части УФСБ России по Приморскому краю, полковник, кандидат медицинских наук; Партин А.П., главный врач Центра медицины катастроф Приморского края;...»

«В.Н. Ш кунов Где волны Инзы плещут. Очерки истории Инзенского района Ульяновской области Ульяновск, 2012 УДК 908 (470) ББК 63.3 (2Рос=Ульян.) Ш 67 Рецензенты: доктор исторических наук, профессор И.А. Чуканов (Ульяновск) доктор исторических наук, профессор А.И. Репинецкий (Самара) Шкунов, В.Н. Ш 67 Где волны Инзы плещут.: Очерки истории Инзенского района Ульяновской области: моногр. / В.Н. Шкунов. - ОАО Первая Образцовая типография, филиал УЛЬЯНОВСКИЙ ДОМ ПЕЧАТИ, 2012. с. ISBN 978-5-98585-07-03...»

«Межрегиональные исследования в общественных науках Министерство образования и науки Российской Федерации ИНОЦЕНТР (Информация. Наука. Образование) Институт имени Кеннана Центра Вудро Вильсона (США) Корпорация Карнеги в Нью-Йорке (США) Фонд Джона Д. и Кэтрин Т. МакАртуров (США) Данное издание осуществлено в рамках программы Межрегиональные исследования в общественных науках, реализуемой совместно Министерством образования и науки РФ, ИНОЦЕНТРом (Информация. Наука. Образование) и Институтом имени...»

«КАЗАХСТАНСКИЙ ИНСТИТУТ СТРАТЕГИЧЕСКИХ ИССЛЕДОВАНИЙ ПРИ ПРЕЗИДЕНТЕ РЕСПУБЛИКИ КАЗАХСТАН МУРАТ ЛАУМУЛИН ЦЕНТРАЛЬНАЯ АЗИЯ В ЗАРУБЕЖНОЙ ПОЛИТОЛОГИИ И МИРОВОЙ ГЕОПОЛИТИКЕ Том V Центральная Азия в XXI столетии Алматы – 2009 УДК 327 ББК 66.4 (0) Л 28 Рекомендовано к печати Ученым Советом Казахстанского института стратегических исследований при Президенте Республики Казахстан Научное издание Рецензенты: Доктор исторических наук, профессор Байзакова К.И. Доктор политических наук, профессор Сыроежкин...»

«Министерство образования науки Российской Федерации Российский университет дружбы народов А. В. ГАГАРИН ПРИРОДООРИЕНТИРОВАННАЯ ДЕЯТЕЛЬНОСТЬ УЧАЩИХСЯ КАК ВЕДУЩЕЕ УСЛОВИЕ ФОРМИРОВАНИЯ ЭКОЛОГИЧЕСКОГО СОЗНАНИЯ Монография Издание второе, доработанное и дополненное Москва Издательство Российского университета дружбы народов 2005 Утверждено ББК 74.58 РИС Ученого совета Г 12 Российского университета дружбы народов Работа выполнена при финансовой поддержке РГНФ (проект № 05-06-06214а) Н а у ч н ы е р е...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РЕСПУБЛИКИ КАЗАХСТАН КОМИТЕТ НАУКИ ИНСТИТУТ ФИЛОСОФИИ И ПОЛИТОЛОГИИ КАЗАХСТАН В ГЛОБАЛЬНОМ МИРЕ: ВЫЗОВЫ И СОХРАНЕНИЕ ИДЕНТИЧНОСТИ Посвящается 20-летию независимости Республики Казахстан Алматы, 2011 1 УДК1/14(574) ББК 87.3 (5каз) К 14 К 14 Казахстан в глобальном мире: вызовы и сохранение идентичности. – Алматы: Институт философии и политологии КН МОН РК, 2011. – 422 с. ISBN – 978-601-7082-50-5 Коллективная монография обобщает результаты комплексного исследования...»

«Vinogradov_book.qxd 12.03.2008 22:02 Page 1 Одна из лучших книг по модернизации Китая в мировой синологии. Особенно привлекательно то обстоятельство, что автор рассматривает про цесс развития КНР в широком историческом и цивилизационном контексте В.Я. Портяков, доктор экономических наук, профессор, заместитель директора Института Дальнего Востока РАН Монография – первый опыт ответа на научный и интеллектуальный (а не политический) вызов краха коммунизма, чем принято считать пре кращение СССР...»

«ISSN 2075-6836 Фе дера льное гос уд арс твенное бюджетное у чреж дение науки ИнстИтут космИческИх ИсследованИй РоссИйской академИИ наук (ИкИ Ран) А. И. НАзАреНко МоделИровАНИе космического мусора серия механИка, упРавленИе И ИнфоРматИка Москва 2013 УДК 519.7 ISSN 2075-6839 Н19 Р е ц е н з е н т ы: д-р физ.-мат. наук, проф. механико-мат. ф-та МГУ имени М. В. Ломоносова А. Б. Киселев; д-р техн. наук, ведущий науч. сотр. Института астрономии РАН С. К. Татевян Назаренко А. И. Моделирование...»

«Министерство сельского хозяйства Российской Федерации Федеральное государственное образовательное учреждение высшего профессионального образования Мичуринский государственный аграрный университет А.Г. КУДРИН ФЕРМЕНТЫ КРОВИ И ПРОГНОЗИРОВАНИЕ ПРОДУКТИВНОСТИ МОЛОЧНОГО СКОТА Мичуринск - наукоград РФ 2006 PDF created with FinePrint pdfFactory Pro trial version www.pdffactory.com УДК 636.2. 082.24 : 591.111.05 Печатается по решению редакционно-издательского ББК 46.0–3:28.672 совета Мичуринского...»

«РОССИЙСКАЯ АКАДЕМИЯ НАУК ИНСТИТУТ ЛИНГВИСТИЧЕСКИХ ИССЛЕДОВАНИЙ Л. З. Сова АФРИКАНИСТИКА И ЭВОЛЮЦИОННАЯ ЛИНГВИСТИКА САНКТ-ПЕТЕРБУРГ 2008 Л. З. Сова. 1994 г. L. Z. Sova AFRICANISTICS AND EVOLUTIONAL LINGUISTICS ST.-PETERSBURG 2008 УДК ББК Л. З. Сова. Африканистика и эволюционная лингвистика // Отв. редактор В. А. Лившиц. СПб.: Издательство Политехнического университета, 2008. 397 с. ISBN В книге собраны опубликованные в разные годы статьи автора по африканскому языкознанию, которые являются...»

«Е.А. Урецкий Ресурсосберегающие технологии в водном хозяйстве промышленных предприятий 1 г. Брест ББК 38.761.2 В 62 УДК.628.3(075.5). Р е ц е н з е н т ы:. Директор ЦИИКИВР д.т.н. М.Ю. Калинин., Директор РУП Брестский центр научно-технической информации и инноваций Государственного комитета по науке и технологиям РБ Мартынюк В.Н Под редакцией Зам. директора по научной работе Полесского аграрно-экологического института НАН Беларуси д.г.н. Волчека А.А Ресурсосберегающие технологии в водном...»

«УА0600900 А. А. Ключников, Э. М. Ю. М. Шигера, В. Ю. Шигера РАДИОАКТИВНЫЕ ОТХОДЫ АЭС И МЕТОДЫ ОБРАЩЕНИЯ С НИМИ Чернобыль 2005 А. А. Ключников, Э. М. Пазухин, Ю. М. Шигера, В. Ю. Шигера РАДИОАКТИВНЫЕ ОТХОДЫ АЭС И МЕТОДЫ ОБРАЩЕНИЯ С НИМИ Монография Под редакцией Ю. М. Шигеры Чернобыль ИПБ АЭС НАН Украины 2005 УДК 621.039.7 ББК31.4 Р15 Радиоактивные отходы АЭС и методы обращения с ними / Ключников А.А., Пазухин Э. М., Шигера Ю. М., Шигера В. Ю. - К.: Институт проблем безопасности АЭС НАН Украины,...»

«Камчатский государственный технический университет Профессорский клуб ЮНЕСКО (г. Владивосток) Е.К. Борисов, С.Г. Алимов, А.Г. Усов Л.Г. Лысак, Т.В. Крылова, Е.А. Степанова ЭКСПЕРИМЕНТАЛЬНАЯ ДИНАМИКА СООРУЖЕНИЙ. МОНИТОРИНГ ТРАНСПОРТНОЙ ВИБРАЦИИ Петропавловск-Камчатский 2007 УДК 624.131.551.4+699.841:519.246 ББК 38.58+38.112 Б82 Рецензенты: И.Б. Друзь, доктор технических наук, профессор Н.В. Земляная, доктор технических наук, профессор В.В. Юдин, доктор физико-математических наук, профессор,...»

«Институт биологии моря ДВО РАН В.В. Исаева, Ю.А. Каретин, А.В. Чернышев, Д.Ю. Шкуратов ФРАКТАЛЫ И ХАОС В БИОЛОГИЧЕСКОМ МОРФОГЕНЕЗЕ Владивосток 2004 2 ББК Монография состоит из двух частей, первая представляет собой адаптированное для биологов и иллюстрированное изложение основных идей нелинейной науки (нередко называемой синергетикой), включающее фрактальную геометрию, теории детерминированного (динамического) хаоса, бифуркаций и катастроф, а также теорию самоорганизации. Во второй части эти...»

«МИНИСТЕРСТВО ОБРАЗОВАНИЯ И НАУКИ РОССИЙСКОЙ ФЕДЕРАЦИИ ФЕДЕРАЛЬНОЕ ГОСУДАРСТВЕННОЕ БЮДЖЕТНОЕ ОБРАЗОВАТЕЛЬНОЕ УЧРЕЖДЕНИЕ ВЫСШЕГО ПРОФЕССИОНАЛЬНОГО ОБРАЗОВАНИЯ САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ ЭКОНОМИКИ И ФИНАНСОВ КАФЕДРА ЦЕНООБРАЗОВАНИЯ И ОЦЕНОЧНОЙ ДЕЯТЕЛЬНОСТИ Т.Г. КАСЬЯНЕНКО СОВРЕМЕННЫЕ ПРОБЛЕМЫ ТЕОРИИ ОЦЕНКИ БИЗНЕСА ИЗДАТЕЛЬСТВО САНКТ-ПЕТЕРБУРГСКОГО ГОСУДАРСТВЕННОГО УНИВЕРСИТЕТА ЭКОНОМИКИ И ФИНАНСОВ ББК 65. К Касьяненко Т.Г. К 28 Современные проблемы теории оценки бизнеса / Т.Г....»






 
© 2013 www.diss.seluk.ru - «Бесплатная электронная библиотека - Авторефераты, Диссертации, Монографии, Методички, учебные программы»

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