ГИДРОДИНАМИКА, ТЕПЛО-И МАССООБМЕННЫЕ ПРОЦЕССЫ, ЭНЕРГЕТИКА
УДК 662.237.5; 544.032
Л. Х. Бадретдинова, О. В. Костиным, Е. Б. Смирнов, А. В. Станкевич, К. А. Тем, Б. П. Толочко
ИЗОТЕРМИЧЕСКОЕ СЖАТИЕ 1,3,5ТРИАМИНО 2,4,6ТРИНИТРОБЕНЗОЛА
Ключевые слова: энергонасыщенный материал, изотермическое сжатие, алмазные наковальни, ТАТБ.
В работе предложен метод определения кривой изотермического сжатия на примере энергонасыщенного вещества 1,3,5триамино 2,4,6тринитробензол (ТАТБ).
Key words: explosives, isothermal compression, diamond anvils, TATB.
In this work a method for determining the isothermal compression curve on the example of explosives TATB is present.
Энергонасыщенный материал (ЭМ) ТАТБ чрезвычайно привлекательно для научных исследований, благодаря своей сложной кристаллической структуре [1-2]. Обладая триклинной сингонией, кристаллы ТАТБ характеризуются высокой анизотропией и низкой симметрией, представляющей значительные трудности для рентгеноструктурного анализа. Вместе с тем, ТАТБ обладает уникально низкой чувствительности к внешним тепловым и механическим воздействиям [3]. Высокая стойкость по сравнению с большинством известных ЭМ позволяет расширить диапазон изменения термодинамических параметров (давления, температуры и т.д.) при проведении экспериментальных исследований, необходимых для построения уравнения состояния.
Параметры уравнения состояния
определялись на основании рентгеноструктурных исследований ТАТБ в изотермических условиях. Кривая изотермического сжатия до давлений ~6,5 ГПа при температуре T0=293 К была получена на ускорительном комплексе Института Ядерной Физики (ИЯФ) Сибирского Отделения РАН с использованием метода порошковой дифракции при сжатии ТАТБ в ячейке с алмазными наковальнями DAC (diamond anvil cell). Схема используемой в данной работе ячейки с алмазными наковальнями, конструкции Меррилла-Бассета [4], приведена на рисунке 1.
ElectromagnetJc Radiation
— Ruby
t_> Sample
Backing Plate
Рис. 1 - Схема ячейки с алмазными наковальнями
В этих устройствах два алмаза ювелирного качества, ограненные специальным способом, сжимаются с помощью несложного рычажно-винтового устройства. Исследуемый образец вместе со смесью метилового и этилового спиртов и кристаллом рубина помещается в отверстие в прокладке (гаскете) из прочного пластичного материала. Снаряженная гаскета устанавливается непосредственно между рабочими плоскостями алмазов. Подробное описание работы таких камер дано в обзорах [5-7]. Смесь спиртов обеспечивает гидростатичность давления, а по смещению линии Я-люминесценции рубина измеряется давление [8], приложенное к исследуемому образцу. Для изменения давления и его измерения на оптическом спектрометре алмазные наковальни снимали с синхротронного пучка и после этих процедур снова на нем юстировали.
Исследования проводились на источнике синхротронного излучения (СИ) ускорительного комплекса ВЭПП-3, расположенном в ИЯФ СО РАН. Источником СИ на данном ускорителе является излучение релятивистских электронов с энергией 2 ГэВ из вигглера с полем 2 Тл. Синхротронное излучение вигглера имеет непрерывный спектр с критической длиной волны Хс= 2.33 А и интегральной мощностью 103 Вт. Схема проведения экспериментов приведена на рис.2. Пучок рентгеновских лучей по вакуумированному каналу, через коллиматор (размер щели не превышал 100 мкм), подводится к алмазной камере высокого давления. С помощью коллиматора можно менять размер падающего пучка и избавляться, таким образом, от рассеяния на материале металлической прокладки (гаскете). Регистрация дифракционных сигналов от сжатого образца ТАТБ осуществлялась при помощи детектирующей системы на основе запоминающего экрана МАЯ345 фирмы МагЯ^еагсИ. Энергия рентгеновских квантов составляла E=33,7 кэВ
[9, 10]. Получаемая картина двумерной дифракции
интегрированием по азимутальному углу переводилась в одномерный спектр "Интенсивность-29", который анализировался в дальнейшем. Полученные таким образом дифрактограммы приведены на рис 2.
Для сравнения на графике приведены данные, полученные другими исследователями [11-13]. Данные по изотермическому сжатию аппроксимировались уравнением полученным из (2) в предположении, что давление складывается из потенциальной и тепловой составляющих
IOC-
Э
Ш 50
й с Ш с
--М-
3L
_л.
jL
P (V,T) = Px (V) + Pt (V,T):
4 S e 7 8 9 10
angle, 20
Рис. 2 - Дифрактограммы ТАТБ при сжатии в алмазных наковальнях температуре T=293K
При атмосферном давлении
регистрировалось до десяти дифракционных пиков, что вполне достаточно для анализа триклинной системы. Однако, с увеличением давления до ~4 ГПа, количество регистрируемых рефлексов сокращалось до четырех. Для расчета параметров элементарной ячейки ТАТБ, необходимо регистрировать минимум шесть дифракционных пиков. На основе дифрактограмм методами рентгеноструктурного анализа, учитывающего реальную
кристаллографическую модель ТАТБ, были определены параметры кристаллической структуры и элементарный объем молекулярного кристалла ТАТБ. Используя данную модель, удалось рассчитать объем элементарной ячейки на основании дифрактограмм с малым (менее 6) количеством рефлексов.
Взаимосвязь гидростатического давления
Px (о)в и относительного объема кристалла ТАТБ приведена на рис.3.
1 2 bjjir' Ф " 1_ и 4 * н^-н. / н-Я(-н ' : 0^ -/ н н
1,05 1, ■ ■ ■ ■ , ф 0 1,15
3 Stevens -В Olinger . ф Pastine О Данная работа ^ ^ ^ аппроксимация •
Рис. 3 - Изотерма ТАТБ при температуре T=293K
\до )T ц
9в+3ТЦ-
8 It
(1)
Тепловой член в уравнении (1) учитывает начальную температуру проведения
Потенциальная составляющая задавалась в виде потенциала
экспериментов. давления Px(6) Борна-Майера:
3PokCok
Px (О) =
b - 3s-1
( 2 О3 exp
(
b
1-5
\
-5
s+1
(2)
а температура Дебая 6(6) - выражением:
Г dPx (5)
t-2
Ф) "
t -2
d5
5-Px (5)k
л-2
5 ~lpokClk
(3)
V У
Функция Грюнайзена у(6) описывалась приближением Ландау-Слейтера. Учитывая, что силами притяжения в молекулярных кристаллах являются силы Ван-дер-Ваальса, значение в в выражении (2) фиксировалось и принималось равным 5 = 2. В ходе аппроксимации методом нелинейной регрессии были найдены коэффициенты Ь и Сок. Полученное уравнение изотермы (1) хорошо согласуется с данными работы [11], полученными для ТАТБ в более широком (до 13 ГПа) диапазоне изменения давления. Полученное в результате аппроксимации значение Сок - скорости звука при Т=0 К оказалось завышенным по сравнению с данными работ [13, 14].
Таким образом, получено уравнение изотермического сжатия взрывчатого вещества 1,3,5триамино 2,4,6тринитробензол, хорошо описывающее статические эксперименты по изотермическому сжатию. В ходе аппроксимации экспериментальных данных по изотермическому сжатию определены все параметры уравнения состояния вида:
Р (У,Т) = Рх (V) + Ра (V)* Рт(У,7),
где Рх(Ю - потенциальная (холодная) часть свободной энергии, которая зависит только от объема V; Ро- свободная энергия нулевых колебаний кристалла при 7=0; РТ(^Т) - тепловая (квазигармоническая) часть свободной энергии, зависящая от объема и температуры.
Можно ожидать, что использование полученного уравнения состояния (1) позволит повысить точность описания термодинамических
V /V
параметров непрореагировавшего энергонасыщенного
материала ТАТБ при численном моделировании
ударно-волновых и детонационных процессов.
Литература
1. Badretdinova L.Ch., Ten K.A., Pruuel E.R., Lukyanchikov L.A., Tolochko B.P., Sharafutdinov M.R., Shmakov A.N., Muzyrya A.K., Kostitsyn O.V., Smirnov E.B., Kuper K.E. Structure of condensed heterogeneous explosives. Conference Interaction of Intense Energy Fluxes with Matter. IIEFM-2013, March 1-6, 2013. Elbrus, Kabardino-Balkaria, Russia.
2. Бадретдинова Л.Х. Квантово-химические модели про-гнорзирования чувствительности к удару полинитросо-единений / Л.Х. Бадретдинова, А.Н. Анисимов, Д.А. Ха-диева, В.Я. Базотов // Вестник КТУ - Казань, 2014. - Т.17, № 7, С. 26-27.
3. Бадретдинова Л.Х. Прогноз параметров чувствительности энергонасыщенных веществ с использованием кластерных моделей / Л.Х. Бадретдинова, А.Н. Анисимов, Д.А. Ха-диева, В.Я. Базотов // Вестник КТУ - Казань, 2014. - Т.17, № 8, С. 25-27.
4. L. Merrill & W.A. Basset, Miniature Diamond Anvil Pressure Cell for Single Crystal X-ray Diffraction Studies, Review of Scientific Instruments 1974, 45, 290—294
5. Eremets M. High Pressure Experimental Methods - Oxford University Press, 1996.
6. Jayaraman A. Diamond anvil cell and high-pressure physical investigations // Reviews of Modern Physics. -1983. -V.55. -N1. P. 65-108.
7. Bohler R., Musshoff H. G., Ditz R., Aquilanti G., Trapananti A. Portable laser-heating stand for synchrotron applications. // Rev. Sci. Instrum. -2009. - V.80. -P.045103-045106.
8. Forman R.A., Piermarini G.J., Barnett J.D. & Block S. Pressure Measurement Made by the Utilization of Ruby Sharp-Line Luminescence // Science. -1972. -V.176. -P.284—285.
9. Тен К.А., Прууэл Э.Р., Лукьянчиков Л.А., Толочко Б.П., Шарафутдинов М.Р., Шмаков А.Н., Аминов Ю.А., Музыря А.К., Костицын О.В., Смирнов Е.Б. Исследование ТАТБ дифракционными методиками. XI Забахинские научные чтения. Снежинск. 2012.
10. БадретдиноваЛ.Х., Тен К.А., Прууэл Э.Р., Лукьянчиков Л.А., Толочко Б.П., Шарафутдинов М.Р., Шмаков А.Н., Музыря А.К., Костицын О.В., Смирнов Е.Б., Купер К.Э. Исследование структуры конденсированных гетерогенных ВВ. Всероссийская конференция "Взрыв в физическом эксперименте" Новосибирск, Академгородок, 16-20 сентября 2013.
11. Stevens L. L., Velisavljevic N., Hooks D. E., Dattelbaum D.M.. Hydrostatic Compression Curve for Triamino-Trinitrobenzene Determined to 13.0 GPa with Powder X-Ray Diffraction. Propellants, Explosives, Pyrotechnics 33, No. 4, 2008.
12. Pastine D. J., Bernecker R. R., P,v,E,T Equation of State for 1,3,5-Triamino-2,4,6-Trinitrobenzene, J. Appl. Phys. 1974, 45.
13 Olinger B., Cady H., The Hydrostatic Compression of Explosives and Detonation Products to 10 GPa (100 kbars) and Their Calculated Shock Compression: Results for PETN, TATB, CO2 and H2O, 6th Symposium (International) on Detonation, Coronado, California, USA, 1976.
14. Grebyonkin K.F., Gerebtsov A.L., Popova V.V., Taranik M.V. P,V,E,T Equation of State for TATB Based Explosives. Shock Compression of Condensed Matter - 2003
© Л. Х. Бадретдинова - к.т.н., асс. кафедры ТТХВ КНИТУ, salamandra_1985@mail.ru, О.В. Костицын - к.т.н., нач. отдела РФЯЦ ВНИИТФ; Е.Б. Смирнов - нач. группы РФЯЦ ВНИИТФ, ewgeny_smirnov@mail.ru; А. В. Станкевич - асп. кафедры ТТХВ КНИТУ, К. А. Тен - к.т.н., с.н.с. ИГиЛ СО РАН, Б. П. Толочко - д.х.н., рук. лаб. ИХТТМ СО РАН.
© L. Kh. Badretdinova - associate professor, KNRTU, salamandra_1985@mail.ru, O. V. Kostitsyn - Ph.D. of technical science; E.B. Smirnov - Ph.D. of technical science, ewgeny_smirnov@mail.ru; A.V. Stankevich - undergraduate KNRTU; K. A. Ten - Ph.D. of technical science, Russian of the academy of sciences; B. P. Tolochko - Doctor of Chemical Science, Russian of the academy of sciences.