Метеорит
Страница 3
На месте катастрофы последующие экспедиции обнаружили вывал леса общей площадью 2000 км2, наблюдались светлые ночи. В районе катастрофы начался пожар и были обнаружены следы радиационного повреждения веток деревьев.
Таким образом над тайгой произошло явление взрывного типа, энергия взрыва была больше ,чем энергия взрыва 1 млн. т. троти-ла.
Работа по математическому моделированию началась в 1969 г. К этому времени уже были собраны данные о характере катастрофы.
Сейчас это исследование проводится В.П.Коробейниковым, П.И.Чушкиным и Л.В.Шуршаловым.
В дальнейшем будем придерживаться двух рабочих гипотез.
1.В атмосферу влетел фрагмент ядра кометы, окружённый пылегазовой атмосферой (комой).
2. Вторгся большой рыхлый метеорит типа углистого хондрита.
Несколько слов о головах комет и углистых хондритах. Голова кометы состоит из ядра и сильно разряжённой атмосферы (около 100 частиц/см3). Ядро кометы - это конгломерат кусков льда, газа и пыли. Средняя плотность вещества ядра не превышает 1 г/см3, давление внутри ядра размером около 1 км 1000 дин/см2. Фрагменты ядра могут соединяться в нём лишь некоторыми частями, поэтому скреплены слабо, возможно отрывание отдельных частей под действием солнечной радиации. Так ,например, ядро кометы Веста в 1976 г. разделилось на четыре фрагмента. Фрагменты могут существовать как малые кометы. По химическому составу кометы в основном содержат воду, метан, ацетилен, углекислоту, водород, соединения углерода и азота с другими элементами.
Углистые хондриты - это весьма редкий тип метеоритов, обнаруженых на Земле. Это каменные метеориты, содержащие повышенное количество углерода как свободного, так и связанного в угеводородах. В них ,как правило, имеются газовые включения и гидросодержащие минералы. Цвет - угольно-чёрный или серочёрный. Содержание воды в них может доходить до 20% (связанная вода), плотность этих метеоритов не более 3 г/см3.Только наиболее плотные и крупные из них достигают поверхности Земли, большинство же рассеивается в атмосфере. Так произошло в 1965 г. с метеоритом Ривелсток, упавшим над Канадой. Общая масса его оценивается в 4 тыс. т ,скорость входа около 12 км/с .Воздушные волны были зарегистрированы барографами за несколько тысяч километров от места падения, и общая энергия возмущения атмосферы оценена в 10-20 тыс. т тротила. Явление по мощности равно атомному взрыву над Херосимой.
Воздушные волны были зарегистрированы на ближайшей сейсмостанции и организованы поиски вещества. Однако было найдено всего около грамма вещества на льду одного озера.
Если бы космическое тело было гораздо больших размеров, чем метеорит Ривелсток, и было углистым хондритом, оно проникло бы гораздо глубже в атмосферу, и могла бы произойти катастрофа, аналогичная Тунгусской в смысле воздействия на земную поверхность.
Как кометная, так и углисто-хондритная гипотезы удовлетворяют основному свойству Тунгусского космического тела: взрывной распад над поверхностью Земли при отсутствии выпадания значительныхмасс вещества. Как кометная, так и углисто-хондритная гипотеза характерна тем, что в состав этих тел входит вода в состоянии льда, углерод и углеводороды. Все эти вещества могут либо испариться, либо сгореть в атмосфере. Кометная гипотеза более полно объясняет помутнение (запыление) атмосферы в период падения и после него, но зато падение углистых хондритов есть явление сравнительно обычное ,а столкновение с ядром малой кометы - явление уникальное.
Приведём пример численного решения задачи входа в атмос-феру Тунгусского космического тела, выполненную конечноразностным методом Л.В.Соколовской.
Газообразное тело в форме циллиндра, высота которого равна диаметру L (L=0.6 км), с начальной скоростью 40 км/с движется в атмосфере, и при t=0, ZE=36 км давление в теле равно атмосферному, плотность rme=0.1225 г/см3; gm=5/3; g=1.4, v=90 (вертикальный вход).На рис.5 показана форма тела для различных высот за время около 0.5 с. Видно, что тело начинает резко расширяться при Z< 20 км. Причём поперечный его размер меняется так:
b»b0+3×10-1vet
Заметим, что скорости бокового разлёта вещества значительны и в конце пути превышают 500 м/с .Тело тормозится до скорости 2 км/с на высоте около 10 км. Таким образом, за время порядка 1с в столбе газа длиной 29 км и толщиной 1 км выделяется энергия около 1026 эрг. По небу проходит гигантская “молния”, от которой расходится гром - след баллистической волны. В нижних слоях атмосферы при Z=13 км температура воздуха за фронтом головной волны достигает 15000 К ,давление около 30 атм. При резком торможении в концевоий части давление на фронте ударной волны падает, но и газ внутри тела, обладающий запасом внутренней энергии и значительным давлением Dp>0, начинает расширятся в окружающую среду, посылая вперед ударную волну взрывного типа, которая будет двигаться в атмосфере в направлении Земли, обгоняя частицы среды метеороида. В действительности, конечно процесс гораздо сложнее, но некоторые общие качественные черты уже улавливающая в этой простой модели.
Пользуясь изложенными выше законами, можно выполнить решение задачи о входе в атмосферу газообразного тела других размеров и энергий.
Вернемся, однако, к Тунгусскому телу. Моделирование процесса его взаимодействия с атмосферой и земной поверхностью проводилось в рамках математической модели, описанной в конце предыдущего раздела. Сначала были проанализированы результаты расчёта обтекаемых твёрдых недеформируемых тел совершенным газом при больших числах Маха M=v/a1 где a12=gp1/s1, p1, s1 - параметры окружающего воздуха. Были так же проведены специальные расчёты такого обтекания при M > 5. В результате этих расчётов определилась как форма ударных волн, так и всё поле течений газа при стационарных условиях обтекания. Оказалось, что для M > 10 картина течения слабо зависит от этого праметра и при x>5L (x - расстояние от лобовой точки вдоль траектории) поле течения выходит на некоторую асимптотику, существенно зависящую лишь от величин r1, g и
(rm,=b/2)
Пример такого расчёта дан на рис.6,а. Здесь изображено стационарное состояние баллистической ударной волны при обтекании гиперзвуковым потоком ( M > 10 ) тела, составленного из сферического затупления радиуса и примыкающего к нему цилиндра толщиной 2rm. Вдоль траектории указаны безразмерные давления `p=p/v2r1 за фронтом баллистической волны для случая rm=70 м, Qe=35°, когда передняя часть волны находиться на высоте 7 км над Землёй. Нестационарность процесса обтекания приближенно можно учитывать лишь меняя p1, r1 и скорость движения тела, которые определяются из тракторных расчётов (например типа представленных на рис. 2 ).
На рис 6,а схематически даны волны для четырёх последовательных моментов времени. В момент времени t отмечен приход волн к земной поверхности и их отражение как в окрестности конечной точки траектории, так и в её балистической части. Оказывается, что в плоскостях, перпендикулярных к движению тела (см. сечение S на рис.16,б ), течение газа аналогично таковому при взрыве шнурового заряда с удельной энергией E0. Это обстоятельство использовалось для приближения расчёта баллистических волн. Задавалось значение E0 в соответствии с (4.21) и затем по теории циллиндрического взрыва определялись параметры баллистических волн при их прохождении в атмосфере. Давления в лобо-вой точке тела за головной ударной волной могут быть вычислены по условиям на ударной волне и по законам сохранения для течения в окрестности критической точки. Оказывается, что давление в лобовой части тела. Параметры баллистических волн вдоль траектории можно расчитать с помощью ЭВМ для широкого набора значений E0(s) вдоль пути s по траектории. Процессы в конечной части траектории (момент t4 на рис. 6,б) моделировались расширением газового шара (раскалённые остатки тела плюс воздух) с давлением pm*. Полная энергия этого шара принималась равной E (объёмный сферический взрыв).