Стартовая Предметный указатель Новости науки и техники
Новости науки и техники
ПРОГНОЗ СОЛНЕЧНОЙ НЕПОГОДЫ
В будущем исследователи будут следить за рентгеновскими лучами от Юпитера, чтобы выяснить, что происходит на дальней стороне Солнца, невидимой с Земли, сообщает New Scientist. Далее...

Солнечная активность

устойчивость упругих систем

УСТОЙЧИВОСТЬ УПРУГИХ СИСТЕМ -свойство упругих систем возвращаться к состоянию равновесия после малых отклонений их из этого состояния. Понятие У. у. с. тесно связано с общими понятиями устойчивости движения и равновесия. Устойчивость является необходимым условием для любой конструкции. Потеря устойчивости

5050-104.jpg

где 5050-105.jpg-матрица 3-поворотов, 5050-106.jpg , 5050-107.jpg Подчеркнём, что частота w в множестве (9) не фиксирована.

Если возмущённый солитон описывать полем

5050-108.jpg то возмущение x определим как 5050-109.jpg Метрики r0, r выберем в виде

5050-110.jpg

где || ||-норма в 5050-111.jpg , значок С обозначает совместную норму в 5050-112.jpg

5050-113.jpg

Изучим Q-yстойчивость солитонных решений (8), наложив условие фиксации заряда, уже предполагающееся в определении (10):

5050-114.jpg

Введем удобные для дальнейшего обозначения:

5050-115.jpg

5050-116.jpg

Выберем в качестве функционала Ляпунова интеграл движения

5050-117.jpg

где 5050-118.jpg-энергия поля. Его вторая вариация может быть представлена в виде

5050-119.jpg

где введены самосопряжённые операторы

5050-120.jpg

Из (12) следует, что для положительной определённости d2V необходимо выполнение неравенств Fp > 0, h > 0. Оказывается, что безузловые солитоны (u>0) могут быть устойчивыми, тогда как узловые солитоны (для к-рых на нек-рых поверхностях u = 0) всегда неустойчивы. Заметим, что для безузловых солитонов спектр оператора 5050-121.jpg неотрицателен, т. к. 5050-122.jpg и>0, и поэтому и - первая собственная ф-ция оператора L^2, тогда как нулевая мода x2=u исключается выбором метрики r. Анализируя структуру второй вариации (12), можно установить справедливость следующей теоремы (Q-теоре-мы): безузловые стационарные решения (8) Q-устойчивы по Ляпунову в области

5050-123.jpg

если в ней оператор 5050-124.jpg имеет единств. отрицат. собств. значение, а собств. ф-ция y- удовлетворяет условию

5050-125.jpg

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

5050-126.jpg

где 5050-127.jpg

Вычисляя его производную W., находим:

5050-128.jpg

Отсюда следует, что в области 5050-129.jpg т. е. имеет место неустойчивость солитонов.

Чтобы убедиться в неустойчивости узловых солитонов, заметим, что в этом случае возмущение x2 всегда содержит решение однородного ур-ния 5050-130.jpg допускающего знакопеременный интеграл "энергии"

5050-131.jpg

т. к. оператор 5050-132.jpg имеет отрицат. собств. значения. Это видно из ур-ния 5050-133.jpg и наличия узлов у ф-ции u(r). Неустойчивость доказывается существованием функционала Четаева 5050-134.jpg для к-рого W>0 в области 5050-135.jpg

Рассмотрим примеры применения Q-теоремы для анализа устойчивости солитонов в D-мерном пространстве.

1) С т е п е н н а я м о д е л ь. В этом случае5050-136.jpg и ф-ция и(х)удовлетворяет ур-нию

5050-137.jpg

к-рое имеет безузловое решение и(r)при условиях |w|<1 0<1-1/n<=2/D. Выполнив в (15) замену переменных: x=r(1-w2)-1/2, u=u(1-w2)s, s-1=2(n-1), находим заряд Q(w) невозмущённого солитона:

5050-138.jpg

Из (16) следует, что условие (13) выполнено для частот

5050-139.jpg

Условие (14) также выполнено, т. к. 5050-140.jpg а ф-ция 5050-141.jpg как первая собств. ф-ция оператора 5050-142.jpg. Поэтому неравенство (17) определяет область устойчивости безузловых солитонов.

2) Л о г а р и ф м и ч е с к а я м о д е л ь задаётся ф-цией F=p+s(1-lns) и допускает решения вида

5050-143.jpg

Отсюда находим зависимость заряда от частоты:

5050-144.jpg

определяющую, согласно (13), область устойчивости:

5050-145.jpg

3) Шрёдингера уравнение нелинейное5050-146.jpg5050-147.jpg , допускает решения (8) с амплитудой u, подчиняющейся ур-нию (15) с переобозначением 5050-148.jpg . Замена переменных x=r|w|-1/2, u=u|w|s, s-1=2(n-1) позволяет найти заряд как ф-цию от w:

5050-149.jpg

Отсюда следует, что в области устойчивости 1 < п < 1+2/D, а при п> 1 +2/D солитоны неустойчивы. Это устанавливается с помощью функционала Четаева5050-150.jpg

3. Метод Захарова - Кузнецова (1974). Метод состоит в доказательстве ограниченности снизу энергии консервативной системы при условии фиксации нек-рых дополнит. интегралов движения. Проиллюстрируем метод на последнем примере, показав, что интеграл энергии 5050-151.jpgв5050-152.jpg оценивается снизу через заряд Q. В самом деле,

5050-153.jpg

Вводя обозначение 5050-154.jpg ..., и используя неравенства 5050-155.jpg приходим к оценке

5050-156.jpg

Если 5>3n, то правая часть этого неравенства имеет минимум при

5050-157.jpg

Поэтому энергия5050-158.jpg при фиксированном I2 = Q также имеет минимум, к-рый и реализуется на нек-рой стабильной конфигурации.

Используем метод Захарова - Кузнецова для доказательства существования стабильных солитонов ещё в двух распространённых моделях.

1) Кортевега - де Фриса уравнение (D = 1)5050-159.jpg5050-160.jpg описывает волны на мелкой воде и допускает законы сохранения энергии

5050-161.jpg

и импульса5050-162.jpg Используя неравенство Гальяр-до - Ниренберга-Ладыженской 5050-163.jpgполучаем оценку для энергии снизу:

5050-164.jpg

Минимизируя правую часть этого неравенства по ||дxj ||, находим5050-165.jpg 5050-166.jpg Т. о., при фиксированном импульсе Р=I2 энергия ограничена снизу и имеет минимум, к-рый реализуется на нек-рой устойчивой конфигурации.

2) Кадомцева - Петвиашвили уравнение (D = 2)

5050-167.jpg

рассматривается как двумерное обобщение ур-ния Кортевега- де Фриса и также допускает законы сохранения энергии

5050-168.jpg

и импульса5050-169.jpg Воспользуемся неравенством Гёльдера5050-170.jpg, а также очевидными неравенствами

5050-171.jpg

объединяя к-pые, приходим к соотношению 5050-172.jpg , позволяющему получить оценку для энергии снизу:

5050-173.jpg

Минимизируя правую часть в (18) по ||дxj|| и ||дyw||, получаем неравенство5050-174.jpg означающее, что при фиксированном импульсе Р=I2 минимум энергии реализуется на нек-рой стабильной солитонной конфигурации.

4. Пример применения прямого метода в кинетической теории плазменных солитонов. Рассмотрим эл--статич. приближение Власова - Пуассона в одномерном случае (D = 1). Ур-ния для ф-ции распределения электронов f(t, x, u) и напряжённости электрич. поля в плазме E(t, x)в приближении тяжёлых ионов имеют вид (распределение ионов не зависит от времени)

5050-175.jpg

С учётом граничных условий

5050-176.jpg

в системе отсчёта, связанной с центром распределения 5050-177.jpg электрич. поле исключается:

5050-178.jpg

Пусть невозмущённое решение ур-ний (19) стационарно:

5050-179.jpg

где 5050-180.jpg -энергия электрона, m = sign u. Т. к. f>0 , полагаем 5050-181.jpg считая c0 решением ур-ния

5050-182.jpg

где 5050-183.jpg При этом возмущение x=c-c0 с учётом (20) и линеаризованного условия нормировки 5050-184.jpg удобно представить в виде

5050-185.jpg

считая, что j удовлетворяет линеаризованному ур-нию

5050-186.jpg

где введены операторы

5050-187.jpg

Из ур-ния (21) следует, что существует интеграл движения

5050-188.jpg

В случае e= -1 функционал (22) положительно определён, что говорит об устойчивости монотонных по энергии w распределений - теорема Ньюкомба - Гарднера (клас-сич. пример: распределение Максвелла - Больцмана f0=Ae-w). Покажем, что монотонные распределения гло-бально устойчивы, выбрав функционал Ляпунова

5050-189.jpg

где l-множитель Лагранжа, G(f) - нек-рая вспомогательная ф-ция, определяемая из условия стационарности V1. Из условия dV1(f0)=0 находим5050-190.jpg или, после дифференцирования по w, 5050-191.jpg. Т. о., V1 - глобально выпуклый функционал. В частности, полагая 5050-192.jpg убеждаемся, что5050-193.jpg Однако если распределение f0 немонотонно по энергии, то функционал (22) знакопеременный, что говорит о неустойчивости. В самом деле, для функционала Четаева

5050-194.jpg

где F-решение вспомогат. ур-ния

5050-195.jpg

найдём, что 5050-196.jpg в области V<0. Т. о., немонотонные распределения неустойчивы по метрикам r0, r где 5050-197.jpg

(Подробное изложение теории прямого метода Ляпуно-ва и его приложений смотри в прилагаемом списке литературы.)

Лит.: Ляпунов А. М., Общая задача об устойчивости движения, 2 изд., Л.- М., 1935; Зубов В. И., Методы А. М. Ляпунова и их применение, Л., 1957; Мовчан А. А., Устойчивость процессов по двум метрикам, "Прикл. матем. и мех.", 1960, т. 24, в. 6, с. 988; Жидков Е. П., Кирчев И. П., Устойчивость решений вида уединенных волн некоторых нелинейных уравнений математической физики, "ЭЧАЯ", 1985, т. 16, в. 3, с. 597; Рыбаков Ю. П., Устой-чивость многомерных солитонов в киральных моделях и гравитации, в кн.: Итоги науки и техники, сер. Классическая теория поля и теория гравитации, т. 2, М., 1991, с. 56; Benjamin Т. В., Stability of solitary waves, "Proc. Roy. Soc.", 1972, v. 328A, p. 153; Makhan-kov V. G., Dynamics of classical solitons (in non-integrable systems), "Phys. Repts", 1978, v. 35, № 1, p. 1; Holm D. D. [a.o.], Nonlinear stability of fluid and plasma equilibrium, "Phys. Repts", 1985, v. 123, № 1 -2, p. l;ShatahJ., Strauss W., Instability of nonlinear bound states, "Comm. Math. Phys.", 1985, v. 100, № 2, p. 173; Kuzne- tsov E. A., Rubenchik A. M., Zakharov V. E., Soliton stability in plasmas and hydrodynamics, "Phys. Repts", 1986, v. 142, № 3, p. 103; Grillakis M., Shatah J., Strauss W., Stability theory of solitary waves in the presence of symmetry. I, II, "J. Funct. Anal.", 1987, v. 74, № 1, p. 160; 1990, v. 94, № 2, p. 308. Ю. П. Рыбаков.

  Предметный указатель