Страница: 6/9
Приведём пример численного решения задачи входа в атмос-феру Тунгусского космического тела, выполненную конечноразностным методом Л.В.Соколовской.
Газообразное тело в форме циллиндра, высота которого равна диаметру 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 и
|
(4.21) |
(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 (объёмный сферический взрыв).
Реферат опубликован: 11/10/2009