CC BY
38
УДК 534.222.2 + 536.46 + 661.215.1
ОСЛАБЛЕНИЕ И ПОДАВЛЕНИЕ ВОЛНЫ ГАЗОВОЙ ДЕТОНАЦИИ ОБЛАКОМ ХИМИЧЕСКИ ИНЕРТНЫХ МИКРОЧАСТИЦ
CLOUD CHEMICALLY INERT MICROPARTICLES ATTENUATION AND SUPPRESSION OF GAS DETONATION
WAVE
П. А. Фомин - канд. физ.-мат. наук, старший научный сотрудник ФГБУН «Институт гидродинамики им. М.А. Лаврентьева СО РАН»
А. А. Васильев - д-р физ.-мат. наук, профессор, заведующий лабораторией ФГБУН «Институт гидродинамики им. М.А. Лаврентьева СО РАН»
P. A. Fomin - Cand. Sci. Sciences, Senior Researcher FGBUN "M.A. Lavrentyev Institute of Hydrodynamics SB of RAS "
A. A. Vasiliev - Dr. Sci. Sciences, Professor, Head of the Laboratory FGBUN "M.A. Lavrentyev Institute of Hydrodynamics SB of RAS "
Проанализирован способ ослабления и гашения волны газовой детонации облаком химически инертных твердых частиц. Рассчитаны параметры детонационной волны Чепмена - Жуге в гетерогенных смесях метан - воздух - микрочастицы при различных массовых долях конденсированной фазы и соотношениях между горючим и окислителем. Оценено количественно влияние физико-химических свойств частиц на ослабление детонационной волны. Найдена минимальная концентрация частиц, при которых происходит гашение волны.
Analyzed way to attenuation and extinction of gas detonation wave by a cloud of chemically inert solid particles. The parameters of the Chapman - Jouget detonation wave in heterogeneous mixtures of methane - air - microparticles with various mass fractions of the condensed phase and the ratio between fuel and oxidizer. Quantified the impact of physical and chemical properties of the particles on the weakening of the detonation wave. Found a minimum concentration of particles, at which the wave extinction occurs.
Ключевые слова: ШАХТНЫЕ ВЗРЫВЫ, ГОРЕНИЕ, ДЕТОНАЦИЯ, ПОДАВЛЕНИЕ, ПЫЛЕВЫЕ ЗАСЛОНЫ
Key words: MINE EXPLOSIONS, BURNING, DETONATION, SUPPRESSION, DUST BARRIERS
Введение. Один из способов ослабления и гашения газовой детонации состоит в создании облака химически инертных частиц перед фронтом волны. Потери энергии на разгон и нагрев частиц выше некоторого предела приводят к гашению детонационной волны (ДВ). Один из физически понятных критериев гашения волны состоит в том, что скорость волны должна уменьшиться до такой величины, когда самовоспламенение газа за фронтом волны становится невозможным. В рамках идеализированной модели ДВ с плоским фронтом в настоящей работе рассмотрено также более сильное условие, когда скорость волны уменьшится до звуковой, а температура в пике Неймана будет стремиться к комнатной. Для проведения подобных оценок прежде всего необходимо уметь рассчитывать параметры волны Чепмена - Жуге (Ч-Ж) в смеси газа с частицами, ко-
торые будут забирать на свое ускорение и нагрев энергию химической реакции. Основной вопрос - каков максимально возможный уровень потерь, а, следовательно, и критическая концентрация частиц, при которых обеспечивается успешное подавление волны?
Обозначения:
а - массовая доля частиц;
р, р„ рС - плотность газа, объемная плотность двухфазной смеси и частиц;
P, T, ц - давление, температура и молярная масса газа;
R - универсальная газовая постоянная;
ц , ц , ц - молярные массы газа в
' а ' mm ' max 1
атомарном, предельно диссоциированном и предельно рекомбинированном состояниях;
A, K+ - константы скорости диссоциации и рекомбинации обобщенных продуктов реакции;
научно-технический журнал № 2-2016
ВЕСТНИК
63
в - эффективная температура возбуждения колебательных степеней свободы молекул;
Е - средняя энергия диссоциации обобщенных продуктов реакции;
и, ие, и* - удельные внутренние энергии газа, конденсированной фазы и двухфазной смеси;
Vе, Не - удельная внутренняя энергия и удельная энтальпия конденсированной фазы; у - показатель адиабаты; ив - скорость детонационной волны; с - скорость звука; и - скорость газа. Индексы:
0 - начальное состояние;
1 - состояние в плоскости Чепмена - Жуге; * - двухфазная смесь газа с частицами;
- параметры в пике Неймана. Параметры волны Чепмена - Жуге в смеси газа и микрочастиц
Выпишем уравнения для расчета параметров детонационной волны Ч-Ж в газовой смеси с химически инертными частицами. Хотя в настоящее время хорошо известна сложная неодномерная структура фронта и заметное отличие от классической одномерной модели ДВ, последняя часто используется при оценках интегральных газодинамических параметров ДВ - усредненных волновой и массовой скоростей, давления и температуры, ... Воспользуемся простотой идеализированной модели и примем следующие основные предположения. Волна стационарная, одномерная, без потерь на трение и теплоотвод в стенки трубы. Структура волны соответствует модели Зельдовича - фон Неймана - Деринга. В плоскости Ч-Ж газ находится в состоянии химического равновесия и является идеальным. Конденсированная фаза несжимаема. Парциальное давление и суммарный объем частиц, а также давление насыщенных паров конденсированной фазы пренебрежимо малы по сравнению с газовым давлением. Наибольшую критику вызывает предположение о том, что частицы находятся в механическом и тепловом равновесии с газом (скорость и температура газа и частиц равны). Правомерность последнего предположения для частиц микронного размера достаточно обоснована, а вот для более крупных частиц с размером в сотню микрон вышеуказанные предположения не выполняются, потому результаты расчетов параметров детонационной волны для них могут рассматриваться лишь как оценочные. Предположения о тепловом и механическом равновесии между фазами справедливы для частиц, имеющих относительно малый размер,
когда характерные времена тепловой и механической релаксации между фазами существенно меньше, чем характерные времена изменения параметров двухфазного потока. Например, для частиц размером меньше или порядка десятка микрон скорость и температура частиц успевают сравняться со скоростью и температурой газа уже внутри зоны реакции детонационной волны [1], и, соответственно, предположение о тепловом и механическом равновесии между фазами является допустимым.
Следуя принятым выше предположениям, можно записать:
Р=Р + Рс, Р=(1 - а)Р*, Рс =аР*, Рс =ар/(1 - а) , (1)
и*=(1 - а)и + аиС. (2)
Для описания термодинамических свойств двухфазной системы горючая смесь - инертные частицы в состоянии химического равновесия модифицируем приведенную модель химического равновесия в газах [2-5], которая успешно и широко используется для научных и инженерных расчетов [6-10]. Модель физически обоснована, все константы модели имеют четкий физический смысл, подгоночные коэффициенты отсутствуют. Она проста, обладает высокой точностью, согласована со вторым началом термодинамики и соответствует принципу Ле-Шателье.
В соответствии с [2-5] удельную внутреннюю энергию и молярную массу газа рассчитываем по формулам:
в/т
2
.1) - ]*?+ ехр(0 /Т)-\ я
I
I
Итт
(3)
^ (1--—)гсхр(£/- I).
[1 «и,¡г.
А Г
(1-схр(-0 :'Т)У .
(4)
Из уравнения состояния идеального газа и уравнений (1) - (4) легко получить формулы для расчета скорости звука с* и показателя адиабаты (изоэнтропы) у* рассматриваемой двухфазной смеси в состоянии химического равновесия:
Р р.
0-сОССц цр-ДГ/рц)
¿-у.
м Г ц.
(5)
Производные и, ит, цр, цт являются явными алгебраическими функциями температуры и плотности. Они могут быть легко получены из уравнения состояния идеального газа и уравнений (3) - (4) и здесь не приводятся.
Таким образом, если известны т, р (или р*) и а, то из уравнения состояния идеального газа и (1) - (5), можно рассчитать молярную массу и давление газа, внутреннюю энергию, скорость
64
звука и показатель изоэнтропы гетерогенной смеси химически равновесного газа с твердыми химически инертными частицами. Отметим, что если а равно нулю, формулы (5) трансформируются в соответствующие выражения модели [25].
Таким образом, система алгебраических уравнений для расчета параметров волны может быть записана в виде:
р1и1=р0ив (6)
I
(7)
1 +: Р|"|
1 -н
1-й
14) + — + + -
р, 2 1-СС 2
раи0{иа ^ + -+ (Нс(Та) +
Ро 2 ]-а 2
М-гаад
Р|
I1!
¡и\ й
(8) (9)
А Т.
4 А'
(10)
■■¡|-и | ■/;.!! ^ (11)
Уравнения (6) - (8) являются законами сохранения массы, импульса и энергии, записанными в системе координат, связанной с фронтом волны. Формула (9) представляет собой уравнение состояния идеального газа. Формула (10) описывает состояние химического равновесия. Уравнение (11) выражает условие Чепмена -Жуге.
Величины и(Т1,ц1) и у*(Т1,ц1,а) могут быть рассчитаны по формулам (3) и (5) при Т = Т1, ц = ц. Отметим, что зависимость НС(Т) учитывает теплоемкость частицы, а при наличии плавления - и теплоту фазового перехода конденсированной фазы. Величина ид может быть рассчитана по формуле:
где т и и(Т) - массовая доля и внутренняя энергия - го реагирующего компонента при Т = Тд. Величины и.(Т), НС(Т) и НС(Т) берутся из соответствующих справочников. Более того, для водородно-воздушных смесей величина ид может быть рассчитана по формуле (3) при Т=Тд,
Ц = Цд.
Итак, уравнения (6) - (11), совместно с формулами (3), (5) и формулами для расчета НС(Т), НС(Т) и ио представляют собой замкнутую систему алгебраических уравнений для нахождения р,, и 1, Р, Т1, ц1, ыв. Отметим, что при а = 0 система уравнений (6) - (11) переходит в систему уравнений, которая была успешно использована ранее (например, в [11]) для расчета параметров ДВ в газовых смесях.
Параметры течения за передним ударным
фронтом ДВ («химпик») могут быть рассчитаны по следующей системе алгебраических уравнений:
I 1 П . 1 ..2
Я™ +
' , I1-.: : ^ ' ' ■ (13)
] -а I -а 4 '
Ц8ж=Цд (16)
Уравнения (12) - (14) являются законами
сохранения массы, импульса и энергии, записанными в системе координат, связанной с фронтом волны. Формула (15) представляет собой уравнение состояния идеального газа. Уравнение (16) учитывает отсутствие химической реакции за передним ударным фронтом волны.
Внутренняя энергия газовой фазы может быть рассчитана по Формуле:
¿V = V (Г;|Н), где и.(Т5№) - внутренняя энергия / - го реагирующего компонента при Т = Т^. Величины и/Т!,№) и НС(Т8№) могут быть взяты из соответствующих справочников. Более того, для водородно-воз-душных смесей величина и^ может быть рассчитана по формуле (3) при Т = Т8№, ц = ц^.
Итак, уравнения (12) - (16) совместно с формулами для расчета энтальпий конденсированной фазы и внутренней энергии газа представляют собой замкнутую систему алгебраических уравнений для расчета р8ЙГ,иш, , Т^, ц^.
Влияние термодинамических параметров частиц на ослабление ДВ
Характерные результаты расчетов скорости детонации и параметров потока в плоскости Ч-Ж (скорости - м/с, температуры - в градусах Кельвина, давления - в атмосферах) представлены на рисунке 1 (а-с). Параметры волны за передним ударным фронтом (с индексом SW) показаны на рисунке 1 (4 е). Расчеты проведены при различных массовых долях конденсированной фазы для нескольких материалов частиц с различными термодинамическими свойствами.
Рисунок 1 иллюстрирует влияние термодинамических параметров материала частиц на эффективность подавления ДВ. Потери энергии на нагрев частиц определяются соответствующим увеличением энтальпии конденсированной фазы, которое, в свою очередь, определяется интегралом от теплоемкости материала частиц по температуре и в случае плавления включает
65
2400- 3600
СН4+202 + 1\Г2
; * т
- ^Х V V V V V дигк
2000 : □ □ □ □ □ яЮт 3200
_ Д Д Д Д Д КО
- Ч^чХ 0 0 0 0 0 №/0
1 600-
1200-
2800
2400-
у у у у у А1203 БЮг
д д д д д
wc
Рисунок 1 - Параметры детонации Ч. - Ж. в зависимости от массовой доли твердой фазы для различного вида материалов частиц. Смесь СН4 + 202 + Ы2
66
Ь
110 :
400
ООО
600
2000.0
СН4-а1г М
3000
Т
2600
2200
0.2 0.4
1400
1000
6000.0
СН4-а1Г (з1)
0.2 0.4 0.
а
"0^8 Ь
СН4—а!г (э!)
0.8 а
Рисунок 2 - Параметры идеализированной ДВ в стехиометрической метано-воздушной смеси в за -
висимости от массовой доли частиц БЮ2
а
в себя также скрытую теплоту фазового перехода. Естественно, увеличение энтальпии при нагреве тем выше, чем выше теплоемкость материала частиц и, при наличии плавления, - чем выше теплота фазового перехода. Вследствие большой теплоемкости увеличение энтальпии AlO3 при нагреве до одной и той же температуры больше, чем соответствующее изменение энтальпии SiO2, KCl и WC.
В результате, как видно на рисунке 1 (a-e), параметры ДВ в смесях с SiO2, KCl и WC (то есть скорость волны, давление и температура позади переднего ударного фронта и в плоскости Ч-Ж) выше по сравнению с соответствующими параметрами смесей, содержащих частицы Al2O3. Таким образом, из рассмотренных сортов частиц частицы Al2O3 наиболее эффективны для подавления ДВ. Это, в частности, означает, что суммарная масса частиц Al2O3, которая необходима
для подавления детонации, меньше, чем аналогичная масса частиц SiO2, KCl или WC.
Дальнейшие результаты расчетов приведены для частиц SiO2, поскольку песок является вполне доступным материалом, а его эффективность при подавлении сравнима с Al2O3, хотя и несколькониже.
Ослабление и подавление ДВ в стехио-метрическойметано-воздушнойсмеси
Результаты расчета параметров идеализированной ДВ в стехиометрической метано-воздушной смеси с частицами SiO2 приведены на рисунке 2. Смесь рассматривалась при нормальных начальных условиях.
Температура самовоспламенения метано-воздушной смеси согласно [12] равна 810 K. Из рисунка 2 (d) видно, что минимальная массовая доля частиц в газе, при которых температура смеси за передним ударным фронтом (падаю-
u
D -
1760
I 440
1 120
Т?СЫ4 + (1-77)
air
7)CH4+(l-7?)air
0.04
b
20-
Рт
г]СЯ4 + (1-г]]
air
0.04
1200 U
1040
720-
560-
400-
7?CH4 + (l-?7)air
a = 0
7]
.04
' d
Рисунок 3 - Скорость ДВ и параметры волны в плоскости Ч - Ж в зависимости от начального состава смеси и массовой доли частиц SiO2. Смесь цСЕ4 + (1 - n) Air, стехиометрическому составу соот-
ветствует п = 0.0951, P = 1 атм, T = 293 K
68
c
35 п :
30^ 25: 20 -
?7СН4+(1— ?7)а1г
15-
104
1 600 Тэ¥ =|
1 400 1 200 1 ООО
6000.04
Г] СН4+Ц
V V V V V = 0
□_□_□_□_□ о. 1
Д А А А А 0 . 2
'з
4
5
О.С
0.12
V
Рисунок 4 - Параметры ДВ за передним ударным фронтом («пик Неймана») в зависимости от начального составасмесиимассовойдоличастиц &'ОгСмесь цСЕ4 + (1-ц)Л^, стехиометрическому составу соответствует п = 0.0951, Рд = 1 атм, Т0 = 293 К
щая волна) становится ниже температуры самовоспламенения, равна 0.4. Соответствующая объемная плотность частиц равна 0.67рд, т.е. 0.8 кг/м3. Массовая доля частиц, при которой скорость волны становится звуковой, а температура за передним ударным фронтом близка к комнатной (в силу чего воспламенения газа за фронтом отраженной волны становится заведомо невозможным), равна 0.9. Объемная плотность частиц приэтомравна 9рд,т.е. 10.4 кг/м3.
Следует еще раз заметить, что данные величины получены при условии отсутствия релаксационных процессов по скорости и температуре, что оправдано для достаточно мелких (микронных) частиц. Для более крупных частиц приведенные величины являются оценкой снизу, т.е. для заметного гашения необходимо не менее указанного только что количества частиц.
Параметры ДВ в метано-воздушной смеси различной стехиометрии с добавками микрочастиц песка
На рисунках 3 и 4 представлены расчеты параметров ДВ в метано-воздушной смеси при различных концентрациях горючего и массовых долях конденсированной фазы. Видно, что наиболее трудно подавить ДВ в стехиометрической смеси.
Заключение
В рамках достаточно простого алгоритма проведена оценка минимальной концентрации микронных частиц, обеспечивающей гашение волны газовой детонации в метановых смесях химически инертными частицами. В стехиометрической метано-воздушной смеси минимальная концентрация частиц песка, при которой
температура газа за передним ударным фронтом волны Ч-Ж меньше температуры самовоспламенения, равна 0.8 кг/м3. Концентрация частиц песка для той же смеси, при которой воспламенение газа за фронтом отраженной волны становится заведомо невозможным, равна 10.4 кг/м3.
Из расчета параметров детонации Ч-Ж в метано-воздушной смеси различной стехиометрии с добавками микрочастиц песка следует, что наиболее трудно подавить ДВ в стехиоме-трической смеси.
Проанализировано влияние физико-химических свойств частиц на ослабление детонационной волны. Показано, что ослабление и гашение ДВ более эффективно, если частицы имеют высокую теплоемкость и большую теплоту плавления.
Если частицы имеют относительно большой размер, газ и конденсированная фаза не будут находиться в тепловом и механическом равновесии внутри зоны реакции ДВ [1]. Таким образом, потери энергии на нагрев и ускорение частиц будут меньше, чем в случае мелких частиц, а параметры волны, соответственно, выше. Поэтому эффективность подавления ДВ будет меньше и масса конденсированной фазы, необходимой для гашения волны, должна быть больше, чем соответствующая величина, рассчитанная в рамках предлагаемой в настоящей работе модели. Для увеличения эффективности подавления детонации желательно использовать настолько малые частицы, насколько это возможно.
69
ь
а
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Борисов, А. А. Влияние твердых инертных частиц на детонацию горючей газовой смеси / А. А. Борисов, Б. Е. Гельфанд, С. А. Губин, С. М. Когарко // Физика горения и взрыва. - 1975. - Т. 11. - № 6.
- С. 909-914.
2. Николаев, Ю. А. Модель кинетики химических реакций при высоких температурах / Ю. А. Николаев // Физика горения и взрыва. - 1978. - Т. 14. - № 4. - С. 73-76.
3. Николаев, Ю. А. О расчете равновесных течений химически реагирующих газов / Ю. А. Николаев, П. А. Фомин // Физика горения и взрыва. - 1982. - Т. 18. - № 1. - С. 66-72.
4. Николаев, Ю. А. Приближенное уравнение кинетики в гетерогенных системах типа газ - конденсированная фаза / Ю. А. Николаев, П. А. Фомин // Физика горения и взрыва. - 1983. - Т. 19. - № 6. - С. 49-58.
5. Николаев, Ю.А. Согласование моделей химических реакций в газах со вторым началом термодинамики / Ю. А. Николаев, Д. В. Зак // Физика горения и взрыва. - 1988. - Т. 24. - № 4. - С. 87-90.
6. Васильев, А. А. Экспериментальное исследование и численное моделирование расширяющейся многофронтовой детонационной волны / А. А. Васильев, А. В. Троцюк // Физика горения и взрыва.
- 2003. - Т. 39. - № 1. - С. 92-103.
7. Николаев, Ю. А. Модель стационарной гетерогенной детонации в газокапельной среде / Ю. А. Николаев, П. А. Фомин // Физика горения и взрыва. - 1984. - Т. 20. - № 4. - С. 97-105.
8. Gavrilenko, T. P. Acceleration of solid particles by gaseous detonation products / T. P. Gavrilenko, V. V. Grigoriev, S. A. Zhdan [et al.] // Combustion and Flame. - 1986. - Vol. 66. - С. 121-128.
9. Ждан, С. А. Величина реактивного импульса от взрыва газовой смеси в полуограниченном пространстве / С. А. Ждан, В. В. Митрофанов, А. И. Сычев // Физика горения и взрыва. - 1994. - Т. 30. - № 5. - С. 90-97.
10. Fomin, P. A. Modeling of detonation processes in chemically active bubble systems at normal and elevated initial pressures / P. A. Fomin, K. S. Mitropetros, H. Hieronymus // Journal of Loss Prevention in the Process Industries. - 2003. - Vol. 16. - N 4. - Pp. 323-331.
11. Mitropetros, K. Behavior of the surface of a bubbly liquid after detonation wave impact / K. Mitropetros, P. A. Fomin, H. Hieronymus // Experiments in Fluids. - 2005. - Vol. 40. - I. 3. - Pp. 431-441.
12. Пожаро-взрывоопасность веществ и материалов и средства их тушения. Справочное издание: В 2 кн. / А. Н. Баратов, А. Я. Корольченко, Г. Н. Кравчук [и др.]. - М.: Химия, 1990.
REFERENCES
1. Borisov, A. A., Gel'fand, B. E., Gubin, S. A., & Kogarko, S. M. (1975). Vlijanie tverdyh inertnyh chastic na detonaciju gorjuchej gazovoj smesi [Effect of inert solid particles on the combustible gas mixture detonation]. Fizika Gorenija I Vzryva, 11(6), 909-914.
2. Nikolaev, J. A. (1978). Model' kinetiki himicheskih reakcij pri vysokih temperaturah [Model of kinetics of chemical reactions at high temperatures]. Fizika Gorenija I Vzryva, 14(4), 73-76.
3. Nikolaev, J. A., & Fomin, P. A. (1982). O raschete ravnovesnyh techenij himicheski reagirujushhih gazov [About computation of the equilibrium flows of chemically reacting gas]. Fizika Gorenija I Vzryva, 18(1), 66-72.
4. Nikolaev, J. A., & Fomin, P. A. (1983). Priblizhennoe uravnenie kinetiki v geterogennyh sistemah tipa gaz - kondensirovannaja faza [An approximate equation of kinetics in heterogeneous systems such as gas -condensed phase]. Fizika Gorenija I Vzryva, 19(6), 49-58.
5. Nikolaev, J. A., & Zak, D. V. (1988). Soglasovanie modelej himicheskih reakcij v gazah so vtorym nachalom termodinamiki [Coordination of models of chemical reactions in gases with the second law of thermodynamics]. Fizika Gorenija I Vzryva, 24(4), 87-90.
6. Vasil'ev, A. A., & Trocjuk, A. V. (2003). Jeksperimental'noe issledovanie i chislennoe modelirovanie rasshirjajushhejsja mnogofrontovoj detonacionnoj volny [Experimental study and numerical simulation of expanding the multi-detonation wave]. Fizika Gorenija I Vzryva, 39(1), 92-103.
7. Nikolaev, J. A., & Fomin, P. A. (1984). Model' stacionarnoj geterogennoj detonacii v gazokapel'noj srede [Model of stationary heterogeneous detonation in gas-drop environment]. Fizika Gorenija I Vzryva, 20(4), 97-105.
8. Gavrilenko, T., Grigoriev, V., Zhdan, S., Nikolaev, Y., Boiko, V., & Papyrin, A. (1986). Acceleration of solid particles by gaseous detonation products. Combustion and Flame, 66(2), 121-128. doi:10.1016/0010-2180(86)90084-2
научно-технический журнал № 2-2016
ВЕСТНИК
9. Zhdan, S. A., Mitrofanov, V. V., & Sychev, A. I. (1994). Velichina reaktivnogo impul'sa ot vzryva gazovoj smesi v poluogranichennom prostranstve [The size of the jet pulse from the explosion of the gas mixture in the semi-infinite space]. Fizika Gorenija I Vzryva, 30(5), 90-97.
10. Fomin, P., Mitropetros, K., & Hieronymus, H. (2003). Modeling of detonation processes in chemically active bubble systems at normal and elevated initial pressures. Journal of Loss Prevention in the Process Industries, 16(4), 323-331. doi:10.1016/s0950-4230(03)00018-4
11. Mitropetros, K., Fomin, P. A., & Hieronymus, H. (2005). Behavior of the surface of a bubbly liquid after detonation wave impact. Exp Fluids Experiments in Fluids, 40(3), 431-441. doi:10.1007/s00348-005-0081-x
12. Baratov, A. N., Korol'chenko, A. J., Kravchuk, G. N., & At al. (1990). Pozharovzryvoopasnost substances and materials and their extinguishing agents. Reference book: In 2 kn. [Fire and explosion hazard of substances and materials and their extinguishing agents. Reference book: In 2 books]. Moscow: Himija.
