Расчет равновесного состава и температуры рабочего тела в камере сгорания и в выходном сечении сопла

Расчёт ведём для 4 температур– 3200, 3400, 3600, 3800 К. При 3200 К справочные значения констант равновесия определяются из таблиц:

К₁ = 11,4025, К₂ = 9,9108, К₃ = 12,7616, К₄ = 19,3274.

Алгоритм метода Ньютона состоит в последовательном логарифмировании и линеаризации всех уравнений системы.

Логарифмирование и линеаризация 1-го уравнения.

;

,

где на i -ом итерационном шаге

;

; ;

.

Второе уравнение:

;

,

где ;

.

Третье уравнение:

;

, где

.

Четвёртое уравнение:

;

, где

.

Пятое уравнение:

логарифмирование

– .

Выражения для частных производных от левой части этого уравнения по логарифмам всех независимых переменных:

;

;

;

;

.

После линеаризации пятое уравнение:

,

где

.

Шестое уравнение:

.

Линеаризация проводится аналогично пятому уравнению:

,

где

.

Седьмым замыкающим является уравнение закона Дальтона

.

линеаризация:

,

где

.

Таким образом, проводя линеаризацию всех исходных уравнений, получим линейную систему алгебраических уравнений, в которой в качестве неизвестных фигурируют параметры (V – парциальные давления и число молей исходного топлива).

Таким образом, на i -ом итерационном шаге для определения поправок надо решить систему линейных алгебраических уравнений (СЛАУ) , где матрица коэффициентов СЛАУ имеет вид:

.

Вектор неизвестных :

.

Вектор правых частей :

Для использования метода Ньютона нужно задать начальные приближения и . Начальные приближения парциальных давлений выберем из условия их равенства, т. е:

Начальное приближение количества молей топлива находим по формуле: .

Для начального приближения найдём коэффициенты матрицы А и вектор правых частей В:

Т.о., при начальном приближении на первом итерационном шаге требуется решить СЛАУ:

Итерационный процесс продолжается до достижения заданной точности решения. В результате для четырех значений температур 3200К, 3400К, 3600К, 3800К находим состав ПС, который задается парциальными давлениями.

По рассчитанному составу ПС для четырех температур находим энтальпии ПС по формуле:

,

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

Для того, чтобы перевести ккал в Кдж, надо полученное значение умножить на коэффициент 4,187.

Энтальпия ПС должна равняться энтальпии топлива,. Поэтому температура в камере сгорания лежит в интервале от 3600 К до 3800 К. Действительное значение температуры можно определить графически или, используя формулу линейной интерполяции. По формуле линейной интерполяции на основании уравнения сохранения энергии в форме баланса энтальпий:

.

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

;

Определим суммарное значение давления в камере сгорания, т.к. использование линейной интерполяции при определении состава ПС могло привести к погрешностям, связанным с нелинейными функциональными связями :

,

т.е. расчет состава ПС можно считать удовлетворительным.

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

.

Для каждого пробного находятся последовательно , , , , , , . Истинное критическое сечение, соответствующее ему давление и другие параметры будут те, при которых удельная площадь получится минимальной.

Для термодинамического расчета в выходном и критическом сечениях сопла необходимо вычислить значение энтропии на входе в сопло:

.

Здесь размерности ,

Значения стандартных энтропий являются табличными значениями. Используя формулу линейной интерполяции, определим значения стандартных энтропий продуктов сгорания для температуры 3690 К:

.

Используя эти значения, находим суммарное значение энтропии ПС .

Срез сопла.

Составим систему уравнений для среза сопла. Считаем, что на срезе сопла в смеси ПС присутствуют только H₂O и H₂, т.е. система будет состоять только из уравнения сохранения вещества и закона Дальтона.

.

Т.е. состав на срезе от Т_а не меняется.

Найдём значение энтропии ПС на срезе в зависимости от Т_а.

Начнём с Т=1200 К, определим

, и

Аналогично для Т=1300 К

Т.к. значение энтропии ПС на срезе сопла должно равняться энтропии на входе в сопло, то истинное значение температуры Т_а лежит между 1200 К и 1300 К. Используя формулу линейной интерполяции, найдем температуру на срезе Т_а.

Найдем энтальпию ПС на срезе сопла, по найденному значению температуры:

, и .

После нахождения термодинамических параметров в характерных сечениях вычисляют важнейшие теоретические параметры истечения ПС:

Удельный импульс на расчетном режиме:

.

Плотность ПС в КС и на срезе:

, , , , ,

Средний показатель изоэнтропы расширения.

.

Удельная площадь среза сопла.

.

Удельный импульс в пустоте.

.

2. Расчет критического сечения

Проведем термодинамический расчет состава и температуры ПС для четырех предполагаемых значений давления в критическом сечении: 110 Атм, 115 Атм, 120 Атм, 125 Атм. Расчет проводится по системе уравнений для КС, но температура и действительный состав ПС определяется из условия . Запишем результаты проведенных расчетов.

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

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

;

;

.

Здесь – это непрерывный аналог дискретной функции .

Экстремум этой функции ищется обычным способом: находится производная и решается уравнение . В результате получается экстремальное значение (т.е. ), причем в данном случае это будет точка минимума функции . Таким образом, давление в критическом сечении найдено. Далее необходимо провести термодинамический расчет для этого давления, в результате которого находим:

; ; ; ; ; ; ; .

3. Теоретические параметры истечения

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

Характеристическая скорость:

.

Коэффициент тяги в пустоте:

.

Относительная площадь среза

.

Найдём средние показатели изоэнтропы течения n до критики и после:

КС – критика:

Критика - срез:

Общий: .