Александр Филиппов - Многоликий солитон

Трудно сказать, интересны ли были бы такие опыты Фарадею или Эрстеду. Скорее всего, они показались бы им чересчур абстрактными, лишенными жизни. Но вот Эйлер или Гayсс наверняка пришли бы в восторг. Они сами были феноменальными вычислителями, а идея численного эксперимента была им близка и понятна. В сущности, Эйлер и ввел понятие о численном эксперименте в математике: «Покажется немало парадоксальным приписывать большое значение наблюдениям даже в той части математических наук, которую обычно называют чистой математикой, ибо существует распространенное убеждение, что наблюдения ограничиваются физическими объектами, которые воздействуют на наши чувства. Поскольку мы уверены, что числа принадлежат одному лишь чистому разуму, нам очень трудно представить себе пользу наблюдений и квазиэкспериментов при изучении природы чисел. В действительности же... известные сегодня свойства чисел были, по большей части, открыты путем наблюдения...»

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

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

Мы, люди XX в., не любим беспорядка, и слово «хаос» для нас почти ругательное. Примерно такой смысл вкладывает в него Борис Пастернак. Поэт же прошлого века, Федор Иванович Тютчев, вкладывает в это слово совсем иной смысл. Чуть раньше он называет его «древним» и «родимым». Он, конечно, имеет в виду хаос древних греков — некое первичное состояние мира, все порождающее и все поглощающее, или «первоматерию». На современном языке такой хаос больше всего похож на состояние нашей Вселенной в первые мгновения после «Большого взрыва». Одна из глубоких проблем современной космологии — понять, как из этого хаотического состояния образовалась современная упорядоченная Вселенная с ее галактиками и звездами. На некоторых этапах развития вселенной солитоны, видимо, помогали ей упорядочиваться. Мы уже вскользь упоминали об этом, когда познакомились со спиральными рукавами галактик.

Другой тип беспорядочного движения хорошо известен каждому из наблюдений над течением воды. Проще всего увидеть рождение хаотического движения воды, если двигать какой-нибудь предмет плохо обтекаемой формы. Уже при небольшой скорости движения возникают вихри. При очень большой скорости может возникнуть так называемый турбулентный след, подобный тому, который наблюдается за кормой быстро движущегося корабля. В области следа частицы воды движутся совершенно беспорядочно, хаотически. Такие движения жидкости впервые начали изучать Кельвин, Буссинеск, Рейнольдс и Рэлей. Термин «турбулентность» ввел в обиход Кельвин, произведя его от латинского «turbulentus» (беспокойный, беспорядочный). Первые опыты по изучению турбулентного движения воды в обыкновенных водопроводных трубах выполнил Рейнольдс в 1883 г.

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

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

Представьте себе ученую «блоху», владеющую всеми тремя действиями арифметики и прыгающую не просто так, а по определенному закону. Если она в момент времени tn = n · Δt (n = 0, 1, 2, 3, ...) сидит в точке xn на оси x, то в следующий момент tn+1 = (n + 1) Δt она перепрыгивает в точку xn+1 = b - xn2, где b — некоторое выбранное число, свое для каждой блохи (назовем ее, скажем, «постоянной блохи»). Пусть блоха начинает движение из некоторой точки отрезка -2  x  +2. Наша задача — определить, куда она может убежать за большое время, т. е. представить себе, каким может быть xn при больших значениях n.

Как ни проста эта задача на вид, вам едва ли удастся найти ее решение без помощи микрокалькулятора. Однако прежде чем приступить к экспериментам, стоит немного подумать, чего от них можно ожидать. Рассмотрим кривую AA0A1 (рис. 7.1), соответствующую уравнению x2 + x = b. Если xn стремится к некоторому пределу, то предельное значение будет лежать на этой кривой.

Численные эксперименты, однако, показывают, что только блохи, для которых точки с координатами (x0, b) лежат внутри фигуры A0B1B1'A0' (кривая A0'A0 получается зеркальным отражением относительно оси Ob), приближаются к кривой A0А1. При этом на ветвь A0A они никогда не попадают, а если их постоянные b и начальные координаты x0 таковы, что точка (x0, b) лежит в заштрихованной области, то такие блохи убегают на бесконечность. Судьба наших ученых блох с постоянной b 0,75 полностью определена. Что бы они ни делали, они либо погибают «в бесконечности», либо притягиваются к точкам на кривой A0А1, где мы их легко ловим. Почти столь же печальна судьба блох, у которых (x0, b) лежит в области B1B2B2'B1'. Они в конце концов попадают на кривую A2A1A2', и при достаточно большом n перескакивают на каждом шаге с A1A2' на A1A2 и обратно. Уравнение этой кривой тоже легко определить, воспользовавшись тем, что в пределе больших n после двух последовательных скачков, блоха оказывается вблизи той же ветви, т. е. x2+n  xn. Отсюда находим xn+2 = b - (b - xn2)2 и в пределе, когда xn+2  xn  x, находим для x уравнение (b - x2)2 - (b - x) = 0. Легко проверить, что левая часть этого уравнения равна произведению двух множителей: (x2 + x - b) и (x2 - x + 1 - b). Обращение в нуль первого множителя дает кривую AA0A1, а второго — кривую A2A1A2'.

До сих пор мы могли бы более или менее точно предсказать, что будет происходить. При увеличении «блошиной» постоянной все, однако, быстро усложняется. Кривые, которые их притягивают, продолжают раздваиваться, и при b 1,5 скачки становятся почти непредсказуемыми, беспорядочными. На рисунке этому соответствует зачерненная область. Например, если блоха начинает движение на отрезке B3B3', то она притянется к отрезку C3C3'. Предельные значения ее координат плотно заполняют этот отрезок (на самом деле, как показывают дальнейшие эксперименты, устройство зачерненного притягивающего множества гораздо сложнее, читатель может попробовать изучить его в экспериментах).