Исследование параметров регистрации теплового потока при горении и детонации в канале газовой смеси

Размещено на http: //www. allbest. ru/

Исследование параметров регистрации теплового потока при горении и детонации в канале газовой смеси

Пай Владимир Васильевич, Гулевич Максим Александрович

Яковлев Игорь Валентинович, Пинаев Александр Владимирович

Васильев Анатолий Александрович, Игнатенко Александр Георгиевич

Саяпин Виталий Викторович, Доманов Виктор Петрович

Трубицына Дарья Анатольевна, Батраков Дмитрий Николаевич

Аннотация

В.В. Пай (д-р физ.-мат. наук, ведущий научный сотрудник Института гидродинамики им. М.А. Лаврентьева СО РАН)

М.А. Гулевич (младший научный сотрудник Института гидродинамики им. М.А. Лаврентьева СО РАН)

И.В. Яковлев (д-р техн. наук, заведующий лабораторией Института гидродинамики им. М.А. Лаврентьева СО РАН)

А.В. Пинаев (д-р физ.-мат. наук, ведущий научный сотрудник Института гидродинамики им. М.А. Лаврентьева СО РАН)

А.А. Васильев (д-р физ.-мат. наук, директор Института гидродинамики им. М.А. Лаврентьева СО РАН)

А.Г. Игнатенко (ведущий инженер ОАО «Новосибирский механический завод «Искра»)

В.В. Саяпин (заведующий лабораторией ОАО «Новосибирский институт программных систем»)

В.П. Доманов (канд. техн. наук, заведующий лабораторией ОАО «НЦ ВостНИИ»)

Д.А. Трубицына (директор ООО «ЦОТ Горный»)

Д.Н. Батраков (научный сотрудник ОАО «НЦ ВостНИИ»)

Представлен экспериментальный метод исследования параметров регистрации тепловых потоков при горении или детонации в канале газовой смеси термопарными датчиками. Подробно рассматриваются тепловые и электромагнитные процессы, происходящие в планарных термопарах при их импульсном нагреве продуктами реакции газовых смесей. Приводятся экспериментальные данные, полученные с помощью таких термопар при детонации смеси 2O2+N2+CH4 в металлической трубе диаметром 70 мм при исходных давлениях - 0,2; 0,4 атм, а также при распространении быстрого горения в трубе диаметром 1800 мм в метано-воздушной смеси при атмосферном давлении с объемной концентрацией метана 9,3 - 9,8 %. Выполнены численные расчеты тепловых потоков из реагирующей газовой смеси на стенки канала. Термобатареи, состоящие из таких планарных термопар, могут быть применены в шахтах, опасных по газу и пыли, в качестве датчиков возгорания и взрыва. Характерное время формирования термобатареей электрического импульса составляет примерно 1 мс, что позволяет использовать его в качестве стартового сигнала для запуска систем взрыво-пожаротушения, имеющих время срабатывания примерно 100 мс.

Работа поддержана грантом РФФИ № 13-03-00791 и программой Президиума РАН № 2.8, 2012-2014.

Ключевые слова: ТЕПЛОВОЙ ПОТОК, ГАЗОВАЯ СМЕСЬ, МЕТАН, ТЕРМОБАТАРЕЯ, СТЕНКИ КАНАЛА, ЭЛЕКТРИЧЕСКИЙ ИМПУЛЬС, СИСТЕМА ПОЖАРОТУШЕНИЯ

тепловой детонация газовый датчик

Annotatіon

RESEARCH OF HEAT STREAM REGISTRATION PARAMETERS DURING BURNING AND DETONATION IN THE CHANNEL OF GAS MIXTURE

V.V. Pai, M.A. Gulevich, I.V. Yakovlev, A.V. Pinaev, A.A. Vasiliev, A.G. Ignatenko, V.V. Sayapin, V.P. Domanov, D.A. Trubitsina, D.N. Batrakov

Experimental method of research of heat stream on the channel walls with thermocouple sensors during burning and detonation of gas mixture in it is presented. The heat and electromagnetic processes happening in planar thermocouples at their impulse heating with the reaction products of gas mixtures are reviewed in details.

The experimental data obtained by such thermocouples in the detonation process of 2O2+N2+CH4 mixture in 70 mm metal pipe with the initial pressure 0,2; 0,4 atm, and also during fast spread of burning in the 1800 mm diameter pipe in methane-air mixture at atmospheric pressure with volumetric methane concentration about 9,3 - 9,8% are presented. Numerical calculations of heat streams from reacting gas mixture on the channel walls are done. Thermal batteries consisting of such planar thermocouples can be used as combustion and explosion sensors in the mines dangerous for gas and dust. Characteristic time of electric impulse formation by thermal battery is about 1 ms, which allows to use it as the start signal to initiate the system of fire-explosion extinguishing which usually take 100 ms to start working.

The work is supported by RFFI grant No. 13-03-00791 and by RASc Presidium program No. 2.8, 2012-2014.

Key words: HEAT STREAM, GAS MIXTURE, METHANE, THERMAL BATTERY, CHANNEL WALLS, ELECTRIC IMPULSE, FIRE EXTINGUISHING SYSTEM

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

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

Для понимания сути предлагаемого метода рассмотрим подробно конструкцию элемента планарной термобатареи, схематически изображенной на рисунке 1.

Рисунок 1 Элемент термобатареи 1 - теплоизолирующая подложка; 2 - константановая лента; 3 - слой меди; 4 - защитный слой

На поверхности теплоизолирующей подложки (1) расположена константановая лента (2), частично покрытая электролитическим способом слоем меди (3). На эту ленту, состоящую из чередующихся монометаллических (константановых) и биметаллических (медь-константановых) частей, нанесен защитный слой (4) так, чтобы каждый второй стык монометалла с биметаллом был изолирован от контакта с горячим газом. Толщины константановой ленты - 5 мкм, слоя меди - 2,75 мкм, защитного слоя - 50 мкм, подложки - 1 мм, ширина ленты - 350 мкм. Часть ленты, состоящая из одной монометаллической и биметаллической частей, представляет собой одну термопару. Длины монометаллической и биметаллической частей, а также участков термопары, покрытых защитным слоем, были одинаковыми и равными 500 мкм.

1 Теория метода

Пусть на открытую часть термопары в результате ее теплообмена с контактирующим горячим газом идет поток тепла с плотностью , а , , - коэффициенты теплопроводности; , , - удельные теплоемкости; , , - плотности меди, константана, подложки и защитного слоя, соответственно [1]; , - толщина медного слоя и константановой ленты. Найдем распределение температуры в термопаре в произвольный момент времени. Вначале заметим, что характерное время выравнивания температуры по толщине ленты как в ее монометаллической, так и биметаллической частях составляет: мкс, в то время как выравнивание температуры по длине происходит в течение мкс, где - длина термопары. Поэтому при характерных временах изменения величины плотности теплового потока , которые реализуются в описываемых ниже экспериментах, можно считать температуру в поперечном сечении термопары одинаковой, функцией лишь одной координаты вдоль ленты. Другими словами, достаточно искать лишь зависимость , где - время. Необходимо также учитывать теплообмен термопары как с подложкой, так и с защитным слоем, имеющими практически одинаковые теплофизические характеристики. Это обусловлено тем, что хотя коэффициенты их теплопроводности примерно в 103 раз меньше, чем у металлов термопары, но поверхность теплообмена больше поперечного сечения ленты примерно в 102 раз.

Для биметаллического участка термопары с защитным слоем уравнение теплопроводности будет:

, (1)

где - плотность теплового потока от термопары как в подложку, так и в защитный слой.

Для незащищенного биметаллического участка:

, (2)

где - плотность теплового потока от горячего газа к термопаре.

Для защищенного монометаллического участка:

. (3)

Для незащищенного монометаллического участка:

. (4)

На стыке монометаллической и биметаллической частей при выполняются условия сшивки:

, .

Поскольку распределение температуры в каждой из термопар одинаково, то условие цикличности будет следующим:

, .

Начальное условие .

Вычисление плотности теплого потока от термопары в подложку и защитный слой проводилось из следующих соображений. Так как глубина прогрева подложки и защитного слоя за время ~ 1 мс составляет мкм, что существенно меньше длины термопары, то в уравнении для расчета теплопроводности в подложке и защитном слое, пренебрегая членами, содержащими производные по , имеем:

.

При этом граничные условия (считаем толщину подложки и толщину защитного слоя бесконечными, так как характерное время их прогрева 25 мс): ; , начальное условие: . Решением этой задачи для произвольной является функция [2]:

.

Из этой формулы определяем

.

После подстановки в уравнения (1), (2), (3), (4) численно может быть найдено распределение температуры в термопаре в зависимости от времени при произвольно заданной функции .

Теперь рассмотрим зависимость термо-ЭДС такой термопары как функцию времени.

В результате неоднородного нагрева термопары в ней возникают электрическое поле и замкнутые электрические токи, в основном сосредоточенные в области существенных изменений температуры границы контакта металлов [3]. Запишем дифференциальный закон Ома с учетом термо-ЭДС для каждого из металлов биметаллической части термопары [4]:

,

где ; - векторы напряженности электрического поля и плотности тока, соответственно;

, - термоэлектрический коэффициент и электрическая проводимость соответствующего металла [1].

Из условия непрерывности тангенциальной составляющей электрического поля к границе контакта металлов (индексы 1 и 2 относятся к меди и константану, соответственно) имеем:

; . (5)

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

; .

Тогда, поскольку полный ток, протекающий через поперечное сечение термопары, равен нулю: , из формулы (5) следует:

.

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

.

Для монометаллической части термопары:

; .

Суммируя, получим выражение для напряжения на одной термопаре:

, т.к. .

Наконец, из найденного выше решения о распределении температуры в термопаре, можно вычислить зависимость термо-ЭДС от времени для произвольно заданной функции . При незначительной разности температур незащищенных монометаллической и биметаллической частей термопары (~ 10-100 °C) в сравнении с температурой газа (~ 2000-3000 °C) [5] можно считать, что плотность теплового потока не зависит от , так как его величина является монотонно растущей функцией разности температур газа и поверхности термопары, то есть .

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

. (6)

Это интегральное уравнение Вольтерра 2-го рода, корректное при определении , ядро которого находится численным расчетом с произвольной точностью. Функция при экспериментально определенной вычислялась путем алгебраизации уравнения (6). Для расчета зависимости необходимо выполнение некорректной процедуры - дифференцирования функции по времени. Это возможно без существенной потери точности, при наличии априорной информации о гладкости искомого решении. Ниже при рассмотрении результатов экспериментов будет обсуждаться такая априорная информация относительно характера поведения .

2 Схемы и результаты экспериментов

Учитывая значительный диапазон изменения скоростей распространения волн горения и детонации, а также принципиальное их различие по структуре фронта и величинам газодинамических параметров течения, тестирование термобатареи проводилось как для детонационного режима, так и для режима горения. В лабораторных экспериментах с детонационным режимом использовалась металлическая труба длиной 7 м и диаметром 70 мм, в которую для измерения тепловых потоков от продуктов детонации помещалась термобатарея площадью 1,25 см2, содержащая термопар. При этом плоскость термобатареи была перпендикулярна диаметральному сечению трубы (рисунок 2). Эксперименты с детонацией проведены на смеси CH4+2O2+N2. Выбор данной смеси и размеров детонационной трубы определялся условием сохранения геометрического подобия детонационных явлений в метано-воздушных смесях в условиях реальной шахты и в метановых смесях в лабораторных условиях. Параметром геометрического подобия служит отношение между характерным размером канала (диаметром трубы или шахтного штрека) и характерным размером зоны химической реакции в детонационной волне (размером детонационной ячейки).

После воспламенения смеси электрическим разрядом и формирования стационарной детонационной волны измерялись ее скорость и напряжение на термобатарее в зависимости от времени

.

Измеренная пьезодатчиками скорость детонационного процесса оказалась равной 2220 м/с, что близко к расчетной скорости детонации [5].

Рисунок 2 Схема эксперимента1 - термобатарея; 2 - газовая смесь CH4+2O2+N2; 3 - детонационный фронт

На рисунке 3 а приведены осциллограммы зависимостей для начальных давлений смеси 0,2 атм и 0,4 атм. Здесь же на рисунке 3 б изображены рассчитанные по уравнению (6) зависимости . Эти зависимости, как видно из сравнения, хорошо аппроксимируются функциями вида с показателем экспоненты мс-1. Вид функциональной зависимости выбирался из следующих соображений. Горячий газ с удельной тепловой энергией, пропорциональной температуре, двигаясь вдоль металлической трубы, непрерывно остывает в процессе теплообмена со стенками, с плотностью теплового потока . Тогда , откуда следует, что

; , .

Значения и были найдены из условия минимума , где и - рассчитанные и экспериментально определенные величины напряжения на термопаре в момент времени , соответственно.

На рисунке 3 в изображены зависимости плотностей тепловых потоков от времени.

Рисунок 3 Результаты экспериментов и расчетов при начальных давлениях газовой смеси 0,2 атм (I) и 0,4 атм (II) а - экспериментальные (1) и рассчитанные (2) зависимости напряжения на термопаре от времени; б - рассчитанные по экспериментальным данным зависимости Q 0(t) (1) и Q(t) (2); в - зависимости плотности теплового потока от времени

Следующая серия крупномасштабных экспериментов, приближенных к шахтным условиям, проводилась в металлической трубе диаметром 1,8 м, длиной 15 м, частично заполненной метано-воздушной смесью, с объемным содержанием метана от 9,3 до 9,8 %. Схема эксперимента изображена на рисунке 4 с указанием мест размещения трех одинаковых термобатарей, их плоскости были ориентированы вдоль стенок трубы. Датчики №1 и 2 регистрировали волну в смеси, датчик №3 - вне смеси, в воздухе. Для инициирования газовой смеси применялся заряд детонита массой 100 г.

Рисунок 4 Схема эксперимента1 - термобатареи; 2 - метано-воздушная смесь; 3 - заряд детонита

На рисунке 5 (верхний график) приводятся осциллограммы зависимости измеренных напряжений от времени на термопарах в трех одинаковых экспериментах (кривые 1, 2, 3). Для каждого из экспериментов рассчитывалась зависимость , показанная на рисунке 5 (нижний график).

Рисунок 5 Осциллограммы зависимости измеренных напряжений от времени на термопарах в трех экспериментах верхний график - зависимости напряжения на термобатареях от времени (1, 2, 3 - номера экспериментов); нижний график - зависимости количества тепла от времени

Как видно из осциллограмм, скорость распространения процесса в метано-воздушной смеси составляла 480 м/с, причем головная часть представляет собой ударно-сжатую метано-воздушную смесь, а за ней с отставанием 5-7 мс движется фронт горения этой смеси. Тепловые потоки в головной части и в зоне горения смеси отличаются приблизительно в 5-8 раз. Такой процесс является существенно неодномерным и нестационарным, чем объясняется плохая повторяемость регистрируемой зависимости . В этом случае имеет смысл указать лишь средние значения плотности теплового потока в областях сжатой смеси и за фронтом горения. Для головной части Вт/см2, для области горения смеси Вт/см2.

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

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

Библиографический список

1 Таблицы физических величин / под ред. акад. И.К. Кикоина. - М.: Атомиздат, 1976. - 1008 с.

2 Карслоу, Г. Теплопроводность твердых тел / Г. Карслоу, Д. Егер. - М.: Наука, 1964. - 488 с.

3 Экспериментально-аналитические методы в задачах динамического нагружения материалов / Г.Е.Кузьмин, В.В.Пай, И.В. Яковлев. - Новосибирск: Изд-во СО РАН, 2002. - 312 с.

4 Ландау, Л.Д. Электродинамика сплошных сред / Л.Д. Ландау, Е.М. Лифшиц. - М.: Наука, 1982. - 624 с.

5 Characteristics of combustion and detonation of methane-coal mixtures / A.A.Vasil'ev, A.A.Trubitsyn, A.V.Pinaev // 9-th International Symposium on Hazards, Prevention, and Mitigation of Industrial Explosions, Krakow, Poland, July 22-27, 2012. 9th ISHPMIE Proceedings (USB, ish063.pdf) - 9 p.

Размещено на Allbest.ru