Среди комментариев к заметке про столкновительные траектории Земли с другими планетами было два вопроса, ответ на которые я решил вынести в отдельный пост.1.Вы же тут пишете что какой бы малой ни была бы ошибка на начальный момент, уже через несколько сотен миллионов лет будет совсем хаос. Но ошибки порядка машинного эпсилон не устранимы в принципе ... Т.е. поведение на дольших масштабах времени обусловлено не физикой, а ошибками округлений на начальный момент вычислений. Тогда не ясно какая научная ценность в таких рассчетах, если нет возможности оценить влияние ошибки вычислений: уже за пол миллиарда лет ошибки набегут настолько, что вся дальнейшая траектория станет одним сплошным "артефактом". Хотя я в комментариях уже пояснял это, но наверно полезно еще раз подчеркнуть. Вначале терминологическое пояснение. Под словом "траектория" тут понимается не траектория одной планеты, а как бы совокупная траектория всех планет. Если следить только за 8 планетами, то в каждый момент времени состояние Солнечной системы будет описываться 24 числами -- пространственными координатами этих планет, т.е. некоторой точкой в 24-мерном пространстве состояний Солнечной системы. Под "траекторией" понимается траектория этой точки в 24-мерном пространстве.Теперь ответ на вопрос. Да, из-за экспоненциального разбегания изначально близких траекторий с ляпуновским показателем около 5 млн. лет нет никаких шансов рассчитать истинную эволюцию Солнечной системы на миллиардо-летнем масштабе. Траектория начинает хаотично диффундировать в доступной ей области многомерного пространства, и в какую сторону она удиффундирует -- неизвестно. Что можно сделать в этой ситуации, так это примерно попытаться очертить эту доступную для хаотического диффундирования область. Ведь наличие хаоса еще не означает, что планеты смогут двигаться как угодно -- в конце концов, законы сохранения энергии и углового момента никто не отменял. Поэтому доступная для хаотического диффундирования область будет низкоразмерной и скорее всего какой-то хитрой формы. Эту область желательно очертить, выяснить, какие траектории для нее типичны, а какие -- нет. Выяснить, в какую сторону диффундирование идет легче, а в какую -- труднее. Именно с этой целью и запускается сразу тысячи сеансов моделирования, которые вначале отличаются на совершенно мизерную величину -- они как бы "зондируют" доступную для диффундирования область. И эта область -- это уже довольно достоверная, не вероятностая характеристика Солнечной системы, а не артефакт численных расчетов. Её форма уже должна быть нечувствительна к погрешностям округления, и такие вещи действительно проверяют.В последние годы эту область только-только начали ощупывать, поэтому никакой общей картинки привести пока не могут. Могут лишь "рассказать словами" о том или ином уголке этой области. Заголовок оригинальной статьи в Nature Physics, если его перевести на этот язык, будет звучать примерно так: "Доступная для 5-млрд.-летнего диффундирования область включает в себя так же и далекий уголок со столкновениями".2....показал ссылку знакомому, вот его комментарий. А что тут комментировать, резонанса между Меркурием и Юпитером не может быть по определению, Солнце рядом... Что ж, пришлось копнуть немножко вглубь.Речь тут идёт не об орбитальных резонансах между Юпитером и Меркурием, а о резонансе между частотами прецессии перигелия орбит этих двух планет.Цитирую по статье того же Жака Ласкара, "Chaotic diffusion in the Solar System", Icarus 185, 312–330 (2008); статья свободно доступна как епринт arXiv:0802.3371.В нынешнюю эпоху перигелий Меркурия прецессирует со скоростью 5,59" за год. Из них 0,43"/год -- вклад ОТО. Перигелий Юпитера прецессирует со скоростью 4,257"/год, вклад ОТО пренебрежимо мал. Т.е. обе орбиты вращаются, совершая оборот примерно за 200 тыс.лет.Скорости прецессии этих двух планет довольно близки уже сейчас, и поскольку с течением времени параметры всех орбит немного меняются, вполне может наступить ситуация, когда довольно долго эти остаются почти одинаковыми. Это будет означать, что в течение длительного времени (очень много периодов вращения) эллипс юпитерианской орбиты будет находиться под почти постоянным углом к эллипсу меркурианской орбиты, и значит слабое влияние Юпитера будет накапливаться, а не усредняться. Так и получается раскачка эксцентриситета Меркурия.Вот одно конкретное моделирование (т.е. одна конкретная траектория) из этой же статьи:
Здесь показан угол между направлениями перигелия Меркурия и Юпитера (угол показан по модулю 6π, а не 2π, просто для визуального удобства). Моделирование делалось в прошлое. Видно, что в период времени между -1284 и -1286 млн. лет (т.е. в течение двух миллионов лет!) этот угол оставался примерно постоянным. А вот что в это время происходило с эксцентриситетом Меркурия:
Видно, что именно в течение этих двух млн. лет его эксцентриситет раскачался с 0,4 до 0,8.Кстати, из приведенных выше чисел по скорости прецессии видна роль ОТО. Если не учитывать поправки ОТО, то частоты прецессии становятся еще ближе, т.е. резонанс наступает гораздо охотнее. Это, по-видимому, и объясняет, почему без ОТО примерно половина траекторий вырабатывает очень большой эксцентриситет Меркурия, а с поправками -- всего около 1%.
Комментарии (0)