книги из ГПНТБ / Мелькумов Т.М. Теория двигателей. I. Теория ракетных двигателей. II. Применение ядерной энергии в силовых установках [учебник]
.pdfЭнтальпия продуктов сгорания зависит как от температуры, так и от состава их; состав же продуктов сгорания в свою очередь зависит от температуры. Для веществ, состоящих из элементов Н, С, N и О, исключение составляют продукты сгорания при низких температурах, если а > 1 ; в этом случае температура не влияет на состав продуктов сгорания.
В общем случае определение температуры сгорания следует . вести одновременно с расчетом состава продуктов сгорания в ка мере. Последовательность расчета состава и температуры продук тов сгорания следующая. Исходными данными для расчета явля ются давление р к* в камере сгорания и химический состав топли ва. Зная состав топлива, определяют его элементарный состав и энтальпию /о при начальной температуре. Далее задаются тремя значениями температур ( Т / , Tz" н Т г'") продуктов сгорания и для каждой температуры производят определение состава про дуктов сгорания. Зная состав продуктов сгорания, определяют их энтальпию для трех выбранных значений температур Г, I" и V" и строят график зависимости / от температуры'(фиг. 8.3).
<1> и г. 8.3. К определению температу |
Фиг. 8.4. К определению молекуляр |
ры сгорания |
ного веса |
Пользуясь соотношением (8.31) и графиком фиг. 8.3, нетруд но найти Tz i Для определения состава продуктов сгорания при найденном значении температуры нет необходимости производить вновь решение системы уравнений, описывающей состав продук тов сгорания, поскольку он может быть найден графической ин терполяцией. Таким же образом найти средний (кажущийся) мо лекулярный вес газов в камере сгорания р.к и газовую постоян ную (фиг. 8.4).
Определение состава и температуры продуктов сгорания в об щем случае связано с решением системы из десяти и более урав нений, что представляет собой довольно трудоемкий процесс. Ис пользование электронных вычислительных машин для этих целей существенно упрощает задачу.
201
§ 8.9. ТЕРМОДИНАМИЧЕСКИЙ РАСЧЕТ ПРОЦЕССА
|
ИСТЕЧЕНИЯ ИЗ СОПЛА |
Как отмечалось |
выше, расчет процесса в сопле проводят |
в предположении его |
изоэнтропичности и затем вводят необходи |
мые поправки для определения действительных величин. |
газа |
При изоэнтропическом процессе в сопле энтропия |
|
в каждом сечении сопла одинакова. |
|
Поэтому |
(8.36) |
5 К= 5С, |
здесь 5К и 5С — энтропия газа в конце камеры сгорания и в вы ходном сечении сопла.
Уравнение (8.36) является исходным для расчета.
Для расчета изоэнтропического процесса в сопле необходимо
знать значения энтропии |
различных соединений, входящих в со |
|
||||||||
став продуктов сгорания. Значения энтропий приводятся в специ |
|
|||||||||
альных таблицах при стандартном давлении. |
|
|
|
|||||||
|
Из термодинамики для |
идеального |
газа известно |
|
|
|||||
|
|
|
|
( d S \ |
AR |
|
|
|
|
|
|
|
|
|
др /т |
У |
|
|
|
||
откуда |
|
|
р |
|
|
|
|
|||
|
|
|
|
ккал |
|
|
|
|||
|
|
|
5 = 5° - |
1,986 In — |
|
|
|
|||
|
|
|
моль град |
|
|
|||||
|
|
|
|
|
р0 |
|
|
|||
где 5° — энтропия при стандартном давлении ро; |
|
|
||||||||
|
S — энтропия при данном давлении. |
|
|
|
|
|||||
Поскольку р0= 1, то |
|
|
|
|
|
|
|
|||
|
|
|
5 = S ° ~ 1,986 Inр. |
|
|
|
||||
ных |
Если известен состав |
продуктов сгорания и энтропии |
отдель |
|
||||||
газов, |
то энтропия |
смеси газов определится следующим |
|
|||||||
образом: |
|
|
|
|
|
|
|
|
||
|
|
5 = ---- -----У. р; (5г° - 1,986 In |
|
ККйЛ■ . |
(8.37) |
|
||||
|
|
|
Немцем |
|
|
|
|
кг град |
|
|
|
Обычно |
заданными при |
расчете являются |
состав топлива, |
|
|||||
давление в |
камере рк* и давление в выходном |
сечении сопла Не |
|
|||||||
основной задачей расчета является определение скорости исте |
|
|||||||||
чения, |
а также проходных сечений сопла. |
|
|
|
|
|||||
|
Расчет изоэнтропического истечения производится в следую |
|
||||||||
щей последовательности: |
|
|
|
|
|
давле |
|
|||
|
1) |
для определенного состава топлива при заданном |
, |
|||||||
нии находят состав и температуру газов в камере сгорания |
(§8.8); |
|||||||||
|
2) |
по известному составу и температуре вычисляют |
энтро- |
/ |
||||||
пию газов в камере: |
|
|
|
|
|
|
|
|||
|
|
|
SK= |
I £ |
Pi (5/° - |
1’986 In P i ) , |
|
|
||
|
|
|
Y-rP * |
|
|
|
|
|
|
|
где |
S t берется из таблиц по температуре |
Тг ; |
|
|
|
2 0 2
|
3) для заданного давления /^ выбирают несколько (обычно |
|||||||||||
три) |
значений |
температур (Ти Т2 и Т3) в |
области |
ожидаемой |
||||||||
температуры газов в выходном |
сечении сопла и для каждого зна |
|||||||||||
чения находят состав продуктов сгорания методами, |
изложенны |
|||||||||||
ми в данной главе; |
далее |
для каждого значения Т по |
формуле- |
|||||||||
(8.37) |
находят энтропию и строят график S=>f(T) |
(фиг. 8.5); |
||||||||||
|
4) |
из условия |
(8.36) по известной величине SK находят тем |
|||||||||
пературу |
Тс„д |
— температу |
|
|
|
|
|
|||||
ру |
газов |
идеального процесса |
|
|
|
|
|
|||||
в |
выходном |
сечении |
|
сопла |
|
|
|
|
|
|||
(фиг. 8.5).; для найденного |
|
|
|
|
|
|||||||
значения |
температуры |
опреде |
|
|
|
|
|
|||||
ляют |
состав газов и |
энталь |
|
|
|
|
|
|||||
пию смеси |
|
|
|
|
|
|
|
|
|
|||
|
/с |
|
1 |
|
ккал |
|
|
|
|
|
||
|
|
|
|
|
кг |
|
|
|
|
Г |
||
|
|
|
|
|
|
|
|
|
|
|||
|
|
, |
, |
|
|
|
(8.38) |
|
Фиг. 8.5. |
К определению |
темпера- |
|
|
|
|
|
|
|
|
||||||
здесь / 1 кал/кг для каждого га- |
|
ТурЫ газа |
в выходном |
сечении сопла |
||||||||
за |
берется при |
температуре |
|
|
|
|
|
|||||
Т |
• |
|
|
|
|
|
|
|
|
|
|
|
1 с ид> |
определяют |
идеальную скорость истечения; эта величина |
||||||||||
|
5) |
|||||||||||
может быть найдена из уравнения сохранения энергии, |
написан |
|||||||||||
ного для входного и выходного сечений сопла: |
|
|
||||||||||
|
|
|
|
|
/к * = h + |
А |
|
|
|
|||
|
|
|
|
|
|
|
|
|
ig |
|
|
|
|
|
|
|
|
®«д |
|
|
|
|
М |
|
|
|
|
|
|
|
|
|
|
|
сек ’ |
|
|
|
|
Учитывая, что |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
||||
. в итоге получим |
|
|
‘ к |
— |
■'O’ |
|
|
|
||||
|
|
|
_____ |
|
|
|
||||||
|
|
|
|
|
wHR= |
|
91,531//„ — /с м/сек. |
|
(8.39)- |
|||
|
- Идеальная удельная тяга на расчетном режиме работы сопла. |
|||||||||||
(Рс=Рн) равна |
Рул нд = |
|
_____ |
|
|
(8.40) |
||||||
|
|
|
|
9,33 1//0 — /с кг сек/кг. |
|
|||||||
|
Площадь выходного сечения сопла для изоэнтропического ис |
|||||||||||
течения найдется, если известен расход газа: |
|
|
||||||||||
|
|
|
|
|
|
|
FС И Д |
|
С?а сек |
|
|
(8.41) |
|
|
|
|
|
|
|
|
Тс ид ® нд |
|
|
||
Здесь |
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
= |
|
Рс |
|
|
|
|
|
|
|
|
|
|
|
|
Т’сндЯ |
|
|
|
203
Для определения площади сечения в произвольном месте по длине сопла необходимо знание р и Т вдоль сопла. Точное опреде ление этих величин возможно по изложенной выше методике. За дача облегчается, если заменить действительный процесс расши рения газов в сопле адиабатой, имеющей средний постоянный по казатель, хотя решение в этом случае получается приближенным. Для этой цели можно воспользоваться известным уравнением для определения среднего показателя k:
k |
lgjOc — IgA,* |
(8.42) |
|
Ig T 'cu -lg T 'z |
|
|
|
Средний показатель изоэнтропического процесса, идущего с изменением состава и теплоемкости газов, меньше показателя адиабаты с постоянными составом и теплоемкостью газов. Это ес тественно, так как процесс равновесного расширения продуктов сгорания вследствие рекомбинации сопровождается превращени ем химической энергии в тепловую, что эквивалентно подводу те
пла к газу постоянного состава. Чем сильнее диссоциирован |
газ |
||
перед соплом, |
тем при прочих |
равных условиях меньше средний |
|
показатель. |
|
|
так: |
Критическое сечение сопла в данном случае найдется |
|||
|
Fкр нд |
(8.43) |
|
где |
определяется по данным термодинамического |
расчета.
Перейдем к определению действительных величин реального процесса течения.
Скорость истечения и удельная тяга на расчетном режиме со шла согласно вышеуказанному найдутся так:
Р у , = ?С фк Р у , ид-
Средний показатель действительного процесса должен отли чаться от среднего показателя изоэнтропы k вследствие влияния гидравлических потерь, теплообмена со стенками и неоднородно сти состава по сечению потока. Однако практически это влияние до сих пор уловить не удалось. Отдельные из указанных факто ров влияют в разную сторону на величину показателя процесса (например, теплоотдача в стенки и гидравлические потери). По этому при определении геометрических размеров будем считать, что показатель действительного процесса равен среднему пока зателю изоэнтропы k.
.204
Температура газа в выходном сечении будет равна:
следовательно, |
|
= |
|
ф 2 Т |
|
|
||
Т |
|
|
(8.45) . |
|||||
1с |
|
|
|
тк |
лс ид* |
|||
Площадь выходного сечения сопла |
|
|||||||
|
|
|
|
|
О и сек |
|
(8.46) |
|
|
|
|
|
|
7с |
|
||
где |
|
|
|
|
|
|
||
|
|
|
|
|
Рс |
|
|
|
|
|
1с |
1 |
|
|
|||
|
|
|
|
г> |
|
|
||
или |
|
|
|
|
|
|
|
|
|
|
|
F |
Чк |
|
|
||
|
|
|
|
|
(8.460 |
|||
|
|
|
|
* С ид |
|
|||
|
|
|
|
|
|
Тс |
|
|
Площадь критического сечения |
сопла найдется из форму- |
|||||||
лы (7.5). |
|
|
|
|
|
|
|
|
р |
|
|
|
|
Os сек Тк Рид |
(8.47) |
||
ГКР |
|
|
|
_ * |
||||
или |
|
|
|
|
|
Рк |
|
|
|
|
— F |
|
ш |
|
|||
F |
кр |
|
(8.47') |
|||||
1 |
|
|
л |
кр ид тк* |
||||
На основании многочисленных расчетов установлены значе |
||||||||
ния показателя изоэнтропы к |
для |
многих топлив |
(см. ниже, |
§ 8.10). Используя эти значения показателя к, можно произве
сти |
приближенный расчет процесса истечения, если известны |
|||
параметры газа |
(р К*, Тг, R) перед соплом. Скорость истечения |
|||
в этом случае |
находят по формуле (2.28), а температуру |
газа |
||
Тс |
— по |
формуле (8.44). Проходные сечения определяют |
по |
|
формулам |
(8.46) я (8.47). |
|
||
|
При пользовании / —5-диаграммой для диссоциированных |
продуктов сгорания объем расчетов значительно сокращается. Од нако необходимо иметь заранее подготовленную диаграмму для данного топлива, с вполне определенным составом, т. е. опреде ленным коэффициентом избытка окислителя. Это предполагает проведение серии предварительных расчетов, целесообразность которых может быть оправдана лишь для стандартных топлив. Одна из таких диаграмм приведена в приложении. На I—5-диа граммах наносят изобары, изохоры,- изотермы и другие кривые. Максимальное значение энтальпии ограничено ее значением для исходных веществ и представляет собой исходную энтальпию то
плива /о- Если давление в камере А ,* известно (фиг. 8 .6 ), то исходная
точка k, соответствующая Состоянию продуктов сгорания в камере-
205
при wK= 0 в предположении отсутствия потерь тепла, |
опреде |
ляется пересечением линии р к*=>const и /0 = const В |
точке к |
с помощью диаграммы можно определить температуру |
газов и |
удельный объем. Линия к—с будет графическим изображением процесса изоэнтропического течения; точка с находится на изоба ре р с, соответствующей давлению на срезе сопла. В точке с мож
но определить параметры |
/с, 7’с ид . Скорость истечения |
оп |
ределяется из уравнения |
(8.39) или снимается с / —S-диаграммы |
|
(см. приложение). |
|
|
Ф и г. 8.6. К использованию / —S -диаграммы
Проходные сечения на выходе из сопла и по его длине опре деляются из уравнения расхода или с помощью построения вспо могательных линий постоянной удельной плотности тока яг/у или
постоянного удельного диаметра й уж= |
—^ R = — (кривые Фанно), |
|
где |
V G i сек |
|
|
|
|
± F |
4 |
1 |
|
(8.48) |
|
Gv сек |
|
ТОУ |
Пользуясь линиями w у= const или а?уд=-const, можно опре делить также критическое сечение сопла (точка «кр») и парамет ры состояния газа в нем. В критическом сечении плотность то ка w у достигает максимального значения, a dyl — минимально-
.го; поэтому кривые dyll = const будут касаться изоэнтропы к ^ с в точке «кр» (фиг. 8.6).
§8.10. РЕЗУЛЬТАТЫ ТЕРМОДИНАМИЧЕСКИХ РАСЧЕТОВ ПРОЦЕССОВ ГОРЕНИЯ И ИСТЕЧЕНИЯ И ИХ АНАЛИЗ
На фиг. 8.7 — 8.16 и в табл. 8.6 представлены результаты тер. модинамических расчетов для ряда ракетных топлив.
На фиг. 8.7 и 8.8 представлены результаты расчета состава продуктов сгорания для двух жидких топлив раздельной подачи.
206
Видно, что с уменьшением коэффициента избытка окислителя уменьшаются доля продуктов полного окисления (из-за недостат ка кислорода) и доля продуктов диссоциации. Уменьшение доли продуктов диссоциации объясняется, с одной стороны, уменьше нием температуры сгорания ввиду понижения теплотворности топлива, а с другой стороны, — увеличением доли продуктов не полного окисления; последние диссоциируют в меньшей степени, чем продукты полного окисления. Рост давления несколько уменьшает долю продуктов диссоциации.
Фиг. 8.7. Зависимость состава про дуктов сгорания от а для топлива:
азотная кислота — керосин:
---------- 2 р .=40 а т а ; --------- |
I р . ~20 а т а |
> |
Фиг. 8.8. Зависимость со става продуктов сгорания. от з для топлива: кисло род—керосин:
-------Е ^ л — 40 а т а \ ------ |
2 ^ = 2 0 а т а |
На фиг. 8.9 — 8.15 показаны зависимости Р удид, Тг, р.к и- k от коэффициента избытка окислителя и давления в камере для некоторых жидких топлив раздельной подачи. При отсутствии диссоциации максимум*температуры Тz должен совпадать с мак
симумом |
теплового эффекта (т. е. * = 1 ) . |
Примером |
является |
||
топливо, |
состоящее из |
окислителя 80 % Нг02 и |
горючего 50% |
||
N2H4H2O -f- 50% СН3ОН |
(фиг.8.9). При |
наличии |
диссоциации |
||
максимум температуры |
смещается обычно в сторону |
а < 1. Не |
которое смещение максимального Значения температуры в сторо ну богатых смесей (фиг. 8.10—8.12) объясняется понижением степени диссоциации и, следовательно, уменьшением доли тепла, идущей на диссоциацию при уменьшении коэффициента избытка окислителя.
207
PyJ.tffl. « г е # / м
= Рн — ' a-та. Окислитель — R0% HsOs; горючее — 50% N2 H4 H.O +
+ 50% СНзОН
h
28
?Ч
20
10
Фиг. 8.10. Зависимость Т2, РуДШ1, рк и А от а при
рк* = 40 ата и рс = /> я = 1 ата.
Окислитель — HNOs; горючее — керосин .
Фиг. 8.11. Зависимость |
Т., Р удпд, |
рк и k от |
|
а при р к* = |
40 ата |
и р с — р н = ] ата. |
|
Окислитель — |
жидкий |
кислород; |
горючее — |
керосин
«
Як
|
|
|
2k |
|
|
|
п |
|
|
|
20 |
|
|
|
(3 |
|
|
|
16 |
0 Л |
0,6 |
0,8 |
1,0 Ы |
Фиг. 8.12. Зависимость Тг, РудИД1 рк |
|||
от а при рк* |
ч0 |
ата и р с = р ц *=-1 ата. |
Окислитель — жидкий кислород; горючее — несимметричный диметилгидразид
14. т. М. Мелькумов, Н. И. Мелик-Пашаев |
20Э |
Увеличение давления несколько повышает температуру вслед
ствие |
понижения степени диссоциации продуктов сгорания |
|
(фиг. |
8.13). Чем |
менее диссоциированы продукты сгорания, тем |
слабее влияние давления на температуру. |
||
Влияние а. |
на средний молекулярный вес определяется изме |
нением состава продуктов сгорания. С уменьшением а растет до-
.ля более легких продуктов неполного окисления, вследствие чего
падает. |
С ростом давления молекулярный вес несколько |
ра |
||||||
стет (фиг. 8.14) |
из^за уменьшения доли продуктов диссоциации. |
|||||||
7 > с |
|
|
|
|
|
|
|
|
|
|
|
|
Г* |
4 |
|
2 |
|
|
|
|
|
2 |
|
4 |
|
|
|
|
|
|
23 |
|
|
3 |
|
|
|
|
|
П |
|
|
|
|
|
|
|
|
21 |
|
|
|
|
|
|
|
|
П |
|
<v> |
1 |
|
|
|
|
|
21 |
|
|
|
|
|
20 |
ОД 60 |
Р*ата |
20 |
40 60 |
Pfvma |
|
|
Фиг. |
8.13. Зависимость |
Фиг. |
8.14. |
Зависи |
|
|||
Т. |
от давления в камере |
мость iJ;t от |
давления |
|
||||
|
|
(<*=0,8): |
|
|
в камере (<*=0,8): |
|
||
1-80% H..O.J —50% |
N.H,H..O -т- |
/ —80%Н3О—50%Ы..Н,Н.,О+ |
|
|||||
+50% СН'зОН; 2 —HNOa —ке- |
-г 50%СН3ОН; 2 - HNO,— |
|
||||||
росик; 3—жидкий кислород— |
—керосин; 3—жидкий кис |
|
||||||
—С^Н^ОН; ‘/—жидкий кисло |
лород—C.jH-iOH; ^-жидкий |
|
||||||
|
|
род—керосин |
|
кислород—керосин |
|
|||
Средний |
показатель изоэнтропы k |
зависит от а, так |
как |
меняются состав и температура продуктов сгорания. С уменьше нием й (при й < 1) уменьшается температура, а следовательно, и степень диссоциации, что ведет к уменьшению теплоемкости и поэтому к увеличению показателя k. Уменьшение коэффициента избытка окислителя приводит к росту доли газов с' Меньшим чис лом атомов, что также способствует увеличению k. С ростом дав ления в камере величина k несколько растет ввиду уменьшения степени диссоциации. При данном/?к*с увеличением степени пони-
жения давления в сопле —51— (т. е. с уменьшением рс) величи не
на k также несколько растет, так как при этом уменьшается вли яние диссоциации на процесс расширения газа в сопле.
Изменение удельной тяги (или, что то же, скорости истече ния) с изменением а определяется влиянием коэффициента из
2 1 0