CC BY
34
ПРИКЛАДНАЯ МАТЕМАТИКА И МЕТОДЫ МАТЕМАТИЧЕСКОГО МОДЕЛИРОВАНИЯ
УДК 536.2+534.222.2
А. В. Аттетков, Е. В. Пилявская
оценка критических условий
удАрно-волнового инициирования
химической реакции в двухфазном пористом энергетическом материале
Исследован процесс диссипативного разогрева ударно-сжатого энергетического материала, содержащего сферические поры с термически тонким покрытием их поверхности. Полученные теоретические результаты использованы для определения порогового давления инициирования очаговой химической реакции в изучаемой системе.
E-mail: e.pilyavskaya@mail.ru
Ключевые слова: двухфазный пористый энергетический материал,
ударно-волновое нагружение, критические условия инициирования химической реакции.
В исследованиях по ударно-волновому инициированию пористых энергетических материалов (ЭМ) важное место занимает вязкопла-стическая модель горячих точек — локализованных зон динамического перегрева ударно-сжатого ЭМ. Несмотря на достигнутые результаты по рассматриваемой проблеме (см., например, обзор [1] и цитируемую в нем литературу), ряд вопросов является еще открытым и требует дальнейшего изучения. В частности, актуальным остается вопрос о влиянии тонкого слоя пластификатора (в дальнейшем — термически тонкого покрытия) на процесс диссипативного разогрева и критические условия инициирования очаговой химической реакции в ударно-сжатом пористом ЭМ. Его изучение является основной целью проведенного исследования.
Рассмотрим задачу [1-7] о стационарной ударной волне (УВ), распространяющейся со скоростью D в двухфазном пористом ЭМ — несжимаемой вязкопластической среде (фаза 1), содержащей сферические поры радиуса а с покрытием (фаза 2; несжимаемая вязкая среда) их поверхности (регулярная ячеистая схема [8]; наличием газа в порах пренебрегается). При этом приняты следующие предположения:
1) характерная длина УВ много больше размера пор и расстояния между ними [2-7];
2) тепловой контакт на границе фаз 1 и 2 идеальный, покрытие является термически тонким и допустима реализация модели "сосредоточенная емкость" [9, 10];
3) объемным содержанием фазы 2 в единице объема пористого ЭМ можно пренебречь;
4) реализуется сильновязкий режим [1] затекания пор с сохранением их сферической формы.
При сделанных предположениях интегралы уравнений сохранения массы и импульса двухфазного пористого материала в системе координат, связанной с УВ, можно представить в виде [5]
Р - Ро = Р1 П2(а0 — а)а-2, (1)
справедливом для всех промежуточных состояний во фронте волны (индекс "0" характеризует параметры состояния перед фронтом УВ), где
2У ао
Ро = "Т" 1п ^---(2)
3 о + а0 — 1
— амплитуда упругого предвестника; У — предел текучести материала фазы 1. Давление р и плотность р двухфазного пористого ЭМ в (1), (2) связаны с фазовыми значениями величин выражениями [5-7]
р = а-1 [0р2 + (1 — 0)рх]; р = а-1 [0р2 + (1 — 0)р1] « а-1 р1;
(3)
с3 — a3
а = —-3; о =
^3
Ь3 — а3' Ь3 — а3' где Ь — радиус сферического объема характерного (представительного) элемента двухфазного материала; с — радиус контактной границы фаз; а, о — концентрационные симплексы подобия изучаемого ЭМ.
Отметим, что реализация рассматриваемого (сильновязкого) режима затекания пор приводит к существованию волны с монотонным профилем (равновесное состояние позади фронта УВ в дальнейшем будет обозначено индексом е), при этом скорость а деформации двухфазного пористого ЭМ внутри фронта УВ определяется соотношением [5, 7]
, N 3а(а — 1)(0 + а — 1) ~
а(а) = —~г—т-т^—-г—-ттг piD х
v ! 4[^а0 + п(а — 1)(1 — 0)]
х
а0 — а 2Y ао(0 + а — 1) +--i^r In
(4)
а2 3р1 V2 а(0 + а0 — 1) где п, Р — коэффициенты вязкости фаз 1 и 2 соответственно.
Процесс формирования температурного поля в ударно-сжатом пористом ЭМ определяется конкуренцией мезоскопических процессов тепловой диссипации (вследствие пластических деформаций материала фазы 1 и работы вязких сил фаз 1, 2) и теплопереноса. На стадии,
предшествующей возбуждению очаговой химической реакции в рассматриваемой двухфазной системе с термически тонким покрытием поверхности затекаемых пор, математическая модель изучаемого процесса имеет вид [5, 6]
ClPl
dT (r,t) dT (r,t)
+ vr
dt
dr
V / V \ 2
- + 12n iyv^j , c(t) < r < b(t), t > 0;
Ai d
r2 dr 'V-N 2 r
>ßT(r, t) dr
r(t)lt=o = ro; v(r,t)lt=o =0; T(r,t)lt=o = T0;
dT (r, t)
dr dT(r, t)
= 0;
r=b(t)
dr
c(t)
(5)
r=c(t) 3A1 5 + a - 1
dT (r, t)
X
X
c2P2
dt
а2
r=c(t) 3 (а - 1)(5 + а - 1)
Vr =
3 .
a'Oa
3(ao - 1)r2
= r0 + a3(t) - a3,
где Т — температура; ск, Ак,к Е {1, 2}, — коэффициенты удельной теплоемкости и теплопроводности фаз 1 и 2 соответственно; скорость а деформации двухфазного пористого ЭМ во фронте УВ определяется равенством (4).
Математическая модель (4), (5) представляет собой смешанную задачу нестационарной теплопроводности, в которой наличие термически тонкого покрытия на поверхности пор фактически учтено граничным условием на подвижной границе фаз при г = с(£), явно содержащим производную от температуры по времени. По принятой терминологии [5, 6, 9, 10] данную математическую модель принято называть моделью "сосредоточенная емкость".
Для достижения основной цели исследования предположим, что процесс теплопереноса не оказывает существенного влияния на формируемое температурное поле в ударно-сжатом пористом ЭМ, формально полагая в (5) А1 = 0. Это позволяет представить скорость изменения температуры на подвижной границе фаз г = с(£) в виде
clPl
dT (r, t)
dt
4^ cipi
а
r=c(t)
3 C2P2 (а - 1)(5 + а - 1)'
(6)
а приращение удельной (на единицу объема) диссипированной энергии на границе фаз в* = с1р1[Т(с, ¿) — Т0] записать как
r
5
r
„ _ 3^cipi Ре
бе - .
* = 4^С1Р1 Г а (а') ¿а'
6 = 3 сгр2} (а' - 1)(8 + а' - 1)' ( )
ао
В УВ с амплитудой ре ^ У, достаточной для практически полного затекатия пор (ае ^ 1), основной вклад в мезоскопический процесс тепловой диссипации анализируемого (сильновязкого) режима их затекания вносят вязкие свойства фаз изучаемого пористого ЭМ [7]. В пренебрежении объемным содержанием фазы 2 (8 ~ 0) это позволяет записать приращение удельной диссипированной энергии контактной границы фаз за фронтом УВ в виде
2
4 С2Р2 П2'П'
где тп = 4п/ре — характерное время затекания поры [1, 11-13]. Далее предполагая, что необходимым условием возбуждения очаговой химической реакции в ударно-сжатом двухфазном пористом ЭМ является условие достижения заданной температуры Тгдп на контактной границе фаз, приращение удельной (на единицу объема) диссипиро-ванной энергии е*дп = с1р1(Тгдп — Т0) при ударном переходе можно представить как
* _ 3С *
6гдп ^ Ргдп ' (8)
где
^ = ^ £ = С2Р2 п' С1Р1
— симлексы подобия вязких и физических свойств изучаемого материала соответственно. Здесь р*дп — предел возбуждения очаговой химической реакции в двухфазном пористом ЭМ при ударно-волновом
нагружении, определяемый, согласно (8), выражением
* ^ *
Ргдп 6гдп'
(9)
Для типичных ЭМ (таких, например, как тротил, ТЭН, гексоген), принимая с1р1 ~ 2 • 106 Дж/(м3^К)и Тгдп — Т0 = 400 ... 600 К, находим е*дп = 0,8''' 1,2 ГПа [4]. Отметим, что оценка (9) предела ударно-волнового инициирования очаговой химической реакции в двухфазном пористом ЭМ (нижнего предела зажигания по давлению [1]) отмечена верхним индексом *, поскольку она может отличаться от фактического порогового давления ргдп возбуждения очаговой реакции в изучаемой системе вследствие проявления эффектов теплопроводности в процессе затекания пор (обычно ргдп > р*дп). Для однофазного (при отсутствии покрытия на поверхности пор) пористого ЭМ подробное
изучение данного вопроса проведено в работах [1,4], для двухфазного — в работе [5].
Теоретическая оценка (9) физически правильно отражает существующие представления о процессе ударно-волнового инициирования очаговой реакции в пористых ЭМ при наличии тонкого слоя пластификатора на поверхности затекаемых пор. В частности, как следует из равенства (9), при Z < 1 термически тонкое покрытие оказывает флегматизирующее действие на стадии генерации очагов разложения в ударно-сжатом пористом ЭМ [14, 15]. Отметим, что симплекс подобия е физических свойств фаз ЭМ в этом случае близок к единице. Практический интерес представляет также случай Z > 1, при реализации которого термически тонкое покрытие оказывает сенсибилизирующее действие на процесс возбуждения очаговой химической реакции в пористом ЭМ при ударно-волновом нагружении. Его физическая интерпретация ассоциируется, например, с анализом проблемы ударно-волновой чувствительности водоэмульсионных взрывчатых веществ, сенсибилизируемых полыми стеклянными микросферами [13, 16].
Работа выполнена при частичной финансовой поддержке гранта Президента РФ по государственной поддержке ведущих научных школ (грант НШ-255.2012.8).
СПИСОК ЛИТЕРАТУРЫ
1. Х а с а и н о в Б. А., Аттетков А. В., Борисов А. А. Ударно-волновое инициирование пористых энергетических материалов и вязкопластическая модель горячих точек // Химическая физика. - 1996. - Т. 15, № 7. - С. 53-125.
2. Дунин С. З., Сурков В. В. Динамика закрытия пор во фронте ударной волны // Прикладная математика и механика. - 1979. - Т. 43, № 3. - С. 511-518.
3. Дунин С. З., Сурков В. В. Эффекты диссипации энергии и влияние плавления на ударное сжатие пористых тел // Журнал прикладной механики и технической физики. - 1982. - № 1. - С. 131-142.
4. Аттетков А. В., Соловьев В. С. О возможности разложения гетерогенных ВВ во фронте слабой ударной волны // Физика горения и взрыва. - 1987. -Т. 23, №4.-С. 113-125.
5. A 11 e t k o v A. V., Golovina E. V., Ermolaev B. S. Mathematical simulation of mesoscopic processes of heat dissipation and heat transfer in a two-phase porous material subjected to shock compression // Journal of Heat Transfer Research. - 2008. - Vol. 39, no. 6. - P. 479-487.
6. Аттетков А. В., Головина Е. В., Ермолаев Б. С. Иерархия моделей процесса теплопереноса в двухфазном пористом материале при ударном сжатии // Труды V Всерос. национальной конф. по теплообмену. - М., 2010. - Т. 8. -С. 50-53.
7. Аттетков А. В., Пилявская Е. В. Эффекты тепловой диссипации при распространении ударной волны в двухфазном пористом материале // Вестник МГТУ им. Н.Э. Баумана. Сер. Естественные науки. - 2011. - № 3. - С. 53-58.
8. НигматулинР. И. Динамика многофазных сред / В 2 ч. - М.: Наука, 1987.
9. Пудовкин М. А., Волков И. К. Краевые задачи математической теории теплопроводности в приложении к расчетам температурных полей в нефтяных пластах при заводнении. - Казань: Изд-во Казанского ун-та, 1978. - 188 с.
10. А т т е т к о в А. В., Волков И. К., П и л я в с к и й С. С. Иерархия математических моделей процесса теплопереноса в твердом теле со сферическим очагом разогрева, обладающим покрытием // Проблемы газодинамики и тепломассообмена в аэрокосмических технологиях: Труды XVII Школы-семинара молодых ученых и специалистов под руководством академика РАН А.И. Леонтьева. - М., 2009. - Т. 1. - С. 166-169.
11. Two-phase viscoplastic model of shock initiation of detonation in high density pressed explosives / B.A. Khasainov, A.A. Borisov et al. // VII Symp. (Int.) on Detonation. - Annapolis, 1981. - P. 435-448.
12. Maiden D., N u 11 G. A hot-spot model for calculating the threshold for shock initiation of pyrotechnic mixtures // XI International Pyrotechnics Seminar Vail. -Colorado, 1986. - P. 813-826.
13. ХасаиновБ. А., Ермолаев Б. С. Возбуждение химической реакции при ударно-волновом сжатии жидких ВВ, содержащих стеклянные микросферы // Химическая физика. - 1992. - Т. 11, № 11. - С. 1588-1600.
14. Ингибирование разложения гексогена невзрывчатыми добавками при механических и ударно-волновых воздействиях / С.Г. Андреев, В.В. Зюзин и др. // Химическая физика. - 1998. - Т. 17, № 1. - С. 45-54.
15. Физика взрыва / Под ред. Л.П. Орленко. - В 2 т. T. 1. - М.: Физматлит, 2002. - 832 c.
15. K h a s a i n o v B. A., E r m o l a e v B. S., P l e s l e s H. N. Shockwave initiation of chemical reaction in liquid high explosives sensitized by glass microballoons // X Symp. (Int.) on Detonation. - Boston, 1993. - P. 40-43.
Статья поступила в редакцию 13.10.2011
Александр Владимирович Аттетков родился в 1955 г., окончил МВТУ им. Н.Э. Баумана в 1979 г. Канд. техн. наук, старший научный сотрудник, доцент кафедры "Прикладная математика" МГТУ им. Н.Э. Баумана. Лауреат премии Правительства РФ в области науки и техники. Автор более 100 научных работ в области физики горения и взрыва, химической физики и математического моделирования.
A.V. Attetkov (b.1955) graduated from Bauman Moscow Higher Technical School in 1979. Ph. D. (Eng.), senior researcher, assoc. professor of "Applied Mathematics" department of the Bauman Moscow State Technical University. Author of more than 100 publications in the field of physics of combustion and explosion, chemical physics and mathematical simulation.
Елена Владимировна Пилявская окончила МГТУ им. Н.Э.Баумана в 2009 г. Аспирантка кафедры "Прикладная математика" МГТУ им. Н.Э. Баумана. Награждена медалью "'За лучшую научную студенческую работу"' по итогам открытого конкурса на лучшую научную работу студентов по естественным, техническим и гуманитарным наукам в вузах РФ. Автор 10 научных работ в области математического моделирования ударно-волновых процессов в многофазных средах.
Ye.V. Pilyavskaya graduated from Bauman Moscow State Technical University in 2009. Post-graduate of "Applied Mathematics" department of the Bauman Moscow State Technical University. Author of 10 publications in the field of mathematical simulation of shock-wave processes in multiphase materials.
