Экспериментальное и численное исследование газодинамики разлета продуктов детонации конденсированных взрывчатых веществ с твердыми частицами тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

>Г6

- 5

ОД

ИЮН 19Я5

МОСКОВСКИЙ ГОСУДАРСТВЕННЫЙ УНИВЕРСИТЕТ им. М.В. ЛОМОНОСОВА

на правах рукописи

Теверовский Михаил Анатольевич

ЭКСПЕРИМЕНТАЛЬНОЕ И ЧИСЛЕННОЕ ИССЛЕДОВАНИЕ ГАЗОДИНАМИКИ РАЗЛЕТА ПРОДУКТОВ ДЕТОНАЦИИ КОНДЕНСИРОВАННЫХ ВЗРЫВЧАТЫХ ВЕЩЕСТВ С ТВЕРДЫНИ ЧАСТИЦАМИ

Специальность: 01.02.05 -механика жидкости, газа и плазмы

Автореферат диссертации на соискание ученой степени кандидата физико-математических наук

Москва 1995

Р.чбота выполнена в Институте механики МГУ им. М.В. Ломоносова.

Научный руководитель: доктор физико-математических наук, профессор А.Л. Гонор

Официальные оппоненты: доктор физико-математических наук, профессор H.H. Смирнов кандидат физико-математических наук, старший научный сотрудник В.А. Куликовский

Ведущая организация: Институт химической физики РАН

30

г. в \ ' часов

Защита состоится " 23 " МлХНиЛ 1995 г. в В

на заседании Диссертационного совета Д.053.05.02

при МГУ им. 11.13. Ломоносова в ауд. Адрес: 119899, Москва, Ленинские горы, главное здание МГУ

С диссертацией можно ознакомиться в библиотеке механико-математического факультета МГУ

Автореферат разослан \М <U£ 19Э5 г.

Ученый секретарь Диссертационного совета

профессор В.П. Карликов

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Среди разнообразных взрывчатых веществ особое место занимают гетерогенные смеси последних с дисперсными, иедетопирующими добавками в виде твердых частиц. Введением частиц удается изменять энергетические характеристики заряда, работоспособность расширяющихся продуктов детонации, параметры детонационной волны. Взрывчатые смеси с дисперсными частицами широко используются в технике и встречаются в природе: в пиротехнических устройствах, при детонационном напылении материалов, при горении металлизированных топлив ракетных двигателей, при смешении угольный пыли с рудничным газом в шахтах н т.д.

При распространении детонации в таких смесях за фронтом волны образуется газовзвесь продуктов детонации с частицами, разлет которой сопровождается образованием мощных ударных волн в окружающей среде и, в случае химически активных частиц, их воспламенением и горением. В настоящее время подробно изучены влияние частиц на параметры детонационной волны и газодинамика одномерного движения газовзвеси продуктов детонации с твердыми частицами. Характер распространения детонационной волны в цилиндрическом заряде взрывчатого вещества с твердыми частицами (гладкость фронта, стационарность) и газодинамические особенности двумерного разлета газовзвеси продуктов детонации изучены недостаточно.

Таким образом, исследование детонации и разлета продуктов взрывчатых веществ с твердыми частицами, проведенное в диссертационной работе, является актуальным в теоретических и практических приложениях.

Цель работы: Экспериментальное исследование процесса разлета цилиндрических образцов, заполненных взрывчатым веществом п металлическими частицами. Построение математической модели и численное исследование газодинамики двумерного разлета продуктов детонации взрывчатого вещества с твердыми, инертными частицами за волной, распространяющейся вдоль оси цилиндрического заряда.

Научная новизна. В диссертации экспериментально исследован характер распространения детонационной волны в цилиндрическом заряде металлизированного взрывчатого вещества. Получены немонотонные зависимости скорости детонации и угла разлета продуктов в окрестности фронта волны в зависимости от массовой концентрации частиц. Экспериментально исследовано влияние горения частиц на растекание продуктов детонации прореагировавшего заряда в объеме взрывной камеры.

На основе результатов экспериментов построена математическая модель осескмметричиого разлета газовзвеси продуктов детонации с твердыми, инертными частицами за плоскостью Чепмена-Жуге в постановках "запыленного" и эффективного газов. Численным моделированием выявлены основные газодинамические особенности разлета, обусловленные влиянием частиц и бокового растекания газовзвесп.

Выполнен расчет параметров сжатого ударной волной слоя воздуха в окрестности фронта детонационной волны с учетом диссоциации за скачком. Показано, что ятя плотных взрывчатых веществ температура сжатого ударно!) волной слоя воздуха оказывается выше температуры продуктов детонации в плоскости Чепмена-Жуге.

Практическая ценность работы. Полученные в диссертации результаты могут быть использованы при разработке технологий с применением детонации зарядов взрывчатых веществ с твердыми частицами, при проектировании устройств, создающих высоконерасчетные двухфазные струн. Проведенный анализ температур в сжатой слое воздуха и плоскости Чепмена-Жуге позволяет по-новому интерпретировать результаты измерения температур при детонации конденсированных взрывчатых веществ.

Апробация работы. Результаты диссертационной работы докладывались и обсуждались ка семинарах в Институте механики (рук. академик Г.Г. Черный) и кафедры "Волковой и газовой динамики" механико-математического факультета (рук. академик В.И. Шемякин) Московского государственного университета им. М.В, Ломоносова. Результаты работы получили положительную оценку.

Структура и объем работы. Диссертация состоит из введения, четырех глав и заключения. Диссертация содержит 162 страницы машинописного текста, включая 33 рисунка н список цитированной литературы из 136 наименований.

СОДЕРЖАНИЕ РАБОТЫ

Во введении отмечена актуальность темы исследования, кратко приведены основные физические эффекты и теоретические результаты в исследованиях детонации конденсированных взрывчатых веществ (ВВ) с твердыми частицами. Сформулированы цели работы н кратко излагается ее содержание.

В первой гладе даны общая классификация взрывчатых смесей с частицами, обзор теоретических и экспериментальных

нсследовашш по детонации ВВ с твердыми частицами.

В § 1.1 рассмотрены две основные группы взрывчатых смесей с частицами:

A. Газовые гетерогенные взрывчатые смеси: взвеси твердых или жидких частиц в газообразной взрывчатой или горючей смеси.

B. Конденсированные ВВ с твердыми частицами.

Внутри каждой группы взрывчатые смеси классифицируются по характеру взаимодействия фаз при детонации.

В § 1.2 описаны газодинамические модели ("запыленного" газа, детонационных волн), используемые в численном моделировании движения газовзвесей с частицами за детонационными волнами. С помощью характерных ударных и детонационных адиабат рассмотрены различные состояния фаз при детонации, структура гетерогенных детонационных волн.

В § 1.3 дан обзор результатов основных экспериментов по исследованию влияния твердых частиц (химически активных пли инертных) на разлет продуктов и параметры детонации конденсн-р о ванных ВВ.

Вторая глава посвящена экспериментальному исследованию процесса разлета образцов, заполненных конденсированным ВВ и мелкодисперсными алюминиевыми частицами. Цель эксперимента — выявление характерных особенностей процесса необходимых для построения математической модели разлета продуктов детонации ВВ с твердыми частицами.

В § 2.1 приведена постановка и методика эксперимента. Экспериментальная установка на базе воздушной взрывной камеры включала в себя высокоскоростную фотокамеру ВСК-5, систему автоматизированного измерения давления в воздухе датчиками ЛИД и дистанционный пульт подрыва образца с детопп-

рующей смесью. Фотокамера применялась для визуализации процесса детонации образца. Высокая частота покадровой сьемкл (282-103 кадров/с) позволила получить подробные кинограммы распространения детонационной волны в образце.

Датчики давления регистрировали разлет продуктов продетонировавшего образца: падающую и отраженную от стенок взрывной камеры ударные волки, растекание облака продуктов детонации. Запуск датчиков происходил при падении ударной волны взрыва на чувствительный элемент. Запись информации с датчиков в запоминающее устройство производилась последова-. тельно. Время между двумя измерениями составляло 31,25 мке; общее время работы датчика от момента запуска - 8 мс.

Схема эксперимента показана на рис. 1. Образец 1 вывешивался в центре взрывной камеры 2 в поле объектива фотокамеры 3. Сигнал на подрыв образца поступал с выносного дистанционного пульта. Оптическая ось объектива ВСК-5 совмещалась с осью смотровых окон взрывной камеры, снабженных обоймами с оптическими стеклами 5. Фотокамера находилась в ждущем режиме. Сигналом для запуска фотокамеры служил им пульс от фотоэлемента 4, реагирующего на вспышку при подрыве детонатора. Т.к. процесс разлета образца сопровождался ярким свечением продуктов детонации, то дополнительной подсветки не требовалось. Датчики давления б помещались в ящик с песком 7 у стенки камеры на расстоянии примерно 1,7 м от образца. Поверхность чувствительного элемента совпадала с уровнем песка в ящике. Блок памяти датчиков располагался в камере и защищался от повреждений специальной крышкой (на рисунке не показан).

Образец представлял собой тонкостенный алюминиевый цилиндр, заполненный взрывчатым веществом (взрывчатый иапол-

Рис. I Схема экспериментальной установки. I - модель, 2 - взрывная камера, 3 - фотокамера ВСК-5, 4 - фотоэлемент, 5 - оптические диски, б - датчики измерения давления АВД, 7 - ящик с песком. .

нитель) и алюминиевыми частицами. В экспериментах использовались образцы с различной массовой концентрацией а5 частиц алюминия в наполннтеле.

В § 2.2 рассматривается процесс развития и распространения детонационной волны в металлизированном ВВ, представленный на кинограммах разлета. Продукты детонации видны б виде ярко светящейся области. Фронт детонационной волны идентифицировался с передней поверхностью светящейся зоны в образце. Скорость детонационной волны определялась по кинограмме измерением перемещения фронта. Окончательный результат получался осреднением по данным нескольких опытов.

Анализ и обработка кинограмм разлета показали, что детонационная волна в металлизированной смеси имеет негладкий фронт, продукты детонации за волной занимают симметричную конусообразную область, процесс распространения детонации в образце является квазистацнонарным.

Получена немонотонная зависимость скорости детонации О от массовой концентрации <х металлических частиц (рис. 2, где пунктирными линиями показан доверительный интервал с учетом 10%—ной ошибки измерения). Явление уменьшения скорости детонации при добавлении в заряд ВВ металлических частиц хорошо известно из экспериментов. Полученный результат — немонотонное изменение скорости детонации при увеличении содержания металлических частиц — раннее не отмечался. *

Для объяснения немонотонного поведения скорости детонации металлизированных ВВ предложен механизм, связанный с одновременным действием двух факторов: искривлением фронта детонационной волны в окрестности частицы и горением алюминиевых частиц за детонационной волной. Уменьшение скорости

Б*10

ч N \ \ ч X У / ' /

Ч \ Ч_ Ч ч ч. ч. ч "»ч У /

0.1

0.2

0.3

0.4

0.5

Рис. 2. Зависимость скорости детонации от массовой концентрации частиц.

Погрешность измерения - 760 м/с.

80

Г(Ч9 70

60

50

40

30

, " ~ ^ ___ ч

/ у / / ' / ' у / ч ч

у / у » И"'' / У \\ \

у

«

0.1 0.2 0.3 0.4 0.5

а.

Рис. 3. Зависимость угла разлета продуктов детонации от массовой концентрации частиц.

Погрешность измерения — 4°'

а

детонации, по-видимому, обусловлено искривлением фронта волны, отмеченное повышение — пороговым характером эффекта горения алюминиевых частиц: при достаточно интенсивном тепловыделении формируется импульс сжатия, который усиливает детонационную волну.

Квазистационарность процесса детонации образцов подтверждают и результаты измерения угла разлета у продуктов детонации в окрестное™ фронта волны. С удовлетворительной точностью этот угол остается постоянным все время движения детонационной волны по образцу. Зависимость угла разлета от массовой концентрация частиц приведена на рис. 3. Как и в случае Дс^) полученные точки зависимости образуют немонотонную кривую.

§ 2.3 посвящен анализу движения продуктов продетониро-вавшнх образцов в объеме взрывной камеры. По результатам измерения давления (рис. 4) сравниваются картины разлета продуктов детонации алюмтшзированного (рис. 4 а, б, в) и "обычного" (без частиц) зарядов ВВ (рис. 4 г).

Ударно—волновая стадия разлета продуктов детонации в камере отмечена характерными всплесками, которые соответствуют падающей (точки 1 на рнс. 4) и отраженной от стенки камеры (точки 2) ударным волнам. Отличительной особенностью детонации металлизированного ВВ является образование газовзвеси продуктов детонации с частицами, в которой протекает реакция горения с интенсивным выделением тепла. Газовзвесь движется за ударной волной в виде высокотемпературного облака с высоким внутренним давлением. Расширяясь, облако создает длительный импульс давления в камере, который фиксируется датчиком (точки 3 на рис. 4 а-в). Отсутствие на графике резких всплесков

6 t, мс 8

кг/см'г

л?< л

кг/см j

8 t ,МС

Рис Ь Зависимость избыточного давления в камере от

а. 1с/ = 0,14, б. -*s=0,23, в. -ot{» 0,44, г. - <*Е

времени,

= О

6 t, мс

указывает на беэволнопой характер этой стадии. Протяженная во времени, по сравнению с взрывом обычных ВВ, фаза сжатия является следствием горения газовзвсск. Детонация ВВ без частиц происходит без образования реагирующей взвеси; по завершении ударно—волновой стадии процесс быстро затухает (рис. 4 г).

Глава 3 поевгацена построению математической модели разлета продуктов детонации (ПД) конденсированного ВВ с твердыми инертными частицами за плоскостью Чепмена-Жуге (ЧЖ) при распространении стационарной детонационной волны (ДВ) вдоль оси цилиндрического заряда. В основу модели положены результаты экспериментальных исследований распространения детонации в металлизированном ВВ.

В § 3.1 дана постановка задачи о разлете. Рассматривается полубесконечный цилиндр радиуса Г!., заполненный ВВ и твердыми частицами. С торца инициируется ДВ. Пренебрегая участком формирования, детонация считается стационарной, характеризуемой постоянной скоростью волны О и описываемой моделью ЧЖ. За плоскостью ЧЖ имеется монодисперсная смесь ПД (несущая фаза) с частицами (дисперсная фаза), движение которой определяется в рамках модели "запыленного" газа. Объемной долей частиц прене-брегается. Частицы химически инертны.

Осесимметричная схема разлета продуктов при детонации ЧЖ цилиндрического заряда ВВ в воздухе в системе координат, связанной с волной, показана иа рис. 5. Возможность применения такой схемы подтверждается расчетом.

Задача о разлете газовзвесн ПД рассматривается кап задача об истечении стационарной, осесимметрмчноп. неравновесной по скорости и температуре двухфазной струп с иерасчет-ностью " = Р]/Рх в сверхзвуковой спутный воздушный поток с

У в

^—

13

--"""

А /¿С — — = = - __ Б'

м,

И 1

Рис. 5. Схема разлета продуктов детонации цилиндрического заряда ВВ с твердыми частицами.

ОА— детонационная волна; АВ~ ударная волна в воздухе; АС- контактная поверхность; АС— сепаратриса потока частиц; БЗ'- висячий скачок.

0.4

0.3

0.2

0.1

0.0 0.01

0

\

\

0.1

Рис. 6. Давление на оси струи Р„(х; /?), а,= 0,34.

0- /3= 0 14- 0= 24 5-

0.2 2- /3= 1 3- /3= 4 100 6- р= 400 7- 0= „

числом Маха М^ = D/a^ (здесь pj — давление и плоскости ЧЖ; рв, —давление, скорость звука в набегающем спутном потоке воздуха, соответственно). Проведенные оценки положения диска Маха показали, что на расстояниях около трех калибров от плоскости ЧЖ течение газа является чисто сверхзвз'ковым.

Стационарное движение неравновесной по скорости и температуре газовзвеси описывается системой уравнений "запыленного" газа (Marble F.E., Ann. Rev. Fluid Mech. V 2. 1970):

где |0£ — "размазанная" плотность среды частиц; и>— пол-

пая энтальпия; /г д. — компоненты силы /8 и тепловой поток в межфазном взаимодействии, соответственно; ге—удельная теплоемкость материала частицы; Я —газовая постоянная струи; остальные обозначения общепринятые. Уравнения написаны в безразмерном виде. В качестве основных масштабов взяты значения параметров в плоскости ЧЖ, давление отнесено к величине

Здесь н ниже индексом отмечены значения параметров в плоскости Ч/К: индексом "о" — начальные параметры дисперсной смеси; индексом "э" - параметры среды частиц; индексом "ос" -

Iрпю) + |Jpvw) + + Ps(vJx+vJy) + = 0

* »2

парамстры спутного потока; "*" - размерные величины.

Система замыкается заданием силы /4, теплового потока gs, действующих на одиночную частицу со стороны несущей фазы и уравнений состояния фаз.

Сила / задается, как сила аэродинамического сопротивления с поправками на инерционность и сжимаемость несущей фазы (Carlson D.E. fc Hoglund R.F., AJAA. V6. N2. 1964). Лоток тепла <¡s определяется конвективным теплообменом частицы с газом с поправкой на инерционность несущей фазы (Drake R.M., J. Heat Transfer. 1961. V. 83):

O co/cs

G = (1+ 0,15- Äe°>687)-( 1 + exp(-^Z---?_)],

s Re¡

D = 1 + 0,23 • .ñejj'55 • iV'33,

fíes = ñe,.Mä = M¡-Iblsi; Re. =

rr n

Т^(и>0), Pr = 47/(97-5)

ai

где G, D — поправки, учитывающие реальные свойства

несущей фазы; Äes, Ms - числа Рейнольдса н Маха относительного

* *

обтекания частицы, соответственно; /?= R /iv;

* * * * *

/v = ms-ii j/(6jTos-/ij) — характерная длина релаксации скорости * * *

частицы; ms-масса частицы; —число Прандтля;

* *

ср — удельная теплоемкость несущей фазы; <rs — радиус частицы; * *

Л , /i — коэффициенты теплопроводности и вязкости газа, соот-*

ветственно; а — скорость звука в газе. Параметр ß описывает интенсивность релаксационных процессов в газовзвеси: при /?-> О частицы сохраняют значения параметров, заданные граничными

условиями ("замороженное" течение), прн /?->ш смесь описывается моделью равновесного (эффективного) газа.

Уравнения состояния фаз. Уравнение состояния дисперсной фазы, в силу условия несжимаемости, имеет вид: = const (здесь р°~ плотность материала частицы). Для ПД (? = 3) и воздуха (7 = 1,4) принимаются уравнения состояния идеального совершенного газа:

ж ♦ А ^ ф

где я = р, /рв, г = lp9, t = Tj /Г..

Граничные условия. Границами расчетной области в задаче о разлете газовзвеси ПД с твердыми частицами являются (рис. 5): детонационная волна (левая граница), ударная волна АВ (внешняя граница), контактная поверхность АС (внутренняя граница), ось симметрии потока, плоскость х = const (правая граница).

Для несущей фазы (ПД) условиями на левой границе служат соотношения на сильной ДВ; для дисперсной — условие непрерывности параметров частиц при переходе через волну детонации. В

системе координат ДВ имеем:

* 1*1*ч+1** ~i D2 * -у

Pj =WiPot>,Pi = V"' ri = (tFITTFIJ'v "j

* * *

«sj = A r6j = ги

*

Скорость детонационной волны jD и начальная плотность ВВ р0

считаются известными. В спутном потоке — невозмущенный воз* * * *

дух: р р^, Т\ = О. В безразмерном виде при i = 0 для ПД, воздуха и частнц имеем

р = l/',.H/j.p= 1, Т= 1, ч = 1. с= 0; р = !/-;», р = 1 /г. Т = 1//. и = iii. г = О; /V = Tf = lft, ,,s = 2±1, rf = 0 где - начальная массовая концентрация частнц.

Положсния внутренней ii внешней границ определяются в ходе решения задачи. На каждом отрезке (ij, элементы внутренней и внешней границы аппроксимируются секущими, на которых записываются условия на контактном разрыве и косом скачке, соответственно. На плоскости х~ const (в расчетах i=5) ввиду сверхзвукового характера течения газа дополнительных условий не требуется. На оси симметрии — ;> = 0.

Уравнения движения, законы взаимодействия фаз, уравнения состояния вместе с граничными условиями полностью определяют осесимметричный стационарный разлет газовзвесп ПД гетерогенного заряда ВВ в воздухе.

В § 3.2 приведено описание метода расчета движения газовзвеси с частицами за плоскостью ЧЖ. Для расчета несущей фазы выбран стационарный аналог метода С.К. Годунова, для расчета дисперсной — дискретно^граекторный метод. Кроме описания стандартных процедур даны модификации метода С.К. Годунова, связанные с численной реализацией граничных условий на оси симметрии, контактной поверхности и ударной волне. Приведен алгоритм вычисления "размазанной" плотности среды частиц.

Использование стационарного метода С.К. Годунова возможно только для сверхзвуковых, в маршевом направлении, течений. Поэтому начальный слоя в численных расчетах располагается не на плоскости ЧЖ, где Mj = 1, а на некотором расстоянии дг = где имеется развитое сверхзвуковое течение. Построение такого течения подробно описано в пункте 3.2.3. Там же приведен алгоритм расчета углов наклона скачка в воздухе и контактной поверхности в окрестности фронта детонационной волны.

В § 3.3 описана методика выбора оптимальных параметров разностной схемы (число ячеек в газе, количество траекторий

частиц, коэффициент устойчивости схемы), обеспечивающих вычисления с заданной точностью и равную величину погрешностей, вносимых разбиением расчетной области эйлеровой и лаг-ранжевой сетками. Приведены результаты методических расчетов. Получены формулы, связывающие параметры разностной схемы с требуемой точностью расчета.

В § 3.4 дака постановка задачи о разлете газовзвеси ПД с твердыми частицами в приближении эффективного газа. Такая постановка оправдана для взрывчатых смесей с малоинерцнонны-ми частицами и позволяет рассчитывать разлет по "обычным" уравнениям газовой динамики с видсизменными свойствами газа. Кроме того, модель эффективного газа является хорош™ тестом для проверки корректности программ, составленных по алгоритму с взаимным влияюте.ч фаз.

Особенность задачи заключается в том, что в плоскости ЧЖ смесь неравновесна: vf 'Г/ Т$. Необходимо вычислить параметры смеси после перехода к равновесному состоянию н принять их в качестве граничных условий для уравнений эффективного газа. Это достигается с помощью законов сохранения массы, импульса и энергии, связывающих начальное "1" (неравновесное в плоскости ЧЖ) и конечное "е" (равновесное) состояния смеси:

Ръ А + Р| = Р^Ч + Ре + />=е4

^.«.•(ср7\ + и?/2) + + и%/2) =

= РеЛ-(сРТ* + '4/2) + + "е/2)

Р = Ръпт

где индекс "»" относится к параметрам газа, индекс 'У к параметрам частиц; остальные обозначения общепринятые.

Получены формулы, дающие решение задачи о вычислении

значений и параметров (газовой постоянной Ле, показателя адиабаты 7е, числа Маха Л^е) эффективного газа в плоскости ЧЖ.

Указаны области значений исходных параметров смеси (Л^,

т, —, 7), для которых стационарное решение задачи о перехо-ср

де к равновесному состоянию несправедливо (здесь а-6,

II ^

5 = т — уД — начальные массовая концентрация частиц, отношения скоростей и температур фаз, соответственно). Приведены формулы пересчета для сравнения результатов численных решений в постановках "запыленного" и эффективного газов.

Глава 4 посвящена описанию газодинамики разлета ПД с частицами по результатам численных расчетов и расчету параметров сжатого ударной волной слоя.воздуха в окрестности фронта ДВ.

В § 4.1 приведены параметры разлета ПД для четырех вариантов, различающихся содержанием твердых частиц. Скорости детонации взяты из данных экспериментов. Газодинамика разлета рассматривается на примере смеси с массовой концентрации дисперсной фазы о6 = 0,34. В расчетах варьировалась величина параметра релаксации Д отвечающая разным размерам частиц. Рассмотрены случаи разлета ПД без частиц (0 = 0) и в приближении эффективного газа (/?=«>). Расчеты при 0—0 охватывают также предельный случай разлета ЯД с крупными частицами, когда отсутствует взаимное влияние фаз. Влияние а5 на газодинамику разлета показано для взвесей ПД с мелкими (¡3 = 400) частицами.

Некоторые характерные особенности разлета, связанные с наличием частиц показаны на рис. 6 - 8. Рис. 6 иллюстрирует сходимость численного решения в постановке "запыленного" газа к решению для эффективного (кривая 7) при возрастании

1.6

1.4

1

4

( /

О

X 6

Рис. 7. Скорость газа на оси струи, а,= 0,34. 1- <3= 0 2- 0= 0,2 3- 0= 24 4- 0= ю

Л

/■ \

_ 1 .1 1-1 ИИ 1 I 1 Г 1.111 1 3 14 1 I |1Лпи

10^ 10 Ю-" 10"' 1 т

Рис. 6. Температура газа а, =0,34

Х= 2: 1- 0= О 2-/5= 400; Х= 5: 3- /3= О 4- /3= 400.

Наличие частиц приводит к формированию поступательно движущегося ядра за плоскостью ЧЖ в газовзвеси ПД с постоянными параметрами под начальной характеристикой эффективного газа с числом Иаха М- ("полки" на кривых 2 — 6). Размеры ядра уменьшаются с возрастанием радиуса частиц. Образование поступательного ядра в потоке и два этапа ускорения ПД с высокими значениями поперечных градиентов скорости газа показаны на рис. 7. Начальное ускорение происходит на участке релаксации фаз до образования ядра, вторичное - после его разрушения. Разлет газовзвеси за плоскостью ЧЖ сопровождается интенсивным охлаждением ПД (рис. 8, где показаны температуры ПД в двух поперечных сечениях расчетной области) и частиц с резким перепадом температур фаз в поперечном направлении. Расширяющиеся ПД окружены тонкий высокотемпературным слоем сжатого в скачке уплотнения воздуха (верхние участки кривых на рис. 8).

Для конденсированных ВВ характерны высокие скорости распространения детонации. Разлет ПД таких веществ сопровождается образованием мощных ударных волн, в которых воздух разогревается до температур, когда становятся существенными процессы диссоциации.

В § 4.2 выполнен расчет параметров сжатого слоя в окрестности фронта ДВ с учетом мгновенной диссоциации воздуха в скачке. Диссоциированный газ за скачком уплотнения описывался равновесным уравнением состояния с измененным показателем адиабаты. Показано, что учет диссоциации приводит к существенному изменению параметров в сжатом слое, уменьшению угла наклона скачка и незначительному увеличению угла разлета ПД. Кроме того, для плотных ВВ температура сжатого слоя воздуха оказывается выше температуры ПД в плоскости ЧЖ.

В заключении перечислены основные результаты диссертационной работы.

1. Разлет образца сопровождается двумя характерными стадиями: а) детонацией металлизированного взрывчатого наполнителя; Ь) растеканием и горением взвеси ПД с частицами в объеме взрывной камеры.

2. Детонация металлизированного ВВ протекает квазиста-ционарио. Для ДВ в металлизированном ВВ характерен искаженный (негладкий) фронт. ПД занимают симметричную конусообразную область за фронтом ДВ.

3. Измеренные значения скорости ДВ и угла разлета ПД в окрестности фронта от массовой концентрации частиц зависят немонотонно.

4. Растекание газовзвеси ПД с частицами в объеме взрывной камеры — процесс нестационарный. Взвесь в виде облака движется за ударной волной и создает протяженный по времени импульс давления в камере. Суммарный импульс давления в замкнутом объеме, создаваемый алюмшшзнровлнпьгм ВВ, превосходит импульс при детонации заряда ВВ без частиц.

5. Численные расчеты подтвердили возможность использования принятой осесиммегричной схемы разлета ПД. Структура течения в соответствии с этой схемой включает в себя: а) присоединенный к плоскости ЧЖ скачок уплотнения в воздухе; Ь) контактную границу области ПД; с) сепаратрису потока частиц, которая может пересекать границу ПД: (I) висячий скачок уплотненна, образующийся вблизи границы ПД. Падение висячего скачка на ось симметрии н появление диска Маха можно ожидать на расстояниях не менее трех калибров заряда от плоскости ЧЖ.

0. Расширяющиеся ЛД окружены тонким высокотемператур-

ним слоем сжатого в скачке уплотнения воздуха. Температура в слое сжатого воздуха меняется при удалении от плоскости ЧЖ заметно медленнее, чем в области ЛД.

7. Взаимодействие несущей я дисперсной фаз за плоскостью ЧЖ приводит к образованию поступательно движущегося ядра в потоке ПД. После зоны релаксации течение перестраивается и выходит на режим, характеризуемый постоянными параметрами некоторого эффективного газа. Размеры ядра уменьшаются с увеличением радиуса частиц. Область ядра ограничивается, как и следовало ожидать, характеристикой, угол наклона которой определяется числом Маха в эффективном газе.

8. Ускорение ПД протекает в два этапа: а) начальное ускорение непосредственно за плоскостью ЧЖ до образования ядра потока; Ь) вторичное - после разрушения ядра. Для этих этапов характерны высокие значения поперечных градиентов параметров.

9. Наличие частиц приводит к уменьшению давленая, плотности, температуры ПД в области формирования и существования ядра потока по сравнению с разлетом ПД заряда ВВ без частиц. После разрушения ядра давление, плотность и температура в ПД с частицами выше, чем в "чистом" газе.

10. Разлет газовзвеси сопровождается интенсивным охлаждением ПД за счет расширения и теплообмена с частицами. Охлаждение настолько велико, что температура ПД становится много меньше температуры окружающего ударно сжатого слоя воздуха. Нагрев частиц происходит в ядре потока ПД. После разрушения ядра частицы интенсивно охлаждаются, сохраняя максимум температуры на оси симметрии. Для дисперсной фазы характерен значительный перепад температуры, в поперечном направлении обрат-

»0 пропорциональный размеру частиц.

11. Учет диссоциации воздуха приводит к существенному изменению параметров в сжатом слое, уменьшению угла наклона скачка и незначительному увеличению угла разлета ПД.

12. Отношение температур воздуха в сжатом слое и ПД в плоскости ЧЖ зависит от начальной плотности ВВ: в плотных ВВ температура ПД оказывается ниже температуры воздуха в сжатом слое.

ПУБЛ11КАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ.

1. Гонор А.Л., Теверовский И.А., Тирская P.I1. Численные исследования разлета капли вязкой жидкости под действием импульса давления Ц Механика реагирующих сред и ее приложения. Сб. науч. тр. СО АН СССР. Новосибирск, 1989. С. 156 - 162.

2. Гилинский М.М., Теверовский U.A. и др. Методика численного решения некоторых задач для двух- и однофазных течений в соплах и струях, истекающих в затопленное пространство или спутный поток. М.: Отчет Института механики МГУ. К 3874. 1989.

3. Гилинский М.М:, Гонор А.Л., Теверовский М.А. Разлет цилиндрических объемов гетерогенных металлизированных сред // Математ. моделирование нестационарных процессов в механике сплошных сред. Сб. науч. тр. Mlffl, 1990. С. 69 - 83.