Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги из ГПНТБ / Горение пористых конденсированных систем и порохов (препринт)

..pdf
Скачиваний:
8
Добавлен:
19.10.2023
Размер:
773.2 Кб
Скачать

АКАДЕМИЯ НАУК СССР

ОРДЕНА ЛЕНИНА ИНСТИТУТ ХИМИЧЕСКОЙ ФИЗИКИ ОТДЕЛЕНИЕ ОРДЕНА ЛЕНИНА ИНСТИТУТА ХИМИЧЕСКОЙ ФИЗИКИ

В. Ф. Дубовицкий, В. Г. Коростелев, А. И. Коротков, Ю. В. Фролов, А. Н. Фирсов, К. Г. Шкадинский,

С. В. Хомик

ГОРЕНИЕ ПОРИСТЫХ КОНДЕНСИРОВАННЫХ

СИСТЕМ И ПОРОХОВ

(Препринт)

ЧЕРНОГОЛОВКА 1974

Послана в журнал «Физика горения и взрыва»

S

Гос. п блинная

I

1

научно-техническая

!

!

библио : ска < OCP

і

I

аизЕмпляр

j

ЧИТАЛЬНОГО ЗАЛА

j

 

 

УДК 536.463

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

пературных продуктов горения или предварительно разогретого в зоне реакции газофазного реагента, т. е. путем конвекции.

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

в различных областях теплофизики и теплотехники,

в химии и

т. д. Принципы

конвекции используются

и в ствольной артил­

лерии: быстрое

воспламенение и горение

пороховых зарядов

обеспечивается

проникновением в поры заряда высокотемпера­

турных продуктов сгорания воспламенителя [1].

В последнее

время начинают активно развиваться основанные на конвектив­ ном горении методы синтеза новых конструктивных и жаропроч­ ных материалов.

Но несмотря на столь широкое использование и большой ин­ терес, проявляемый к теории конвективного горения, его режи­ мы и закономерности изучены весьма слабо, в основном качест­ венно. Исследования в этой интересной области горения были начаты А. Ф. Беляеым и К. К. Андреевым и в дальнейшем про­ должены их учениками и последователями [2—11]. Основное внимание в этих работах было уделено условиям перехода нор­ мального послойного горения пористых систем к ускоренному режиму конвективного горения и детонации. Логическим про­ должением этих работ явились исследования по выявлению ус­

ловий горения в щелевых зарядах [5]. Но многие стороны этого явления оставались и остаются далеко не изученными.

Задачей данного исследования явилось дальнейшее изучение процесса развития конвективного горения пористых гетероген­ ных конденсированных систем и порохов на образцах ограничен­

1

ной длины с различными граничными условиями и на базе это­ го— разработка математической модели их горения [11].

Задача о конвективном горении пористых систем, как пока­ зывает анализ, может быть разбита на две части: а) собственно горение пористых систем при конвективном механизме теплопе­ редачи, б) диспергирование из зоны горения дисперсных частиц и их дореагирование в объеме. Между двумя этими процесса­ ми, как правило, существует тесная связь.

Если условия применения пористых систем таковы, что до­ горание продуктов диспергирования существенно влияет на тем­

пературу или давление в зоне

горения,

что,

например,

имеет

место при сжигании в жесткой оболочке

(малый замкнутый или

полузамкнутый объем), то процесс, как это

уже неоднократно

отмечалось

[3—5], является

неустойчивым

и развивается до

скоростей в

сотни—тысячи метров в секунду

вплоть до

дето­

нации. При этом развитие конвективного горения зависит от ре­ жима диспергирования и догорания указанных частиц. Если же горение происходит в большом объеме и пористые образцы за­

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

характер горения может быть совершенно иным. В частности,

здесь может быть поставлен вопрос о стационарном режиме конвективного горения.

Поскольку практически все предыдущие исследования прове­ дены в жестких оболочках, в данной работе основное внимание

было уделено режиму горения пористых систем при постоянных внешних (за зоной горения) условиях. Рассмотрен режим «вы­ нужденного» конвективного горения, когда его развитие начи­

нается с фильтрации высокотемпературных газов под действием созданного инертным газом (азотом) перепада давлений в по­ рах образца и внешнем объеме — объеме установки постоянного

давления БПД-150.

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

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

ции скорости распростанения

фронта горения по длине образ­

ца, давления в зоне горения,

оценки температурного градиента

в зоне горения в зависимости

от структуры, размеров, состава

2

и условий воспламенения образца, т. е. от характерных пара­ метров образца.

Для исключения влияния эффектов догорания диспергиру­ емой части вещества эксперименты проводили в установке пос­ тоянного давления с образцами ограниченной длины, заброниро­ ванными для предотвращения заполнения пор инертным газом

весьма тонкой разрушающейся (слой эпоксидной смолы) или сгорающей (слой баллиститного пороха) оболочкой.

Регистрация скорости распространения процесса осуществ­ ляется с помощью серпи термопар или тонких проволочек, за­ делываемых па равных расстоян иях в ооразец перпендикулярно

 

его оси.

В первом

случае

ре­

 

 

гистрировалось

и

изменение

 

 

температуры

в зоне

прогрева

 

 

образца.

 

полноты

реагирова­

 

 

Оценку

 

 

ния в

зоне горения

проводили

 

 

па основании сравнения осцил­

 

 

лограмм

записи

давления

у

 

 

донного торца

п

в объеме

 

ус­

 

 

тановки.

 

 

 

 

 

 

 

 

 

IIa

рис. 1

приведены

 

ре­

 

 

зультаты

 

измерения скорости

 

 

распространения

фронта

кон­

 

 

вективного

горения

но длине

 

J.

пористого образца из смесевых

X, MM

составов

гексоген—Г1ХА (ГГ—

 

■'

ПХА)

[11].

Для рассматрива­

Рис. 1. Скорость распространения

 

емых условий

можно выделить

фронта конвективного горения по дли­

 

три характерные зоны.Ila

уча­

не пористого образца системы IlXA-

 

стке I

определяющим является

ГГ с различной температурой горе­

 

ния:

 

установление фронта фильтра­

7_ 7'r = 3200oK; 2 — Tr=SlOO0K

 

ции, разогрева

и воспламене­

приблизительным ПОСТОЯ HCT-

 

ния. Участок II характеризуется

 

вом конвективного горения. IIa участке IIl скорость горения па­

 

дает по мере приближения фронта к закрытому концу.

 

Картина качественно не меняется при изменении дисперсно­

 

сти компонентов.

Однако

с ростом размера частиц при одной и

 

той же плотности

образца скорость горения возрастает.

>Весьма существенное влияние на среднюю скорость конвек­

тивного горения V оказывает энергетика состава или температу­

ра продутов его горения. Так, снижение Tr

с 3200 до 2IOOoK за

счет изменения соотношений между гексогеном и HXA привело

к падению скорости в области H с 200 до 60 c.M cen.

При этом,

если для состава с высокой Tr на первом

участке

мы имели

3

спад

скорости υ, то для состава

с низкой Tr при одной

и той

же температуре газов наблюдается постепенный рост ѵ, анало­

гичный наблюдаемому

в случае

горения

единичной поры

[5].

Таким образом, на первой стадии развития конвективного горе­

ния

наблюдается

два

характерных

нестационарных

режима,

связанных с установлением фронта фильтрации

и тепловыделе­

ния.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

В таблице приведены результаты измерения

максимального

давления, регистрируемого в

различных сечениях горящего об­

разца,

забронированного тонким слоем эпоксидной смолы.

Здесь •

 

 

 

 

 

 

 

 

же приведены значения

 

 

 

 

 

 

 

 

конечного

давления

в

Размер

час­

 

Давление

 

 

объеме установки.

 

Из

Длина

Давление

таблицы

следует,

что

тиц перхлора­

образца.

в порах

в бомбе,

давление

в

зоне

горе­

та аммония,

MM

образца,

 

атм

 

MKM

 

 

атм

 

 

ния в 1,5—2 раза пре­

 

 

 

19

 

144

 

78

вышает

внешнее

 

дав­

 

 

 

 

 

ление и

несколько уве­

 

 

 

38

 

168

 

85

10012.6

 

 

личивается

 

по

 

мере

*38

 

258

 

 

 

 

 

 

 

 

увеличения

длины

об­

 

 

 

57

 

195

 

97

разца. Давление в объ­

 

 

 

19

 

95

 

74

еме также увеличивает­

200—250

 

 

ся

па

15—20%.

C

38

 

130

 

81

 

 

уменьшением

дисперс­

 

 

 

 

 

57

180

89

* Эксперименты с образцами в прочной метал­ лической оболочке. Относительная плотность образцов A = 0,8.

ности кристаллов ПХА

максимальное давление уменьшается.

Скорость нараста­

ния давления в измеря­ емом сечении составляет (5÷7)∙106 атм/сек. Это на 2—3 поряд­

ка выше, чем изменение давления в условиях чистой фильтрации при аналогичном перепаде давлений Ap на торце пористой сре­ ды.

Сочетание повышенного давления в зоне горения с большим градиентом dρ dt способствует увеличению скорости фильтра­ ции, а следовательно, и скорости распространения фронта вос­ пламенения и горения. Сравнение данных графика на рис. 1 и таблицы показывает, что изменение скорости распространения υ ,

V происходит пропорционально изменению давления в порах ве- ¡

щества.

Оценка полноты горения по измерениям давления (одновре­

менно в порах вещества и в объеме) показала, что для рассмат­ риваемых систем во фронте горения сгорает не более 20—30% вещества. Остальное же вещество в случае сгорающих и легко

4

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

Рассмотрим далее вопрос о развитии конвективного горения.

Экспериментальные результаты, представленные на рис. 1, от­ мечают постоянство либо слабый рост скорости горения, сопро­

вождаемый постепенным увеличением давления во внешнем объеме. Однако они не дают определенного ответа на вопрос, \/ может ли быть конвективный режим горения стационарным. Ответ па этот вопрос может быть получен либо при детальном

изучении фронта конвективного горения с одновременным заме­ ром давления, температуры и т. п., либо при проведении экспери­ ментов на образцах, размеры которых значительно больше ве­ личины зон фильтрации и горения. Действительно, если справед­ лива идея о спонтанном развитии конвективного горения за счет наложения тепловых воли прогрева [5—8], то увеличение длины

образца и при отсутствии оболочки

должно

приводить к росту

скорости горения.

 

 

 

Эксперименты па образцах длиной 200—150 мм (в установке

постоянного давления объемом более 20 л)

показали (рис. 2),

что режим горения является

пульсирующим

со средней скоро­

стью, близкой к полученной

ранее

на укороченных образцах.

Рис. 2. Скорость распространения фронта конвективного горе­

ния на удлиненных образцах пороха при постоянном внешнем давлении р=70 атм

5

Однако пульсации в этих опытах связаны не с особенностями процесса конвективного горения, а с периодическим разрушени­

ем хотя и тонкой, по все-таки существующей оболочки. Эти экс­ перименты показали, что стационарный (или квазистацнопарпый) режим конвективного горения пористых систем возможен. Но поскольку стационарный режим конвективного горения воз­

можен, ему должна соответствовать и определенная

структура

фронта горения.

 

 

 

что

Из самого определения конвективного горения следует,

фронт горения должен включать в себя,

по крайней мере,

три

характерные зоны: зону фильтрации, зону прогрева

вещества и

зону химического превращения (сгорания) вещества.

 

 

О величине зоны фильтрации дает представление область III

(рис. I). Существование зоны спада скорости

горения связано

с тем, что фронт фильтрации, опережая фронт

воспламенения,

приводит к заполнению пор продуктами

горения, снижению пе­

репада давления ∆p, а в силу этого и скорости горения [1 1].

По­

этому расстояние, начиная с которого чувствуется эффект сни­

жения ѵ, дает возможность судить

о глубине зоны фильтрации.

Для исследуемых систем оно составляет приблизительно

15—

20 мм, причем с увеличением Tr

намечается тенденция

к его

уменьшению.

Оценка зоны прогрева проведена по осциллограммам записи температуры T (/) в различных сечениях образца и по фотогра­ фиям образцов, погашенных резким сбросом давления. Ее вели­ чина составляет 4± 1 мм.

Зона химической реакции (т. е. расстояние от сечения с тем­ пературой воспламенения до верхнего торца образца) находится

в пределах 1—2 мм, что на порядок меньше, чем зона фильтра­ ции. Ее протяженность, так же как и протяженность двух пер­ вых зон, зависит от физико-химических свойств среды и усло­ вий сжигания (ро, Ta и т. и.).

Дополнительно остановимся на участке III графика на рис. 1.

Как уже отмечалось, если торец образца глухой, скорость сго­

рания при подходе к нему падает.

Картина изменяется, если нижний торец образца открыт. В

этом случае скорость ѵ по мере подхода к торцу возрастает.

Этот результат аналогичен полученному в работе 16], где отме­ чалось повышение скорости горения пористого заряда у донного торца, если последний контактирует с дополнительным и сво­ бодным объемом (аналогичный эффект, по-видимому, должен быть и при распространении горения в единичной поре). Уско­ рение горения в этом случае объясняется увеличением градиен­ та dp dx при подходе к открытому торцу, сообщающемуся . с

внешней атмосферой, т. е. увеличением скорости фильтрации. Но

6

с последней, как известно, тесно связаны условия теплопередачи от газа к конденсированному веществу, т. е. прогрев и воспламе­ нение вещества.

Рассмотренные физические аспекты

горения конденсирован-

,

ных пористых систем

(см.

также [12]) позволяют сформулиро-

'

вать математическую модель горения.

 

однородная

конденси­

 

Имеется

пористая,

газопроницаемая,

 

рованная среда, представляющая собой совокупность сферичес­

 

ких частиц. Длина пористой среды в направлении

оси х ограни­

 

чена

и равна L.

Параметры

пористой среды:

относительная

 

ПЛОТНОСТЬ

/X = p/p,,,aʃ,

HOpHCTOCTb Ш,

рЭДИѴС ЧЗСТИЦ R0.

Поры

 

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

 

нием р„; начальная температура газа и пористой среды T0.

 

 

Теплофизические параметры пористой среды и газа считают­

 

ся известными. Сечение X = О конденсированной пористой среды

 

граничит с объемом, давление р0 и температура Тв„,

в котором

 

считаются заданными и постоянными.

 

 

 

 

газонепрони­

 

Правая граница

(х = L)

пористой среды либо

 

цаема, либо открыта. В последнем случае давление и температу­

 

ра па пей аналогичны начальным давлению и температуре в об­

 

разце

(/?

= /?„; 7’

= T0).

t

= О газ,

нагретый

до температуры

 

В начальный

момент

 

Tatt > T0,

из внешнего объема под действием перепада давлений

 

в объеме и порах ∆p = р0 — pll начинает проникать в поры в со­

 

ответствии с законом фильтрации.

 

 

 

 

между

газом п

 

В

результате

конвективного

теплообмена

 

частицами последние

прогреваются газовым

потоком по закону

 

прогрева индивидуальных

сферических частиц.

Считается,

что

 

теплообмен между частицами теплопроводностью отсутствует.

 

По мере прогрева частиц

па их поверхности начинается хи­

 

мическая

реакция,

сопровождающаяся

газовыделением;

ско­

 

рость газообразования

равна w.

Образующиеся

газообразные

 

продукты разложения вступают в реакцию горения, скорость ко­

 

торой считается бесконечно

большой

по сравнению

с

w.

При

 

этом температура продуктов реакции равна адиабатической тем­

 

пературе горения 7’ад.

 

 

 

 

 

 

 

 

 

к интенсив­

 

Газовыделение и подъем температуры приводят

 

ному росту давления

в зоне реакции.

Поскольку

давление во

 

внешнем объеме считается постоянным, отток газа из зоны реак­

 

ции будет проходить

 

как в

поры вещества,

так и

во внешний

 

объем, вызывая диспергирование к-фазы.

Сопротивление тече­

 

нию газа в образец создается пористой средой.

Сопротивление

 

оттоку газа из зоны горения

во вне зависит

от течения газа в

 

зоне реакции

и

от наличия перед

ней пористого

каркаса.

По­

 

следнее определяется полнотой химических реакций разложения конденсированной фазы пороха.

7

При расчете используется предположение, что величина кар­ каса (считается заданной) по мере сгорания образца остается все время постоянной. Его газопроницаемость определяется полнотой выгорания η. За каркасом диспергируемые во внешний объем частицы догорают, не повышая в нем р0 и Tnu.

Пористая среда является изотропной и состоит из сферичес­ ких частиц радиуса Rfj.

Внутренние силы взаимодействия между отдельными части­ цами пористой среды и эрозионные эффекты на повехпости час­ тиц не учитываются. Суммарная полнота выгорания ηκ частиц

считается конечной (задается в пределах 0,05—0,20).

Фронты фильтрации, прогрева и горения принимаются плос­

кими.

Развитие горения пористых систем, как это следует из выше­

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

матичесая модель включает уравнения,

 

описывающие все эти

процессы.

 

 

 

1.

Уравнение сохранения массы

 

 

 

дртKJUIfl

. I/д (Ulllllрти)I

 

 

 

—------ 1----- -- —- = W.

 

 

∂t

дх

 

 

2.

Уравнение сохранения импульса

 

 

du

др

— uw.

 

рт • —

- — т • —— К

 

dt

дх

 

 

3.Уравнение сохранения энергии

4.Уравнение состояния газа

Эти уравнения составлены для газовой фазы в пористом ве­ ществе.

5.

Уравнение теплопроводности для конденсированных час­

тиц

∂Tκ

 

{∂2T

 

 

 

—- = аJ---

 

∂t

κ

дг-

6.

Уравнение скорости

химической реакции (газовыделения)

W= Spκko exp (-E∕RTκ s).

7.Уравнение глубины превращения к-фазы

de _ W dt pκ∆

8

Система уравнений решается при следующих начальных и граличных условиях:

t = O: T = Tκ= T0,

и = О,

η = О, р = р„\

X = Xk

P = Po,

dn = О;

 

дх

λ'k(0) = 0]

'Γ = Tβ",

X = L: u(t, L) = О — торец глухой;

дх- = O —торец

открыт;

г = 0:

∂Tκ

 

дх

r = rs-. lκ.^ = ^τ-τκs).

дг

Здесь: t — время; х—пространственная координата, 0≤.v<L;

/¿0, rs = R0(l—η)13 — начальный и

текущий радиус частиц;

T —

температура

газа;

Tκs — температура па

поверхности

частиц;

T0— начальная температура

газа

и

пористой

среды;

ро —

начальная

плотность

пористой

среды;

 

рк—плотность

частиц

(к-фазы);

р — давление;

и — скорость газа;

/г —удельная эн­

тальпия;

s — удельная

поверхность;

ск—теплоемкость;

λ —

теплопроводность;

μ — динамическая

 

 

вязкость

газа;

 

Q —

тепловой эффект химической реакции;

k0 — предэкспонент;

E —

энергия активации;

η — глубина превращения;

т — коэффици­

ент пористости;

— относительная плотность

ʌ = p pmaχ',

ɑ—

коэффициент

теплоотдачи

 

Nuk

ак—температуроировод-

α = 9-;

 

 

 

 

 

 

 

Z/

s∙

 

 

 

 

 

 

 

 

пость, ακ- — ; g— ускорение земного

 

притяжения;

А— меха-

 

CpK

 

 

 

L — длина образца;

хк

текущая

ническпй эквивалент работы;

координата каркаса;

R — газовая постоянная.

 

 

 

 

 

Закон сопротивления

при движении

 

газа

в пористой среде

взят в форме [13]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

pMS (4,7/г

4-

1,8

і n 7Γ

 

+ 0,321

и |

),

 

 

 

 

 

 

т

 

 

 

2

 

 

 

 

 

 

 

 

 

где

 

 

SOl

3ʌ/* Л

m =—

»

4 г і

 

\

 

 

 

п = _S-——,

 

I — ʌ (

1 — и),

 

 

 

 

 

m?

 

R3a

 

 

 

 

 

 

 

 

 

 

 

--= ɪ ( ɪ - η) eχp (-e∕* t× s)= ⅞

 

 

 

 

 

∂t

R11

 

 

 

 

 

 

 

Pkδ

 

 

 

 

9

Соседние файлы в папке книги из ГПНТБ