книги / Термодинамическое проектирование баллиститных артиллерийских порохов и твердых ракетных топлив
..pdfMu =39506-05 |
Lt= 11235+00 |
Lt"= 38037+00 |
MM = 22973+02 |
Ср.г=4055+00 |
k.r= 12734+01 |
MM.r =2282+02 |
R.r = 37145+02 |
Z= 89943-02 |
П л = 0 0 0 0 -0 0 |
Bm = 11198+00 |
n = 12078+01 |
W= 18486+04 |
WA= 28565+01 |
FIF*= 52897+01 |
F"= 16675-02 |
/удп= 2083+03
Содержание компонентов - объем, (масс.) доли
Н 10712-08 |
Н2 29648+00 |
Н20 75489-01 |
N2 99088-02 |
NH3 11183-04 |
СО 37704+00 |
С02 13924+00 |
СНЗ 48199-09 |
СН4 72736-02 |
С2Н4 31491-07 |
С2Н6 31003-07 |
СН20 74232-07 |
СН202 40708-07 |
N2C 94555-01 |
HCN 27882-06 |
k*Cu 39943-02 |
1с*СаСОЗ 5 0 0 0 0 -0 2 |
|
3.3.1. |
Размерности выходных данных |
Размерности выходных величин, помещаемых на экран дисплея, и
выводимые документы программного комплекса определяются введенны ми или предполагаемыми по умолчанию директивами: <PrSI, <Prte, <РгМ, <PrP, <PrR, <РгС.
|
Размерности общих характеристик равновесия приведены в таблице. |
||
|
Размерности выходных величин |
|
|
Пара |
Название параметра |
с и |
Техническая |
метр |
|
|
система |
|
|
|
единиц |
1 |
2 |
3 |
4 |
Р |
Давление |
МПа |
кгс/см2 |
Т |
Температура |
к |
к |
V |
Удельный объем |
м3/кг |
м3/кг |
S |
Энтропия |
кДж/(кг*К) |
ккал/(кг*К) |
1 |
2 |
/Энтальпия
иВнутренняя энергия
мОбщее число молей компонентов
Ср |
Удельная теплоемкость при постоян |
|
ном давлении (замороженная) |
кПараметр к= Cp/Cv
АРавновесная скорость звука
Ми |
Коэффициент динамической вязкости |
и |
Коэффициент теплопроводности |
Lt” |
Полный коэффициент теплопровод |
|
ности |
ММСредняя молярная масса
Ср.г |
Удельная теплоемкость газовой фазы |
к.г |
Параметр кг = Cpr/Cvr |
ММ.г |
Средняя молярная масса газовой фазы |
R s |
Газовая постоянная |
Z |
Массовая доля всех конденс. фаз |
Пл |
Плотность исходной смеси |
В т |
Окислительный потенциал |
N |
Показатель процесса расширения |
W |
Скорость потока |
WIA |
Число Маха |
1удп |
Удельная тяга (импульс) в пустоте |
ВРасходный комплекс
Окончание таблицы 3 4
кДж/кг ккал/кг кДж/кг ккал/кг моль/кг моль/кг кДж/(кг-К) ккал/(кг-К)
1 |
1 |
м/с |
м/с |
Нх/м2 |
кгсс/м2 |
Вт/(мК) |
ккал/(м*час-К) |
Вт/(м*К) ккал/(м-час-К) |
|
г/моль |
г/моль |
кДж/(кг-К) |
ккал/(кг-К) |
1 |
1 |
г/моль |
г/моль |
кДж/(кгК) |
кгсм/(кг-К) |
1 |
1 |
кг/куб.м |
г/куб.м |
1 |
1 |
1 |
1 |
м/с |
м/с |
1 |
1 |
м/с |
с |
Печать равновесных концентраций компонентов всех фаз иницииру ется директивами <PrM, <РгР, <PrR или <РгС. Каждая из них обеспечи вает вывод состава в определенных единицах измерения. По умолчанию
предполагается директива <РгМ.
Директива <Prdom обеспечивает вывод концентраций тех индивиду
альных веществ, содержание которых в рабочем теле не менее
0,000001 моль/кг. По директиве <Prfull выводятся концентрации всех ве ществ, участвовавших в расчетах.
В данной работе выходные величины R.г, Z, W и ММ.г далее исполь зуются для расчета следующих характеристик исследуемых систем:
• Rpp - газовой постоянной для смеси продуктов сгорания;
Rnp= Rs (1-Z), [кДж/ (кг-К)] в СИ и [кгс-м/ (кг-К)] в технической системе единиц;
•/ - силы пороха, / = Дпр• Т\
• 1\ - единичного расчетного импульса (удельной тяги при Ра = Рнар), /, = 079,80665 (с);
• Мюпр - среднего молекулярного веса смеси продуктов сгорания (отношения полного веса продуктов сгорания к числу молей газовой фазы),
Мюпр= MM.r/(l-Z) (г/моль);
где ММ.г —отношение веса газовой фазы продуктов к числу молей газовой фазы; ММ - отношение полного веса продуктов к общему числу молей М (газы+конденсаты).
3.3.2. Вывод на экран дисплея Результаты, полученные для каждого равновесного состояния, ото
бражаются на экране дисплея в верхнем окне. После завершения первого этапа расчета пользователю предоставляется возможность просматривать полученные результаты в асинхронном режиме, т. е. не прерывая парал
лельно продолжающихся вычислений (если они предусмотрены входными данными). При переходе к следующему заданию все ранее полученные ре зультаты (но не более 300 последних строк) продолжают оставаться дос тупными для обозрения.
3.3.3. Вывод в дисковый файл
Вывод результатов в дисковый файл, если при настройке программ ного комплекса он был запрошен, осуществляется после определения ха рактеристик равновесия в каждой расчетной точке.
Перед выводом результатов в файл записываются входные данные к расчету в том виде, в каком они были подготовлены в окне встроенного текстового редактора. Файл вывода не изменяется в течение всего сеанса работы с программным комплексом, и в него записываются результаты всех рассчитанных вариантов.
Запись результатов в дисковый файл производится средствами опе рационной системы, которая, в свою очередь, использует режим буфериза ции. Из-за этого может происходить частичная потеря информации на дис ке при прерывании работы программного комплекса до завершения цикла вычислений. Для обеспечения вывода всех полученных данных рекоменду ется завершать работу после возврата в окно текстового редактора.
При необходимости дисковый файл может быть распечатан.4
4. РАСЧЕТНО-ТЕОРЕТИЧЕСКИЕ ИССЛЕДОВАНИЯ ГОРЕНИЯ
ДВОЙНЫХ СМЕСЕЙ ЭЛЕМЕНТАРНЫХ ВЕЩЕСТВ
Основным фактором, влияющим на термодинамические и энергети ческие характеристики продуктов горения, является соотношение между окислительными и окисляющимися элементами в составе пороха. Соотно-
шение между окислительными и окисляющимися элементами - коэффи
циент |
избытка |
окислительных элементов а, определяют как |
|
|
a = (Z(- 4 v i)/(2(+)wiv<), |
здесь |
Y.(~)wivi - |
суммарное число окислительных валентностей в услов |
ной молекуле топлива, равное сумме произведений числа атомов окисли тельных элементов v* на их высшую валентность (w*). Аналогично для
горючего |
- суммарное число окисляющихся валентностей. |
|
В зависимости от величины коэффициента избытка |
окислительных |
|
элементов а изменяется состав продуктов горения. Для |
топлива, состоя |
|
щего из |
элементов H,C,0,N, при низких температурах в составе про |
|
дуктов горения преобладают следующие газы: Н2О, СО2, N2 при а=1, |
||
Н20, С 02, 0 2, N2 при а > 1, Н20, С 02, СО, N2 при а < 1. При значитель |
ном недостатке кислорода могут образовываться, кроме того, низшие уг
леводороды (СН4), а также свободный углерод. |
|
|
|
Следующим важным фактором, |
определяющим состав продуктов |
||
горения, является температура. Состав продуктов горения |
при высоких |
||
температурах отличен от состава при низких температурах из-за |
диссо |
||
циации и ионизации компонентов продуктов горения (рисунок). |
|||
Диссоциация начинается при температуре - 2000 К. |
Для |
топлива, |
|
состоящего из Н, С, О, N, при температуре выше 2000 К в продуктах го |
|||
рения дополнительно появляются газы |
ОН, О, Н, N0, N. |
При очень вы |
соких температурах возможно появление свободного углерода. Появление этих газов является результатом диссоциации сложных соединений на бо
лее простые: |
СО <=> СО + 1/2 О2, Н2 О <=> Н2 + О, Н20 <=> НО + Н. |
||
С |
повышением температуры в результате увеличения внутренней |
||
энергии |
молекул (вращательной, колебательной) |
колебательные движе |
|
ния атомов |
в некоторых из них развиваются |
настолько интенсивно, |
что энергия этих колебаний оказывается достаточной для преодоления действующих между атомами сил притяжения и молекула разрушается. Процесс диссоциации требует затраты части энергии теплового движения молекул, что приводит к уменьшению температуры рабочего тела.
Рис. Зависимость состава продуктов горения от температуры при а = 1
|
Процессу диссоциации сопутствует |
процесс рекомбинации - |
вос |
|||||
соединения атомов и |
осколков молекул |
в |
молекулы. Этот процесс обу |
|||||
словлен эффективными столкновениями, |
в |
которых участвуют необхо |
||||||
димые |
осколки |
молекул, и сопровождается выделением того же количе |
||||||
ства тепла, которое было поглощено при диссоциации. Реакции |
диссоциа |
|||||||
ции |
и рекомбинации обратимы. В определенных условиях скорости |
пря |
||||||
мой |
и |
обратной реакций будут равны, т.е. устанавливается состояние |
||||||
равновесия. Химический состав рабочего |
тела, отвечающий |
этому |
со |
|||||
стоянию, называется |
химически равновесным. |
|
|
|||||
|
Другим важнейшим фактором, определяющим положение рав |
|||||||
новесия, |
является |
давление. Согласно известному принципу |
смещения |
|||||
равновесия, увеличение давления вызывает |
процессы, способствующие |
относительному снижению давления, т.е. процессы рекомбинации, веду щие к снижению числа молекул смеси.
Результатом сильного возбуждения электронов при высокой тем пературе может быть отрыв электрона и превращение атомов в положи тельно заряженный ион. Ионизация продуктов горения химических
топлив ничтожна и ее влияние целесообразно учитывать лишь при расчете некоторых свойств продуктов горения.
Влияние соотношения окислителя с горючим на термодинамические характеристики продуктов их горения исследуют на примере двойных сме сей элементарных веществ: С+О2, Н2+О2, Al+Ог, Mg+Ог, при а =0,3; 0,6;
1,0; 1,3; 1,5. Величина коэффициента а задается следующим образом. На
пример, кислородо-водородная смесь взята в стехиометрическом отноше нии: 2Нг + Ог = 2Н20 , т.е. в реакции принимают участие 2 моль водорода и
1 моль кислорода, в результате образуется 2 моль воды. Суммарное |
число |
|
окислительных валентностей в данном случае |
fv|. = 2 - 2 1 = 4, где 2 - |
|
валентность кислорода и число атомов в молекуле; 1 - число молей. |
Сум |
|
марное число окисляющихся валентностей |
также равно 4 |
(1-валентный водород, 2 атома в молекуле, 2 моль). Таким образом, а = 1 Обратная задача: задан коэффициент а = 0,3. Определяем количество мо
лей кислорода, приходящихся на 1 моль водорода при заданном а:
N 2 2 |
т.е. N = 0,15. |
0,3 = ■—— |
|
1- 2-1 |
|
Пример подготовки исходных данных для случая С+О2 при а =1: <lNSI<PRSI<PRR<PRDOM<NOION>
/*=300, и=о,
1С|=1, [01=1.
Здесь содержание компонентов задано в мольных долях согласно ре акции С+02 =С02: на 1 моль углерода приходится 1 моль кислорода.
При а =0,6 последняя строка будет выглядеть так: [С]=1,[О]=0,6. Для проведения расчетов необходимо выполнить следующие дейст
вия:
•Вызвать программу расчета термодинамических характеристик. •Ввести директивы к расчету.
•Подготовить и ввести исходные данные: давление и термодинамиские параметры, определяющие условия равновесия.
•Ввести содержание химических элементов в исследуемой системе. Самый простой способ описания химического состава в данной задаче - прямое задание мольного содержания каждого элемента. Названия элемен тов изображаются общепринятыми символами Периодической системы элементов Менделеева.
П р и м е р . Исследуется двойная смесь: углерод - кислород при
а = 1.
Исходные данные: <INSI<PRSI<PRR<PRDOM<NOION>
/*=300, и=о,
[С| = 1, [О] = 1- Результаты расчета:
Характеристики равновесия СИ:
Р = 3 0 0 0 0 + 0 3 Т=44221+04 V=43705-02 S=79150+01 1=13112+04
U = 0 0 0 0 - 0 0 М=35660+02 Ср=13517+01 k=12810+01 Ср"=13661+01
k"=12748+01 A=12924+04 M u=1030-03 Lt=1976 + 0 0 |
L t"= 20012+ 00 |
||||
MM=2804+02 |
Cp.r=1351+01 |
k.i-1281+Ol |
MM.p=2804+02 |
||
R.r=2965+03 |
Z= 0 0 0 0 0 - 0 0 |
Пл= 0 0 0 0 0 - 0 0 |
Bm= 63150-03 |
||
Содержание компонентов - мае. доли: |
|
|
|||
О 60390-01 |
02 |
17746+00 |
С 39387-06 |
СО |
41533+00 |
С02 34681+00 |
С 20 |
52819-06 |
С302 33758-08 |
03 |
86131-05 |
По результатам расчетов построить графические зависимости темпе ратуры, состава, среднего молекулярного веса продуктов горения, силы пороха применительно к условиям артиллерийского выстрела и удельного импульса ТРТ применительно к условиям работы ракетного двигателя от величины а. Сила и удельный импульс композиции рассчитываются по формулам, приведенным в разделе 3.3.1.
5. РАСЧЕТНО-ТЕОРЕТИЧЕСКОЕ ИССЛЕДОВАНИЕ ГОРЕНИЯ
ТРОЙНЫХ СМЕСЕЙ ЭЛЕМЕНТАРНЫХ ВЕЩЕСТВ
Расчетно-теоретическое исследование горения тройных смесей эле
ментарных веществ проводят после выполнения исследования горения
двойных смесей. Тройная смесь представляет собой двухкомпонентное го рючее, третьим компонентом является окислитель - кислород. Исследова ния проводят при различном соотношении содержания различных компо нентов горючего и при постоянном коэффициенте избытка окисляющих
элементов а. |
|
|
|
|
|
Задают следующие соотношения компонентов: |
|
|
|||
1-й компонент(моль) |
1 |
1 |
1 |
2 |
3 |
2-й компонент (моль) |
3 |
2 |
1 |
1 |
1 |
В ходе выполнения данной работы рассматривают следующие пары компонентов горючего: Н - С, Н - А1, С - Al, Н - Mg, С - Mg, А1 - Mg.
Мольное содержание окислителя при заданном значении а рассчи тывают по формуле а = (Е(”)w,v,)/(Е( + ) ).
П р и м е р . Для соотношения Н/С=3/1 имеем 3 моль одновалентного водорода по два атома в молекуле 1-3-2 и 1 моль четырехвалентного углерода, т.е. (2(+)v*'I.v,.) = l-3-2 + 4 1 1 = 10.
N моль двухвалентного кислорода (по два атома в молекуле) дают суммарное число окислительных валентностей в условной молекуле топ лива: (L(~)wivi) = N ’2-2 = 4N При а = 1 суммарное число окислитель ных валентностей 4N должно быть равно 10, отсюда N= 2,5. Задавая таким образом мольное содержание каждого компонента, при помощи программ ного комплекса определяют термодинамические и энергетические характе ристики тройных смесей. Для этого с экрана дисплея выписывают сле дующие данные: температуру горения Т, средний молекулярный вес про дуктов горения ММ.г, газовую постоянную для смеси продуктов сгорания
R.r, массовую долю всех конденсированных фаз Z, а также массовую долю основных продуктов горения.
По результатам расчетов строят графические зависимости темпера туры, среднего молекулярного веса и состава продуктов горения, силы и удельного импульса от соотношения компонентов в горючем.
6. ИССЛЕДОВАНИЕ ВЛИЯНИЯ
СТЕПЕНИ ЭТЕРИФИКАЦИИ НИТРОЦЕЛЛЮЛОЗЫ
НА ЭНЕРГЕТИЧЕСКИЕ ХАРАКТЕРИСТИКИ, ТЕМПЕРАТУРУ
И СОСТАВ ПРОДУКТОВ ГОРЕНИЯ
Используя исходные данные табл. 3 приложения, рассчитывают тер модинамические и энергетические характеристики не менее чем для пяти образцов НЦ с различным содержанием азота и с экрана дисплея выписы вают следующие данные: темпералуру горения 71, средний молекулярный вес продуктов горения ММ.г, газовую постоянную для смеси продуктов сгорания R.г, массовую долю всех конденсированных фаз Z, а также мас совую долю СО, СО2, Н2О и N2.
По результатам расчетов строят графические зависимости темпера туры, среднего молекулярного веса и состава продуктов горения, силы или
зо