ПРИБЛИЖЕННЫЙ РАСЧЕТ СОСТАВА ПРОДУКТОВ СГОРАНИЯ

Щипская В.В, Пятышкин Г.Г.
Донецкий национальный технический университет


Источник: Металлургия глазами молодых - 2011 / Матеріали науково-технічної конференції молодих учених та студентів. - Донецьк, ДонНТУ - 2011, с. **-**.


ПРИБЛИЖЕННЫЙ РАСЧЕТ СОСТАВА ПРОДУКТОВ СГОРАНИЯ

Многие процессы, протекающие внутри рабочего пространства теплового агрегата существенно зависят от состава и свойств продуктов горения топлива, поэтому определение состава продуктов горения является важной задачей. В первом приближении расчет состава продуктов горения сводится к отысканию и решению системы уравнений относительно коэффициентов реакции сгорания топлива, которая в общем виде записывается так:

СaРbOcNd = xCO2 + yCO + zH2 + uH2O + (d/2)*N2 (1)

где СaРbOcNd - так называемая условная формула топлива;

a,b,c,d - числа грамм-атомов соответствующих элементов в условной молекуле топлива;

x,y,z,u,d/2 - числа грамм-молекул соответствующих веществ в продуктах сгорания топлива (коэффициенты реакции).

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

2CO = CO2 + C (2)

H2O+CO = CO2 + H2(3) и реакции диссоциации.

Для давлений и температур, характерных для тепловых агрегатов константа равновесия реакции (2) настолько велика, что эту реакцию можно считать предельно сдвинутой влево, т.е. не рассматривать при практических расчетах. Реакции диссоциации также не играют существенной роли, так что остается принять во внимание лишь реакцию (3), константа равновесия этой реакции дает дополнительное уравнение, с учетом которого система уравнений для определения коэффициентов реакции горения становится замкнутой:

a = x + y; (баланс по углероду)

b = 2z + 2u; (баланс по водороду)

c = 2x + y+ u; (баланс по кислороду)

K = uy/xz; (константа равновесия)

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

F = |x+y-a| + | 2z+2u-b| + |2x+y+u-c| + |Kxz-uy|. (4)

Алгоритм решения данных задач был реализован на языке программирования Pascal. Прежде чем приступить к интересующему расчету сгорания на более простом примере был изучен способ составления и решения системы уравнений для определения состава продуктов сгорания. Были рассмотрен случай горения водорода при заданном постоянном давлении, при изменяющейся температуре и при изменяющемся давлении, [табл. 1, табл. 2]

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

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

Литература

  1. Синярев Г.Б. Добровольский М.В. Жидкостные ракетные двигатели. - М.: Государственное издательство оборонной промышленности, 1955, 479с
  2. На Цунг-Йен. Вычислительные методы решения прикладных граничных задач. - М.: Мир, 1982. 295 с
  3. ДьяконовВ.П.Справочник по алгоритмам и программам на языке БЕЙСИК для персональных ЭВМ. - М.: Наука,1989,240с