Движение твердого метеороида в атмосфере - Метеоритное вещество и метеориты
Как уже говорилось выше, всю область полета метеороида можно разбить на две зоны. Первая зона будет соответствовать большим числам Кнудсена KN 0.1 ,а вторая зона - малым числам Кнудсена KN < 0.1. Эффектами вращения тела принебрегаем, форму его будем считать сферической с радиусом r. Будем предполагать тело однородным.
Сначала построим модель для первой зоны. В этой зоне изменением массы метеороида можно приберечь, так как абляции и разрушения тела практически нет. Уравнения движения следуют из законов ньютоновской механики:
Здесь
M - масса метеороида,
V - скорость,
- угол наклона вектора скорости к поверхности Земли,
G - ускорение силы тяжести,
- плотность атмосферы в точке,
A=rE2 - площадь поперечного сечения метеороида (площадь миделя),
Z - высота, отсчитываемая от уровня моря,
T - время,
CD - коэффициент сопротивления воздуха,
R3 - радиус Земли.
Изменение плотности воздух с высотой будем находить по барометрической формуле:
Где - плотность на уровне моря. Коэффициент CD Можно считать зависящим от числа Кнудсена, причем он убывает с высотой и меняется в пределах 2>CD>0.92 при изменении KN от 10 до 0.1.
Систему (4.1)-(4.3) нужно решать в предположении, что начальный момент времени при t=0 заданы zE=z, E=, vE=v, mE=m, то есть параметры входа метероида. За координату z, можно принять ту высоту, где согласно (4.1) сила тяготения Земли выравнивается с сопротивлением, то есть когда уравнение (4.5) при заданных m=mE, v=vE, можно считать за определение. Пренебрежем также изменением угла, то есть примем E= (это не внесет погрешностей, ибо есть малая величина для диапазона скоростей от 11 до < 70 км/с
(< 0.001 c-1).
После интегрирования уравнения (4.1) при условии пренебрежения силой mg sin и для z<zE получаем
Где B - баллистический коэффициент.
Приближенную формулу (4.6) можно использовать для оценки поведения решения при больших v. Видно, что vvE При zH. Это означает, что скорость тела практически не меняется.
Используя несложную компьютерную технику, систему (4.1)-(4.3) можно проинтегрировать с помощью любого подходящего численного метода, например метода Эйдлера с пересчетом. Сущность этого метода состоит в том, что для уравнения y'=f(x, y) сначала мы находим значение y'1=f(x0,y0) x+y0 где x0, y0 - начальная точка, а x - шаг интегрирования, затем берем
И находим уточненное значение y1=y'x+y0+O(x2)
Аналогичная процедура используется в случае системы уравнений.
Этот метод весьма прост для реализации даже с помощью программируемых микрокалькуляторов (вследствие простоты правых частей системы (4.1)-(4.3)).
Для расчета движения метеорита в нижних слоях атмосферы система (4.1)-(4.4)не годится, т. к. она не учитывает абляцию (изменение массы), поэтому перейдем к описанию более сложной модели, пригодной для низких высот, т.е. для второй зоны.
Систему уравнений так называемой физической теории метеоритов (KN<0.1) запишем в предположении движения тела в плоскости, проходящей через ось z:
Здесь
F - коэффициент реактивной отдачи, -1<f<1;
CL - коэффициент подъемной силы,
I* - эффективная энтальпия разрушения
(характерная теплота сублимации или парообразования),
CH=CH(r, v,) - коэффициент теплопередачи;
Остальные обозначения такие же, как и в системе (4.1)-(4.3). Реактивной силой в уравнении (4.7) можно принебречь, если I*>1000 кал/г. Площадь А в общем случае - величина переменная, ибо масса тела меняется, причем для для случая шара:
Уравнению (4.10) можно придать следующий физический смысл: изменение массы,- m, за время пропорциональное кинетической энергии газа в объеме,"охваченном" за это время миделем вдоль траектории (tvA), и обратно пропорциональное энергии разрушения, то есть
Приведем теперь численные значения констант. Для высоты H=7.16 км; 0=1.29*10-3 г/см; g=9.8 м/с; R3=6371.7 км. Коэффициенты CD и CH зависят от v,,r и находятся специальными расчетами, однако коэффициент CD можно приближенно считать равным 0.9; CH как функция v,,r приводится в руководствах по метеоритике и аэродинамике. Коэффициент теплообмена обычно состоит из двух частей:
- -конвективного теплообмена. -радиационного.
Для крупных тел главную роль будет играть радиационный теплообмен. Для тел размером около 0.5 м при скорости входа vE=20 км/с и массе mE=200 кг оценки показывают, что
0.01 CH0.1; v1 км/c
Коэффициент подъемной силы CL, как правило, мал, и его обычно не учитывают в приближенных теориях, т. к. силы, действующие поперек траектории, малы. Эти силы могут возникать из-за неоднородности среды, реактивного эффекта, сильного ветра, угла наклона тела к направлению движения (угла атаки). На рис. 1 дан график изменения скорости движения тела в зависимости от высоты для фрагмента каменно-железного метеорита Лост-Сити, полет которого был зафиксирован оптической камерой сети наблюдений. Найденная часть метеорита имела массу 15 кг, его скорость входа была vE=14.2 км/с, плотность M=3.6 г/см3, I*=1300 кал/г, E=43 (рис.1). Кружки на графике соответствуют данным наблюдений до скорости 3 км/с, когда метеорит перестал светиться. Потеря массы составила около 3 кг. Видно, что представленная модель для такого случая вполне удовлетворительна. Здесь же на рис.1 дана зависимость z(v) для случая vE=14.2 км/с, mE=490 кг, M=3.6 г/см3, I*=500 кал/г, E=43 (штриховая линия). Видно, что траектории отличаются не так уж сильно, хотя абляция должна быть весьма интенсивной. Здесь могут быть и такие случаи, когда практически вся масса метеороида испарится и снесется в спутный поток, то есть (m/mE)1.
Американский астрофизик Д. О.Ре-Вилл выполнил расчеты для системы (4.7)-(4.10) при CL=0, vE=30 км/с, E=45, M=3.7 г/см3, I*=2000 кал/г, mE=10000 кг. Оказалось, что mmE на высотах, где v=3.5 км/с.
Таким образом, практически все вещество распылилось в виде пара и мелких частичек в следе метеороида. Космическое тело "сгорело" до касания поверхности Земли.
Здесь механизм испарения обусловлен сильными лучистыми потоками к поверхности воздуха, прошедшего через баллистическую волну при высоких скоростях до (до 5 км/с.). Различные исследователи проводили опыты по деформации и разрушению водяных капель в потоках воздуха.
По Дж. Ханту (Англия), при временах порядка tB происходит струйное "пробивание" в центре эллипсоидального тела и образование объема в форме тора, который уже потом разрушается на более мелкие капли.
Расчеты показали существенную роль процессов абляции и изменения формы при взаимодействии метеорита с атмосферой.
Так как влетающий в атмосферу метеороид холодный (температура его внутренних частей ниже температуры окружающей Среды), то можно считать, что энергия состоит только из кинетической. Углерод метеорита может гореть в атмосфере при соответствующих условиях.
Но недостаток кислорода не позволит сгореть большому количеству углерода, и выделившаяся энергия не будет превышать кинетической энергии тела. Как же расходуется кинетическая энергия тела.
Пусть тело затормозилось от скорости vE до скорости vC на пути. Это означает, что на этом пути энергия передалась окружающей среде за время tC. Время tC около 1-10 с, S порядка 80-40 км.
Отсюда получаем, что с точки зрения воздействия на атмосферу метеориты подобны молнии: за малое время вдоль траектории выделяется энергия E, на единицу длины приходится E/S.
Рассмотрим пример. Для метеорита типа Лост-Сити mE=18 кг, mC=15 кг, vE=14 км/с, vC=3 км/с
EmEVE2/2
S=50 км, E0=E/S=360 дж/см. метеороид подобен весьма длинной молнии с удельной энергией E0= mEVE2/2S. Для "сгорающих" метеороидов есть и внешнее сходство: они сверкают в небесной выси, как молнии. Очевидцы, наблюдавшие падение метеороидов, слышали и раскаты грома; баллистическая волна распространялась в атмосфере, подобно грозовой ударной волне.
Сформулированная выше упрощенная модель движения метеороида, объединенная с теорией линейных взрывов (грозовых разрядов), дает возможность создать модель движения и взаимодействия метеороидов с атмосферой.
В заключении этого раздела коснемся вопроса о характере и многообразии траекторий метеороидов. Не будем учитывать изменения массы, т. е. положим dm/dt=0, но CL0; отношение (CL/CD)=k называется аэродинамическим качеством движущегося тела.
Будем считать, что |k|1, причем отрицательные значения k соответствуют наличию поперечной силы, действующей на тело "вниз" - в отрицательном направлении оси y местной системы координат, где ось x направлена вдоль вектора скорости, а ось y к ней перпендикулярна. Характерную величину m/CDA обозначим через. За величину примем значение 1515 кг/м2, что будет соответствовать входу в атмосферу сферического тела радиуса rE=97.8 м и плотностью M=0.03 г/см. Обозначим через S расстояние вдоль поверхности Земли от проекции условной точки входа в атмосферу на эту поверхность. Пусть угол входа равен 20, ZE=60км, vE=30 км/с. Меняя значения k, мы получим разные траектории и скорости тела при значениях аэродинамического качества k=0.5;-0.125;0;0.125;0.5 (S - расстояние от поверхности Земли) (рис. 2). При k=0.5 наблюдается явление рикошета.
При значениях k<0 траектории могут иметь вертикальную касательную, а при k<-1 пролетную - образную траекторию.
Из рис. 2 видно, что скорость тела остается практически постоянной до высоты 40 км.
Кроме описанных выше параметров вычисляется интенсивность свечения I по формуле
(4.11)
Где 0 - коэффициент эффективности свечения (опытный параметр).
Опишем вкратце более общую модель входа метеороида в атмосферу. Уравнения (4.7)-(4.10) описывают движение центра масс метеороида. Кроме этого следовало бы описать движение метеороида около центра масс. Довольно трудной задачей является определение параметров тела и окружающего воздуха, включая след за телом. Для этой задачи следует на определенных этапах (для дискретного набора времени t=tJ) проводить расчет обтекания и абляции, а так же механической деструкции тел, с учетом эффекта теплопередачи и излучения, а так же высвечивание метеороидов в различных спектральных диапазонах). Нужно рассчитывать распространение атмосферных возмущений в пространстве и времени. Следует изучить вопросы, связанные с моделированием воздействия удара метеороидов и балистических волн о поверхность Земли.
Похожие статьи
-
Тунгусское космическое тело - Математическое моделирование движения метеоритов в атмосфере
30 июня 1908 г. произошло столкновение с атмосферой Земли космического тела, нижняя часть траектории которого проходила над Вост. Сибирью. Траектория...
-
Метеоритное вещество и метеориты - Математическое моделирование движения метеоритов в атмосфере
Каменные и железные тела, упавшие на Землю из межпланетного пространства, называются метеоритами, а наука, их изучающая-метеоритикой. В околоземном...
-
Приведем наглядный пример с автомобилем, который то увеличивает, то уменьшает скорость. Это позволяет нам провести некоторую аналогию с изменением...
-
Сила сопротивления атмосферы - Движение искусственного спутника Земли в нецентральном поле тяготения
При изучении поступательного движения спутника принимают во внимание лишь ту компоненту аэродинамических сил, направление которой противоположно вектору...
-
Плотность атмосферы - Движение искусственного спутника Земли в нецентральном поле тяготения
Плотность воздуха верхней атмосферы определяется как непосредственно, при помощи аппаратуры, установленной на спутниках и ракетах, так и косвенным путем,...
-
2.1 Дифференциальные уравнения движения искусственного спутника Мы ввели подвижную, жестко связанную с Землей, систему координат и соответствующие ей...
-
Атмосферные явления и феномены, Движение атмосферы - Юпитер среди планет Солнечной системы
Движение атмосферы Анимация вращения Юпитера, созданная по фотографиям с "Вояджера-1", 1979 г. Скорость ветров на Юпитере может превышать 600 км/ч. В...
-
Межслоевые процессы, Движение атмосферы - Планеты Солнечной системы: Юпитер
Локализовать все процессы внутри независимых слоев невозможно: необходимо объяснять недостаток химических элементов в атмосфере, избыточное излучение и...
-
Введение В предыдущих главах были рассмотрены возмущения элементов орбиты спутника, вызываемые зональными и тессеральными гармониками геопотенциала и...
-
Введение Помимо несферичности Земли, притяжения Луны и Солнца, сопротивления атмосферы и светового давления на движение спутника действует целый ряд...
-
Влияние электромагнитных сил - Движение искусственного спутника Земли в нецентральном поле тяготения
Пусть спутник обладает электрическим зарядом, равным Q. Тогда при движении в магнитном поле Земли на него будет действовать сила F, определяемая формулой...
-
Метеориты - Малые тела Солнечной системы
В околоземном космическом пространстве движутся самые различные метеороиды (космические осколки больших астероидов и комет). Их скорости лежат в...
-
Астероиды и метеориты - Астероиды и метеориты
Довольно часто на небе появляются космические пришельцы. Их размеры исчисляются от нескольких сотен метров до тысячи километров. Это астероиды и кометы....
-
Параболическое движение, Прямолинейное движение - Математическое моделирование движения небесных тел
(p = 0, e = 1) Уравнение параболической орбиты записывают в видеp r = 1 + cos v (1.80) Где величина определяет расстояние от центра притяжения M0 до...
-
Покажем, что для любой классической системы, обладающей центральной симметрией и заданной энергией, существует такая метрика, что действие системы будет...
-
Как образуется кратер - Астероиды и метеориты и их влияние на поверхность Земли
Размер, радиус R Кратера, который образуется при сверхзвуковом столкновении метеорита с поверхностью, можно приближенно установить из подсчета того, на...
-
РЕАКТИВНОЕ ДВИЖЕНИЕ - К. Э. Циолковский, его вклад в развитие ракетостроения
"Практические дела делаются только исходя из общих начал, только при знакомстве с абстрактами, до них относящимися". (Д. И. Менделеев) Среди великих...
-
Нецентральность гравитационного поля Земли - Возмущенное движение космического аппарата
Возмущенный движение гравитационный орбита При решении ограниченной задачи двух тел Земля представляется шаром со сферическим распределением плотности. В...
-
Строение Солнечной системы - Законы движения небесных тел и строение Солнечной системы
Хорошо известно, что основная масса Солнечной системы (около 99.8%) приходится на ее единственную звезду - Солнце. Суммарная масса планет составляет...
-
Возмущенное движение - Возмущенное движение космического аппарата
Описание и изучение орбит КА на основе решения ограниченной задачи двух тел является лишь первым этапом при определении реальных движений тел любой...
-
Безопорное движение&;nbsp; - Двигатель космолета на эффекте гравитационного самоускорения
Такое "вытягивание" потенциала, зависящего от скорости удаления Солнца, наводит на интересную мысль. А что, если "получатель потенциала", скажем,...
-
Вращательное движение - Звездная аберрация против релятивистской астрономии
Посмотрим, как можно объяснить парадокс Эренфеста, если считать, что пространство является общим для всех инерциальных систем, а время для них едино....
-
ЗАКЛЮЧЕНИЕ - Три закона движения планет
Законы Кеплера замечательны и тем, что они, если можно так выразиться, более точны, чем сама действительность. Они представляют собой точные...
-
Метеориты - Астероиды, метеориты, кометы, мелкие осколки
В околоземном космическом пространстве движутся самые различные метеороиды (космические осколки больших астероидов и комет). Их скорости лежат в...
-
Круговое движение, Эллиптическое движение - Математическое моделирование движения небесных тел
(p = 0, e = 0) Круговое движение представляет собой наиболее простой случай движения в задаче двух тел. Только для кругового движения (и прямолинейного...
-
Введение - Математическое моделирование движения небесных тел
В небесной механике для описания движений небесных тел в зависимости от конкретных условий используются различные физические модели - идеализированные...
-
Движение тела, Первая космическая скорость - Космические скорости
Для определения значений космических скоростей сначала рассмотрим движение тела, брошенного на расстоянии h от поверхности Земли, с начальной скоростью v...
-
Сущность метеоритов и комет - Чрезвычайные ситуации космического характера
Метеорит -- твердое тело космического происхождения, упавшее на поверхность Земли. Большинство найденных метеоритов имеют вес от нескольких граммов до...
-
Уравнения движения МКА Рассмотрим невозмущенное движение материальных точек М и m в некоторой инерциальной системе координат. Движение совершается под...
-
Происхождение метеоритов, Состав метеоритов и их вещества - Метеориты
В настоящее время во многих музеях мира хранится не менее 500 тонн метеоритного вещества. Расчет показывает, что в виде метеоритов и метеорной пыли за...
-
Начало метеоритных исследований - Метеориты
Метеориты известны человеку уже многие тысячи лет. Обнаружены орудия первобытных людей, сделанные из метеоритного железа. Случайно находя метеориты, люди...
-
Тринадцать движений Земли - Планета Земля
Прежде чем подробно рассмотреть те движения нашей планеты, которые имеют непосредственное отношение к ее недрам, представим общую картину очень сложно...
-
Орбиту можно получить как линию пересечения двух поверхностей. Уравнение одной поверхности - это уравнение плоскости орбиты. Уравнение второй поверхности...
-
Гиперболическое движение - Математическое моделирование движения небесных тел
(p=0,e>1) Каноническое уравнение гиперболы в центральных прямоугольных координатах O?!?! представляется в виде (1.68) Где a - действительная, а b...
-
Астроклимат, это совокупность атмосферных условий, влияющих на качество астрономических наблюдений. Важнейшие из них - прозрачность воздуха, степень его...
-
Уравнение Бине - Математическое моделирование движения небесных тел
Другой способ получения траектории движения в задаче двух тел связан с широко известным уравнением Бине. Это уравнение записывается в цилиндрической...
-
Мелкие осколки - Астероиды, метеориты, кометы, мелкие осколки
Околосолнечное пространство заполнено мелкими частицами, источниками которых служат разрушающиеся ядра комет и столкновения тел, в основном, в поясе...
-
Общее решение уравнений относительного движения - Математическое моделирование движения небесных тел
Рассмотрим на небесной сфере сферический треугольник N M x, где M - проекция текущего положения тела M на небесную сферу. Сторонами этого треугольника...
-
Основные слои в атмосфере Солнца - Солнце и Солнечная система
Солнце - центральное тело Солнечной системы - представляет собой очень горячий плазменный шар. Солнце - ближайшая к земле звезда. Свет от него доходит до...
-
Метеориты - вестники космоса - Строение Солнечной системы
Метеориты - каменные или железные тела, падающие на Землю из межпланетного пространства. Падение метеоритов на Землю сопровождается звуковым, световым и...
Движение твердого метеороида в атмосфере - Метеоритное вещество и метеориты