11 (1013742), страница 5
Текст из файла (страница 5)
где m4=m3=m2 [кг] в соответствии с законом сохранения массы; H2=H3 [Дж] в соответствии с законом сохранения энергии.
Определяется число Маха на выходе из сопла:
где [м/с] скорость звука в сечении 4-4 сопла.
Тяга двигателя на расчетной высоте у(р4=ру):
где [кг/с] - секундный расход топлива (заданная величина).
Интерполяция табличных данных.
При расчете параметров в сечениях 3-3 и 4-4 двигателя табличные значения температуры могут не совпадать с рассчитанными значениями температуры в рассматриваемых сечениях, поэтому возникает необходимость проводить линейную интерполяцию при определении из таблиц, например, избыточной энтальпии , теплоемкости
, энтропии
.
Интерполяция может проводиться либо графически, либо расчетным путем. Например, при определении некоторого параметра Yx при температуре Tx строится график искомой функции Y=Y(T) по двум значениям температуры T1 и Т2 :
Тогда искомое значение функции Yx=Y(Tx) определяется либо графически по пересечению прямой T=Tx=const с зависимостью Y=Y(T), либо расчетным путем из соотношения:
Пример термодинамического расчета ракетного двигателя с учетом одной химической реакции
Исходные данные:
Высота полета y=5км; давление в камере сгорания ; горючее – метан СН4; окислитель – смесь дикислорода
и диазота
; коэффициент избытка окислителя
, секундный расход топлива
кг/с.
Уравнение химической реакции
CH4+2O2=CO2+2H2O
или CO2+2H2O-CH4-2O2-ON2=0,
Результаты расчета ЖРД с учетом одной химической реакции представлены в следующей таблице.
обозначения | ед. измерения | сечение 2-2 | сечение 3-3 | сечение 4-4 |
T | K | 298,15 | 2921,36 | 1038 |
p | Пa | 11145750 | 11145750 | 54048,3 |
V | м3 | 0,0008383 | 0,0082108 | 0,6018740 |
моль | 1 | 0,0000001 | 0,0000001 | |
моль | 2,225806 | 2,225806 | 2,225806 | |
моль | 0,543478 | 0,043478 | 0,43478 | |
моль | 0 | 0,49999995 | 0,49999995 | |
моль | 0 | 0,9999999 | 0,9999999 | |
n | моль | 3,769284 | 3,769284 | 3,769284 |
кг | 0,032 | 0,00000000032 | 0,00000000032 | |
кг | 0,062323 | 0,062323 | 0,062323 | |
кг | 0,0086944 | 0,0006956 | 0,0006956 | |
кг | 0 | 0,021999997 | 0,021999997 | |
кг | 0 | 0,017999998 | 0,017999998 | |
m | кг | 0,103017 | 0,103017 | 0,103017 |
кг/моль | 0,02733 | 0,02733 | 0,02733 | |
| кг/м3 | 122,8880 | 12,5465 | 0,17116 |
R’ | Дж/(кгK) | 304,272 | 304,272 | 304,272 |
H | Дж | -2790,7 | -2790,7 | -306278,7 |
S | Дж/K | 614,728 | 941,942 | 941,942 |
Cp | Дж/K | 113,602 | 175,173 | 145,895 |
Cv | Дж/K | 82,263 | 143,834 | 114,556 |
к | - | 1,381 | 1,218 | 1,274 |
a | м/c | 353,949 | 1040,278 | 634,148 |
Wc | м/c | - | - | 2427,330 |
Mc | - | - | - | 3,828 |
P | Н | - | - | 84956,554 |
Изложенный метод термодинамического расчета позволяет наряду с расчетом реального процесса сгорания в ЖРД использовать модели смесей идеальных газов для характерных точек цикла для простой, закрытой термодинамической системы, а использование таблиц индивидуальных веществ позволяет учесть зависимость термодинамических свойств реагентов от температуры, что увеличивает точность термодинамических расчетов ракетных двигателей.
11.5. Решение СУХР для систем с одной химической реакцией методом Ньютона
Для решения СУХР системы с одной химической реакцией методом Ньютона будем использовать условие равенства нулю химического сродства в состоянии химического равновесия:
Введем безразмерные величины для химического сродства и свободной энергии Гиббса
Тогда в состоянии химического равновесия имеем: ,
и СУХР в канонической форме записи примет вид:
где - непрерывные функции от 3-х независимых переменных
Введем новые переменные:
. Тогда
Система уравнений (2) решается приближенным итерационным методом Ньютона. При этом нелинейные уравнения (2) заменяются подходящим линейным отображением, используя линеаризацию через разложение функций по формуле Тейлора.
При решении системы линейных уравнений (3) в (к+1) - ом приближении определяются поправки:
Если при к-ом приближении значения известны, то при
(к+1) -ом приближении значения этих параметров будут равны:
Решение (4) оказывается подходящим приближением к решению системы уравнений (3) при условии, что поправки в (к+1) –ом приближении будут менее заданной величины погрешности , т.е.
где погрешность - малое число (10-3,….., 10-6). Итерации продолжаются до выполнения условия (5).
Последовательность расчета:
1) Выбираются величины к-го приближения: .
2) Вычисляется матрица коэффициентов уравнений системы (3) (матрица Якоби): . Эти коэффициенты являются частными производными величин
и
.
-
Решается система уравнений (3) для определения величин
в (к+1)-ом приближении. Процесс итераций продолжается до тех пор, пока изменение величин
(к+1)-го приближения не окажется в пределах погрешности
.
Коэффициенты системы уравнений (3) представлены в таблице для ряда задач, решаемых в авиационно-космической технике.
Таблица коэффициентов системы (3)
Величины в таблице приведены к безразмерному виду путем деления на одну из следующих величин :
,
,
,
.