УДК 541.126.011.2
КРИТЕРИИ ЗАЖИГАНИЯ ЭНЕРГЕТИЧЕСКИХ МАТЕРИАЛОВ КОРОТКИМ ЛАЗЕРНЫМ И ЭЛЕКТРОННЫМ ИМПУЛЬСАМИ
А. В. Ханефт, В. А. Долгачев, Е. В. Дугинов, Г. А. Иванов
CRITERIA OF ENERGETIC MATERIALS IGNITION BY SHORT LASER AND ELECTRON PULSES
A. V. Khaneft, V. A. Dolgachev, E. V. Duginov, G. A. Ivanov
В работе рассмотрены критерии зажигания конденсированных взрывчатых веществ коротким лазерным и электронным импульсами. Рассмотрено влияние радиуса светового пучка, зависимость коэффициента поглощения от температуры, многократное отражение светового потока и плавление. Проведено сравнение результатов расчета с экспериментом.
This paper discusses the criteria of condensed explosives ignition by short laser and electronic pulses. The influence of the light beam radius, the dependence of the absorption coefficient on the temperature, the multiple reflection of a light beam and melting are investigated. A comparison of the calculations with the experimental data is shown.
Ключевые слова: критерий зажигания, взрывчатое вещество, лазерный импульс, электронный импульс, размерный эффект, коэффициент поглощения, многократное отражение, плавление.
Keywords: ignition criterion, explosive, laser pulse, electron pulse, size effect, absorption coefficient, multiple reflections, melting.
Введение
Два последних десятилетия начались интенсивно проводиться эксперименты по инициированию конденсированных взрывчатых веществ (ВВ) лазерными импульсами и пучками электронов наносекундной длительности [1 - 4]. Данные работы интересны как с точки зрения теории, так и практики. Например, интерес к лазерному инициированию ВВ обусловлен не только изучением механизма инициирования, а также возможностью их применения в качестве светодето-наторов, помехозащищенность которых значительно выше, чем электродетонаторов [5].
Для оценки пороговой плотности энергии W * импульсного излучения, приводящей к инициированию ВВ, необходимо иметь теорию, позволяющую получать простые и физически наглядные критерии. К такой теории относится работа А. А. Ковальского с соавторами [6]. В данной работе предложен критерий зажигания конденсированных ВВ, нагретых тепловым импульсом. Суть данного критерия заключается в том, что по окончании действия внешнего источника тепла скорость теплоприхода в результате химической реакции в твердом теле Q+ должна превысить
скорость теплоотвода в глубь вещества Q_ :
Q+s Q-. (1)
При этом наибольшая скорость выделения тепла вследствие экзотермической реакции происходит в узком приповерхностном слое - реакционной зоне.
Ширина приповерхностного слоя 21 определяется таким образом, чтобы скорость химической реакции в плоскости 2 = 21 была в е раз меньше, чем на поверхности (рис. 1). Вне слоя тепловыделением химической реакцией пренебрегают. Скорость теплоотвода определяется градиентом температуры на границе реакционной зоны 2 = 21. В этом случае условие зажигания конденсированного вещества запишется в виде [6]:
z1pQZ exp
E
RTm
= -x dT
dz
-1, (2)
где р - плотность вещества; Q, Е - тепловой эффект химической реакции на единицу массы вещества и энергия активации; X - частотный фактор;
Я - универсальная газовая постоянная; X - коэффициент теплопроводности; Тт - температура в плоскости 2 = 0; Т1 - температура в плоскости
2 = 2, .
Рис. 1. К расчету ширины реакционной зоны химической реакции в твердом теле
Из условия, что Q+ (0) / Q+ (21) = е следует [6], что температура
T
-. (3)
1 + ЯТт / Е
Распределение температуры определяется из решения инертной задачи.
В данной работе приведен обзор работ авторов, а также приведены новые результаты по численному моделированию инициирования ВВ.
1. Критерии зажигания взрывчатых веществ коротким лазерным импульсом
1.1. Одномерная задача.
Для плоского светового пучка уравнение теплопроводности имеет вид:
дТ д 2Т
рс = X —— + а(1 - А)10 (і) ехр(-аг) +
ді
дг1
+ рQZ ехрі -
Е
КТ
(4)
Уравнение (4) имеет следующие начальные и граничные условия:
Т (^,0) = То
дТ
дг
г=0
дТ
дг
= 0 .
(5)
г=Ь
W
10(і) = 7---(4/Тт ) ехР(-4 Ч 1т ),
6т
(6)
где тт - длительность переднего фронта импульса. Длительность лазерного импульса хг- = 1,19тт. При-
чем 110 (і)Л = W,
АТ„ =
а(1 - A)W
сР '
(9)
Из сравнения выражений для Т1 (3) и (8) следует, что ширина реакционной зоны равна [6]:
г1 = — 1п Е а
(1 + у)
(1 - УТ0 / ДТт )
Производная ёТ1
—1 = -аАГт ехР(-аг1) . аг
(10)
(11)
Исходя из выражений (10) и (11) критерий зажигания конденсированного ВВ коротким световым импульсом (2) примет вид [6]:
z1рQZ ехр
Е
кт
ХаАТт
Е (у)
где
у = КТт /Е , Е(у) =-
1 + у
(13)
1 -УТ)/ АТт
Критическая плотность энергии короткого лазерного импульса вычисляется по формуле:
W * = _срАТт
а(1 - А)
(14)
Здесь Ь - толщина образца; Т0 = 300 К - начальная температура образца; с - коэффициент теплопроводности; А - коэффициент отражения; 10(і) -плотность потока энергии лазерного импульса. Внешний теплоотвод в задаче не учитывается, т. к. длительность лазерного импульса и время задержки зажигания значительно меньше характерного времени внешнего теплоотвода.
Плотность потока энергии определяется выражением:
Согласно работе [7], критерий (12) не был апробирован сравнением с численными расчетами на ЭВМ. В связи с этим в работе [8] численно решено уравнение теплопроводности (4) с граничными условиями (5), моделирующим зажигание азида свинца коротким импульсом электронов с экспоненциальным профилем поглощения по образцу. В [8] учитывалось также выгорание образца. Результаты численных расчетов критической температуры зажигания ДТт и
*
критической плотности энергии W совпали с точностью ~ 2 % с результатами, полученными из критерия (12). Кроме того, результаты данной работы показали, что выгоранием ВВ к моменту зажигания можно пренебречь.
1.2. Влияние радиуса светового пучка лазерного импульса на условие зажигания конденсированного ВВ.
Запишем уравнение теплопроводности для гауссова пучка в цилиндрической системе координат:
дТ х рс эТ = х
дТ
г дг I дг
д 2Т ^ дг 2
(
+ а(1 - А)10 (і) ехр
где W - плотность энергии лазерного импульса.
Согласно [6], после воздействия короткого светового импульса распределение температуры в твердом теле определяется формулой:
ДТ = Т - Т0 =а(1 - А ехр(-а2). (7)
ср
Температура в плоскости 21 определяется выражением:
ДТ1 = Т1 - Т0 = ДТт ехр(- а21) , (8)
где
2 Л
- аг -
г0 у
+рQz ехр|- КТ
(15)
Уравнение (15) имеет следующие начальные и граничные условия:
Т (г, г,0) = Т0
дТ_
дг
дТ_
дг
г=0
дТ_
дг
дТ_
дг
= 0.
г = К
= 0 .
(16)
Здесь К0 и Ь - радиус и длина образца в виде цилиндра. Распределение интенсивности по сечению лазерного пучка в уравнении (13) имеет гауссову форму:
I(г, і) = 10(і)ехр(-г2/ г02), где г0 - характерный радиус пучка. Интенсивность 10(і) определяется формулой (6).
Распределение температуры по окончании короткого лазерного импульса определяется распределением поглощенной энергии:
Т(г, г) = Т0 + а(1 А')1¥ ехр сР
(
2
- аг -
(17)
'0 У
Критерий зажигания (1) в этом случае примет вид
(12) [9]:
т
+
+
г=0
г=Ь
(
= -Хлг,
пг1 г1рQZ ехр 21 дТ
E
м
дг
г=0
V КТт У
-
(18)
г=0
г1рQZ ехр
Е
V КТт у
= -ХАТт
а ехі
:р(-аг ) + -2-
4 г1 -^ехр
( „2 Л
V г0 У
(19)
При 2 = 0, согласно (16), температура на периметре окружности радиусом г равна:
( ..2 Л
Т1 = Т0 + АТт ехр
0У
Из сравнения выражений (3) и (20) получим:
1
1п
1 + у
(20)
(21)
.1 -УТ^ДТт
Подставив выражения для 21 и г1 (10) и (21) в
(19) получим критерий зажигания ВВ коротким лазерным импульсом в зависимости от радиуса светового пучка [9]:
z1рQZ ехр
Е
КТт
ХАТт
Е (у)
а+
4г1
(22)
0У
В статье [9] в данной формуле допущена опечатка. Критическая плотность энергии лазерного импульса вычисляется по формуле (14).
Для типичных ВВ параметр у << 1. В этом случае
УТт
УТт
сРКТт
дК0 Е
ехр
(ЕЛ
VКТт У
(
і, 1 = а
4 1
Л
2 + 2 г12 г12 У
то есть скорость тепловыделения в реакционном объеме V = пг121 равняется потоку тепла из этого объема в глубь вещества. Здесь г1 - радиус реакционного объема. Уравнение (18) с учетом (17) преобразуется к виду:
1.3. Влияние зависимости коэффициента поглощения от температуры на условие зажигания конденсированных ВВ лазерным импульсом.
Уравнение теплопроводности в цилиндрической системе координат с учетом зависимости коэффициента поглощения от температуры имеет вид:
дТ
ср--------= X
ді
1 А( дТ
г дг V дг
д2Т Л дг 2
(
+ а(1 - А)10 (і )ехр
+ PQZ ехр| -
Е
КТ
(25)
Решение уравнения (25) при поглощении короткого лазерного импульса ВВ при условии, что
а(Т) = а 0 + Р(Т - Т0),
(26)
а время задержки зажигания значительно превышает длительность лазерного импульса, имеет вид [10]:
Т(г,г) = Т0 +ОТт)(ЬА^ехр(- Гайг' ср 0
Л
(27)
Исходя из уравнений (18) и (27) критерий зажигания ВВ примет вид [10]:
=Х
г1рQZ ехр
АТт
Е
КТт
Е(у)[1 -вАТт /(а(Тт)Е(у))]
срАТт
а+
4г1
&
W =
а(Тт)(1 - А)
(28)
(29)
(23)
1 ДТт 1 аДТт '
а функция Е(у) и 1. При у << 1 критерий зажигания (22) с учетом выражений для гх и 21 (23) можно преобразовать к виду [10]:
tad = , (24)
где ^ - адиабатическое время разогрева и характерное время тепловой релаксации реакционного объема, определяемые выражениями
Здесь а = X / ср - коэффициент температуропроводности.
Следовательно, условию зажигания (1) можно придать новый физический смысл: критическая температура зажигания определяется из равенства адиабатического периода индукции экзотермической реакции характерному времени тепловой релаксации реакционного объема.
В статье [10] в выражении для Е (у) имеется опечатка.
Результаты расчета зависимости критической плотности энергии зажигания ВВ от радиуса пучка лазерного импульса при решении уравнения теплопроводности (25) с параметрами РЕТО (С5И8М4012-
ТЭН) при а0 = 150 см-1 приведены на рис. 2 (кривые 1 - 3) [10]. Из рис. 2 видно, что если коэффициент поглощения зависит от температуры, то в случае в > 0 (кривая 3) критическая плотность энергии лазерного импульса уменьшается, а зависимость критической плотности энергии W * от радиуса пучка становится более пологой относительно кривой (1), которая рассчитана при в = 0. В случае в < 0 (кривая 2) критическая плотность энергии лазерного импульса, наоборот, увеличивается, а зависимость W * от радиуса пучка становится более резкой по сравнению с кривой (1).
+
2
г
+
г,
г
г
г, = г
г
Рис. 2. Влияние радиуса светового пучка лазерного импульса на критическую плотность энергии зажигания взрывчатого вещества при различных зависимостях коэффициента поглощения от температуры: кривые 1, 2, 3 - решение уравнения (2), точки 4, 5, 6 - решение по критерию (2, 4 - в = -0,1
см
1 К-1; 1, 5 - в = 0 ; 3, 6 - в = 0,1 см-1-К-1)
ж * =-
2
р[с + Я/8(Т - Т/)]
+ а(1 - А)/0(і)ехр
дТ
ді
(
= Х
г
г0
1_д_
г дг
2 Л
дТ Л д 2Т
2Л
+
дг У дг
Е
(30)
W =-
а(1 - А) а(1 - А)
(31)
На рис. 3 представлены результаты расчета изменения температуры АТ5 в центре поверхности РБТМ вблизи порога зажигания при воздействии лазерного импульса длительностью 40 нс и радиусе светового пучка г0 = 10- см. Как видно из данного рисунка, при АТ5 = 113 К на зависимости АТ5 = /(і) наблюдается «полочка», которая обусловлена плавлением ВВ. Слева от полочки ВВ не расплавлено, а справа от нее приповерхностный слой расплавлен.
На рис. 2 изображены также результаты расчетов критической плотности энергии лазерного импульса, вычисленные по формулам (28), (29) (точки 4, 5 и 6). Наилучшее совпадение дают результаты расчета при в = 0 (кривая 2, точка 5). Увеличение ошибки при в Ф 0 обусловлено тем, что при выводе критерия (28) полагалось, что коэффициент поглощения изменяется мгновенно. При численном решении уравнения теплопроводности а изменяется вследствие нагрева при поглощении энергии.
Ошибка в вычислении критической плотности энергии по формуле (14) уменьшается, если брать среднее значение критической плотности энергии
W\а(Тт)) + W *(а 0)
Рис. 3. Зависимость температуры ДТ5 на поверхности образца в центре пучка от времени при W = 16,3 (1), 16,4 (2) и 16,5 Дж/см2 (3)
На рис. 4 приведены результаты расчета зависимости критической плотности энергии инициирования ВВ лазерным импульсом длительностью хг- = 40 нс от радиуса пучка без учета плавления (кривая 1) и с учетом плавления (кривая 2). Как видно из сравнения двух кривых, плавление повышает критическую плотность энергии инициирования. Кривая 2 сдвинута вверх относительно кривой 1 на одну и ту же величину Д W ~ 6 Дж/см2.
1.4. Влияние плавления на условие зажигания конденсированных ВВ лазерным импульсом.
В случае плавления уравнение теплопроводности запишется в виде:
где Ну - скрытая теплота плавления; 8(х) - дельт-функция; Ту - температура плавления.
Процесс плавления не влияет на критическую температуру зажигания при воздействии короткого лазерного импульса. Процесс плавления влияет только на время задержки инициирования. Выражение для критической температуры зажигания ВВ коротким лазерным импульсом в зависимости от радиуса светового пучка в случае независимости коэффициента поглощения от температуры определяется формулой (22), а критическая плотность энергии определяется формулой [11]:
СРДТтах + рНУ
Рис. 4. Влияние радиуса светового пучка лазерного импульса на критическую плотность энергии зажигания РЕТ№ кривая 1 - расчет без учета плавления; кривая 2 - расчет с учетом плавления;
■ - расчет по критерию (формулы (22) и (31))
1.5. Влияние многократного отражения светового потока на условие зажигания конденсированных ВВ лазерным импульсом.
Уравнение теплопроводности с учетом плавления и многократного отражения светового потока от противоположных сторон образца имеет вид:
+
г
- аг
р[с + Иґ 5(Т - Тґ)] дТ = Х
(1 д ( дТЛ д2ТЛ
г дг V дг
дг 2
(
+ а(1 - А1) 10 (і) ехр
„2 Л
х [1 + А2 ехр(2а(г - щ +рQz ехр(^_£.
[1 - А! А2ехр(-2аЬ)] V КТ
(32)
ZlPQZ ехр1 -
Е_
КТ
Е (У)
1 2 аЬ --1п( А2 Е 2(у))
4г1
(33)
Выражение для критической плотности энергии лазерного импульса с учетом многократного отражения светового пучка имеет вид:
W * = [1 - А^2 ехр(-2а£)] х
[1 + А1 ехр(-2аЬ)]
срАТт
рН/
а(1 - А1) а(1 - А1)
(34)
Рис. 5. Зависимость пороговой энергии инициирования РЕТМ от радиуса светового пучка: линия - численный расчет при а = 0,065 см-1,
■ - эксперимент [12], ▲ - расчет по критерию
где А1 - коэффициент отражения поверхности на входе светового пучка в образец; А2 - коэффициент отражения тыльной поверхности на выходе светового пучка из образца.
Критерий зажигания ВВ с температурой плавления ниже температуры зажигания коротким лазерным импульсом гауссовой формы по сечению в области прозрачности образца с учетом многократного отражения определяется формулой:
По формулам (33), (34) рассчитана зависимость критической плотности энергии зажигания РЕТМ от радиуса светового пучка, а также энергия лазерного импульса:
Е * = ж& *. (35)
Результаты расчетов приведены на рис. 5. Как видно из данного рисунка, результаты расчетов по критерию достаточно хорошо совпадают как с результатами численного решения уравнения теплопроводности (32), так и с экспериментом.
Рис. 6. Зависимость критической энергии зажигания РЕТМ от коэффициента отражения А2 тыльной поверхности образца при г0 = 0,25 мм и А1 = 0,036 (1), 0,3 (2) и 0,6 (3)
На рис. 6 приведены результаты расчета критической энергии зажигания РЕТМ на основе формул (33)
- (35) в зависимости от коэффициента отражения тыльной поверхности образца. Расчеты проведены при радиусе светового пучка г0 = 0,25 мм и коэффициенте отражения внешней поверхности А1 = 0,036 (кривая 1), 0,3 (кривая 2) и 0,6 (кривая 3). Естественно, что чем больше А1 , тем больше Е* . Коэффициент отражения тыльной поверхности образца менялся от френелевского до коэффициента отражения металла. Как видно из рисунка, что чем больше коэффициент отражения А1 , тем в более широких пределах можно
регулировать энергию зажигания Е* , изменяя коэффициент отражения А2.
2. Критерии зажигания взрывчатых веществ коротким электронным импульсом
2.1. Одномерная задача.
Распределение поглощенной энергии пучка электронов в твердом теле имеет вид кривой с максимумом (рис. 7). На рис. 7 Д21 есть ширина реакционного слоя, Яеу и Яех есть эффективная и экстраполированная длины линейного пробега электрона.
Рис. 7. Качественная кривая распределения плотности поглощенной энергии пучка электронов
Экспериментальную кривую распределения плотности поглощенной энергии твердым телом обычно
+
+
-аг -
х
г,
г
X
аппроксимируют полиномом третьей степени. В случае PETN:
WI Wm - Л(|) - 0,7 +1,57| - 2,З1|2 + 0,61|З , (З6)
где | - z I Ref. При I>Iex Л(|) - 0, а интеграл
Ie
ex
-1,
где ^ = 1,44, Яу = 173,6 -10-4 см,
Кех = 1,44 Яе/ = 250 • 10-4 см.
Уравнение теплопроводности с учетом плавления в случае нагрева ВВ электронным пучком имеет вид:
P[c + Hf 5(T - Tf)]
дт
=к-д!Г+
dt dz2
+!<tm+pqz expf-E. 1. Ref І RTу
(З7)
При малом времени воздействия импульса электронов распределение температуры в образце по окончании импульса определяется зависимостью:
T (I) - To +
W
cpRef
Л(|).
(ЗВ)
Критическая температура зажигания определяется из условия:
AzjpQZ exp
Г E Л
RT»
К дТ ' Ref д1
При I -1» согласно (ЗВ) температура
Tm - To +
W
cPRef
(З9)
(40)
Г E Л
V RTrn у а ширина реакционного слоя
KAT» д2Л»
R* dI2
(42)
Az1 - R~ef I-
2T1TmЛ rnR
at»eд2Л»I dI2
(4З)
(44), Теплофизические и кинетические параметры ВВ приведены в табл. 2.
Таблица 1
W, Дж^м2 PETN RDX HMX
По критерию (40, 42) 14,ЗВ З2,99 21,1В
Численный расчет 15,0 З4,0 22,0
Таблица 2
Параметр PETN RDX HMX
E, 196,6 197,З 220,В
кДжЛмоль [14] [14] [14]
Z , c- 6,3.1019 [14] 2,02 •101В [14] 5 •Ю19 [14]
Q, МДж7кг 1,26 [14] 2,1 [14] 2,1 [14]
c, Дж!(кг-К) 1255,2 2092 1250
[15] [1В] [20]
К, ВтI(м•К) 0,250В 0,105 0,29З
[14] [14] [14]
p, кг!мЗ 1,77 •10З 1,В2 •10З 1,9 •10З
[16] [19] [19]
Hf, кДж!кг 19З [17] 2З5,5 [17] 192,46 [21]
Tf , К 41З [14] 476 [14] 55В [14]
Экстраполированная длина пробега электронов с начальной энергией Е0 = 250 кэВ для КЭХ и ИМХ вычислялась по эмпирической формуле, предложенной в работе [18]:
Получим выражение для критической температуры зажигания ВВ. Для этого разложим в ряд Тейлора функцию Л(§) в окрестности |m [13]:
Л^-и = Л(^) + ^ 2m )2 д ^ m ) . (41)
2 дс,
Исходя из выражений (3), (38), (40) и (41) критерий зажигания (39) примет вид:
pQZ exp
1 i/i ч a3e0
—ln(1+a2eo) - -----------ar
a2 1 + a4e0a5
(45)
Здесь р - плотность среды в г/см ; е0 = Е0/тс2 (тс2 - энергия покоя электрона, равная 511 кэВ); постоянные:
0,2335Ае
а =-
1,209 f , a2 - 1,7В •10-4 Zef ,
Zef
a3 - 0,9В9 - З • 10-4 Zef , a4 - 1,46В - 1,1В • 10-2 Zef ,
1,2З2
Критическая плотность энергии зажигания ВВ электронным пучком с учетом плавления:
W * =рКе/ (сАТт + Нг). (44)
Численное решение уравнения (37) с длительностью импульса 15 нс позволило определить критическое значение плотности энергии электронного пучка для РБТК Кроме РБТМ расчеты были проведены для ИБХ (С3И6М606-гексоген) и НМХ (С4И8М808 -октогена). Данные результаты расчетов, как это видно из табл. 1, достаточно хорошо согласуются с результатами, вычисленными по критерию (формулы (42) и
5 7 0,109 ’
Zef
где Ае^ - эффективная атомная масса; Zef - эффективный атомный номер. Эффективная атомная масса и эффективный атомный номер вычисляются по формулам:
п z
Zef =£^ , А^ е
І = 1
ef
(X / А)еу
(X /А)еу = Х (УХ /А), У = А / .
г=1 г=1
Здесь Хг, Аг - атомный номер и атомная масса
г -го элемента; /г - его весовая доля; п - число
весовых долей. Результаты расчетов Zf и Aef для PETN, RDX и HMX приведены в табл. З.
Таблица 3
Параметр PETN RDX HMX
Zef 69,2 39,7 53
Aef 135 77,37 103,2
Рис. 8. Зависимость эффективной длины пробега электронов от начальной энергии электронов:
1 - РЕТЩ 2 - НМХ, 3 - КОХ
Рис. 9. Зависимость критической плотности энергии зажигания ВВ от Яеу : 1 - РЕТМ,
2 - НМХ, 3 - ЯВХ
На рис. 10 представлены результаты расчета распределения температуры в РЕТЫ в области поглощения пучка электронов с плотностью энергии
W = 15,5 Дж/см2. Появления полочки на кривых (1)
- (3) обусловлено плавлением. Справа от полочки РБТЫ не расплавлен, а слева от полочки РБТЫ расплавлен.
Оценка Яеу для РЕТЫ с использованием формулы
(45) при |ех = 1,44 показал, что расчетное значение примерно на 14 % меньше Яеу , полученного из обработки эксперимента. Поэтому расчетные значения Яеу при моделировании зажигания гексогена и окто-
гена электронным пучком были увеличены соответственно на 14 %. Расчеты показали, что Яеу линейно
возрастает с увеличением начальной энергии электронов (рис.8). Это приводит к возрастанию критической энергии зажигания ВВ в области поглощения пучка электронов (рис. 9).
Рис. 10. Динамика распределения температуры в области поглощения пучка электронов при зажигании РЕТМ: / = 10-4 (1), 6 -10-4 (2) и 6,42 -10-4 с.
На рис. 11 представлены результаты расчетов за*
висимости времени задержки і инициирования ВВ в области поглощения пучка электронов. Критическая плотность энергии электронного пучка, согласно табл. 1, согласуется с экспериментом [4]. Однако рас*
четное значение времени задержки і для РБТЫ превышает экспериментальное примерно на 2,5 порядка. Таким образом «классический тепловой взрыв» не позволяет описать полностью эксперименты по инициированию РЕТЫ пучком электронов.
Рис. 11. Зависимость времени задержки зажигания ВВ от плотности энергии импульса электронов (линии расчет: 1 - РЕТЩ 2 - НМХ, 3 - КОХ,
точка - эксперимент [4]
Данный эксперимент достаточно хорошо объясняется при учете термоупругих напряжений, которые, как известно, влияют на энергию активации химической реакции [23], а также радиационно-термическим механизмом [24; 25]. Следует все таки заметить, что экзотермическая реакция вносит основной вклад в процесс зажигания ВВ электронным пучком наносе-кундной длительности.
Литература
1. Адуев, Б. П. Взрывное разложение азидов тяжелых металлов / Б. П. Адуев, Э. Д. Алукер, Г. М. Белокуров [и др.] // ЖЭТФ. - 1999. - Т. 116. - Вып. 5(11). - С. 1676 - 1693.
2. Медведев, В. В. Влияние неоднородного облучения на пороги зажигания двухосновного пористого топлива / В. В. Медведев // Химическая физика. - 2009. - Т. 28. - № 6. - С. 74 - 76.
3. Oleshko, V. I. Explosive Decomposition of FTDO Initiated by Laser Radiation Pulse and pulse of Accelerated Electrons / V. I. Oleshko, V. P. Tsipilev, V. V. Lysyk [et al.] // Известия вузов. Физика. - 2012. - Т. 55. - № 11/3. -С. 158 - 161.
4. Адуев, Б. П. Исследование ранних стадий взрывного разложения кристаллов тетранитропентаэритрита при инициировании импульсными электронными пучками / Б. П. Адуев, Г. М. Белокуров, С. С. Гречин [и др.], // Известия вузов. Физика. - 2007. - № 2. - С. 3 - 9.
5. Илющин, М. А. Влияние добавок ультрадисперсных частиц углерода на порог лазерного инициирования полимерсодержащего светочувствительного взрывчатого состава / М. А. Илющин, И. А. Угрюмов [и др.] // Химическая физика. - 2005. - Т. 24. - № 10. - С. 49 - 56.
6. Ковальский, А. А. К вопросу о зажигании баллиститных порохов / А. А. Ковальский, С. С. Хлевной,
В. В. Михеев // Физика горения и взрыва. - 1967. - Т. 3. - № 4. - С. 527 - 541.
7. Вилюнов, В. Н. Теория зажигания конденсированных веществ / В. Н. Вилюнов. - Новосибирск: Наука, 1984. - 190 с.
8. Ханефт, А. В. К инициированию азида свинца электронным импульсом / А. В. Ханефт // Физика горения и взрыва. - 1993. - Т. 29. - № 5. - С. 63 - 66.
9. Ханефт, А. В. Влияние распределения светового потока в лазерном пучке на критическую энергию зажигания конденсированного вещества / А. В. Ханефт // Химическая физика. - 1998. - Т. 17. - № 10. - С. 67 - 70.
10. Дугинов, Е. В. Влияние зависимости коэффициента поглощения от температуры на критическую энергию зажигания конденсированного вещества лазерным импульсом / Е. В. Дугинов, А. В. Ханефт // Физика горения и взрыва. - 2011. - Т. 47. - № 4. - С. 127 - 136.
11. Ханефт, А. В. Влияние плавления на критическую энергию зажигания конденсированного взрывчатого вещества коротким лазерным импульсом / А. В. Ханефт, Е. В. Дугинов // Физика горения и взрыва. - 2012. -Т. 48. - № 6. - С. 47 - 53.
12. Быхало, А. И. Инициирование ТЭНа мощным лазерным излучением / А. И. Быхало, Е. В. Жужукало, Н. Г. Ковальский [и др.] // Физика горения и взрыва. - 1985. - Т. 21. - № 4. - С. 110 - 113.
13. Ханефт, А. В. Критерий зажигания конденсированных веществ электронным импульсом / А. В. Ханефт. // Химическая физика. - 1998. - Т. 17. - № 8. - С. 132 - 137.
14. Физика взрыва; под ред. Л. П. Орленко. - М.: Физматлит, 2004. - Т. 1. - 824 с.
15. Беляев, А. В. Переход горения конденсированных систем во взрыв / А. В. Беляев, В. К. Боболев, А. И. Коротков [и др.]. - М.: Наука, 1973. - 292 с.
16. Детонация и взрывчатые вещества: сб. ст.; под ред. А. А. Борисова. - М.: Мир, 1981. - 392 с.
17. Баум, Ф. А. Термостойкие взрывчатые вещества и их действие в глубоких скважинах / Ф. А. Баум, А. С. Державец, Н. Н. Санасарян. - M.: Недра, 1969. - 160 с.
18. Мейдер, Ч. Численное моделирование детонации / Ч. Мейдер. - М: Мир, 1985. - 384 с.
19. Шалл, Р. Физика детонации / Р. Шалл // Физика высоких плотностей энергии. - М.: Мир, 1974. - 488 с.
20. Ассовский, И. Г. Физика горения и внутренняя баллистика / И. Г. Ассовский. - М.: Наука, 2005. - 358 с.
21. Струнин, В. А. Моделирование горения октогена / В. А. Струнин, Л. И. Николаева, Г. Б. Манелис // Химическая физика. - 2010. - Т. 29. - № 7. - С. 63 - 70.
22. Tatsuo Tabata and Rinsuke Ito. An Algoritm for the Energy Deposition by Fast Electrons // Nuclear Science and Engineering. - 1974. - V. 53. - P. 226 - 239.
23. Ханефт, А. В. Моделирование инициирования ТЭНа пучком электронов наносекундной длительности / А. В. Ханефт, Е. В. Дугинов, Г. А. Иванов // Химическая физика и мезоскопия. - 2012. - Т. 14. - № 1. - С. 28 -39.
24. Khaneft, A. V. Radiation-thermal mechanism of PETN initiation in the absorption / A. V. Khaneft, G. A. Ivanov // Известия вузов. Физика. - 2012. - Т. 55. - № 11/3. - С. 71 - 75.
25. Khaneft, A. V. Radiation-thermal mechanism of initiation of PETN in the absorption region of the electron beam. Energetic Materials: 43th International Annual Conference of ICT, FRG, Karlsruhe, 26 - 29 iune / A. V. Khaneft, G. A. Ivanov. - 2012. - Р. 17-1 - 17-11.
Информация об авторах:
Ханефт Александр Вилливич - доктор физико-математических наук, профессор кафедры теоретической физики КемГУ, 8-9132947480, [email protected].
Alexander V. Khaneft - Doctor of Physics and Mathematics, Professor at the Department of Theoretical Physics, Kemerovo State University.
Долгачев Вадим Александрович - аспирант кафедры теоретической физики КемГУ, 89505843359, [email protected].
Vadim A. Dolgachev - post-graduate student at the Department of Theoretical Physics, Kemerovo State University.
Дугинов Евгений Владимирович - кандидат физико-математических наук, доцент кафедры общенаучных дисциплин Кемеровского государственного сельскохозяйственного института, 89502719570,
Evgeniy V. Duginov - Candidate of Physics and Mathematics, Associate Professor, Department of General Scientific Disciplines, Kemerovo State Agricultural Institute.
Иванов Георгий Анатольевич - аспирант кафедры теоретической физики КемГУ, 89236094492, [email protected].
Georgiy A. Ivanov - post-graduate student at the Department of Theoretical Physics, Kemerovo State University.