Взаимодействие космических тел с атмосферой и поверхностью Земли

Светцов, Владимир Владимирович

Бесплатный автореферат и диссертация по наукам о земле на тему
Взаимодействие космических тел с атмосферой и поверхностью Земли
ВАК РФ 25.00.29, Физика атмосферы и гидросферы

Автореферат диссертации по теме "Взаимодействие космических тел с атмосферой и поверхностью Земли"

Институт динамики геосфер Российской академии наук

ии3445459

На правах рукописи

УДК 550 3+518 5

Светцов Владимир Владимирович

ВЗАИМОДЕЙСТВИЕ КОСМИЧЕСКИХ ТЕЛ С АТМОСФЕРОЙ И ПОВЕРХНОСТЬЮ ЗЕМЛИ

Специальность 25 00 29 - Физика атмосферы и гидросферы

АВТОРЕФЕРАТ

диссертации на соискание ученой степени доктора физико-математических наук

РГБ ОД

2 В АВГ 2008

Москва 2008

003445459

Работа выполнена в Институте динамики геосфер РАН

Официальные оппоненты:

Доктор физико-математических наук

Витязев А В

Член-корр РАН,

доктор физико-математических наук

Лобковский JIИ

Доктор физико-математических наук

Чарахчьян А А

Ведущая организация: Институт математического моделирования РАН

Защита диссертации состоится 16 октября 2008 г в 11 часов на заседании диссертационного совета Д002 050 01 при Институте динамики геосфер Российской академии наук (ИДГ РАН) по адресу 119334 Москва, Ленинский проспект 38, корпус 1

С диссертацией можно ознакомиться в библиотеке ИДГ РАН Автореферат разослан и^Л- 2008 г

Ученый секретарь диссертационного совета кандидат физико-

математических наук

В А Рыбаков

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность проблемы. Падения космических тел на Землю представляют большой интерес в геологии, геофизике, планетологии и астрономии, а также и с общей естественнонаучной точки зрения Земля, как и другие планеты, рождалась в процессе многочисленных ударов планетезималей о ее поверхность Удары крупных тел влияли на образование и эволюцию геосфер, в том числе атмосферы Кометы и астероиды, занося органическое вещество на поверхность Земли, могли способствовать зарождению жизни С другой стороны, удары крупных тел приводят к гибели живых организмов Астероиды и кометы размером порядка десяти километров могли приводить к глобальным катастрофам и массовым вымираниям видов в фанерозое Угрожающая человечеству астероидно-кометная опасность все больше осознается с углублением наших знаний В США развернута программа астрономических наблюдений за потенциально опасными космическими объектами, которые могут приблизиться к Земле Наблюдения и каталогизация этих тел ведутся и в России Реальную опасность для людей представляют объекты, начиная с диаметров 50-100 м (предполагаемый размер тела, вызвавшего Тунгусскую катастрофу 1908 г) Поэтому была поставлена задача выявления космических тел размером до ~100 м, число которых составляет многие десятки тысяч Один из известных объектов - астероид Апофис диаметром около 350 м, который сблизится с Землей в 2029 г

Проблема астероидно-кометной опасности выдвинула много проблем В их числе изучение процессов, происходящих при столкновении космических тел с планетой, механизмов их воздействия на окружающую среду Падающие космические тела взаимодействуют прежде всего с атмосферой, причем в большинстве случаев метеороиды не долетают до поверхности Изучение взаимодействия космических тел с атмосферой актуально не только с точки зрения определения состояния атмосферы после удара, но и с точки зрения определения параметров падающих на Землю тел С целью выяснения их свойств проводятся дорогостоящие космические эксперименты, такие как Stardust (захват частиц кометы Вильде 2), Deep Impact (удар по комете Темпеля 1) или планируемый Европейским космическим агентством эксперимент Дон Кихот, в котором предполагается осуществить удар по астероиду Но и поток небольших метеороидов, часто входящих в атмосферу Земли, несет в себе огромную информацию о составе и природе малых, а также и больших космических тел, которая требует расшифровки

Научная новизна. Удары космических тел по поверхности планет изучаются и моделируются уже довольно давно Вычислители, обратившиеся к

задачам об ударах космических тел, и геологи, использующие имеющиеся компьютерные программы, проделали большую работу, которая позволила понять основные процессы при образовании ударных структур и геологические последствия ударов, хотя многие вопросы, связанные с ударным кратерообразованием и интерпретацией геологических данных, все еще остаются нерешенными Процессы же, связанные с движением космических тел в атмосферах планет, изучены в гораздо меньшей степени. Новым общим итогом диссертации стало обнаружение общих закономерностей при падении и взаимодействии с атмосферой тел разного масштаба Значительная часть конкретных результатов работы была получена впервые Новизна результатов, связанная как с методами расчетов, так и с постановкой задач, заключается в следующем (1) Построена лагранжева полностью консервативная численная схема, обладающая определенными преимуществами при решении газодинамических задач, связанных с ударами космических тел (2) Разработан новый экономичный метод численного решения одномерных уравнений переноса излучения для расчетов остывания областей нагретого газа (3) Получены и исследованы новые решения уравнений газовой динамики при обтекании тел гиперзвуковым потоком газа (4) Проведены исследования нескольких падений космических тел Тунгусского, Сихотэ-Алинского, болида Шумава, кометы Шумейкер-Леви 9, в которых выявлены не изученные ранее особенности этих событий (5) Впервые полноценно исследованы касательные удары космических тел (под очень острыми углами к поверхности планеты), показано, какие геологические следы они могли оставить и определена их возможная роль в снабжении Земли органическим веществом (6) Впервые исследованы сильные акустико-гравитационные волны после ударов крупных космических тел, которые могут производить существенное механическое воздействие в глобальном масштабе (7) Построена модель эволюции атмосферы в период завершения аккумуляции Земли (8) Проведено математическое моделирование ударов очень крупных космических тел, приводящих к образованию временных силикатных атмосфер на ранней стадии существования планеты.

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

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

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

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

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

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

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

2 Исследовать атмосферные эффекты, вызываемые ударами космических тел при выделении их энергии в воздухе, определить параметры ударных волн, сильных акустико-гравитационных волн, потоки излучения Существенные атмосферные эффекты возникают при падении тел размером более 50-100 м Тела размером 50 - 100 м могут создать кратер, подобный кратеру Бэрринжера в Аризоне, если достигают поверхности, или произвести мощный взрыв в воздухе, как в Тунгусском явлении 1908 г Они представляют прямую угрозу для людей при падении в населенных районах

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

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

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

Основные научные положения, защищаемые в диссертации

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

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

3 Согласно проведенным расчетам в типичных случаях падений каменных тел размером 10-100 м на Землю эти тела дробятся на большое число мелких фрагментов, а потоки излучения могут приводить к полному испарению осколков в атмосфере Излучение из огненного шара, создаваемого в воздухе космическим телом, кометой или астероидом, размером 100 м и более метров, может расплавить грунт с образованием стекол и инициировать пожары в протяженных областях поверхности Земли

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

5 Моделирование ударов тел размером 100-1000 км по ранней Земле продемонстрировало, что в результате этих ударов образовывались силикатные атмосферы с коротким временем жизни Показано, что удары космических тел на завершающей стадии аккумуляции Земли играли существенную роль в образовании и эрозии ранних газовых атмосфер, которые были более массивны, чем современная атмосфера, с давлениями у поверхности 10-100 бар

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

неразрушенной кометной органики ~4109 лет назад был сравнимым с эндогенными источниками и потоком, приносимым межпланетной пылью

Апробация и внедрение результатов. Результаты, полученные в ходе выполнения диссертационной работы, вошли в отчеты Института динамики геосфер по программам Президиума РАН №13 и №16 "Изменения окружающей среды и климата природные катастрофы" и №25 "Проблемы зарождения биосферы Земли и ее эволюции", в отчет по астероидной опасности для МЧС по теме "Анализ и оценка последствий падения на Землю опасных космических объектов и разработка предложений по их смягчению" и монографию "Катастрофические воздействия космических тел", написанную коллективом авторов ИДГ РАН (Вариант этой монографии на английском языке "Catastrophic events caused by cosmic objects" был издан в 2008 г в издательстве Springer) Ряд работ вошел в отчеты, выполненные по грантам РФФИ, контрактам с Национальной лабораторией Сандиа (США) и по гранту Американского фонда исследований и развития (CRDF)

Результаты исследований докладывались на международных конференциях по астероидной опасности (Tucson, AZ, USA, 1993, Livermore, CA, USA, 1995), на лунно-планетных конференциях (Lunar Planetary Science Conference, Houston, TX, USA, 1994, 1996, 2002, 2005, 2006, 2007), на европейском совещании по столкновению кометы Шумейкер-Леви 9 с Юпитером (Garching, Germany, 1995), на международной конференции по космосу "Space-96" (Albuquerque, NM, USA, 1996), на международных совещаниях по ударному кратерообразованию, эволюции Земли и глобальным катастрофам (Ancona, Italy, 1995, Postojna, Slovenia, 1996, Quillan, France, 1999, Vienna, Austna, 2000, Granada, Spain, 2001, Mora, Sweden, 2002, Nordhngen, Germany, 2003), на совещаниях Метеоритного общества (Annual Meetings of the Meteontical Society, Munster, Germany, 2003, Zunch, Switzerland, 2006), на конференциях, посвященных Тунгусскому явлению 1908 г (Bologna, Italy, 1996, Красноярск, 1998, Москва, 2003, 2008), на всероссийской конференции "Астероидно-кометная опасность 2005" (Санкт-Петербург, 2005), на международной конференции "Околоземная астрономия 2007" (Терскол 2007)

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

Структура и объем диссертации. Диссертация состоит из введения, семи глав и заключения Объем диссертации - 335 страниц, включая 92 иллюстрации Список литературы содержит 549 наименований

Благодарности. Ряд разделов настоящей работы выполнялся совместно с сотрудниками ИДГ РАН И В Немчиновым, В М Хазинсом, В В Шуваловым, автор благодарит их за полезное сотрудничество Автор благодарен ИВ Немчинову и В В Шувалову за внимание к работе, поддержку и полезные консультации, И Б Косареву за предоставление таблиц оптических и термодинамических свойств веществ, Н А Артемьевой, Д Т Вассону, А.ГТ Голубю, Б А Иванову, Г В Печерниковой, О П Поповой, Т И Орловой и А В Тетереву за полезные обсуждения Автор искренне признателен Л М Белецкой за помощь в оформлении результатов

СОДЕРЖАНИЕ РАБОТЫ

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

Глава 1. Развитие методов математического моделирования течений, возникающих при ударах космических тел

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

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

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

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

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

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

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

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

В расчетах обтекания тел использовались несколько известных методов численного решения газодинамических уравнений крупных частиц, Годунова, TVD, РРМ (метод Годунова с высоким порядком точности) П 1 3 посвящен деталям применения этих методов и введенным дополнениям Численные схемы проверяются на тесте Квирка на четно-нечетные возмущения сетки, определяющем карбункулярную неустойчивость

В п 14 описывается экономичный метод расчета переноса излучения, который был разработан автором для решения одномерных радиационно-газодинамических задач Уравнения осредняются по углу и получается система уравнений для двух направлений Способ осреднения позволяет преодолеть трудности, связанные с появлением членов, обратно пропорциональных радиусу в задачах со сферической и цилиндрической симметрией Это достигается путем специального осреднения с введением дополнительного уравнения для интенсивности излучения вдоль радиуса Неявные разностные уравнения устойчивы без ограничений на шаг по времени Зависимость потоков излучения от энергии квантов учитывается в многогрупповом приближении Разностная схема обеспечивает точное решение в предельных случаях объемного высвечивания и лучистой теплопроводности

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

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

В диссертации рассматривается задача, в которой на пути движущегося с высокой скоростью тела возникает ограниченная в пространстве область газа пониженной плотности, что вносит кратковременные возмущения в течение газа перед телом Когда тело пролетает эту область, вихревые движения полностью не затухают, и оказывается, что при малых у существуют два режима стандартное стационарное обтекание тела гиперзвуковым потоком и нестационарное с вихрями, размер которых порядка отхода ударной волны от тела, и искаженным фронтом ударной волны Задача исследовалась с помощью численных решений уравнений Эйлера и уравнений Навье-Стокса различными методами Получаемая картина течения для вихревого режима при постоянной величине у в отдельные моменты времени показана на рис 1 Численные нестационарные решения зависят от метода расчета и разностной сетки, причем

они не проявляют сходимости с измельчением сетки. В связи с этим в работе рассматриваются вопросы применимости численных методов.

Рис. !. Вихревой режим обтекания бруска в безразмерных координатах. Ширина тела - две безразмерных единицы. Скорость набегающего потока направлена слева направо. Черным цветом показаны положения ударных волн. Серым цветом показаны области дозвукового течения. Результаты расчетов для у=1.1 методом РРМ на сетке с шагом h= 0.01 показаны для двух моментов безразмерного времени 1=680 (а) и /=720 (б). Результаты расчетов для у=1.05 методом Годунова на сетке с шагом /г=0.005 показаны для ¿=765 (в) и ?=820 (г).

Рис. 2. Результаты расчета обтекания метеороида в виде бруска шириной 2 м, падающего вертикально со скоростью 30 км/с на высоте 80 км. (а) - без учета излучения, (б) - с учетом излучения, картина течения блика к стандартному режиму обтекания, (в) и (г) - с излучением, вихревой режим обтекания (1=25.9 мс (в) и 1=27.4 мс (г) от начала расчета). Серым цветом показаны области дозвукового течения, а наклонной штриховкой - обтекаемое тело.

Моделирование обтекания ступеньки с использованием уравнений Навье-Стокса показало, что при у=1.1 и числе Рейнольдса 7?е=500 (определяемом по

параметрам набегающего потока) после окончания действия искусственного возмущения потока вихревой режим не развивается, а течение приходит к обычному стационарному режиму При Ле=1000 вихревой режим возникает, но отход ударной волны колеблется около среднего значения с относительно небольшой амплитудой При йе=2000 вихревой режим существует бесконечно долго Уменьшение у приводит к усилению вихревого режима, он образуется при большей вязкости К такому же эффекту приводит и излучение, которое охлаждает газ за фронтом ударной волны, обеспечивая большее сжатие и меньшее расстояние между ударной волной и телом Это проиллюстрировано на рис 2, где приведены результаты моделирования с использованием табличных уравнения состояния и коэффициентов поглощения воздуха

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

О 100 200 300

Время

Рис 3 Изменение максимального давления (отнесенного к давлению торможения) на передней поверхности обтекаемой ступеньки, вычисленного на сетке с шагами h= 0 01 (а) и h=0 005 (б) Расчет проведен методом РРМ для у=1 1

Данные, получаемые при регистрации падений метеорных тел в атмосфере, традиционно интерпретируются наблюдателями таким образом, что тела, разрушающиеся и тормозящиеся на больших высотах, считаются малоплотным кометным материалом Среди всех падающих метеороидов размером до 1 м, для которых получена статистика, таких тел насчитывается около 40% При такой интерпретации и экстраполяции наблюдательных данных на тела большего размера оказывается, что 70% всех метеороидов размером порядка 10 м составляют кометы Между тем, оценки, основанные на вычислении количества комет и астероидов, сталкивающихся с Землей, показывают, что кометы составляют не более 20% крупных ударников Из существования вихревого режима следует, что определенная группа каменных тел должна разрушаться примерно на 20 км выше, чем остальные, и поэтому, вероятно, что именно из-за вихревых эффектов, а не из-за малой плотности и прочности вещества значительное количество метеороидов разрушается на больших высотах

Глава 3. Падения разрушенных метеороидов в атмосфере Земли и их особенности (тела размером 1-100 м)

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

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

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

Рис.4. Вертикальное падение ледяного метеороида радиусом 100 м в атмосфере. Предполагалось, что в начальный момент времени сферическое тело находится на высоте 50 км, обладая скоростью 50 км/с. Частицы, расположеные внутри тела, не показаны. В моменты 0.5 и 0.6 с показаны положения частиц воздуха. В моменты времени £>0.7 с кружки представляют собой частицы метеороида, которые отделились от основного тела.

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

Исследуются несколько конкретных случаев падений тел - Сихотэ-Алинского метеоритного дождя, болида Шумава, падения фрагментов кометы Шумейкер-Леви 9 в атмосфере Юпитера и Тунгусского явления 1908 г. Сопоставление моделей с данными, полученными при натурных исследованиях последствий падения Сихоте-Алинского метеорита, показывает, что при

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

Анализ зависимости интенсивности излучения болида Шумава от высоты привел к выводу, что это тело размером порядка 1 м, имеющее скорость входа 27 км/с, стало разрушаться на высотах 80-85 км. Если учесть зависимость прочности от размера и то обстоятельство, что обтекание космического тела могло происходить в вихревом режиме с увеличенными осциллирующими нагрузками, оказывается, что прочность тела не выходила за рамки прочности типичных каменных тел. Фрагментация на определенных высотах, где достигается наибольшая эффективность теплопередачи излучением, привела к резкому росту коэффициентов абляции и усилению интенсивности излучения. В отличие от других работ, где предполагалась кометная природа этого тела, построенная модель объясняет поведение болида как каменного тела.

Численное моделирование падения фрагмента кометы Шумейкер-Леви 9 радиусом 1 км в атмосфере Юпитера продемонстрировало, что аналитические модели расплющивающегося тела вполне применимы к этому случаю. Вместе с тем, оказалось, что имеется различие по сравнению с аналогичными результатами моделирования других авторов. Все расчеты предсказывают основное энерговыделение ниже облачного покрова Юпитера, однако при этом отличаются как темпы и глубина выделения энергии, так и детали процесса разлета фрагментов. Это объясняется тем, что неустойчивые возмущения по-разному вносятся и усиливаются в зависимости от используемого алгоритма.

Рис. 5. Рой осколков Тунгусского метеороида на стадии его полной дезинтеграции. Т -высота от поверхности Земли, У горизонтальная координата. Черными кружками изображены частицы метеороида, точками - частицы воздуха. Каменный метеороид, 58 м в диаметре, в расчетах состоял из 650 частиц, 340 из них остались в области, охваченной рисунком. Рой осколков в этот момент имеет среднюю скорость 10 км/с, в то время как начальная скорость метеороида предполагалась равной 15 км/с. Впереди роя летят более массивные частицы.

Моделирование взаимодействия Тунгусского космического тела с атмосферой лагранжевым методом показало, что на стадии сильного

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

Л Вт/стер Р1Р.

I, с tit.

Рис 6 Мощность излучения в единицу телесного угла в зависимости от времени при наблюдении вспышки гипотетического Тунгусского падения с большой высоты (а), б -сравнение радиационного импульса в Тунгусском событии со световыми импульсами, зарегистрированными со спутников в событиях 1 10 1990, 5-8 кт, (пунктир) и 1 02 1994, 3040 кт, (штриховая линия) в безразмерных переменных Мощность излучения отнесена к ее максимальному значению Р>, а время - к величине /« = £»/Р>, где Е> - интеграл от мощности излучения по времени от начала импульса до его максимума

Была вычислена мощность излучения (световой импульс) Тунгусского метеороида в предположении каменного тела диаметром 58 м, вошедшего в атмосферу под углом 45° со скоростью 15 км/с (рис 6) Полная доля энергии, уносимой световым излучением, составляет около 20% от начальной кинетической энергии метеороида На рис 6 приводятся также световые кривые, зарегистрированные датчиком со спутника при падении довольно крупных метеорных тел, но существенно меньшего размера, чем Тунгусский астероид Участок роста интенсивности излучения в зарегистрированных световых импульсах неплохо согласуется с расчетным для Тунгусского явления, но спад мощности по времени в зарегистрированных падениях

происходит резче, поскольку разрушение и максимальное энерговыделение происходили на больших высотах 30 км в событии 1 10 1990 и 21 км в событии 1 02 1994 Вычисление потоков излучения на поверхности Земли при падении Тунгусского космического тела показало, что они существенно зависят от состояния атмосферы, и при видимости 20-40 км согласуются с данными натурных исследований

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

Глава 4. Касательные удары километровых тел

Рассматриваются такие падения космических тел, когда они входят в атмосферу под очень острыми углами, либо их траектории вообще не касаются поверхности Земли, но проходят через достаточно плотные слои атмосферы Вычисленная вероятность таких событий не является пренебрежимо малой - на протяжении фанерозоя произошло 25-50 касательных ударов километровых тел В этих случаях тело проходит большой путь в атмосфере и поэтому может воздействовать тепловым излучением на протяженную область земной поверхности Касательный удар аналогичен вертикальному удару тела с меньшим размером В п 4 2 на основе аналитической модели расплющивающегося метеороида вычисляются потоки излучения на поверхности Земли и исследуется возможность инициирования пожаров в случае касательных ударов Оценки показали, что в фанерозое несколько касательных ударников километровых размеров пролетали над континентальной поверхностью Земли и могли поджечь леса на площадях от 106 до 2 10б км2 Длина зоны начального воспламенения могла составлять от 1000 до 3000 км, а ширина - несколько сотен километров Но относительно узкая область начального воспламенения могла затем развиться в пожары континентального масштаба

В и 4 3 численное моделирование применяется к тем случаям касательных ударов, когда излучение, испускаемое из объема нагретого воздуха и пара, плавит и испаряет грунт Вычисляются толщины расплавленного слоя кварцевого песка

-20 -15 -10

-50-45-40-35-30-25-20-15-10 -5 0 5 Горизонтальная координата, км

Горизонтальная координата, км

Рис. 7. Положение вещества астероида (частицы изображены черными кружками) и изолинии плотности в последовательные моменты времени при движении астероида в атмосфере. Первоначально астероид имел кубическую форму с ребром, равным 3 км. Скорость тела на бесконечности 1 км/с (в перигее 11.2 км/с), а траектория проходит близко от поверхности Земли. Начало горизонтальной координаты совпадает с передним краем тела.

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

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

Расстояние вдоль траектории, км

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

С помощью гипотезы воздушного взрыва объясняется возникновение стекол Ливийской пустыни, площадь рассеяния которых составляет лишь около 7000 км2. Ударных кратеров в окрестностях этого поля рассеяния не найдено. Как показали расчеты, ливийские стекла могли быть образованы при ударах тел размером в сотни метров под острыми углами к поверхности (рис. 9).

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

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

Расстояние вдоль проекции траектории, км

Рис. 9. Изолинии энергии излучения, приходящейся на единицу площади поверхности, после падений кометоподобных тел размером 300 м, входящих в атмосферу со скоростью 50 км/с под углами 10е (а) и 20° (б). Значения радиационной энергии в кДж/см2 указаны у изолиний. Начало координат помещено в точку с максимальной плотностью поверхностной энергии. Тела двигались справа налево.

От 1.5 до 3.5% сталкивающихся с Землей комет, размером 5—10 км, имели орбиты с перигейными расстояниями близкими к радиусу Земли, при которых они существенно тормозились в атмосфере (при р= 10 бар) без значительных потерь массы. Мягкое приземление этих крупных комет создавало высокие

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

ц, кг/год

10|0г

Рис 10 Средний поток кометного ОВ на Землю как функция времени (отсчитываемого в прошлое) для гипотетических атмосфер с давлениями от 3 до 100 бар Поток кометного ОВ 4 109 лет назад мог составлять от 3 107 до 2 108 кг/год при атмосферных давлениях от 10 до 100 бар, что вполне сравнимо с эндогенными источниками и потоком, приносимым межпланетной пылью Основной вклад в средний поток кометного ОВ на Землю осуществлялся за счет тед размером порядка километра Стометровые тела полностью испарялись в плотной атмосфере, а большинство тел размером 10 км достигали поверхности со слишком высокой скоростью

38 39 40 41 42 43 44 Л млрд лет

Глава 5. Сильные акустико-гравитационные волны в атмосфере, инициируемые ударами крупных (~10 км) космических тел

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

Результаты расчетов проиллюстрированы на рис 11 Движение вверх приводит к выбросу некоторой части газа на высоты более 100 км, а затем падению этого газа сверху на атмосферу За время около 1 часа формируется внутренняя волна - аналог волны Лэмба, которая распространяется от точки удара вдоль поверхности Земли в нижней атмосфере Примерно за 5 часов возмущения, связанные с волнами, распространяющимися с большей скоростью в средней атмосфере, достигают точки, диаметрально противоположной точке удара Но эти возмущения не вызывают существенного движения в нижних слоях атмосферы перед фронтом внутренней волны, движущейся вдоль поверхности Через сутки после удара изменения плотности в нижних слоях атмосферы уже невелики по сравнению с невозмущенной атмосферой На высотах более 40 км плотность воздуха оказывается в

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

<Л

8 3

1 со

30 минут

Горизонтальная координата, км Горизонтальная координата, км Гориэонтальная координата, км

22 часа

5 часов

7 часов

10 часов

Рис. 11. Контуры относительной плотности воздуха в определенные моменты времени в процессе распространения волн после выделения энергии 2.5- 1022 Дж (5% от энергии каменного астероида диаметром 15 км при скорости 15 км/с) в нижних слоях атмосферы в области с характерным горизонтальным масштабом ¡000 км. В зонах белого цвета относительная плотность воздуха ниже 0.5, в серых - от 0.5 до 1.5 (включая нормальную плотность), а в темно-серых зонах - выше 1.5 (плотность отнесена к значениямплотности воздуха в невозмущенной атмосфере). Изолинии построены в интервале относительных плотностей ог 0.5 до 1.5 с шагом 0.05. Горизонтальная координата отсчитывается вдоль поверхности Земли от точки удара до диаметрально противоположной точки.

Волна Лэмба имеет большую амплитуду. Скорости газа у поверхности Земли показаны на рис. 12. При выделении энергии более 2- 10" Дж в области нижних слоев атмосферы с характерным горизонтальным масштабом порядка 1000 км скорости частиц воздуха за фронтом больше 30 м/с. Такие скорости ветра могут привести к частичному вывалу и обламыванию ветвей деревьев по всему Земному шару. Указанное выделение энергии может произойти при падении на Землю астероида диаметром 15 км с типичной скоростью 15 км/с. Кратер Чиксулуб в Мексике возрастом 65 млн. лет мог быть создан как раз при таком ударе. По-видимому, для поражающего действия волны Лэмба это критический размер тела, так как при падении меньших тел достаточно

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

Горизонтальная координата, 108 м

Рис. 12. Зависимость скорости воздуха у поверхности Земли от расстояния от места удара в отдельные моменты времени, которые в часах указаны у соответствующих кривых. Вариант расчета соответствует рис. 11.

Глава 6. Эрозия и пополнение атмосферы за счет ударов космических

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

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

В диссертационной работе проводилось численное моделирование вертикальных ударов космических тел размером от 100 м до 10 км и рассчитывались массы атмосферного газа и вещества ударника, которые преодолевают притяжение планеты На основании расчетов и оценок были построены аппроксимационные формулы доли теряемой массы как функции диаметра тела, скорости, угла входа и плотности атмосферы Показано, что тела в несколько сотен метров, сильно фрагментирующие при падении в атмосфере, более эффективны с точки зрения уноса атмосферы, приводя к потерям атмосферы из ее нижних слоев Тем не менее, согласно разработанной модели конечные параметры атмосферы существенно зависят от того периода роста планеты, когда дефицит слоя аккумулируемого Землей вещества ЗЬ изменялся от -100 до -10 км и, следовательно, еще происходили удары очень крупных тел (дЪ - толщина слоя поверхности растущей планеты, который ей предстоит нарастить, чтобы достигнуть окончательного размера) Более ранний период роста планеты, когда <5/г>100 км (т е Земля набрала менее 97 5% своей массы), практически не влияет на конечное состояние атмосферы

С использованием характерных распределений ударников по массам и по углам были построены зависимости масс ранних атмосфер Земли и Марса от времени При этом делались различные предположения о составе ударяющих тел Варьировались содержание атмофильных элементов в ударяющих телах, показатель степени ц в распределении ударяющих тел по массам (в большинстве вариантов предполагалось <7=1 8) и время начала кометпой бомбардировки Соотношение массы астероидов, короткопериодических и долгопериодических комет принималось равным 0 82 0 12 0 Об на Земле и 0 94 0 04 0 02 на Марсе

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

При моделировании ударов было определено количество воды, выбрасываемой в верхние слои современной атмосферы при ударах тел размером порядка 1 км в глубоководные бассейны. Разложение воды под действием солнечного излучения достаточно быстро приводит к образованию водорода, который затем составляет верхнюю часть атмосферы и постепенно проникает сверху в стратосферу, разрушая озоновый слой. Изменение концентрации водорода в стратосфере, где находится озоновый слой, со временем оценивалось путем решения уравнения диффузии для малой компоненты; учитывалась утечка Джинса в верхней разреженной части атмосферы и сток водорода в поверхность планеты. Была оценена связанная с выбросом воды в верхнюю атмосферу опасность для современной Земли. В случае удара тела диаметром три километра заметное понижение концентрации озона в стратосфере может продолжаться в течение нескольких десятилетий.

Глава 7. Удары очень крупных (до -1000 км) астероидов и комет по ранней Земле, временные силикатные атмосферы

,15!-<-.-1-.->- .-¡.51-г-1-;-----1-----.---- -1.51--' •'■'» ---^--'-----

-1.5 -1.0 -0.5 00 05 1.0 1 5 -1.5 -1.0 -0.5 00 0.5 1.0 1.5 -1.5 -1 0 -0 5 0.0 0.5 1.0 1.5

х/КЕ Х/Ке Х/Я*

Рис. 13. Удар астероида диаметром 3000 км со скоростью 15 км/с. В правых частях рисунков показаны изолинии постоянной плотности, которые построены с интервалом 0.5 г/см3, изолиния с минимальной плотностью находится у поверхности Земли. На границе мантии и ядра, где происходит скачок плотности от 5.9 г/см3 до 10.3 г/см3, девять изолиний совпадают. В левых частях рисунков показано движение вещества Земли и астероида. Линии являются границами между определенными слоями внутри Земли, сферическими до удара. Частицы астероида показаны мелкими темными кружками. ЯЕ - радиус Земли.

Моделировались удары тел диаметром от 50 до 3000 км по ранней Земле, которые приводят к перестройке и частичному испарению коры и мантии, испарению океана и образованию временных силикатных и водяных атмосфер. Земля рассматривается как тело сферической формы с определенным строением - корой, мантией и ядром. Пример расчета удара наиболее крупного тела приведен на рис. 13. Астероид проникает в Землю до глубины, приближенно равной его диаметру, образуя через 10 минут переходный кратер. Затем сжатые слои Земли поднимаются вверх внутрь полости переходного кратера, который заполняется расширяющимся веществом мантии и коры. Плотное вещество поднимается до высоты 1500 км над первоначальной поверхностью Земли через 30 мин, а затем в течение часа падает вниз. Рассчитывались распределения энергии, параметры слоя выбросов и теряемые после удара массы.

Рис. 14. Образование силикатной атмосферы после удара астероида диаметром 3000 км со скоростью 15 км/с. В правых частях рисунков показаны изолинии постоянной плотности, которые построены с интервалом 10°3 (в логарифмической шкале) от 10"8 г/см3 (внешние изолинии) до 10 г/см3 (внутренние изолинии, расположенные на поверхности внутреннего ядра Земли). Цифры на кривых указывают десятичные логарифмы плотности. В левых частях рисунков показаны изолинии температуры, измеряемой в тысячах градусов Кельвина

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

атмосферы показан на рис. 14. Время жизни атмосферы определяется излучением верхнего слоя пара в космическое пространство, которое приводит к охлаждению и конденсации пара. Предполагалось, что атмосфера приходит в состояние равновесия, и на этой стадии решалась одномерная сферически симметричная гидродинамическая задача с излучением. Сконденсировавшееся вещество на внешней границе атмосферы медленно падает к Земле со скоростью около нескольких метров в секунду. Атмосфера из силикатного пара, образованная после удара астероида диаметром 3000 км, продолжает существовать в течение примерно 10 лет, а после удара тела размером 500 км -всего три дня. Вероятность дополнительных ударов крупных тел за этот период времени очень мала. Ударники размером более 500 км нагревают поверхность Земли до температур в несколько сотен градусов Кельвина и более и испаряют водяной океан. Время жизни водяной атмосферы определяется излучением водяного пара и примерно одинаково для всех ударов, испаряющих океан. Оно составляет около 3000 лет для океана глубиной 3 км.

Рис. 15. Сечение Земли вертикальной плоскостью, проходящей через место удара. Заштрихованы области расплавленного вещества после ударов астероидов диаметрами 200 км (а) и 400 км (б). Яе - радиус Земли.

В работе моделировались также удары космических тел размером 50-400 км, которые могли произойти в архее в тот период, когда кора Земли уже полностью сформировалась. Геологическая активность и эрозия быстрее уничтожают мелкие кратеры, поэтому в архее более вероятно найти следы самых крупных из возможных ударов. Вычислялись размеры кратеров, количество испаренного и расплавленного вещества коры и мантии, толщины слоя сконденсировавшегося вещества. Объемы, занимаемые расплавленной породой через 1 час после удара, показаны на рис. 15. Максимальная глубина,

где мантия оказывается полностью расплавленной, изменяется от ~100 км при диаметре ударника 50 км до -450 км при размере тела 400 км Результаты расчетов аппроксимировались формулами Полученные данные можно использовать при поиске следов ударов крупных космических тел на древних платформах Основными признаками крупных ударов в архее могут служить большие объемы переплавленного вещества, особенно ультраосновного состава, и относительно толстые поверхностные слои сконденсированного вещества, включающие сферулы высокоскоростных выбросов Усредненная по всей поверхности Земли толщина слоя выбросов составляет около 10 см при ударе тела размером 50 км и 1 м при размере тела 100 км

Основные выводы по результатам исследований

1 Развиты численные методы решения уравнений газовой динамики для моделирования ударов и атмосферных течений, вызываемых вторжением космических тел в атмосферу Разработан оригинальный свободно-лагранжев метод расчета Эйлеров метод БОУА модифицирован для двумерных расчетов в сферической системе координат Разработан экономичный метод расчета одномерных уравнений переноса излучения, который позволяет учитывать эффекты переноса тепла и потери энергии на излучение в геофизических задачах, связанных с возникновением областей нагретого газа

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

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

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

4 Проведены расчеты деформации и фрагментации тел размером 50100 м в атмосфере Земли, которые, будучи сильно раздробленными, ведут себя подобно несжимаемой жидкости С помощью численного моделирования проверена теория увеличения поперечного сечения метеороида при падении в атмосфере, которая была выдвинута ранее другими авторами на основании простых оценок Показано, что существенную роль в фрагментации раздробленного метеороида диаметром 10-100 м играют неустойчивости Релея-Тейлора и Кельвина-Гельмгольца, развивающиеся на поверхности тела Вследствие развития крупномасштабных неустойчивостей тело в полете может принимать как расплющенную форму с отверстиями, так и вытянутую хорошо обтекаемую форму Высота конечного торможения тел с фиксированными массой и составом не является вполне определенной величиной вследствие случайного характера развития неустойчивостей, а также различной формы падающих тел

5 Проведены расчеты потоков излучения внутри и вне светящейся области, создаваемой раздробленным каменным космическим телом, параметры которого соответствуют Тунгусскому событию 1908 года (размер 50-70 м при скорости 15-20 км/с) Показано, что при дроблении тела на осколки размером менее 10 см, эти фрагменты могут полностью испариться Дробление на столь малые фрагменты полностью согласуется с аэродинамическими нагрузками и типичной прочностью метеоритов Подтверждена выдвинутая ранее гипотеза о том, что взрывы в атмосфере, подобные Тунгусскому, с отсутствием кратеров и фрагментов метеоритов на поверхности грунта, могут происходить при падении каменных тел размером в десятки метров

6 Исследованы касательные удары, при которых космическое тело подходит к Земле по гиперболической траектории с перигеем, отстоящим от 0 до нескольких десятков километров от поверхности планеты Определено, что кометы диаметром 1-10 км, движущиеся по касательным траекториям, могут,

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

7 Исследованы сценарии ударов космических тел, при которых могут образоваться природные стекла за счет плавления груша излучением Показано, что слоистые тектиты на площади 700000 км2 в Юго-Восточной Азии могли образоваться в результате касательного удара астероида диаметром порядка 1 км, который дробится и рассеивается в атмосфере Природные стекла с меньшим полем рассеяния, типа ливийских стекол (площадь 7000 км2), могут образоваться при ударах тел размером в сотни метров под острыми углами к поверхности Земли Излучение из образующегося воздушного взрыва производит плавление сантиметровых слоев кремнезема

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

9 Показано, что удары тел размером более 10 км генерируют в атмосфере Земли сильные акустико-гравитационные волны Идущая вдоль поверхности акустико-гравитационная волна типа волны Лэмба имеет столь большую амплитуду, что скорость газа за ее фронтом соответствует ураганному или штормовому ветру, причем эта волна огибает всю поверхность Земли Обосновано предположение, что действие поверхностной акустико-гравитационной волны после удара, образовавшего кратер Чиксулуб 65 млн лет назад, могло вызвать вывал леса или обламывание ветвей деревьев по всей Земле Большое количество высохшего древесного материала могло послужить причиной распространения глобальных пожаров, которые усугубили

стрессовую ситуацию в биосфере и способствовали гибели биоты на границе мелового и палеогенового периодов

10 Исследованы вертикальные удары крупных (50-1000 км) астероидов и комет по ранней Земле Показано, что значительная доля энергии самых крупных ударников тратится на перемещение мантийного материала, который поднимается вверх и растекается по поверхности планеты Предложены аппроксимационные формулы для массы расплавленного и испаренного вещества мантии, коры и ударника Оценена толщина слоя вещества, выпадающего на Землю после удара

Необходимо отметить, что ряд результатов был получен автором диссертации совместно с сотрудниками лаборатории математического моделирования геофизических процессов ИДГ РАН Разработка свободно-лагранжева метода (глава 1) проводилась совместно с В М Хазинсом Исследование падений метеороидов в атмосфере (глава 3) проводилось в тесном взаимодействии с И В Немчиновым, В В Шуваловым, А В Тетеревым, О П Поповой, А П Голубем, Н А Артемьевой Результаты по сильным акустико-гравитационным волнам (глава 5) получены совместно с В В Шуваловым При построении аппроксимационных формул потерь атмосферного газа после ударов (глава 6) использовались результаты расчетов косых ударов Н А Артемьевой и В В Шувалова Исследование по предсказанию возможных следов ударов в архее (глава7) проведено вместе с В В Шуваловым

Основные результаты диссертации опубликованы в следующих работах

1 Светцов В В Экономичный численный метод для одномерных задач радиационной газодинамики // Журнал вычислительной математики и математической физики 1994 Т 34 №3 С 432-445

2 Светцов В В Взрывы в нижней и средней атмосфере - сферически симметричная стадия//Физика горения и взрыва 1994 Т 30 № 5 С 129142

3 Светгрв В В Падение кометы в атмосфере Юпитера // Астрономический вестник 1995 Т 29 №4 С 331-340

4 Немчинов ИВ, Попова ОП, Светцов ВВ, Шувалов В В О фотометрической массе и радиационном размере крупных метеороидов //Астрономический вестник 1995 Т 29 № 2 С 155-173

5 Светцов В В Куда делись осколки Тунгусского метеороида"? // Астрономический вестник 1996 Т 30 №5 С 427—441

6 Светцов В В Взрывы метеороидов и оценка их параметров по световому излучению//Физика горения и взрыва 1998 Т 34 №4 С 117-128

7 Светг/ов В В Загадки кратерного поля Сихотэ-Алиня // Астрономический вестник 1998 Т 32 № 1 С 76-88

8 Светцов В В Нестационарный режим сверхзвукового обтекания // Журнал технической физики 1999 Т 69 №12 С 95-97

9 Светцов В В Нарушение стандартного режима обтекания и болид Шумава//Астрономический вестник 2000 Т 34 №4 С 331-347

10 Светцов В В Об эффективности ударного механизма эрозии атмосферы // Астрономический вестник 2000 Т 34 № 5 С 441-453

11 Светг/ов В В Оценка привноса органического вещества кометами на Землю на начальных стадиях ее эволюции // Астрономический вестник 2002 Т 36 № 1 С 55-67

12 Светцов В В Пожары и ураганные акустико-гравитационные волны от ударов крупных космических тел // Динамика взаимодействующих геосфер Сборник научных трудов ИДГ РАН Москва, 2004 С 258-267

13 Светцов В В Тепловое излучение и пожары при падениях космических тел И Катастрофические воздействия космических тел / Ред Адушкин В В , Немчинов И В Москва ИКЦ "Академкнига", 2005 С 151-165

14 Светцов В В Выброс воды в верхнюю атмосферу при ударах космических тел в океан возможность разрушения озонного слоя // Всероссийская конференция "Астероидно-кометная опасность - 2005" (АКО-2005) Материалы конференции СПб ИПА РАН, 2005 С 290293

15 Светцов В В, Шувалов В В Сильные акустико-гравитационные волны от ударов крупных космических тел // Известия РАН, Физика атмосферы и океана 2005 Т 41 № 5 С 602-612

16 Светцов В В, Шувалов В В Тунгусская катастрофа 30 июня 1908 г // Катастрофические воздействия космических тел / Ред Адушкин В В, Немчинов И В Москва ИКЦ "Академкнига", 2005 С 167-199

17 Светцов В В, Шувалов В В Воздействие ударов крупных космических тел на поверхность ранней Земли // Доклады РАН 2005 Т 403 № 6 С 812-814

18 Светцов В В Выброс воды в атмосферу при ударах крупных космических тел в океан возможность разрушения озонного слоя // Динамические процессы во взаимодействующих геосферах Сборник научных трудов ИДГ РАН Москва ГЕОС, 2006 С 247-250

19 Светцов В В Потерн и пополнение атмосферы за счет ударов космических тел на Земле и Марсе // Астрономический вестник 2007 Т 40 Ш С 33-46

20 Светцов В В Оценки энергии поверхностных волн при взрывах в атмосфере и параметров источника Тунгусского события // Физика Земли 2007 Т 43 № 7 С 57-66

21 Nemchmov I V, Svetsov VV Global consequences of radiation impulse caused by comet impact // Advances in Space Research 1991 V 11 No 6 P (6)95-(6)97

22 Ivanov В A, Nemchinov IV, Svetsov VV, Provalov AA, Khazins VM Impact cratermg on Venus Physical and Mechanical Models // Journal of Geophysical Research 1992 V 97 No ЕЮ P 16167-16181

23 Hazins VM, Svetsov VV A conservative stable smoothness-enhancing free-Lagrangian method // Journal of Computational Physics 1993 V 105 No 2 P 187-198

24 Svetsov V V Radiation emitted during the flight application to assessment of bolide parameters from the satellite recorded light flashes // Lunar and Planetary Science Conf XXV Houston LPSI, 1994 P 1365-1366

25 Nemtchmov I V, Popova OP, Shuvalov VV, Svetsov VV Radiation emitted during the flight of asteroids and comets through the atmosphere // Planetary and Space Science 1994 V 42 No 6 P 491-506

26 Svetsov V V, Nemtchmov I V, Teterev A V Disintegration of large meteoroids in Earth's atmosphere Theoretical models//Icarus 1995 V 116 No 1 P 131-153

21 Nemtchmov IV, Kosarev IB, Popova OP, Shuvalov VV, Svetsov VV, Spaulding RE, Jacobs С, Shavez J, Tagliaferry E Historical evidence of recent impacts on the Earth // Proceedings of the Planetary Defense Workshop Lawrence Livermore National Laboratory, Livermore, California, May 22-26, 1995 P 31-49

28 Orlova TI, Svetsov V V Deep penetration of projectile material into the planet during the impact//Earth, Moon, and Planets 1995 V 71 P 225-263

29 Svetsov VV Total ablation of the debns from the 1908 Tunguska explosion II Nature 1996 V 383 No 6602 P 697-699

3d Nemtchmov IV, Svetsov VV, Kosarev IB, Golub' AP, Popova OP, Shuvalov VV, Spalding RE, Jacobs С, Tagliaferri E Assessment of kinetic energy of meteoroids detected by satellite-based light sensors // Icarus 1997 V 130 No 2 P 259-274

31 Nemtchmov IV, Artem'eva NA, Kosarev IB, Shuvalov VV, Svetsov VV, Neukum G, Hahn G, De Niem D Luminosity of the bolides created by SL-9 comet fragments in the Jovian atmosphere // International Journal of Impact Engineering 1997 V 20 No 6-10 P 591-599

32 Nemtchmov IV, Shuvalov VV, Kosarev IB, Artem'eva NA, Trubetskaya IA , Svetsov V V, Ivanov В A , Loseva T V, Neukum G, Hahn G, de Niem D

Assessment of Comet Shoemaker-Levy 9 fragment sizes using light curves measured by Galileo spacecraft instruments // Planetary and Space Science 1997 V 45 No 3 P 311-326

33 Svetsov V V Could the Tunguska debris survive the terminal flare7 // Planetary and Space Science 1998 V 46 No 2/3 P 261-268

34 Svetsov V V Vortical regime of the flow behind the bow shock wave // Shock Waves 2001 V 11 No 3 P 229-244

35 Svetsov V V Grazing meteoroids could ignite contmental-scale fires // Catastrophic events and mass extinctions Impacts and beyond / Eds C Koeberl, K G MacLeod Bolder, Colorado, Geological Society of America Special Paper 356 2002 P 685-694

36 Svetsov VV Assessment of organics delivery by comets to the early Earth // Lunar and Planetary Science Conf XXXIII (Houston, Texas), 2002 Abstract #1451

37 Svetsov V V Numerical modelling of large asteroidal impacts on the Earth // International Journal of Impact Engineering 2003 V 29 P 671-682.

38 Svetsov VV Numerical simulations of very large impacts on the Earth // Planetary and Space Science 2005 V 53 No 12 P 1205-1220

39 Svetsov V V Impact erosion of atmosphere Some results of numerical simulations for vertical impacts // Lunar and Planetary Science Conf XXXVI (League City, Texas), 2005 Abstract #1675

40 Svetsov V V Thermal radiation on the ground from large aerial bursts caused by Tunguska-hke impacts // Lunar and Planetary Science Conf XXXVII (League City, Texas), 2006 Abstract #1553

41 Svetsov V V, Wasson J T Melting of soil rich in quartz by radiation from aerial bursts - a possible cause of formation of Libyan Desert Glass and layered tektites // Lunar and Planetary Science Conf XXXVIII (League City, Texas), 2007. Abstract # 1499

Подписано в печать 14 июля 2008 г. Объем 1,8 п л Тираж 120 экз Заказ № 624 Отпечатано в Центре оперативной полиграфии ООО «Ол Би Принт» Москва, Ленинский пр-т, д 37

Содержание диссертации, доктора физико-математических наук, Светцов, Владимир Владимирович

Введение.

1 Развитие методов математического моделирования течений, возникающих при ударах космических тел.

1.1 Свободно-лагранжев метод.

1.1.1 Построение сетки.

1.1.1.1 Триангуляция.

1.1.1.2 Дополнительные связи.

1.1.1.3 Процедура преобразования (меппинг).

1.2.1 Конечно-разностная схема.

1.2.1.1 Разностные уравнения.

1.2.1.2 Процедура обратного преобразования (ремеппинг).

1.2.1.3 Итерационная процедура.

1.2 Метод 80УА.

1.2.1 Применение метода 80УА.

1.2.2 Метод БОУА в сферических координатах.

1.3 Другие методы решения уравнений газовой динамики.

1.3.1 Методы, используемые в задаче обтекания метеороида.

1.3.2 Тест на карбункулярную неустойчивость.

1.4 Экономичный метод расчета переноса излучения.

1.4.1 Приближенные уравнения.

1.4.2 Разностная схема.

1.4.3 Многогрупповое приближение.

Выводы по главе

2 Вихревой режим обтекания космических тел в атмосфере.

2.1 Неустойчивость ударной волны и течения.

2.2 Тестовая задача.

2.3 Численные решения тестовой задачи в рамках уравнений Эйлера . 67 2.4. Решения на основе уравнения Навье-Стокса.

2.5 Применимость численных методов.

2.6 Изменение показателя адиабаты и геометрии.

2.7 Учет переноса излучения.

2.8 Геофизические аспекты.

Выводы по главе 2.

3 Падения разрушенных метеороидов в атмосфере Земли и их особенности (тела размером 1—100 м).

3.1 Модели движения разрушенных тел.

3.1.1 Роль абляции при движении метеороидов в атмосфере.

3.1.2 Прочность и начало разрушения.

3.1.3 Аналитические и полуаналитические модели разрушения и движения метеороидов.

3.1.4 Гидродинамические численные модели.

3.2 Особенности Сихоте-Алинского метеоритного дождя.

3.3 Вспышки в атмосфере.

3.4 Болид Шумава.

3.5 Комета Шумейкер-Леви 9 и Тунгусское космическое тело.

3.6 Воздействие ударных волн на поверхность Земли - сейсмический . эффект при взрыве метероида в воздухе.

Выводы по главе 3.

4 Касательные удары километровых тел.

4.1 Касательные удары и их специфические черты.

4.2 Излучение и возникновение пожаров.

4.2.1 Проблема глобальных пожаров.

4.2.2 Модель расплющивающегося тела при касательном ударе.

4.2.3 Результаты вычислений.

4.3 Плавление поверхности и образование природных стекол.

4.3.1 Плавление и испарение поверхности при касательных ударах

4.3.2 Образование тектитов при ударах космических тел.

4.3.3 Плавление кремнезема.

4.3.4 Прямое численное моделирование касательных ударов.

4.3.5 Прямое численное моделирование ударов под острыми углами.

4.4 Доставка органического вещества и возникновение жизни.

4.4.1 Проблема доставки органического вещества (ОВ) на Землю.

4.4.2 Приближенные формулы для потока массы ОВ на Землю.

4.4.3 Торможение и абляция тела в атмосфере.

4.4.4 Приземление.

4.4.5 Результаты вычислений.

Выводы по главе 4.

5 Сильные акустико-гравитационные волны в атмосфере, инициируемые ударами крупных (~10 км) космических тел.

5.1 Последствия удара космического тела на границе мел-палеоген

5.2 Постановка задачи о генерации акустико-гравитационных волн при ударе крупного космического тела.

5.3 Особенности моделирования сильных АГВ.

5.4 Результаты моделирования.

Выводы по главе 5.

6 Эрозия и пополнение атмосферы за счет ударов космических тел

6.1 Ранняя атмосфера и механизмы ее эрозии.

6.2 Модели ударной эрозии атмосферы.

6.3 Оценки потерь атмосферы при ударах тел диаметром 100 м - 1 км

6.4 Численное моделирование вертикальных ударов и потери атмосферы для современной Земли.

6.5 Численное моделирование ударов по Марсу.

6.6 Аппроксимационные формулы.

6.7 Эволюция атмосфер планет земной группы под действием ударов

6.8 Выброс воды при ударах в океан: возможность разрушения озонного слоя.

Выводы по главе 6.

7 Удары очень крупных (до —1000 км) астероидов и комет по ранней Земле, временные силикатные атмосферы.

7.1 Размеры наиболее крупных ударников на ранней стадии эволюции Земли.

7.2 Особенности моделирования ударов крупных тел.

7.3 Воздействие крупных ударов на Землю.

7.4 Распределение энергии.

7.5 Образование и исчезновение силикатных атмосфер.

7.6 Испарение океана.

7.7 Ударная переработка поверхности Земли и поиск следов ударов в архее.АОУ

Выводы по главе 7.^

Введение Диссертация по наукам о земле, на тему "Взаимодействие космических тел с атмосферой и поверхностью Земли"

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

С другой стороны, удары крупных тел истребляли живые существа. Астероиды и кометы размером порядка десяти километров могли приводить к глобальным катастрофам и массовым вымираниям видов в фанерозое. Нобелевский лауреат Луис Альварес с коллегами в 1980 г. связал падение крупного космического тела с глобальной катастрофой на границе мел-палеоген. Позже был найден и соответствующий кратер Чиксулуб в Мексике, созданный, как показали исследования, ударом тела диаметром 10-15 км. Хотя прямая связь этого события с гибелью биоты не доказана, результаты многочисленных работ по последствиям падений космических тел показывают, что такой удар мог привести к глобальной катастрофе.

Угрожающая человечеству астероидно-кометная опасность все больше осознается с течением времени. В 1981 г. НАС А начало изучение космической угрозы. Проведенные с тех пор многочисленные международные и российские конференции позволили заострить внимание на этой проблеме и выработать ряд рекомендаций. В США была развернута программа астрономических наблюдений за потенциально опасными космическими объектами, которые могут приблизиться к Земле, с помощью телескопов, радаров и космических аппаратов. Наблюдения и каталогизация ведутся и в России. Учет потенциально опасных объектов диаметром более километра предполагается завершить к 2010 г. Но реальную опасность для людей представляют и меньшие объекты, начиная с диаметров 50-100 м (предполагаемый размер тела, вызвавшего Тунгусскую катастрофу 1908 г.). Поэтому было поставлено целью выявление космических тел размером до ~100 м, число которых составляет многие десятки тысяч. Один из известных уже выявленных объектов - это астероид Апофис диаметром около 350 м, который сблизится с Землей в 2029 г.

Проблема астероидно-кометной опасности выдвинула много задач, которые предстоит решить [207]. В числе этих задач определение и изучение процессов, происходящих при столкновении космических тел с планетой. Исследования последствий ударов космических тел и изучение механизмов воздействия играют важную роль при разработке мер снижения опасности. Для предотвращения угрозы, при применении тех или иных средств защиты от опасных объектов очень важными являются также сведения о строении и составе космических тел. Но знания в этой области недостаточны, и особенно скудны сведения о строении кометных ядер.

Падающие космические тела взаимодействуют прежде всего с атмосферным воздухом, причем большинство метеороидов не долетает до поверхности. Они оказывают воздействие на атмосферу, тем более сильное, чем больше размер и скорость ударника. Изучение взаимодействия космических тел с атмосферой актуально не только с точки зрения состояния атмосферы после удара, но и с точки зрения определения параметров падающих на землю тел. С целью выяснения их свойств проводятся дорогостоящие космические эксперименты, такие как Star Dust (захват частиц кометы Вильде 2), Deep Impact (удар по комете Темпеля 1) или планируемый Европейским космическим агентством эксперимент Дон Кихот, в котором предполагается осуществить удар по астероиду. Но по-прежнему поток метеороидов, входящих в атмосферу Земли, несет в себе огромную информацию о составе и природе малых космических тел, которая требует расшифровки.

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

Новизна результатов, связанная как с методами расчетов, так и с постановкой задач, заключается в следующем. (1) Построена лагранжева полностью консервативная численная схема, обладающая определенными преимуществами при решении газодинамических задач, связанных с ударами космических тел. (2) Разработан новый экономичный метод численного решения одномерных уравнения переноса излучения для расчета остывания областей нагретого газа. (3) Получены и иследованы новые решения уравнений газовой динамики при обтекании тел гиперзвуковым потоком газа. (4) Проведены исследования нескольких падений космических тел: Тунгусского, Сихотэ-Алинского, болида Шумава, кометы Шумейкер-Леви 9, в которых выявлены не изученные ранее особенности этих событий. (5) Впервые полноценно исследованы касательные удары космических тел (под очень острыми углами к горизонтали), показано, какие геологические следы они могли оставить и определена их роль в снабжении Земли органическим веществом. (6) Впервые исследованы сильные акустико-гравитационные волны после ударов крупных космических тел, которые могут производить существенное механическое воздействие в глобальном масштабе. (7) Построена модель эволюции атмосферы в период завершения аккумуляции Земли. (8) Проведено математическое моделирование ударов очень крупных космических тел, приводящих к образованию временных силикатных атмосфер на ранней стадии существования планеты.

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

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

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

Методы исследования. Задачи обтекания космических аппаратов, которые решались до сих пор в связи с развитием космической техники, представляют хорошую базу для понимания процессов взаимодействия космических тел с атмосферой [176], но значительно более высокие скорости входа астероидов и комет в атмосферу Земли, фрагментация и разрушение вносят новые элементы и создают существенные трудности для создания математических моделей. При наблюдениях падений метеорных тел удается зарегистрировать, в лучшем случае, зависимость скорости тела от высоты, мощность излучения и спектр, и поэтому несмотря на обширный наблюдательный материал," поставляемый наземными болидными сетями и регистрацией со спутников, интерпретации падений метеорных тел вызывают много противоречий. Так до сих пор не удается согласовать динамические и фотометрические массы, определяемые для метеоров, а следовательно и определить их истинные массы [29, 262]. Чем крупнее тела, тем реже они сталкиваются с Землей, и тем меньше наблюдательных данных, которые могли бы служить основой для исследования взаимодействия крупных метеороидов (более 10 м) с атмосферой. Определенные сходства с выделением энергии астероидов и комет в атмосфере имеют воздушные ядерные взрывы. Но процесс торможения, наличие следа и материал космического тела вносят в эти явления существенные различия.

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

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

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

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

2. Исследовать атмосферные эффекты, вызываемые ударами космических тел при выделении их энергии в воздухе, определить параметры ударных волн, сильных акустико-гравитационных волн, потоки излучения. Существенные атмосферные эффекты возникают при падении тел размером более 50-100 м. Тела размером 50 - 100 м могут создать кратер, подобный кратеру Бэрринжера в Аризоне, если достигают поверхности, или произвести мощный взрыв в воздухе, как в Тунгусском явлении 1908 г. Они представляют прямую угрозу для людей при падении в населенных районах.

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

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

Структура и объем диссертации. Диссертация состоит из введения, семи глав и заключения. Объем диссертации - 335 страниц, включая 92 иллюстрации. Список литературы содержит 549 наименований.

Заключение Диссертация по теме "Физика атмосферы и гидросферы", Светцов, Владимир Владимирович

Основные выводы по результатам исследований

1. Развиты численные методы решения уравнений газовой динамики для моделирования ударов и атмосферных течений, вызываемых вторжением 1 космических тел в атмосферу. Разработан оригинальный свободно-лагранжев метод расчета. Эйлеров метод БОУА, предназначенный для моделирования движения нескольких веществ в разных фазовых состояниях, модифицирован для двумерных расчетов в сферической системе координат. Разработан экономичный метод расчета одномерных уравнений переноса излучения, который позволяет учитывать эффекты переноса тепла и потери энергии на излучение в геофизических задачах, связанных с возникновением областей нагретого газа.

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

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

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

4. Проведены расчеты деформации и фрагментации тел размером 50100 м в атмосфере Земли, которые, будучи сильно раздробленными, ведут себя подобно несжимаемой жидкости. С помощью численного моделирования проверена теория увеличения поперечного сечения метеороида при падении в атмосфере, которая была выдвинута ранее другими авторами на основании простых оценок. Показано, что существенную роль в фрагментации раздробленного метеороида диаметром 10-100 м играют неустойчивости Релея-Тейлора и Кельвина-Гельмгольца, развивающиеся на поверхности тела. Вследствие развития крупномасштабных неустойчивостей тело в полете может принимать как расплющенную форму с отверстиями, так и вытянутую хорошо обтекаемую форму. Таким образом, высота конечного торможения тел с фиксированными массой и составом не является вполне определенной величиной вследствие случайного характера развития неустойчивостей, а также различной формы падающих тел.

5. Проведены расчеты потоков излучения внутри и вне светящейся области, создаваемой раздробленным каменным космическим телом, параметры которого соответствуют Тунгусскому событию 1908 года (размер 50-70 м при скорости 15-20 км/с). Показано, что при дроблении тела на осколки размером менее 10 см, эти фрагменты могут полностью испариться. Дробление на столь малые фрагменты полностью согласуется с аэродинамическими нагрузками и типичной прочностью метеоритов. Подтверждена выдвинутая ранее гипотеза о том, что взрывы в атмосфере, подобные Тунгусскому, с отсутствием кратеров и фрагментов метеоритов на поверхности грунта могут происходить при падении каменных тел размером в десятки метров.

6. Исследованы касательные удары, при которых космическое тело подходит к Земле по гиперболической траектории с перигеем, отстоящим от О до нескольких десятков километров от поверхности планеты. Такие удары редки, но, как показано в работе, неоднократно происходили в истории Земли. Определено, что кометы диаметром 1-10 км, движущиеся по касательным траекториям, могут, постепенно тормозясь, совершить мягкую посадку на Землю без разрушения содержащегося в космическом теле органического материала. Этот механизм привноса органики мог послужить основой возникновения жизни на Земле. Показано, что излучение крупного метеороида, движущегося по касательной траектории может вызвать поджог горючих (древесных) материалов в полосе, длина которой сравнима с размером континента, а ширина составляет сотни километров.

7. Предложены и исследованы сценарии ударов космических тел, при которых могут образоваться природные стекла. Показано, что слоистые тектиты на площади 700000 км (где они найдены в Юго-Восточной Азии) могли образоваться в результате касательного удара астероида диаметром порядка 1 км, который дробится и рассеивается в атмосфере. Фрагменты такого тела тормозятся и выделяют свою энергию в воздухе над обширной поверхностью суши. Излучение из нагретой области обеспечивает плавление слоя кварцевого песка толщиной в сантиметры. Природные стекла с меньшим полем рассеяния, типа ливийских стекол (площадь 7000 км2), могут образоваться при ударах тел размером в сотни метров под острыми углами к поверхности Земли. Излучение из образующегося воздушного взрыва в этом случае также производит плавление сантиметровых слоев кремнезема.

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

9. Показано, что удары тел размером более 10 км генерируют в атмосфере Земли сильные акустико-гравитационные волны. Идущая вдоль поверхности акустико-гравитационная волна типа волны Лэмба имеет' столь большую амплитуду, что скорость газа за ее фронтом соответствует ураганному или штормовому ветру, причем эта волна огибает всю поверхность Земли. Обосновано предположение, что действие поверхностной акустико-гравитационной волны после удара, образовавшего кратер Чиксулуб 65 млн. лет назад, могло вызвать вывал леса или обламывание ветвей деревьев по всей Земле. Большое количество высохшего древесного материала могло послужить причиной распространения глобальных пожаров, которые усугубили стрессовую ситуацию в биосфере и способствовали гибели биоты на границе мелового и палеогенового периодов.

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

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

Необходимо отметить, что ряд результатов был получен автором диссертации совместно с сотрудниками лаборатории математического моделирования геофизических процессов ИДГ РАН. Разработка свободно-лагранжева метода (глава 1) проводилась совместно с В.М.Хазинсом. Исследование падений метеороидов в атмосфере (глава 3) проводилось в тесном взаимодействии с И.В.Немчиновым, В.В.Шуваловым, А.В.Тетеревым, О.П.Поповой, А.П.Голубем, Н.А.Артемьевой. Результаты по сильным акустико-гравитационным волнам (глава 5) получены совместно с В.В.Шуваловым. При построении аппроксимационных формул потерь атмосферного газа после ударов (глава 6) использовались результаты расчетов косых ударов Н.А.Артемьевой и В.В.Шувалова. Исследование по предсказанию возможных следов ударов в архее (глава7) проведено вместе с В.В.Шуваловым.

Апробация и внедрение результатов

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

13 и №16 "Изменения окружающей среды и климата: природные катастрофы" и №25 "Проблемы зарождения биосферы Земли и ее эволюции", в отчет по астероидной опасности для МЧС по теме "Анализ и оценка

300 последствий падения на Землю опасных космических объектов и разработка предложений по их смягчению" и монографию "Катастрофические воздействия космических тел", написанную коллективом авторов ИДГ РАН. (Вариант этой монографии на английском языке "Catastrophic events caused by cosmic objects" был издан в 2008 г. в издательстве Springer). Ряд работ вошел в отчеты, выполненные по грантам РФФИ, контрактам с Национальной лабораторией Сандиа (США) и по гранту Американского фонда исследований и развития (CRDF).

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

Результаты исследований докладывались на международных и российских конференциях: на международных конференциях по астероидной опасности (Tucson, AZ, USA, 1993; Livermore, СА, USA, 1995); на лунно-планетных конференциях (Lunar Planetary Science Conference, Houston, TX, USA, 1994, 1996, 2002, 2005, 2006, 2007); на европейском совещании по столкновению кометы SL-9 с Юпитером (Garching, Germany, 1995); на международной конференции по космосу "Space -96" (Albuquerque, NM, USA, 1996); на международных совещаниях по ударному кратерообразованию, эволюции Земли и глобальным катастрофам (Ancona, Italy, 1995; Postojna, Slovenia, 1996; Quillan, France, 1999; Vienna, Austria, 2000; Granada, Spain, 2001; Mora, Sweden, 2002; Nordlingen, Germany, 2003); на совещаниях Метеоритного общества (Annual Meetings of the Meteoritical Society, Munster, Germany, 2003; Zurich, Switzerland, 2006); на конференциях, посвященных Тунгусскому явлению 1908 г., (Bologna, Italy, 1996; Красноярск, 1998; Москва, 2003, 2008); на всероссийской конференции "Астероидно-кометная опасность 2005" (Санкт-Петербург, 2005); на международной конференции "Околоземная астрономия 2007" (Терскол, 2007).

Заключение

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

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

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

Исследование конкретных падений космических тел: Сихоте-Алинского метеоритного дождя, болида Шумава, Тунгусского космического тела, кометы Шумейкер-Леви 9, вспышек излучения, зарегистрированных со спутников, изучение возможности образования слоистых тектитов и ливийских стекол под действием излучения воздушных взрывов позволило проиллюстрировать основные процессы взаимодействия космических тел с атмосферой и в каждом случае выявить своеобразные стороны отдельных событий.

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

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

Библиография Диссертация по наукам о земле, доктора физико-математических наук, Светцов, Владимир Владимирович, Москва

1. Авдуевский B.C., Иванов A.B. Течение разреженного газа вблизи передней критической точки затупленного тела при гиперзвуковых скоростях // Изв. АН СССР. Механика жидкости и газа. 1968. № 3. С. 26-34.

2. Авилова И.В., Биберман U.M., Воробьев B.C. и др. Оптические свойства горячего воздуха. М.: Наука, 1970. 320 с.

3. Адушкин В.В., Кудрявцев В.П., Журавлева JI.A. Влияние дегазации Земли на ее озоновый слой // Динамика взаимодействующих геосфер. Сборник научных трудов. М: ИДГ РАН, 2004. С. 191-201.

4. Азовский А.Н., Шувалов В.В. Влияние геометрической формы на движение фрагментов разрушенного метеороида // Астрономический вестник. 2002. Т. 36. № 4. С. 1-7.

5. Артемьев В.И., Бергельсон В.И., Немчинов И.В., Орлова Т.И., Смирнов В.А., Хазинс В.М. Изменение режима сверхзвукового обтекания препятствия при возникновении перед ним тонкого разреженного канала //Изв. Ан СССР, МЖГ. 1989. № 5. С. 146-151.

6. Аки К., Ричарде П. Количественная сейсмология. М.: Мир, 1983. Т.1. 520 с.

7. Апштейн Э.З., Вартанян Н.В., Сахаров В.И. О распределении лучистых тепловых потоков по поверхности пространственных и осесимметричных. тел при сверхзвуковом обтекании их идеальным газом // Изв. АН СССР. Механика жидкости и газа. 1986. № 1. С. 92-97.

8. Апштейн Э.З., Пилюгин H.H., Севастъяненко В.Г., Тирский Г.А. Радиационный теплообмен при входе тел в атмосферу Земли и планет со сверхорбитальными скоростями // Итоги науки и техники. Сер. Механ. жидкости и газа. М.: ВИНИТИ, 1989. Т. 23. С. 116-236.

9. Артемьева H.A. Высокоскоростные выбросы из кратеров и обмен веществом между планетами // Катастрофические воздействия космических тел / Ред. Адушкин В.В., Немчинов И.В. М: ИКЦ "Академкнига", 2005. С. 201-212.

10. Асланов С.К. О гидродинамическом моделировании процесса абляции поверхностного слоя метеороида // Астрономический вестник. 2000. Т. 34. № 4. С. 348-356.

11. Астапович КС. Новые материалы по полету большого метеорита 30 июня 1908 г. в Центральной Сибири // Астрономический журнал. 1933. Т. 10. № 4. С. 465-486.

12. Атмосфера стандартная. ГОСТ 4401-81. М: Изд. Стандартов. 1981. 180 с.

13. Бармин A.A., Егорушкин С.А. Устойчивость поверхностей сильного разрыва в газах (обзор) // Изв. АН СССР. Механика жидкости и газа. 1996. № 2. С. 3-22.14

Информация о работе

Светцов, Владимир Владимирович
доктора физико-математических наук
Москва, 2008
ВАК 25.00.29

Диссертация

Взаимодействие космических тел с атмосферой и поверхностью Земли - тема диссертации по наукам о земле, скачайте бесплатно

Автореферат

Похожие работы