Studying the influence of relative motion of suspended inert particles on the rate of the gas mixture combustion front.pdf Газодисперсные системы, к которым относятся различные газовзвеси инертных частиц, а также частиц органических горючих, часто используются в химической, пищевой промышленности, инертные порошки используются для пожаротушения [1-5]. Исследование и моделирование процессов, происходящих при воспламенении и горении газовзвесей, изучение качественных и количественных характеристик - скоростей горения, температур, плотностей и т.д., является необходимым для эффективного развития промышленных технологий. В данной работе представлены результаты теоретического исследования горения газовой смеси со взвешенными в ней мелкими инертными частицами с учетом теплового расширения газа. Проведено сравнение результатов с результатами вычислений по модели, в которой не учитывается тепловое расширение газовой фазы и движение частиц. Физико-математическая постановка задачи Горючая газовая смесь с взвесью инертной пыли занимает область 0 < х < да. Считается, что реакция в газовой фазе проходит по второму порядку. Частицы одного размера, сферической формы, с радиусом rk, в начальный момент времени распределены равномерно в газовой смеси. В координате х = 0 размещена горячая непроницаемая стенка с температурой равной адиабатической температуре сгорания газа в изобарических условиях. Такое граничное условие моделирует зажигание у «закрытого торца» [6]. Горение протекает в ламинарном режиме при постоянном давлении [7-9]. Предполагается, что число Льюиса и коэффициент теплопроводности не зависят от температуры; теплообмен между газовой фазой и частицами описывается законом Ньютона; при движении на частицы со стороны газа 2013 № 2(22) Математика и механика действует сила сопротивления, описываемая законом Стокса, из-за малой объемной концентрации частиц пренебрегаем их обратным влиянием на движение газа. Система уравнений, описывающая процессы в горючем газе с инертными частицами, в безразмерной форме запишется в виде (1) (2) (3) (4) (5) (6) (7) (8) 1 + [ .SCj = Le д2Cj _ 56г 59g 1 529г pk в, ч — + й— = ---gB( -9k) + p„CjC2 exp dC _ _ -ГГ + u ^ = ^ —Г - YpgC1C2 exp дт 1 + в дС2 Le д2С2 ( v2^2 Y Yw \ YpgC1C2 exp pg ^Г 5C, - + u дт д£, 1 +( pg v1h 20 д_g дРgU n —- +—— = 0 ; дт д£, gpk =0. дт д£, д9к _д9к —- + w—дт д£, (9g-9k). 2 3 PrC pk(u - w); дт д£, 2 yNu pg =(1 -p/y)/(1+ P9 g). Граничные условия: д9g (oo,т) SCj (0,т) 9g (0, т) = 0 дC2 (0, т) = 0, дC1 (o, т) = 0 дC2 (о, т) д^ =0 =0, =0, д9k (0,т) д_k (0,т) =0, = 0, U (0, т) = 0, w (0, т) = 0. (9) д£, д£, Начальные условия: 9g (|,0) = -Vy , 9k (1,0) = -1/y, C, (|,0) = C|Q, C2 (|,0) = 1, U(|,0) = 0, w(|,0) = 0, pk (|,0) = 1. (10) Уравнение (1) описывает перенос тепла в газовой смеси за счет теплопроводности, конвекции, с учетом теплообмена между газовой фазой и частицами, тепловыделения за счет химической реакции; (2) - уравнение концентрации горючего; (3) - уравнение концентрации окислителя; (4), (5) - уравнения сохранения массы газа и частиц; (6) - уравнение энергии частицы, с учетом ее движения и теплообмена с газом; (7) - уравнение движения частиц. (8) - уравнение состояния идеального газа, записанное при постоянном давлении. При записи системы уравнений выбраны следующие масштабы: времени t* = RT2cg/ (QEpg0Yi0Y20k0exp(-E/(RTa))); координаты x* = (^t»/(cgpg0))1/2; скорости V* = х*И*; плотности газа pg0 - начальная плотность газа; приведенной плотности частиц рк0 - начальная приведенная плотность частиц; температуры T* = RTa2/E, где Ta - адиабатическая температура продуктов сгорания газа в изобарических условиях. Безразмерные переменные: т = t/t* - время; | = х/х* - координата; Qg,k = (E/RTa )/(Tg,k - Ta) - температура; C1 = Y1 / Y10, C2 = Y2 / Y20 - относительные концентрации горючего и окислителя в газовой фазе, где Y1 =р1 / pg , Y2 =р2 / pg, pg = р / р 0 - плотность газа; pk = pk / pk0 - приведенная плотность частиц; u = u / V* - скорость газа; w = w / V* - скорость частиц. Индексы: g - газ, k - частицы, 0 - начальные условия, 1 - горючее, 2 - окислитель, * - масштаб, черта над обозначением - безразмерная величина. Безразмерные параметры (критерии): Pr = n/g/X - число Прандтля, Nu = ark/X -число Нуссельта, Le = Dcgpg0/X - число Льюиса, C = ck /cg - отношение теплоемко-стей, B = 4nrk3ckpk0n/(3cgpg0) - начальная массовая концентрация частиц, X = rk2ckpk0/2NuXt» - размер частиц, у = cgRTa2/(QEYi0), в = RTJE - параметры, характеризующие тепловой эффект и энергию активации химической реакции. Обозначения: t - время, х - координата, T - температура, Y - отношение плотности окислителя/горючего к плотности газовой смеси, р - плотность, u - скорость газа, w - скорость частиц, c - удельная теплоемкость, Х - коэффициент теплопроводности, D - коэффициент диффузии, ng - коэффициент вязкости, д - молекулярная масса, v - стехиометрический коэффициент, a - коэффициент теплоотдачи, rk - радиус частиц, n - количество частиц в единице объема, R - универсальная газовая постоянная, Q - тепловой эффект реакции, k0 - предэкспоненци-альный множитель, E - энергия активации химической реакции. Математическая модель (1) - (10) описывает распространение пламени в газовой смеси с инертными частицами с учетом теплового расширения. При u = 0, w = 0, pg = 1, модель переходит в теплодиффузионную модель горения газодисперсных сред [10-13]. Система уравнений (1) - (10) решалась численно. Уравнения (1) - (3) решались по неявной разностной схеме, конвективное слагаемое аппроксимировалось разностями против потока. Значения температуры газа на новом временном слое использовались для нахождения плотности газа из уравнения (8). Затем из уравнения сохранения массы (4) определяли скорость газа. Уравнения (5) - (7) решались по явной разностной схеме. Расчеты проводились при значениях параметров: у = 0,1, в = 0,07, Le = 1, Nu = 1, vi д = 16, v2^2 = 64, Y10 = 0,05, Y20 = 0,2, Pr = 1, C = 1, B = 0,15. В расчетах варьировались значения C1 в интервале от 1 до 0,7, параметр х в интервале от 2000 до 500000, B в интервале от 0 до 0,25. В численных экспериментах горение газовзвеси инициировали передачей тепла от горячей стенки. В момент времени, когда фронт пламени достигал безразмерной координаты | = 500, граничное условие 6g (0, т) = 0 заменялось граничным условием 56g (0, тЦ/д^ = 0. Скорость фронта определялась скоростью перемещения координаты с концентрацией горючего C} = 0,5. На рис. 1 представлены распределения параметров среды в последовательные моменты времени с интервалом безразмерного времени Ат = 120000 для модели без учета теплового расширения, на рис. 2 - для модели с учетом расширения при значениях х = 50000, C1 = 0,8, B = 0,15. Температура частиц во фронте горения ниже, чем у газовой смеси, за фронтом их температуры выравниваются (рис. 1, а, 2, а). Концентрация горючего, взятого в недостатке, за фронтом пламени обращается в ноль. Окислитель, взятый в избытке, расходуется до концентрации C2 = 0,2 (рис. 1, б, 2, б). е C1,C2 е C1,C2 Рис. 2. Распределения температуры (а), концентрации (б), скорости (в), плотности (г); (а, в, г: сплошная линия - газ, пунктир - инертные частицы), (б: сплошная линия - горючее, пунктир - окислитель). Модель с учетом расширения Рис. 1. Распределения температуры (а) (сплошная линия - газ, пунктир -инертные частицы), концентраций (б) (сплошная линия - горючее, пунктир - окислитель). Модель без учета расширения Перед установившимся фронтом горения скорость газа и частиц имеет одинаковое значение, позади фронта скорость газа и частиц равна нулю, во фронте прогрева газовой фазы частицы имеют скорость выше скорости газовой фазы из-за их механической инерции (рис. 2, в). Частицы в области химических превращений имеют увеличенную приведенную плотность (рис. 2, г). В силу этого эффекта в зоне интенсивных химических реакций локально присутствует большее количество частиц, которые локально понижают температуру в зоне реакции сильнее, чем это получается в расчетах на основе модели без учета расширения газа и относительного движения частиц. Подогретые в зоне химических реакций частицы и имеющие скорость большую, чем газ, переносят тепло в зону прогрева перед фронтом горения. Это способствует растяжению волны горения в газе. На рис. 3 представлены зависимости скорости горения от начальной концентрации горючего, полученные на основе модели без учета теплового расширения (рис. 3, а) и с его учетом (рис. 3, б) при массовой концентрации частиц B = 0,15 и трех значениях параметра х, характеризующего размер частиц. Видно, что тепловое расширение увеличивает скорость распространения фронта пламени в чистом газе на величину от 20 до 30 % (кривые 1 на рис. 3, а, б). Присутствие крупных инертных частиц в газе не сказывается на скорости распространения пламени в модели без учета теплового расширения газа (кривая 2, рис. 3, б). Расчет на основе модели с учетом теплового расширения показывает, что скорость пламени уменьшается, несмотря на то, что частицы во фронте пламени не успевают нагреться, забрав тепло из зоны химической реакции. Этот эффект обусловлен меньшей температурой газа позади фронта реакции и в связи с этим меньшим объемным расширением газа. С уменьшением размера частицы начинают влиять на скорость химической реакции во фронте горения, уменьшая в нем температуру (кривые 3, 4, рис. 3, а, б). Существует интервал значений х, при которых частицы слабо влияют на скорость пламени при стехиометрическом и близком к нему соотношении окислителя и горючего и сильно влияют на скорость при значительном отклонении состава газовой смеси от стехиометрической (кривые 3, рис. 3, а, б). Рис. 3. Зависимости значений скоростей фронта пламени от концентрации горючего: а - модель с учетом расширения (1 - незапыленная газовая смесь, 2 - х = 500000, 3 - х = 50000, 4 - х = 2000); б - модель без учета расширения (1 - незапыленная газовая смесь, 2 -х = 300000, 3 - х = 50000, 4 - х = 2000) На рис. 4 представлены зависимости относительной безразмерной скорости распространения фронта пламени от параметров дисперсной фазы - B и х при двух значениях начальной концентрации горючего в газовой смеси (С1 = 0,85 и C1 = 1), рассчитанные по модели, учитывающей тепловое расширение. V0 - значение скорости фронта горения для чистой газовой смеси, которое имеет значение при C1 = 1 - V0 = 0,02, при C1 = 0,85 - V0 = 0,01. Параметр х изменялся в следующем интервале: 2000 < х < 500 000, параметр B : 0 < B < 0,25. Из рис. 4 видно, что при фиксированном х с увеличением массовой концентрации скорость пламени уменьшается. Чем мельче частицы, тем их влияние на скорость пламени значительнее. Рис. 4. Зависимости значений относительных скоростей фронта пламени от параметров дисперсной инертной фазы, пунктир - концентрация C1 = 0,85, сплошная линия - C1 = 1: а - зависимость от величины, характеризующей размер частиц х (1 - B = 0,05, 2 - B = 0,1, 3 -B = 0,15, 4 - B = 0,2, 5 - B = 0,25,); б - зависимость от массовой концентрации B (1 - х = 500 000, 2 - х = 120 000, 3 - х = 50 000, 4 х=2000) Сформулирована математическая модель горения газовой смеси с взвешенными в ней инертными частицами, учитывающая тепловое расширение газа и относительное движение фаз. Получены зависимости скорости распространения пламени от массовой концентрации, размера инертных частиц и концентрации горючего в газовой смеси.

Скачать электронную версию публикации