Систему уравнений так называемой физической теории метеоритов (Kn<0.1) запишем в предположении движения тела в плоскости, проходящей через ось z:

Здесь

f - коэффициент реактивной отдачи, -1<f<1;

CL - коэффициент подъёмной силы,

i* - эффективная энтальпия разрушения

(характерная теплота сублимации или парообразования),

CH=CH(r,v,r) - коэффициент теплопередачи;

остальные обозначения такие же, как и в системе (4.1)-(4.3). Реактивной силой в уравнении (4.7) можно принебречь, если i*>1000 кал/г. Площадь А в общем случае - величина переменная, ибо масса тела меняется, причём для для случая шара:

Уравнению (4.10) можно придать следующий физический смысл: изменение массы,- Dm, за время пропорциональное кинетической энергии газа в объеме ,”охваченном” за это время миделем вдоль траектории (DtvA), и обратно пропорциональное энергии разрушения, то есть

Приведём теперь численные значения констант. Для высоты H=7.16 км; r0=1.29*10-3 г/см; g=9.8 м/с; R3=6371.7 км. Коэффициенты CD и CH зависят от v,r,r и находятся специальными расчётами, однако коэффициент CD можно приближённо считать равным 0.9; CH как функция v,r,r приводится в руководствах по метеоритике и аэродинамике. Коэффициент теплообмена обычно состоит из двух частей:

-конвективного теплообмена.

-радиационного.

Для крупных тел главную роль будет играть радиационный теплообмен. Для тел размером около 0.5 м при скорости входа ve=20 км/с и массе me=200 кг оценки показывают, что

0.01<CH<0.1; v>1 км/c

Коэффициент подъёмной силы CL, как правило, мал, и его обычно не учитывают в приближённых теориях, т.к. силы, действующие поперёк траектории ,малы. Эти силы могут возникать из-за неоднородности среды, реактивного эффекта, сильного ветра, угла наклона тела к направлению движения (угла атаки). На рис. 1 дан график изменения скорости движения тела в зависимости от высоты для фрагмента каменно-железного метеорита Лост-Сити, полёт которого был зафиксирован оптической камерой сети наблюдений. Найденная часть метеорита имела массу 15 кг, его скорость входа была ve=14.2 км/с, плотность rm=3.6 г/см3, i*=1300 кал/г, qe=43° (рис.1). Кружки на графике соответствуют данным наблюлений до скорости 3 км/с, когда метеорит перестал светиться. Потеря массы составила около 3 кг. Видно, что представленная модель для такого случая вполне удовлетворительна. Здесь же на рис.1 дана зависимость z(v) для случая ve=14.2 км/с, me=490 кг, rm=3.6 г/см3, i*=500 кал/г, qe=43° (штриховая линия). Видно, что траектории отличаются не так уж сильно, хотя абляция должна должна быть весьма интенсивной. Здесь могут быть и такие случаи случаи, когда практически вся масса метеороида испарится и снесётся в спутныйё поток, то есть (Dm/me)»1.

Американский астрофизик Д.О.Ре-Вилл выполнил расчёты для системы (4.7)-(4.10) при CL=0, ve=30 км/с, qe=45°, rm=3.7 г/см3, i*=2000 кал/г, me=10000 кг. Оказалось, что Dm»me на высотах, где v=3.5 км/с.

Таким образом, практически всё вещество распылилось в виде пара и мелких частичек в следе метеороида. Космическое тело “сгорело” до касания поверхности Земли. Здесь механизм испарения обусловлен сильными лучистыми потоками к поверхности воздуха, прошедшего через баллистическую волну при высоких скоростях до (до 5 км/с.)

Различные исследователи проводили опыты по деформации и разрушению водяных капель в потоках воздуха. По Дж. Ханту (Англия), при временах порядка tb происхрдит струйное “пробивание” в центре эллипсоидального тела и образование объёма в форме тора, который уже потом разрушается на более мелкие капли. Расчёты показали существенную роль процессов абляции и изменения формы при взаимодействии метеорита с атмосферой.

Так как влетающий в атмосферу метеороид холодный (температура его внутренних частей ниже температуры окружающей Среды), то можно считать, что энергия состоит только из кинетической. Углерод метеорита может гореть в атмосфере при сответствующих условиях. Но недостаток кислорода не позволит сгореть большому количеству углерода, и выделившаяся энергия не будет превышать кинетической энергии тела. Как же расходуется кинетическая энергия тела. Пусть тело затормозилось от скорости ve до скорости vc на пути. Это означает, что на этом пути энергия передалась окружающей среде за время tc. Время tc около 1-10 с, S порядка 80-40 км. Отсюда получаем, что с точки зрения воздействия на атмосферу метеориты подобны молнии: за малое время вдоль траектории выделяется энергия DE, на единицу длины приходится DE/S. Рассмотрим пример. Для метеорита типа Лост-Сити me=18 кг, mc=15 кг, ve=14 км/с, vc=3 км/с

DE»meve2/2

s=50 км, E0=DE/S=360 дж/см. метеороид подобен весьма длинной молнии с удельной энергией E0= meve2/2S. Для “сгорающих” метеороидов есть и внешнее сходство: они сверкают в небесной выси, как молнии. Очевидцы, наблюдавшие падение метеороидов, слышали и раскаты грома; баллистическая волна распространялась в атмосфере, подобно грозовой ударной волне.

Сформулированная выше упрощённая модель движения метеороида, объединённая с теорией линейных взрывов (грозовых разрядов), даёт возможность создать модель движения и взаимодействия метеороидов с атмосферой.

В заключении этого раздела коснёмся вопроса о характере и многообразии траекторий метеороидов. Не будем учитывать изменения массы, т.е. положим dm/dt=0, но CL¹0; отношение (CL/CD)=k называется аэродинамическим качеством движущегося тела.

Будем считать, что |k|£1, причём отрицательные значения k соответствуют наличию поперечной силы, действующей на тело “вниз” - в отрицательном направлении оси y местной системы координат, где ось x направлена вдоль вектора скорости, а ось y к ней перпендекулярна. Характерную величину m/CDA обозначим через b. За величину b примем значение 1515 кг/м2, что будет соответствовать входу в атмосферу сферического тела радиуса rE=97.8 м и плотностью rm=0.03 г/см. Обозначим через S расстояние вдоль поверхности Земли от проекции условной точки входа в атмосферу на эту поверхность. Пусть угол входа равен 20°, ZE=60км, vE=30 км/с. Меняя значения k, мы получим разные траектории и скорости тела при значениях аэродинамического качества k=0.5;-0.125;0;0.125;0.5 (S -расстояние от поверхности Земли) (рис. 2). При k=0.5 наблюдается явление рикошета .