Моделирование образования планет на ЭВМ
Еще французский математик Лаплас и немецкий философ Иммануил Кант предполагали, что Солнце и планеты образовались из вращающейся туманности. Сегодня подобный процесс можно попытаться смоделировать на компьютере. В дальнейшем я буду основываться на результатах расчетов, которые осуществили калифорнийский астрофизик Питер Боденгеймер и Вернер Чарнутер частью по отдельности, частью совместно в Мюнхене. Поначалу они намеревались объяснить происхождение Солнца и планет. Но дело обернулось совсем иначе.
Насколько просто моделировать на компьютере процессы, обладающие сферической симметрией, понимаешь лишь тогда, когда берешься за следующую по сложности задачу. В сферически-симметричной задаче в любой момент времени все параметры зависят только от расстояния до центра. Когда, например, в модели Ларсона вещество нагревается, то одновременно нагревается весь слой, расположенный на одном и том же расстоянии от центра, т. е. слой, лежащий на поверхности сферы определенного радиуса. Если вещество не вращается, сферическая симметрия является хорошим приближением; все частицы, участвующие в коллапсе, имеют одинаковую судьбу независимо от того, по какому направлению они движутся.
Вращение нарушает сферическую симметрию. Силы, действующие на частицы, движущиеся со стороны полюса, отличаются от сил, действующих на частицы, приходящие со стороны экватора. Сферической симметрии больше нет. Но это не значит, что задача приобретает невероятную сложность. Она сохраняет определенную степень симметрии. Например, в экваториальной плоскости частицы движутся к центру по различным направлениям, но все эти направления равнозначны. В таком случае говорят, что процесс обладает осевой (аксиальной) симметрией. Осесимметричные процессы рассчитывать на ЭВМ уже гораздо труднее, но и к ним можно найти подход. Боденгеймер и Чарнутер построили компьютерную модель коллапсирующего вращающегося облака (рис. 13.2). Вначале все идет по Ларсону: облако сжимается, и в центре образуется уплотнение. Чем сильнее сжимается облако, тем больше дает о себе знать центробежная сила: облако сплющивается. В конце концов образуется плоский диск. Теперь в коллапсе участвует лишь вещество, находящееся вблизи оси вращения, в экваториальной же плоскости газ движется к центру медленно и в какой-то момент прекращает движение. Вместо ядра, на которое со всех сторон падает вещество, мы имеем диск, на который вещество падает лишь по оси. Диск, экваториальный радиус которого в восемь раз больше его толщины, занимает огромное пространство с поперечником около 120 радиусов орбиты Плутона, самой далекой планеты Солнечной системы. Один оборот вокруг центра совершается за 300000 лет.