Математическая модель расчёта термодинамических параметров гомогенной смеси продуктов сгорания углеводородного топлива в термодинамическом цикле газотурбинных двигателей

Тип работы:
Реферат
Предмет:
ТЕХНИЧЕСКИЕ НАУКИ


Узнать стоимость

Детальная информация о работе

Выдержка из работы

МАТЕМАТИЧЕСКАЯ МОДЕЛЬ РАСЧЁТА ТЕРМОДИНАМИЧЕСКИХ ПАРАМЕТРОВ ГОМОГЕННОЙ СМЕСИ ПРОДУКТОВ СГОРАНИЯ УГЛЕВОДОРОДНОГО ТОПЛИВА В ТЕРМОДИНАМИЧЕСКОМ ЦИКЛЕ ГАЗОТУРБИННЫХ ДВИГАТЕЛЕЙ
Болдырев О. И.
ОАО «Научно-производственное предприятие «Мотор»
Уфа, Россия (450 039, г. Уфа, ул. Сельская-Богородская, 2) boldyrevoi@rambler. ru
На основе метода расчётов ракетных двигателей автором предложена математическая модель расчёта термодинамических параметров гомогенной смеси продуктов сгорания углеводородного топлива с учётом термической диссоциации и рекомбинации. Модель построена применительно к рабочим процессам высокотемпературных газотурбинных двигателей на основе констант химического равновесия рассматриваемых реакций. Дополнительно в модели учтён ряд факторов, необходимых для расчёта термодинамического цикла ГТД, а также фактор неравновесности состава смеси в зависимости от времени пребывания рабочего тела в объёме каждого узла двигателя. В таком виде модель применена для расчёта рабочих процессов сжатия, расширения и смешения рабочего тела в узлах двигателя, а также процессов сгорания топлива в основной и форсажной камере сгорания. Применение математической модели в составе системы моделирования рабочих процессов ГТД позволяет проводить термодинамический расчёт ГТД и его узлов с учётом термической диссоциации и рекомбинации.
Ключевые слова: термическая диссоциация, равновесное состояние, продукты сгорания, газотурбинные двигатели.
MATHEMATICAL MODEL OF CALCULATION OF THERMODYNAMIC PARAMETERS OF THE HOMOGENEOUS MIX OF HYDROCARBONIC FUEL COMBUSTION PRODUCTS IN THE THERMODYNAMIC CYCLE OF GAS-TURBINE ENGINES
Boldyrev O.I.
Joint-stock company «Scientific-and-production enterprise& quot-MOTOR"-
Ufa, Russia (450 039, Ufa, Selskaya-Bogorodskayastr, 2) boldyrevoi@rambler. ru
On the basis of a method of calculations of rocket engines the author offers mathematical model of calculation of thermodynamic parametres of a homogeneous mix of combustion products of hydrocarbonic fuel with the account of thermal dissociation and recombination. The model is constructed with reference to working processes of high-temperature gas-turbine engines on the basis of constants of chemical equilibrium of considered reactions. In addition in model a number of the factors necessary for calculation of a thermodynamic cycle of a gas-turbine engine, and also the non-equilibrium mix structure factor depending on propulsive mass residence time in volume of each node of the engine is considered. In such kind the model is applied to calculation of working compression processes, expansion and propulsive mass mixture in engine nodes, and also processes of combustion of fuel combustion chamber and the afterburner. Application of mathematical model as a part of system of modelling of working processes of a gas-turbine engine allows to conduct thermodynamic calculation of a gas-turbine engine and its nodes with the account of thermal dissociation and recombination.
Key words: thermal dissociation, equilibrium condition, combustion products, gas-turbine engines.
Качество математической модели во многом определяет успешность проекта, сроки и материальные затраты на его реализацию. Подходы к созданию математических моделей и принятые в них допущения подвержены изменениям в процессе развития газотурбинных двигателей (ГТД), как в связи с повышением уровня параметров, так и возможно-
стей вычислительных средств. Проектирование эксплуатируемых в настоящее время авиационных ГТД выполнялось с применением математических моделей, основанных на допущении отсутствия диссоциации. Однако, в связи с разработкой проектов перспективных двигателей со значительно увеличенной температурой продуктов сгорания, приближающихся к 2000 К в основной камере сгорания и более 2100 К в форсажной, возникла необходимость учёта факторов диссоциации, рекомбинации и неравновесности на параметры рабочего процесса двигателя.
Одним из направлений совершенствования современных математических моделей ГТД, используемых при проектировочных термодинамических расчётах, является введение расчёта температуры и термодинамических свойств гомогенной смеси продуктов сгорания с учётом термической диссоциации и рекомбинации компонентов.
Основные положения
Для расчёта свойств рабочего тела в настоящее время применяется методика ЦИАМ [3], основывающаяся на ряде основных допущений:
— не учитывается влияние термической диссоциации газов (далее — диссоциации) —
— сгорание топлива — полное при коэффициенте избытка воздуха а& gt-1,0, продукты сгорания — нереагирующая смесь СО2, Н2О (г), О2 и атмосферного азота, объёмный состав которых зависит только от величины, а и состава топлива. При этом критерий Дамкёлера Da, представляющий собой отношение характерного времени физического процесса к характерному времени химической реакции, равен нулю-
— рабочие тела (воздух и продукты сгорания) представляют собой смесь компонентов, обладающую свойствами идеального газа, с неизменными термодинамическими свойствами, зависящими только от температуры.
Методика учёта фактора диссоциации была создана в 1970-х гг. [4] для расчётов ракетных двигателей. Теоретические основы методики базируются на более строгих, чем в [3], допущениях:
— наиболее вероятному, т. е. устойчивому равновесному состоянию изолированной термодинамической системы, соответствует состав, для которого при заданных исходных данных (суммарном давлении p, температуре T и полной энтальпии /) энтропия? будет максимальной. Такому состоянию отвечает соотношение парциальных давлений газообразных компонентов, задаваемое значениями констант равновесия обратимых химических реакций, идущих равновесно между всеми входящими в продукты сгорания индивидуальными веществами. При этом условиями устойчивости будут условия ё?& lt-0,ср=0-
— продукты сгорания — смесь химически реагирующих газов, состав и объёмное содержание которых определяется с учётом диссоциации и рекомбинации при постоянном давлении по уравнениям химического равновесия и баланса масс химических элементов (при этом значение Da стремится к бесконечности) —
— компоненты продуктов сгорания обладают свойствами идеального газа, их изобарная теплоёмкость зависит только от температуры, однако изобарная теплоёмкость смеси в целом зависит также и от уровня давления и тем значительнее, чем интенсивнее процесс диссоциации.
Для адаптации методики к расчёту рабочих процессов ГТД необходимо внести следующие дополнения:
— набор учитываемых химических элементов целесообразно расширить до следующего: учитываются шесть химических элементов С, Н, К, О, Аг, Б с увеличением количества учитываемых индивидуальных веществ до двадцати: Н, Н2, ОН, Н20, НО2, Н2О2, О,
О2, С, СО, С02, К, N, N0, К02, К20, Аг, Б, БО, Б02, наиболее вероятно образующихся при сгорании углеводородного топлива в воздушной среде (вещества Б, Б0 и Б02 введены для возможности расчёта с топливом, содержащим серу). При этом не учитываются вещества групп (КН)х, (НК0)х, (СН0)х, обладающие незначительными объёмными долями (менее 1×10−7) и образующиеся при а& lt-1-
— в определение энтальпии топлива необходимо внести коэффициент, учитывающий потери располагаемого тепла топлива (степень неполноты сгорания топлива) —
— для расчёта форсажной камеры сгорания необходимо внести учёт продуктов сгорания в окислителе-
— для реализации расчёта процессов сжатия, расширения и смешения рабочего тела методику необходимо дополнить алгоритмами нахождения термодинамических параметров смеси по температуре, энтальпии или энтропии смеси (при известном составе топлива, коэффициенте избытка воздуха, а и давлении р).
Для предлагаемой методики выбор химических реакций является условным и произвольным в том смысле, что компоненты продуктов сгорания могут получаться в результате различных химических реакций, в том числе и той, которая выбрана в качестве единственной, тогда как в камере сгорания значительная часть химических превращений происходит в виде многостадийных разветвлённых реакций. Каждой выбранной реакции соответствует константа равновесия с известной зависимостью от температуры [5]. Причём предполагается, что зависимости сохраняются и в случаях, когда реакция идёт в присутствии реагентов, участвующих в других реакциях.
Это предположение соответствует правдоподобной гипотезе о том, что динамическое равновесие смеси всех газов достигается при установлении динамического химического равновесия для каждой из перечисленных реакций. Кроме того, согласно закону Г есса, тепловой эффект реакции (в нашем случае — реакции горения) не зависит от промежуточных стадий и определяется лишь начальным и конечным состоянием системы. Т. е. если из заданных исходных продуктов можно различными путями получить заданные конечные продукты, то независимо от вида промежуточных реакций суммарный тепловой эффект будет одним и тем же (при условии протекания процесса при постоянном объёме или давлении).
Система уравнений для расчёта термодинамических параметров
Система уравнений, описывающая термодинамическое состояние продуктов сгорания, состоит из следующих групп:
1. Уравнения химического равновесия для наиболее вероятных при 1800… 2600 К химических реакций сгорания углеводородного топлива в воздушной среде с учётом диссоциации и рекомбинации. При этом для получения величин парциальных давлений веществ используются константы химического равновесия (четырнадцать уравнений по количеству входящих в систему химических элементов и индивидуальных веществ). Приняты следующие уравнения реакций равновесия:
С02~С0+0,502 К, (тп*с)= РсоР°2 — (1)
Р.
СО 2
р р0,5
Н2°~Н2+0,5°2 к2 (тп*с)= и^°2 — (2)
7* ГН^О
-IV ПС
2
р Н2°
Н2°~0Н+0,5Н2 к3 (тп*с) = р°НрН2 — (3)
р Н2°
0,5Н2+°2^Н0?
Н2+°2^Н2°2
к 4 (с)=(4& gt-
рН р0,
к5 (тп*с)= рН202 — (5)
р н р
0,5К2+02~№ 2 к 6 (тп*с) = ^ - (6)
р ^ р0,
К2+0,502^К20 N2+02^^ 2N0 С0~С+0
К 7 ТЛе) =
К, (С)_
р0−5 р N2 р N0
Р02 РN
*)_ рс р0 •
9 т ПС у
РС0
к9 Тс)_
во^в+о К, 0 (гп'-с)_
в02"в+02 К"(тП*С)_ РР0'
К, 2 (С)
*)_ рв р0.
ПС /
Рв0
И2^2И
0220
рв02
)_ рн •
— ПС ,
р н,
2
К, 3 (гп*с)_
Ро,
К14 (С)
р N
Р N2
(7)
(8) (9)
(10)
(11)
(12)
(13)
(14)
где Кг (тПс) — константы равновесия обратимых химических реакций, значения которых в различных формах представления приведены в [5].
В дальнейшем с целью сокращения записи уравнений заменим К ((г1*С) на К
2. Уравнения свойств смеси идеальных газов:
— уравнение состояния:
ж п + п
(15)
0 П0К + пт *
Xт, р, _-^Ъ-тРх
/_1 Кг
ЬПС
где КПС — количество молей продуктов сгорания-
ПОК — расход окислителя в кг/с на входе в камеру сгорания, для основной камеры сгорания он равен расходу воздуха, для форсажной — расходу газа-
ПТ — расход топлива в кг/с,
р- - парциальное давление /-го компонента, МПа-
*
р х — суммарное давление смеси продуктов сгорания, МПа-
— уравнение Дальтона (как следствие закона сохранения импульса):
20
(16)
/ _1
3. Уравнения материального баланса химических элементов (шесть уравнений по количеству входящих в систему химических элементов):
— баланс углерода
п0к + сПт *- (17)
РС0 +РС0 +рС _ ~ Рх — ()
тс К ПС
— баланс кислорода
р0н + рн20 + 2 рн02 + 2 рн202 + р0 + 2 р02 + РС0 + 2 РС02 + рж + 2 р N02 + РN20 +
п0К+0Пт
т0 К ПС
Рх-
N0
— баланс водорода
Рн + 2 рн2 + Рая + 2 рн2о + рно2 + 2 рн2о2
— баланс азота
О^ + НОт
тн К ПС
.. оок + коТ
РN + 2рN + рN0 + р№о2 + 2рК2О = Тг
2 2 2
Рх —
Рх-
ПС
— баланс серы
— баланс аргона
Рз + Рзо + Р80, =
о Ок + їо т
тз К ПС
Р х-
Р Аг =
О,
Аг
ок
т Аг К ПС
Рх,
(19)
(20)
(21)
(22)
где ГС гн г0 /~* N ^ Лг
где п0К, п0К, п0К, п0К, п0К, п0К
серы и аргона в составе окислителя, кг/с,
массы углерода, водорода, кислорода, азота,
Оок = О оск + О ок + О ок + О ок + О ок + О оАк
ве,
С, О, Н, N, S — массовые доли углерода, кислорода, водорода, азота и серы в топли-
С + О + Н + N + Ї = 1,0
тС, т0, тн, т^ тв, тЛг — массовые числа углерода, кислорода, водорода, азота, серы и аргона.
Поскольку окислителем может быть как атмосферный воздух, так и смесь воздуха с продуктами сгорания, то последний вариант принимается как более общий, для формирования системы уравнений.
4. Уравнение первого закона термодинамики, выраженное через полные энтальпии топлива, окислителя и продуктов сгорания (уравнение сохранения энергии):
Пт • /Т Тт)+ П0К • /0к Тк МПж + Пт)• /Пс Тс)
(23)
где
і0 (Т0) =
*ПС Vі ПС /
20
І I (тПс)р
І =1
І т іРі
і=1
В этом уравнении /(Т*ПС) — полная энтальпия /-го компонента (с учётом энтальпии образования вещества), кДж/кмоль,
/0Т (ГТ), /00К (Т 0К), /0ПС (Т ПС) — полные энтальпии топлива, окислителя и продуктов сгорания, взятые для соответствующих температур, кДж/кг.
При этом из физического смысла следуют ограничения на знак неизвестных величин — парциальных давлений и температуры:
Р/& gt-0 и Т*пс& gt-0
Задача определения равновесного состояния системы сводится к решению системы нелинейных алгебраических уравнений (1) — (23), содержащих в качестве неизвестных значения парциальных давлений двадцати компонентов продуктов сгорания, полной температуры, энтальпии и кажущейся молекулярной массы смеси продуктов сгорания. Таким образом, число уравнений равно 23, число неизвестных равно 23, при задании условий на входе в камеру сгорания (температур топлива и окислителя и суммарного давления) и соответствующих граничных условий решение системы, если оно существует, является единственным.
Необходимая для расчёта полная энтальпия топлива при температуре топлива 298,15 К определена следующим образом.
При температуре стандартного состояния Т0 полная энтальпия топлива 1 Т (Т0) равна теплоте его образования. Согласно закону Гесса, тепловой эффект реакции Ни^ КС при неполном сгорании (^ КС & lt- 1,0)и соответствующая ему полная энтальпия топлива вычисляется из уравнения
/Т (298,15) = Ни Лкс +
ДНСо ±^~АНно ±^-АН±^АН80о Л V 12,01 2 2,016 н20 14,007 2 32,066 802 у
где АНСо = -39 351С, 0 кДж/кг, АН0 =-241 814,0 кДж/кг, АН^ = 34 192,0
кДж/кг, АН80о = -296 810,0 кДж/кг — энтальпии образования СО2, н20(г), К02, в02 при
Т0 = 298,15 К по данным [3] с учётом содержания в топливе азота и серы.
Для определения термодинамических параметров смеси добавим вычисление изобарной теплоёмкости, энтропии и газовой постоянной по известным формулам, что дополняет систему несколькими уравнениями и таким же количеством формальных переменных, что не усложняет систему в целом.
Таким образом, после решения системы уравнений мы получаем все термодинамические параметры рабочего тела, необходимые для дальнейшего термодинамического расчёта ГТД.
Необходимость учёта неравновесности состава в зависимости от времени пребывания рабочего тела в объёме узла
При расширении рабочего тела в турбине и реактивном сопле ГТД происходит снижение температуры потока. При этом, если в конечной точке расширения моделировать равновесное состояние для конечной температуры, то необходимо понимать, что для этой точки полностью завершены все рассматриваемые реакции диссоциации и рекомбинации. Принимая во внимание очень малое время пребывания рабочего тела в объёме узла, значительно меньшее, чем скорость протекания многих из рассматриваемых реакций (1−14), необходимо параметры турбины и сопла моделировать по «замороженному» состоянию на входе, т. к. на выходе узла равновесное состояние не успеет установиться, и его моделирование по конечной температуре расширения приведёт к значительным погрешностям. А для реакций с наибольшим временем протекания (1, 2, 3, 6) «замороженное» состояние необходимо принимать по параметрам на выходе из камеры сгорания.
Это означает, что расчёт равновесного состояния происходит для процессов сгорания топлива в основной и форсажной камерах сгорания, а далее для расчёта процессов, происходящих в турбине, смесителе и сопле, расчёт происходит по частично «замороженным» химическим реакциям.
Введение учёта кинетики химических реакций проведём следующим способом.
Пусть за камерой сгорания стоит узел ГТД с индексом У1. При известных размерах узла У1 и скорости потока в осевом направлении определяется время пребывания рабочего тела в пределах узлаУ1, следующего за камерой сгорания. Также при данной температуре известно время прохождения каждой из рассматриваемых реакций ^ или скорости
реакции -^7. Если выполняется условие
*У1 & lt-р 1,
то реакция считается незавершённой и «замороженной» по условиям на входе в узел У1. Далее за узлом У1 следует узел У2. Для него также проверяется условие
(Су1 + Су2) & lt- Ьр1
при температуре на выходе узла У1. Если условие выполняется, реакция продолжает считаться «замороженной», если же условие не выполняется, реакция считается завершённой и включается в число рассматриваемых в равновесном состоянии реакций. Общий
принцип проиллюстрирован на рисунке 1. Сведения по скоростям протекания реакций приведены в справочной литературе.
Равновесное состо-
К
Неравновесное состояние, часть реакций «замо-
Неравновесное состояние, часть реакций перешла в равновесное состояние, часть самых «медленных» реакций продолжает оста-
Рисунок 1. Принцип учёта неравновесности состава смеси в зависимости от скорости реакций и времени пребывания рабочего тела в объёме узла
После определения параметров равновесного состояния (1−23) происходит поправка термодинамических параметров с учётом неравновесности по указанным выше условиям.
Поскольку уравнения химических реакций, протекающих равновесно, учитывают как диссоциацию, так и рекомбинацию, то в результате расчёта равновесного состояния учитываются оба этих фактора, а фактор неравновесности состава учитывается «замораживанием» отдельных реакций. При этом степень влияния термической диссоциации и рекомбинации будет выше для двигателей большей размерности, для которого объём узлов сравнительно велик и реакции диссоциации и рекомбинации протекают более полно.
Верификация математической модели
Алгоритм решения предложенной системы уравнений реализован в виде отдельной программы БАвво, позволяющей проводить тепловые расчёты основной и форсажной камер сгорания. Методика верифицирована по результатам расчёта процесса горения в камерах сгорания ракетных двигателей для сжигания керосина в воздушной среде [1] и показала хорошую сходимость результатов по определению объёмных долей и энтальпии смеси с погрешностью, не превышающей 0,3% в диапазоне температур от 1800 до 2800 К [2] (рисунок 2). При этом общее время расчёта термодинамических параметров ГТД в целом увеличивается несущественно.
С целью обеспечения сравнения полученных результатов расчёта по рассматриваемой методике с опубликованными в справочной и научной литературе данными, термодинамические свойства индивидуальных веществ и сведения по константам равновесия для химических реакций приняты по [5]. Однако расчёт может быть выполнен с любыми более современными данными.
Рисунок 2. Сравнение значений истинной изобарной «замороженной» (пунктир) и равновесной (сплошная линия) теплоёмкости, рассчитанной по методике [3] (точки) и Б1б-
БО
Выводы
На основе проверенной методики В. Е. Алемасова создана адаптированная к расчётам рабочих процессов ГТД математическая модель расчёта термодинамических параметров рабочего тела с учётом фактора термической диссоциации и рекомбинации в гомогенной смеси продуктов сгорания углеводородного топлива.
Созданная математическая модель позволяет проводить расчёт процесса сгорания топлива с определением величин температуры, энтальпии, энтропии и теплоёмкости смеси с использованием констант химического равновесия. При реализации возможности расчёта термодинамических параметров смеси по заданной температуре, энтальпии или энтропии смеси обеспечивается реализация набора необходимых функций для расчёта термодинамических процессов сжатия, расширения и смешения рабочего тела в узлах ГТД.
При реализации математической модели учитывается фактор неравновесности, который заключается в возможности «замораживания» отдельных химических реакций для расчётов узлов ГТД с малым, по сравнению со скоростью протекания реакции, временем пребывания рабочего тела в объёме узла. Таким образом, степень влияния термической диссоциации и рекомбинации будет различной для двигателей разной размерности.
Применение математической модели в составе системы моделирования рабочих процессов ГТД позволяет проводить термодинамический расчёт ГТД и его узлов с учётом термической диссоциации и рекомбинации.
Список литературы
1. Алемасов В. Е. и др. Термодинамические и теплофизические свойства продуктов сгорания: В 5 т. — М.: ВИНИТИ, 1971−1974.
2. Болдырев О. И. Оценка влияния термической диссоциации на температуру продуктов сгорания углеводородного топлива в основной и форсажной камерах сгорания авиационных ГТД. Научно-технический отчёт о НИР ОАО «НПП «Мотор» № 199ДО-017. -Уфа, ОАО «НПП «Мотор», 2010.
3. Руководящий технический материал авиационной техники РТМ 1677−83: Двигатели авиационные газотурбинные: Методы и подпрограммы расчёта термодинамических параметров воздуха и продуктов сгорания углеводородных топлив. — М.: ЦИАМ, 1983. — 92 с.
4. Теория ракетных двигателей: Учебник для вузов / В. Е. Алемасов, А. Ф. Дрегалин, А.П. Тишин- / Под ред. В. П. Глушко. — М.: Машиностроение, 1989. — 464 с.
5. Термодинамические свойства индивидуальных веществ: Справочник: В 4 т. / Под ред. В. П. Глушко. — М.: Наука, 1979−1982.
Рецензенты:
Гишваров А. С., д.т.н., профессор, зав. кафедрой авиационных двигателей ФБГОУ ВПО Уфимский государственный авиационный технический университет, г. Уфа.
Кривошеев И. А., д.т.н., профессор, декан факультета авиационных двигателей ФБГОУ ВПО Уфимский государственный авиационный технический университет, г. Уфа.

ПоказатьСвернуть
Заполнить форму текущей работой