Моделирование режимов работы низкотемпературного газогенератора

Расчет и конструирование

УДК 519.6:536.24

МОДЕЛИРОВАНИЕ РЕЖИМОВ РАБОТЫ НИЗКОТЕМПЕРАТУРНОГО ГАЗОГЕНЕРАТОРА

В.В. Кириллов

Южно-Уральский государственный университет, г. Челябинск

Низкотемпературный газогенератор (НТГГ) предназначен для получения газа с температурой 350...450 К. Он состоит из двух основных частей - камеры сгорания с зарядом твёрдого топлива и камеры охлаждения с гранулами твёрдого охладителя. Продукты сгорания твёрдого топлива имеют высокую температуру, что не позволяет использовать их в целом ряде технических устройств. Охлаждение газа происходит в результате его взаимодействия с гранулами твёрдого охладителя, на термическое разложение которого расходуется часть внутренней энергии газа. Для описания процессов тепло- и массообмена в НТГГ разработана математическая модель на основе одномерных уравнений неразрывности, импульса, энергии, баланса массы компонентов газовой смеси. Газовая смесь состоит из продуктов сгорания твёрдого топлива, продуктов разложения охладителя и воздуха. Математическая модель учитывает воспламенение и горение топлива, течение газовой смеси в НТГГ, теплообмен с элементами конструкции НТГГ и кинетику разложения гранул охладителя. Для реализации математической модели разработан численный метод решения краевой задачи на основе метода конечных разностей. Система уравнений аппроксимирована неявными разностными схемами. Для решения нелинейных разностных уравнений применён метод Ньютона. Краевая задача решалась ортогональной прогонкой. Результаты расчётов сравнены с экспериментальными данными. Выяснено влияние конструктивных и режимных параметров на изменение температуры, давления, расхода, скорости разложения охладителя.

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

НТГГ применяются в системах пожаротушения [1-4], для раскрутки турбин турбонасосных агрегатов [5-7], производства низкотемпературного газа для наддува ёмкостей различного назначения [8-12]. Для описания процессов тепло- и массообмена в НТГГ применялись математические модели различного уровня сложности [11, 12]. Математическое моделирование НТГГ позволяет сократить сроки и стоимость проектирования, а также выявить влияние различных факторов на ход рабочих процессов, что достаточно сложно и дорого выполнить экспериментальными исследованиями.

На рис. 1 представлена схема экспериментального низкотемпературного газогенератора [11].

Экспериментальные исследования [11, 12] показали наличие заметных отклонений в характере изменения температуры (рис. 2) и давления (рис. 3) на выходе из НТГГ. Данное обстоятельство может быть вызвано как погрешностями проведения экспериментов, так и отклонениями размеров конструктивных элементов НТГГ.

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

14 15 16 17 18 19 20

Рис. 1. Схема низкотемпературного газогенератора: 1 - воспламенитель; 2 - передняя крышка; 3 - теплоизоляция; 4 - передний объём; 5 - передняя решётка; 6 - стенка камеры сгорания; 7 - теплоизоляция; 8 - заряд топлива; 9 - задняя решётка; 10 - задний объём; 11 - теплоизоляция; 12 - переходник; 13 - патрубок; 14 - коллектор; 15 - газоход; 16 - верхняя стенка кассет; 17 - боковая стенка кассеты; 18 - кольцевой канал; 19 -наружная стенка камеры охлаждения; 20 - кассета с охладителем; 21 - задняя стенка;

I - камера сгорания; II - камера охлаждения

400 ев 350

Ц 300 250

0 0,4 0,8 1,2 1,6 2,0 2,4 2,8 3,2 3,6 4,0 время, с

Рис. 2. Изменение температуры на выходе из НТГГ

0'

0 0,4 0,8 1,2 1,6 2,0 2,4 2,8 3,2 3,6 4,0

? ? 1 5 ? ?

время, с 1—2—3

Рис. 3. Изменение давления на выходе из НТГГ

1. Математическая модель НТГГ

Для анализа рабочих процессов в НТГГ применялась одномерная математическая модель [13]. Математическая модель камеры сгорания (КС) включает в себя уравнения рабочих процессов в воспламенителе, переднем объёме, канальной части с шашками топлива и в выходном объёме. Рабочие процессы в воспламенителе, переднем и выходном объёмах камеры сгорания

описываются в рамках термодинамических математических моделей, включающих в себя уравнения баланса массы, энергии смеси газов, уравнения баланса массы компонентов газовой смеси, уравнения теплопроводности стенок. Газовую смесь в КС можно разделить на неконденсирующиеся (НК) продукты сгорания топлива (ПСТ) и воспламенителя (ПСВ), такие, как углекислый газ, окись углерода и т.п., воздух и водяной пар. Предполагается, что все компоненты газовой смеси с достаточной точностью подчиняются уравнению состояния совершенного газа.

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

дФ + дТ _

дх дх

В камере сгорания

Ркс

Ф _

кс кс

Окс

п

Рксекс ^кс

О

^ ' + ркс ^кс ркс Лкс

°ксекс - Ркс^кс^к

Е _

ЫтРт Пт

О1П

8^к2срк

РтЫтКПт - О,

V2

_с Т + —^•

укс кс 2 ' "кс"кс£>кс,п + дОкс<?кс,п

+-

_дркс^кс?к

дх

дркс ^кс gкс,в

дХ

др кс ^кс gкс,т

дХ

Ркс РксRксTкс •

+

дx

д°ксЁк дx

дОкс g

Окс акс (Tкс

_ ЫтРтПтРп,т ;

_ 0;

кс кс,т

дx

_ ы,

Рт Пт (1 -рп,т );

дТ„,

дх

дАЯ^ г

г дг

дГ ^ ^

г

дг

д % г + д[((-2^Лr-aITU

дх АЯ ду2

ду

+ ЧуГАЯ;

у _•

-АМ; АЯ _Я2(х)-Я1(х). АЯ

В (1)-(9) р - плотность; О - расход; р - давление; Т - температура; S - площадь свободного сечения; V - скорость; К - энтальпия; ы - скорость горения; П - периметр; О - тепловой поток; g - массовая доля; Рп,т - массовая доля водяного пара в продуктах сгорания топлива; Я - газовая

постоянная; а - коэффициент теплоотдачи; су - удельная теплоёмкость при постоянном объёме; а - коэффициент температуропроводности; Ях, Я2 - внутренний и наружный радиусы шашки топлива; ^ - коэффициент гидравлического сопротивления; х - время; х - продольная координата; г - радиальная координата; индексы: кс - камера сгорания; т - топливо; п - водяной пар; в - воспламенитель; w - стенка; м - металл.

Процессы тепло- и массообмена в кассетах с охладителем описываются одномерными уравнениями

дФ дТ

+-_ Е .

дх дг

кс

е

а

м

а

т

г

Векторы Ф, F имеют следующий вид:

PiSi

Ф _ Gi ;

PieiSi'

Gl , gV, + рД'

Giei — P,SV,

; f _

Пэкв Jm Пэкв Jc ;

Si Pi s Gi Пэфф , „ ' „ '

--S -+ S V'Jm - Si VJ

Г 8S Pi

Q

/

— П j h + П i r

экв./ m s экв./ c s

SPiSigi,n , SGigi,n _ S '  P S '

" +---_ Si JmPx — Si Jc ;

Si

ФДх,

Si

ФД^Х Si

dTwi a,.

• + -

+ -

Sr

sg^

Sr

sg^

Sr

_ 0;

_ Si Jm (1 — P х );

1 (T-1 — Tw,i )Fw,i , ai (T- — Tw,i )Fw,i .

• +

cmpm^w

^х смрУ,

мг м Ш

£ = 2%г1г.

В уравнениях (10)-(15) / - номер кассеты; в - пористость слоя гранул; Рх - массовая доля водяного пара в продуктах разложения охладителя; ]т - плотность потока массы продуктов разложения охладителя; ]с - плотность потока массы конденсации водяного пара; г - теплота парообразования; I - ширина кассеты; V - скорость; Пэфф - эффективный гидравлический периметр кассеты; V - объём; х - охладитель.

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

В газоходе и кольцевом канале векторы Ф, Т, ^ имеют следующий вид:

0

Ф _

PS

G

peS

G

Gw + pS Ge + pwS

F _

— S

pw2 П

— П ^

SPSgj SGgj _ .

+--:— _ 0 , J = в, т, п, х.

дх дх

Температурное поле в стенке описывается уравнение теплопроводности вида (8) с граничными условиями третьего рода.

Начальные условия системы уравнений определяют значения температуры, расхода и давления в начальный момент времени. Температура всех элементов одинакова и равна некоторой начальной температуре Т0 . Давление во всех точках НТГГ равно начальному р0, расход равен нулю. Массовые доли продуктов сгорания воспламенителя, топлива и продуктов разложения охладителя равны нулю.

Рассмотрим постановку граничных условий. На входе в канальную часть КС

с (х,0) .

Т^ксЛ. (l)_- ^

Ткс1 — 1

Ркс (i,0) + Pre wL

Ткс — 1 Pкс1 (i) 2

Ркс1(i) _ Ркс(i,0)+P™ (1,0)2(i,0) ;

§кс,п (i,0) _ gкс1,п (i) ; gкс,в (i,0) _ gкс1,в (i) .

В (18)-(20) к = ср!еу; индекс кс1 обозначает параметры переднего объёма КС. На входе в газоход

Ркс2 (0=Рг (,0)+Арю; ккс2 (0=К М; °кс2 = Ог М;

&кс2л

В (21) индекс кс2 обозначает параметры выходного объёма КС.

На входе в кассеты из газохода поток газа разветвляется. Часть газа входит в кассету, часть газа проходит дальше в газоход. Схема потоков представлена на рис. 4. В соответствии со схемой

аг+1 = -О,; Кг,г+1 = К,; Рг,г+1 = Рг,г -4р;

&гн,,'+1 &гн,,'; &гп,,'+1 &гп,,' .

Здесь О1 - расход газа в ,-ю кассету; Ог,, О - расходы газа в газоходе до и после точки разветвления; Ку+1, Ку - энтальпия газа до и после точки разветвления; £гп,,'+1, &гп,, , gгн,,+1, &гн,, - массовые доли водяного пара и НК газов в газоходе до и после точки раз- Рис. 4. Расчётная схема потоков в газоходе ветвления; рг,, рг,+1, Ар, - давление газа до и после

точки разветвления, а также потери давления при переходе через точку разветвления. Индекс г означает параметры газохода. На входе и выходе ,'-й кассеты

к О г,,'-1 О -1 О г, г \О ' 1 О г,,'+1 ^ О ,+1 О г,'+2

¡г ,' -1 - ' Г* , ,'+1

О, =

2кг

кг -1

Рг,г'Рг,

2 кг+1 1

[ Рг (0)' 7 - Г Р, (0)' Т

V Рг,г V V Рг,г V

Р,

(

>

2

Л'

V кг + 1 V

кг-1

МЛ

V кг + 1 V

кг-1 2кг

1кг +1

Рг,г'Рг,г ,

Рг

<

V кг + 1 V

кг-1

К (0)=;

&,,п ( ^вн ) = &гп,,; &',т ( ^вн ) = &гн,,; &',х ^ ^вн ) = 0.

О вых _

IIвых Л

м вых Л

2к,

к , -1

Р,' Рг

с Р Х-

Рк к'

кг +1

Рк к'

Л. \к, -1

к +1

Г Л

кк, +1

1

к,-1

2к Рк

—+7Р,Р,, — < 'к, +1 Р ,

Л.

А к,-1

V кг + 1

В (26), (27) ^вн - внутренний радиус кассеты; Р,, Рк - давление на выходе из г-ой кассеты и давление в кольцевом канале в точке слияния, соответственно; к,, р, - отношение теплоёмкостей и плотность на выходе из ,-ой кассеты. В кольцевом канале потоки газа из кассет смешиваются с газом в канале. Схема потоков газа изображена на рис. 5. В соответствии со схемой

Ок, ,+1 = Ок, , + 0; Кк, ,+1Ок,,+1 = Кк, А, , + ; Рк, ,+1 = Рк, , -АР,;

&ф'+1Ок,,+1 = £кы°к,, + &|,Д- , ] = т, п, х.

Здесь К, О, &т,,, &х,,', &п, ,' - энтальпия, расход, массовые доли НК газов ПСТ и ПСВ, НК газов продуктов разложения охладителя, водяного пара на выходе из г-й кассеты. Индекс к означает параметры кольцевого канала.

Расход из кольцевого канала к потребителю определяется по формулам

0к,,-2 ,'+1 0к,,-1

О

г+1

¡г

г

Ок, -1 О

к, +1

М-

О

О -1

Рис. 5. Схема потоков в кольцевом канале

к

г

г,г

к

г

2

2

г,

2

2

к

>

2

G —

и S

"вых вых.

2k

k -1

Рк Р

2 k+1 1

Г Po ^ k Г Po ^ PSL

l Рк J l Рк J Рк

с

ивых Sвых

2 1 k-i 2k

-Р Р — <

чк +1) \к + 1 к к рк ^к +1

Здесь рк, рк - плотность и давление на выходе из кольцевого канала, р0 - давление у потребителя. 2. Численный метод решения уравнений математической модели

Система уравнений состоит из обыкновенных дифференциальных уравнений и уравнений в частных производных. Для численного интегрирования обыкновенных дифференциальных уравнений применялся метод Рунге-Кутта 4-го порядка. Для решения уравнений в частных производных применялся метод конечных разностей. Краевая задача решалась ортогональной прогонкой [16].

Уравнения неразрывности, количества движения и энергии аппроксимируются неявными разностными схемами и приводятся к векторным уравнениям вида [15]

>

k

2 1 k-i

k

2 1 k-i

k +1

А ут+1 _ в у™-1 = с

Здесь Ап, Вп—1 матрицы 3 х 3 , У = О, р, Т\', Сп - вектор размерности 3. Граничные условия на входе в КС преобразуются к виду

А1У1™+1 = С,,

где А1 - матрица размерности 2 х 3 ; С1 - вектор размерности 2.

В КС и газоходе дифференциальные уравнения аппроксимировались на разностной сетке, конструкция которой представлена на рис. 6. Узел М - последний узел КС, узел М +1 - первый узел газохода. Газоход разделён на отрезки, количество которых равно числу кассет плюс один, каждый из которых содержит по несколько узлов.

1 КС

N N + 1

N

M M + 1

Рис. 6. Разностная сетка в КС и газоходе

Первые три условия перехода от КС к газоходу (22) в разностной аппроксимации имеют следующий вид:

Ym+1 —nym+i + „ 1M ~ U-1M+1 ^ A 5

m+1

где а - квадратная матрица 3 х 3 , а х - вектор размерности 3. Условия перехода в газоходе (23) в разностной аппроксимации также могут быть представлены в виде

У^+1 = а^+1 +Хг.

В конце газохода граничное условие записывается в виде

Оы = 0.

Краевая задача для КС и газохода решается ортогональной прогонкой.

Уравнения (10)-(12) аппроксимируются аналогично уравнениям КС и газохода и приводятся к виду (31). Граничные условия (24), (25) преобразуются к виду (32), условие (27) преобразуется к виду [15]

ву+1 = Сы

где Вн - матрица 1 х 3 , Сн - скаляр. Решение уравнений кассет выполняется последовательными приближениями, аналогично решению задачи для КС и газохода, после чего решаются уравнения баланса массы компонентов газовой смеси. Затем рассчитываются расходы на выходе из кассет по формулам (27).

После расчёта кассет с охладителем решаются уравнения математической модели кольцевого канала. Конструкция разностной сетки для кольцевого канала аналогична конструкции сетки в газоходе. В узлах стыковки кольцевого канала с кассетами граничные условия (28) в разностной аппроксимации преобразуются к виду (32), а условия (29) записываются следующим образом.

т+1 _ + § у

§4 N+1 _

+1

, ] = т,п,х.

Граничное условие на выходе из кольцевого канала (30) приводится к виду (36). Краевая задача решается ортогональной прогонкой.

На рис. 7-14 представлены результаты расчётов режимов 1-3.

250

0 0,4 0,8 1,2 1,6 2,0 2,4 2,8 3,2 3,6 4,0

время, с — эксперимент — расчёт Рис. 7. Изменение температуры на выходе из НТГГ

время, с — эксперимент — расчёт Рис. 8. Изменение давления на выходе из НТГГ

время, с — эксперимент — расчёт

Рис. 9. Изменение температуры на выходе из НТГГ

В расчётах выявлен рост температуры в начальный период работы длительностью 0,2 с. В экспериментах использовалась термопара с постоянной времени 0,8 с, что не позволило зафиксировать подъём температуры в начальный момент, вызванный горением воспламенителя.

время, с

-эксперимент — расчёт

Рис. 10. Изменение давления на выходе из НТГГ

время, с — эксперимент — расчёт Рис. 11. Изменение температуры на выходе из НТГГ

время, с

-эксперимент — расчёт

Рис. 12. Изменение давления на выходе из НТГГ

Во всех вариантах после первой секунды процесса разброс значений температуры находится в пределах 15 К. Качественные и количественные различия изменения температуры наблюдаются на начальной стадии работы НТГГ.

На рис. 13, 14 показано изменение температуры и расхода газовой смеси по радиусу кассет в различные моменты времени. Расчёты показывают, что различие в параметрах кассет проявляется в пределах 0,1 с после начала. После этого значения расходов, температур, давлений в кассетах практически одинаковы. Данное обстоятельство позволяет уменьшить объём вычислительной работы, считая кассеты идентичными.

В расчётах исследовалось влияние теплообмена газа с элементами конструкции НТГГ, диаметров отверстий в стенке газохода, диаметров входа в газоход и патрубка. Расчёты показывают, что пренебрежение теплообменом увеличивает температуру газа на 10...12 К. Время расчёта при этом уменьшается незначительно.

800 W 700

а 600

а 500

(D

g 400

<D

Н 300 200

1 3 5 7 9 11 13 15

узлы сетки по радиусу

Рис. 13. Изменение температуры по радиусу кассет в момент времени 0,01 с

0,07 -------

- 0q0Q6; Г____

0-------

1 3 5 7 9 11 13 15

узлы сетки по радиусу Рис. 14. Изменение расхода по радиусу кассет в момент времени 0,01 с

Установлено, что изменение диаметров отверстий в стенке газохода, на входе в газоход и выходе из коллектора практически не влияют на величину температуры. Однако эти изменения очень сильно влияют на характер и величину изменения давления как в КС, так и на выходе из НТГГ. Увеличение диаметра отверстий на 25 % приводят к уменьшению давления на 45 %. Изменение давления в КС существенно зависит от диаметров отверстий в стенке газохода и диаметра входа в газоход, что приводит к увеличению времени горения топлива на 70 %. Изменение диаметра выходного патрубка практически не влияет на давление в КС и время работы НТГГ.

1400 1200

<а 1000

а

fr

5 800

а

s 600

(D

400 200

123

456

7

5 6 7 t—t—1-

--^

12 3 4

время, с — 0,1И8 ; — 0,5И ; — И Рис. 15. Влияние энтальпии разложения на температуру на выходе из КО

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

Данная ситуация была смоделирована путём применения охладителя, скорость разложения которого равна скорости разложения карбоната аммония (КА), но с энтальпией разложения составляющей 50 и 10 % от энтальпии разложения КА.

Уменьшение энтальпии разложения в два раза не привело к заметному изменению температуры на выходе из КО. Однако изменение энтальпии разложения в 10 раз приводит к росту температуры в начале КО уже через 0,2 с от начала процесса. В результате повышается разность потока и температуры разложения по длине КО (рис. 16), что вызывает увеличение притока массы охладителя за счёт ускорения процесса разложения (рис. 17). Как следствие ускоряется уменьшение размеров гранул, что приводит в уменьшению массопритока в начале КО и дальнейшему росту температуры. Приток массы продуктов разложения по длине КО растёт, вследствие этого растёт и расход потока по длине КО в гораздо большей степени, чем при высокой энтальпии разложения. Количество разложившегося охладителя увеличивается, а температура на выходе из КО не только не уменьшается, а увеличивается. Между тем в случае 50 % уменьшения энтальпии разложения температура на выходе из КО практически не отличается от случая 100 % величины энтальпии разложения. Таким образом, возможно применение охладителей с высокой скоростью разложения (как у КА) и с меньшей в два раза энтальпией разложения. При этом следует соответствующим образом увеличить начальную массу охладителя.

1600

1400

1200

200

5 7 9

узлы сетки по радиусу

—0,1HS

11

H.

м., —0,5HS ; — ns Рис. 16. Изменение температуры по радиусу КО

13

15

0,12

0,10

« 0,08

0,06

0,04

0,02

0,00

7 9

узлы сетки по радиусу

11

13

15

0,1H.

0,5HS ;

Hs

Рис. 17. Изменение массопритока по радиусу КО

1

3

1

3

5

Выводы

1. Разработанные математические модели на основе одномерных уравнений неразрывности, количества движения, энергии качественно и количественно правильно описывают процессы в НТГГ с камерой охлаждения.

2. Разработан экономичный численный метод расчёта, обладающий высокой вычислительной производительностью.

3. Результаты расчёта показали, что значения температуры на выходе из КО качественно и количественно согласуются с имеющимися экспериментальными данными. Исследовано влияние геометрических параметров НТГГ на характер протекания процессов тепло- и массообмена. Отмечена высокая чувствительность результатов расчёта к изменению геометрических характеристик. Особенно велико влияние геометрических параметров на изменение давления.

4. Получена информация о характере изменения температуры, расхода, давления по длине элементов НТГГ и во времени, что позволяет глубже оценить влияние конструктивных и режимных факторов на процессы течения и теплообмена.

5. Выявлено слабое влияние учёта теплообмена газа с элементами конструкции в математической модели, что объясняется малым временем работы.

Литература

1. Production of the fire extinguishing mixture by solid propellant propulsion / A.I. Karpov, A.Y. Lesthev, A.M. Lipanov, G.A. Lesthev // Journal of less prevention in the process industries. - 2013. -Vol. 26. - P. 338-343.

2. Srinivasan, R. Experiments on thermal response of low aspect ratio packed beds at high Reynolds numbers with varying inflow temperatures / R. Srinivasan, B.N. Raghunandan // Experimental thermal and fluid science. - 2012. - Vol. 44. - P. 323-333.

3. Krishnan, S. Experimental investigation of erosive burning of composite propellants under supersonic crossflows / S. Krishnan, K.K. Rajesh // International Journal of Energetic Materials and Chemical Propulsion. - 2002. - Vol. 5. - Iss. 1-6. - P. 316-325.

4. Cheng, H.G. The design and numerical simulation study of gas generator sealing device / H.G. Cheng, J.S. Zhou, T.B. Wang, C.G. Dong // Proceedings of the third international conference on mechanical engineering and mechanics. - 2009. - Vol. 1, 2. - P. 1213-1217.

5. Thurbochaged solid propellant ramjet for actical missile / S. Yang, G.Q. He, Y. Liu, J. Li // Mechanical engineering and materials. - 2012. - Vol. 152-154. - P. 204-209.

6. Hong, M. Experimental correction of combustion gas properties ofAN-based composite solid propellant used for turbo-pump starter / M. Hong // Aerospace science and technology. - 2012. - Vol. 16. -P. 56-60.

7. Oserov, A. Analytical modelling of the gas generator frequency response in hybtid rocket boosters / A. Ozerov, B. Natan, A. Gany //ACTA Astronautics. - 1986. - Vol. 39, No. 8. - P. 589-598.

8. Kim, A. Explosion suppression with hybrid gas generator system / A. Kim, G. Crampton //Progress in safety science and technjljgy. - 2008. - Vol. 7. - P. 891-895.

9. A functional on-clip pressure generator using solod chemical propellant lab-on-a-chip / G.C. Hong, S. Murugesan, S. Kim et al. //Lab on a Chip. - 2003. - Vol. 3. - P. 281-286.

10. Engelen, K. Chemical formulation of solod propellant for specific gas generators / K. Engelen, M.H. Lefebvre, J. De Ruyek // Twenty-fourth international pyrotechnics seminar, 1998. - P. 203-216.

11. Низкотемпературные газогенераторы на твёрдом топливе / С.Д. Ваулин, А.М. Калинкин, С.Г. Ковин и др.; под общ. ред. акад. А.М. Липанова. - Ижевск: ИПМ УрО РАН, 2006. - 236 с.

12. Теоретическое и экспериментальное исследование низкотемпературных газогенераторов /Д.Д. Аксёненко, С.Д. Ваулин, В.Г. Зезин и др.; под общ. ред. акад. А.М. Липанова. - Ижевск: ИПМ УрО РАН, 2008. - 264 с.

13. Численный эксперимент в теории РДТТ/А.М. Липанов, В.П. Бобрышев А.В. Алиев и др. -Екатеринбург: ИПМ УИФ Наука, 1994. - 303 с.

14. Вукалович, М.П. Таблицы теплофизических свойств воды и водяного пара /М.П. Вукало-вич, С.Л. Ривкин, А.А. Александров. - М. : Изд-во стандартов, 1969. - 408 с.

15. Кириллов, В.В. Расчётно-теоретическое исследование процессов тепло- и массообмена в низкотемпературных газогенераторах / В.В. Кириллов // Хим. физика и мезоскопия. - 2008. -Т. 10, № 4. - С. 428-435.

16. Самарский, А.А. Методы решения сеточных уравнений / А.А. Самарский, Е.С. Николаев. -М. : Наука, 1978. - 592 с.

Кириллов Валерий Владимирович. Доктор технических наук, профессор кафедры «Двигатели летательных аппаратов», Южно-Уральский государственный университет, г. Челябинск, valery.v.kirillov@gmail.com.

Поступила в редакцию 7 мая 2015 г.

SIMULATION OF LOW-TEMPERATURE OPERATION MODE GAS GENERATOR

V.V. Kirillov, South Ural State University, Chelyabinsk, Russian Federation, valery.v.kirillov@gmail.com

Low-temperature gas generator (LTGG) for producing gas with a temperature of 350 ... 450 K. It consists of two main parts - a combustion chamber with a charge of solid propellant and cooling chamber with granules of the solid coolant. The products of combustion of solid propellant have a high temperature that does not allow to use them in a variety of technical devices. Cooling gas is the result of its interaction with the solid granules cooler on the thermal decomposition of which consumes part of the internal energy of the gas. To describe the processes of heat and mass transfer in LTGG developed a mathematical model based on the one-dimensional equations of continuity, momentum, energy and mass balance of the components of the gas mixture. The gas mixture is composed of products of combustion of solid propellant, degradation products and the air cooler. The mathematical model takes into account the ignition and combustion of the fuel within the gas mixture in LTGG, heat exchange with the structural elements LTGG and kinetics of degradation of pellets cooler. To implement the mathematical model developed numerical method for solving the boundary value problem based on the finite difference method. The system of equations is approximated by implicit difference schemes. For the solution of nonlinear differential equations applied Newton's method. A boundary value problem was solved orthogonal factorization. The calculation results are compared with experimental data. The influence of constructive and regime parameters on the change in temperature, pressure, flow, speed expansion cooler.

Keywords: low-temperature gas generator, mathematical model, numerical method.

References

1. Karpov A.I., Lesthev A.Y., Lipanov A.M., Lesthev G.A. Production of the Fire Extinguishing Mixture by Solid Propellant Propulsion. Journal of less prevention in the process industries, 2013, vol. 26, pp. 338-343. DOI: 10.1016/j.jlp.2011.10.007

2. Srinivasan R., Raghunandan B.N. Experiments on Thermal Response of Low Aspect Ratio Packed Beds at High Reynolds Numbers with Varying Inflow Temperatures. Experimental thermal and fluid science, 2012, vol. 44, pp. 323-333. DOI:10.1016/j.expthermflusci.2012.07.004

3. Krishnan S., Rajesh K.K. Experimental Investigation of Erosive Burning of Composite Propel-lants Under Supersonic Crossflows. International Journal of Energetic Materials and Chemical Propulsion, 2002, vol. 5, iss. 1-6, pp. 316-325. DOI: 10.1615/IntJEnergeticMaterialsChemProp.v5.i1-6.340

4. Cheng H.G., Zhou J.S., Wang T.B., Dong C.G. The Design and Numerical Simulation Study of Gas Generator Sealing Device. Proceedings of the third international conference on mechanical engineering and mechanics, 2009, vol. 1, 2, pp. 1213-1217.

5. Yang S., He G.Q., LiuY., Li J. Thurbochaged Solid Propellant Ramjet for actical Missile. Mechanical engineering and materials, 2012, vol. 152-154, pp. 204-209.

6. Hong M. Experimental Correction of Combustion Gas Properties of AN-based Composite Solid Propellant Used for Turbo-Pump Starter. Aerospace science and technology, 2012, vol. 16, pp. 56-60. D01:10.1016/j.ast.2011.03.001

7. Oserov A., Natan B., Gany A. Analytical Modelling of the Gas Generator Frequency Response in Hybtid Rocket Boosters. ACTA Astronautics, 1886, vol. 39, no. 8, pp. 589-598.

8. Kim A., Crampton G. Explosion Suppression with Hybrid Gas Generator System. Progress in safety science and technology, 2008, vol. 7, pp. 891-895.

9. Hong G.C., Murugesan S., Kim S., Beacage G., Choi G., Ahn C.H. A Functional On-Clip Pressure Generator Using Solod Chemical Propellant Lab-on-a-chip. Lab on a Chip, 2003, vol. 3, pp. 281-286.

10. Engelen K., Lefebvre M.H., De Ruyek J. Chemical Formulation of Solod Propellant for Specific Gas Generators. Twenty-Fourth international pyrotechnics seminar, 1998, pp. 203-216.

11. Vaulin S.D., Kalinkin A.M., Kovin S.G., Mariash V.I., Simonov E.A., Feofilaktov V.I. Niz-kotemperaturnye gazogeneratory na tverdom toplive [Low-Temperature Solid Fuel Gas Generator]. Izshevsk, IAMScience, 2006. 236 p.

12. Aksionenko D.D., Vaulin S.D., Zezin V.G., Kirillov V.V., Mazitov R.Z., Mariash V.I., Sav-chenkov V.I., Feofilaktov V.I. Teoreticheskoe i eksperimental'noe issledovanie nizkotemperaturnykh gazogeneratorov [Theoretical and Experimental Study of Low-Temperature Gas Generators]. Izshevsk, IAM Science, 2008. 264 p.

13. Lipanov A.M., Bobryshev V.P., Aliev A.V., Spiridonov F.F., Lisitsa V.D. Chislenny experiment v teorii RDTT [Calculation Experiment in the RESF Theory]. Ecaterinburg, IAM Science, 1994. 303 p.

14. Vukalovich M.P., Rivkin S.L., Aleksandrov A.A. Tablitsy teplofizicheskikh svoystv vody i vo-dyanogopara [Table of the Water and Water Steam Property]. Moscow, 1969. 408 p.

15. Kirillov V.V. [Settlement and Theoretical Study of the Processes of Heat and Mass Transfer in Low-Temperature Gas Generators]. Khimicheskaya fizika i mezoskopia [Chemistry physics and mesos-copy], 2008, vol. 10, no. 4, pp. 428-435.

16. Samarsky A.A., Nikolaev E.S. Metody resheniya setochnykh uravneny [Calculations Methods for Difference Equations]. Moscow, 1978. 601 p.

Received 7 May 2015

БИБЛИОГРАФИЧЕСКОЕ ОПИСАНИЕ СТАТЬИ

REFERENCE TO ARTICLE

Кириллов, В.В. Моделирование режимов работы низкотемпературного газогенератора / В.В. Кириллов // Вестник ЮУрГУ. Серия «Машиностроение». - 2015. -Т. 15, № 3. - С. 5-17.

Kirillov V.V. Simulation of Low-Temperature Operation Mode Gas Generator. Bulletin of the South Ural State University. Ser. Mechanical Engineering Industry, 2015, vol. 15, no. 3, pp. 5-17.