Борис Шустов - Астероидно-кометная опасность: вчера, сегодня, завтра

Избыточное и динамическое давления. Эффекты разрушения, вызванные взрывной волной, обычно соотносят с максимумом избыточного давления. В Нагасаки жилые дома разрушились на расстоянии вплоть до 2 км от эпицентра взрыва, где максимум избыточного давления Δp по оценкам составлял 20 кПа [Glasstone and Dolan, 1977]. Результаты ядерных испытаний показывают, что почти полное разрушение двухэтажного деревянного каркасного дома и неукрепленного кирпичного дома происходит при Δp = 30–35 кПа, прочный дом с каркасом из стали с алюминиевыми панелями обваливался при Δp = 20 кПа. Максимум избыточного давления, вызванного ударной волной, зависит не только от расстояния до эпицентра, но также от высоты взрыва. Величина Δp = 30 кПа достигалась на расстоянии 5 км от эпицентра для ядерного взрыва с энергией 1 Мт, когда он происходил на поверхности, и на расстоянии 9 км при взрыве на высоте 3,6 км [Glasstone and Dolan, 1977]. Для взрыва с энергией 30 Мт простые оценки, основанные на гидродинамическом подобии, дают радиус поражения R = 25 км при взрыве на высоте 10 км. Таким образом, размеры областей, которые могли бы быть опустошены падением метеороида с энергией около 30 Мт, сравнимы с радиусом крупного города. Существует другой важный фактор — динамическое давление, которое определяет повреждения, вызванные сильными ветрами, возникающими за фронтом взрывной волны. Для Δp = 35 кПа максимальное динамическое давление q = 4 кПа и максимальная скорость ветра составляет 260 км/ч. Но мы не будем углубляться в такие детали. Лишь упомянем, что полученная с помощью теоретических газодинамических расчетов [Коробейников и др., 1991] зона разрушения лесного массива неплохо согласуется с данными натурных исследований зоны воздействия при Тунгусском событии 1908 г.

Оценки величины зоны разрушений ударными волнами, вызванными падениями астероидов и комет, были даны в работе [Chapman and Morrison, 1994]. Авторы принимали во внимание избыточное давление на фронте ударной волны и динамическое давление, которое вызвало падение деревьев. Для калибровки они использовали Тунгусское событие [Зоткин, Цикулин, 1966], где ударные волны вызвали падение деревьев на площади 2000 км 2, которая эквивалентна площади круга с радиусом R s= 25 км. Предполагалось, что площадь A sсильных повреждений строений приблизительно равна площади вывала леса. Для оценок может быть использовано следующее выражение:

A s= πR 2s = 200E 2/3k, Rs = 8E k 1/3,

где радиус R sвыражен в километрах, площадь A s — в квадратных километрах, энергия ударяющего тела E k — в мегатоннах. Коэффициенты в формуле изменяются с высотой взрыва. «Оптимальная» высота h [км] составляет 6,4E 1/3[Glasstone and Dolan, 1977; Hills and Goda, 1993]. Если высота взрыва уменьшается до нуля, радиус R sуменьшается приблизительно в 1,4 раза. Для E = 30 Мт получается радиус R s= 18 км.

Прорыв атмосферы. Законы подобия, которыми мы пользовались для оценок, основаны на теоретических исследованиях гидродинамической задачи о распространении ударной волны после точечного взрыва в однородной атмосфере и на исследованиях с помощью численного эксперимента распространения взрывных волн, образованных химическими взрывами с высокой энергией и ядерными взрывами с энергиями менее 10 Мт. Но для больших энергий радиус R sоказывается порядка характерной высоты атмосферы или превышает ее (для E = 10 3Мт мы получаем R s= 80 км), поэтому оценки следует уточнять численными расчетами.

Двумерные численные расчеты взрыва в неоднородной атмосфере и различные теоретические оценки прорыва атмосферы показывают, что ударная волна вследствие уменьшения плотности воздуха с высотой над поверхностью Земли движется вверх быстрее, чем в радиальном направлении. Численные расчеты [Jones and Sanford, 1977; Jones and Kodis, 1982] взрыва с энергией 500 Мт, произведенного на поверхности, показали, что динамическое давление превосходит порог вывала леса на расстояниях в 27,5 км вместо 45 км, как это следует из закона подобия. В действительности давление и скоростной напор снижаются еще больше за счет влияния следа, не учтенного в этих расчетах.

Результаты расчета для большого тела (диаметром 10 км) приведены на рис. 8.3. Картина распределения плотности и температуры в более поздние моменты времени показана на рис. 8.4. Хотя след здесь также присутствует, но выброс вверх в основном облегчен за счет быстрого падения плотности атмосферы с высотой.

8.1.3. Световой импульс и пожары.Падение Тунгусского космического тела 30 июня 1908 г. вызвало пожар на площади около 500 км 2[Vasilyev, 1998], что в 4 раза меньше, чем площадь опустошения леса взрывными волнами (2000 км 2). Этот пожар наглядно демонстрирует роль светового излучения. В Хиросиме и Нагасаки 20–30 % всех жертв были ранены за счет ожога от прямого действия теплового излучения вспышки. Используя эти данные, можно ожидать, что лучистое воздействие при энергии взрыва в 10–30 Мт могло бы быть причиной ожогов незащищенной кожи первой степени (обратимое повреждение) для 82 % населения, а 15 % получили бы ожоги второй степени (которые можно вылечить за одну или две недели) [Glasstone and Dolan, 1977]. Конечно, жертвы прямого действия теплового излучения вне зоны «огненного шара» могут быть сокращены простыми способами гражданской обороны (убежище и другие меры защиты) при условии предупреждения о возможном падении космического объекта. Мы должны упомянуть также глазные травмы, вызывающие слепоту и ожоги сетчатки, но они также могут быть сокращены адекватной тренировкой использования специальных фильтров для глаз, опять же при условии предупреждения об опасности.

Рис. 8.3. Распределение относительной плотности ( а ) и температуры ( б ) непосредственно перед столкновением каменного тела диаметром 10 км с поверхностью Земли

Рис. 8.4. Распределение плотности (верхняя панель) и температуры (нижняя панель) через 5 с ( а ) и 15 с ( б ) после удара о Землю каменного тела диаметром 10 км со скоростью 20 км/с

Используя «экспериментальное» значение площади пожара для энергии E = 30 Мт (Тунгусское событие), получим следующие соотношения:

A r= 30E r, R f= 3E r 1/2 N = 2,9 10 3E r,

где A r — площадь воздействия теплового излучения в км 2, R f — радиус зоны пожара в км, N — число жертв, E r — энергия теплового излучения в Мт, которая для больших взрывов примерно составляет 30–40 % от энергии взрыва E [Glasstone and Dolan, 1977] или энергии ударника E k.

При ударе метеороида и его фрагментов о поверхность твердого тела со скоростью более 15 км/с происходит испарение ударника и поверхности мишени. При расширении в вакуум характерный размер излучающей области R fсоставляет примерно 10–15 размеров ударника [Melosh et al., 1993; Немчинов и др., 1998]. Время излучения порядка R f/V, где V — скорость удара о Землю. Излучательная эффективность в вакууме η зависит от скорости ударника, но в среднем очень мала (η ≈ 10 -4–10 -2). Такие значения применимы при ударах об астероиды, планеты и их спутники, лишенные атмосферы, например о Луну [Ortiz et al., 2000]. Именно эти значения приняты в программе СММ для излучения факела. Однако в условиях весьма плотной атмосферы Земли они сильно преуменьшают реальную излучательную эффективность. Во-первых, наличие атмосферы сильно сдерживает разлет паров и скорость их охлаждения, увеличивает их плотность и оптическую толщину. Во-вторых, воздух нагревается в ударной волне, генерируемой при расширении паров, и сам излучает.